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

    基于離散縫洞網(wǎng)絡(luò)模型的縫洞型油藏混合模型

    2017-10-12 08:41:28嚴(yán)俠黃朝琴李陽姚軍樊冬艷
    關(guān)鍵詞:縫洞溶洞滲透率

    嚴(yán)俠,黃朝琴,李陽,姚軍,樊冬艷

    ?

    基于離散縫洞網(wǎng)絡(luò)模型的縫洞型油藏混合模型

    嚴(yán)俠1,黃朝琴1,李陽2,姚軍1,樊冬艷1

    (1. 中國石油大學(xué)(華東) 石油工程學(xué)院,山東青島,266580;2. 中國石化油田勘探開發(fā)事業(yè)部,北京,100728)

    針對存在大尺度導(dǎo)流裂縫的復(fù)雜縫洞型油藏,為了克服目前常用數(shù)學(xué)模型不能準(zhǔn)確高效對其進(jìn)行流動(dòng)模擬的缺陷,提出一種新的混合模型。首先根據(jù)縫洞的結(jié)構(gòu)特征劃分合適的正交網(wǎng)格,在每個(gè)網(wǎng)格單元內(nèi)建立離散縫洞網(wǎng)絡(luò)模型(DFVN),基于超樣本技術(shù)和體積平均方法求得網(wǎng)格單元的等效滲透率張量,對小尺度的縫洞結(jié)構(gòu)進(jìn)行等效處理;以形成的正交網(wǎng)格為基礎(chǔ)建立嵌入式離散裂縫模型,對大尺度導(dǎo)流裂縫進(jìn)行顯示處理,并基于模擬有限差分方法建立能準(zhǔn)確處理全張量形式滲透率的數(shù)值計(jì)算格式;最后通過單相穩(wěn)定滲流數(shù)值算例驗(yàn)證本文方法的正確性和優(yōu)越性。研究結(jié)果表明:對于縫洞結(jié)構(gòu)易發(fā)生變化的油藏,只需根據(jù)變化后的縫洞結(jié)構(gòu)進(jìn)行局部的參數(shù)修正,避免了類似于離散縫洞網(wǎng)絡(luò)模型的整體網(wǎng)格重構(gòu),有效提高了計(jì)算效率。

    縫洞型油藏;離散縫洞網(wǎng)絡(luò)模型;超樣本技術(shù);嵌入式離散裂縫模型;模擬有限差分

    縫洞型碳酸鹽巖油氣藏屬于非常規(guī)油氣藏,其儲(chǔ)量規(guī)模較大。在世界已探明儲(chǔ)量中,碳酸鹽巖油氣藏儲(chǔ)量約占50%,產(chǎn)量占60%以上。我國碳酸鹽巖油氣藏儲(chǔ)量也十分可觀,是今后增儲(chǔ)的主要領(lǐng)域[1]。然而,由于這類油氣藏構(gòu)造及流動(dòng)模式十分復(fù)雜,因此,如何準(zhǔn)確模擬縫洞型油藏中的流體流動(dòng)成為一項(xiàng)具有挑戰(zhàn)性的研究工作。目前用于縫洞型油藏流動(dòng)數(shù)值模擬的數(shù)學(xué)模型主要有等效連續(xù)介質(zhì)模型[2?3]、多重介質(zhì)模型[4?7]以及離散縫洞網(wǎng)絡(luò)模型[8?9]。等效連續(xù)介質(zhì)模型是以等效滲透率為基礎(chǔ),用連續(xù)介質(zhì)方法描述滲流問題的數(shù)學(xué)模型。其突出優(yōu)點(diǎn)是可以利用各向異性連續(xù)介質(zhì)理論進(jìn)行分析。該模型的關(guān)鍵問題是等效滲透率的求取以及全張量的數(shù)值模擬方法。目前,縫洞型介質(zhì)等效滲透率的計(jì)算方法主要是基于流動(dòng)模擬的方法,可以分為2類:一類是基于宏觀達(dá)西定律的流量等效方法[10?12],另一類是基于均化理論的尺度升級方法[13?14]。這些方法在處理小尺度縫洞時(shí)可以得到較準(zhǔn)確的結(jié)果,但處理大尺度導(dǎo)流裂縫時(shí)往往會(huì)低估其強(qiáng)導(dǎo)流能力,致使得到的結(jié)果誤差較大[15]。多重介質(zhì)模型適用于油田規(guī)模數(shù)值模擬,在一定程度上反映了縫洞型油藏中的優(yōu)先流現(xiàn)象,并考慮裂縫系統(tǒng)、基巖系統(tǒng)和溶洞系統(tǒng)間的流體交換,比較符合實(shí)際。但被裂隙分割的基巖和溶洞系統(tǒng)被假定具有相同的大小和形狀,過于簡化而不能充分表現(xiàn)縫洞型油藏的強(qiáng)各向異性、不連續(xù)性等特征。同時(shí)該模型竄流系數(shù)難以確定,并且不適用于含有若干不連通大尺度導(dǎo)流縫洞的油藏。離散縫洞網(wǎng)絡(luò)模型將縫洞顯式表征,基巖和裂縫系統(tǒng)視為滲流區(qū)域,溶洞系統(tǒng)視為自由流區(qū)域,是離散裂縫模型的有效延伸,屬于宏觀尺度的精細(xì)模型,能準(zhǔn)確刻畫出縫洞型油藏中的流動(dòng)特征,但受目前計(jì)算機(jī)硬件和數(shù)值模擬技術(shù)的制約,只能處理小規(guī)模的縫洞型油藏。鑒于以上各模型對存在大尺度裂縫的縫洞型油藏難以進(jìn)行高效準(zhǔn)確的流動(dòng)模擬,本文作者將離散介質(zhì)模型與等效連續(xù)介質(zhì)模型相結(jié)合:對于小尺度裂縫和溶洞,基于離散縫洞網(wǎng)絡(luò)模型運(yùn)用超樣本技術(shù)和體積平均方法計(jì)算等效滲透率張量,對其進(jìn)行等效處理;對于大尺度裂縫,采用嵌入式離散裂縫模 型[15?17]進(jìn)行顯示處理,在減小計(jì)算復(fù)雜度的同時(shí)保證了計(jì)算精度。

    1 離散縫洞網(wǎng)絡(luò)模型

    巖和裂縫系統(tǒng)視為滲流區(qū)域,其運(yùn)動(dòng)方程符合Darcy定律,即:

    式中:為流體黏度,mPa?s;為滲透率,μm2;為滲流速度,m/s;對于巖塊,=m,對于裂縫,=f;為多孔介質(zhì)中的平均壓力,Pa;為流體密度,kg/m3;為單位質(zhì)量力,m/s2。

    溶洞的空間尺度較大,其流動(dòng)模式為黏性流體的自由流動(dòng),為簡單起見本文僅考慮牛頓流體的低雷諾數(shù)運(yùn)動(dòng)。為了避免引入復(fù)雜的耦合邊界條件,降低自由流?滲流耦合問題的復(fù)雜性,本文采用Brinkman方程描述溶洞內(nèi)的自由流:

    式中:為流體真實(shí)速度,m/s;為溶洞內(nèi)滲透率(無充填時(shí)為無窮大),μm2;為溶洞內(nèi)孔隙度(無充填時(shí)為1);p為流體壓力,Pa。

    考慮等溫過程的單相穩(wěn)態(tài)流動(dòng),上述運(yùn)動(dòng)方程的統(tǒng)一形式結(jié)合質(zhì)量守恒方程得到離散縫洞網(wǎng)絡(luò)模型的控制微分方程組:

    式中:

    ;

    結(jié)合合適的邊界條件可得到完整的離散縫洞網(wǎng)絡(luò)單相穩(wěn)態(tài)滲流數(shù)學(xué)模型。

    2 混合模型原理

    由于實(shí)際的縫洞型油藏中不僅存在小尺度的裂縫和溶洞,還存在貫穿大部分油藏區(qū)域的大尺度高導(dǎo)流裂縫,對于這類實(shí)際油藏,單獨(dú)采用等效連續(xù)介質(zhì)模型和多重介質(zhì)模型往往得不到準(zhǔn)確的計(jì)算結(jié)果,離散縫洞網(wǎng)絡(luò)模型雖然可以準(zhǔn)確刻畫出縫洞型油藏中的流動(dòng)特征,但其計(jì)算量較大,僅適用于小范圍的油藏 模擬。

    本文對油藏中的小尺度裂縫和溶洞(裂縫長度和溶洞大小與所劃分的正交網(wǎng)格尺寸相差不大)通過超樣本技術(shù)和體積平均法結(jié)合離散縫洞網(wǎng)絡(luò)模型進(jìn)行等效處理,得到整個(gè)油藏的等效滲透率場;對于大尺度的裂縫(裂縫長度遠(yuǎn)大于正交網(wǎng)格尺寸),采用嵌入式離散裂縫模型將裂縫直接嵌入到前面通過等效計(jì)算得到的非均質(zhì)各向異性油藏中?;旌夏P偷恼w計(jì)算過程如圖1所示。

    2.1 小尺度裂縫及溶洞等效處理

    2.1.1 超樣本技術(shù)

    目前,基于流量等效計(jì)算等效滲透率的方法,都將各個(gè)網(wǎng)格單獨(dú)考慮而忽略網(wǎng)格周圍介質(zhì)對它的影響,這樣得到的等效滲透率往往不準(zhǔn)確。對此,本文對網(wǎng)格進(jìn)行超樣本技術(shù)處理,首先獲取每個(gè)網(wǎng)格的超樣本單元;然后在超樣本單元內(nèi)建立離散縫洞網(wǎng)絡(luò)模型并進(jìn)行有限元數(shù)值求解,從而獲取超樣本單元內(nèi)的壓力和速度場;最后,基于體積平均法計(jì)算目標(biāo)網(wǎng)格內(nèi)的速度、壓力梯度的體積平均值,求取目標(biāo)網(wǎng)格的等效滲透率張量。

    超樣本處理是將目標(biāo)原始網(wǎng)格及其周圍網(wǎng)格組成一個(gè)大的超樣本單元,如圖2所示,以二維系統(tǒng)為例:當(dāng)網(wǎng)格在角落時(shí),超樣本網(wǎng)格中包含4個(gè)網(wǎng)格;當(dāng)網(wǎng)格在邊界上時(shí),超樣本網(wǎng)格中有6個(gè)網(wǎng)格;當(dāng)網(wǎng)格在內(nèi)部時(shí),超樣本網(wǎng)格包含9個(gè)網(wǎng)格。

    圖1 混合模型計(jì)算過程示意圖

    (a) 含有1條裂縫的裂縫性介質(zhì)網(wǎng)格示意圖;(b) 網(wǎng)格在角落超樣本處理圖;(c) 網(wǎng)格在內(nèi)部超樣本處理圖;(d) 網(wǎng)格在邊界超樣本處理圖

    2.1.2 等效滲透率計(jì)算

    在超樣本單元內(nèi)結(jié)合封閉定壓邊界條件建立離散縫洞網(wǎng)絡(luò)模型,采用有限單元法進(jìn)行求解(具體的求解過程見參考文獻(xiàn)[8]),根據(jù)求得的壓力及速度場,運(yùn)用體積平均法求取目標(biāo)原始網(wǎng)格上的速度跟壓力梯度的體積平均值[18]:

    (5)

    (6)

    根據(jù)以上公式可以求得等效滲透率張量。由于定壓封閉邊界條件下得到的等效滲透率并不一定是對稱的,這時(shí)可以取取。

    2.2 大尺度裂縫顯示處理

    通過等效處理將只包含小尺度縫洞的油藏轉(zhuǎn)化為非均質(zhì)各向異性油藏。由于之前計(jì)算油藏等效滲透率場時(shí)已劃分了正交網(wǎng)格,采用嵌入式離散裂縫模型將大尺度裂縫直接嵌入到正交網(wǎng)格系統(tǒng)中,以裂縫與正交網(wǎng)格邊界相交的交點(diǎn)對裂縫進(jìn)行網(wǎng)格劃分。此時(shí)數(shù)學(xué)模型為:

    等效基質(zhì)系統(tǒng):

    (8)

    大尺度裂縫系統(tǒng):

    式中:m為基質(zhì)中滲流速度;m為基質(zhì)的等效滲透率張量;f為大尺度裂縫滲透率;m和f分別為基質(zhì)和裂縫的流體壓力;f和m分別為裂縫單元和基質(zhì)單元的體積;m和f分別為基質(zhì)和裂縫的源匯項(xiàng);為沿裂縫方向的局部坐標(biāo)系;mf為基巖與裂縫之間的竄流量;ff為相交裂縫單元之間的竄流量;

    ;

    借鑒常規(guī)油藏?cái)?shù)值模擬方法中生產(chǎn)井和網(wǎng)格單元間流量交換的計(jì)算格式,可以將基巖與裂縫之間的竄流量表達(dá)式定義為

    (10)

    與網(wǎng)格尺度相比,裂縫開度很小,且裂縫滲透率遠(yuǎn)大于基巖滲透率,因此可以認(rèn)為裂縫兩側(cè)的壓力連續(xù),根據(jù)達(dá)西定律,基巖與裂縫間的竄流量為[15]

    (11)

    本文采用Karimi-Fard[19]計(jì)算相交裂縫段間傳遞系數(shù)的方法計(jì)算裂縫與裂縫間竄流量[16]:

    3 混合模型求解

    由于求得的各基質(zhì)網(wǎng)格等效滲透率為全張量形式,因此,本文采用模擬有限差分方法對基質(zhì)部分的方程進(jìn)行離散。作為一種新型數(shù)值計(jì)算方法,模擬有限差分具有良好的局部守恒性和對復(fù)雜網(wǎng)格的適用性,并且能有效地處理全張量滲透率問題[20]。而大尺度裂縫滲透率考慮為標(biāo)量,因此,采用有限差分隱式離散。

    3.1 大尺度裂縫顯示處理

    基質(zhì)部分被正交網(wǎng)格所剖分,如圖3所示,取任意單元Ω進(jìn)行分析,其中Ω為其相鄰單元,交界面,為交界面A的面積加權(quán)法向量,為單位外法線向量。在單元中心點(diǎn)和邊界面中心點(diǎn)上分別定義單元平均壓力和邊界面平局壓力。

    圖3 模擬有限差分網(wǎng)格單元分析示意圖

    其中:mi為基質(zhì)網(wǎng)格的傳導(dǎo)矩陣,,為單元Ω的邊界面數(shù),,因此,模擬有限差分方法的關(guān)鍵問題是獲取矩陣mi。在此,假設(shè)壓力在基質(zhì)網(wǎng)格內(nèi)呈線性變化,即,顯然,由方程(7)可得:

    (14)

    (16)

    式中:

    (18)

    對于方程(8),直接在基質(zhì)網(wǎng)格單元上積分并運(yùn)用散度定理:

    考慮單元邊界面上的速度連續(xù)條件,結(jié)合方程(13)和(19),得到基質(zhì)部分模擬有限差分的數(shù)值計(jì)算格式:

    (20)

    其中:N為網(wǎng)格單元總數(shù);i=。

    3.2 裂縫部分有限差分隱式離散

    對于一維裂縫系統(tǒng)采用隱式差分,方程左右兩端同乘網(wǎng)格單元體積,則式(9)的差分方程為

    3.3 整體計(jì)算格式

    以存在2條交叉裂縫為例,將基質(zhì)部分和裂縫部分的數(shù)值計(jì)算格式組裝到一起,得到基于模擬有限差分的嵌入式離散裂縫模型計(jì)算方式:

    (23)

    其中:mfi=[mfi],表示第條裂縫與基巖竄流系數(shù)矩陣;ff=[ff],表示裂縫之間的竄流系數(shù)矩陣;fi和fi分別表示第條裂縫的有限差分傳導(dǎo)系數(shù)矩陣和裂縫單元壓力列陣。多條裂縫的計(jì)算格式與上類似。

    4 算例分析

    通過分析單相穩(wěn)定滲流模擬算例,并離散介質(zhì)模型結(jié)果對比,說明本文方法的正確性和優(yōu)越性。離散介質(zhì)模型(離散裂縫模型和離散縫洞網(wǎng)絡(luò)模型)作為一種宏觀精細(xì)模型,國內(nèi)外學(xué)者通過物理實(shí)驗(yàn)[16, 22?23]以及數(shù)值實(shí)驗(yàn)[8?9, 19]等方法對其準(zhǔn)確性進(jìn)行了驗(yàn)證,結(jié)果表明離散介質(zhì)模型能準(zhǔn)確刻畫縫洞型油藏中的流動(dòng)特征。因此,本文中將離散介質(zhì)模型作為對比驗(yàn)證的 標(biāo)準(zhǔn)。

    4.1 二維復(fù)雜裂縫性油藏

    100 m×100 m裂縫性油藏幾何模型如圖4(a)所示,其中存在2條長度分別為82.5 m和60.8 m的大尺度導(dǎo)流裂縫,裂縫開度和滲透率如表1所示。圖4(b)和圖4(c)所示分別為10×10正交網(wǎng)格和20×20正交網(wǎng)格的劃分結(jié)果。模型其他參數(shù)如表1所示。

    表1 裂縫性油藏模型基本參數(shù)

    不考慮大尺度裂縫時(shí)的裂縫性油藏如圖5(a)所示,通過劃分粗網(wǎng)格,建立離散裂縫模型并結(jié)合超樣本技術(shù)和體積平均法,得到該裂縫性油藏的等效滲透率場,見圖5(b)和圖5(c)。

    分別運(yùn)用離散裂縫模型(圖6(a))和本文提出的混合模型結(jié)合封閉定壓邊界條件(上、下邊界封閉,左、右邊界定壓,壓力值分別為1 MPa和2 MPa)對該裂縫性油藏進(jìn)行單相穩(wěn)定滲流模擬,其中混合模型分別10×10網(wǎng)格系統(tǒng)(圖6(b))和20×20網(wǎng)格系統(tǒng)(圖6(c))求解。圖7所示為=50 m直線上各方法求得的壓力分布曲線。通過式(24)可進(jìn)一步計(jì)算各方法相應(yīng)的誤差范數(shù)0。

    從圖7可以看出:2種方法的模型結(jié)果都能夠反映出裂縫性油藏基本流動(dòng)特征:裂縫內(nèi)壓力傳播快,大裂縫及周圍相連的小裂縫構(gòu)成明顯的等勢體,同時(shí)其他小裂縫也在一定程度上增加了油藏的非均質(zhì)性。采用混合模型時(shí),10×10的網(wǎng)格系統(tǒng)得到的壓力場分布與離散裂縫模型有一定偏差,誤差范數(shù)為2.2%,而采用20×20網(wǎng)格系統(tǒng)時(shí)混合模型的壓力場分布與離散裂縫模型的基本一致,誤差范數(shù)為0.85%。因此,本文提出的混合模型可以有效地模擬復(fù)雜裂縫性油藏中單相滲流壓力場分布,同時(shí),可以看出網(wǎng)格劃分對于模型精度有一定影響。

    (a) 裂縫性油藏幾何模型;(b) 10×10的正交網(wǎng)格系統(tǒng);(c) 20×20的正交網(wǎng)格系統(tǒng)

    圖4 裂縫性油藏幾何模型及網(wǎng)格劃分結(jié)果

    Fig. 4 Geometrical model and grids of fractured reservoir

    (a) 裂縫性油藏幾何模型;(b) 10×10的正交網(wǎng)格系統(tǒng);(c) 20×20的正交網(wǎng)格系統(tǒng)

    (a) 離散裂縫模型;(b) 混合模型(10×10);(c) 混合模型(20×20)

    1—離散裂縫模型(參考解);2—混合模型(20×20);3—混合模型(10×10)。

    4.2 二維復(fù)雜縫洞型油藏

    100 m×100 m縫洞型油藏幾何模型如圖8(a)所示,其中存在3條長度分別為70.7,99.0和70.7 m的大尺度導(dǎo)流裂縫,裂縫開度和滲透率如表1所示。模型其他參數(shù)如表1所示,其中溶洞部分視為無充填的自由流。圖8(b)所示為20×20正交網(wǎng)格劃分結(jié)果。圖8(c)所示為求得的不考慮大尺度裂縫時(shí)的等效滲透率場圖。

    分別運(yùn)用離散縫洞網(wǎng)絡(luò)模型(圖9(a))和本文提出的混合模型(圖9(b))結(jié)合封閉定壓邊界條件(上、下邊界封閉,左、右邊界定壓,壓力分別為1 MPa和2 MPa)對縫洞型油藏進(jìn)行單相穩(wěn)定滲流模擬。圖10所示為=47.5 m和=52.5 m直線上2種方法求得的壓力分布曲線。

    從圖10可以看出:模擬結(jié)果能夠反映出縫洞型油藏基本流動(dòng)特征。裂縫和溶洞內(nèi)壓力傳播快,大裂縫及周圍相連的小縫洞構(gòu)成明顯的等勢體,同時(shí)其他小縫洞也在一定程度上增加了油藏的非均質(zhì)性。采用混合模型20×20的網(wǎng)格系統(tǒng)得到的壓力場分布與離散縫洞網(wǎng)絡(luò)模型的結(jié)果基本一致,其誤差范數(shù)分別為1.2%和1.1%。因此,本文提出的混合模型同樣適用于二維復(fù)雜縫洞型油藏單相流模擬。并且本文模型對于縫洞參數(shù)易發(fā)生變化的油藏,當(dāng)局部小尺度縫洞結(jié)構(gòu)發(fā)生變化,只需要更新這部分的等效參數(shù),而當(dāng)大尺度裂縫結(jié)構(gòu)發(fā)生變化時(shí),只需要重新劃分該大尺度裂縫網(wǎng)格并更新裂縫與基質(zhì)之間的竄流信息,從而避免了離散縫洞網(wǎng)絡(luò)模型那樣的整體網(wǎng)格重新劃分,可以有效提高計(jì)算效率,圖9(c)所示為當(dāng)大尺度裂縫結(jié)構(gòu)發(fā)生變化后的壓力場計(jì)算結(jié)果。

    (a) 縫洞型油藏幾何模型;(b) 20×20的正交網(wǎng)格系統(tǒng);(c) 等效滲透率場圖

    (a) DFVN;(b) 混合模型;(c) 大尺度裂縫結(jié)構(gòu)發(fā)生改變

    圖10 不同位置壓力對比曲線

    4.3 三維復(fù)雜縫洞型油藏

    根據(jù)某一實(shí)際縫洞型油藏的裂縫、溶洞信息統(tǒng)計(jì)資料,生成相應(yīng)的三維復(fù)雜縫洞型油藏幾何模型,其長×寬×高為200 m×100 m×20 m,如圖11所示。其中裂縫在垂直方向上都貫穿地層,溶洞簡化成半徑為4 m的規(guī)則球體,溶洞內(nèi)不考慮充填,其流動(dòng)視為自由流。模型參數(shù)如表2所示。

    表2 三維縫洞型油藏模型基本參數(shù)

    分別運(yùn)用離散縫洞模型(圖12(a))和本文方法(圖12(b)),結(jié)合封閉定壓邊界條件(左邊界1和右邊界2定壓,壓力分別為1 MPa和2 MPa,其余邊界均為封閉邊界)對該復(fù)雜縫洞型油藏進(jìn)行單相穩(wěn)定滲流模擬。圖13所示為地層上表面=50 m直線上2種方法求得的壓力曲線,可以看出兩者結(jié)果吻合很好。因此,本文方法對于三維復(fù)雜縫洞型油藏同樣適用。

    (a) 縫洞型油藏幾何模型;(b) 等效滲透率場

    (a) DFVN;(b) 混合模型

    圖13 上表面y=50 m線上的壓力分布比較

    5 結(jié)論

    1) 將離散介質(zhì)模型與等效連續(xù)介質(zhì)模型相結(jié)合:對于小尺度裂縫和溶洞,基于離散縫洞網(wǎng)絡(luò)模型采用超樣本技術(shù)和體積平均方法對其進(jìn)行等效處理;對于大尺度裂縫,采用嵌入式離散裂縫模型進(jìn)行顯示處理;并基于模擬有限差分方法構(gòu)造了可以準(zhǔn)確處理全張量形式等效滲透率的數(shù)值計(jì)算格式。

    2) 本文方法對于復(fù)雜縫洞型油藏單相穩(wěn)定滲流模擬有較高的準(zhǔn)確性和靈活性。

    [1] 李陽. 塔河油田碳酸鹽巖縫洞型油藏開發(fā)理論及方法[J]. 石油學(xué)報(bào), 2013, 34(1): 115?121. LI Yang. The theory and method for development of carbonate fractured-cavity reservoirs in Tahe oilfield[J]. Acta Petrolei Sinica, 2013, 34(1): 115?121.

    [2] 李亞軍. 縫洞型介質(zhì)等效連續(xù)模型油水兩相流動(dòng)模擬理論研究[D]. 青島: 中國石油大學(xué)石油工程學(xué)院, 2011: 17?64. LI Yajun. Theoretical study of the equivalent continuum model for simulating oil-water two-phase flow in fractured-vuggy porous media[D]. Qingdao: China University of Petroleum. School of Petroleum Engineering, 2011: 17?64.

    [3] HUANG Zhaoqin, YAO Jun, WANG Yueying. An efficient numerical model for immiscible two-phase flow in fractured Karst reservoirs[J]. Communications in Computational Physics, 2013, 13(2): 540?558.

    [4] 吳玉樹, 葛家理. 三重介質(zhì)裂?隙油藏中的滲流問題[J]. 力學(xué)學(xué)報(bào), 1983, 19(1): 81?85. WU Yushu, GE Jiali. The transient flow in naturally fractured reservoirs with three-porosity systems[J]. Chinese Journal of Theoretical and Applied Mechanics, 1983, 19(1): 81?85.

    [5] ABDASSAH D, ERSHAGHI I. Triple-porosity systems for representing naturally fractured reservoirs[J]. SPE Formation Evaluation, 1986, 1(2): 113?127.

    [6] 楊堅(jiān), 姚軍, 王子勝. 三重介質(zhì)復(fù)合油藏壓力動(dòng)態(tài)特征研究[J]. 水動(dòng)力學(xué)研究與進(jìn)展(A輯), 2005, 20(4): 418?425. YANG Jian, YAO Jun, WANG Zisheng. Study of pressure-transient characteristic for triple-medium composite reservoirs[J]. Journal of Hydrodynamics (Series A), 2005, 20(4): 418?425.

    [7] WU Yushu, QIN Guan, EWING R, et al. A multiple-continuum approach for modeling multiphase flow in naturally fractured vuggy petroleum reservoirs[C]// SPE International oil & Gas Conference and Exhibition. Beijing: SPE, 2006: 1?12.

    [8] 姚軍, 黃朝琴, 王子勝, 等. 縫洞型油藏的離散縫洞網(wǎng)絡(luò)流動(dòng)數(shù)學(xué)模型[J]. 石油學(xué)報(bào), 2010, 31 (5): 15?20. YAO Jun, HUANG Zhaoqin, WANG Zisheng, et al. Mathematical model of fluid flow in fractured vuggy reservoirs based on discrete fracture-vug network[J]. Acta Petrolei Sinica, 2010, 31(5): 815?819.

    [9] YAO Jun, HUANG Zhaoqin, LI Yajun, et al. Discrete fracture-vug network model for modeling fluid flow in fractured vuggy porous media[C]// The CPS/SPE International Oil & Gas Conference and Exhibition. Beijing, China: SPE 2010: 1?14.

    [10] 李亞軍, 姚軍, 黃朝琴, 等. 基于Darcy-Stokes 耦合模型的縫洞型介質(zhì)等效滲透率分析[J]. 中國石油大學(xué)學(xué)報(bào)(自然科學(xué)版), 2011, 35(2): 91?95. LI Yajun, YAO Jun, HUANG Zhaoqin, et al. Estimating equivalent permeability of fractured-vuggy media based on coupled Darcy-Stokes model[J]. Journal of China University of Petroleum (Edition of Natural Science), 2011, 35(2): 91?95.

    [11] 李亞軍, 姚軍, 黃朝琴, 等. 裂縫性油藏等效滲透率張量計(jì)算及表征單元體積研究[J]. 水動(dòng)力學(xué)研究與進(jìn)展(A輯), 2010, 25(1): 1?7. LI Yajun, YAO Jun, HUANG Zhaoqin, et al. Calculation of equivalent permeability tensor and study on representative element volume for modeling fractured reservoirs[J]. Journal of Hydrodynamics (Series A), 2010, 25(1): 1?7.

    [12] 張世明, 嚴(yán)俠, 孫紅霞, 等. 基于離散裂縫模型的裂縫性介質(zhì)等效滲透率求解新方法[J]. 科學(xué)技術(shù)與工程, 2014, 14(16): 36?40. ZHANG Shiming, YAN Xia, SUN Hongxia, et al. An efficient calculation of equivalent permeability of fractured porous media based on discrete fracture model[J]. Science Technology and Engineering, 2014, 14(16): 36?40.

    [13] HUANG Zhaoqin, YAO Jun, LI Yajun, et al. Permeability analysis of fractured vuggy porous media based on homogenization theory[J]. Science China Technological Sciences, 2010, 53(3): 839?847.

    [14] HUANG Zhaoqin, YAO Jun, LI Yajun, et al. Numerical calculation of equivalent permeability tensor for fractured vuggy porous media based on homogenization theory[J]. Communications in Computational Physics, 2011, 9(1): 180?204.

    [15] LI L, LEE S H. Efficient field-scale simulation of black oil in a naturally fractured reservoir through discrete fracture networks and homogenized media[J]. SPE Reservoir Evaluation & Engineering, 2008, 11(04): 750?758.

    [16] 嚴(yán)俠, 黃朝琴, 姚軍, 等. 基于模擬有限差分的嵌入式離散裂縫數(shù)學(xué)模型[J]. 中國科學(xué):技術(shù)科學(xué), 2014, 44(12): 1333?1342. YAN Xia, HUANG Zhaoqin, YAO Jun, et al. The embeded discrete fracture model based on mimetic finite difference method[J]. Science China: Technological Sciences, 2014, 44(12): 1333?1342.

    [17] ZHOU Fangqi, SHI Anfeng, WANG Xiaohong. An efficient finite difference model for multiphase flow in fractured reservoirs[J]. Petroleum Exploration and Development, 2014, 41(2): 262?266.

    [18] DURLOFSKY L J. Upscaling and gridding of fine scale geological models for flow simulation[C]// 8th International Forum on Reservoir Simulation Iles Borromees, Stresa, Italy. 2005: 20?24.

    [19] KARIMI F M, DURLOFSKY L J, AZIZ K. An efficient discrete-fracture model applicable for general-purpose reservoir simulators[J]. SPE Journal, 2004, 9(2): 227?236.

    [20] LIPNIKOV K, MANZINI G, SHASHKOV M. Mimetic finite difference method[J]. Journal of Computational Physics, 2014, 257: 1163?1227.

    [21] BREZZI F, LIPNIKOV K, SIMONCINI V. A family of mimetic finite difference methods on polygonal and polyhedral meshes[J]. Mathematical Models and Methods in Applied Sciences, 2005, 15(10): 1533?1551.

    [22] 黃朝琴, 高博, 王月英, 等. 基于模擬有限差分法的離散裂縫模型兩相流動(dòng)模擬[J]. 中國石油大學(xué)學(xué)報(bào)(自然科學(xué)版), 2014, 38(6): 97?105. HUANG Zhaoqin, GAO Bo, WANG Yueying, et al. Two-phase flow simulation of discrete fracture model using a novel mimetic finite difference method[J]. Journal of China University of Petroleum (Edition of Natural Science), 2014, 38(6): 97?105.

    [23] YAN Xia, HUANG Zhaoqin, YAO Jun, et al. An efficient embedded discrete fracture model based on mimetic finite difference method[J]. Journal of Petroleum Science and Engineering, 2016, 145: 11?21.

    (編輯 趙俊)

    An efficient hybrid model for fractured-vuggy reservoir based on discrete fracture-vug network model

    YAN Xia1, HUANG Zhaoqin1, LI Yang2, YAO Jun1, FAN Dongyan1

    (1. School of Petroleum Engineering, China University of Petroleum, Qingdao 266580, China;2. Department of Oilfield Exploration & Development, SINOPEC, Beijing 100728, China)

    The existing numerical simulation methods are inaccurate or inefficient when modeling the flow and transport processes in fractured-vuggy reservoir with long fractures. A numerically efficient hybrid model was proposed to solve this problem, in which small fractures and vugs were modeled by equivalent permeability tensor, and long fractures were modeled by embedded discrete fracture model. Firstly, the orthogonal grid system was made according to the structure characteristics of fractures and vugs. Then, the discrete fracture-vug model (DFVN) was implemented in each orthogonal grid to calculate equivalent permeability tensor based on volume average method and oversampling technology, and the long fractures were embedded in equivalent media system. After that, an efficient numerical simulator was devised to solve the coupled system of long fracture and equivalent media based on mimetic finite difference method. At last, several single-phase steady numerical examples were shown to verify the validity and accuracy of the hybrid model. The results show that the proposed model is suitable for the fractured-vuggy reservoirs with dynamically changing fractures and vugs, because they can be simulated with an element effectively adapting itself to changes. Therefore the challenges associated with the remeshing for whole reservoirs, which has to be done in the DFVN, are bypassed entirely, and the computational costs can be reduced greatly.

    fractured-vuggy reservoir; discrete fracture-vug network model; oversampling technology; embedded discrete fracture model; mimetic finite difference

    10.11817/j.issn.1672-7207.2017.09.028

    TE319;O241.82

    A

    1672?7207(2017)09?2474?10

    2016?09?01;

    2016?12?22

    國家自然科學(xué)基金資助項(xiàng)目(51404292,51234007,51504277);山東省自然科學(xué)基金資助項(xiàng)目(ZR2014EEQ010);中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金資助項(xiàng)目(15CX05037A,14CX06091A);中國石油大學(xué)(華東)研究生創(chuàng)新工程資助項(xiàng)目(YCXJ2016015) (Projects(51404292, 51234007, 51504277) supported by the National Natural Science Foundation of China; Project(ZR2014EEQ010) supported by Shandong Provincial Natural Science Foundation, China; Projects(15CX05037A, 14CX06091A) supported by the Fundamental Research Funds for the Central Universities; Project(YCXJ2016015) supported by the Innovative Project of China University of Petroleum)

    黃朝琴, 博士,從事復(fù)雜介質(zhì)滲流理論及油藏?cái)?shù)值模擬;E-mail: huangzhqin@upc.edu.cn

    猜你喜歡
    縫洞溶洞滲透率
    別有洞天
    碳酸鹽巖縫洞儲(chǔ)集體分尺度量化表征
    出發(fā)吧,去溶洞
    哈拉哈塘奧陶系縫洞型成巖圈閉及其成因
    中煤階煤層氣井排采階段劃分及滲透率變化
    不同滲透率巖芯孔徑分布與可動(dòng)流體研究
    SAGD井微壓裂儲(chǔ)層滲透率變化規(guī)律研究
    神秘的溶洞
    幼兒100(2017年31期)2017-11-27 02:37:45
    高滲透率風(fēng)電并網(wǎng)對電力系統(tǒng)失步振蕩的影響
    縫洞型介質(zhì)結(jié)構(gòu)對非混相氣驅(qū)油采收率的影響
    桃红色精品国产亚洲av| 淫妇啪啪啪对白视频| 久久热在线av| 国内精品久久久久精免费| 黑人操中国人逼视频| 动漫黄色视频在线观看| 亚洲自拍偷在线| 又粗又爽又猛毛片免费看| 国产成人av教育| 国产欧美日韩一区二区精品| 国产精品av视频在线免费观看| 曰老女人黄片| 手机成人av网站| 国产成+人综合+亚洲专区| 可以在线观看毛片的网站| 不卡av一区二区三区| 亚洲最大成人中文| 韩国av一区二区三区四区| 欧美日韩瑟瑟在线播放| 欧美性猛交黑人性爽| av欧美777| 亚洲av日韩精品久久久久久密| 精品久久久久久久毛片微露脸| 我要搜黄色片| 亚洲av日韩精品久久久久久密| 悠悠久久av| 巨乳人妻的诱惑在线观看| 婷婷六月久久综合丁香| 免费观看人在逋| 别揉我奶头~嗯~啊~动态视频| 欧美成人免费av一区二区三区| 久久热在线av| 真人做人爱边吃奶动态| 国产亚洲av高清不卡| 国产成人系列免费观看| 久久香蕉国产精品| 亚洲五月天丁香| 女人被狂操c到高潮| 亚洲 欧美一区二区三区| 曰老女人黄片| 999精品在线视频| 中文字幕人成人乱码亚洲影| 亚洲精品粉嫩美女一区| 欧美xxxx黑人xx丫x性爽| av片东京热男人的天堂| 啦啦啦观看免费观看视频高清| 日日摸夜夜添夜夜添小说| 国产精品1区2区在线观看.| 日韩国内少妇激情av| 日韩国内少妇激情av| www国产在线视频色| 免费观看精品视频网站| 麻豆成人av在线观看| 麻豆久久精品国产亚洲av| 欧美一区二区精品小视频在线| 男女下面进入的视频免费午夜| 国产精品爽爽va在线观看网站| 十八禁人妻一区二区| 欧美一级毛片孕妇| 亚洲人成网站高清观看| 成人18禁在线播放| 日本三级黄在线观看| 亚洲无线观看免费| 国产精品一区二区精品视频观看| 国产一级毛片七仙女欲春2| 日本免费一区二区三区高清不卡| 黄色日韩在线| 黄色 视频免费看| 日本熟妇午夜| 亚洲成av人片在线播放无| 天堂网av新在线| 狂野欧美激情性xxxx| 99热精品在线国产| 老司机福利观看| 脱女人内裤的视频| 日韩精品中文字幕看吧| 亚洲人成网站在线播放欧美日韩| 亚洲自拍偷在线| 国产成人福利小说| or卡值多少钱| 1024香蕉在线观看| 免费观看人在逋| 女警被强在线播放| 在线观看日韩欧美| 精品99又大又爽又粗少妇毛片 | 高清毛片免费观看视频网站| 很黄的视频免费| 一级作爱视频免费观看| 国产 一区 欧美 日韩| 午夜久久久久精精品| 两性午夜刺激爽爽歪歪视频在线观看| 中文资源天堂在线| 精品一区二区三区视频在线观看免费| 神马国产精品三级电影在线观看| 深夜精品福利| 亚洲精品中文字幕一二三四区| 久久午夜亚洲精品久久| 老司机福利观看| 欧美日韩国产亚洲二区| 欧美一级a爱片免费观看看| 啦啦啦韩国在线观看视频| 又紧又爽又黄一区二区| 亚洲成人精品中文字幕电影| 美女午夜性视频免费| 久久中文字幕人妻熟女| 99久久无色码亚洲精品果冻| 91麻豆av在线| 亚洲第一欧美日韩一区二区三区| 国产 一区 欧美 日韩| 可以在线观看毛片的网站| 亚洲 欧美 日韩 在线 免费| 国产美女午夜福利| 久久国产精品人妻蜜桃| 亚洲一区高清亚洲精品| 麻豆一二三区av精品| 亚洲最大成人中文| 日韩人妻高清精品专区| 女人高潮潮喷娇喘18禁视频| 99在线人妻在线中文字幕| 高清毛片免费观看视频网站| 热99re8久久精品国产| 国内精品久久久久久久电影| 国内精品久久久久精免费| 香蕉久久夜色| 亚洲精品在线观看二区| 欧美日韩黄片免| 丁香欧美五月| 国产午夜精品久久久久久| 久久香蕉国产精品| 日韩av在线大香蕉| 久久性视频一级片| 成人三级做爰电影| 国产成人精品久久二区二区免费| 欧美日韩精品网址| 搡老岳熟女国产| 天天添夜夜摸| 一进一出抽搐动态| 99热精品在线国产| 欧美三级亚洲精品| 亚洲片人在线观看| 精品国产乱码久久久久久男人| 成年免费大片在线观看| 99热精品在线国产| 欧美在线黄色| tocl精华| 欧美日韩瑟瑟在线播放| 亚洲精品国产精品久久久不卡| 国产69精品久久久久777片 | 亚洲精华国产精华精| 国产aⅴ精品一区二区三区波| 1000部很黄的大片| 成人午夜高清在线视频| 男女那种视频在线观看| 在线视频色国产色| 成人av在线播放网站| 波多野结衣巨乳人妻| 亚洲欧美日韩无卡精品| 午夜精品在线福利| 欧美激情久久久久久爽电影| 国产午夜精品久久久久久| 欧美中文日本在线观看视频| 国产99白浆流出| 亚洲av日韩精品久久久久久密| 婷婷丁香在线五月| 久久久久久人人人人人| 午夜福利在线观看吧| 中文字幕久久专区| 男人舔女人下体高潮全视频| 18禁黄网站禁片午夜丰满| 在线十欧美十亚洲十日本专区| 黄片小视频在线播放| 一进一出抽搐动态| 亚洲熟女毛片儿| 老汉色av国产亚洲站长工具| 国产久久久一区二区三区| 我的老师免费观看完整版| 黄色视频,在线免费观看| 亚洲av电影不卡..在线观看| 日韩人妻高清精品专区| 麻豆av在线久日| 可以在线观看毛片的网站| 巨乳人妻的诱惑在线观看| 给我免费播放毛片高清在线观看| 国产 一区 欧美 日韩| 中文资源天堂在线| 又爽又黄无遮挡网站| 欧美大码av| 婷婷精品国产亚洲av| 97人妻精品一区二区三区麻豆| 欧美一区二区国产精品久久精品| 国产淫片久久久久久久久 | 国产精品综合久久久久久久免费| 日本精品一区二区三区蜜桃| 最近最新中文字幕大全免费视频| 在线观看免费午夜福利视频| 宅男免费午夜| 性色av乱码一区二区三区2| 亚洲人成伊人成综合网2020| 色综合婷婷激情| 日韩人妻高清精品专区| 国产午夜精品久久久久久| 国产淫片久久久久久久久 | 美女免费视频网站| 可以在线观看的亚洲视频| 色老头精品视频在线观看| 国产高清视频在线观看网站| 美女高潮的动态| 亚洲国产高清在线一区二区三| 欧美一区二区国产精品久久精品| 国产精品一区二区免费欧美| 在线观看一区二区三区| 亚洲成人久久爱视频| 一个人看的www免费观看视频| 久久精品综合一区二区三区| 最近最新中文字幕大全电影3| av天堂在线播放| 在线观看午夜福利视频| 在线国产一区二区在线| 长腿黑丝高跟| 国产人伦9x9x在线观看| 亚洲乱码一区二区免费版| 十八禁人妻一区二区| 亚洲无线观看免费| 国产精品日韩av在线免费观看| 日本三级黄在线观看| 欧美又色又爽又黄视频| 久久久色成人| 成人高潮视频无遮挡免费网站| 天堂动漫精品| 国内久久婷婷六月综合欲色啪| 波多野结衣高清无吗| 中国美女看黄片| 欧美日韩国产亚洲二区| 亚洲无线在线观看| 国产精品一区二区免费欧美| 久久99热这里只有精品18| 99久久精品国产亚洲精品| 久久香蕉精品热| 叶爱在线成人免费视频播放| 一级毛片精品| 久久久成人免费电影| 在线看三级毛片| 精品一区二区三区视频在线观看免费| 麻豆av在线久日| 成人永久免费在线观看视频| 国产精品日韩av在线免费观看| 天堂影院成人在线观看| 啦啦啦免费观看视频1| 欧美日韩一级在线毛片| 一区二区三区激情视频| 日韩av在线大香蕉| 免费看光身美女| 久久久久久久午夜电影| 级片在线观看| 欧美三级亚洲精品| 久久99热这里只有精品18| 身体一侧抽搐| 97超视频在线观看视频| 母亲3免费完整高清在线观看| 欧美极品一区二区三区四区| 欧美乱色亚洲激情| 一区福利在线观看| 叶爱在线成人免费视频播放| 久久伊人香网站| 夜夜看夜夜爽夜夜摸| 一二三四社区在线视频社区8| 亚洲成人久久爱视频| 亚洲无线观看免费| 国产免费av片在线观看野外av| 亚洲欧美日韩高清专用| 亚洲成av人片在线播放无| avwww免费| 国产欧美日韩一区二区精品| 欧美乱码精品一区二区三区| 国产一区在线观看成人免费| 岛国在线免费视频观看| 视频区欧美日本亚洲| 亚洲av熟女| 亚洲欧美一区二区三区黑人| 国产v大片淫在线免费观看| 亚洲国产欧美网| 国产精品av视频在线免费观看| 好男人在线观看高清免费视频| 亚洲av成人一区二区三| 精品国产亚洲在线| a级毛片a级免费在线| 亚洲精品乱码久久久v下载方式 | 亚洲午夜精品一区,二区,三区| 国产私拍福利视频在线观看| 在线观看日韩欧美| 日韩精品中文字幕看吧| 国产欧美日韩精品亚洲av| 亚洲一区高清亚洲精品| 在线观看免费视频日本深夜| 一个人免费在线观看的高清视频| 中出人妻视频一区二区| 国产午夜福利久久久久久| 久久中文字幕人妻熟女| 蜜桃久久精品国产亚洲av| 久久久久久人人人人人| 黄频高清免费视频| 国产人伦9x9x在线观看| 1000部很黄的大片| 亚洲狠狠婷婷综合久久图片| 欧美日韩亚洲国产一区二区在线观看| 精品国产亚洲在线| 制服丝袜大香蕉在线| 久久午夜综合久久蜜桃| 成人av一区二区三区在线看| 国产精品精品国产色婷婷| 亚洲第一电影网av| 精品电影一区二区在线| av在线天堂中文字幕| 亚洲avbb在线观看| 黄色 视频免费看| 一夜夜www| 午夜精品一区二区三区免费看| 国产精品影院久久| 精品国产乱码久久久久久男人| 久久久久精品国产欧美久久久| 最好的美女福利视频网| 久久亚洲真实| 国产综合懂色| 国产爱豆传媒在线观看| 亚洲九九香蕉| 欧洲精品卡2卡3卡4卡5卡区| 又紧又爽又黄一区二区| 成人特级av手机在线观看| 国产97色在线日韩免费| 免费看光身美女| 男女床上黄色一级片免费看| 美女高潮的动态| 国产99白浆流出| 黄色女人牲交| 97超级碰碰碰精品色视频在线观看| 黄色日韩在线| 热99在线观看视频| 老鸭窝网址在线观看| 成人一区二区视频在线观看| 亚洲成人精品中文字幕电影| 欧美在线一区亚洲| 亚洲在线自拍视频| 亚洲精品在线观看二区| 欧美成人性av电影在线观看| av黄色大香蕉| 精品国产超薄肉色丝袜足j| 久久精品综合一区二区三区| 久久精品国产99精品国产亚洲性色| 网址你懂的国产日韩在线| 欧美黑人欧美精品刺激| 久久这里只有精品19| 一二三四在线观看免费中文在| 九九热线精品视视频播放| 国产精品永久免费网站| 大型黄色视频在线免费观看| 国产欧美日韩精品亚洲av| 美女高潮的动态| 一二三四在线观看免费中文在| 99在线人妻在线中文字幕| 亚洲第一欧美日韩一区二区三区| av在线蜜桃| 亚洲精品中文字幕一二三四区| 欧美大码av| 久久国产乱子伦精品免费另类| 男女视频在线观看网站免费| 一级毛片高清免费大全| 日韩欧美国产一区二区入口| 男插女下体视频免费在线播放| 亚洲一区二区三区不卡视频| 91麻豆精品激情在线观看国产| 三级毛片av免费| 欧美日本视频| 婷婷六月久久综合丁香| 欧美丝袜亚洲另类 | 看片在线看免费视频| 国产三级在线视频| 美女cb高潮喷水在线观看 | 九九热线精品视视频播放| 午夜福利成人在线免费观看| 老汉色∧v一级毛片| 日本免费a在线| 欧美一级a爱片免费观看看| 少妇人妻一区二区三区视频| 男插女下体视频免费在线播放| 国产1区2区3区精品| 男女做爰动态图高潮gif福利片| 国产精品美女特级片免费视频播放器 | 18禁美女被吸乳视频| 精品久久久久久久久久久久久| 亚洲国产看品久久| 国内久久婷婷六月综合欲色啪| 国产伦精品一区二区三区视频9 | 亚洲国产日韩欧美精品在线观看 | 首页视频小说图片口味搜索| 日韩欧美精品v在线| e午夜精品久久久久久久| 国产三级在线视频| 特大巨黑吊av在线直播| 色噜噜av男人的天堂激情| 精品久久蜜臀av无| 国产精品香港三级国产av潘金莲| 国产av一区在线观看免费| 久久婷婷人人爽人人干人人爱| 亚洲国产精品久久男人天堂| 亚洲激情在线av| 久久久久精品国产欧美久久久| 免费av毛片视频| 欧美性猛交╳xxx乱大交人| 国产69精品久久久久777片 | 亚洲成av人片免费观看| 国产精品国产高清国产av| 国产精品九九99| 91老司机精品| 美女免费视频网站| xxxwww97欧美| 欧美一区二区精品小视频在线| 中国美女看黄片| 国产精品一区二区三区四区免费观看 | 亚洲专区国产一区二区| 男女之事视频高清在线观看| netflix在线观看网站| 九色成人免费人妻av| 中文字幕最新亚洲高清| 啦啦啦免费观看视频1| 桃色一区二区三区在线观看| 99久久久亚洲精品蜜臀av| 高清毛片免费观看视频网站| 亚洲av五月六月丁香网| 免费看美女性在线毛片视频| 国产1区2区3区精品| 中文字幕熟女人妻在线| 成人国产综合亚洲| 青草久久国产| 好男人电影高清在线观看| 国产视频内射| 国内揄拍国产精品人妻在线| 级片在线观看| 日韩三级视频一区二区三区| 看黄色毛片网站| 亚洲性夜色夜夜综合| 一个人看视频在线观看www免费 | 在线播放国产精品三级| 88av欧美| 国产1区2区3区精品| 在线免费观看的www视频| 欧洲精品卡2卡3卡4卡5卡区| 国产熟女xx| 国产精品一区二区精品视频观看| 国产精品av久久久久免费| 成人亚洲精品av一区二区| av在线蜜桃| 禁无遮挡网站| 又黄又爽又免费观看的视频| 成在线人永久免费视频| 午夜久久久久精精品| 欧美一区二区精品小视频在线| 在线看三级毛片| 亚洲国产看品久久| 国产精品久久电影中文字幕| 久久中文看片网| 99国产精品99久久久久| av视频在线观看入口| 亚洲人成网站在线播放欧美日韩| 一本一本综合久久| 村上凉子中文字幕在线| 久久久久精品国产欧美久久久| 国产私拍福利视频在线观看| 欧美极品一区二区三区四区| 神马国产精品三级电影在线观看| 亚洲av日韩精品久久久久久密| 99国产精品一区二区蜜桃av| 中文资源天堂在线| 久久这里只有精品中国| 久久这里只有精品19| 日韩免费av在线播放| 麻豆成人午夜福利视频| 老熟妇仑乱视频hdxx| 国产精品久久久久久人妻精品电影| 国产成人av激情在线播放| 午夜a级毛片| 热99re8久久精品国产| 欧美乱色亚洲激情| 国产一区二区激情短视频| 精品国内亚洲2022精品成人| 少妇裸体淫交视频免费看高清| 在线观看66精品国产| 亚洲国产色片| 亚洲国产精品sss在线观看| 狂野欧美激情性xxxx| 国产91精品成人一区二区三区| 精品久久久久久久人妻蜜臀av| 日日干狠狠操夜夜爽| 亚洲精品在线美女| 欧美3d第一页| 久久精品国产清高在天天线| 小蜜桃在线观看免费完整版高清| 日日干狠狠操夜夜爽| 国产99白浆流出| 色视频www国产| 国产精品久久久人人做人人爽| 淫妇啪啪啪对白视频| 成在线人永久免费视频| 黑人巨大精品欧美一区二区mp4| 亚洲成人久久爱视频| 日韩 欧美 亚洲 中文字幕| 夜夜躁狠狠躁天天躁| 精品国产乱码久久久久久男人| 99视频精品全部免费 在线 | 午夜福利高清视频| 欧美另类亚洲清纯唯美| 亚洲精品粉嫩美女一区| 久久久久国产精品人妻aⅴ院| 男人和女人高潮做爰伦理| 久久国产乱子伦精品免费另类| 最新美女视频免费是黄的| 听说在线观看完整版免费高清| 亚洲专区国产一区二区| 精品一区二区三区四区五区乱码| 日韩欧美一区二区三区在线观看| 波多野结衣巨乳人妻| 母亲3免费完整高清在线观看| 老熟妇仑乱视频hdxx| 亚洲av电影不卡..在线观看| 国产午夜精品久久久久久| 岛国在线观看网站| 麻豆成人午夜福利视频| 色噜噜av男人的天堂激情| 亚洲 欧美一区二区三区| av视频在线观看入口| 99久久久亚洲精品蜜臀av| 人妻久久中文字幕网| 99在线人妻在线中文字幕| 欧美日韩瑟瑟在线播放| 欧美日韩一级在线毛片| 国产高清激情床上av| 小蜜桃在线观看免费完整版高清| 精品久久久久久久末码| 又粗又爽又猛毛片免费看| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品国产精品久久久不卡| 曰老女人黄片| 美女 人体艺术 gogo| 男女做爰动态图高潮gif福利片| 波多野结衣巨乳人妻| 国产精品久久视频播放| 日本五十路高清| 久久人妻av系列| 黄频高清免费视频| 每晚都被弄得嗷嗷叫到高潮| 国产精品久久电影中文字幕| 九九热线精品视视频播放| 久久国产精品影院| 成年免费大片在线观看| 国产精品一区二区三区四区久久| 免费在线观看视频国产中文字幕亚洲| 禁无遮挡网站| 欧美在线黄色| 波多野结衣巨乳人妻| 久久久久国内视频| 日韩成人在线观看一区二区三区| 国内精品久久久久精免费| 国产成人精品无人区| 18禁观看日本| 色尼玛亚洲综合影院| 国产精品精品国产色婷婷| 精品一区二区三区视频在线观看免费| 欧美一级毛片孕妇| 久久久久久国产a免费观看| 在线观看66精品国产| 成年免费大片在线观看| 亚洲熟妇中文字幕五十中出| 免费无遮挡裸体视频| 亚洲人成伊人成综合网2020| 欧美黄色片欧美黄色片| 欧美日韩瑟瑟在线播放| 狂野欧美白嫩少妇大欣赏| 国产精品一及| 免费看光身美女| 国产美女午夜福利| 香蕉av资源在线| 一个人看视频在线观看www免费 | 一区二区三区国产精品乱码| 99riav亚洲国产免费| 在线视频色国产色| 丁香六月欧美| 18禁裸乳无遮挡免费网站照片| 久久国产乱子伦精品免费另类| 欧美最黄视频在线播放免费| 一个人观看的视频www高清免费观看 | 亚洲精品粉嫩美女一区| 一级黄色大片毛片| 国产97色在线日韩免费| 99在线人妻在线中文字幕| 亚洲专区国产一区二区| 性欧美人与动物交配| 国产一区二区在线av高清观看| 亚洲午夜精品一区,二区,三区| 免费在线观看亚洲国产| 波多野结衣高清作品| 国产精品99久久99久久久不卡| 最新在线观看一区二区三区| 亚洲熟妇中文字幕五十中出| 国语自产精品视频在线第100页| 老司机午夜十八禁免费视频| www.999成人在线观看| 亚洲av五月六月丁香网| 夜夜夜夜夜久久久久| 嫩草影院入口| 色在线成人网| 国产精品野战在线观看| 黄片大片在线免费观看| 99re在线观看精品视频| 99久久99久久久精品蜜桃| 亚洲欧美日韩高清专用| 男女午夜视频在线观看|