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

    基于楊氏模量、泊松比和各向異性梯度的裂縫介質(zhì)AVAZ反演方法

    2015-06-27 05:54:44杜炳毅楊午陽王恩利張廣智高建虎
    石油物探 2015年2期
    關(guān)鍵詞:楊氏模量縱波泊松比

    杜炳毅,楊午陽,王恩利,張廣智,高建虎

    (1.中國石油勘探開發(fā)研究院西北分院,甘肅蘭州730020;2.中國石油大學(xué)(華東)地球科學(xué)與技術(shù)學(xué)院,山東青島266580)

    基于楊氏模量、泊松比和各向異性梯度的裂縫介質(zhì)AVAZ反演方法

    杜炳毅1,楊午陽1,王恩利1,張廣智2,高建虎1

    (1.中國石油勘探開發(fā)研究院西北分院,甘肅蘭州730020;2.中國石油大學(xué)(華東)地球科學(xué)與技術(shù)學(xué)院,山東青島266580)

    為了全面地描述裂縫介質(zhì)的儲層特征、流體性質(zhì)和各向異性特征,提出了由寬方位疊前道集數(shù)據(jù)直接獲取裂縫介質(zhì)彈性參數(shù)和各向異性梯度參數(shù)的振幅隨方位角變化(Amplitude Versus Azimuth,AVAZ)反演方法。首先推導(dǎo)了基于楊氏模量、泊松比和各向異性梯度的各向異性AVO方程,通過與Ruger近似進行對比,分析了該方程的近似精度。然后利用實際測井資料和二維逆掩斷層模型建立了不同方位角的合成疊前角度道集,對未加噪聲和信噪比為3的寬方位疊前角度道集進行了AVAZ反演方法測試,結(jié)果表明,未加噪聲和信噪比為3的合成寬方位疊前道集均能反演得到符合地震反演精度的楊氏模量(E),泊松比(σ)和各向異性梯度(Γ),噪聲對楊氏模量的反演影響較小,不含噪聲時泊松比和各向異性梯度的反演精度更高,與真實值吻合程度較高。

    楊氏模量;泊松比;各向異性梯度;裂縫介質(zhì);AVAZ反演

    Thomsen[1]提出了弱各向異性介質(zhì)理論,引入了弱各向異性參數(shù),為各向異性的研究奠定了基礎(chǔ);Ruger[2]對HTI介質(zhì)的反射與透射特征進行了詳細(xì)的研究,推導(dǎo)了縱波反射系數(shù)近似公式,并分析了方位各向異性與入射角和方位角的關(guān)系;Bakulin等[3]和Sayers等[4]利用帶有線性滑動邊界條件的有限平行裂縫和硬幣型裂縫理論來估計各向異性參數(shù)與裂縫結(jié)構(gòu)之間的關(guān)系,進而估測裂縫的發(fā)育方向和裂縫密度,利用各向異性差值判斷裂縫的含油氣性。AVAZ反演中彈性參數(shù)和各向異性參數(shù)的反演具有很多的不確定性,不同的學(xué)者對此進行了相關(guān)的研究。Subhashis等[5]利用縱波和轉(zhuǎn)換波隨入射角和方位角的變化特征進行AVAZ反演,得到穩(wěn)健的各向異性參數(shù)估計值;Downtown等[6-7]運用貝葉斯理論估計HTI介質(zhì)AVAZ反演的彈性參數(shù)預(yù)測值的不確定性,提出了同步彈性參數(shù)反演,加入巖石物理模型的約束,克服了AVAZ反演的固有缺陷,從而定量地解釋AVAZ反演的穩(wěn)定性。

    為了更好地識別裂縫儲層特征,孔麗云等[8]深入分析了裂縫孔隙、裂縫弱度、裂縫滲透率等參數(shù)對雙孔隙裂縫介質(zhì)地震各向異性的影響。P波反射振幅與彈性參數(shù)之間的關(guān)系是裂縫性儲層反演研究的主要問題之一,劉軍迎等[9]認(rèn)為裂縫儲層的反射振幅隨著觀測方位的變化而變化,表現(xiàn)出明顯的方位各向異性特征,可以應(yīng)用該特征來預(yù)測碳酸鹽巖儲層中的裂縫密度及發(fā)育方向。張廣智等[10]研究了裂縫彈性參數(shù)的各向異性AVO反演方法,準(zhǔn)確地得到了反映儲層特征的縱波阻抗、橫波阻抗以及各向異性梯度。陳懷震等[11-12]研究了縱波動力學(xué)特征與裂縫流體因子、縫隙流體因子、裂縫密度、入射角和方位角等因素之間的關(guān)系特征,通過合成的模型數(shù)據(jù)對該方法的適用性進行了詳細(xì)的分析。宗兆云等[13]分析了縱波模量、橫波模量、楊氏模量和泊松比對儲層反射特征的影響,提出了利用疊前AVO反演方法直接估算縱波模量、橫波模量、楊氏模量和泊松比的算法,避免了間接求取帶來的誤差累積,為裂縫型儲層的識別提供了更加可靠的理論依據(jù);宗兆云等[14]通過建立縱橫波模量、密度以及各向異性差值與縱波反射系數(shù)之間的關(guān)系,進一步估測地層的地應(yīng)力,這對頁巖氣儲層的評價是非常有利的。裂縫儲層的彈性阻抗也表現(xiàn)為強方位各向異性特征,通過方位各向異性彈性阻抗反演可以獲取表征裂縫性質(zhì)的彈性參數(shù)和各向異性參數(shù)[15-16]。

    楊氏模量是表示巖石抗壓縮能力的物理量,表征了儲層的巖性特性;泊松比是橫向應(yīng)變與縱向應(yīng)變的比值絕對值,與儲層流體性質(zhì)有關(guān)。為了全面地描述裂縫介質(zhì)的儲層特征、流體性質(zhì)和各向異性特征,本文在上述研究的基礎(chǔ)上,建立了裂縫介質(zhì)的縱波反射系數(shù)與楊氏模量、泊松比和各向異性梯度之間的關(guān)系,利用模型數(shù)據(jù)對新的AVAZ反演算法進行測試,取得了較好的應(yīng)用效果。

    1 方法原理

    裂縫介質(zhì)的各向異性特征可以由Thomsen提出的弱各向異性參數(shù)ε(v),δ(v),γ來表征,它們與彈性矩陣之間的關(guān)系為[1]:

    (1)

    Ruger研究了各向異性半空間界面的地震波反射特征,推導(dǎo)了縱波反射系數(shù)隨入射角、方位角變化的關(guān)系式[2],如(2)式所示:

    (2)

    按照擾動理論,分為背景場和擾動場,可以將縱波的AVO近似方程表示為:

    (3)

    (4)

    (5)

    宗兆云推導(dǎo)了縱波近似反射系數(shù)與楊氏模量(E),泊松比(σ)和密度反射系數(shù)之間的變化關(guān)系式[14]:

    (6)

    (7)

    式中:C是一個常數(shù);g為縱波速度的冪指數(shù)。

    對其求一階偏導(dǎo),求出密度反射系數(shù)與縱波速度反射系數(shù)的關(guān)系:

    (8)

    最后,可以求得密度反射系數(shù)與楊氏模量反射系數(shù)、泊松比反射系數(shù)的函數(shù)表示式:

    (9)

    將其代入到(6)式中,進一步得到:

    (10)

    Chen等[11]指出,在入射角小于30°的情況下,可以令sin2θtan2θ≈0,(5)式可以進一步簡化為:

    (11)

    Ruger定義各向異性梯度Γ為:

    (12)

    最終,可以得到基于楊氏模量、泊松比和各向異性梯度的P波AVAZ近似方程:

    (13)

    在上述推導(dǎo)的基礎(chǔ)上,選取兩層的裂縫介質(zhì)模型,相關(guān)參數(shù)如表1所示,對新推導(dǎo)的縱波AVAZ近似方程進行精度分析。由Ruger近似得到曲線和新推導(dǎo)公式的近似曲線對比結(jié)果(圖1),對比發(fā)現(xiàn),在入射角較大的情況下,新推導(dǎo)方程的近似精度與Ruger近似吻合程度較高,因此,利用該方程進行裂縫介質(zhì)AVAZ反演可以提高裂縫預(yù)測的精度。

    表1 裂縫模型參數(shù)

    在不同方位角和入射角的情況下,可以將(13)式表示為矩陣的形式:

    (14)

    圖1 新推導(dǎo)方程與Ruger近似方程的精度對比

    令:

    (15)

    (16)

    (17)

    (14)式可以表示為:

    d=Gm

    (18)

    利用最小二乘算法求取出模型參數(shù)的最佳估算值:

    (19)

    2 實際模型測試

    2.1 實際工區(qū)測井資料合成疊前角度道集的AVAZ反演方法測試

    為了驗證反演算法的有效性,本文選取某實際工區(qū)的縱波速度、橫波速度和各向異性參數(shù)等測井曲線(圖2),利用新推導(dǎo)的公式求取不同入射角(0~30°)和方位角(0,30°,60°,90°)的P波AVO反射系數(shù),并且與40Hz的雷克子波褶積得到不同方位角的疊前合成角度道集,對其添加信噪比為3的隨機噪聲。分別對不含噪聲和信噪比為3的合成數(shù)據(jù)進行算法測試。圖3是密度和縱波速度的冪指數(shù)擬合,擬合點基本分布在公式(7)所表示的曲線附近,說明密度和縱波速度之間滿足關(guān)系式(7)。

    圖4是不含噪聲時合成的不同方位角(0,30°,60°,90°)的疊前角度道集;圖5是未加噪聲合成數(shù)據(jù)經(jīng)過基于楊氏模量、泊松比和各向異性梯度反演得到的反演結(jié)果與真實結(jié)果的對比;圖6是信噪比為3時合成的不同方位角(0,30°,60°,90°)的疊前角度道集;圖7是信噪比為3的合成數(shù)據(jù)經(jīng)過基于楊氏模量、泊松比和各向異性梯度反演得到的反演結(jié)果與真實結(jié)果的對比。對比圖7與圖5可以看出,噪聲對楊氏模量的反演影響較?。徊缓肼晻r泊松比和各向異性梯度的反演精度更高,與真實值吻合程度較高。

    圖2 某實際工區(qū)的測井曲線

    圖3 密度和縱波速度的冪指數(shù)擬合

    圖4 不加噪聲時合成的不同方位角的疊前角度道集

    圖5 未加噪聲合成數(shù)據(jù)的楊氏模量、泊松比和各向異性梯度反演值(紅色)與真實值(藍(lán)色)對比

    圖6 信噪比為3時合成的不同方位角的疊前角度道集

    圖7 信噪比為3合成數(shù)據(jù)的楊氏模量、泊松比和各向異性梯度反演值(紅色)與真實值(藍(lán)色)對比

    2.2 二維逆掩斷層模型測試

    再應(yīng)用上述AVAZ反演方法對二維逆掩斷層模型進行算法測試,進一步驗證該算法的合理性。圖8給出了二維逆掩斷層的原始模型剖面,分別是楊氏模量剖面(圖8a)、泊松比剖面(圖8b)和各向異性梯度剖面(圖8c)。通過正演模擬得到二維逆掩斷層模型不同方位角(0,30°,60°,90°)的合成地震角度道集,并且加入信噪比為3的隨機噪聲。利用含噪聲合成地震記錄直接通過AVAZ反演獲取二維逆掩斷層模型的楊氏模量、泊松比和各向異性梯度,反演結(jié)果如圖9所示。由圖9可以看到,楊氏模量、泊松比和各向異性梯度均能準(zhǔn)確反映該模型的構(gòu)造特征,清晰地展示了地層界面,具有很強的橫向連續(xù)性。測試結(jié)果進一步驗證了新的P波AVAZ近似公式的正確性和反演算法的有效性。

    圖8 二維逆掩斷層模型的楊氏模量(a)、泊松比(b)和各向異性梯度(c)的原始剖面

    圖9 二維逆掩斷層模型合成道集(信噪比為3)的楊氏模量(a)、泊松比(b)和各向異性梯度(c)的反演結(jié)果

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

    本文推導(dǎo)了基于楊氏模量、泊松比和各向異性梯度的方位各向異性AVO方程,提出了一種AVAZ反演算法;選用某實際工區(qū)的測井?dāng)?shù)據(jù)和二維逆掩斷層模型合成不同方位角的疊前角度道集測試了該AVAZ反演算法,為裂縫介質(zhì)參數(shù)反演的發(fā)展提供了新的思路。通過以上研究工作可以得出以下結(jié)論與認(rèn)識。

    1) 對新推導(dǎo)的方位各向異性AVO方程進行精度分析,通過與Ruger近似方程對比,可以看出新推導(dǎo)的方程能夠滿足裂縫介質(zhì)AVAZ反演的精度要求。

    2) 選擇實際井資料模型和二維逆掩斷層模型,合成不同方位角的疊前角度道集對AVAZ反演算法進行測試,結(jié)果表明,楊氏模量、泊松比和各向異性梯度的反演結(jié)果與真實值吻合程度較高,能夠反演裂縫儲層的構(gòu)造特征和巖性特征等信息。

    3) 值得指出的是,裂縫介質(zhì)中流體的識別是裂縫介質(zhì)研究的主要目的之一,因此,下一步需要對基于上述AVAZ反演結(jié)果的裂縫流體識別進行深入的研究。

    [1] Thomsen L.Weak elastic anisotropy[J].Geophysics,1986,51(10):1954-1966

    [2] Ruger A.Reflection coefficients and azimuthal AVO analysis in anisotropic media[D].Colorado:Colorado School of Mines Golden,1996

    [3] Bakulin A,Grechka V,Tsvankin I.Estimation of fracture parameters from reflection seismic data—part I:HTI model due to a single fracture set[J].Geophysics,2000,65(6):1788-1802

    [4] Sayers C M,Simon D.Azimuth-dependent AVO in reservoirs containing non-orthogonal fracture sets[J].Geophysical Prospecting,2001,49(1):100-106

    [5] Subhashis M,Diana G.Azimuthal reflectivity and quantitative evaluation of anisotropic parameters from seismic data a feasibility study[J].Expanded Abstracts of 75thAnnual Internat SEG Mtg,2005,182-185

    [6] Downton J,Gray D.AVAZ parameter uncertainty estimation[J].Expanded Abstracts of 76thAnnual Internat SEG Mtg,2006,234-238

    [7] Downton J,Roure R.Azimuthal simultaneous elastic inversion for fracture detection[J].Expanded Abstracts of 80thAnnual Internat SEG Mtg,2010,263-267

    [8] 孔麗云,王一博,楊慧珠.裂縫誘導(dǎo)HTI雙孔隙介質(zhì)中的裂縫參數(shù)分析[J].地球物理學(xué)報,2012,55(1):189-196 Kong L Y,Wang Y B,Yang H Z.Fracture parameters analyses in fracture -induced HTI double-porosity medium[J].Chinese Journal of Geophysics,2012,55(1):189-196

    [9] 劉軍迎,雍學(xué)善,楊午陽,等.基于疊前方位振幅的大港—埕海地區(qū)奧陶系風(fēng)化殼裂縫儲層的疊前預(yù)測[J].地球物理學(xué)進展,2012,27(4):1588-1597 Liu J Y,Yong X S,Yang W Y,et al.Pre-stack fracture reservoir bed estimation of Ordovician weather layer of Dagang-Chenghai area based on azimuthal pre-stack seismic amplitude[J].Progress in Geophysics,2012,27(4):1588-1597

    [10] 張廣智,陳懷震,印興耀,等.基于各向異性AVO的裂縫彈性參數(shù)疊前反演方法[J].吉林大學(xué)學(xué)報(地球科學(xué)版),2012,42(3):845-871 Zhang G Z,Chen H Z,Yin X Y,et al.Method of fracture elastic parameter inversion based on anisotropy AVO[J].Journal of Jilin University(Earth Science Edition),2012,42(3):845-871

    [11] Chen H Z,Zhang G,Yin X.AVAZ inversion for elastic parameter and fracture fluid factor[J].Expanded Abstracts of 82ndAnnual Internat SEG Mtg,2012,1-5

    [12] 陳懷震,印興耀,高成國,等.基于各向異性巖石物理的縫隙流體因子AVAZ反演[J].地球物理學(xué)報,2014,57(3):968-978 Chen H Z,Yin X Y,Gao C G,et al.AVAZ inversion for fluid factor based on fracture anisotropic rock physics theory[J].Chinese Journal of Geophysics,2014,57(3):968-978

    [13] 宗兆云,印興耀,吳國忱.基于疊前地震縱橫波模量直接反演的流體檢測方法[J].地球物理學(xué)報,2011,55(1):284-292 Zong Z Y,Yin X Y,Wu G C.Reflection coefficient equation and pre-stack inversion with Young’s modulus and poisson ratio[J].Chinese Journal of Geophysics,2011,55(1):284-292

    [14] 宗兆云,印興耀,張峰,等.楊氏模量和泊松比反射系數(shù)近似方程及疊前地震反演[J].地球物理學(xué)報,2012,55(11):3786-3794 Zong Z Y,Yin X Y,Zhang F,et al.Reflection coefficient equation and pre-stack seismic inversion with Young’s modulus and poisson ratio[J].Chinese Journal of Geophysics,2012,55(11):3786-3794

    [15] Zong Y Z,Yin X Y,Wu G C.AVAZ inversion and stress evaluation in heterogeneous medium[J].Expanded Abstracts of 83rdAnnual Internat SEG Mtg,2013,428-431

    [16] 陳懷震,印興耀,杜炳毅,等.裂縫型碳酸鹽巖儲層方位各向異性彈性阻抗反演[J].地球物理學(xué)進展,2013,28(6):3073-3079 Chen H Z,Yin X Y,Du B Y,et al.Azimuthal anisotropic elastic impedance inversion in fractured layered carbonate rock reservoir[J].Progress in Geophysics,2013,28(6):3073-3079

    [17] 竇喜英,韓立國,劉春成,等.裂縫型HTI介質(zhì)中的彈性阻抗[J].吉林大學(xué)學(xué)報地球科學(xué)版,2012,42(4):1192-1197 Dou X Y,Han L G,Liu C C,et al.Elastic impedance in fracture medai[J].Journal of Jilin University(Earth Science Edition),2012,42(4):1192-1197

    (編輯:陳 杰)

    AVAZ inversion based on Young’s modulus,Poisson’s ratio and anisotropy gradient in fractured media

    Du Bingyi1,Yang Wuyang1,Wang Enli1,Zhang Guangzhi2,Gao Jianhu1

    (1.NorthwestBranch,ResearchInstituteofPetroleumExploration&Development,Petrochina,Lanzhou730020,China; 2.SchoolofGeosciences,ChinaUniversityofPetroleum,Qingdao266580,China)

    In order to fully describe the reserroir characteristics,fluid properties and anisotropy in fractured media,we proposed an algorithm of Amplitude Versus Azimuth (AVAZ) inversion to obtain elastic parameters and anisotropy parameters of fractured media by using wide azimuthally pre-stack seismic data.Firstly,a new P-wave azimuth AVO approximate equation based on Young’s modulus,Poisson’s ratio and anisotropy gradient is deduced.The precision of this new equation was analyzed and compared with Ruger’s equation.Then,we utilized actual log data and 2D overthrust model to generate synthetic pre-stack angle gathers of wide-azimuth data.The algorithm is tested through the wide-azimuth synthetic angle gathers without noise and with random noise (SNR=3) respectively.The result shows that the inverted Young’s modulus,Poisson’s ratio and anisotropy gradient without noise is estimated exactly and agreed with the original data.Meanwhile,the inversion result of the data with random noise (SNR=3) is stable and can meet the demand of seismic inversion.Noise causes less influence on the inversion of Young’s modulus,the invesion accuracy of Poisson’s ratio and anisotropy gradient is higher without noise which is well coinciding with the actual data.

    Young’s modulus,Poisson’s ratio,anisotropy gradient,fractured media

    2014-09-13;改回日期:2014-12-19。

    杜炳毅(1985—),男,碩士,主要從事疊前地震反演和儲層流體識別方面的研究。

    國家科技重大專項(2011ZX05019-008-005)資助。

    P631

    A

    1000-1441(2015)02-0218-08

    10.3969/j.issn.1000-1441.2015.02.014

    猜你喜歡
    楊氏模量縱波泊松比
    武漢大學(xué)研究團隊發(fā)現(xiàn)迄今“最剛強”物質(zhì)
    河南科技(2023年10期)2023-06-07 13:33:44
    具有負(fù)泊松比效應(yīng)的紗線研發(fā)
    負(fù)泊松比功能的結(jié)構(gòu)復(fù)合紡紗技術(shù)進展
    毛紡科技(2020年6期)2021-01-06 03:41:48
    考慮粘彈性泊松比的固體推進劑蠕變型本構(gòu)模型①
    固體推進劑粘彈性泊松比應(yīng)變率-溫度等效關(guān)系
    黃257井區(qū)疊前縱波方位各向異性裂縫分布預(yù)測
    近距二次反射式楊氏模量測量儀簡介
    物理實驗(2017年2期)2017-03-21 07:20:41
    變截面階梯桿中的縱波傳播特性實驗
    拉伸法測楊氏模量中的橫梁形變對實驗的影響
    關(guān)于縱波的波峰和波谷的探討
    物理與工程(2010年3期)2010-03-25 10:02:04
    成人国产一区最新在线观看| 欧美日韩瑟瑟在线播放| av网站在线播放免费| 狠狠狠狠99中文字幕| 亚洲美女黄片视频| 精品久久久精品久久久| av片东京热男人的天堂| 国产成人精品久久二区二区91| 精品一区二区三卡| 又黄又爽又免费观看的视频| 免费日韩欧美在线观看| 丰满迷人的少妇在线观看| 日韩大码丰满熟妇| 在线免费观看的www视频| 老熟妇仑乱视频hdxx| 欧美日韩国产mv在线观看视频| 婷婷丁香在线五月| 18禁美女被吸乳视频| 丰满的人妻完整版| 丝袜在线中文字幕| 国产成人精品久久二区二区免费| 女人被狂操c到高潮| 欧美精品人与动牲交sv欧美| 黄片小视频在线播放| 亚洲专区中文字幕在线| 欧美国产精品va在线观看不卡| 免费av中文字幕在线| 69精品国产乱码久久久| 国产精品一区二区在线观看99| 日韩欧美一区二区三区在线观看 | 国产精品98久久久久久宅男小说| 国产aⅴ精品一区二区三区波| 99热国产这里只有精品6| 狠狠婷婷综合久久久久久88av| 国产精品国产av在线观看| 婷婷成人精品国产| 一本综合久久免费| 国产高清videossex| 欧美日韩一级在线毛片| 51午夜福利影视在线观看| 欧美老熟妇乱子伦牲交| 黑人巨大精品欧美一区二区蜜桃| 女性被躁到高潮视频| 亚洲在线自拍视频| e午夜精品久久久久久久| 中亚洲国语对白在线视频| 黄色视频不卡| 巨乳人妻的诱惑在线观看| 一二三四在线观看免费中文在| 一级毛片精品| 老司机福利观看| 成年人黄色毛片网站| 中文字幕人妻丝袜制服| 亚洲第一青青草原| 91成人精品电影| 黑人巨大精品欧美一区二区mp4| 在线观看免费高清a一片| 高清av免费在线| 国产熟女午夜一区二区三区| 国产亚洲精品久久久久5区| a级片在线免费高清观看视频| av天堂在线播放| 精品少妇一区二区三区视频日本电影| 亚洲色图综合在线观看| 99久久综合精品五月天人人| 制服诱惑二区| 久久精品国产亚洲av高清一级| 国产精品二区激情视频| tube8黄色片| 丝袜在线中文字幕| 亚洲午夜精品一区,二区,三区| 日日夜夜操网爽| 法律面前人人平等表现在哪些方面| 麻豆国产av国片精品| 亚洲av日韩在线播放| 99在线人妻在线中文字幕 | 99久久综合精品五月天人人| 黄色丝袜av网址大全| 久9热在线精品视频| 51午夜福利影视在线观看| 久久久水蜜桃国产精品网| 欧美黑人精品巨大| 天天躁夜夜躁狠狠躁躁| 少妇 在线观看| 女人爽到高潮嗷嗷叫在线视频| 另类亚洲欧美激情| 精品一区二区三区视频在线观看免费 | 国产精品亚洲一级av第二区| 美女福利国产在线| 看片在线看免费视频| 欧美中文综合在线视频| 极品人妻少妇av视频| 99riav亚洲国产免费| 亚洲九九香蕉| 十八禁网站免费在线| 欧美精品高潮呻吟av久久| 12—13女人毛片做爰片一| 亚洲精品一卡2卡三卡4卡5卡| 欧美乱码精品一区二区三区| 黄色a级毛片大全视频| www.精华液| 国产亚洲欧美98| 69精品国产乱码久久久| 午夜福利在线观看吧| 黄频高清免费视频| 视频区欧美日本亚洲| 1024香蕉在线观看| 国产淫语在线视频| 精品亚洲成a人片在线观看| 成人av一区二区三区在线看| 一二三四在线观看免费中文在| 亚洲自偷自拍图片 自拍| 一区二区三区国产精品乱码| 两个人看的免费小视频| 成熟少妇高潮喷水视频| 亚洲欧美激情综合另类| 午夜精品久久久久久毛片777| 久久精品91无色码中文字幕| 51午夜福利影视在线观看| 女人久久www免费人成看片| 午夜免费观看网址| av欧美777| 久久国产乱子伦精品免费另类| 亚洲全国av大片| 超碰97精品在线观看| 麻豆成人av在线观看| 久久精品亚洲熟妇少妇任你| 中文字幕高清在线视频| 怎么达到女性高潮| 国产精品一区二区在线不卡| 欧美另类亚洲清纯唯美| 国产成人免费观看mmmm| 18禁黄网站禁片午夜丰满| 午夜精品在线福利| a在线观看视频网站| 久久天躁狠狠躁夜夜2o2o| 亚洲av熟女| 中文字幕另类日韩欧美亚洲嫩草| 国产精品.久久久| 十八禁网站免费在线| 欧美黄色片欧美黄色片| av在线播放免费不卡| 亚洲美女黄片视频| 免费看a级黄色片| 伊人久久大香线蕉亚洲五| 成人av一区二区三区在线看| 国产精品久久久久久精品古装| 一边摸一边做爽爽视频免费| 天堂俺去俺来也www色官网| 女人高潮潮喷娇喘18禁视频| 自线自在国产av| 国产高清视频在线播放一区| 一区二区三区激情视频| 亚洲成人免费电影在线观看| 一级黄色大片毛片| 少妇 在线观看| 久热爱精品视频在线9| 国产精品99久久99久久久不卡| 下体分泌物呈黄色| 岛国在线观看网站| 欧美日韩瑟瑟在线播放| 色婷婷久久久亚洲欧美| 色婷婷久久久亚洲欧美| 国产精品乱码一区二三区的特点 | 久久久国产精品麻豆| 日韩熟女老妇一区二区性免费视频| 亚洲熟妇熟女久久| 在线天堂中文资源库| 国产午夜精品久久久久久| 欧美日本中文国产一区发布| 日韩三级视频一区二区三区| 日本黄色日本黄色录像| 大香蕉久久网| 国产亚洲欧美在线一区二区| 国产三级黄色录像| 欧美日韩中文字幕国产精品一区二区三区 | 国产成人av教育| 两性夫妻黄色片| 两性夫妻黄色片| 波多野结衣一区麻豆| 男人操女人黄网站| 精品熟女少妇八av免费久了| 久久香蕉激情| 国产欧美日韩综合在线一区二区| av福利片在线| 国产成人av教育| 两人在一起打扑克的视频| 国产区一区二久久| 久久久久视频综合| 国产亚洲欧美精品永久| 亚洲国产中文字幕在线视频| 午夜成年电影在线免费观看| 国产精品永久免费网站| av天堂久久9| 国产在线精品亚洲第一网站| 久久人人爽av亚洲精品天堂| 人人妻人人添人人爽欧美一区卜| av国产精品久久久久影院| 在线观看免费高清a一片| 亚洲午夜理论影院| 高清欧美精品videossex| 一边摸一边抽搐一进一出视频| 午夜激情av网站| 一级黄色大片毛片| 热re99久久精品国产66热6| 午夜视频精品福利| 嫩草影视91久久| 国产在线一区二区三区精| 精品亚洲成国产av| 啦啦啦免费观看视频1| 国产蜜桃级精品一区二区三区 | 老司机福利观看| 国产片内射在线| 精品少妇久久久久久888优播| 久久国产亚洲av麻豆专区| 精品国产国语对白av| 欧美日韩瑟瑟在线播放| 欧美成人午夜精品| 免费在线观看影片大全网站| 丝瓜视频免费看黄片| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品aⅴ一区二区三区四区| 午夜福利,免费看| 国产精品av久久久久免费| 精品熟女少妇八av免费久了| 桃红色精品国产亚洲av| 欧美日韩一级在线毛片| 十分钟在线观看高清视频www| av网站在线播放免费| 丰满迷人的少妇在线观看| 在线观看免费高清a一片| 超碰成人久久| 午夜两性在线视频| 好男人电影高清在线观看| 丁香欧美五月| 欧美老熟妇乱子伦牲交| av有码第一页| 亚洲国产精品一区二区三区在线| 亚洲精品自拍成人| av福利片在线| 精品欧美一区二区三区在线| av天堂久久9| 亚洲av片天天在线观看| 9热在线视频观看99| 欧美在线黄色| 99精品欧美一区二区三区四区| 色94色欧美一区二区| 国产成人精品无人区| tocl精华| aaaaa片日本免费| 丰满人妻熟妇乱又伦精品不卡| 亚洲 欧美一区二区三区| 国产免费男女视频| 丰满饥渴人妻一区二区三| 国产精品影院久久| 久久久精品免费免费高清| 19禁男女啪啪无遮挡网站| 亚洲国产精品一区二区三区在线| 国产精品成人在线| 国产97色在线日韩免费| 国产片内射在线| 看黄色毛片网站| av免费在线观看网站| 久久久久久人人人人人| 国产成人系列免费观看| 亚洲av美国av| 精品国产乱码久久久久久男人| 欧美亚洲 丝袜 人妻 在线| 精品一区二区三区四区五区乱码| 日韩欧美在线二视频 | 亚洲视频免费观看视频| 日韩熟女老妇一区二区性免费视频| 国产精品久久久久久人妻精品电影| 午夜精品久久久久久毛片777| 婷婷丁香在线五月| 亚洲 欧美一区二区三区| 久久久久国产精品人妻aⅴ院 | 在线观看午夜福利视频| 老司机午夜十八禁免费视频| 久久ye,这里只有精品| 91在线观看av| 一区二区三区激情视频| 精品久久久久久,| 日韩免费av在线播放| 午夜福利一区二区在线看| 国产精品.久久久| 国产精品 国内视频| 很黄的视频免费| 国产日韩一区二区三区精品不卡| 美国免费a级毛片| 国产无遮挡羞羞视频在线观看| 波多野结衣av一区二区av| 99国产精品免费福利视频| 天堂中文最新版在线下载| 波多野结衣一区麻豆| 国产男女内射视频| 色播在线永久视频| 精品少妇一区二区三区视频日本电影| 婷婷精品国产亚洲av在线 | 亚洲avbb在线观看| 在线国产一区二区在线| 成人特级黄色片久久久久久久| 欧美av亚洲av综合av国产av| 在线天堂中文资源库| 一级,二级,三级黄色视频| 精品无人区乱码1区二区| 黑人欧美特级aaaaaa片| 色婷婷久久久亚洲欧美| 欧美+亚洲+日韩+国产| 欧美日本中文国产一区发布| 日韩有码中文字幕| 中文字幕av电影在线播放| 九色亚洲精品在线播放| 亚洲伊人色综图| 深夜精品福利| 久久久久久久久久久久大奶| 很黄的视频免费| 免费在线观看日本一区| 精品福利观看| 新久久久久国产一级毛片| 老汉色av国产亚洲站长工具| 一级,二级,三级黄色视频| 少妇的丰满在线观看| 国产1区2区3区精品| 国产精品亚洲一级av第二区| 亚洲欧美色中文字幕在线| 日韩熟女老妇一区二区性免费视频| 欧美日韩成人在线一区二区| 欧美日韩亚洲国产一区二区在线观看 | 法律面前人人平等表现在哪些方面| 国产成人精品久久二区二区免费| 精品一区二区三区视频在线观看免费 | x7x7x7水蜜桃| 欧美最黄视频在线播放免费 | svipshipincom国产片| 欧美 亚洲 国产 日韩一| 99国产极品粉嫩在线观看| 啦啦啦 在线观看视频| 男人的好看免费观看在线视频 | 日韩欧美一区视频在线观看| 国产有黄有色有爽视频| 视频区图区小说| 看免费av毛片| 丰满饥渴人妻一区二区三| 黄色女人牲交| 人成视频在线观看免费观看| 亚洲精品乱久久久久久| 亚洲欧美一区二区三区久久| 国产欧美日韩一区二区三区在线| 人人妻人人爽人人添夜夜欢视频| 视频在线观看一区二区三区| 欧美另类亚洲清纯唯美| 免费在线观看视频国产中文字幕亚洲| 国产深夜福利视频在线观看| 国产精品九九99| 免费观看a级毛片全部| 亚洲情色 制服丝袜| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产看品久久| 亚洲欧美一区二区三区黑人| 欧美丝袜亚洲另类 | 天堂中文最新版在线下载| 亚洲欧美激情综合另类| 黄色毛片三级朝国网站| 少妇粗大呻吟视频| 男女床上黄色一级片免费看| 又大又爽又粗| 正在播放国产对白刺激| 国产精品免费视频内射| 在线观看免费视频网站a站| 人成视频在线观看免费观看| 成人三级做爰电影| 老司机亚洲免费影院| 在线观看一区二区三区激情| 欧美日韩乱码在线| 18禁裸乳无遮挡动漫免费视频| 天天躁日日躁夜夜躁夜夜| av中文乱码字幕在线| 国产男女超爽视频在线观看| 美女 人体艺术 gogo| 午夜日韩欧美国产| av天堂在线播放| 亚洲九九香蕉| 国产97色在线日韩免费| 一本大道久久a久久精品| 日韩成人在线观看一区二区三区| 日韩欧美免费精品| 女人被狂操c到高潮| av不卡在线播放| 亚洲国产看品久久| 国产精品久久久久久人妻精品电影| 国产亚洲欧美精品永久| 老司机在亚洲福利影院| av视频免费观看在线观看| 亚洲一区二区三区欧美精品| 男男h啪啪无遮挡| 亚洲,欧美精品.| 一区二区三区精品91| 午夜成年电影在线免费观看| 欧美性长视频在线观看| 精品一区二区三卡| 日韩人妻精品一区2区三区| 国产人伦9x9x在线观看| 精品国产亚洲在线| 捣出白浆h1v1| 久久久精品国产亚洲av高清涩受| 9191精品国产免费久久| 免费黄频网站在线观看国产| 日本a在线网址| 国产成人啪精品午夜网站| av免费在线观看网站| a级片在线免费高清观看视频| 成年人黄色毛片网站| 黄色毛片三级朝国网站| av天堂在线播放| 一夜夜www| 国产精品永久免费网站| 在线观看www视频免费| 亚洲精品粉嫩美女一区| 性少妇av在线| 不卡av一区二区三区| 久久久精品免费免费高清| 亚洲五月婷婷丁香| 国产激情欧美一区二区| av天堂在线播放| 亚洲一码二码三码区别大吗| 久久久久精品国产欧美久久久| 国产亚洲一区二区精品| 国产亚洲欧美在线一区二区| e午夜精品久久久久久久| 欧美精品啪啪一区二区三区| 亚洲片人在线观看| 男人的好看免费观看在线视频 | 亚洲熟妇熟女久久| aaaaa片日本免费| 一边摸一边抽搐一进一出视频| 国产一卡二卡三卡精品| 欧美日韩一级在线毛片| 亚洲一区二区三区不卡视频| 在线观看免费视频网站a站| 一区二区三区精品91| 99国产精品99久久久久| 国产成人系列免费观看| 在线观看午夜福利视频| 亚洲精品av麻豆狂野| 波多野结衣av一区二区av| 久久精品91无色码中文字幕| 亚洲精品久久午夜乱码| 国产男靠女视频免费网站| 啪啪无遮挡十八禁网站| 黄色怎么调成土黄色| 精品人妻1区二区| 高清黄色对白视频在线免费看| xxx96com| 涩涩av久久男人的天堂| 欧美精品人与动牲交sv欧美| 人妻丰满熟妇av一区二区三区 | 无限看片的www在线观看| 在线视频色国产色| 国产黄色免费在线视频| 国产成人免费无遮挡视频| 色在线成人网| 日韩大码丰满熟妇| 看黄色毛片网站| avwww免费| 中文字幕人妻熟女乱码| 国产精品乱码一区二三区的特点 | 伦理电影免费视频| 男人舔女人的私密视频| 男女下面插进去视频免费观看| 久久人人爽av亚洲精品天堂| 免费看a级黄色片| 99riav亚洲国产免费| 色播在线永久视频| 精品人妻在线不人妻| 免费在线观看黄色视频的| www日本在线高清视频| 99热只有精品国产| 最新在线观看一区二区三区| 亚洲欧美日韩另类电影网站| 一个人免费在线观看的高清视频| 午夜成年电影在线免费观看| 国产欧美日韩一区二区三区在线| 最近最新中文字幕大全免费视频| 精品一区二区三区av网在线观看| 色在线成人网| 欧美乱妇无乱码| 国产精品乱码一区二三区的特点 | 久久人人爽av亚洲精品天堂| 精品亚洲成a人片在线观看| 嫁个100分男人电影在线观看| 9191精品国产免费久久| 1024香蕉在线观看| 最新在线观看一区二区三区| 国产亚洲一区二区精品| 男女之事视频高清在线观看| 亚洲欧美日韩高清在线视频| 欧美在线黄色| 18禁黄网站禁片午夜丰满| 国产精品免费视频内射| 亚洲人成77777在线视频| 久久午夜综合久久蜜桃| 波多野结衣av一区二区av| tocl精华| a级毛片黄视频| 国产成人精品在线电影| 淫妇啪啪啪对白视频| 自拍欧美九色日韩亚洲蝌蚪91| 黑人操中国人逼视频| 一个人免费在线观看的高清视频| 91精品国产国语对白视频| 亚洲精品美女久久久久99蜜臀| 精品国产乱子伦一区二区三区| 人人妻人人爽人人添夜夜欢视频| 亚洲精品av麻豆狂野| 90打野战视频偷拍视频| www日本在线高清视频| 久久久久久久精品吃奶| 国产精品1区2区在线观看. | 激情在线观看视频在线高清 | 亚洲伊人色综图| 精品久久久久久电影网| 精品人妻1区二区| 欧美成人免费av一区二区三区 | 欧美另类亚洲清纯唯美| 免费在线观看完整版高清| 国产成人免费无遮挡视频| 精品免费久久久久久久清纯 | 三级毛片av免费| 国产精品久久久av美女十八| 丰满迷人的少妇在线观看| 午夜福利在线免费观看网站| 亚洲精品一卡2卡三卡4卡5卡| 国产黄色免费在线视频| www.自偷自拍.com| 久久久国产一区二区| 91精品三级在线观看| 美国免费a级毛片| 国产黄色免费在线视频| 国产日韩欧美亚洲二区| av天堂久久9| 午夜福利在线免费观看网站| 亚洲精品成人av观看孕妇| 国产成人av教育| 在线观看免费视频网站a站| 精品国内亚洲2022精品成人 | netflix在线观看网站| a级毛片在线看网站| 精品电影一区二区在线| 欧美另类亚洲清纯唯美| 成人免费观看视频高清| 欧美精品亚洲一区二区| 三级毛片av免费| 他把我摸到了高潮在线观看| 亚洲av片天天在线观看| 亚洲一区二区三区欧美精品| 亚洲伊人色综图| 久久人妻av系列| 久久中文字幕人妻熟女| 中文亚洲av片在线观看爽 | 另类亚洲欧美激情| 色综合婷婷激情| 午夜福利在线观看吧| 纯流量卡能插随身wifi吗| 国产亚洲欧美精品永久| 欧美激情极品国产一区二区三区| 在线十欧美十亚洲十日本专区| 欧美日韩视频精品一区| 国产xxxxx性猛交| 国产99久久九九免费精品| 亚洲av熟女| 国产高清视频在线播放一区| 欧美另类亚洲清纯唯美| 一个人免费在线观看的高清视频| 亚洲三区欧美一区| 女人被狂操c到高潮| 人人妻人人爽人人添夜夜欢视频| 好看av亚洲va欧美ⅴa在| 久久天躁狠狠躁夜夜2o2o| 国产精品影院久久| 欧美大码av| 香蕉丝袜av| 国产无遮挡羞羞视频在线观看| 久久草成人影院| 99在线人妻在线中文字幕 | 视频区图区小说| 国产91精品成人一区二区三区| 又黄又粗又硬又大视频| 男人操女人黄网站| 国产淫语在线视频| 精品高清国产在线一区| 久久久久国内视频| 亚洲美女黄片视频| 夜夜躁狠狠躁天天躁| 满18在线观看网站| 男人舔女人的私密视频| av一本久久久久| 俄罗斯特黄特色一大片| 国产99白浆流出| 色综合欧美亚洲国产小说| 成人亚洲精品一区在线观看| 激情视频va一区二区三区| 99久久综合精品五月天人人| 国产精品九九99| 麻豆成人av在线观看| 欧美黄色淫秽网站| 捣出白浆h1v1| 水蜜桃什么品种好| 色婷婷av一区二区三区视频| 深夜精品福利| 激情在线观看视频在线高清 | 欧美午夜高清在线| 无限看片的www在线观看| 久9热在线精品视频| 国产主播在线观看一区二区|