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

    納米流體驅(qū)提高原油采收率的三維孔隙尺度模擬

    2021-11-03 11:54:08MINAKOVAndreyViktorovichGUZEIDmitriyViktorovichPRYAZHNIKOVMaximIvanovichFILIMONOVSergeyAnatolyevichVORONENKOVAYuliaOlegovna
    石油勘探與開發(fā) 2021年4期
    關(guān)鍵詞:驅(qū)油采收率巖心

    MINAKOV Andrey Viktorovich,GUZEI Dmitriy Viktorovich,PRYAZHNIKOV Maxim Ivanovich,FILIMONOV Sergey Anatol'yevich,VORONENKOVA Yulia Olegovna

    (1.西伯利亞聯(lián)邦大學(xué),克拉斯諾亞爾斯克 660041,俄羅斯;2.庫塔捷拉澤熱物理研究所,俄羅斯科學(xué)院西伯利亞分院,新西伯利亞 630090,俄羅斯)

    0 引言

    通常采用水或含表面活性劑或聚合物添加劑的溶液驅(qū)替原油的方式開發(fā)油藏。研究發(fā)現(xiàn)納米流體應(yīng)用于水驅(qū)可以提高采收率[1-3]。一般認(rèn)為納米流體是添加了納米顆粒的懸浮液(通常為水基懸浮液),納米顆粒的尺寸、濃度和類型不同使得納米顆粒懸浮液具有不同的物理性質(zhì)[1-5]。在驅(qū)油過程中,最重要的影響因素包括驅(qū)替液黏度、接觸角和界面張力,均可以通過納米顆粒進(jìn)行控制。近年來學(xué)者通過室內(nèi)實驗[6-8]和數(shù)值模型[9-25]進(jìn)行了大量關(guān)于水驅(qū)過程中添加各種納米顆粒提高原油采收率的研究。

    模擬納米流體驅(qū)油的數(shù)值模型可分為兩類:將巖石看作均質(zhì)多孔介質(zhì)[9-16]并基于達(dá)西方程求解的模型,以及基于Navier-Stokes方程描述孔隙空間流體流動的直接數(shù)值模型(孔隙尺度模型)[17-25]。Ju等[9-10]首次采用基于達(dá)西定律的一維模型模擬納米流體在多孔介質(zhì)中的流動,但未考慮納米顆粒引起的潤濕性變化;基于此,El-Amin等[11-12]建立了各向異性多孔介質(zhì)中考慮毛管壓力作用和布朗運(yùn)動的二維模型,同樣忽略了潤濕效應(yīng);Feng等[13]基于El-Amin等的模型,研究了納米流體的注入時間、注入量等參數(shù)對采油效率的影響,但未考慮毛管壓力和潤濕性的變化;Sepehri等[14]通過達(dá)西模型研究了納米顆粒引起的巖石潤濕性和地層穩(wěn)定性的變化;Yu等[15]模擬了納米流體在白云巖巖心中的流動,發(fā)現(xiàn)水的礦化度對納米顆粒的運(yùn)移有顯著影響;Rahmani等[16]介紹了一種用于示蹤劑診斷的超順磁性納米顆粒在地層中運(yùn)移的模型。

    在孔隙尺度下對納米流體驅(qū)油進(jìn)行建模的研究很少。該方式采用網(wǎng)格描述多孔介質(zhì)結(jié)構(gòu),需要大量的計算資源,因此主要采用二維模型進(jìn)行相關(guān)研究?;煜嗄P蚚17-20]以及計算相界面的非混相流體模型[21-25]都可以用來描述孔隙尺度模型的多相流,且不需要介質(zhì)孔隙度、相滲透率等經(jīng)驗常數(shù)。求解混相模型中整個計算體的質(zhì)量守恒方程、動量守恒方程和能量守恒方程,基于已知的納米顆粒濃度,根據(jù)混相規(guī)則可以確定每個控制體中液體的物理性質(zhì)。Gharibshahi等[17]通過混相模型研究了二維微觀模型中孔隙的形狀和分布對驅(qū)油過程的影響,基于此探討了納米顆粒的類型、濃度和尺寸以及注入流體的溫度對驅(qū)油效果的影響[18],但仍舊沒有考慮潤濕性的影響。

    與混相模型相反,基于流體體積法(VOF)的非混相流體模型直接計算非混相液體界面,從而準(zhǔn)確得到兩種液體之間的表面張力和接觸角,而這兩種性質(zhì)在注納米流體的過程中對提高原油采收率起著關(guān)鍵作用[6-8]。VOF在提高原油采收率相關(guān)的研究中(如注表面活性劑、聚合物、二氧化碳、蒸汽)應(yīng)用非常廣泛,然而,關(guān)于采用VOF對納米流體驅(qū)油建模的研究很少,一般使用二維模型[21]或單個孔道的簡化模型進(jìn)行研究,Zhao和Wen[22]基于此研究了納米流體從兩個通過狹窄通道相連的圓形孔隙中驅(qū)油的過程。

    本研究基于巖心三維微觀模型對納米流體驅(qū)油過程進(jìn)行直接數(shù)值模擬,該模型在結(jié)構(gòu)上與真實巖心非常接近。在三維微觀模型中采用VOF研究納米顆粒的質(zhì)量分?jǐn)?shù)和粒徑、驅(qū)替流體流速、油黏度和巖心滲透率等因素對納米流體驅(qū)油效率的影響。

    1 實驗研究

    為了進(jìn)行準(zhǔn)確的數(shù)值模擬,對納米流體-油-巖石體系中的界面張力和潤濕性進(jìn)行了系統(tǒng)的實驗研究獲取相關(guān)數(shù)據(jù);關(guān)于納米顆粒對潤濕性影響的實驗結(jié)果詳見文獻(xiàn)[26],實驗中使用的輕質(zhì)油密度為0.831 g/cm3,黏度為7.8 mPa·s。

    將SiO2納米粉末添加到蒸餾水中制備納米流體,然后進(jìn)行充分機(jī)械混合。其中納米顆粒的粒徑為5~50 nm,顆粒質(zhì)量分?jǐn)?shù)為0~1%。將懸浮液進(jìn)行超聲分散處理分解納米顆粒團(tuán)塊,采用TurbiscanLAB分析儀研究發(fā)現(xiàn)納米流體在制備后10 d內(nèi)是穩(wěn)定的。使用聲學(xué)和電聲分析儀 DT1202直接測量懸浮液中的顆粒平均粒徑和Zeta電位。由于納米顆粒易于聚集,液體中的顆粒平均粒徑與原始粉末的粒徑不同,且液體中顆粒的平均粒徑與原始粒徑成正比(見表 1)。大部分懸浮液的Zeta電位絕對值超過20 mV,說明其膠體穩(wěn)定性高。

    表1 納米顆粒的實驗測量數(shù)據(jù)

    用IFT-820-P自動張力儀測量納米流體-油-巖石界面張力和接觸角,該張力儀的工作原理基于懸滴法,即通過測量液滴的幾何參數(shù)來確定界面張力。使用DropImage Advanced軟件處理測量結(jié)果,對4個測量值取平均值,實驗中的數(shù)據(jù)方差不超過5%。用該表面張力測量方法對水和乙二醇進(jìn)行了測試,測試結(jié)果與參考數(shù)據(jù)吻合較好[27-28],證明方法可靠。

    使用 OFITE高溫高壓黏度計測定了納米流體黏度,測量誤差約2%。采用Mettler Toledo精密天平測量樣品質(zhì)量、采用刻度燒杯測量樣品體積,以計算納米流體的密度,密度測量誤差約1%。

    含平均粒徑為5 nm的SiO2顆粒的納米流體測量結(jié)果如圖1和圖2所示。納米顆粒質(zhì)量分?jǐn)?shù)對巖石表面親油特性影響顯著(見圖 1),油接觸角隨著納米顆粒質(zhì)量分?jǐn)?shù)的增加,從73°增至153°,因此,向液體中添加納米顆粒會削弱巖石的親油性;當(dāng)納米顆粒質(zhì)量分?jǐn)?shù)為1%時,巖心不親油。同時,隨著納米顆粒質(zhì)量分?jǐn)?shù)的增加,界面張力變化很?。ㄒ妶D2),最大降幅約為7%。

    圖1 不同納米顆粒質(zhì)量分?jǐn)?shù)下納米流體中的油滴在巖心表面的照片

    圖2 納米流體-油的界面張力、接觸角與顆粒質(zhì)量分?jǐn)?shù)的關(guān)系

    通過實驗研究納米顆粒粒徑對納米流體中巖石親油性的影響。向水中添加粒徑為5~50 nm的SiO2納米顆粒,納米顆粒質(zhì)量分?jǐn)?shù)為 0.5%。研究了 25 ℃下納米顆粒粒徑對納米流體-油-巖石接觸角和界面張力的影響(見表2)??梢钥闯觯S著納米顆粒粒徑的減小,油-納米流體-巖石邊界處的接觸角增大,因此,隨著納米顆粒粒徑的減小,納米流體的洗油能力提高;油和納米流體的界面張力隨著納米顆粒粒徑的增加而降低約25%。

    表2 含不同粒徑SiO2納米顆粒的懸浮液測量結(jié)果

    2 數(shù)學(xué)模型與數(shù)值方法

    采用三維多孔介質(zhì)微觀模型對納米流體驅(qū)油過程進(jìn)行了系統(tǒng)的數(shù)值研究。在僅考慮層流的條件下,使用流體體積法(VOF)[30]模擬多孔介質(zhì)中的兩相流。該模型通過求解動量守恒方程,利用各相的體積分?jǐn)?shù)模擬非混相的液-液多組分流動。引入驅(qū)替液體積分?jǐn)?shù)α和油體積分?jǐn)?shù) β,α+β=1。首先,用巖石的非潤濕相(水)驅(qū)替潤濕相(油)。Navier-Stokes方程中混合物的密度和黏度用下式確定:

    質(zhì)量守恒方程標(biāo)準(zhǔn)形式如下:

    v是混合物速度矢量,通過求解動量方程來定義:

    通過求解輸運(yùn)方程來模擬驅(qū)油過程中的界面運(yùn)動:

    Brackbill等[31]用連續(xù)表面力模型模擬毛管壓力,該方法通過在動量方程中添加源項來模擬表面張力:

    在計算域壁面上,使用以下表達(dá)式確定法向量:

    在計算中,納米流體的黏度、密度、界面張力和接觸角的實驗測量值取決于顆粒質(zhì)量分?jǐn)?shù)和粒徑。本文使用的數(shù)值技術(shù)及其測試結(jié)果在文獻(xiàn)[32-33]中有詳細(xì)的描述。使用有限體積法(FVM)求解非線性微分方程組(3)式—(5)式;采用壓力耦合關(guān)聯(lián)方程的半隱式方法(SIMPLEC)實現(xiàn)了速度場和壓力場的耦合;使用預(yù)壓交錯選項(PRESTO)進(jìn)行連續(xù)平衡計算;使用二階隱式法和二階中心差分法分別估算輸運(yùn)方程的非穩(wěn)態(tài)項和對流項;使用帶有高分辨率界面捕捉(HRIC)方案的全變差遞減法(TVD)求解(5)式。

    3 數(shù)值研究結(jié)果

    3.1 模型描述

    本研究對納米流體驅(qū)油過程開展了數(shù)值模擬,并提出采用三維數(shù)字巖心模型進(jìn)行建模的方法,結(jié)構(gòu)接近真實巖心。采用孔隙度分別為 25%,20%,7%的數(shù)字巖心,滲透率為(3.1~52.0)×10?3μm2。由于最小流動孔徑明顯大于最小納米顆粒粒徑,因此,本研究忽略納米顆粒對巖心滲透率的影響。計算域是尺寸為150 μm的立方體,圖3所示為不同滲透率的巖心模型。

    圖3 不同滲透率的數(shù)字巖心模型

    通過一系列方法學(xué)計算選擇了 350萬個單元的細(xì)化網(wǎng)格進(jìn)行模擬,保證采收率的誤差為零?;谇拔慕o出的黏度、界面張力和接觸角相關(guān)實驗測量數(shù)據(jù)進(jìn)行計算,選擇含粒徑為5 nm的SiO2顆粒的納米流體,顆粒質(zhì)量分?jǐn)?shù)為0~1%?;A(chǔ)油模型的黏度設(shè)定為7.8 mPa·s,密度為 0.831 g/cm3,油水界面張力為 22.5 mN/m。模型的邊界條件為:計算域上邊緣的驅(qū)替液流速為 7.8×10?6~2.3×10?1m/s,計算域下表面設(shè)置為自由出口,計算域的側(cè)壁設(shè)置黏附條件。

    在非穩(wěn)態(tài)下進(jìn)行計算,計算開始時孔隙空間被油飽和。初始含水飽和度為零,模擬研究了計算域內(nèi)的驅(qū)替前緣動態(tài)、油突破時間、剩余油飽和度、原油采收率的瞬時值和最終值以及入口到出口的壓降。

    3.2 納米顆粒質(zhì)量分?jǐn)?shù)對驅(qū)油效率的影響

    為了研究納米流體中顆粒質(zhì)量分?jǐn)?shù)對驅(qū)油效率的影響,設(shè)計了基本計算條件:巖心滲透率為3.1×10?3μm2,油黏度為 7.8 mPa·s,驅(qū)替液流速為 7.05×10?4m/s。圖4是水和顆粒質(zhì)量分?jǐn)?shù)為1%的納米流體的驅(qū)油過程,用驅(qū)替液在多孔介質(zhì)壁面上的體積分?jǐn)?shù)的等值面來表征。納米流體驅(qū)油時,隨著顆粒質(zhì)量分?jǐn)?shù)的增加,巖石表面由疏水性轉(zhuǎn)變?yōu)橛H水性。對孔隙通道內(nèi)局部驅(qū)替前緣的分析表明,在水濕系統(tǒng)中(w=1%),水相沿著多孔介質(zhì)孔壁運(yùn)移,此時水驅(qū)更容易洗油;在油濕系統(tǒng)中(w=0),水相沿著孔隙通道的中心運(yùn)移,孔壁被油浸潤。與水驅(qū)相比,納米流體填充在孔隙通道橫截面中,納米流體的驅(qū)替前緣更接近活塞驅(qū)替且占據(jù)的孔道體積比例顯著增加,在整個巖心界面上的前緣分布更均勻;而水主要沿一個特定的方向流動,大部分油不會被驅(qū)掃(見圖 4)。另外,水突破時間為0.013 s,而納米流體的突破時間為0.027 s。

    圖4 不同驅(qū)替時間水和納米流體在多孔介質(zhì)中的分布(驅(qū)替方向從上到下)

    圖 5為具有不同顆粒質(zhì)量分?jǐn)?shù)的納米流體的采收率,其中 V/Vp是注入流體體積與孔隙體積之比??梢钥闯?,隨著驅(qū)替液中顆粒質(zhì)量分?jǐn)?shù)的增加采收率增大,油停止流出巖心的時間也隨之延長,表明油幾乎完全洗凈。最終采收率隨顆粒質(zhì)量分?jǐn)?shù)的變化如圖6a所示,采收率隨著顆粒質(zhì)量分?jǐn)?shù)的增加而增大。與水相比,顆粒質(zhì)量分?jǐn)?shù)為 0.25%的 SiO2納米流體可以提高采收率約10%;顆粒質(zhì)量分?jǐn)?shù)增加到0.5%時,采收率可提高約19%。當(dāng)質(zhì)量分?jǐn)?shù)從0.5%變?yōu)?.0%時,接觸角和界面張力的變化都不明顯,因此增加顆粒質(zhì)量分?jǐn)?shù)不會進(jìn)一步增加采收率。納米流體均勻的驅(qū)油剖面會顯著延長水的突破時間(見圖6b),即通過納米流體可以從多孔介質(zhì)中洗出更多的油。

    圖5 不同顆粒質(zhì)量分?jǐn)?shù)下SiO2納米流體的采收率

    圖6 原油采收率(a)和水的突破時間(b)與SiO2顆粒質(zhì)量分?jǐn)?shù)的關(guān)系

    圖 7顯示了不同顆粒質(zhì)量分?jǐn)?shù)下納米流體相對體積與計算域入口與出口間壓降的關(guān)系??梢钥闯?,驅(qū)替液滲透到飽和油多孔介質(zhì)中的過程伴隨壓力脈沖,這與油從巖壁中的微分離有關(guān)。在驅(qū)替液完全突破計算域的出口后,壓降減小,壓力脈沖停止,流動呈擬穩(wěn)態(tài)。盡管納米流體的黏度比純水高約 10%,但其壓降通常低于純水,這是由于添加納米顆粒可以降低界面張力、改善巖石的潤濕性,從而降低毛管壓力,但隨著納米顆粒質(zhì)量分?jǐn)?shù)的進(jìn)一步增加,由于界面張力和驅(qū)替液黏度的增加,總壓力增加。

    圖7 不同顆粒質(zhì)量分?jǐn)?shù)下SiO2納米流體在驅(qū)替過程中的壓降

    3.3 納米顆粒粒徑對驅(qū)替效果的影響

    對4種顆粒粒徑分別為5,18,22,50 nm的SiO2納米流體進(jìn)行模擬計算,顆粒質(zhì)量分?jǐn)?shù)為0.5%,納米流體的接觸角、界面張力和黏度的實驗數(shù)據(jù)如表 2所示。隨納米顆粒粒徑的減小,黏度和親水性增加,界面張力增大,接觸角增大幅度最大。增加驅(qū)替液的黏度和接觸角有利于提高采收率,因此納米顆粒粒徑是控制驅(qū)油效率的主要因素。通過減小納米流體中的顆粒粒徑,巖石表面從疏水轉(zhuǎn)變?yōu)橛H水,數(shù)值模擬結(jié)果表明驅(qū)油量增加。SiO2納米顆粒粒徑從50 nm減小到5 nm,采收率提高了約29個百分點(diǎn)。從納米流體中顆粒粒徑與采收率的關(guān)系可以看出(見圖8),使用較小粒徑的納米顆??梢蕴岣咴筒墒章?,但減小納米顆粒粒徑的同時會削弱界面張力對采收率的影響。在本研究中,隨著顆粒粒徑的減小,界面張力增大了約30%,低界面張力更有利于采油,此時降低顆粒尺寸增加接觸角對采收率的影響比增大界面張力的影響更為明顯。一般情況下,存在最優(yōu)納米顆粒粒徑,可以在給定的顆粒濃度下使采收率最大化。

    圖8 納米流體驅(qū)采收率與顆粒粒徑的關(guān)系

    3.4 驅(qū)替液流速對驅(qū)替效果的影響

    通過研究驅(qū)替液流速對驅(qū)替特性的影響[34-37],發(fā)現(xiàn)剩余油飽和度取決于毛管數(shù),即毛管數(shù)對原油采收率有很大影響。本文使用定義式Nco=vμ/σ確定毛管數(shù),先前的研究表明,如果Nco超過臨界值,剩余油飽和度會單調(diào)減小,采收率單調(diào)遞增。在巖心實驗中發(fā)現(xiàn),毛管數(shù)的臨界值與潤濕性有關(guān),為 1×10?5~1×10?3。利用毛管數(shù)針對納米流體驅(qū)油進(jìn)行了一系列模擬計算,結(jié)果表明,驅(qū)替液的流速為 7.8×10?6~2.3×10?1m/s時,對應(yīng)的毛管數(shù)為 3.6×10?7~1.1×10?2。不同驅(qū)替速度下水驅(qū)和顆粒質(zhì)量分?jǐn)?shù)為1%的納米流體驅(qū)注入3倍孔隙體積后巖心中部截面的剩余油分布如圖 9所示。納米流體驅(qū)替時,因為孔隙壁面呈水濕,納米流體沿孔隙通道驅(qū)替更均勻,從而洗出了更多的油,此時剩余油分布于孔隙通道的中心;水驅(qū)替時,剩余油主要分布于大通道壁面的滯留腔中。隨著流速的增加,毛管力的作用減弱,多數(shù)孔隙通道參與流體流動,從而提高了采油效率。

    圖9 不同毛管數(shù)下巖心中部截面剩余油的最終分布(驅(qū)替方向從上到下)

    本研究中低毛管數(shù)條件下,毛管壓力滯留了多孔介質(zhì)中的大部分油,因此該模式下的剩余油飽和度和采收率與驅(qū)替速度無關(guān);當(dāng)毛管數(shù)達(dá)到臨界值時,剩余油飽和度隨驅(qū)替速度的增加而降低,而采收率隨驅(qū)替速度的增加而增加(見圖10),在此模式下黏滯力作用超過毛管壓力。在本實驗驅(qū)替速度范圍內(nèi),納米流體驅(qū)油的采收率均高于水驅(qū)采收率(見圖10a),同時,納米顆粒對采收率的影響隨著驅(qū)替速度的增加而減弱;Nco=3.6×10?7時屬于毛管壓力驅(qū)替模式,納米流體驅(qū)油采收率是水驅(qū)采收率的1.4倍;Nco=1.1×10?2時屬于壓力驅(qū)替模式,納米流體驅(qū)油采收率是水驅(qū)采收率的1.14倍。因此,在毛管壓力驅(qū)替模式下,納米顆粒對驅(qū)油效率的影響更顯著。此外,驅(qū)替液中添加納米顆粒時,毛管數(shù)臨界值從約 1×10?4減小到約 1×10?5,降低了幾乎一個數(shù)量級(見圖 10b),說明注入少量的驅(qū)替液就可以提高采收率。

    圖10 采收率與驅(qū)替速度的關(guān)系(a)及殘余油飽和度與毛管數(shù)的關(guān)系(b)

    3.5 原油黏度的影響

    采用納米流體驅(qū)替含不同黏度原油的巖心并計算其采收率。原油黏度為0.78~780 mPa·s,原油黏度與驅(qū)替液黏度之比為 0.8~800.0,懸浮液中的 SiO2顆粒質(zhì)量分?jǐn)?shù)為1%,驅(qū)替液流速對應(yīng)的毛管數(shù)為3.6×10?7(毛管壓力驅(qū)替模式),巖心滲透率為 3.1×10?3μm2。圖11所示為注入3倍孔隙體積的水和納米流體驅(qū)替不同黏度原油后驅(qū)替流體的分布,可以看出納米流體驅(qū)替前緣擴(kuò)散更寬,而水則形成穿透巖心的分散射流。研究表明,當(dāng)水驅(qū)替高黏度原油時,水在驅(qū)替前緣形成局部突進(jìn)并迅速到達(dá)巖心出口,導(dǎo)致采收率降低;而納米流體的驅(qū)替前緣移動明顯不同,由于納米流體具有親水性,因此主要沿巖石壁面運(yùn)移。納米流體的潤濕性良好,可以洗掉黏度更大的油,可用于提高稠油采收率。

    圖11 不同黏度比下巖心中部縱向剖面中驅(qū)替流體的最終分布(驅(qū)替方向從上到下)

    如圖12所示,采收率隨著原油與驅(qū)替液黏度比的增加而降低,這與前人的研究一致[32-36]。在任何黏度比下,向驅(qū)替液中添加質(zhì)量分?jǐn)?shù)為1%的納米顆粒都可以提高采收率,且隨著黏度比增加提高采收率幅度增加。黏度比為0.8時,納米流體驅(qū)油的采收率是水驅(qū)的1.2倍;而黏度比為800.0時,納米流體驅(qū)的采收率是水驅(qū)的2.2倍。通常當(dāng)油與驅(qū)替液的黏度比超過10時,采收率開始急劇下降;但從圖 12a看出,納米流體驅(qū)油采收率的急劇降低始于黏度比大于 100時,說明納米流體可以在更大的黏度比范圍內(nèi)有效驅(qū)油。從驅(qū)替液突破時間與黏度比的關(guān)系曲線可以看出(見圖12b),突破時間隨著黏度比的增加而減少;納米流體突破時間比水更長,可以從巖心中驅(qū)出更多的油。因此,納米流體比水的波及效率更高,其提高采收率幅度隨油與驅(qū)替介質(zhì)的黏度比的增加而增加。

    圖12 采收率(a)和突破時間(b)與黏度比的關(guān)系

    迄今為止,已經(jīng)進(jìn)行了大量關(guān)于砂巖驅(qū)油的研究,但由于殘余油飽和度與驅(qū)替速度、流體黏度、潤濕性、孔隙度等多種因素有關(guān),目前還沒有統(tǒng)一的毛管數(shù)定義[34-37]。毛管數(shù)的臨界值決定了毛細(xì)管滯留油的能力,本研究使用Abrams[35]提出的毛管數(shù)確定法:

    將歸一化的殘余油飽和度與已知的實驗數(shù)據(jù)進(jìn)行比較(見圖13)??梢钥闯?,無論是毛管數(shù)的臨界值,還是壓力驅(qū)替階段曲線的斜率,計算結(jié)果和實驗結(jié)果都具有相當(dāng)好的一致性,證明了數(shù)值模擬的可靠性。

    圖13 歸一化殘余油飽和度與毛管數(shù)的關(guān)系

    3.6 巖心滲透率的影響

    通過對納米流體驅(qū)替不同滲透率的巖心進(jìn)行數(shù)值模擬研究,考慮低滲和高滲儲集層模型,巖心模型滲透率為(3.1~52.0)×10?3μm2。對水和顆粒質(zhì)量分?jǐn)?shù)為 1%的 SiO2納米流體驅(qū)油進(jìn)行模擬計算,毛管數(shù)為3.6×10?7。巖心模型中納米顆粒運(yùn)動軌跡和驅(qū)替后油的三維分布如圖14所示。

    圖14 納米流體驅(qū)替不同滲透率巖心納米顆粒運(yùn)動軌跡和驅(qū)替后油的最終分布(驅(qū)替方向從上到下)

    儲集層的采收率主要由儲集層性質(zhì)決定,通常采收率隨滲透率的減小而降低(見圖15)。這是因為低滲儲集層模型中的孔隙通道比高滲儲集層中的窄,狹窄孔隙通道越多,被毛細(xì)管滯留的剩余油越多,驅(qū)替效率越低(見圖16)。

    圖15 采收率與巖心滲透率的關(guān)系

    圖16 不同滲透率巖心的中部縱向剖面中驅(qū)替液的最終分布(驅(qū)替方向從上到下)

    納米流體對采收率的影響程度取決于儲集層滲透率,添加納米顆粒的效果隨滲透率的增加而降低。研究表明相對于水驅(qū),低滲儲集層和高滲儲集層中納米流體驅(qū)油的采收率分別提高了40和30個百分點(diǎn)。在較高滲透率的巖心中,驅(qū)替液呈較寬的驅(qū)替前緣展開,占據(jù)了大量的孔隙空間,大孔道的油被水沖洗得相對較好,納米流體對表面潤濕性的影響與驅(qū)替效率沒有太大相關(guān)性。對于低滲儲集層,毛管壓力將大量油滯留在狹窄的孔隙通道中,在這種情況下,添加納米顆粒至驅(qū)替液中時,表面潤濕性的改善起重要作用,納米流體可以洗掉更多的毛細(xì)管滯留油(見圖 16)。因此,結(jié)果表明,納米流體可以顯著提高原油采收率;當(dāng)巖石滲透率降低時,使用納米流體的效果增強(qiáng)。

    4 結(jié)論

    向驅(qū)替液中添加納米顆粒會顯著影響多孔介質(zhì)驅(qū)油效果,采收率隨著顆粒質(zhì)量分?jǐn)?shù)的增加而增加,納米顆粒質(zhì)量分?jǐn)?shù)增加到0.5%時,與水驅(qū)相比采收率可提高約 19%;納米流體驅(qū)油效率取決于納米顆粒的尺寸,采收率隨著納米顆粒粒徑的減小而增加;將納米顆粒添加到驅(qū)替液中會使得毛管數(shù)的臨界值降低一個數(shù)量級,在毛管數(shù)接近臨界值的注入模式下使用納米流體驅(qū)油對高黏低滲儲集層提高采收率最為有效,且提高采收率幅度隨著驅(qū)替速度的增加而減??;原油黏度越大,巖石滲透率越小,納米流體驅(qū)油提高采收率的效果越明顯。

    致謝:實驗研究中的潤濕性測定由俄羅斯科學(xué)基金會(Russian Science Foundation(17-79-20218-P))支持,數(shù)值模擬由國家科學(xué)任務(wù)體系支持(西伯利亞聯(lián)邦大學(xué)(Siberian Federal University,F(xiàn)SRZ-2020-0012)),納米顆粒表征由克拉斯諾亞爾斯克地區(qū)共享研究中心(俄羅斯科學(xué)院西伯利亞分支克拉斯諾亞爾斯克科學(xué)中心)和西伯利亞聯(lián)邦大學(xué)共享研究中心完成,在此一并致謝。

    符號注釋:

    Nco——毛管數(shù),無因次;Nco*——根據(jù) Abrams提出的方法計算的毛管數(shù),無因次;Fs——由毛細(xì)管力引起的體積力矢量,N/m3;k——界面的曲率,m?1;n——相之間界面的法向量,無因次;nw——與壁向量正交的向量,無因次;p——混合物的靜壓力,Pa;t——時間,s;v——混合速度矢量,m/s;V/Vp——注入流體體積與孔隙體積之比,無因次;v——多孔介質(zhì)中的流體流速,m/s;w——納米顆粒質(zhì)量分?jǐn)?shù),%;α,β——驅(qū)替液體積分?jǐn)?shù)和油體積分?jǐn)?shù),無因次;Δp——計算域入口與出口間的壓降,Pa;θ(φ)——壁面處的接觸角,(°);μ——黏度,Pa·s;μo——油的黏度,Pa·s;μr——原油黏度與驅(qū)替流體黏度之比,無因次;μw——驅(qū)替液的黏度,Pa·s;μ0——基液黏度,Pa·s;μ1(φ)——納米流體的黏度,Pa·s;ρ——密度,kg/m3;ρ1(φ)——納米流體的密度,kg/m3;ρo——原油密度,kg/m3;σ(φ)——界面張力,N/m;wτ——與壁向量相切的向量,無因次;φ——納米顆粒體積分?jǐn)?shù),%。

    猜你喜歡
    驅(qū)油采收率巖心
    《油氣地質(zhì)與采收率》征稿簡則
    《油氣地質(zhì)與采收率》征稿簡則
    《油氣地質(zhì)與采收率》第六屆編委會
    《油氣地質(zhì)與采收率》征稿簡則
    Osteotomized folded scapular tip free flap for complex midfacial reconstruction
    注氣驅(qū)油技術(shù)發(fā)展應(yīng)用及海上油田啟示
    一種頁巖巖心資料的保存方法
    化工管理(2017年23期)2017-09-11 14:14:22
    Acellular allogeneic nerve grafting combined with bone marrow mesenchymal stem cell transplantation for the repair of long-segment sciatic nerve defects: biomechanics and validation of mathematical models
    CO2驅(qū)油與埋存對低碳經(jīng)濟(jì)的意義
    長巖心注CO2氣水交替驅(qū)試驗?zāi)M研究
    少妇精品久久久久久久| 一边亲一边摸免费视频| 91精品伊人久久大香线蕉| 少妇熟女欧美另类| 国产午夜精品一二区理论片| 国产黄色免费在线视频| 黑人高潮一二区| 在线观看免费高清a一片| 国产一级毛片在线| 午夜免费观看性视频| 老司机影院毛片| 3wmmmm亚洲av在线观看| 亚洲,一卡二卡三卡| 天堂俺去俺来也www色官网| 久久婷婷青草| 日韩人妻高清精品专区| 亚洲欧洲日产国产| 少妇人妻精品综合一区二区| 欧美高清成人免费视频www| 中文天堂在线官网| 午夜激情久久久久久久| 日韩 亚洲 欧美在线| 成人综合一区亚洲| a级毛片在线看网站| 国产亚洲91精品色在线| av女优亚洲男人天堂| 涩涩av久久男人的天堂| 寂寞人妻少妇视频99o| 亚洲伊人久久精品综合| 我要看黄色一级片免费的| av线在线观看网站| 亚洲欧美一区二区三区黑人 | 最近2019中文字幕mv第一页| 国产精品麻豆人妻色哟哟久久| 熟女av电影| 97精品久久久久久久久久精品| 成人18禁高潮啪啪吃奶动态图 | 国产av一区二区精品久久| 亚洲精品一二三| 亚洲成人手机| 国产亚洲精品久久久com| 2021少妇久久久久久久久久久| 亚洲丝袜综合中文字幕| 一级毛片aaaaaa免费看小| 在线播放无遮挡| 少妇高潮的动态图| 亚洲国产精品一区二区三区在线| 自拍偷自拍亚洲精品老妇| 国产极品粉嫩免费观看在线 | 2022亚洲国产成人精品| 亚洲伊人久久精品综合| 亚洲国产精品一区三区| 少妇 在线观看| 精品久久久噜噜| 免费看光身美女| 91aial.com中文字幕在线观看| 久久影院123| 狂野欧美激情性bbbbbb| 亚洲国产成人一精品久久久| av国产久精品久网站免费入址| 久久久亚洲精品成人影院| 大话2 男鬼变身卡| 亚洲av二区三区四区| av有码第一页| 麻豆成人av视频| 成人亚洲精品一区在线观看| 国产成人精品久久久久久| 亚洲激情五月婷婷啪啪| 日本午夜av视频| 亚洲av.av天堂| 精品久久国产蜜桃| 91在线精品国自产拍蜜月| 一本大道久久a久久精品| 成年人午夜在线观看视频| 久久精品久久久久久久性| 国产乱人偷精品视频| 午夜免费观看性视频| 黄色一级大片看看| 桃花免费在线播放| 亚洲国产精品成人久久小说| 人人妻人人爽人人添夜夜欢视频 | 亚洲精品久久午夜乱码| 日日啪夜夜爽| 中文天堂在线官网| 日本av手机在线免费观看| 极品教师在线视频| 高清午夜精品一区二区三区| 免费av中文字幕在线| 少妇 在线观看| 日韩成人伦理影院| 婷婷色综合大香蕉| 精品国产国语对白av| 精品人妻一区二区三区麻豆| 两个人的视频大全免费| 免费大片18禁| 国产深夜福利视频在线观看| 成年美女黄网站色视频大全免费 | 亚洲国产色片| 国产精品不卡视频一区二区| 国产 一区精品| 91在线精品国自产拍蜜月| 成年av动漫网址| 久久久久久久久久人人人人人人| 欧美精品国产亚洲| 韩国高清视频一区二区三区| 国产白丝娇喘喷水9色精品| 亚洲国产av新网站| av女优亚洲男人天堂| av视频免费观看在线观看| 蜜臀久久99精品久久宅男| 欧美日韩视频高清一区二区三区二| 另类精品久久| 国产av码专区亚洲av| 97超视频在线观看视频| 国产日韩欧美视频二区| 国产精品.久久久| 七月丁香在线播放| 亚洲精品成人av观看孕妇| 99re6热这里在线精品视频| 3wmmmm亚洲av在线观看| 国产精品熟女久久久久浪| 两个人免费观看高清视频 | 中文字幕亚洲精品专区| 九九在线视频观看精品| 丰满乱子伦码专区| 18禁动态无遮挡网站| 久久久久久久久久人人人人人人| 国产精品三级大全| 亚洲av不卡在线观看| 黄色配什么色好看| 亚洲精品色激情综合| 在现免费观看毛片| 欧美人与善性xxx| 日韩中字成人| 亚洲欧美中文字幕日韩二区| 黄色一级大片看看| 啦啦啦在线观看免费高清www| 亚洲人与动物交配视频| av有码第一页| 国国产精品蜜臀av免费| 中文字幕制服av| 亚洲成人av在线免费| 大片电影免费在线观看免费| 建设人人有责人人尽责人人享有的| 国产成人freesex在线| 搡女人真爽免费视频火全软件| 一本大道久久a久久精品| 亚洲丝袜综合中文字幕| 视频区图区小说| 人妻人人澡人人爽人人| 蜜桃在线观看..| 五月伊人婷婷丁香| 亚洲国产精品成人久久小说| 女人久久www免费人成看片| 久久久久视频综合| 免费看日本二区| 99视频精品全部免费 在线| 国产男人的电影天堂91| 久久精品久久久久久久性| 久久久久国产网址| 一边亲一边摸免费视频| 九九爱精品视频在线观看| 成人亚洲欧美一区二区av| 亚洲伊人久久精品综合| 色视频www国产| 国产欧美亚洲国产| 在线观看三级黄色| 亚洲精品国产av蜜桃| 国产免费福利视频在线观看| 多毛熟女@视频| 国产无遮挡羞羞视频在线观看| 五月开心婷婷网| 国产精品国产三级国产av玫瑰| 亚洲人与动物交配视频| 国产av码专区亚洲av| 亚洲自偷自拍三级| 国产一级毛片在线| 日韩人妻高清精品专区| 午夜91福利影院| 久久国产精品男人的天堂亚洲 | 99久久精品一区二区三区| 日本与韩国留学比较| 亚洲国产最新在线播放| 亚洲,欧美,日韩| 中国三级夫妇交换| 久久精品久久久久久久性| 精品国产露脸久久av麻豆| av一本久久久久| 成年av动漫网址| 大话2 男鬼变身卡| 国产成人免费无遮挡视频| 国产色婷婷99| 国内揄拍国产精品人妻在线| videossex国产| 精品99又大又爽又粗少妇毛片| 成人美女网站在线观看视频| 久久亚洲国产成人精品v| 亚洲av日韩在线播放| 成人国产av品久久久| 午夜免费观看性视频| 欧美精品一区二区大全| av又黄又爽大尺度在线免费看| 亚洲欧美清纯卡通| 一个人免费看片子| 久久国产精品大桥未久av | 国产精品麻豆人妻色哟哟久久| 最新的欧美精品一区二区| 王馨瑶露胸无遮挡在线观看| 91精品一卡2卡3卡4卡| 草草在线视频免费看| 亚洲国产精品成人久久小说| 国产男女超爽视频在线观看| 国产综合精华液| 国产色爽女视频免费观看| 国产极品天堂在线| 日韩亚洲欧美综合| 黑丝袜美女国产一区| 寂寞人妻少妇视频99o| 中国国产av一级| 久久久国产欧美日韩av| 丝袜喷水一区| 国产精品福利在线免费观看| av不卡在线播放| 蜜桃在线观看..| 色婷婷av一区二区三区视频| 女的被弄到高潮叫床怎么办| 亚洲精品,欧美精品| 成人无遮挡网站| 国产成人午夜福利电影在线观看| 亚洲精品色激情综合| 香蕉精品网在线| 人妻夜夜爽99麻豆av| 日日爽夜夜爽网站| 赤兔流量卡办理| 少妇人妻久久综合中文| 国产一级毛片在线| 国产亚洲5aaaaa淫片| 午夜老司机福利剧场| 尾随美女入室| 亚洲av男天堂| 精品99又大又爽又粗少妇毛片| 久久久亚洲精品成人影院| 国产免费一区二区三区四区乱码| 十八禁网站网址无遮挡 | a级片在线免费高清观看视频| 色婷婷av一区二区三区视频| 欧美 日韩 精品 国产| 日韩欧美 国产精品| 在线观看免费高清a一片| 成人午夜精彩视频在线观看| 一本色道久久久久久精品综合| 亚洲国产欧美日韩在线播放 | 亚洲精品456在线播放app| 成人毛片60女人毛片免费| 啦啦啦视频在线资源免费观看| av专区在线播放| 精品一区二区免费观看| 18禁动态无遮挡网站| a级毛片在线看网站| 欧美日韩国产mv在线观看视频| 人人妻人人爽人人添夜夜欢视频 | 乱人伦中国视频| 国产精品偷伦视频观看了| 高清欧美精品videossex| 国产一区二区三区综合在线观看 | 我的女老师完整版在线观看| 久久久久久久国产电影| freevideosex欧美| 黄色视频在线播放观看不卡| 亚洲国产精品一区二区三区在线| 精品卡一卡二卡四卡免费| 亚洲性久久影院| av又黄又爽大尺度在线免费看| 蜜桃久久精品国产亚洲av| 亚洲真实伦在线观看| 国产精品嫩草影院av在线观看| 国产精品一区二区三区四区免费观看| 国产综合精华液| a级一级毛片免费在线观看| 最近最新中文字幕免费大全7| 免费看日本二区| 日韩一本色道免费dvd| 夜夜骑夜夜射夜夜干| www.av在线官网国产| 亚洲欧洲精品一区二区精品久久久 | videossex国产| 亚洲在久久综合| 黑人猛操日本美女一级片| 亚洲国产精品国产精品| 久久精品国产自在天天线| 亚洲欧美日韩另类电影网站| 亚洲av日韩在线播放| 精品卡一卡二卡四卡免费| 不卡视频在线观看欧美| 免费大片18禁| 午夜福利,免费看| 2018国产大陆天天弄谢| 大又大粗又爽又黄少妇毛片口| 九九久久精品国产亚洲av麻豆| 免费看日本二区| 日日啪夜夜撸| 午夜激情福利司机影院| 成人亚洲欧美一区二区av| 久久青草综合色| 熟女av电影| 色哟哟·www| 午夜影院在线不卡| 亚洲自偷自拍三级| 精品久久久噜噜| 亚洲av男天堂| 精品亚洲成a人片在线观看| 日韩精品有码人妻一区| 狂野欧美白嫩少妇大欣赏| 夜夜看夜夜爽夜夜摸| 六月丁香七月| 夫妻性生交免费视频一级片| 好男人视频免费观看在线| 亚洲欧美日韩东京热| 久久久久人妻精品一区果冻| 少妇的逼水好多| 九九久久精品国产亚洲av麻豆| 亚洲av不卡在线观看| 亚洲真实伦在线观看| 久久精品夜色国产| 日韩av在线免费看完整版不卡| 国产亚洲精品久久久com| 久久久久久久久久人人人人人人| 在线观看美女被高潮喷水网站| 97超碰精品成人国产| 欧美国产精品一级二级三级 | 视频中文字幕在线观看| 十八禁网站网址无遮挡 | 国产av精品麻豆| 亚洲成人手机| 国产高清不卡午夜福利| 国产成人a∨麻豆精品| 亚洲国产日韩一区二区| 久久国产亚洲av麻豆专区| 伦精品一区二区三区| 欧美国产精品一级二级三级 | 最近中文字幕高清免费大全6| 一级毛片电影观看| 精品亚洲成a人片在线观看| 日韩av不卡免费在线播放| 极品人妻少妇av视频| 我要看黄色一级片免费的| 精品视频人人做人人爽| 肉色欧美久久久久久久蜜桃| 亚洲欧洲精品一区二区精品久久久 | 精品卡一卡二卡四卡免费| 久久婷婷青草| 国产熟女午夜一区二区三区 | 国产精品一二三区在线看| 久久久国产精品麻豆| 亚洲内射少妇av| 欧美日韩亚洲高清精品| 高清黄色对白视频在线免费看 | 久久久久久久久久成人| 国产成人91sexporn| 热99国产精品久久久久久7| 国产成人一区二区在线| 爱豆传媒免费全集在线观看| 妹子高潮喷水视频| 午夜福利网站1000一区二区三区| videossex国产| av福利片在线观看| 美女cb高潮喷水在线观看| av天堂中文字幕网| 亚洲欧美清纯卡通| 国产女主播在线喷水免费视频网站| 在现免费观看毛片| a级毛片在线看网站| 51国产日韩欧美| 久久久久久久久久久久大奶| 国产成人aa在线观看| 亚洲精品aⅴ在线观看| 在线 av 中文字幕| 国产永久视频网站| 最近的中文字幕免费完整| 亚洲激情五月婷婷啪啪| 国产亚洲一区二区精品| av免费观看日本| 妹子高潮喷水视频| 美女主播在线视频| 欧美成人午夜免费资源| av免费在线看不卡| 国产成人精品一,二区| 五月玫瑰六月丁香| 秋霞伦理黄片| 久久久久人妻精品一区果冻| 精品人妻熟女毛片av久久网站| 最近中文字幕高清免费大全6| 九草在线视频观看| 一区二区三区免费毛片| 欧美日韩一区二区视频在线观看视频在线| 男人爽女人下面视频在线观看| 天天操日日干夜夜撸| 国产精品99久久99久久久不卡 | 亚洲欧洲日产国产| 久久亚洲国产成人精品v| 人人澡人人妻人| 内射极品少妇av片p| 视频中文字幕在线观看| 日韩熟女老妇一区二区性免费视频| av卡一久久| 亚洲av二区三区四区| 在线观看免费高清a一片| 波野结衣二区三区在线| 亚洲av免费高清在线观看| 午夜福利网站1000一区二区三区| 国产中年淑女户外野战色| 国产精品偷伦视频观看了| 亚洲精品456在线播放app| 26uuu在线亚洲综合色| 九九爱精品视频在线观看| 色婷婷av一区二区三区视频| 青春草视频在线免费观看| 色哟哟·www| 免费看不卡的av| 免费大片18禁| 亚洲成色77777| 精品少妇内射三级| 久久人妻熟女aⅴ| 老司机亚洲免费影院| 一级,二级,三级黄色视频| 一级黄片播放器| 欧美日韩在线观看h| 午夜精品国产一区二区电影| 国产欧美日韩综合在线一区二区 | 免费观看a级毛片全部| 中文字幕免费在线视频6| 精品国产一区二区久久| 精品一区二区三卡| 国产精品一区二区在线不卡| 久久鲁丝午夜福利片| 精品国产国语对白av| av.在线天堂| 黄色视频在线播放观看不卡| 欧美精品一区二区免费开放| 亚洲美女搞黄在线观看| 亚洲国产欧美在线一区| 国产成人午夜福利电影在线观看| 久久6这里有精品| 色94色欧美一区二区| 亚洲欧美清纯卡通| 午夜日本视频在线| 久久久久久久大尺度免费视频| 少妇人妻久久综合中文| 日韩精品免费视频一区二区三区 | 欧美三级亚洲精品| 最近手机中文字幕大全| 久久久久久久久久人人人人人人| 美女福利国产在线| 久久久久久人妻| 青春草国产在线视频| 亚洲激情五月婷婷啪啪| 亚洲国产精品成人久久小说| 插逼视频在线观看| 久久97久久精品| 国产成人freesex在线| 偷拍熟女少妇极品色| 亚洲无线观看免费| 一级毛片电影观看| 久久精品久久精品一区二区三区| 亚洲精华国产精华液的使用体验| 日本猛色少妇xxxxx猛交久久| 一级爰片在线观看| 免费播放大片免费观看视频在线观看| 亚洲综合精品二区| 乱码一卡2卡4卡精品| 亚洲av.av天堂| 亚洲三级黄色毛片| 久久影院123| 性高湖久久久久久久久免费观看| 热99国产精品久久久久久7| 久久免费观看电影| 亚洲精品中文字幕在线视频 | 人妻人人澡人人爽人人| 日韩人妻高清精品专区| 中文字幕免费在线视频6| 99久久中文字幕三级久久日本| a级一级毛片免费在线观看| 久久99热这里只频精品6学生| .国产精品久久| 亚洲欧美中文字幕日韩二区| 精华霜和精华液先用哪个| 十八禁高潮呻吟视频 | 国产精品99久久久久久久久| 国模一区二区三区四区视频| 日韩av不卡免费在线播放| 18禁动态无遮挡网站| 爱豆传媒免费全集在线观看| 久久av网站| 国产高清有码在线观看视频| 九九在线视频观看精品| 人妻系列 视频| 99热国产这里只有精品6| 免费少妇av软件| 一级二级三级毛片免费看| 久久国产精品大桥未久av | 婷婷色综合大香蕉| 精品久久久噜噜| 一级毛片久久久久久久久女| 欧美区成人在线视频| 狂野欧美白嫩少妇大欣赏| 婷婷色综合大香蕉| 成人综合一区亚洲| 亚洲色图综合在线观看| 国产精品一区二区性色av| 男人狂女人下面高潮的视频| 亚洲自偷自拍三级| 99国产精品免费福利视频| 午夜福利影视在线免费观看| 亚洲精品日韩在线中文字幕| 久久精品国产鲁丝片午夜精品| 高清av免费在线| 一个人看视频在线观看www免费| 最近手机中文字幕大全| 国产精品成人在线| av在线app专区| 男女国产视频网站| 国产黄片美女视频| 黄色视频在线播放观看不卡| 日韩精品有码人妻一区| 午夜视频国产福利| 91精品伊人久久大香线蕉| 边亲边吃奶的免费视频| 国产欧美日韩一区二区三区在线 | 国产片特级美女逼逼视频| 少妇人妻一区二区三区视频| 欧美精品高潮呻吟av久久| 青青草视频在线视频观看| 免费大片黄手机在线观看| 热re99久久国产66热| 九九在线视频观看精品| 亚洲,一卡二卡三卡| 欧美 亚洲 国产 日韩一| 欧美激情极品国产一区二区三区 | 久久婷婷青草| av免费在线看不卡| 观看免费一级毛片| 亚洲精品乱久久久久久| 夜夜爽夜夜爽视频| 五月玫瑰六月丁香| 午夜免费观看性视频| 亚洲精品国产色婷婷电影| 蜜桃在线观看..| 少妇裸体淫交视频免费看高清| 国产在线视频一区二区| 亚洲,一卡二卡三卡| av福利片在线观看| 日韩欧美精品免费久久| 你懂的网址亚洲精品在线观看| 欧美 亚洲 国产 日韩一| 22中文网久久字幕| 丰满饥渴人妻一区二区三| 黄片无遮挡物在线观看| 国产无遮挡羞羞视频在线观看| 精品久久久久久久久亚洲| 成年av动漫网址| 一级毛片 在线播放| 国产淫语在线视频| 日韩强制内射视频| 欧美日本中文国产一区发布| 一级毛片我不卡| 在线看a的网站| 国产精品99久久99久久久不卡 | 亚洲欧美日韩卡通动漫| 国产女主播在线喷水免费视频网站| 狠狠精品人妻久久久久久综合| 多毛熟女@视频| 亚洲成色77777| 国产精品一二三区在线看| 一本—道久久a久久精品蜜桃钙片| 欧美激情国产日韩精品一区| av女优亚洲男人天堂| 国产伦精品一区二区三区视频9| 蜜桃在线观看..| 久久国内精品自在自线图片| 如日韩欧美国产精品一区二区三区 | 国产一级毛片在线| 免费久久久久久久精品成人欧美视频 | 精品卡一卡二卡四卡免费| 日韩av在线免费看完整版不卡| 精品人妻一区二区三区麻豆| 搡老乐熟女国产| 制服丝袜香蕉在线| 亚洲人成网站在线播| 国产亚洲最大av| 日本免费在线观看一区| 晚上一个人看的免费电影| 国产亚洲91精品色在线| 色94色欧美一区二区| 丰满饥渴人妻一区二区三| 大片电影免费在线观看免费| 夫妻性生交免费视频一级片| 色视频www国产| 免费久久久久久久精品成人欧美视频 | 一本大道久久a久久精品| 日本午夜av视频| 另类亚洲欧美激情| 国语对白做爰xxxⅹ性视频网站| 国产精品无大码| 国产色婷婷99| av国产久精品久网站免费入址| 街头女战士在线观看网站| 成年女人在线观看亚洲视频| 一区二区三区免费毛片| 日本av免费视频播放| 国产 一区精品| 亚洲美女搞黄在线观看| 久久综合国产亚洲精品| 中文天堂在线官网| 老女人水多毛片| av在线播放精品|