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

    捕集數(shù)對(duì)碳酸鹽巖油藏生物聚合物驅(qū)采收率的影響

    2022-08-25 08:50:28ALSHALABIEmad
    石油勘探與開(kāi)發(fā) 2022年4期
    關(guān)鍵詞:采收率巖心滲透率

    AL-SHALABI Emad W.

    (阿聯(lián)酋哈利法科技大學(xué)石油工程學(xué)院,阿布扎比 127788,阿聯(lián)酋)

    0 引言

    在一次采油和二次采油后,約 60%的原始石油地質(zhì)儲(chǔ)量仍滯留在儲(chǔ)集層中,需要應(yīng)用溶劑法、熱力法、化學(xué)驅(qū)法等方法進(jìn)行三次開(kāi)發(fā),此類(lèi)提高采收率方法主要開(kāi)采未波及的原油以及被毛細(xì)管捕集的殘余油[1],其中聚合物驅(qū)屬于成熟的化學(xué)驅(qū)方法。全球40多個(gè)成功的聚合物驅(qū)項(xiàng)目大多應(yīng)用在溫度低于60 ℃、地層水礦化度低于 1×105mg/L、滲透率大于40×10-3μm2的中高滲透砂巖儲(chǔ)集層中[2]。中東地區(qū)碳酸鹽巖儲(chǔ)集層具有地層溫度高(高于85 ℃)、地層水礦化度高(大于 1×105mg/L)、非均質(zhì)性強(qiáng)、滲透率低(小于100×10-3μm2)的特點(diǎn)[3-5],儲(chǔ)集層為混合潤(rùn)濕—油濕,原油采收率低,大部分提高采收率技術(shù)均難以實(shí)施[6-7]。因此,研究人員一直致力于擴(kuò)大聚合物驅(qū)的應(yīng)用范圍,使其能適應(yīng)惡劣條件下的碳酸鹽巖,主要研究方向包括開(kāi)發(fā)新型聚合物,以及將聚合物驅(qū)與其他提高采收率技術(shù)相結(jié)合[8-9]。

    聚合物根據(jù)其來(lái)源分為合成聚合物和生物聚合物。由于生物聚合物含剛性分子鏈而合成聚合物含柔性分子鏈,高礦化度溶液中生物聚合物性質(zhì)優(yōu)于合成聚合物[10-11]。生物聚合物包括黃原膠、硬葡聚糖和裂褶菌多糖等,其中裂褶菌多糖具有與多糖親水性生物聚合物相同的化學(xué)特性和流變特性,且能耐受高礦化度和硬度,抗機(jī)械降解;裂褶菌多糖由于具有三螺旋分子構(gòu)象,熱穩(wěn)定性更好[12-13];裂褶菌多糖在厭氧條件下能在長(zhǎng)達(dá)8個(gè)月內(nèi)保持黏度。Quadri[14]研究指出,裂褶菌多糖在溫度為120 ℃、礦化度為2×105mg/L的儲(chǔ)集層中具有應(yīng)用前景。但生物聚合物在惡劣儲(chǔ)集層條件下的應(yīng)用一直集中于實(shí)驗(yàn)研究方面,關(guān)于生物聚合物行為的模擬研究很少。

    Zeito[15]開(kāi)發(fā)了三維油藏模擬器,研究了均質(zhì)和非均質(zhì)油藏中聚合物驅(qū)對(duì)提高垂向波及效率的影響;Bondor等[16]建立了聚合物的流變行為模型,研究了近井地帶對(duì)聚合物注入性的影響;Goudarzi等[17]總結(jié)了數(shù)值模擬中聚合物驅(qū)模塊常用的主要函數(shù),包括聚合物黏度、聚合物吸附量、滲透率降低系數(shù)、不可及孔隙體積、礦化度和硬度。近年來(lái)鮮有分析聚合物對(duì)提高惡劣條件下碳酸鹽巖儲(chǔ)集層采收率影響的數(shù)值模擬研究。筆者[18]使用德克薩斯大學(xué)奧斯汀分校開(kāi)發(fā)的化學(xué)驅(qū)數(shù)值模擬器 UTCHEM[19]研究了裂褶菌多糖生物聚合物對(duì)碳酸鹽巖巖心采收率的影響,結(jié)果表明,與水驅(qū)相比,生物聚合物驅(qū)可提高碳酸鹽巖采收率。

    本文基于文獻(xiàn)[18]的生物聚合物驅(qū)油巖心數(shù)值模型,對(duì)碳酸鹽巖油藏裂褶菌多糖生物聚合物驅(qū)進(jìn)行模擬研究,對(duì)巖心驅(qū)替實(shí)驗(yàn)數(shù)據(jù)進(jìn)行歷史擬合,分析降低殘余油飽和度的主要影響因素。

    1 數(shù)值模型與數(shù)據(jù)

    本文采用文獻(xiàn)[18]的碳酸鹽巖巖心數(shù)值模型(見(jiàn)圖1),對(duì)裂褶菌多糖生物聚合物驅(qū)進(jìn)行數(shù)值模擬研究。

    圖1 模擬所用的碳酸鹽巖巖心數(shù)值模型[18]

    裂褶菌多糖生物聚合物平均相對(duì)分子質(zhì)量為(2~6)×106,由重復(fù)的結(jié)構(gòu)單元組成,每個(gè)結(jié)構(gòu)單元由3個(gè)葡萄糖單元和1個(gè)支鏈葡萄糖單元組成(見(jiàn)圖2)[20]。裂褶菌多糖生物聚合物在120 ℃及高達(dá)2×105mg/L的礦化度條件下表現(xiàn)出剪切稀化行為和強(qiáng)熱穩(wěn)定性[14]。裂褶菌多糖對(duì)滲透率高于 30×10-3μm2的巖心注入能力良好,其在滲透率為(3~163)×10-3μm2的巖心上的動(dòng)態(tài)吸附量為7~48 μg/g。Li[21]通過(guò)碳酸鹽巖巖心驅(qū)替實(shí)驗(yàn)得出,在三次采油模式下,質(zhì)量濃度為200 mg/L的裂褶菌多糖生物聚合物提高采收率幅度達(dá) 7~10個(gè)百分點(diǎn)。本文數(shù)值模擬研究采用文獻(xiàn)[14]中的裂褶菌多糖生物聚合物溶液特性參數(shù)和文獻(xiàn)[21]的三次采油驅(qū)替方式。

    圖2 裂褶菌多糖生物聚合物的化學(xué)結(jié)構(gòu)[20]

    Li[21]的巖心驅(qū)替實(shí)驗(yàn)在水驅(qū)后進(jìn)行生物聚合物驅(qū),實(shí)驗(yàn)步驟為:巖心飽和地層水,然后飽和地層原油造束縛水,而后在油藏條件(120 ℃,20 684 kPa)下注入地層水進(jìn)行水驅(qū),最后采用生物聚合物驅(qū)進(jìn)一步提高巖心采收率。驅(qū)替實(shí)驗(yàn)中巖心和流體的物性參數(shù)如表1所示,更多數(shù)據(jù)信息參考文獻(xiàn)[14]和[21]。文獻(xiàn)[18]對(duì)該驅(qū)替實(shí)驗(yàn)的水驅(qū)階段進(jìn)行了歷史擬合(見(jiàn)圖3),本文主要對(duì)該實(shí)驗(yàn)的三次采油階段進(jìn)行歷史擬合。

    表1 驅(qū)替實(shí)驗(yàn)中巖心和流體的物性參數(shù)

    圖3 巖心驅(qū)替實(shí)驗(yàn)水驅(qū)階段歷史擬合曲線[18](PV—孔隙體積)

    2 生物聚合物特性模擬

    2.1 生物聚合物黏度

    生物聚合物黏度是生物聚合物質(zhì)量濃度和注入水礦化度的函數(shù),利用Flory-Huggins方程描述[22]:

    該式描述了生物聚合物黏度隨礦化度增加及質(zhì)量濃度降低而降低的性質(zhì)。溫度對(duì)生物聚合物黏度沒(méi)有明顯影響,當(dāng)溫度升至150 ℃時(shí)該生物聚合物仍具有熱穩(wěn)定性[14,21],本文采用(1)式對(duì) 120 ℃下的生物聚合物黏度進(jìn)行了計(jì)算。用Csep描述注入水礦化度和硬度對(duì)生物聚合物黏度的影響,其定義為:

    圖4a和圖4b顯示了生物聚合物質(zhì)量濃度和注入水礦化度對(duì)生物聚合物黏度的影響?;A(chǔ)案例中生物聚合物質(zhì)量濃度為200 mg/L,地層水礦化度為168 948 mg/L,地層水中 Na+,Cl-,Ca2+,Mg2+,SO42-,HCO3-,K+的質(zhì)量濃度分別為53 394,102 989,9 260,1 450,0,1 107,748 mg/L。

    圖4 120 ℃下生物聚合物的模擬數(shù)據(jù)(綠色點(diǎn)代表基礎(chǔ)案例的相應(yīng)數(shù)據(jù))

    該方程描述了生物聚合物溶液的剪切稀化行為,生物聚合物黏度隨著剪切速率的增加而降低。在(3)式中,μp,0為低剪切極限(接近零)時(shí)的生物聚合物黏度,μw為高剪切極限(接近無(wú)窮大)時(shí)的水黏度。等效剪切速率用Cannella方程定義[24]:

    圖4c顯示了剪切速率對(duì)生物聚合物黏度的影響,當(dāng)剪切速率從10 s-1增加到1 000 s-1時(shí),生物聚合物黏度從 20 mPa·s減小到 0.4 mPa·s。其中,10 s-1和1 000 s-1分別代表油藏和注入井附近的剪切速率。

    2.2 剪切效應(yīng)

    剪切對(duì)生物聚合物黏度的影響使用 Meter方程描述[23]:

    2.3 生物聚合物吸附量

    與機(jī)械捕集和流體動(dòng)力學(xué)捕集相比,生物聚合物吸附是最主要的生物聚合物滯留機(jī)制。聚合物吸附取決于巖石表面、注入水溶液和聚合物本身的性質(zhì)[11]。本文研究中生物聚合物吸附用Langmuir等溫吸附曲線描述[25]:

    2.4 滲透率降低系數(shù)

    生物聚合物吸附可能導(dǎo)致滲透率顯著降低,采用滲透率降低系數(shù)(FKr)表示:

    模擬器中該系數(shù)使用以下方程計(jì)算:

    其中,假定FKr,max經(jīng)驗(yàn)值為10。

    2.5 不可及孔隙體積

    在生物聚合物驅(qū)過(guò)程中,生物聚合物無(wú)法波及儲(chǔ)集層中的所有孔隙,因此引入不可及孔隙體積參數(shù)Vip,定義為半徑小于生物聚合物顆粒半徑的孔隙體積占總孔隙體積的比例。尤其使用高相對(duì)分子質(zhì)量生物聚合物時(shí),由于聚合物單體分子尺寸更大,不可及孔隙體積也更大。不可及孔隙體積對(duì)生物聚合物波及效率既有積極影響,也有消極影響:積極方面包括促進(jìn)生物聚合物溶液驅(qū)替前緣的推進(jìn)并減少聚合物滯留量,消極方面在于不能波及不可及孔隙中的油[26]。模擬器中通過(guò)將生物聚合物守恒方程中的孔隙度乘以輸入?yún)?shù)Vap(關(guān)鍵字EPHI4)考慮Vip的影響:

    2.6 捕集數(shù)

    3 數(shù)值模擬結(jié)果與討論

    Li[21]的巖心驅(qū)替實(shí)驗(yàn)表明,裂褶菌多糖生物聚合物驅(qū)與水驅(qū)相比,采收率增加了 7~10個(gè)百分點(diǎn)。在巖心驅(qū)替過(guò)程中,先注入3 PV(孔隙體積倍數(shù))水,然后注入3 PV生物聚合物。文獻(xiàn)[18]對(duì)該驅(qū)替實(shí)驗(yàn)的水驅(qū)階段進(jìn)行了歷史擬合(見(jiàn)圖3),然后對(duì)三次采油進(jìn)行模擬預(yù)測(cè),結(jié)果顯示注入生物聚合物后采收率增加了 2個(gè)百分點(diǎn),實(shí)驗(yàn)結(jié)果與數(shù)值模擬結(jié)果差異大。本文基于敏感性分析,對(duì)生物聚合物驅(qū)階段實(shí)驗(yàn)數(shù)據(jù)進(jìn)行歷史擬合。

    對(duì)不可及孔隙體積、生物聚合物吸附量、滲透率降低系數(shù)、剪切速率系數(shù)和注入水的硬度等變量進(jìn)行了敏感性分析,以了解其對(duì)生物聚合物驅(qū)提高采收率的影響程度。表2為各變量對(duì)應(yīng)的 UTCHEM模擬器輸入?yún)?shù)及其不同取值,各設(shè)定值均基于實(shí)驗(yàn)和以往文獻(xiàn)的典型值確定。本文模擬研究中溫度為120 ℃,生物聚合物質(zhì)量濃度為 200 mg/L,注入水礦化度為168 948 mg/L、硬度為10 710 mg/L。

    表2 生物聚合物驅(qū)三次采油的敏感性分析參數(shù)

    3.1 不可及孔隙體積

    在基礎(chǔ)案例中,假設(shè)不可及孔隙體積分?jǐn)?shù)為零(即Vip=0),可以通過(guò)將UTCHEM模擬器中的波及孔隙體積分?jǐn)?shù)(Vap)設(shè)置為 1來(lái)實(shí)現(xiàn)。根據(jù)表3調(diào)整Vap參數(shù)值,對(duì)不可及孔隙體積進(jìn)行敏感性分析。不可及孔隙體積對(duì)生物聚合物驅(qū)三次采油的影響如表3和圖5所示。

    表3 不可及孔隙體積敏感性分析結(jié)果

    圖5 不可及孔隙體積對(duì)生物聚合物驅(qū)采收率的影響

    表3中Vip不同取值下的采收率與基礎(chǔ)案例采收率差異微小,圖5中所有模擬曲線幾乎重疊,表明生物聚合物驅(qū)采收率對(duì)不可及孔隙體積不敏感。總體來(lái)看,不可及孔隙體積的影響不大,但隨著不可及孔隙體積的增加,采油速度提高,而生物聚合物驅(qū)采收率略有下降。這是因?yàn)椴豢杉翱紫扼w積的增大使生物聚合物驅(qū)替前緣推動(dòng)加快,減少了生物聚合物滯留,但同時(shí)有大量石油殘留在不可及孔隙中。

    3.2 生物聚合物吸附量

    在UTCHEM模擬器中采用Langmuir模型描述裂褶菌多糖生物聚合物的吸附,調(diào)整參數(shù)ap1以描述吸附的生物聚合物質(zhì)量濃度。在基礎(chǔ)案例中,ap1=1.857,生物聚合物在120 ℃下的吸附量為6.9 μg/g。根據(jù)表4調(diào)整ap1,對(duì)生物聚合物吸附量進(jìn)行敏感性分析。

    表4 生物聚合物吸附量的敏感性分析結(jié)果

    表4和圖6都表明生物聚合物吸附對(duì)提高采收率略有負(fù)面影響。生物聚合物吸附量的增加導(dǎo)致提高采收率幅度的降低,這是因?yàn)楸晃降纳锞酆衔锩撾x水相,從而不能發(fā)揮生物聚合物的增黏作用。但這種負(fù)面效果并不明顯,吸附量從零增加至20 μg/g時(shí)相對(duì)水驅(qū)的采收率增幅降低0.33個(gè)百分點(diǎn);吸附量增加至200 μg/g時(shí)相對(duì)水驅(qū)的采收率增幅降低1.36個(gè)百分點(diǎn)。從圖6可以看出無(wú)吸附時(shí)為最佳方案,而如果實(shí)驗(yàn)中吸附量高達(dá) 200 μg/g,則說(shuō)明該生物聚合物驅(qū)油表現(xiàn)不佳且經(jīng)濟(jì)性差,不能用于油田規(guī)模的驅(qū)替。

    圖6 生物聚合物吸附量對(duì)采收率的影響

    3.3 滲透率降低系數(shù)

    在UTCHEM模擬器中調(diào)整輸入?yún)?shù)ckr模擬生物聚合物導(dǎo)致的水相有效滲透率降低對(duì)采收率的影響。對(duì)于基礎(chǔ)案例,ckr=0.008 83 m,滲透率降低系數(shù)FKr為3.9。根據(jù)表5調(diào)整ckr,模擬結(jié)果見(jiàn)表5和圖7。

    表5 滲透率降低系數(shù)的敏感性分析結(jié)果

    圖7 滲透率降低系數(shù)對(duì)生物聚合物驅(qū)采收率的影響

    結(jié)果表明,隨著滲透率降低系數(shù)的增加,生物聚合物驅(qū)提高采收率幅度增加。這是因?yàn)樗挠行B透率降低,導(dǎo)致含水率降低、采收率提高。如表5和圖7所示,不考慮滲透率降低(FKr=1)情況下的采收率與其他FKr取值情況下的采收率存在顯著差異。需要注意的是,滲透率降低系數(shù)通常不會(huì)超過(guò)10。

    3.4 剪切速率系數(shù)

    在UTCHEM模擬器中通過(guò)輸入?yún)?shù)cγ˙描述剪切速率對(duì)生物聚合物黏度的影響,cγ˙取決于剪切速率系數(shù)(C)。根據(jù) Kulawardana等[13]的研究,基礎(chǔ)案例中假設(shè)剪切速率系數(shù)為2.55。根據(jù)表6調(diào)整參數(shù)cγ˙,對(duì)剪切速率系數(shù)進(jìn)行敏感性分析。剪切速率系數(shù)對(duì)生物聚合物驅(qū)采收率的影響如表6和圖8所示。

    表6 剪切速率系數(shù)的敏感性分析結(jié)果

    圖8 剪切速率系數(shù)對(duì)生物聚合物驅(qū)采收率的影響

    表6中的采收率差異及圖8幾乎重疊的曲線表明生物聚合物驅(qū)對(duì)剪切速率不敏感,當(dāng)剪切速率系數(shù)增加到 1以上時(shí),采收率不會(huì)受到太大影響。另外,模擬水驅(qū)持續(xù)注水6 PV,對(duì)比分析高剪切速率系數(shù)下生物聚合物驅(qū)對(duì)提高采收率是否仍有效果。結(jié)果表明即使在高剪切速率系數(shù)下,與連續(xù)注水6 PV相比,生物聚合物驅(qū)仍可以提高采收率,證明在本研究的剪切系數(shù)取值下,生物聚合物沒(méi)有完全失去增黏能力。

    當(dāng)剪切速率系數(shù)為零時(shí),即忽略剪切速率效應(yīng)時(shí),生物聚合物驅(qū)采收率最高,與實(shí)驗(yàn)數(shù)據(jù)歷史擬合相對(duì)良好。然而這種情況不合理且不現(xiàn)實(shí),由于生物聚合物分子在多孔介質(zhì)中的滑移效應(yīng),一定會(huì)產(chǎn)生有效剪切,尤其是在非均質(zhì)性相對(duì)較強(qiáng)的地層中剪切更為明顯。此外,隨著剪切速率系數(shù)的增大,生物聚合物驅(qū)的采收率降低。因?yàn)殡S著剪切速率的增加,剪切稀化效應(yīng)導(dǎo)致生物聚合物在高剪切速率下的增黏能力下降,因此降低了流度控制能力和采收率。應(yīng)注意的是,生物聚合物黏度的這種損失通常是暫時(shí)的,其剪切作用可逆。

    3.5 注入水的硬度

    注入水的礦化度和硬度采用聚合物有效礦化度(Csep)進(jìn)行描述。水硬度參數(shù)βp用于確定二價(jià)陽(yáng)離子對(duì)生物聚合物有效礦化度的影響,在文獻(xiàn)中的典型值為10。當(dāng)βp等于1時(shí),一價(jià)和二價(jià)陽(yáng)離子的作用均等;然而,當(dāng)βp大于1時(shí),二價(jià)陽(yáng)離子對(duì)聚合物性能的影響更為明顯[11]?;A(chǔ)案例中取βp=1,代表忽略注入水硬度的影響。根據(jù)表7調(diào)整 UTCHEM 模擬器中的βp值[30-32],進(jìn)行水硬度的敏感性分析。

    表7 注入水硬度的敏感性分析結(jié)果

    在該分析中,總體水礦化度保持恒定,僅調(diào)整了硬度參數(shù)。從表7中極小的采收率差異以及圖9中幾乎重疊的曲線可以看出,隨著注入水硬度的增加,采收率略有下降,然而這種效果并不明顯。這是因?yàn)榱疡蘧嗵巧锞酆衔锸且环N非離子聚合物,對(duì)注入水的礦化度和硬度不敏感。

    圖9 注入水硬度對(duì)生物聚合物驅(qū)采收率的影響

    3.6 捕集數(shù)

    從上述不同參數(shù)的敏感性分析結(jié)果可以看出,通過(guò)調(diào)整這些參數(shù)無(wú)法擬合巖心驅(qū)替實(shí)驗(yàn)的采收率數(shù)據(jù),未能對(duì)生物聚合物驅(qū)提高采收率的控制因素進(jìn)行解釋?zhuān)时疚囊氩都瘮?shù),分析其對(duì)生物聚合物驅(qū)采收率的影響。

    在本文研究中,將捕集數(shù)簡(jiǎn)化為毛管數(shù),在注入3 PV水后注入3 PV生物聚合物,分別計(jì)算了水驅(qū)和生物聚合物驅(qū)的捕集數(shù)(見(jiàn)表8)。與水驅(qū)相比,生物聚合物驅(qū)期間的捕集數(shù)增加了約兩個(gè)數(shù)量級(jí)。捕集數(shù)的增加主要是由于溶液中加入生物聚合物導(dǎo)致水相黏度增加,捕集數(shù)增加很可能導(dǎo)致殘余油飽和度降低,尤其是在碳酸鹽巖中。因此,對(duì)捕集數(shù)與殘余油飽和度之間的關(guān)系進(jìn)行了研究。

    表8 水驅(qū)和生物聚合物驅(qū)階段的捕集數(shù)

    不同研究人員使用與Li[21]實(shí)驗(yàn)中相似的碳酸鹽巖巖心對(duì)捕集數(shù)與殘余油飽和度的關(guān)系進(jìn)行了研究[33-35],其中 Kamath等[34]研究發(fā)現(xiàn),注水期間的捕集數(shù)為3.77×10-8,水驅(qū)后殘余油飽和度約 25.1%;生物聚合物驅(qū)期間的捕集數(shù)為9.80×10-7,聚合物驅(qū)后殘余油飽和度下降到約10.0%。因此,生物聚合物驅(qū)的提高采收率幅度與捕集數(shù)的大小有關(guān)。

    本文研究中采用CDC模型描述捕集數(shù)對(duì)殘余油飽和度的影響(見(jiàn)(14)式)。對(duì)于碳酸鹽巖,通常假設(shè)Sor,high為零,Sor,low取水驅(qū)后殘余油飽和度25.1%;τ取文獻(xiàn)中碳酸鹽巖的典型值0.8;捕集參數(shù)To為100 000,通過(guò)擬合生物聚合物驅(qū)期間殘余油飽和度與捕集數(shù)的關(guān)系曲線獲得。通常對(duì)于砂巖,τ取值為1,捕集參數(shù)To為25 000,Sor,high為零,Sor,low為25.1%。

    CDC模型計(jì)算結(jié)果如圖10所示,本文研究中生物聚合物驅(qū)階段的捕集數(shù)為9.80×10-7,小于砂巖的臨界捕集數(shù)NT*,即殘余油飽和度曲線開(kāi)始下降時(shí)的捕集數(shù)。通常,砂巖的臨界捕集數(shù)為1×10-6~1×10-5,超過(guò)碳酸鹽巖的臨界捕集數(shù)至少1~2個(gè)數(shù)量級(jí)。通常碳酸鹽巖孔徑分布范圍比砂巖更大,非均質(zhì)性相對(duì)更強(qiáng),碳酸鹽巖的CDC曲線更平緩。因而,從碳酸鹽巖儲(chǔ)集層中開(kāi)采部分原油較容易,但剩余油很難開(kāi)采;而對(duì)于砂巖儲(chǔ)集層,只要超過(guò)臨界捕集數(shù)就可以開(kāi)采出大部分原油。因此,為了顯著降低碳酸鹽巖殘余油飽和度,捕集數(shù)需達(dá)到 1×10-5~1×10-4。

    圖10 砂巖和碳酸鹽巖的CDC曲線

    表9所示為捕集數(shù)小于臨界捕集數(shù)和大于臨界捕集數(shù)情況下計(jì)算的相對(duì)滲透率參數(shù)(見(jiàn)(15)—(18)式),圖11為兩種情況下的相對(duì)滲透率曲線。對(duì)于水相,假定為0.3,注水過(guò)程中為 0.063,nw,high假定為 1。此外,假定束縛水飽和度恒定,超過(guò)臨界捕集數(shù)后的油相端點(diǎn)相對(duì)滲透率和Corey指數(shù)no保持不變。

    表9 捕集數(shù)小于臨界捕集數(shù)和大于臨界捕集數(shù)情況下的相對(duì)滲透率參數(shù)

    圖11 捕集數(shù)小于臨界捕集數(shù)和大于臨界捕集數(shù)情況下的相對(duì)滲透率曲線

    將上述CDC模型模擬結(jié)果用于巖心驅(qū)替實(shí)驗(yàn)數(shù)據(jù)的歷史擬合中,擬合結(jié)果如圖12所示,可見(jiàn)考慮捕集數(shù)影響的模擬結(jié)果與Li[21]實(shí)驗(yàn)數(shù)據(jù)擬合更好。不考慮捕集數(shù)影響時(shí)生物聚合物驅(qū)相對(duì)水驅(qū)提高采收率幅度為3.90個(gè)百分點(diǎn),考慮捕集數(shù)影響時(shí)提高采收率幅度為8.85個(gè)百分點(diǎn),與Li[21]實(shí)驗(yàn)中得到的水驅(qū)后生物聚合物驅(qū)可以提高采收率 7~10個(gè)百分點(diǎn)的觀點(diǎn)吻合。圖13所示為巖心生物聚合物驅(qū)實(shí)驗(yàn)數(shù)據(jù)的歷史擬合效果,水驅(qū)階段和生物聚合物驅(qū)階段均成功擬合。

    圖12 捕集數(shù)對(duì)生物聚合物驅(qū)采收率的影響

    圖13 巖心驅(qū)替實(shí)驗(yàn)數(shù)據(jù)歷史擬合效果

    從表8看出,殘余油飽和度由低捕集數(shù)下的25.1%降低到高捕集數(shù)下的10.0%。捕集數(shù)的顯著影響是由于聚合物驅(qū)后水的黏度比聚合物驅(qū)前增加了大約兩個(gè)數(shù)量級(jí)。生物聚合物驅(qū)期間忽略捕集數(shù)影響可能會(huì)出現(xiàn)誤導(dǎo),尤其在以平緩CDC曲線為特征的碳酸鹽巖中更為明顯。捕集數(shù)對(duì)生物聚合物驅(qū)采收率的影響程度取決于目標(biāo)巖石類(lèi)型和生物聚合物黏度。需要說(shuō)明的是,本文重點(diǎn)研究生物聚合物提高驅(qū)替效率和波及效率的能力,僅考慮聚合物的彈性效應(yīng),不考慮黏彈性效應(yīng),因?yàn)榕c合成聚合物相比,生物聚合物通常不表現(xiàn)出黏彈性效應(yīng)。

    4 結(jié)論

    采用 UTCHEM 模擬器成功評(píng)價(jià)了碳酸鹽巖儲(chǔ)集層水驅(qū)后注入裂褶菌多糖生物聚合物驅(qū)油階段捕集數(shù)對(duì)采收率的影響,對(duì)巖心驅(qū)替實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了良好的歷史擬合。

    敏感性分析表明,在本文研究范圍內(nèi),不可及孔隙體積、生物聚合物吸附量、滲透率降低系數(shù)、剪切速率系數(shù)和注入水硬度均無(wú)法準(zhǔn)確描述生物聚合物驅(qū)對(duì)碳酸鹽巖的提高采收率效果;而采用毛細(xì)管去飽和曲線模型考慮捕集數(shù)影響時(shí),殘余油飽和度從不考慮其影響或低捕集數(shù)條件下的25.1%降低到10.0%,采收率擬合效果較好。在碳酸鹽巖儲(chǔ)集層生物聚合物驅(qū)模擬中,不能忽略捕集數(shù)的影響。為了顯著降低殘余油飽和度,本文研究條件下捕集數(shù)需達(dá)到1×10-5~1×10-4的數(shù)量級(jí)。

    致謝:感謝哈利法科技大學(xué)提供的基金資助(FSU-2018-26);感謝德克薩斯大學(xué)奧斯汀分校Kamy Sepehrnoori教授提供UTCHEM模擬器的使用權(quán)限。

    符號(hào)注釋?zhuān)?/p>

    ap——經(jīng)驗(yàn)常數(shù),無(wú)因次;ap1——擬合參數(shù),無(wú)因次;ap2——擬合參數(shù),L/mg;Ap1——擬合常數(shù),(L/mg)Sp+1;Ap2——擬合常數(shù),(L/mg)Sp+2;Ap3——擬合常數(shù),(L/mg)Sp+3;Ap4——擬合常數(shù),(L/mg)Sp;bKr——輸入?yún)?shù),L/mg;bp——經(jīng)驗(yàn)常數(shù),L/mg;cKr——UTCHEM模擬器中滲透率降低系數(shù)輸入?yún)?shù),m;C——剪切速率系數(shù),無(wú)因次;C11——溶液中水的體積分?jǐn)?shù),無(wú)因次;C51,C61——溶液中陰離子和二價(jià)陽(yáng)離子的質(zhì)量濃度,mg/L;Cp——生物聚合物質(zhì)量濃度,mg/L;p——被吸附的生物聚合物濃度,mg/L;Csep——聚合物有效礦化度,mg/L;D——深度,m;FKr——滲透率降低系數(shù),最大值取10,無(wú)因次;g——重力加速度,m/s2;K——滲透率,m2;K——滲透率張量,m2;Keff,p——聚合物驅(qū)時(shí)巖石的有效滲透率,m2;Keff,w——水驅(qū)時(shí)巖石的有效滲透率,m2;Kro,Krw——油相、水相相對(duì)滲透率,無(wú)因次;——油相、水相相對(duì)滲透率端點(diǎn)值,無(wú)因次;Kref——吸附測(cè)量實(shí)驗(yàn)中巖石的參考滲透率,m2;K——平均滲透率,m2;m——裂褶菌多糖生物聚合物結(jié)構(gòu)單元數(shù)量;n——Corey指數(shù),無(wú)因次;NB——邦德數(shù),無(wú)因次;Nc——毛管數(shù),無(wú)因次;NT——捕集數(shù),無(wú)因次;NT*——臨界捕集數(shù),無(wú)因次;S——飽和度,%;Sor——?dú)堄嘤惋柡投龋?;Sp——(μp,0-μw)/μw與Csep雙對(duì)數(shù)圖的斜率,本文取-0.33;Swi——束縛水飽和度,%;T——捕集參數(shù),無(wú)因次;v——巖心中的達(dá)西流動(dòng)速率,m/s;Vap——可及孔隙體積分?jǐn)?shù),無(wú)因次;Vip——不可及孔隙體積分?jǐn)?shù),無(wú)因次;x,y,z——坐標(biāo)軸,cm;α——經(jīng)驗(yàn)參數(shù),通過(guò)擬合黏度實(shí)驗(yàn)測(cè)量數(shù)據(jù)求取,無(wú)因次;βp——硬度參數(shù),通常約為10,無(wú)因次;γ˙c——UTCHEM模擬器中等效剪切速率輸入?yún)?shù),無(wú)因次;γ˙eq——等效剪切速率,s-1;1/2——聚合物溶液黏度為其初始值一半時(shí)的剪切速率,s-1;θ——接觸角,(°);?——驅(qū)替實(shí)驗(yàn)中巖心傾角,(°);μp——生物聚合物的黏度,mPa·s;μp,0——零剪切速率下的生物聚合物黏度,Pa·s;μw——水的黏度,Pa·s;ρ——密度,kg/m3;σow——水相和油相之間界面張力,N/m;τ——與巖石孔徑分布相關(guān)的參數(shù),無(wú)因次;φ——孔隙度,%;?pw——水相的壓力梯度,Pa;?Φw——水相的壓力勢(shì)梯度,Pa。下標(biāo):high——高捕集數(shù);Kr——滲透率降低;low——低捕集數(shù);max——最大值;o——油相;w——水相。

    猜你喜歡
    采收率巖心滲透率
    《油氣地質(zhì)與采收率》征稿簡(jiǎn)則
    《油氣地質(zhì)與采收率》征稿簡(jiǎn)則
    《油氣地質(zhì)與采收率》第六屆編委會(huì)
    《油氣地質(zhì)與采收率》征稿簡(jiǎn)則
    Osteotomized folded scapular tip free flap for complex midfacial reconstruction
    中煤階煤層氣井排采階段劃分及滲透率變化
    不同滲透率巖芯孔徑分布與可動(dòng)流體研究
    SAGD井微壓裂儲(chǔ)層滲透率變化規(guī)律研究
    一種頁(yè)巖巖心資料的保存方法
    化工管理(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
    一个人观看的视频www高清免费观看| 2021天堂中文幕一二区在线观| 97超碰精品成人国产| 三级国产精品片| 欧美日韩视频高清一区二区三区二| 日本色播在线视频| 成人高潮视频无遮挡免费网站| 18禁动态无遮挡网站| 少妇猛男粗大的猛烈进出视频 | 51国产日韩欧美| 欧美一区二区亚洲| 少妇人妻 视频| 日本-黄色视频高清免费观看| 久久久久久久亚洲中文字幕| 国产人妻一区二区三区在| 成人特级av手机在线观看| av在线观看视频网站免费| 九草在线视频观看| 国产亚洲av片在线观看秒播厂| 日本欧美国产在线视频| 亚洲国产欧美人成| 黄片无遮挡物在线观看| 麻豆成人av视频| 精品亚洲乱码少妇综合久久| 国产在视频线精品| 有码 亚洲区| 一区二区三区精品91| 网址你懂的国产日韩在线| videos熟女内射| 国产成人aa在线观看| 久久99精品国语久久久| 国产毛片在线视频| 91精品国产九色| 久久久久九九精品影院| 亚洲内射少妇av| 国产精品国产三级专区第一集| 免费大片黄手机在线观看| videos熟女内射| 免费黄色在线免费观看| 天美传媒精品一区二区| 又大又黄又爽视频免费| 欧美日韩综合久久久久久| 久热这里只有精品99| 国产亚洲一区二区精品| 少妇猛男粗大的猛烈进出视频 | 日本黄色片子视频| 午夜福利视频精品| 日韩三级伦理在线观看| 中文欧美无线码| 老司机影院毛片| 日韩一区二区视频免费看| 国产精品蜜桃在线观看| 麻豆成人午夜福利视频| 嫩草影院新地址| 久久人人爽人人爽人人片va| 精品人妻视频免费看| 熟女电影av网| 人妻一区二区av| av在线老鸭窝| 啦啦啦啦在线视频资源| 日韩免费高清中文字幕av| 国产精品精品国产色婷婷| 美女国产视频在线观看| 婷婷色麻豆天堂久久| 老司机影院毛片| 麻豆成人av视频| 真实男女啪啪啪动态图| 天堂俺去俺来也www色官网| 国语对白做爰xxxⅹ性视频网站| 欧美一区二区亚洲| 国产精品人妻久久久影院| 亚洲天堂av无毛| 日本黄大片高清| 天天一区二区日本电影三级| 一级爰片在线观看| 亚洲国产成人一精品久久久| 中文乱码字字幕精品一区二区三区| 少妇裸体淫交视频免费看高清| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲欧美日韩无卡精品| 全区人妻精品视频| 久久这里有精品视频免费| 高清在线视频一区二区三区| 99热全是精品| 最后的刺客免费高清国语| 青青草视频在线视频观看| 亚洲av福利一区| 亚洲欧洲国产日韩| 亚洲精品久久久久久婷婷小说| 欧美日韩综合久久久久久| 伦精品一区二区三区| 大又大粗又爽又黄少妇毛片口| 嘟嘟电影网在线观看| 国产精品久久久久久精品电影小说 | 国产探花在线观看一区二区| 蜜桃久久精品国产亚洲av| a级毛色黄片| 国产精品一区www在线观看| 男插女下体视频免费在线播放| 男男h啪啪无遮挡| 欧美激情在线99| 亚洲精品久久午夜乱码| 欧美性感艳星| 18禁在线播放成人免费| 亚洲av电影在线观看一区二区三区 | 王馨瑶露胸无遮挡在线观看| 一区二区三区四区激情视频| av在线天堂中文字幕| 我要看日韩黄色一级片| 在线观看国产h片| 欧美国产精品一级二级三级 | 成年女人在线观看亚洲视频 | 国产伦精品一区二区三区视频9| 亚洲自拍偷在线| 91午夜精品亚洲一区二区三区| 亚洲天堂国产精品一区在线| 男人爽女人下面视频在线观看| 一级毛片我不卡| 少妇猛男粗大的猛烈进出视频 | 亚洲av中文字字幕乱码综合| 黄色欧美视频在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产欧美日韩一区二区三区在线 | 色综合色国产| 可以在线观看毛片的网站| 黄色欧美视频在线观看| 日韩视频在线欧美| 亚洲三级黄色毛片| 精品一区二区三卡| 成人欧美大片| 联通29元200g的流量卡| 免费少妇av软件| 成人美女网站在线观看视频| av在线app专区| 国产久久久一区二区三区| 欧美日韩亚洲高清精品| 韩国av在线不卡| 亚洲人与动物交配视频| 国产精品久久久久久精品电影小说 | 搡女人真爽免费视频火全软件| 精品国产乱码久久久久久小说| 亚洲真实伦在线观看| 久久99精品国语久久久| 国产精品国产三级专区第一集| 亚洲激情五月婷婷啪啪| 精品人妻偷拍中文字幕| 国产高潮美女av| 精品人妻偷拍中文字幕| 欧美精品人与动牲交sv欧美| 22中文网久久字幕| 免费观看av网站的网址| 国产av码专区亚洲av| 插逼视频在线观看| 国产免费福利视频在线观看| 中文欧美无线码| 日本午夜av视频| 亚洲欧美日韩无卡精品| 美女主播在线视频| 九九久久精品国产亚洲av麻豆| 久久久成人免费电影| 亚洲自偷自拍三级| 欧美少妇被猛烈插入视频| 我要看日韩黄色一级片| 最近中文字幕高清免费大全6| 日日摸夜夜添夜夜爱| 性色avwww在线观看| 亚洲精品日本国产第一区| 少妇人妻精品综合一区二区| 免费黄色在线免费观看| 男人狂女人下面高潮的视频| 狠狠精品人妻久久久久久综合| 国产综合懂色| av在线app专区| 一级毛片aaaaaa免费看小| a级一级毛片免费在线观看| 亚洲性久久影院| 精品酒店卫生间| 亚洲精品一二三| 国产伦在线观看视频一区| 日韩,欧美,国产一区二区三区| 国产高清不卡午夜福利| 午夜福利网站1000一区二区三区| 婷婷色麻豆天堂久久| 欧美日韩精品成人综合77777| 亚洲国产精品专区欧美| 中文精品一卡2卡3卡4更新| 亚洲精品国产成人久久av| 综合色丁香网| 免费观看无遮挡的男女| 国产淫片久久久久久久久| 国产美女午夜福利| 国产乱来视频区| 久热久热在线精品观看| 伊人久久国产一区二区| 成年女人在线观看亚洲视频 | 日韩强制内射视频| 亚洲,一卡二卡三卡| 男人和女人高潮做爰伦理| 麻豆精品久久久久久蜜桃| 亚洲精品国产av蜜桃| 国产精品麻豆人妻色哟哟久久| 另类亚洲欧美激情| 少妇高潮的动态图| 日本-黄色视频高清免费观看| freevideosex欧美| 欧美日韩精品成人综合77777| 一本色道久久久久久精品综合| 一级毛片电影观看| 少妇猛男粗大的猛烈进出视频 | 国产有黄有色有爽视频| 久久热精品热| 免费播放大片免费观看视频在线观看| 日韩三级伦理在线观看| 搡老乐熟女国产| 性插视频无遮挡在线免费观看| videossex国产| 久久久久精品性色| 免费看日本二区| 亚洲人成网站高清观看| 大陆偷拍与自拍| 国产精品一及| 国产精品爽爽va在线观看网站| 热re99久久精品国产66热6| 少妇的逼水好多| 男女啪啪激烈高潮av片| 国产精品秋霞免费鲁丝片| 久久精品综合一区二区三区| 大片免费播放器 马上看| 亚洲精品国产av蜜桃| 亚洲精品一二三| www.色视频.com| 下体分泌物呈黄色| 国产又色又爽无遮挡免| 亚洲成人精品中文字幕电影| 日本av手机在线免费观看| 男人狂女人下面高潮的视频| 亚洲欧美日韩东京热| 少妇人妻 视频| 精品久久久久久久末码| 中文字幕人妻熟人妻熟丝袜美| 国产伦在线观看视频一区| 一级毛片 在线播放| 久久午夜福利片| 国产一区有黄有色的免费视频| 大话2 男鬼变身卡| 麻豆国产97在线/欧美| 蜜桃亚洲精品一区二区三区| 三级经典国产精品| 七月丁香在线播放| 亚洲国产精品成人久久小说| 夜夜爽夜夜爽视频| 欧美精品人与动牲交sv欧美| 国产国拍精品亚洲av在线观看| 亚洲人与动物交配视频| 国产又色又爽无遮挡免| 国产精品偷伦视频观看了| 日韩av免费高清视频| 少妇人妻精品综合一区二区| 特级一级黄色大片| 欧美日韩精品成人综合77777| 亚洲高清免费不卡视频| 国产精品久久久久久久久免| 在线看a的网站| 麻豆国产97在线/欧美| 麻豆乱淫一区二区| 欧美+日韩+精品| 黄片无遮挡物在线观看| 国产免费视频播放在线视频| 日本三级黄在线观看| 全区人妻精品视频| 亚洲精品乱久久久久久| 亚洲婷婷狠狠爱综合网| 麻豆乱淫一区二区| 国产色爽女视频免费观看| 国产毛片在线视频| 欧美精品一区二区大全| 国产成人91sexporn| 69人妻影院| 中文字幕亚洲精品专区| 在线精品无人区一区二区三 | 亚洲精华国产精华液的使用体验| 三级男女做爰猛烈吃奶摸视频| 亚洲精品日韩av片在线观看| 国产高清有码在线观看视频| 免费看日本二区| 成人欧美大片| 亚洲av二区三区四区| 日韩av在线免费看完整版不卡| 亚洲精品日韩av片在线观看| 国产成人福利小说| 欧美成人a在线观看| 人妻系列 视频| 又黄又爽又刺激的免费视频.| 亚洲欧美日韩东京热| 在线观看av片永久免费下载| 亚洲伊人久久精品综合| 午夜亚洲福利在线播放| 好男人在线观看高清免费视频| 欧美 日韩 精品 国产| 99久久精品一区二区三区| 免费av不卡在线播放| 免费看a级黄色片| 中文精品一卡2卡3卡4更新| 午夜福利视频精品| 亚洲精品第二区| 天堂网av新在线| 亚洲经典国产精华液单| 神马国产精品三级电影在线观看| 免费电影在线观看免费观看| 久久久久久久精品精品| 日韩亚洲欧美综合| 国产久久久一区二区三区| 国产伦理片在线播放av一区| 国产成人精品久久久久久| 亚洲av成人精品一二三区| 成人漫画全彩无遮挡| 国产精品嫩草影院av在线观看| 亚洲图色成人| 黄色配什么色好看| 国产精品偷伦视频观看了| 国产一区二区亚洲精品在线观看| 日韩免费高清中文字幕av| 蜜桃久久精品国产亚洲av| 国产毛片在线视频| 国产永久视频网站| av又黄又爽大尺度在线免费看| 大片免费播放器 马上看| 亚洲精品中文字幕在线视频 | 亚洲av日韩在线播放| 日日撸夜夜添| 午夜福利视频精品| 国产一级毛片在线| 国产成人午夜福利电影在线观看| 制服丝袜香蕉在线| 一本色道久久久久久精品综合| 六月丁香七月| 在线亚洲精品国产二区图片欧美 | 亚洲人成网站在线观看播放| 国产午夜精品一二区理论片| 国产免费一级a男人的天堂| av免费在线看不卡| 亚洲精品影视一区二区三区av| 女的被弄到高潮叫床怎么办| 国产精品人妻久久久久久| 午夜福利高清视频| 啦啦啦在线观看免费高清www| 最近手机中文字幕大全| 老女人水多毛片| 国产有黄有色有爽视频| 一本色道久久久久久精品综合| 亚洲va在线va天堂va国产| 久久久久久久久久久丰满| 日韩一本色道免费dvd| av在线天堂中文字幕| 成人鲁丝片一二三区免费| 熟女电影av网| 少妇的逼好多水| 成年人午夜在线观看视频| 日韩成人伦理影院| 日韩精品有码人妻一区| 亚洲人与动物交配视频| 三级经典国产精品| 五月玫瑰六月丁香| 日本黄色片子视频| 久久久久久久久久久免费av| 免费播放大片免费观看视频在线观看| 2021少妇久久久久久久久久久| 精品视频人人做人人爽| 神马国产精品三级电影在线观看| 在线观看国产h片| av线在线观看网站| 亚洲成人av在线免费| 亚洲国产欧美在线一区| 黄色视频在线播放观看不卡| 日韩欧美一区视频在线观看 | 女人被狂操c到高潮| av福利片在线观看| 高清av免费在线| 欧美3d第一页| 99精国产麻豆久久婷婷| 麻豆国产97在线/欧美| 国产精品一二三区在线看| 久久精品国产亚洲av涩爱| 亚洲av成人精品一二三区| 色婷婷久久久亚洲欧美| 成人午夜精彩视频在线观看| 久久99热这里只有精品18| 久久久久九九精品影院| 成人欧美大片| 亚洲欧美日韩无卡精品| 一个人观看的视频www高清免费观看| 极品少妇高潮喷水抽搐| 亚洲av国产av综合av卡| 超碰av人人做人人爽久久| 国产久久久一区二区三区| 免费观看av网站的网址| 搡老乐熟女国产| 亚洲av不卡在线观看| 亚洲欧美一区二区三区黑人 | 听说在线观看完整版免费高清| 涩涩av久久男人的天堂| 亚洲精品久久久久久婷婷小说| 欧美激情在线99| 亚洲人成网站在线观看播放| 亚洲精品色激情综合| 国产精品嫩草影院av在线观看| 国产熟女欧美一区二区| av卡一久久| 不卡视频在线观看欧美| 亚洲欧洲日产国产| 韩国av在线不卡| 久久这里有精品视频免费| 国产黄片美女视频| 黄片wwwwww| 亚洲怡红院男人天堂| 丰满少妇做爰视频| 尤物成人国产欧美一区二区三区| 在线a可以看的网站| 我的女老师完整版在线观看| 97热精品久久久久久| 日韩视频在线欧美| 亚洲成人中文字幕在线播放| 国产黄片视频在线免费观看| 久久这里有精品视频免费| 免费av毛片视频| 青春草国产在线视频| 国产成人精品久久久久久| 日本爱情动作片www.在线观看| 中文字幕制服av| 日日啪夜夜撸| 国内精品宾馆在线| 国产精品国产av在线观看| 少妇人妻一区二区三区视频| 丝袜美腿在线中文| 免费黄频网站在线观看国产| 两个人的视频大全免费| 听说在线观看完整版免费高清| 女人被狂操c到高潮| 伦理电影大哥的女人| 一区二区av电影网| 午夜福利视频1000在线观看| 女人久久www免费人成看片| 国产精品人妻久久久影院| 国产精品国产三级国产专区5o| 国产精品一区二区性色av| 毛片一级片免费看久久久久| 一级毛片aaaaaa免费看小| 久久精品国产自在天天线| 国精品久久久久久国模美| 熟妇人妻不卡中文字幕| 寂寞人妻少妇视频99o| 国产视频首页在线观看| 亚洲国产成人一精品久久久| 永久网站在线| 欧美日韩视频高清一区二区三区二| 亚洲av在线观看美女高潮| 久久热精品热| 亚洲国产精品成人综合色| 国产探花在线观看一区二区| 国产一区二区三区av在线| 亚洲最大成人手机在线| 成年版毛片免费区| 九色成人免费人妻av| 国产爽快片一区二区三区| 中文字幕av成人在线电影| av卡一久久| 国产精品国产av在线观看| 亚洲伊人久久精品综合| 三级国产精品片| 日本免费在线观看一区| 综合色丁香网| 中文字幕久久专区| 久久精品国产自在天天线| 熟女人妻精品中文字幕| 伦理电影大哥的女人| 精品一区在线观看国产| 国产日韩欧美亚洲二区| 国产精品偷伦视频观看了| 亚洲最大成人手机在线| 成人高潮视频无遮挡免费网站| 日韩av免费高清视频| 欧美成人午夜免费资源| 亚洲av国产av综合av卡| 午夜福利视频1000在线观看| 伊人久久精品亚洲午夜| 亚洲欧洲日产国产| 国产乱人偷精品视频| 91aial.com中文字幕在线观看| 免费观看的影片在线观看| 精品少妇黑人巨大在线播放| 国产成人福利小说| 97在线人人人人妻| 精品人妻视频免费看| 丝袜脚勾引网站| 成人亚洲欧美一区二区av| 少妇裸体淫交视频免费看高清| 精品视频人人做人人爽| 亚洲怡红院男人天堂| 国产亚洲一区二区精品| 纵有疾风起免费观看全集完整版| 欧美成人a在线观看| 国产成人一区二区在线| 777米奇影视久久| 国产久久久一区二区三区| 国产午夜精品久久久久久一区二区三区| 日韩亚洲欧美综合| 人人妻人人澡人人爽人人夜夜| 欧美一区二区亚洲| 亚洲丝袜综合中文字幕| av免费在线看不卡| 欧美极品一区二区三区四区| 国产亚洲5aaaaa淫片| 男人爽女人下面视频在线观看| 亚洲成人av在线免费| 国产av国产精品国产| 人人妻人人看人人澡| 啦啦啦在线观看免费高清www| 高清在线视频一区二区三区| 黄片wwwwww| 干丝袜人妻中文字幕| 亚洲精品亚洲一区二区| 国产亚洲精品久久久com| 亚洲色图av天堂| 亚洲精品久久午夜乱码| 亚洲欧美成人精品一区二区| 一级a做视频免费观看| 一级爰片在线观看| 亚洲美女视频黄频| 特大巨黑吊av在线直播| 免费看a级黄色片| 欧美xxxx性猛交bbbb| 日本黄大片高清| 丰满少妇做爰视频| 看十八女毛片水多多多| 久久精品久久精品一区二区三区| 国产伦理片在线播放av一区| 久久久久久久午夜电影| 青春草视频在线免费观看| 免费不卡的大黄色大毛片视频在线观看| 亚洲图色成人| 又大又黄又爽视频免费| 欧美成人精品欧美一级黄| 国产精品成人在线| 亚洲欧美日韩另类电影网站 | 国产探花在线观看一区二区| 各种免费的搞黄视频| 老师上课跳d突然被开到最大视频| 看非洲黑人一级黄片| 蜜桃亚洲精品一区二区三区| 街头女战士在线观看网站| 一级黄片播放器| 亚洲精品乱码久久久久久按摩| 国产精品嫩草影院av在线观看| 免费观看av网站的网址| 国语对白做爰xxxⅹ性视频网站| 一个人看的www免费观看视频| 亚洲aⅴ乱码一区二区在线播放| 十八禁网站网址无遮挡 | 久久鲁丝午夜福利片| 校园人妻丝袜中文字幕| 欧美日本视频| 狂野欧美激情性xxxx在线观看| 欧美变态另类bdsm刘玥| 国产成人福利小说| 日韩电影二区| 亚洲av福利一区| 精品久久久噜噜| 亚洲自偷自拍三级| 欧美一级a爱片免费观看看| 丰满人妻一区二区三区视频av| 涩涩av久久男人的天堂| 久久久久久伊人网av| 老师上课跳d突然被开到最大视频| 亚洲综合色惰| 老女人水多毛片| 各种免费的搞黄视频| 国产成人精品一,二区| 亚州av有码| 精品一区二区三卡| 麻豆国产97在线/欧美| 日韩欧美一区视频在线观看 | 日韩在线高清观看一区二区三区| 搡老乐熟女国产| 91狼人影院| 啦啦啦中文免费视频观看日本| www.av在线官网国产| 毛片一级片免费看久久久久| 波野结衣二区三区在线| 一区二区三区乱码不卡18| 免费观看在线日韩| 毛片女人毛片| 国产女主播在线喷水免费视频网站| 成年人午夜在线观看视频| 22中文网久久字幕| 99热全是精品| 亚洲欧美日韩无卡精品| 欧美极品一区二区三区四区| 草草在线视频免费看| 国产亚洲av片在线观看秒播厂| 国产精品人妻久久久影院| 成人欧美大片| 嫩草影院精品99| 天天一区二区日本电影三级| 男男h啪啪无遮挡| 一级毛片黄色毛片免费观看视频| 色婷婷久久久亚洲欧美| 男女边摸边吃奶| 久久精品久久久久久久性| 3wmmmm亚洲av在线观看| 久久精品国产鲁丝片午夜精品| 各种免费的搞黄视频| 美女被艹到高潮喷水动态| 在线观看国产h片| 看免费成人av毛片| 嫩草影院入口|