• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于數(shù)值模擬研究瀝青質(zhì)沉積對(duì)油田開發(fā)的影響

    2022-10-27 09:30:52TANANYKHINSTRUCHKOVKHORMALIROSCHIN
    石油勘探與開發(fā) 2022年5期
    關(guān)鍵詞:油藏組分原油

    TANANYKHIN D S,STRUCHKOV I A,KHORMALI A,ROSCHIN P V

    (1.俄羅斯圣彼得堡礦業(yè)大學(xué),圣彼得堡 199106,俄羅斯;2.獨(dú)聯(lián)體秋明石油研究中心,秋明 625000,俄羅斯;3.伊朗岡巴德卡沃斯大學(xué),岡巴德卡沃斯 49717-99151,伊朗;4.獨(dú)聯(lián)體 SamaraNIPIneft公司,薩馬拉 443010,俄羅斯)

    0 引言

    在油田長(zhǎng)期開發(fā)過(guò)程中,原油中瀝青質(zhì)的溶解能力可能會(huì)降低,從而產(chǎn)生絮凝。在滲流過(guò)程中,瀝青質(zhì)沉積在巖石的多孔介質(zhì)表面,減小了孔隙通道水力半徑,并產(chǎn)生了滲流阻力,從而導(dǎo)致地層傷害[1]。通過(guò)在組分模擬器中激活瀝青質(zhì)選項(xiàng),可以模擬上述過(guò)程。在模擬油藏不同生產(chǎn)方式時(shí),使用該模擬器可以考慮原油中熱力學(xué)平衡的破壞,考慮油田生產(chǎn)過(guò)程中由于原油組分和壓力-溫度條件的變化而引起的原油中瀝青質(zhì)動(dòng)態(tài)特征的不確定性。了解原油中瀝青質(zhì)的相變條件,可以合理選擇油藏生產(chǎn)方案以及油井作業(yè)條件,最大限度降低生產(chǎn)風(fēng)險(xiǎn)。先后有多位研究者利用組分模擬器進(jìn)行了瀝青質(zhì)沉積模擬,提出了多個(gè)描述瀝青質(zhì)行為的模型,研究了瀝青質(zhì)沉積造成地層傷害的程度和機(jī)理[2-3]。

    本文使用組分模擬器,基于原油與巖心樣品的實(shí)驗(yàn)室研究數(shù)據(jù),對(duì)俄羅斯某油田多孔介質(zhì)中的瀝青質(zhì)沉積進(jìn)行模擬,并進(jìn)行敏感性分析和不確定性評(píng)價(jià)。進(jìn)行不確定性評(píng)價(jià)時(shí),使用模擬器中瀝青質(zhì)選項(xiàng)的參數(shù),以評(píng)估瀝青質(zhì)選項(xiàng)的所有參數(shù)對(duì)油田產(chǎn)量的綜合影響。

    1 油藏模型

    本文以俄羅斯某油田為例開展研究,該油田已開發(fā)近46年。油田開發(fā)的第1階段為密集鉆井開發(fā),共鉆240口井,其中216口為生產(chǎn)井,24口為注入井。這些井部署在含油厚度最大的位置,沿油田主軸排成3排,井網(wǎng)尺寸為400 m×500 m。油藏模型建立過(guò)程描述如下。

    1.1 靜態(tài)模型

    建立了研究區(qū)的靜態(tài)模型,其沿主軸的長(zhǎng)度為19 km,沿次軸的長(zhǎng)度為5 km。沿著生產(chǎn)層位的頂面和底面建立油田靜態(tài)模型的主要構(gòu)造框架。根據(jù)綜合測(cè)井?dāng)?shù)據(jù)識(shí)別油藏的地層邊界,包括中子測(cè)井、自然電位測(cè)井和電阻率測(cè)井。將測(cè)井圖與相鄰井的數(shù)據(jù)進(jìn)行了比較,根據(jù)測(cè)井曲線的典型結(jié)構(gòu),考慮到研究區(qū)的沉積地層序列和相構(gòu)造特征,對(duì)相關(guān)夾層進(jìn)行了研究。在考慮油田面積和平均井距的情況下,選擇了網(wǎng)格單元的線性尺寸,滿足了在平面上相鄰油井之間至少有5個(gè)網(wǎng)格單元的條件。平面上網(wǎng)格單元尺寸為 50 m×50 m,垂向網(wǎng)格的數(shù)量及尺寸根據(jù)滲透層的最小厚度(0.5 m)來(lái)確定。根據(jù)測(cè)井?dāng)?shù)據(jù)解釋得到的孔隙度曲線得到井點(diǎn)位置孔隙度,使用3D模型對(duì)井?dāng)?shù)據(jù)進(jìn)行插值得到孔隙度的空間分布??紤]油水界面的位置,利用從測(cè)井?dāng)?shù)據(jù)中獲得的含油飽和度曲線計(jì)算含油飽和度空間分布。油藏的平均孔隙度為 19%,平均滲透率0.8 μm2,平均原始含油厚度為30 m。

    1.2 動(dòng)態(tài)模型

    地質(zhì)模型包含了全部油層和部分水層。通過(guò)增大油藏邊界附近飽和水網(wǎng)格的孔隙體積來(lái)模擬水驅(qū)。在建立油田動(dòng)態(tài)模型時(shí),使用了角點(diǎn)法。通過(guò)對(duì)靜態(tài)模型的垂向網(wǎng)格進(jìn)行粗化得到動(dòng)態(tài)模型網(wǎng)格。結(jié)合層的劃分,確定了垂向上層的數(shù)量及大小。由于該油田普遍發(fā)育低滲夾層,因此粗化網(wǎng)格時(shí)將這些夾層劃分為獨(dú)立網(wǎng)格。動(dòng)態(tài)模型網(wǎng)格單元的尺寸為50 m×50 m×3.6 m。全油田動(dòng)態(tài)模型由1 100 000多個(gè)網(wǎng)格組成。通過(guò)粗化賦值得到動(dòng)態(tài)模型網(wǎng)格屬性。從全油田模型中截取了一個(gè)位于油田中心、尺寸為1.5 km×1.5 km的井區(qū)以減少計(jì)算時(shí)間,該井區(qū)模型由33 000個(gè)網(wǎng)格組成,模型僅包括4口生產(chǎn)井(見圖1)。調(diào)整了模型中的含水層,以恢復(fù)油藏動(dòng)態(tài)壓力和產(chǎn)量。

    圖1 動(dòng)態(tài)模型

    研究水錐進(jìn)、瀝青質(zhì)沉積等復(fù)雜動(dòng)態(tài)過(guò)程,需要對(duì)局部區(qū)域網(wǎng)格進(jìn)行大規(guī)模離散??紤]到計(jì)算所需時(shí)間,可能無(wú)法對(duì)整個(gè)模型進(jìn)行網(wǎng)格加密。此外,也沒(méi)必要在受這些過(guò)程影響最小或不受影響的區(qū)域進(jìn)行網(wǎng)格加密。因而采用 Heinemann等[4]提出的局部網(wǎng)格加密方法進(jìn)行加密。從未水侵區(qū)域(原始網(wǎng)格)開始加密,每個(gè)后續(xù)單元對(duì)相鄰前一個(gè)單元的加密倍數(shù)不超過(guò)3倍,直到達(dá)到研究區(qū)所需的網(wǎng)格尺寸[5]。此規(guī)則是為了在計(jì)算原始網(wǎng)格和加密網(wǎng)格之間的流動(dòng)時(shí)達(dá)到所需的精度。

    在每口生產(chǎn)井區(qū)域進(jìn)行局部網(wǎng)格加密(見圖2)。平面上,在距離生產(chǎn)井75~125 m處,新網(wǎng)格大小為25 m×25 m,在距離生產(chǎn)井25~75 m處,新網(wǎng)格大小為12.5 m×12.5 m,在距離生產(chǎn)井0~25 m處,新網(wǎng)格大小為6.25 m×6.25 m。垂向上,將原始網(wǎng)格一分為二,即新網(wǎng)格垂向尺寸為1.8 m。

    圖2 生產(chǎn)井周圍的局部加密網(wǎng)格

    動(dòng)態(tài)模型中使用的相對(duì)滲透率曲線如圖3所示。

    圖3 相對(duì)滲透率曲線

    1.3 PVT模型

    根據(jù) PVT(壓力-體積-溫度)實(shí)驗(yàn)結(jié)果,利用Soave-Redlich-Kwong狀態(tài)方程建立了研究區(qū)流體PVT模型。該模型采用Lorenz-Bray-Clark黏度計(jì)算方法。根據(jù)實(shí)驗(yàn)結(jié)果調(diào)整了臨界溫度、臨界壓力、偏心因子、體積轉(zhuǎn)換系數(shù)、二元交互作用系數(shù)等參數(shù)。PVT模型模擬值與實(shí)驗(yàn)值的誤差小于1%。

    PVT模型中的原油組分被劃分為7個(gè)擬組分:N2和C1歸為組分X1,因?yàn)樗鼈兙哂邢嗨频臒崃W(xué)性質(zhì);CO2、C2和C3歸為組分X2,因?yàn)镃O2表現(xiàn)出與輕質(zhì)油中間組分相似的特征;C4—C6歸為組分 X3;C7—C14歸為組分 X4;C15—C21歸為組分 X5;C22+劃分為非沉積組分 X6nA(高相對(duì)分子質(zhì)量石蠟)和沉積組分 X6A(瀝青質(zhì))。各組分性質(zhì)如表1所示。與非沉積組分相比,瀝青質(zhì)組分具有較高的二元相互作用系數(shù)(見表1),這解釋了其與輕質(zhì)油組分互溶性較低的原因。

    表1 各擬組分性質(zhì)

    恒組成膨脹(CCE)實(shí)驗(yàn)結(jié)果如圖4所示。油藏壓力為20 MPa,泡點(diǎn)壓力為4.7 MPa,油藏條件下的原油黏度為2.6 mPa·s,地面條件下的原油密度為840 kg/m3。

    圖4 恒組成膨脹實(shí)驗(yàn)結(jié)果

    筆者評(píng)估了瀝青質(zhì)初始析出壓力[6],得出其值為6.8 MPa,據(jù)此對(duì)PVT模型中的瀝青質(zhì)絮凝模型進(jìn)行了調(diào)整。壓力低于瀝青質(zhì)初始析出壓力時(shí),原油中的瀝青質(zhì)開始絮凝,隨著壓力進(jìn)一步降低,絮凝的瀝青質(zhì)尺寸增大并沉積在固體表面,從而導(dǎo)致采油困難。因此,建議將油藏壓力保持在瀝青質(zhì)初始析出壓力以上。

    1.4 瀝青質(zhì)選項(xiàng)

    Ghadimi等[7]指出,根據(jù)油田的生產(chǎn)條件、規(guī)模、計(jì)算周期和瀝青質(zhì)絮凝、沉積相關(guān)參數(shù),在油藏組分模擬器中激活瀝青質(zhì)選項(xiàng),會(huì)使預(yù)測(cè)的原油采收率平均降低 20%~35%。本文根據(jù)筆者前期巖心驅(qū)替實(shí)驗(yàn)結(jié)果(見圖 5)和巖心孔喉直徑分布研究結(jié)果(見圖6)[6],對(duì)瀝青質(zhì)在多孔介質(zhì)中的沉積模型進(jìn)行了調(diào)整。將根據(jù)實(shí)驗(yàn)結(jié)果調(diào)整后的代表瀝青質(zhì)在多孔介質(zhì)中沉積的系數(shù)作為基準(zhǔn)參考值。

    圖5 注入孔隙體積倍數(shù)與壓差關(guān)系[6](階段1、階段5為煤油驅(qū)替,階段2—4為原油驅(qū)替)

    圖6 含懸浮瀝青質(zhì)原油驅(qū)替前后的孔喉直徑分布[6]

    在模擬器中,以簡(jiǎn)化的形式設(shè)置原油中瀝青質(zhì)的狀態(tài)。通常情況下,以表格形式列出原油中懸浮瀝青質(zhì)摩爾分?jǐn)?shù)(懸浮瀝青質(zhì)的物質(zhì)的量與原油中瀝青質(zhì)的總物質(zhì)的量之比)與壓力、溫度、相關(guān)原油組分或化學(xué)劑的摩爾分?jǐn)?shù)的相互關(guān)系,這些參數(shù)的相互關(guān)系可以通過(guò)實(shí)驗(yàn)研究獲得。這種方法需要的設(shè)計(jì)能力和時(shí)間成本最小,且計(jì)算精度足夠高。通過(guò)專用程序創(chuàng)建相互關(guān)系表,用于建立PVT模型。根據(jù)特定實(shí)驗(yàn)結(jié)果或動(dòng)態(tài)模型的歷史擬合結(jié)果,在模擬器中設(shè)定瀝青質(zhì)吸附、解吸和堵塞系數(shù)。

    計(jì)算多孔介質(zhì)中瀝青質(zhì)沉積量的流程如下:①計(jì)算每個(gè)網(wǎng)格的壓力。②求解流動(dòng)的質(zhì)量守恒方程組,確定每個(gè)有效網(wǎng)格中各組分的物質(zhì)的量。③將網(wǎng)格壓力與懸浮瀝青質(zhì)摩爾分?jǐn)?shù)-壓力關(guān)系表對(duì)比,確定不同網(wǎng)格中的固相濃度。④計(jì)算各時(shí)間步的瀝青質(zhì)沉積速率。⑤將吸附、解吸和堵塞系數(shù)引入瀝青質(zhì)沉積模型,使用油藏含油飽和度和局部孔隙度。⑥計(jì)算滲流速率。⑦根據(jù)瀝青質(zhì)動(dòng)態(tài)沉積速率計(jì)算多孔介質(zhì)中瀝青質(zhì)沉積物體積。⑧從初始孔隙體積中減去瀝青質(zhì)沉積物體積,得到新的局部孔隙度,然后根據(jù)孔隙度確定滲透率。如果滲流速率超過(guò)臨界速率,則部分沉積物從孔隙表面脫離,并與原油一起通過(guò)孔隙通道流向生產(chǎn)井井底。⑨重新計(jì)算含油飽和度和壓力。⑩重復(fù)上述步驟直至計(jì)算結(jié)束。

    2 不確定性評(píng)估

    進(jìn)行油田開發(fā)時(shí),為了在開始全面投資之前評(píng)估項(xiàng)目風(fēng)險(xiǎn),需對(duì)認(rèn)識(shí)不足的資產(chǎn)進(jìn)行初步評(píng)估。多變量計(jì)算(MVC)方法能夠在不確定性條件下對(duì)投資項(xiàng)目進(jìn)行全面評(píng)估,以做出最佳管理決策并將財(cái)務(wù)損失降至最低。在MVC方法中,不確定性可以理解為由于儀器、方法、研究技術(shù)、覆蓋范圍及其他原因造成的測(cè)量誤差,參數(shù)在較大范圍內(nèi)變化。各種不確定性可能對(duì)最終結(jié)果產(chǎn)生不同的影響[8]。在這種情況下,不確定性評(píng)估是將原始數(shù)據(jù)的不確定性轉(zhuǎn)換為最終結(jié)果的不確定性。MVC方法的關(guān)鍵目標(biāo)之一是研究不確定性范圍,評(píng)估油氣生產(chǎn)的所有可能方案以及不確定性參數(shù)對(duì)產(chǎn)量和累計(jì)產(chǎn)量的影響。油田鉆井時(shí)會(huì)進(jìn)行取心分析,以及地層流體分析、測(cè)井和試井。這些研究資料將構(gòu)成動(dòng)態(tài)模型中的原始數(shù)據(jù),其中每個(gè)參數(shù)都有一定的誤差。將原始數(shù)據(jù)加載到模型后,有些參數(shù)的誤差可能達(dá)到 20%。對(duì)研究區(qū)原始數(shù)據(jù)進(jìn)行了研究,根據(jù)測(cè)量?jī)x器誤差、研究方法精度以及本次研究范圍,確定了不確定性參數(shù),估計(jì)了其變化范圍及分布函數(shù)。

    由于原始數(shù)據(jù)中不確定性參數(shù)較多,需要進(jìn)行敏感性分析,包括優(yōu)化輸入?yún)?shù)的數(shù)量。具體而言,需要確定對(duì)目標(biāo)函數(shù)影響最大的不確定性參數(shù),評(píng)估其相互影響,并通過(guò)特殊設(shè)計(jì)的實(shí)驗(yàn)方法進(jìn)行相關(guān)性計(jì)算,排除冗余的影響較小參數(shù)。使用單變量(OVAT)方法[9]開展敏感性分析,該方法可以評(píng)估輸入?yún)?shù)和響應(yīng)參數(shù)之間的相互作用,并確定影響最大的參數(shù)。

    敏感性分析之后,繼續(xù)進(jìn)行不確定性評(píng)估,直接進(jìn)行多元計(jì)算。設(shè)計(jì)多元計(jì)算實(shí)驗(yàn)的常用方法是McKay等[10]于 2000年提出的拉丁超立方體法。使用該方法可以獲得形式復(fù)雜(通常為多項(xiàng)式)的響應(yīng)面,也稱為代理模型,來(lái)描述整個(gè)系統(tǒng)的變化規(guī)律。該方法可以用較少的實(shí)驗(yàn)次數(shù)計(jì)算出所有參數(shù)空間。在這種情況下,可以統(tǒng)一考慮所有不確定性參數(shù)的可能值,即使用參數(shù)整個(gè)取值范圍,避免出現(xiàn)“空白”區(qū)域。假設(shè)實(shí)驗(yàn)次數(shù)為N,則每個(gè)參數(shù)的取值范圍被劃分為N個(gè)等概率級(jí)別。然后,N次實(shí)驗(yàn)隨機(jī)地用以下限制填充參數(shù)空間:每次實(shí)驗(yàn)在各級(jí)別上占據(jù)一個(gè)隨機(jī)位置,每個(gè)級(jí)別上只進(jìn)行一次實(shí)驗(yàn)。設(shè)置實(shí)驗(yàn)次數(shù)時(shí),應(yīng)足以構(gòu)建一個(gè)完整的代理模型。

    實(shí)施多元計(jì)算的算法包括以下步驟:①創(chuàng)建不確定性參數(shù)矩陣,包含一組輸入?yún)?shù)及其變化范圍和概率分布,并創(chuàng)建目標(biāo)函數(shù)。②使用低層次設(shè)計(jì)方法(OVAT方法)設(shè)計(jì)實(shí)驗(yàn),并進(jìn)行計(jì)算。③構(gòu)建響應(yīng)面,繪制帕累托圖,以確定最重要的輸入?yún)?shù)(主控因素)。④調(diào)整不確定性參數(shù)的范圍,去除對(duì)目標(biāo)函數(shù)影響較小的參數(shù)。⑤根據(jù)上一步定義的輸入?yún)?shù),使用高層次設(shè)計(jì)方法(拉丁超立方體法)設(shè)計(jì)實(shí)驗(yàn),并進(jìn)行計(jì)算。⑥重復(fù)步驟④、⑤,直到目標(biāo)函數(shù)最小化。⑦構(gòu)建響應(yīng)面,即代理模型。⑧盲選輸入?yún)?shù),采用代理模型進(jìn)行計(jì)算,并檢查代理模型的可靠性。⑨使用蒙特卡羅方法創(chuàng)建一組實(shí)驗(yàn)。使用步驟⑦中構(gòu)建的代理模型進(jìn)行計(jì)算。

    動(dòng)態(tài)模型應(yīng)成為實(shí)際油田的數(shù)字孿生模型,并充分模擬生產(chǎn)特征(所有油井動(dòng)態(tài)歷史數(shù)據(jù)),以便可靠地預(yù)測(cè)產(chǎn)量,這是歷史擬合的目標(biāo)。在這種情況下,一個(gè)動(dòng)態(tài)模型與另一個(gè)動(dòng)態(tài)模型的對(duì)比標(biāo)準(zhǔn)是目標(biāo)函數(shù)。因此,可將歷史擬合過(guò)程簡(jiǎn)化為目標(biāo)函數(shù)最小化問(wèn)題,即響應(yīng)參數(shù)的模擬值與實(shí)際值之間的差值最小化。響應(yīng)參數(shù)通常是油氣水產(chǎn)量、累計(jì)產(chǎn)量等。

    由于本文研究的油田生產(chǎn)歷史較長(zhǎng),進(jìn)行了大量巖心和流體研究,因此可以消除油藏地質(zhì)參數(shù)和流體性質(zhì)的不確定性,只考慮瀝青質(zhì)沉積引起的不確定性。首先,在不激活瀝青質(zhì)選項(xiàng)的情況下,使用全油田模型進(jìn)行常規(guī)歷史擬合。在歷史擬合的第 1階段,重新建立油藏壓力動(dòng)態(tài),并將整個(gè)油田的產(chǎn)液量差異(計(jì)算值與實(shí)際值之差)降至最低。為此,對(duì)輸入數(shù)據(jù)進(jìn)行了以下調(diào)整:①修正飽和水網(wǎng)格(含水層)孔隙體積;②校正產(chǎn)液量差異較大的主力產(chǎn)區(qū)的滲透率;③復(fù)查油井射孔層段。歷史擬合的第 2階段調(diào)整含水動(dòng)態(tài),即調(diào)整各相相對(duì)滲透率函數(shù),并得到較好的單井?dāng)M合程度。根據(jù)室內(nèi)驅(qū)油實(shí)驗(yàn)數(shù)據(jù),生成相對(duì)滲透率曲線,表征生產(chǎn)井的含水動(dòng)態(tài)。然后,使用從室內(nèi)實(shí)驗(yàn)結(jié)果中獲得的基本參數(shù)激活瀝青質(zhì)選項(xiàng)。從已擬合的全油田模型中截取一個(gè)井區(qū)(見圖 1),并在此基礎(chǔ)上模擬生產(chǎn)歷史,建立與相鄰井網(wǎng)模型的邊界條件。分析了模型對(duì)瀝青質(zhì)選項(xiàng)參數(shù)的敏感性,根據(jù)敏感性分析結(jié)果選擇了對(duì)目標(biāo)函數(shù)影響最大的參數(shù)。

    僅使用瀝青質(zhì)選項(xiàng)參數(shù)進(jìn)行多變量歷史擬合,以確定瀝青質(zhì)選項(xiàng)參數(shù)的變化范圍。多變量歷史擬合的目標(biāo)函數(shù)如(1)式所示。目標(biāo)函數(shù)中的響應(yīng)參數(shù)不應(yīng)太多,從而優(yōu)化歷史擬合時(shí)間,但響應(yīng)參數(shù)數(shù)量應(yīng)足以獲得滿意的擬合效果。根據(jù)拉丁超立方體法設(shè)計(jì)并進(jìn)行了500次計(jì)算。根據(jù)多變量歷史擬合的結(jié)果,確定參數(shù)模擬值與實(shí)際值的容許偏差在 10%以內(nèi)。劃分瀝青質(zhì)選項(xiàng)參數(shù)的相應(yīng)變化范圍,并在此基礎(chǔ)上通過(guò)拉丁超立方體法生成500個(gè)預(yù)測(cè)值,進(jìn)而得到預(yù)測(cè)剖面。

    通過(guò)多變量歷史擬合獲得了一組調(diào)整后的精度可接受的模型輸入?yún)?shù)。根據(jù)這些輸入?yún)?shù)進(jìn)行計(jì)算可得到一系列預(yù)測(cè)生產(chǎn)剖面。然后進(jìn)行概率分析,確定概率分布曲線上累計(jì)概率10%,50%,90%時(shí)對(duì)應(yīng)的響應(yīng)參數(shù)值,分別用P10,P50,P90表示。可見,使用MVC方法可以在不確定條件下評(píng)估決策風(fēng)險(xiǎn),而傳統(tǒng)方法只得到1個(gè)預(yù)測(cè)生產(chǎn)剖面,缺乏概率分析。

    3 結(jié)果和討論

    3.1 敏感性分析

    本文的不確定性參數(shù)矩陣由11個(gè)變量組成。根據(jù)文獻(xiàn)確定瀝青質(zhì)選項(xiàng)參數(shù)的變化范圍[7]。為了進(jìn)行敏感性分析,使用OVAT方法進(jìn)行了25次實(shí)驗(yàn),以確定對(duì)目標(biāo)函數(shù)影響最大的輸入?yún)?shù),并評(píng)估其相互影響。實(shí)驗(yàn)計(jì)算結(jié)果如圖 7的帕累托圖所示。帕累托圖中頻率低于5%的變量不會(huì)影響目標(biāo)函數(shù),可以排除。在所選變量中,對(duì)目標(biāo)函數(shù)影響最大的變量是 A3、A4、A5、G、N和S,使用這些參數(shù)進(jìn)行了多變量歷史擬合。

    圖7 研究區(qū)瀝青質(zhì)參數(shù)帕累托圖

    3.2 衰竭開采油田產(chǎn)量預(yù)測(cè)

    使用拉丁超立方體方法和敏感性分析階段選定的6個(gè)變量創(chuàng)建 500個(gè)實(shí)驗(yàn),模擬衰竭式開采的生產(chǎn)狀況。圖 8顯示了不考慮瀝青質(zhì)沉積時(shí)的產(chǎn)量預(yù)測(cè)結(jié)果和考慮瀝青質(zhì)沉積時(shí)的產(chǎn)量預(yù)測(cè)結(jié)果概率分布(用P10,P50,P90表示)。將產(chǎn)油量作為控制參數(shù),在歷史擬合階段,累計(jì)產(chǎn)油量與實(shí)際值的偏差為零,累計(jì)產(chǎn)氣量、累計(jì)產(chǎn)水量的P50值與實(shí)際值的偏差均小于5%,累計(jì)產(chǎn)氣量、累計(jì)產(chǎn)水量的P10和P90值與實(shí)際值的偏差均小于10%。

    圖8 衰竭式開采時(shí)的計(jì)算結(jié)果

    P90,P50,P10值對(duì)應(yīng)的瀝青質(zhì)沉積量依次增加。油藏的有效孔隙體積隨著瀝青質(zhì)沉積量的增加而減小,因此,平均油藏壓力的P90,P50,P10值依次增加。由于不考慮瀝青質(zhì)沉積情況下的總產(chǎn)液量低于考慮瀝青質(zhì)沉積情況下總產(chǎn)液量的P10值,不考慮瀝青質(zhì)沉積情況下的平均油藏壓力比考慮瀝青質(zhì)沉積情況下平均油藏壓力的P90,P50,P10值更高。不確定性評(píng)估結(jié)果表明,結(jié)合瀝青質(zhì)選項(xiàng)可以對(duì)油田生產(chǎn)參數(shù)進(jìn)行合理的風(fēng)險(xiǎn)評(píng)估。有瀝青質(zhì)沉積情況下累計(jì)產(chǎn)油量的P90,P50,P10值分別比無(wú)瀝青質(zhì)沉積情況下的累計(jì)產(chǎn)油量低4.0%,8.1%,12.3%。

    在與生產(chǎn)井有一定距離的地方,生產(chǎn)過(guò)程中壓力降至瀝青質(zhì)初始析出壓力以下,原油中的瀝青質(zhì)開始絮凝。當(dāng)瀝青質(zhì)向生產(chǎn)井移動(dòng)時(shí),就會(huì)吸附到孔隙空間表面,隨著壓力繼續(xù)降低,原油中瀝青質(zhì)絮凝速度加快,生產(chǎn)井附近瀝青質(zhì)沉積速率最高(見圖9)。

    圖9 衰竭開采條件下預(yù)測(cè)期末多孔介質(zhì)中吸附瀝青質(zhì)摩爾分?jǐn)?shù)的分布(P50)

    3.3 注水開采油田產(chǎn)量預(yù)測(cè)

    為了保持油藏能量并防止瀝青質(zhì)沉積向油藏深處擴(kuò)散,有必要進(jìn)行注水增壓,保持油藏壓力高于瀝青質(zhì)初始析出壓力,以實(shí)現(xiàn)長(zhǎng)時(shí)間的低傷害開發(fā)。進(jìn)行了不確定性評(píng)估,以分析注水系統(tǒng)對(duì)產(chǎn)量預(yù)測(cè)結(jié)果的影響。該分析基于使用拉丁超立方體方法創(chuàng)建的 500個(gè)實(shí)驗(yàn),得到對(duì)目標(biāo)函數(shù)影響最大的參數(shù)。注水系統(tǒng)以井區(qū)中心的 1口虛擬注水井為代表。該井在油藏的含水層段射孔,歷史擬合階段未注水,僅在預(yù)測(cè)時(shí)注水,注采比設(shè)為1,計(jì)算結(jié)果如圖10所示。歷史擬合階段采用衰竭式開采,隨著注水的開始,P90,P50,P10值對(duì)應(yīng)的瀝青質(zhì)沉積量依次增加,平均油藏壓力的P90,P50,P10值依次增加。考慮瀝青質(zhì)沉積情況下累計(jì)產(chǎn)油量的P90,P50,P10值分別比不考慮瀝青質(zhì)沉積情況下的累計(jì)產(chǎn)油量低4.8%,9.1%,13.3%。

    圖10 注水開采時(shí)的計(jì)算結(jié)果

    計(jì)算結(jié)果表明,多孔介質(zhì)中瀝青質(zhì)的沉積顯著降低了累計(jì)產(chǎn)油量。對(duì)比圖9和圖11可知,注水保持油藏壓力可以有效減緩原油中瀝青質(zhì)的絮凝,與衰竭式開采相比,注水后瀝青質(zhì)沉積體積更小,且沉積部位更接近油藏頂面,因?yàn)樽⑺沟撞坑筒貕毫Φ玫交謴?fù)。向油藏含水層注水可以使油藏壓力保持在瀝青質(zhì)初始析出壓力以上的時(shí)間更長(zhǎng)。壓力從初始油藏壓力降低至瀝青質(zhì)初始析出壓力經(jīng)歷時(shí)間越長(zhǎng),多孔介質(zhì)中瀝青質(zhì)沉積物的形成速度就越慢,最終觀察到的瀝青質(zhì)沉積物就越少。因此,注水使原油從油藏流向油井的滲流阻力減小,從而提高油井產(chǎn)能。需要說(shuō)明的是,本文認(rèn)為注入水不會(huì)與瀝青質(zhì)發(fā)生直接物理或化學(xué)相互作用而直接影響瀝青質(zhì)的沉積,因此模擬時(shí)假設(shè)水對(duì)瀝青質(zhì)具有化學(xué)惰性。

    圖11 注水開采條件下預(yù)測(cè)期末多孔介質(zhì)中吸附瀝青質(zhì)摩爾分?jǐn)?shù)的分布(P50)

    Yaseen等[11]研究表明,當(dāng)原油與乳化水相互作用時(shí),原油中瀝青質(zhì)的溶解度降低,瀝青質(zhì)開始絮凝,這種情況下可能形成強(qiáng)乳液。該研究中考慮的是在油層注水的情況,油層注水的主要目標(biāo)是將油藏壓力保持在初始水平,并沿注入井到生產(chǎn)井的方向驅(qū)油。而本文考慮的是在水層注水,注水只是為了保持油藏壓力,與油層注水相比,水層注水使得油水其他作用更小。此外,本文通過(guò)水質(zhì)控制(礦化度控制、添加表面活性劑、注入前預(yù)熱)消除了對(duì)瀝青質(zhì)沉積的不良影響。本文還做了以下假設(shè):水與油沒(méi)有任何物理或化學(xué)相互作用(形成乳液、使油中瀝青質(zhì)溶解度降低等);巖石的潤(rùn)濕性沒(méi)有因?yàn)r青質(zhì)沉積發(fā)生變化;原油相對(duì)滲透率的變化僅由局部巖石絕對(duì)滲透率的變化引起,原油黏度與瀝青質(zhì)懸浮固體顆粒的濃度無(wú)關(guān)。瀝青質(zhì)沉積對(duì)水驅(qū)機(jī)理的影響較復(fù)雜,需要單獨(dú)研究。Kamath等[12]指出,瀝青質(zhì)沉積提高了水驅(qū)油效率,這是因?yàn)闉r青質(zhì)沉積導(dǎo)致多孔介質(zhì)表面疏水化以及原油相對(duì)滲透率增加,與此同時(shí)油藏的絕對(duì)滲透率降低。Attar等[13]表示,與水層注水相比,油層注水累計(jì)產(chǎn)油量較低,這與油層注水條件下瀝青質(zhì)沉積更嚴(yán)重有關(guān)。

    4 結(jié)論

    通過(guò)對(duì)組分油藏模擬器中瀝青質(zhì)選項(xiàng)參數(shù)的敏感性分析,確定了對(duì)多變量歷史擬合目標(biāo)函數(shù)影響最大的輸入?yún)?shù):決定溶解瀝青質(zhì)摩爾分?jǐn)?shù)的系數(shù)、瀝青質(zhì)解離率、瀝青質(zhì)吸附系數(shù)和油層中原油流動(dòng)臨界速率。

    針對(duì)激活和未激活瀝青質(zhì)選項(xiàng)這兩種情況,預(yù)測(cè)了衰竭式開采和注水開采條件下的生產(chǎn)剖面。在激活瀝青質(zhì)選項(xiàng)的情況下,使用不同的瀝青質(zhì)參數(shù)計(jì)算了500個(gè)不確定性評(píng)估實(shí)例。結(jié)果表明,衰竭式開采條件下,考慮瀝青質(zhì)沉積情況下累計(jì)產(chǎn)油量的P90,P50,P10值分別比不考慮瀝青質(zhì)沉積情況下的累計(jì)產(chǎn)油量低4.0%,8.1%,12.3%。實(shí)施注水保壓開采后,考慮瀝青質(zhì)沉積情況下累計(jì)產(chǎn)油量的P90,P50,P10值分別比不考慮瀝青質(zhì)沉積情況下的累計(jì)產(chǎn)油量低4.8%,9.1%,13.3%。

    在衰竭式開采條件下,油藏壓力顯著下降并形成壓降漏斗,導(dǎo)致生產(chǎn)井井底瀝青質(zhì)沉積,產(chǎn)能降低。水層注水可以使油藏壓力長(zhǎng)期保持在瀝青質(zhì)初始析出壓力以上,延緩多孔介質(zhì)中瀝青質(zhì)沉積物的形成。與油層注水相比,水層注水使得油水其他作用(如原油與乳化水相互作用使得瀝青質(zhì)的溶解度降低,加速瀝青質(zhì)的析出)相對(duì)更小,可抑制瀝青質(zhì)沉積,因而累計(jì)產(chǎn)油量更高。

    符號(hào)注釋:

    F——多變量歷史擬合目標(biāo)函數(shù);i——響應(yīng)參數(shù)序號(hào);j——實(shí)驗(yàn)次數(shù)序號(hào);N——實(shí)驗(yàn)次數(shù);W——權(quán)重因子;y——響應(yīng)參數(shù)的模擬值;z——響應(yīng)參數(shù)的實(shí)際值;σ——根據(jù)可用歷史數(shù)據(jù)確定的響應(yīng)參數(shù)的標(biāo)準(zhǔn)偏差。

    猜你喜歡
    油藏組分原油
    組分分發(fā)管理系統(tǒng)在天然氣計(jì)量的應(yīng)用
    一種難溶難熔未知組分板材的定性分析
    低滲油藏注采結(jié)構(gòu)調(diào)整與挖潛技術(shù)
    云南化工(2020年11期)2021-01-14 00:51:02
    基于模糊數(shù)學(xué)的油藏干層識(shí)別研究
    黑順片不同組分對(duì)正常小鼠的急性毒性
    中成藥(2018年8期)2018-08-29 01:28:26
    金雀花中黃酮苷類組分鑒定及2種成分測(cè)定
    中成藥(2018年2期)2018-05-09 07:20:09
    放開我國(guó)原油進(jìn)口權(quán)的思考
    初探原油二元期權(quán)
    能源(2016年2期)2016-12-01 05:10:43
    “1噸原油等于7.33桶”的來(lái)龍去脈
    注CO2混相驅(qū)候選油藏篩選評(píng)價(jià)新方法
    亚洲国产精品成人综合色| 成人漫画全彩无遮挡| 人人妻人人看人人澡| 性欧美人与动物交配| 欧美性猛交╳xxx乱大交人| 久久99热这里只有精品18| 91午夜精品亚洲一区二区三区| 国产午夜精品论理片| 日本三级黄在线观看| 免费观看人在逋| 日日撸夜夜添| 久久久久久九九精品二区国产| 大又大粗又爽又黄少妇毛片口| 极品教师在线视频| 国产欧美日韩一区二区精品| 成人性生交大片免费视频hd| 欧美xxxx黑人xx丫x性爽| 亚洲五月天丁香| 毛片一级片免费看久久久久| 日日摸夜夜添夜夜爱| 国产精品女同一区二区软件| 国产精品日韩av在线免费观看| 97人妻精品一区二区三区麻豆| 干丝袜人妻中文字幕| 我要看日韩黄色一级片| 热99在线观看视频| 悠悠久久av| 欧美激情久久久久久爽电影| 99久久精品一区二区三区| 校园春色视频在线观看| 最近在线观看免费完整版| 日韩欧美免费精品| 亚洲综合色惰| 在线观看av片永久免费下载| 亚洲av中文字字幕乱码综合| 中文字幕免费在线视频6| 国产麻豆成人av免费视频| 亚洲av一区综合| 久久人人爽人人爽人人片va| 美女大奶头视频| 99久久久亚洲精品蜜臀av| 欧美+日韩+精品| 国产麻豆成人av免费视频| 亚洲国产精品成人综合色| 久久久久久久午夜电影| 亚洲成人中文字幕在线播放| 91久久精品电影网| 婷婷亚洲欧美| 国产 一区 欧美 日韩| 在线免费观看不下载黄p国产| 观看美女的网站| АⅤ资源中文在线天堂| 女同久久另类99精品国产91| 可以在线观看毛片的网站| 亚洲成a人片在线一区二区| 激情 狠狠 欧美| 99久久中文字幕三级久久日本| 女同久久另类99精品国产91| 欧美潮喷喷水| 欧美另类亚洲清纯唯美| 一本一本综合久久| 偷拍熟女少妇极品色| 国产蜜桃级精品一区二区三区| h日本视频在线播放| 亚洲在线观看片| 久久精品国产99精品国产亚洲性色| 久久综合国产亚洲精品| 男女做爰动态图高潮gif福利片| 国产熟女欧美一区二区| 亚洲欧美成人精品一区二区| 99久国产av精品| 嫩草影院入口| 色哟哟·www| 色噜噜av男人的天堂激情| 在线观看66精品国产| 亚洲色图av天堂| 成人一区二区视频在线观看| 欧美性猛交黑人性爽| 国产男靠女视频免费网站| 在线观看66精品国产| 黄色欧美视频在线观看| 色视频www国产| 成人午夜高清在线视频| 日本五十路高清| 国产伦精品一区二区三区四那| 欧美xxxx黑人xx丫x性爽| 国产 一区精品| 欧美一级a爱片免费观看看| 狂野欧美白嫩少妇大欣赏| av在线播放精品| 国产蜜桃级精品一区二区三区| 22中文网久久字幕| 精品久久久噜噜| 性欧美人与动物交配| 免费观看的影片在线观看| 欧美性感艳星| 99热只有精品国产| av天堂在线播放| 12—13女人毛片做爰片一| 又黄又爽又免费观看的视频| 日韩大尺度精品在线看网址| 一级毛片我不卡| 天堂√8在线中文| 午夜a级毛片| 亚洲电影在线观看av| 国产精品不卡视频一区二区| 免费黄网站久久成人精品| 99热全是精品| 午夜免费激情av| 免费大片18禁| 亚洲18禁久久av| 亚洲精品成人久久久久久| av天堂中文字幕网| 成人一区二区视频在线观看| 日韩欧美国产在线观看| 啦啦啦啦在线视频资源| 久久久国产成人免费| 国产成人a∨麻豆精品| 国产精品三级大全| 日韩欧美一区二区三区在线观看| 国内精品美女久久久久久| 欧美人与善性xxx| 天堂av国产一区二区熟女人妻| 美女高潮的动态| 日日干狠狠操夜夜爽| 麻豆国产97在线/欧美| 菩萨蛮人人尽说江南好唐韦庄 | 午夜精品国产一区二区电影 | 日韩中字成人| 色哟哟哟哟哟哟| 小说图片视频综合网站| 看免费成人av毛片| 波多野结衣巨乳人妻| 国产又黄又爽又无遮挡在线| 又爽又黄无遮挡网站| 日韩亚洲欧美综合| 我要搜黄色片| a级毛片a级免费在线| 噜噜噜噜噜久久久久久91| 偷拍熟女少妇极品色| 九九热线精品视视频播放| a级毛片免费高清观看在线播放| 亚洲av免费在线观看| 91麻豆精品激情在线观看国产| 亚洲美女黄片视频| 亚洲乱码一区二区免费版| 两个人的视频大全免费| av福利片在线观看| 波野结衣二区三区在线| 国产久久久一区二区三区| 久久99热6这里只有精品| 国产中年淑女户外野战色| 在线免费十八禁| 亚洲四区av| 久久久久性生活片| 亚洲人与动物交配视频| 在线免费十八禁| 麻豆成人午夜福利视频| 午夜精品在线福利| 男插女下体视频免费在线播放| 精品欧美国产一区二区三| 狠狠狠狠99中文字幕| 丰满乱子伦码专区| 女生性感内裤真人,穿戴方法视频| 亚洲18禁久久av| 精品乱码久久久久久99久播| 男女做爰动态图高潮gif福利片| avwww免费| 久久精品久久久久久噜噜老黄 | 春色校园在线视频观看| 欧美在线一区亚洲| 成人特级av手机在线观看| 草草在线视频免费看| 精品国产三级普通话版| 亚洲自拍偷在线| 国产真实伦视频高清在线观看| 国产高清激情床上av| 久久久久久伊人网av| 99热网站在线观看| 国产精品,欧美在线| 搡女人真爽免费视频火全软件 | 中文亚洲av片在线观看爽| 国产视频一区二区在线看| 两个人的视频大全免费| 此物有八面人人有两片| 男人舔奶头视频| 久久国产乱子免费精品| 久久久久久九九精品二区国产| av在线蜜桃| 不卡一级毛片| 日本免费一区二区三区高清不卡| 国产av在哪里看| 大型黄色视频在线免费观看| 久久久久久九九精品二区国产| 日本-黄色视频高清免费观看| 欧美日韩综合久久久久久| av福利片在线观看| 国产熟女欧美一区二区| 午夜福利18| 国产精品无大码| 神马国产精品三级电影在线观看| 少妇人妻一区二区三区视频| 美女大奶头视频| 国产成人91sexporn| 亚洲综合色惰| 在线免费观看的www视频| 深夜精品福利| 五月玫瑰六月丁香| 有码 亚洲区| 日韩欧美免费精品| 男女视频在线观看网站免费| 成人特级av手机在线观看| 国产成人福利小说| 波多野结衣高清作品| 天堂动漫精品| 变态另类丝袜制服| 麻豆国产97在线/欧美| 长腿黑丝高跟| 男女边吃奶边做爰视频| 91久久精品国产一区二区成人| 女人十人毛片免费观看3o分钟| 国产黄色小视频在线观看| 亚洲av熟女| 久久国内精品自在自线图片| av女优亚洲男人天堂| 精品无人区乱码1区二区| 精品日产1卡2卡| av黄色大香蕉| 又黄又爽又刺激的免费视频.| 亚洲图色成人| 国产精品伦人一区二区| 啦啦啦啦在线视频资源| 在线免费观看不下载黄p国产| 国产v大片淫在线免费观看| 亚洲国产精品国产精品| 久久精品国产亚洲av天美| 亚洲在线自拍视频| 国产在线男女| 久久久久国产网址| 国产视频一区二区在线看| 日本-黄色视频高清免费观看| 国产爱豆传媒在线观看| 国产精品福利在线免费观看| 日韩精品有码人妻一区| 亚洲美女视频黄频| 国产一区亚洲一区在线观看| a级毛片免费高清观看在线播放| 欧美性猛交╳xxx乱大交人| 国产高清有码在线观看视频| 99国产精品一区二区蜜桃av| 欧美xxxx黑人xx丫x性爽| 亚洲欧美日韩高清专用| 国产69精品久久久久777片| 国产精品永久免费网站| 免费一级毛片在线播放高清视频| 国产免费男女视频| 精品不卡国产一区二区三区| 亚洲激情五月婷婷啪啪| 久久人人精品亚洲av| 国内揄拍国产精品人妻在线| 亚洲综合色惰| 国产一区二区亚洲精品在线观看| 搡老妇女老女人老熟妇| 国模一区二区三区四区视频| 国产淫片久久久久久久久| 美女高潮的动态| 亚洲图色成人| 午夜视频国产福利| 国产白丝娇喘喷水9色精品| 精品福利观看| 久久精品人妻少妇| 高清毛片免费观看视频网站| 日韩av不卡免费在线播放| 天堂动漫精品| 一卡2卡三卡四卡精品乱码亚洲| 国产成人福利小说| 一夜夜www| 成人午夜高清在线视频| 国产真实乱freesex| 99热全是精品| av黄色大香蕉| 国产在线男女| 日本熟妇午夜| 午夜福利在线观看免费完整高清在 | 超碰av人人做人人爽久久| 91久久精品国产一区二区成人| 午夜精品在线福利| 免费看日本二区| 级片在线观看| 国产精品亚洲美女久久久| 我的老师免费观看完整版| 村上凉子中文字幕在线| 亚洲国产欧洲综合997久久,| 日韩成人av中文字幕在线观看 | 精品一区二区三区av网在线观看| 亚洲人成网站在线观看播放| 国语自产精品视频在线第100页| 国产精品综合久久久久久久免费| 国内揄拍国产精品人妻在线| 一级毛片电影观看 | 国国产精品蜜臀av免费| 舔av片在线| 天天躁日日操中文字幕| 久久精品国产亚洲网站| 国产精品一区二区免费欧美| 91狼人影院| 99久久久亚洲精品蜜臀av| 美女大奶头视频| 女生性感内裤真人,穿戴方法视频| 男插女下体视频免费在线播放| 热99re8久久精品国产| 色尼玛亚洲综合影院| 日韩欧美在线乱码| 夜夜看夜夜爽夜夜摸| 久久精品人妻少妇| 亚洲自偷自拍三级| 亚洲七黄色美女视频| 中文字幕精品亚洲无线码一区| 国产成人a区在线观看| 91狼人影院| 日韩av不卡免费在线播放| 美女黄网站色视频| 俄罗斯特黄特色一大片| 麻豆国产av国片精品| 精品人妻一区二区三区麻豆 | 国产精品电影一区二区三区| h日本视频在线播放| 亚洲精品日韩av片在线观看| 天堂√8在线中文| 欧美xxxx黑人xx丫x性爽| 色视频www国产| 成人美女网站在线观看视频| 亚洲国产精品合色在线| 日日啪夜夜撸| 深夜精品福利| 天天躁日日操中文字幕| 久久久国产成人精品二区| 国产在线男女| 国产精品一区二区性色av| 国产不卡一卡二| 国产精品日韩av在线免费观看| 亚洲精品国产av成人精品 | 国产白丝娇喘喷水9色精品| 成人毛片a级毛片在线播放| 欧美日韩精品成人综合77777| 国产一区二区在线av高清观看| 午夜a级毛片| 午夜亚洲福利在线播放| 99在线视频只有这里精品首页| 淫妇啪啪啪对白视频| 欧美区成人在线视频| 黄片wwwwww| 国产一区二区三区在线臀色熟女| 99在线视频只有这里精品首页| 日本一本二区三区精品| 中文字幕av在线有码专区| 免费观看在线日韩| 丰满的人妻完整版| 亚洲av第一区精品v没综合| 韩国av在线不卡| 免费一级毛片在线播放高清视频| 你懂的网址亚洲精品在线观看 | 99久久精品热视频| 久久久久久国产a免费观看| 免费观看的影片在线观看| 内地一区二区视频在线| 久久久久久久久中文| 免费看日本二区| 日日摸夜夜添夜夜添av毛片| 国产极品精品免费视频能看的| 国产蜜桃级精品一区二区三区| 国产精品国产三级国产av玫瑰| 蜜桃久久精品国产亚洲av| 成人综合一区亚洲| 少妇高潮的动态图| 婷婷精品国产亚洲av在线| 久久久久久国产a免费观看| 午夜精品在线福利| 久久久久精品国产欧美久久久| 国产一级毛片七仙女欲春2| 91av网一区二区| 久久久久久久午夜电影| 久99久视频精品免费| 亚洲在线自拍视频| 熟女电影av网| 99热网站在线观看| 亚洲五月天丁香| 国产又黄又爽又无遮挡在线| 男人舔奶头视频| 亚洲高清免费不卡视频| 亚洲中文日韩欧美视频| 国产不卡一卡二| 国产片特级美女逼逼视频| 亚洲性久久影院| 一进一出好大好爽视频| 美女 人体艺术 gogo| 在线a可以看的网站| 丰满人妻一区二区三区视频av| 最近2019中文字幕mv第一页| 又黄又爽又免费观看的视频| 久久久久久久久久黄片| 免费看av在线观看网站| 亚洲欧美日韩东京热| 国产成年人精品一区二区| 欧美3d第一页| 高清午夜精品一区二区三区 | av国产免费在线观看| 狠狠狠狠99中文字幕| 欧美成人精品欧美一级黄| 看免费成人av毛片| 色噜噜av男人的天堂激情| 亚洲国产欧美人成| 99热全是精品| 真实男女啪啪啪动态图| 亚洲av五月六月丁香网| 3wmmmm亚洲av在线观看| 岛国在线免费视频观看| 亚洲人成网站在线播| 婷婷精品国产亚洲av| 婷婷六月久久综合丁香| 亚洲精品久久国产高清桃花| 尤物成人国产欧美一区二区三区| 在线a可以看的网站| 亚洲av美国av| 国产真实乱freesex| 女人被狂操c到高潮| 成人综合一区亚洲| 三级国产精品欧美在线观看| 黑人高潮一二区| 精品免费久久久久久久清纯| 可以在线观看的亚洲视频| 啦啦啦观看免费观看视频高清| 99久久久亚洲精品蜜臀av| 国产精品免费一区二区三区在线| 插阴视频在线观看视频| 国产aⅴ精品一区二区三区波| 亚洲人成网站在线观看播放| 久久久久国产网址| 国产一区亚洲一区在线观看| 亚洲五月天丁香| 亚洲国产精品国产精品| 99热精品在线国产| 国产成人精品久久久久久| 床上黄色一级片| 午夜福利视频1000在线观看| 干丝袜人妻中文字幕| 成人精品一区二区免费| 国产男靠女视频免费网站| 人人妻人人澡欧美一区二区| 女生性感内裤真人,穿戴方法视频| 最近在线观看免费完整版| 在线国产一区二区在线| 国产黄片美女视频| or卡值多少钱| 国产精品伦人一区二区| 国产白丝娇喘喷水9色精品| 我的女老师完整版在线观看| 久久九九热精品免费| 午夜福利成人在线免费观看| 中文字幕久久专区| 亚洲专区国产一区二区| 亚洲av美国av| 国产精品美女特级片免费视频播放器| 国产伦一二天堂av在线观看| 人人妻人人澡欧美一区二区| 一级毛片电影观看 | 欧美精品国产亚洲| 美女内射精品一级片tv| 成年av动漫网址| 国产熟女欧美一区二区| 色噜噜av男人的天堂激情| 美女大奶头视频| 老熟妇乱子伦视频在线观看| 搡女人真爽免费视频火全软件 | 91久久精品国产一区二区三区| 国产成人影院久久av| 国产在线精品亚洲第一网站| 国产女主播在线喷水免费视频网站 | 在线国产一区二区在线| 日日摸夜夜添夜夜添小说| 一进一出抽搐gif免费好疼| 国产精品久久久久久精品电影| 国产精品一区www在线观看| 赤兔流量卡办理| 日韩精品有码人妻一区| 一夜夜www| 亚洲av不卡在线观看| 中文亚洲av片在线观看爽| 看非洲黑人一级黄片| 欧美成人一区二区免费高清观看| 青春草视频在线免费观看| 国产一区二区三区av在线 | 中文字幕免费在线视频6| 亚洲成人av在线免费| 成人特级av手机在线观看| 男人狂女人下面高潮的视频| 亚洲成人av在线免费| 久久精品夜夜夜夜夜久久蜜豆| 美女xxoo啪啪120秒动态图| 最新在线观看一区二区三区| 国产老妇女一区| 中国国产av一级| 国产精品久久电影中文字幕| 亚洲成人久久性| av视频在线观看入口| 大又大粗又爽又黄少妇毛片口| 黄片wwwwww| 男人的好看免费观看在线视频| 亚洲人成网站在线观看播放| 看非洲黑人一级黄片| 久久国内精品自在自线图片| 国产伦精品一区二区三区四那| 色噜噜av男人的天堂激情| 国产精品一区www在线观看| 欧美绝顶高潮抽搐喷水| 黑人高潮一二区| 成年免费大片在线观看| 国产黄片美女视频| 久久久久久国产a免费观看| 精品午夜福利视频在线观看一区| 在线观看66精品国产| 最新中文字幕久久久久| 亚洲国产色片| 黄色配什么色好看| 波多野结衣高清作品| 国产v大片淫在线免费观看| 久久精品91蜜桃| 免费看美女性在线毛片视频| 直男gayav资源| 成人亚洲欧美一区二区av| 俄罗斯特黄特色一大片| 精品久久久久久久人妻蜜臀av| 国产一区二区在线av高清观看| 国产精品亚洲美女久久久| 免费黄网站久久成人精品| 97碰自拍视频| 乱系列少妇在线播放| av在线蜜桃| 别揉我奶头~嗯~啊~动态视频| 精品一区二区三区人妻视频| 麻豆成人午夜福利视频| 亚洲欧美精品自产自拍| 日韩亚洲欧美综合| 亚洲一区二区三区色噜噜| 给我免费播放毛片高清在线观看| 免费观看在线日韩| 男人和女人高潮做爰伦理| 欧美区成人在线视频| 日韩高清综合在线| 高清午夜精品一区二区三区 | 男女边吃奶边做爰视频| 人人妻人人看人人澡| 噜噜噜噜噜久久久久久91| 欧美日本视频| 国产精品精品国产色婷婷| 91在线观看av| 美女xxoo啪啪120秒动态图| 欧美人与善性xxx| 亚洲av第一区精品v没综合| 内射极品少妇av片p| 男女之事视频高清在线观看| 国产精品乱码一区二三区的特点| 欧美一区二区亚洲| 久久精品国产亚洲av涩爱 | 精品免费久久久久久久清纯| 久久久久国产网址| 欧美成人一区二区免费高清观看| 亚洲国产高清在线一区二区三| 午夜福利视频1000在线观看| 国产精品一区二区性色av| 麻豆国产av国片精品| 12—13女人毛片做爰片一| 国产欧美日韩一区二区精品| 一进一出抽搐gif免费好疼| 亚洲成人久久爱视频| 日韩人妻高清精品专区| 日本黄色片子视频| 久久人人精品亚洲av| 中出人妻视频一区二区| 日韩中字成人| 精品久久久久久久久久久久久| 日韩欧美三级三区| 亚洲三级黄色毛片| 国产一区二区激情短视频| 中文资源天堂在线| 亚洲av熟女| 在线观看66精品国产| 好男人在线观看高清免费视频| 悠悠久久av| 麻豆乱淫一区二区| 深夜a级毛片| 一卡2卡三卡四卡精品乱码亚洲| www.色视频.com| 青春草视频在线免费观看| 国产一区二区亚洲精品在线观看| 亚洲国产日韩欧美精品在线观看| 少妇裸体淫交视频免费看高清| 国产精品一区二区三区四区免费观看 | a级毛色黄片| 日本色播在线视频| 国产av一区在线观看免费| 国产精品一区二区三区四区久久| 国产av麻豆久久久久久久| 欧美性感艳星| 国产精品一区二区三区四区免费观看 | 一级av片app| 亚洲美女视频黄频| 亚洲在线观看片| 亚洲无线在线观看| 内射极品少妇av片p| 免费人成视频x8x8入口观看| 插逼视频在线观看| 在线免费十八禁| 国产一区二区三区av在线 | 国产又黄又爽又无遮挡在线|