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

    基于擠噴流效應(yīng)的Russell流體因子推廣及應(yīng)用

    2016-11-08 03:04:40張玉潔劉洪崔棟桂生馮玲麗
    地球物理學(xué)報 2016年10期
    關(guān)鍵詞:噴流泊松比裂隙

    張玉潔, 劉洪, 崔棟, 桂生, 馮玲麗

    1 中國科學(xué)院地質(zhì)與地球物理研究所,中國科學(xué)院油氣資源研究重點(diǎn)實驗室, 北京 100029 2 中國科學(xué)院大學(xué), 北京 100049 3 中國石油勘探開發(fā)研究院, 北京 100083 4 中國石油新疆油田分公司勘探開發(fā)研究院, 新疆克拉瑪依 834000

    ?

    基于擠噴流效應(yīng)的Russell流體因子推廣及應(yīng)用

    張玉潔1,2, 劉洪1*, 崔棟3, 桂生1,2, 馮玲麗4

    1 中國科學(xué)院地質(zhì)與地球物理研究所,中國科學(xué)院油氣資源研究重點(diǎn)實驗室, 北京100029 2 中國科學(xué)院大學(xué), 北京100049 3 中國石油勘探開發(fā)研究院, 北京100083 4 中國石油新疆油田分公司勘探開發(fā)研究院, 新疆克拉瑪依834000

    在油氣勘探與開發(fā)過程中,尋求能夠從地震資料中直接識別儲層油氣的流體指示因子至關(guān)重要.由于現(xiàn)存多數(shù)流體指示因子都是在Biot理論假設(shè)前提下建立起來的,因此在雙相孔隙介質(zhì)條件下不能有效地識別流體.為此,本文基于前人提出的雙孔介質(zhì)統(tǒng)一波動理論,考慮巖石裂縫間擠噴流效應(yīng),在經(jīng)典流體指示因子基礎(chǔ)上構(gòu)建了一種新的流體指示因子.該流體因子能夠較好地反映巖石內(nèi)孔隙流體的變化對波傳播產(chǎn)生的影響.在實際資料的應(yīng)用中,采用多種方法分析、對比該流體識別因子與以Biot理論為基礎(chǔ)的傳統(tǒng)流體因子的優(yōu)劣,理論分析與實際資料的驗證表明該流體因子對于儲層中油氣的檢測有較高的精度和靈敏度.

    流體因子; 流體識別; 擠噴流效應(yīng); Biot-Gassmann理論

    1 引言

    流體因子實質(zhì)上是一種用來表征地層含油氣特征的地震屬性,通常建立在巖石物理學(xué)基礎(chǔ)上,與孔隙流體在地震資料中表現(xiàn)出的異常特性有關(guān).而流體識別則是利用流體因子來識別儲層內(nèi)所含流體的類型,因此構(gòu)建流體因子是識別油氣藏流體類型的關(guān)鍵技術(shù)(印興耀等. 2014).流體因子最早由Smith和Gidlow(1987)提出,并利用疊前數(shù)據(jù)計算了流體因子.隨著進(jìn)一步的研究,流體因子被賦予了更多樣、更廣泛的意義.Goodway等(1999)提出了λ-μ-ρ技術(shù)來識別儲層流體,并將其作為流體因子.孔隙模量概念的提出,使得地震波縱橫波阻抗成為流體識別的工具.Russel等(2003)基于Biot-Gassmann理論首先提出了用于描述流體性質(zhì)的普適性方程,該方程比聲波阻抗對飽和孔隙流體更加敏感.Smith和Gidlow(2003)依據(jù)疊前AVO分析,通過計算流體因子和交會圖兩種角度來識別孔隙中的流體.為了增大含流體介質(zhì)的縱橫波速區(qū)別,寧忠華等(2006)基于前人的方法對Russell流體因子進(jìn)行改進(jìn),加入了縱橫波阻抗的高階項,構(gòu)建了高靈敏度流體因子.Quakenbush等(2006)結(jié)合泊松比和密度屬性提出了泊松阻抗的概念.王棟等(2009)則是利用縱橫波阻抗多次冪的組合形式,使高次冪放大需要突出差異的地方,同時使低次量綱壓制噪聲.

    以上方法都是基于Biot-Gassmann理論的孔隙結(jié)構(gòu)假設(shè).研究發(fā)現(xiàn),基于Biot-Gassmann理論的流體因子對于固結(jié)成熟的碎屑巖儲層孔隙流體具有較高的靈敏度.地下巖石介質(zhì)含有大量的孔隙和裂隙,而這些孔隙和裂隙對于油氣的賦存至關(guān)重要.孔隙與裂隙并存的情況下,裂隙的縱橫比要比孔隙小很多,在外力的作用下,裂縫中的流體會被擠壓到與其相連的孔隙中,產(chǎn)生擠噴流現(xiàn)象,因此描述巖石的彈性性質(zhì)會相當(dāng)復(fù)雜.Eshelby(1957)、Walsh(1965)等就針對裂隙對巖石彈性性質(zhì)的影響進(jìn)行了研究.在裂隙彈性力學(xué)方面,Budiansky和O′Connell(1976,1980),O′Connell(1984)指出裂隙介質(zhì)填充流體在外力作用下會產(chǎn)生“擠壓噴流”現(xiàn)象.Biot理論及O′Connel和Budiansky的理論(這里簡稱OB理論)均不能完整地描述孔隙、裂隙并存時介質(zhì)的彈性波動力學(xué)特征.為了尋求一個能夠描述雙相孔隙介質(zhì)及對彈性波動影響的統(tǒng)一理論,White(1975),Johnson(2001)以及Carcione等(2003)用孔隙流體的不均勻分布模擬了流體在裂隙、孔隙中的局部流動效應(yīng).Thomsen(1985)提出了一個使OB理論與Biot 理論相互兼容的理論,簡稱為Biot 相恰(Biot-Consistent)理論.Kuster和Toks?z(1974),Hudson(1981)探索了孔隙分布不均勻性對波傳播造成的影響,并在彈性波的散射理論的基礎(chǔ)上進(jìn)行進(jìn)一步的研究.Mavko和Nur(1979)用噴射流來描述流體在孔隙裂隙間局部流動產(chǎn)生的影響.Dvorkin和Nur(1993)在此基礎(chǔ)上提出BISQ(Biot/squirt)理論,用來解釋測量縱波時發(fā)現(xiàn)的高值頻散現(xiàn)象.孔隙介質(zhì)的研究中比較典型的一類是研究孔隙裂隙并存的介質(zhì),即雙重孔隙介質(zhì)理論和方法(Pride and Berryman, 2003; Berryman and Wang, 1995;Berryman,2006).巴晶(2010)也在Biot理論的基礎(chǔ)上研究了彈性波在雙相孔隙介質(zhì)中的波動效應(yīng).Gurevich等(2010)將Biot理論連通孔隙介質(zhì)條件推廣到不連通或半連通,并構(gòu)建了擠噴流頻散和衰減的模型.

    唐曉明(2011)對雙孔之間的擠噴流效應(yīng)進(jìn)行了具體分析,并通過擠噴流效應(yīng)與描述裂隙介質(zhì)的重要參數(shù)之間的聯(lián)系,得到了統(tǒng)一的波動理論,可以將孔隙與裂隙彈性波動力學(xué)結(jié)合起來.該理論不僅保持Biot理論基本特征及結(jié)構(gòu),還加入因裂隙造成的擠噴流項,提高了預(yù)測和模擬雙孔介質(zhì)中彈性波的精度,從而使得該理論的應(yīng)用前景更為廣泛.

    本文以唐曉明提出的含擠噴流項的巖石物理方程為基礎(chǔ),對比Russell流體因子,推導(dǎo)出新的流體因子.

    2 理論基礎(chǔ)

    本文所構(gòu)造的流體因子是基于Biot-Gassmann理論的Russell流體因子構(gòu)建與含孔隙、裂隙介質(zhì)彈性波動的統(tǒng)一理論,因此首先介紹一下這兩種理論.

    2.1Biot-Gassmann理論

    影響地震反射特性的主要因素是介質(zhì)速度和密度,而介質(zhì)含有孔隙并包含流體時,情況就變得復(fù)雜了.該問題首先被Biot(1941)提出,隨后Gassmann(1951)也對該問題進(jìn)行了研究,Kiref等(1990)指出Biot和Gassmann的結(jié)論是等價的.雖然之后又有很多理論被提出來,但是Biot-Gassmann理論是描述多孔介質(zhì)縱橫波速度穩(wěn)定性最好、使用最廣泛的理論.

    Biot利用拉梅常數(shù)的關(guān)系得出:

    (1)

    其中λsat和λdry分別為飽和巖石與巖石骨架的第一拉梅常數(shù),β為在孔隙內(nèi)流體壓力恒定的情況下流體與巖石體積變化之比,又稱Biot系數(shù),M為在確定巖石體積的情況下,流體擠入巖石骨架時所需的壓力,是一個特定的模量.同理,Gassmann從K、μ、ρ出發(fā),經(jīng)推導(dǎo)得出以下的關(guān)系:

    (2)

    其中Ksat為飽和巖石的體積模量,Kdry為干燥巖石的體積模量.由(1)式和(2)式可得

    (3)

    μsat、μdry分別為飽和巖石和干燥巖石的剪切模量.這表明剪切模量是不受孔隙流體影響的.

    在低頻假設(shè)前提下,Gassmann指出:

    (4)

    (5)

    其中,φ為孔隙度,Km、Kf分別為巖石骨架和流體的體積彈性模量.將β和M代入公式(2)可得Gaussman巖石物理方程具體的形式:

    (6)

    2.2含孔隙、裂隙介質(zhì)彈性波動的統(tǒng)一理論

    唐曉明在Gassmann的基礎(chǔ)上提出新的巖石物理方程(公式(7)),該方程考慮了雙孔介質(zhì)的情況和彈性波動條件下裂隙之間流體的擠噴流效應(yīng),因此它可以有效地模擬與預(yù)測裂隙間流體的瞬間流動對波傳播的影響. 不僅如此, 該方程還保持了Biot理論的基本特征及其結(jié)構(gòu).與Gassmann方程(6)相比,式(7)在分式的分母中多了一項表征擠噴流貢獻(xiàn)的項S(ω).

    (7)

    (8)

    (9)

    其中f為頻率變化因子,γ代表裂隙縱向尺度與橫向尺度之比,η為流體的黏滯系數(shù),ε為裂隙密度,Ks為基質(zhì)的體積模量.J0是零階貝塞爾函數(shù),J1是一階貝塞爾函數(shù).ω是頻率.式(8)中的K等效于式(7)中的K,這種情況下(8)式中關(guān)于K的表達(dá)式轉(zhuǎn)變?yōu)殡[函數(shù)關(guān)系,導(dǎo)致計算復(fù)雜化.在一級近似條件下,可以由不考慮“擠噴流”效應(yīng)影響的(6)式計算得到(8)式中的K值,而μ的取值為式(3)中的干燥巖石剪切模量.類似的,υ可以是干燥介質(zhì)的泊松比.需要說明的是,式(8)的計算結(jié)果主要由ε和γ二者共同控制的,它們是裂隙介質(zhì)描述時最重要的參數(shù).在實際資料的應(yīng)用中這些參數(shù)可以由測井資料計算得到.

    3 構(gòu)建流體因子

    Russell流體因子的構(gòu)建過程是本文流體因子構(gòu)建的基礎(chǔ),本文將從Russell流體因子的推導(dǎo)過程出發(fā)逐步推出新的流體因子.

    均勻彈性無孔介質(zhì)中縱橫波速度的基本方程為:

    (10)

    這里保留用作區(qū)分干燥與飽和兩種不同情況的β2M項,以強(qiáng)調(diào)其并不依賴方程(1)和(2)中等號左側(cè)第一項.重新整理含飽和流體巖石的P波速度公式,在Biot理論框架下縱橫波速度公式寫為:

    (11)

    其中ρsat為巖石含有飽和流體時的密度.

    由于這里考慮的流體為理想液體不存在切應(yīng)力,對于飽和情況下的橫波速度方程僅僅改變了密度項.因此這里僅對縱波速度進(jìn)行討論.

    將公式(11)中受流體影響和不受流體影響的項分別用s和f表示,那么縱波速度表達(dá)式被重新寫為

    (12)

    提取式(12)中的f并整理,得Russell流體因子表達(dá)式:

    (13)

    其中

    (14)

    ZS=ρVS,

    (15)

    (16)

    Russell流體因子是基于Biot理論的流體因子,不考慮裂隙而僅考慮孔隙,且孔隙分布均勻,另外孔隙連通的假設(shè)也保證了波傳播引起的孔隙流體流動是充分均衡的,即忽略了“擠噴流”效應(yīng)對波傳播的影響,使得彈性波傳播經(jīng)過裂隙時產(chǎn)生的相當(dāng)可觀的衰減和頻散現(xiàn)象不能夠得到更加準(zhǔn)確地描述.為更加準(zhǔn)確地描述孔隙、孔隙和裂縫并存的情況對彈性波傳播的影響和含油氣性,構(gòu)造能更準(zhǔn)確地指示流體賦存情況的流體因子,本文將唐曉明提出的擠噴流效應(yīng)與流體因子結(jié)合進(jìn)來.

    將(7)式化為(2)式的形式可得

    (17)

    對比(11)式縱波速度可寫為

    (18)

    此時

    (19)

    (20)

    (21)

    其中S(ω)為K的隱函數(shù)形式,可以對工區(qū)測井?dāng)?shù)據(jù)進(jìn)行統(tǒng)計得到,而M可以用疊前反演方法所得到的縱橫波速度與密度得到.將ρf作為新推導(dǎo)出的流體因子,記為NF.

    該流體因子不僅保留了常規(guī)流體因子原有的特點(diǎn)及結(jié)構(gòu),還反映了擠噴流效應(yīng)對波傳播的影響.擠噴流項是由裂隙介質(zhì)的兩個最重要的參數(shù),裂隙密度和縱橫比所控制,因此擠噴流項能有效地體現(xiàn)裂縫對彈性波傳播的影響.因此加入該項的新流體因子能夠較好地模擬、預(yù)測介質(zhì)含裂隙情況下彈性波在雙孔介質(zhì)中傳播情況,提高流體預(yù)測的精度.新流體因子公式(20)中所需要的裂縫的密度及縱橫比、縱波速度、橫波速度與密度等參數(shù)均由測井資料與反演中得來.

    圖1 基準(zhǔn)井裂縫情況Fig.1 Fractures in stratigraphic well

    4 數(shù)據(jù)試算

    4.1儲層巖石物理敏感參數(shù)分析

    本文選取了某油田油頁巖工區(qū)進(jìn)行分析.主要巖性為泥巖、粉砂巖、油泥巖、泥質(zhì)白云巖、泥質(zhì)砂巖、白云巖、灰質(zhì)泥巖、灰質(zhì)油頁巖、泥質(zhì)灰?guī)r、油頁巖.目的層沉積環(huán)境為半深湖-深湖相.工區(qū)裂隙發(fā)育豐富(如圖1).目標(biāo)層地勢西南高,東北低.在中部靠東北側(cè)有一塊高速體,已知為巖漿侵入體,且該巖漿侵入體區(qū)域內(nèi)不含油氣.儲層巖石物理敏感參數(shù)分析選用所研究工區(qū)內(nèi)有橫波資料的井作為基準(zhǔn)井,對其進(jìn)行了多方面流體敏感度分析.

    首先利用基準(zhǔn)井各測井參數(shù)求出該井的NF,IF,泊松比,縱橫波速比,HSFIF曲線(如圖2).選取這些曲線上典型的含油段、含水段、干井段,然后再將這些流體指示因子的值與泊松比兩兩做成交會圖進(jìn)行比較.

    圖3為各流體因子與泊松比交會圖,其中只有NF流體因子指示的參數(shù)與泊松比的交會圖對含油層、干層和水層區(qū)分較為明顯.我們還分析了Russell流體因子(這里簡稱FI流體因子),該流體因子能夠區(qū)分油層和水層,但無法區(qū)分油層和干層.縱橫波速度比與泊松比的交會圖、HSFIF流體因子與泊松比交會圖都無法將油層區(qū)分出來.

    圖2 基準(zhǔn)井流體因子Fig.2 Fluid factors of stratigraphic well

    圖3 巖石物理參數(shù)分析交會圖藍(lán)色為油層,紅色為干層,黃色為含水層. Fig.3 Cross plots of rock physical parameters Blue points represents oil, red points represents dry rock, yellow points represents water.

    本文計算了該區(qū)有背景巖性影響情況和無背景巖性影響情況下多個流體因子的流體指示系數(shù),并對結(jié)果進(jìn)行綜合分析.無背景巖性影響(Dillon, 2003)情況下(如圖4),Russell流體因子f,泊松比,VP/VS的流體指示系數(shù)較高.再結(jié)合交會圖(圖3),以上流體因子可以在目標(biāo)區(qū)域?qū)⒂退畢^(qū)分開,但是卻不能將含油儲層與背景巖層明確區(qū)分.

    有背景巖性影響情況下(如圖5)NF,泊松比,HSFIF流體指示系數(shù)較高,再與交會圖進(jìn)行綜合分析,可以將油與水和背景巖層區(qū)分開.因此,NF流體因子受工區(qū)背景巖性影響較小,對本工區(qū)流體的指示有著較好的效果.

    4.2有利儲層分析及油氣預(yù)測

    圖6為目標(biāo)層構(gòu)造圖.該區(qū)為油頁巖區(qū)塊,斷層裂隙發(fā)育.在圖7的NF流體因子預(yù)測中可以看到含油范圍與井產(chǎn)油情況基本吻合.從儲層疊前屬性切片(圖7)上分析,NF值較低的位置為良好的油氣儲層位置,與實際鉆井產(chǎn)油情況基本吻合,黃色區(qū)域的分布范圍基本包含了所有的產(chǎn)油井,干井基本都在藍(lán)色區(qū)域內(nèi).這驗證了NF流體因子的有效性,與產(chǎn)油井情況的吻合度為85.7%.對比泊松比切片(圖8)和HSFIF流體因子切片(圖9),泊松比切片受地勢構(gòu)造影響嚴(yán)重.西南方向高地勢位置流體指示過多,而東北部低地勢位置流體指示過少.HSFIF流體因子受背景巖性影響嚴(yán)重,在巖漿侵入體位置有錯誤的油氣指示.

    圖4 無背景巖性影響情況下的流體指示系數(shù)Fig.4 Coefficient of fluid indicator without background lithology influence

    圖5 有背景巖性影響情況下的流體指示系數(shù)Fig.5 Coefficient of fluid indicator with background lithology influence

    圖6 目標(biāo)層構(gòu)造圖Fig.6 Structure of target layer

    圖7 NF流體因子目標(biāo)沿層切片F(xiàn)ig.7 NF slice of target layer

    圖8 泊松比目標(biāo)沿層切片F(xiàn)ig.8 Poisson ratio slice of target layer

    圖9 HSFIF流體因子目標(biāo)沿層切片F(xiàn)ig.9 HSFIF slice of target layer

    將構(gòu)造圖(圖6)覆蓋于流體因子切片之上,也可以看到斷裂發(fā)育的部位同時發(fā)育良好儲層.說明反演結(jié)果不但能夠反映儲層地球物理特征,還對工區(qū)內(nèi)的構(gòu)造特征進(jìn)行了一定程度的反映.

    5 結(jié)論與認(rèn)識

    流體因子能夠較好地指示儲層流體的分布情況,對于儲層的預(yù)測起著至關(guān)重要的作用.通過反演得到的速度密度信息,可以方便地計算多種流體因子.本文結(jié)合擠噴流效應(yīng)對流體指示的影響,提出了一種新的流體指示因子,即NF流體因子.交會分析表明NF流體因子對裂縫、裂隙油氣儲集體中油氣指示效果良好.實際工區(qū)內(nèi)含油氣范圍也表明該流體因子對于儲層的展布情況有較為準(zhǔn)確的預(yù)測.通過理論分析和實際資料驗證,NF流體因子能夠?qū)搸r裂縫發(fā)育區(qū)域的流體進(jìn)行較準(zhǔn)確的預(yù)測,有較好的應(yīng)用前景,為裂縫發(fā)育的泥頁巖探區(qū)的勘探開發(fā)提供可靠的依據(jù).

    Ba J. 2010. Wave propagation theory in double-porosity medium and experimental analysis on seismic responses.ScientiaSinica(Physica,Mechanica&Astronomica) (in Chinese), 40(11): 1398-1409.

    Berryman J G, Wang H F. 1995. The elastic coefficients of double-porosity models for fluid transport in jointed rock.JournalofGeophysicalResearch:SolidEarth(1978—2012), 100(B12): 24611-24627.

    Berryman J G. 2006. Effective medium theories for multicomponent poroelastic composites.JournalofEngineeringMechanics, 132(5): 519-531.

    Biot M A. 1941. General theory of three-dimensional consolidation.JournalofAppliedPhysics, 12(2): 155-164.

    Budiansky B, O′Connell R J. 1976. Elastic moduli of a cracked solid.InternationalJournalofSolidsandStructures, 12(2): 81-97.

    Budiansky B, O′Connell R. 1980. Bulk dissipation in heterogeneous media.SolidEarthGeophysicsandGeotechnology, 42: 1-10.

    Carcione J M, Helle H B, Pham N H. 2003. White's model for wave propagation in partially saturated rocks: Comparison with poroelastic numerical experiments.Geophysics, 68(4): 1389-1398. Dillon L, Schwedersky G, Vásquez G, et al. 2003. A multiscale DHI elastic attributes evaluation.TheLeadingEdge, 22(10): 1024-1029.

    Dvorkin J, Nur A. 1993. Dynamic poroelasticity: A unified model with the squirt and the Biot mechanisms.Geophysics, 58(4): 524-533.

    Eshelby J D. 1957. The determination of the elastic field of an ellipsoidal inclusion, and related problems. Proceedings of the Royal Society of London.SeriesA:MathematicalandPhysicalSciences, 241(1226): 376-396.

    Gassmann F. 1951. Elastic waves through a packing of spheres.Geophysics, 16(4): 673-685.

    Goodway B, Chen T W, Downton J. 1999. Improved AVO fluid detection and lithology discrimination using Lamé petrophysical parameters; “λρ”, “μρ”, & “λ/μfluid stack”, from P and S inversions.∥ 1997 SEG Annual Meeting. Dallas, Texas: SEG.

    Gurevich B, Makarynska D, de Paula O B, et al. 2010. A simple model for squirt-flow dispersion and attenuation in fluid-saturated granular rocks.Geophysics, 75(6): N109-N120.

    Hudson J A. 1981. Wave speeds and attenuation of elastic waves in material containing cracks.GeophysicalJournalInternational, 64(1): 133-150.

    Johnson D L. 2001. Theory of frequency dependent acoustics in patchy-saturated porous media.TheJournaloftheAcousticalSocietyofAmerica, 110(2): 682-694.

    Krief M, Garat J, Stellingwerff, et al. 1990. A petrophysical interpretation using the velocities of P and S waves (full-waveform sonic).LogAnalyst, 31: 355-369.

    Kuster G T, Toks?z M N. 1974. Velocity and attenuation of seismic waves in two-phase media: Part I. Theoretical formulations.Geophysics, 39(5): 587-606. Mavko G M, Nur A. 1979. Wave attenuation in partially saturated rocks.Geophysics, 44(2): 161-178.

    Ning Z H, He Z H, Huang D J. 2006. High sensitive fluid identification based on seismic data.GeophysicalProspectingforPetroleum(in Chinese), 45(3): 239-241.

    O′Connell R J. 1984. A viscoelastic model of anelasticity of fluid saturated porous rocks.∥ Physics and Chemistry of Porous Media. Schlumberger: AIP Publishing, 166-175.

    Pride S R, Berryman J G. 2003. Linear dynamics of double-porosity dual-permeability materials. I. Governing equations and acoustic attenuation.PhysicalReviewE, 68(3): 036603.

    Quakenbush M, Shang B, Tuttle C. 2006. Poisson impedance.TheLeadingEdge, 25(2): 128-138.

    Russell B H, Hedlin K, Hilterman F J, et al. 2003. Fluid-property discrimination with AVO: A Biot-Gassmann perspective.Geophysics, 68(1): 29-39. Smith G C, Gidlow P M. 1987. Weighted stacking for rock property estimation and detection of GAS.GeophysicalProspecting, 35(9): 993-1014.

    Smith G C, Gidlow M. 2003. The fluid factor angle and the crossplot angle. ∥ 2003 SEG Annual Meeting. Dallas, Texas: SEG, 185-188.

    Tang X M. 2011. A unified theory for elastic wave propagation through porous media containing cracks—An extension of Biot′s poroelastic wave theory.ScientiaSinica(EarthScience) (in Chinese), 41(6): 784-795.

    Thomsen L. 1985. Biot-consistent elastic moduli of porous rocks: Low-frequency limit.Geophysics, 50(12): 2797-2807

    Walsh J B. 1965. The effect of cracks on the compressibility of rock.JournalofGeophysicalResearch, 70(2): 381-389.

    Wang D, He Z H, Huang D J. 2009. Construction of a new fluid identification factor and analysis on its application.GeophysicalProspectingforPetroleum(in Chinese), 48(2): 141-145.

    White J E. 1975. Computed seismic speeds and attenuation in rocks with partial gas saturation.Geophysics, 40(2): 224-232.

    Yin X Y, Cao D P, Wang B L, et al. 2014. Research progress of fluid discrimination with prestack seismic inversion.OilGeophysicalProspecting(in Chinese), 49(1): 22-34.

    附中文參考文獻(xiàn)

    巴晶. 2010. 雙重孔隙介質(zhì)波傳播理論與地震響應(yīng)實驗分析. 中國科學(xué): 物理學(xué) 力學(xué) 天文學(xué), 40(11): 1398-1409.

    寧忠華, 賀振華, 黃德濟(jì). 2006. 基于地震資料的高靈敏度流體識別因子. 石油物探, 45(3): 239-241.

    唐曉明. 2011. 含孔隙、裂隙介質(zhì)彈性波動的統(tǒng)一理論——Biot理論的推廣. 中國科學(xué): 地球科學(xué), 41(6): 784-795.

    王棟, 賀振華, 黃德濟(jì). 2009. 新流體識別因子的構(gòu)建與應(yīng)用分析. 石油物探, 48(2): 141-145.

    印興耀, 曹丹平, 王保麗等. 2014. 基于疊前地震反演的流體識別方法研究進(jìn)展. 石油地球物理勘探, 49(1): 22-34.

    (本文編輯何燕)

    Construction and application of the Russell fluid factor with squirt flow effect

    ZHANG Yu-Jie1,2, LIU Hong1*, CUI Dong3, GUI Sheng1,2, FENG Ling-Li4

    1KeyLaboratoryofPetroleumResourcesResearch,InstituteofGeologyandGeophysics,ChineseAcademyofSciences,Beijing100029,China2UniversityofChineseAcademyofSciences,Beijing100049,China3ResearchInstituteofPetroleumExplorationandDevelopment,CNPC,Beijing100083,China4InstituteofExplorationandDevelopment,XinjiangOilfieldCompany,CNPC,XinjiangKaramay834000,China

    In the exploration process of oil and gas, it is important that using a fluid factor to identify gas and oil from seismic data directly. In this work, on the basis of the classical fluid indication factor analysis, we built a new indication factor named as NF fluid factor. Adding the squirt fluid effect between the cracks of rocks, the fluid factor indicates the fluid in the rock pore more possible. It plays a vital role in predicting reservoirs. We compared it with other fluid factors through many methods. We found that it not only indicates the oil but also indicates the dispersion of fractures without the influence of lithology. In oil shale reservoirs, fractures are the main accumulate space, which means that it is the fracture where oil is present. We also compared the production of wells with NF fluid factor time slice and found that the percent of consistence is 85%. In conclusion, the NF fluid factor can indicate the oil from oil shale and ignore the influence of lithology. It is a perfect fluid factor to predict oil in shale. It has been used in real exploration and works well.

    Fluid factor; Fluid identification; Squirt flow effect; Biot-Gassmann theory

    10.6038/cjg20161033.

    中石油新技術(shù)新方法項目“GPU/CPU加速彈性介質(zhì)全波形反演研究”、中國石油天然氣集團(tuán)公司科學(xué)研究與技術(shù)開發(fā)項目“彈性波地震成像技術(shù)合作研發(fā)與應(yīng)用”,國家高技術(shù)研究發(fā)展計劃(863計劃)(2012AA061202),國家重大專項(2016ZX05003003)聯(lián)合資助.作者簡介張玉潔,女,1988年生,博士研究生,主要從事地震波成像和油儲地球物理研究. E-mail: zhangyujie@mail.iggcas.ac.cn

    劉洪,男,1959年生,研究員,主要從事地震波成像和油儲地球物理研究. E-mail: liuhong@mail.iggcas.ac.cn

    10.6038/cjg20161033

    P631

    2015-08-02,2015-12-31收修定稿

    張玉潔,劉洪,崔棟等. 2016. 基于擠噴流效應(yīng)的Russell流體因子推廣及應(yīng)用. 地球物理學(xué)報,59(10):3901-3908,

    Zhang Y J, Liu H, Cui D,et al. 2016. Construction and application of the Russell fluid factor with squirt flow effect.ChineseJ.Geophys. (in Chinese),59(10):3901-3908,doi:10.6038/cjg20161033.

    猜你喜歡
    噴流泊松比裂隙
    具有負(fù)泊松比效應(yīng)的紗線研發(fā)
    裂隙腦室綜合征的診斷治療新進(jìn)展
    負(fù)泊松比功能的結(jié)構(gòu)復(fù)合紡紗技術(shù)進(jìn)展
    毛紡科技(2020年6期)2021-01-06 03:41:48
    “慧眼”發(fā)現(xiàn)迄今距離黑洞最近的高速噴流
    考慮粘彈性泊松比的固體推進(jìn)劑蠕變型本構(gòu)模型①
    固體推進(jìn)劑粘彈性泊松比應(yīng)變率-溫度等效關(guān)系
    裂隙燈檢查的個性化應(yīng)用(下)
    噴流干擾氣動熱數(shù)值模擬的若干影響因素
    《老炮兒》:在時代裂隙中揚(yáng)棄焦慮
    電影新作(2016年1期)2016-02-27 09:16:48
    耀變體噴流高能電子譜的形成機(jī)制
    在线观看免费日韩欧美大片| 亚洲国产欧美一区二区综合| 国产av在哪里看| 亚洲精品一卡2卡三卡4卡5卡| 黄色丝袜av网址大全| 国产日本99.免费观看| 国产真实乱freesex| 两个人视频免费观看高清| 亚洲精品粉嫩美女一区| 搡老岳熟女国产| 亚洲成人免费电影在线观看| 又黄又粗又硬又大视频| 麻豆一二三区av精品| 久久草成人影院| 老司机午夜十八禁免费视频| 久久这里只有精品19| tocl精华| 一边摸一边抽搐一进一小说| 精品国产乱码久久久久久男人| 国产一区二区激情短视频| 国产精品一及| 99re在线观看精品视频| 日韩欧美精品v在线| 精品人妻1区二区| 男人舔奶头视频| 色哟哟哟哟哟哟| 精品一区二区三区视频在线观看免费| 国产亚洲精品久久久久久毛片| 全区人妻精品视频| 欧美日本亚洲视频在线播放| 精品一区二区三区av网在线观看| 最新美女视频免费是黄的| 国产精品永久免费网站| 久久伊人香网站| 久久久久精品国产欧美久久久| 久久午夜综合久久蜜桃| 日本三级黄在线观看| 亚洲av电影不卡..在线观看| www.精华液| 精品日产1卡2卡| 国产成人aa在线观看| 97超级碰碰碰精品色视频在线观看| 久久久久性生活片| 狂野欧美激情性xxxx| 日本三级黄在线观看| xxx96com| 国产精品一区二区免费欧美| 免费观看精品视频网站| 国内精品久久久久精免费| 日本一本二区三区精品| 国产精品野战在线观看| 女人爽到高潮嗷嗷叫在线视频| 久久精品国产综合久久久| 亚洲欧美日韩东京热| 无限看片的www在线观看| 老熟妇乱子伦视频在线观看| 亚洲精品在线美女| 午夜视频精品福利| 精品久久久久久,| av在线播放免费不卡| 亚洲成人国产一区在线观看| 床上黄色一级片| 国产精品av久久久久免费| 叶爱在线成人免费视频播放| 欧美一区二区国产精品久久精品 | 亚洲片人在线观看| 精品乱码久久久久久99久播| 国产精品久久久人人做人人爽| 国产黄a三级三级三级人| 久久久水蜜桃国产精品网| 琪琪午夜伦伦电影理论片6080| 精品国产乱码久久久久久男人| 最近最新免费中文字幕在线| 夜夜爽天天搞| 青草久久国产| 热99re8久久精品国产| x7x7x7水蜜桃| 欧美3d第一页| 校园春色视频在线观看| 观看免费一级毛片| 人妻丰满熟妇av一区二区三区| 国产精品 欧美亚洲| 亚洲男人的天堂狠狠| 精品第一国产精品| 不卡av一区二区三区| 老汉色∧v一级毛片| 伊人久久大香线蕉亚洲五| 成人国产一区最新在线观看| 波多野结衣高清无吗| 国产野战对白在线观看| 在线观看一区二区三区| 亚洲,欧美精品.| 最近视频中文字幕2019在线8| 18禁美女被吸乳视频| 欧美在线一区亚洲| 给我免费播放毛片高清在线观看| 麻豆国产av国片精品| 在线观看免费午夜福利视频| 99精品久久久久人妻精品| 欧美最黄视频在线播放免费| 99热这里只有精品一区 | 国产精品一区二区三区四区久久| 变态另类丝袜制服| 亚洲精品中文字幕一二三四区| 男女视频在线观看网站免费 | 日本熟妇午夜| 欧美 亚洲 国产 日韩一| 国产精品乱码一区二三区的特点| 日本熟妇午夜| 又黄又粗又硬又大视频| 国产精品一区二区精品视频观看| 亚洲色图av天堂| 欧美不卡视频在线免费观看 | 日韩av在线大香蕉| 午夜精品在线福利| 桃色一区二区三区在线观看| 国内精品久久久久久久电影| xxxwww97欧美| svipshipincom国产片| 亚洲国产欧美网| 成人高潮视频无遮挡免费网站| 久久国产精品影院| 亚洲欧美精品综合一区二区三区| 国产99白浆流出| 久久久国产欧美日韩av| 亚洲成人精品中文字幕电影| 精华霜和精华液先用哪个| 久久香蕉国产精品| 欧美日韩中文字幕国产精品一区二区三区| 久久婷婷人人爽人人干人人爱| 1024香蕉在线观看| www.999成人在线观看| 亚洲无线在线观看| 亚洲精品中文字幕一二三四区| 一边摸一边抽搐一进一小说| 好男人在线观看高清免费视频| e午夜精品久久久久久久| 午夜成年电影在线免费观看| 国产在线精品亚洲第一网站| netflix在线观看网站| 国产亚洲精品av在线| 757午夜福利合集在线观看| 欧美三级亚洲精品| 又粗又爽又猛毛片免费看| 国产精品久久久av美女十八| 亚洲精品色激情综合| 日日干狠狠操夜夜爽| 国产熟女午夜一区二区三区| 婷婷丁香在线五月| 最近最新中文字幕大全电影3| 一级毛片高清免费大全| 搡老熟女国产l中国老女人| 黑人巨大精品欧美一区二区mp4| 免费看a级黄色片| 亚洲国产欧洲综合997久久,| 日韩av在线大香蕉| 非洲黑人性xxxx精品又粗又长| www.精华液| 久久 成人 亚洲| 特级一级黄色大片| 我的老师免费观看完整版| 一a级毛片在线观看| 日韩三级视频一区二区三区| 99国产精品一区二区蜜桃av| 麻豆成人午夜福利视频| 精品久久久久久成人av| 亚洲av中文字字幕乱码综合| 亚洲av五月六月丁香网| 久久久久久九九精品二区国产 | 给我免费播放毛片高清在线观看| 99热6这里只有精品| 久久久精品大字幕| 亚洲午夜理论影院| 欧美精品亚洲一区二区| 午夜视频精品福利| 岛国视频午夜一区免费看| 狠狠狠狠99中文字幕| 女警被强在线播放| 国产片内射在线| 国产激情偷乱视频一区二区| 又黄又爽又免费观看的视频| 欧美一区二区精品小视频在线| 久久精品夜夜夜夜夜久久蜜豆 | 一本一本综合久久| 国产单亲对白刺激| 脱女人内裤的视频| 国语自产精品视频在线第100页| 波多野结衣高清无吗| 亚洲av美国av| 久久精品91无色码中文字幕| 久久这里只有精品19| 国产片内射在线| 欧美激情久久久久久爽电影| 婷婷精品国产亚洲av在线| 岛国视频午夜一区免费看| 男女视频在线观看网站免费 | 超碰成人久久| 18禁黄网站禁片午夜丰满| 亚洲欧美日韩高清专用| 国产成人欧美在线观看| 亚洲欧美日韩高清在线视频| 中文资源天堂在线| 亚洲国产高清在线一区二区三| 国产男靠女视频免费网站| 他把我摸到了高潮在线观看| 国产高清视频在线观看网站| 国产一区二区三区在线臀色熟女| 99国产精品一区二区蜜桃av| 精品无人区乱码1区二区| 精品久久久久久久久久久久久| 在线国产一区二区在线| 一本综合久久免费| 亚洲中文字幕一区二区三区有码在线看 | 欧美在线黄色| 天堂动漫精品| 亚洲欧美一区二区三区黑人| 真人一进一出gif抽搐免费| 国产av一区在线观看免费| 天天一区二区日本电影三级| av天堂在线播放| 99热只有精品国产| 午夜福利在线观看吧| 久久久久性生活片| 精品日产1卡2卡| 啦啦啦观看免费观看视频高清| 曰老女人黄片| 欧美日本亚洲视频在线播放| 19禁男女啪啪无遮挡网站| 久久久久国内视频| 国产69精品久久久久777片 | 国产欧美日韩一区二区精品| 老司机午夜十八禁免费视频| 国产av一区在线观看免费| 一级毛片精品| 黄色 视频免费看| 视频区欧美日本亚洲| 亚洲精品粉嫩美女一区| 大型av网站在线播放| 国内精品久久久久久久电影| 精品久久久久久久久久免费视频| 亚洲自偷自拍图片 自拍| 神马国产精品三级电影在线观看 | 久久草成人影院| 在线国产一区二区在线| 日本免费一区二区三区高清不卡| 免费av毛片视频| 两个人视频免费观看高清| 亚洲精品久久国产高清桃花| 99热6这里只有精品| 国产高清videossex| 亚洲一区中文字幕在线| 欧美 亚洲 国产 日韩一| 亚洲一区二区三区色噜噜| 夜夜看夜夜爽夜夜摸| 村上凉子中文字幕在线| 午夜福利免费观看在线| 18禁裸乳无遮挡免费网站照片| 99久久无色码亚洲精品果冻| 国产成人影院久久av| 欧美日韩亚洲国产一区二区在线观看| 久久久精品大字幕| 国产99白浆流出| 两个人免费观看高清视频| 免费观看人在逋| 桃色一区二区三区在线观看| 欧美乱码精品一区二区三区| 国产精品综合久久久久久久免费| 国产免费av片在线观看野外av| 99久久无色码亚洲精品果冻| 熟女少妇亚洲综合色aaa.| 亚洲专区字幕在线| 精品久久久久久,| 免费看美女性在线毛片视频| 别揉我奶头~嗯~啊~动态视频| 久久午夜综合久久蜜桃| 久久久久国产精品人妻aⅴ院| 一区福利在线观看| 午夜亚洲福利在线播放| 99精品久久久久人妻精品| 亚洲 欧美一区二区三区| 亚洲精品中文字幕在线视频| 宅男免费午夜| 亚洲精品久久国产高清桃花| 久久久国产成人精品二区| 特级一级黄色大片| 免费看a级黄色片| 欧美性猛交╳xxx乱大交人| 精品一区二区三区四区五区乱码| 一区二区三区国产精品乱码| 好男人电影高清在线观看| 床上黄色一级片| 成年人黄色毛片网站| 亚洲av成人精品一区久久| 久久婷婷成人综合色麻豆| 丰满人妻一区二区三区视频av | 十八禁人妻一区二区| 正在播放国产对白刺激| 一边摸一边做爽爽视频免费| 99精品欧美一区二区三区四区| 午夜福利视频1000在线观看| 亚洲avbb在线观看| 亚洲最大成人中文| www日本黄色视频网| 国产精品av久久久久免费| 夜夜夜夜夜久久久久| 亚洲av中文字字幕乱码综合| 久久精品人妻少妇| 亚洲av成人av| 国产av又大| 国产精品电影一区二区三区| 国产男靠女视频免费网站| 免费av毛片视频| 国产不卡一卡二| 麻豆成人av在线观看| 国产麻豆成人av免费视频| 成在线人永久免费视频| 精品国产超薄肉色丝袜足j| 黄色a级毛片大全视频| 精品午夜福利视频在线观看一区| 亚洲成a人片在线一区二区| 亚洲国产精品合色在线| 婷婷六月久久综合丁香| 高清毛片免费观看视频网站| 午夜精品一区二区三区免费看| 国产精品影院久久| 嫩草影院精品99| 国产成人av教育| 国产精品亚洲美女久久久| 亚洲美女视频黄频| 两性夫妻黄色片| 50天的宝宝边吃奶边哭怎么回事| 欧美日韩黄片免| 亚洲性夜色夜夜综合| 可以免费在线观看a视频的电影网站| 黄色a级毛片大全视频| 又粗又爽又猛毛片免费看| 99国产综合亚洲精品| 三级毛片av免费| 老司机深夜福利视频在线观看| 日韩精品青青久久久久久| 久久中文字幕人妻熟女| 蜜桃久久精品国产亚洲av| 欧美成人午夜精品| 国产免费av片在线观看野外av| 欧美成狂野欧美在线观看| 好看av亚洲va欧美ⅴa在| 欧美黑人精品巨大| 给我免费播放毛片高清在线观看| 啦啦啦韩国在线观看视频| 精品久久久久久,| a级毛片在线看网站| www.熟女人妻精品国产| 此物有八面人人有两片| 国内精品一区二区在线观看| 日本免费一区二区三区高清不卡| 精品久久久久久久毛片微露脸| 精品一区二区三区四区五区乱码| 国内精品一区二区在线观看| 一进一出好大好爽视频| 在线观看免费午夜福利视频| 中文在线观看免费www的网站 | 日本 av在线| 国产av又大| 伦理电影免费视频| 日本黄色视频三级网站网址| 国产亚洲欧美在线一区二区| 国产亚洲欧美98| 色综合婷婷激情| 亚洲自拍偷在线| 日韩av在线大香蕉| 级片在线观看| 午夜福利高清视频| 一个人免费在线观看电影 | www.熟女人妻精品国产| 国产爱豆传媒在线观看 | 好看av亚洲va欧美ⅴa在| 在线观看一区二区三区| 两性夫妻黄色片| 欧美乱色亚洲激情| 国产精品免费一区二区三区在线| 日韩有码中文字幕| 搞女人的毛片| 91麻豆精品激情在线观看国产| 欧美日韩一级在线毛片| 女人爽到高潮嗷嗷叫在线视频| 国产欧美日韩一区二区三| 午夜激情福利司机影院| 99国产精品一区二区三区| 免费在线观看完整版高清| av片东京热男人的天堂| 久久久久国产精品人妻aⅴ院| 亚洲国产欧洲综合997久久,| 亚洲欧美激情综合另类| 国产精品乱码一区二三区的特点| 日本五十路高清| 欧美黑人巨大hd| 国产av麻豆久久久久久久| 亚洲国产看品久久| 别揉我奶头~嗯~啊~动态视频| 成人国产综合亚洲| 国内揄拍国产精品人妻在线| 成人av在线播放网站| 欧美性长视频在线观看| 国产午夜精品久久久久久| 国产精品久久电影中文字幕| 91字幕亚洲| 精品久久久久久久末码| 免费观看精品视频网站| 美女午夜性视频免费| 麻豆国产av国片精品| 日韩精品中文字幕看吧| 18禁国产床啪视频网站| 禁无遮挡网站| 搡老熟女国产l中国老女人| 亚洲精品久久成人aⅴ小说| 久久久国产成人免费| 天天一区二区日本电影三级| 变态另类丝袜制服| 国产激情久久老熟女| 好男人电影高清在线观看| 日韩av在线大香蕉| 日日夜夜操网爽| 国产精品久久久人人做人人爽| 久久这里只有精品19| 一级黄色大片毛片| 丝袜美腿诱惑在线| 两个人看的免费小视频| 日韩成人在线观看一区二区三区| 国产av又大| 亚洲专区国产一区二区| 国产精品自产拍在线观看55亚洲| 久久久久精品国产欧美久久久| 无遮挡黄片免费观看| 人成视频在线观看免费观看| 成人精品一区二区免费| 精品国产美女av久久久久小说| 两人在一起打扑克的视频| 全区人妻精品视频| 国产精品野战在线观看| 日韩精品青青久久久久久| 高清毛片免费观看视频网站| 99久久99久久久精品蜜桃| 18禁国产床啪视频网站| 桃色一区二区三区在线观看| 国产精品久久久久久久电影 | 亚洲自偷自拍图片 自拍| 99精品欧美一区二区三区四区| 在线观看免费日韩欧美大片| 老司机深夜福利视频在线观看| av有码第一页| 在线播放国产精品三级| 亚洲成人中文字幕在线播放| 非洲黑人性xxxx精品又粗又长| 国产精品久久久久久久电影 | 女人爽到高潮嗷嗷叫在线视频| 色综合亚洲欧美另类图片| e午夜精品久久久久久久| 在线a可以看的网站| 国产精品自产拍在线观看55亚洲| 99热这里只有精品一区 | 校园春色视频在线观看| www.自偷自拍.com| 日韩高清综合在线| 欧美绝顶高潮抽搐喷水| 国产亚洲欧美在线一区二区| 国产在线精品亚洲第一网站| 国产一区二区在线观看日韩 | 亚洲av五月六月丁香网| 听说在线观看完整版免费高清| 日本一本二区三区精品| 少妇人妻一区二区三区视频| 后天国语完整版免费观看| 又紧又爽又黄一区二区| 亚洲七黄色美女视频| 国产激情偷乱视频一区二区| 日韩欧美三级三区| 欧美国产日韩亚洲一区| 午夜福利18| 女警被强在线播放| 久久久久久九九精品二区国产 | 看黄色毛片网站| 国产av在哪里看| 亚洲成av人片免费观看| 亚洲专区国产一区二区| 国产高清有码在线观看视频 | 国产精品久久久久久亚洲av鲁大| 久久九九热精品免费| 欧美日韩瑟瑟在线播放| 嫩草影院精品99| 久久久久久久久免费视频了| 国产成人系列免费观看| 1024手机看黄色片| 欧美精品亚洲一区二区| 18禁美女被吸乳视频| 中文字幕高清在线视频| 丰满人妻一区二区三区视频av | 久久精品91蜜桃| 国产一区二区三区在线臀色熟女| 啦啦啦观看免费观看视频高清| 午夜免费观看网址| 国产精品久久电影中文字幕| 亚洲国产精品999在线| 美女 人体艺术 gogo| 人妻夜夜爽99麻豆av| 午夜免费激情av| 变态另类成人亚洲欧美熟女| 少妇裸体淫交视频免费看高清 | av天堂在线播放| 麻豆成人午夜福利视频| 亚洲人成网站在线播放欧美日韩| 18禁美女被吸乳视频| 香蕉久久夜色| 极品教师在线免费播放| 黑人操中国人逼视频| 久久精品综合一区二区三区| 亚洲人成网站高清观看| 99久久精品热视频| 亚洲av电影不卡..在线观看| av国产免费在线观看| 国产区一区二久久| 亚洲 国产 在线| 后天国语完整版免费观看| av福利片在线观看| а√天堂www在线а√下载| 亚洲va日本ⅴa欧美va伊人久久| 久久中文字幕一级| 一本综合久久免费| 日本精品一区二区三区蜜桃| 亚洲第一电影网av| 久久精品夜夜夜夜夜久久蜜豆 | 757午夜福利合集在线观看| 天天添夜夜摸| 亚洲精品av麻豆狂野| 亚洲国产欧洲综合997久久,| 欧美大码av| 久久天堂一区二区三区四区| 国产区一区二久久| 亚洲精品一区av在线观看| 久久香蕉激情| 一进一出抽搐gif免费好疼| 在线观看免费午夜福利视频| 一进一出抽搐gif免费好疼| 国产精品美女特级片免费视频播放器 | 一个人免费在线观看的高清视频| 国产av麻豆久久久久久久| 女生性感内裤真人,穿戴方法视频| 中文在线观看免费www的网站 | 亚洲国产中文字幕在线视频| 五月伊人婷婷丁香| 精品欧美一区二区三区在线| 国产男靠女视频免费网站| 美女 人体艺术 gogo| 在线a可以看的网站| 男女视频在线观看网站免费 | 日本免费a在线| 欧美日韩国产亚洲二区| 制服丝袜大香蕉在线| 搞女人的毛片| 久久久精品欧美日韩精品| 欧美日本视频| 久久久久久国产a免费观看| 精品一区二区三区av网在线观看| 午夜福利欧美成人| 日本撒尿小便嘘嘘汇集6| 又粗又爽又猛毛片免费看| 在线看三级毛片| 国产野战对白在线观看| 又爽又黄无遮挡网站| 午夜福利欧美成人| 亚洲欧美精品综合久久99| 国产精品久久久av美女十八| 成年女人毛片免费观看观看9| 亚洲自拍偷在线| 97超级碰碰碰精品色视频在线观看| 欧美又色又爽又黄视频| 欧美色欧美亚洲另类二区| 国产精品久久久久久亚洲av鲁大| 无人区码免费观看不卡| 亚洲中文字幕一区二区三区有码在线看 | 亚洲美女视频黄频| 国产又黄又爽又无遮挡在线| 亚洲,欧美精品.| 中文字幕久久专区| 女生性感内裤真人,穿戴方法视频| 欧美国产日韩亚洲一区| 在线免费观看的www视频| 中文在线观看免费www的网站 | 高清在线国产一区| 久久精品国产亚洲av香蕉五月| 国产成人啪精品午夜网站| 国产精品久久久久久精品电影| 日本一本二区三区精品| 天堂√8在线中文| 国产1区2区3区精品| 国内毛片毛片毛片毛片毛片| 久久精品亚洲精品国产色婷小说| 欧美zozozo另类| 熟女电影av网| 欧美色视频一区免费| 午夜福利在线在线| 精品久久蜜臀av无| 婷婷六月久久综合丁香| 亚洲天堂国产精品一区在线| 久久精品影院6| 波多野结衣巨乳人妻| 国产黄色小视频在线观看| 天天一区二区日本电影三级| 狠狠狠狠99中文字幕| 一二三四在线观看免费中文在| 国产精华一区二区三区| 一进一出抽搐gif免费好疼| 久久久精品国产亚洲av高清涩受| 国产在线观看jvid| 法律面前人人平等表现在哪些方面| 精品国产美女av久久久久小说|