• <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在线免费观看 | 一级a爱片免费观看的视频| 国产精品久久久久久人妻精品电影| 看片在线看免费视频| 欧美黑人欧美精品刺激| 午夜成年电影在线免费观看| 亚洲美女黄片视频| 国产精品秋霞免费鲁丝片| 国产精品综合久久久久久久免费 | 丝袜在线中文字幕| 欧美精品一区二区免费开放| 成人特级黄色片久久久久久久| 一级片免费观看大全| 美国免费a级毛片| 嫩草影院精品99| 黄色毛片三级朝国网站| 19禁男女啪啪无遮挡网站| 国产成人av激情在线播放| 美女高潮喷水抽搐中文字幕| 久久午夜综合久久蜜桃| 久热这里只有精品99| 久久久久精品国产欧美久久久| 亚洲精品一二三| 丰满迷人的少妇在线观看| 国产成年人精品一区二区 | 欧美av亚洲av综合av国产av| 色婷婷av一区二区三区视频| 在线视频色国产色| 校园春色视频在线观看| 丰满的人妻完整版| 亚洲精品国产色婷婷电影| 88av欧美| 一级毛片高清免费大全| av在线天堂中文字幕 | 亚洲国产欧美日韩在线播放| 操出白浆在线播放| 成人av一区二区三区在线看| 18禁黄网站禁片午夜丰满| av天堂在线播放| 欧美一级毛片孕妇| 另类亚洲欧美激情| 亚洲欧美日韩高清在线视频| 国产精品日韩av在线免费观看 | 天天添夜夜摸| 80岁老熟妇乱子伦牲交| 神马国产精品三级电影在线观看 | 黄色丝袜av网址大全| 黄色视频不卡| 色婷婷久久久亚洲欧美| 一区二区三区精品91| 免费在线观看影片大全网站| e午夜精品久久久久久久| 久久精品国产亚洲av高清一级| 精品少妇一区二区三区视频日本电影| 国产av又大| 黄色片一级片一级黄色片| 99热国产这里只有精品6| cao死你这个sao货| 19禁男女啪啪无遮挡网站| 午夜a级毛片| 大型黄色视频在线免费观看| 人人妻,人人澡人人爽秒播| 黑人欧美特级aaaaaa片| 免费观看人在逋| 丝袜在线中文字幕| 亚洲avbb在线观看| 国产一区二区激情短视频| 国产精品1区2区在线观看.| 亚洲一卡2卡3卡4卡5卡精品中文| 天天躁狠狠躁夜夜躁狠狠躁| 久久热在线av| 免费观看人在逋| 999久久久精品免费观看国产| 免费不卡黄色视频| 他把我摸到了高潮在线观看| 看黄色毛片网站| 久久久久国内视频| 国产三级在线视频| 成人18禁高潮啪啪吃奶动态图| 一个人免费在线观看的高清视频| 亚洲色图av天堂| 国产欧美日韩综合在线一区二区| 国产不卡一卡二| xxxhd国产人妻xxx| 人成视频在线观看免费观看| 亚洲国产欧美日韩在线播放| 精品欧美一区二区三区在线| 亚洲熟妇中文字幕五十中出 | 欧美激情 高清一区二区三区| 看黄色毛片网站| 嫁个100分男人电影在线观看| 男女午夜视频在线观看| 人人澡人人妻人| 久久精品国产亚洲av高清一级| 12—13女人毛片做爰片一| 午夜久久久在线观看| 日韩视频一区二区在线观看| 黑人巨大精品欧美一区二区mp4| aaaaa片日本免费| 欧美黄色片欧美黄色片| www日本在线高清视频| 欧美日韩一级在线毛片| 国产免费av片在线观看野外av| 午夜91福利影院| 一区二区三区精品91| 男女下面进入的视频免费午夜 | 搡老熟女国产l中国老女人| 欧美乱妇无乱码| 日本 av在线| 午夜福利影视在线免费观看| 真人一进一出gif抽搐免费| 久久狼人影院| 午夜福利影视在线免费观看| 精品熟女少妇八av免费久了| 久久天堂一区二区三区四区| 老司机深夜福利视频在线观看| 别揉我奶头~嗯~啊~动态视频| 久热这里只有精品99| 欧美 亚洲 国产 日韩一| 久久久精品国产亚洲av高清涩受| 亚洲成人国产一区在线观看| 激情在线观看视频在线高清| 色婷婷av一区二区三区视频| 黄色视频不卡| 不卡一级毛片| 老司机福利观看| 91成年电影在线观看| 999久久久精品免费观看国产| 国产精品 欧美亚洲| 精品国产乱码久久久久久男人| 亚洲全国av大片| 淫妇啪啪啪对白视频| 99riav亚洲国产免费| 国产主播在线观看一区二区| 一本综合久久免费| 精品无人区乱码1区二区| 五月开心婷婷网| 欧美日韩黄片免| 搡老乐熟女国产| 久久久久久久久久久久大奶| 亚洲国产精品999在线| 国产精品一区二区精品视频观看| 无遮挡黄片免费观看| 国产亚洲精品第一综合不卡| 欧美一区二区精品小视频在线| 欧美日韩黄片免| 妹子高潮喷水视频| 动漫黄色视频在线观看| xxx96com| 国产欧美日韩一区二区精品| 长腿黑丝高跟| 国产高清视频在线播放一区| 日日摸夜夜添夜夜添小说| 久久午夜亚洲精品久久| 亚洲熟妇中文字幕五十中出 | 亚洲av第一区精品v没综合| 丁香欧美五月| 日韩免费高清中文字幕av| 国产无遮挡羞羞视频在线观看| 日本黄色视频三级网站网址| 一区二区三区激情视频| 丰满饥渴人妻一区二区三| 黄网站色视频无遮挡免费观看| 一边摸一边做爽爽视频免费| 91字幕亚洲| 国产亚洲精品一区二区www| 老司机午夜福利在线观看视频| 欧美在线一区亚洲| 男女做爰动态图高潮gif福利片 | 欧美日韩一级在线毛片| 色综合站精品国产| 久久人人精品亚洲av| 不卡av一区二区三区| 亚洲国产中文字幕在线视频| 亚洲国产欧美一区二区综合| 国产精品国产av在线观看| 动漫黄色视频在线观看| av国产精品久久久久影院| 国产国语露脸激情在线看| 亚洲狠狠婷婷综合久久图片| 亚洲精品在线观看二区| 老司机深夜福利视频在线观看| 人妻久久中文字幕网| 1024视频免费在线观看| 亚洲av美国av| 丝袜在线中文字幕| 91在线观看av| 精品久久蜜臀av无| 国产精品久久视频播放| 亚洲av美国av| 伊人久久大香线蕉亚洲五| 91九色精品人成在线观看| www国产在线视频色| av超薄肉色丝袜交足视频| 欧美午夜高清在线| 亚洲精品国产区一区二| 中文字幕高清在线视频| 久久草成人影院| a级毛片在线看网站| 亚洲av美国av| 99国产极品粉嫩在线观看| a级毛片黄视频| 淫妇啪啪啪对白视频| 老司机午夜十八禁免费视频| 精品免费久久久久久久清纯| 国产蜜桃级精品一区二区三区| 亚洲成人久久性| 又黄又粗又硬又大视频| 国产精品 欧美亚洲| av福利片在线| 俄罗斯特黄特色一大片| 黄色女人牲交| 亚洲全国av大片| 日韩精品中文字幕看吧| 黄色丝袜av网址大全| 女警被强在线播放| 亚洲国产精品一区二区三区在线| 另类亚洲欧美激情| 天天躁夜夜躁狠狠躁躁| 国产精品久久久久成人av| 多毛熟女@视频| 欧美色视频一区免费| 欧美日韩乱码在线| 国产精品国产av在线观看| 亚洲 国产 在线| 成人永久免费在线观看视频| 正在播放国产对白刺激| 黄色片一级片一级黄色片| av网站免费在线观看视频| 十八禁网站免费在线| 欧美日本亚洲视频在线播放| 变态另类成人亚洲欧美熟女 | 在线av久久热| 成人手机av| 香蕉丝袜av| a级片在线免费高清观看视频| 婷婷精品国产亚洲av在线| 国产精品爽爽va在线观看网站 | 99在线人妻在线中文字幕| 美女大奶头视频| 热re99久久精品国产66热6| 色播在线永久视频| 91大片在线观看| 久久精品国产99精品国产亚洲性色 | 免费观看人在逋| 国产蜜桃级精品一区二区三区| 国产精品影院久久| 久久香蕉精品热| 无限看片的www在线观看| 国产精品一区二区免费欧美| 久久草成人影院| 日韩有码中文字幕| 欧美在线一区亚洲| 成人免费观看视频高清| 在线天堂中文资源库| 欧美人与性动交α欧美精品济南到| 老司机福利观看| 欧美最黄视频在线播放免费 | 日韩精品青青久久久久久| 亚洲第一av免费看| 亚洲精品一区av在线观看| 12—13女人毛片做爰片一| 侵犯人妻中文字幕一二三四区| 女人爽到高潮嗷嗷叫在线视频| 免费在线观看日本一区| 国产精品永久免费网站| 狂野欧美激情性xxxx| 国产精品美女特级片免费视频播放器 | 日韩精品免费视频一区二区三区| 国产亚洲精品第一综合不卡| 午夜福利在线观看吧| 亚洲中文日韩欧美视频| 国产真人三级小视频在线观看| 日韩大码丰满熟妇| 欧美乱码精品一区二区三区| 婷婷六月久久综合丁香| 国产精品 国内视频| 午夜免费鲁丝| 国产深夜福利视频在线观看| 中文字幕人妻丝袜制服| 国产人伦9x9x在线观看| 国产高清视频在线播放一区| 亚洲精品中文字幕一二三四区| 亚洲欧美精品综合久久99| 757午夜福利合集在线观看| 久久精品国产清高在天天线| 久久精品亚洲av国产电影网| 中文字幕高清在线视频| 自线自在国产av| 国产av精品麻豆| 国产精品久久电影中文字幕| 老司机靠b影院| 免费高清在线观看日韩| av网站在线播放免费| 999久久久精品免费观看国产| 亚洲精品粉嫩美女一区| 国产激情久久老熟女| 纯流量卡能插随身wifi吗| 亚洲欧美一区二区三区黑人| 美女扒开内裤让男人捅视频| av电影中文网址| 亚洲精品久久成人aⅴ小说| 在线观看午夜福利视频| 99久久99久久久精品蜜桃| 国产精品久久久av美女十八| 成人黄色视频免费在线看| 黄色片一级片一级黄色片| av网站免费在线观看视频| 男女高潮啪啪啪动态图| 亚洲av成人不卡在线观看播放网| 在线看a的网站| 丰满迷人的少妇在线观看| 久久久久亚洲av毛片大全| 18美女黄网站色大片免费观看| 看黄色毛片网站| 好男人电影高清在线观看| 亚洲精品中文字幕一二三四区| 国产aⅴ精品一区二区三区波| 涩涩av久久男人的天堂| 青草久久国产| 欧美国产精品va在线观看不卡| 美女高潮喷水抽搐中文字幕| 日本wwww免费看| 亚洲精品一卡2卡三卡4卡5卡| 久久久久久久午夜电影 | 免费看a级黄色片| 一边摸一边抽搐一进一小说| 咕卡用的链子| 国产精品一区二区三区四区久久 | 欧美日本亚洲视频在线播放| 又大又爽又粗| 国产熟女xx| 99国产精品免费福利视频| www.精华液| 嫩草影视91久久| av网站免费在线观看视频| 成人精品一区二区免费| 国产精品免费一区二区三区在线| 色婷婷av一区二区三区视频| 亚洲专区国产一区二区| 极品教师在线免费播放| xxx96com| 久久人人精品亚洲av| a在线观看视频网站| 久久精品国产综合久久久| 国产欧美日韩一区二区三| 久久香蕉激情| 国产精品国产高清国产av| 成年人免费黄色播放视频| 少妇的丰满在线观看| 亚洲一区二区三区色噜噜 | 在线av久久热| 亚洲三区欧美一区| 最近最新中文字幕大全免费视频| 国产成人免费无遮挡视频| 久久精品国产清高在天天线| 欧美中文日本在线观看视频| 亚洲欧美精品综合一区二区三区| 99精品欧美一区二区三区四区| 人人妻人人澡人人看| 777久久人妻少妇嫩草av网站| 亚洲精品美女久久久久99蜜臀| 老司机午夜福利在线观看视频| 午夜两性在线视频| 久久中文字幕人妻熟女| 黄片播放在线免费| 麻豆久久精品国产亚洲av | netflix在线观看网站| 亚洲,欧美精品.| 熟女少妇亚洲综合色aaa.| 精品少妇一区二区三区视频日本电影| 亚洲中文av在线| 成人影院久久| 亚洲美女黄片视频| tocl精华| 国产精品综合久久久久久久免费 | 大码成人一级视频| 一区福利在线观看| 欧美大码av| 99热只有精品国产| 看免费av毛片| 亚洲av日韩精品久久久久久密| 国产成人免费无遮挡视频| 成人免费观看视频高清| 香蕉久久夜色| 男女下面插进去视频免费观看| 欧美日韩瑟瑟在线播放| 啪啪无遮挡十八禁网站| 亚洲精品成人av观看孕妇| 十八禁人妻一区二区| 亚洲成人免费电影在线观看| 一区在线观看完整版| 99久久久亚洲精品蜜臀av| 在线观看免费午夜福利视频| 国产精品影院久久| 中出人妻视频一区二区| 国产av一区二区精品久久| 黄色丝袜av网址大全| 欧美激情 高清一区二区三区| 精品一区二区三区av网在线观看| 欧美激情极品国产一区二区三区| 在线观看舔阴道视频| 日本免费a在线| 女性被躁到高潮视频| 免费人成视频x8x8入口观看| 亚洲熟女毛片儿| 两性夫妻黄色片| 免费在线观看影片大全网站| 精品熟女少妇八av免费久了| 国产男靠女视频免费网站| 夜夜夜夜夜久久久久| 国产精品久久久av美女十八| 色尼玛亚洲综合影院| 久久精品国产清高在天天线| 国产99久久九九免费精品| av有码第一页| 欧美人与性动交α欧美软件| a级毛片在线看网站| 亚洲成人免费电影在线观看| 欧美日韩福利视频一区二区| 精品人妻1区二区| 少妇的丰满在线观看| 午夜a级毛片| 天堂影院成人在线观看| 国产黄a三级三级三级人| 久久久久久久久中文| 欧美激情高清一区二区三区| 国产精品99久久99久久久不卡| 亚洲人成电影观看| 天天躁夜夜躁狠狠躁躁| 可以免费在线观看a视频的电影网站| 精品国产一区二区三区四区第35| 制服人妻中文乱码| 欧美日本中文国产一区发布| 黑人操中国人逼视频| 日本a在线网址| 久热爱精品视频在线9| 国产色视频综合| 亚洲专区中文字幕在线| 可以免费在线观看a视频的电影网站| 亚洲片人在线观看| 50天的宝宝边吃奶边哭怎么回事| 午夜福利在线免费观看网站| a级毛片黄视频| 午夜免费激情av| 亚洲成人久久性| 日本撒尿小便嘘嘘汇集6| 无遮挡黄片免费观看| 中文字幕色久视频| 18禁国产床啪视频网站| 99国产精品99久久久久| 一区二区日韩欧美中文字幕| 国产精品 欧美亚洲| 99精品欧美一区二区三区四区| 国产成人免费无遮挡视频| 欧美激情高清一区二区三区| 久久国产精品影院| 在线观看免费午夜福利视频| 丁香欧美五月| 精品第一国产精品| 成人亚洲精品一区在线观看| 日韩精品中文字幕看吧| 如日韩欧美国产精品一区二区三区| 97碰自拍视频| 亚洲avbb在线观看| 91精品三级在线观看| 男男h啪啪无遮挡| 免费女性裸体啪啪无遮挡网站| 亚洲专区国产一区二区| 免费在线观看黄色视频的| 亚洲成人免费av在线播放| 久久久国产欧美日韩av| 多毛熟女@视频| 老汉色av国产亚洲站长工具| 欧美乱色亚洲激情| 99riav亚洲国产免费| 三上悠亚av全集在线观看| 性欧美人与动物交配| 他把我摸到了高潮在线观看| 91精品国产国语对白视频| 一夜夜www| 亚洲av五月六月丁香网| 日本a在线网址| 国产熟女xx| 久久中文字幕一级| av在线天堂中文字幕 | 悠悠久久av| 久久精品人人爽人人爽视色| 搡老熟女国产l中国老女人| 午夜福利,免费看| 50天的宝宝边吃奶边哭怎么回事| 十八禁网站免费在线| 91老司机精品| www.999成人在线观看| www国产在线视频色| 日韩欧美一区视频在线观看| 欧洲精品卡2卡3卡4卡5卡区| 男女做爰动态图高潮gif福利片 | 美女 人体艺术 gogo| 亚洲伊人色综图| 午夜老司机福利片| 亚洲全国av大片| 免费不卡黄色视频| 啦啦啦免费观看视频1| 人妻久久中文字幕网| 精品人妻在线不人妻| 精品久久久久久成人av| 国产欧美日韩一区二区三| 久久国产精品人妻蜜桃| www.999成人在线观看| 久久国产亚洲av麻豆专区| 成年人免费黄色播放视频| 在线观看免费高清a一片| 国产一区二区三区在线臀色熟女 | 51午夜福利影视在线观看| 97人妻天天添夜夜摸| 超色免费av| 亚洲国产欧美日韩在线播放| 一进一出抽搐动态| 国产av在哪里看| 久久久久国产一级毛片高清牌| 国产真人三级小视频在线观看| 欧美中文综合在线视频| 国产精品免费一区二区三区在线| 50天的宝宝边吃奶边哭怎么回事| 国产一区二区激情短视频| 韩国av一区二区三区四区| 亚洲七黄色美女视频| 久热爱精品视频在线9| 91国产中文字幕| 18禁观看日本| 久久久久久亚洲精品国产蜜桃av| 欧美在线黄色| 国产高清国产精品国产三级| 国产精品二区激情视频| 在线国产一区二区在线| 国产伦人伦偷精品视频| 777久久人妻少妇嫩草av网站| 桃红色精品国产亚洲av| 老熟妇仑乱视频hdxx| aaaaa片日本免费| 午夜福利欧美成人| 亚洲欧洲精品一区二区精品久久久| 亚洲成av片中文字幕在线观看| 一级毛片高清免费大全| 成人18禁在线播放| 久久中文字幕一级| 99久久99久久久精品蜜桃| 亚洲精品美女久久久久99蜜臀| 久久国产精品人妻蜜桃| 日韩三级视频一区二区三区| 国产精品电影一区二区三区| 天天躁狠狠躁夜夜躁狠狠躁| 精品久久久久久成人av| 色老头精品视频在线观看| 亚洲一码二码三码区别大吗| 色婷婷久久久亚洲欧美| 91av网站免费观看| 国产亚洲精品一区二区www| 香蕉丝袜av| 日韩免费高清中文字幕av| 狠狠狠狠99中文字幕| 9热在线视频观看99| 久久天躁狠狠躁夜夜2o2o| 欧美 亚洲 国产 日韩一| 夜夜爽天天搞| 国产99白浆流出| 欧美色视频一区免费| 国产精品综合久久久久久久免费 | 国产免费现黄频在线看| 国产亚洲精品综合一区在线观看 | 久久精品国产99精品国产亚洲性色 | 国产亚洲精品久久久久5区| www日本在线高清视频| 夜夜爽天天搞| 国产亚洲精品久久久久久毛片| 欧美色视频一区免费| 99久久人妻综合| 大码成人一级视频| 黄色女人牲交| 老司机福利观看| 黄色女人牲交| 精品国产乱子伦一区二区三区| 女性生殖器流出的白浆| 美国免费a级毛片| 免费搜索国产男女视频| 精品少妇一区二区三区视频日本电影| 亚洲色图 男人天堂 中文字幕| 悠悠久久av| 香蕉国产在线看| 韩国精品一区二区三区| 少妇的丰满在线观看| 国产av一区二区精品久久| 人妻丰满熟妇av一区二区三区| 国产亚洲欧美在线一区二区| 中文字幕高清在线视频| 午夜亚洲福利在线播放| 亚洲欧美精品综合一区二区三区| 国产伦人伦偷精品视频| 97人妻天天添夜夜摸| 久久午夜亚洲精品久久| 国产深夜福利视频在线观看| 两人在一起打扑克的视频| 亚洲欧美日韩无卡精品| 在线播放国产精品三级| 欧美黑人欧美精品刺激| 国产日韩一区二区三区精品不卡| 欧美激情极品国产一区二区三区| av电影中文网址| 亚洲一区二区三区欧美精品| 99久久久亚洲精品蜜臀av| 桃色一区二区三区在线观看| 天天影视国产精品| 黑人欧美特级aaaaaa片| 国产精品久久久人人做人人爽|