• <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
    日本vs欧美在线观看视频| 国产成人一区二区在线| 亚洲第一区二区三区不卡| 国产精品一国产av| 校园人妻丝袜中文字幕| 国产成人欧美| a级片在线免费高清观看视频| 大话2 男鬼变身卡| 啦啦啦在线观看免费高清www| 免费观看性生交大片5| 下体分泌物呈黄色| 日本欧美视频一区| 国产av码专区亚洲av| 视频在线观看一区二区三区| 国产永久视频网站| 色吧在线观看| 一本大道久久a久久精品| 亚洲精品中文字幕在线视频| 国产午夜精品一二区理论片| 香蕉丝袜av| 日本与韩国留学比较| 高清av免费在线| 国产在视频线精品| 乱码一卡2卡4卡精品| 精品酒店卫生间| 2022亚洲国产成人精品| 有码 亚洲区| av电影中文网址| 久久久久久人妻| 人人妻人人爽人人添夜夜欢视频| av播播在线观看一区| 精品一区在线观看国产| 久久人人爽人人片av| 乱码一卡2卡4卡精品| av有码第一页| 成人国产麻豆网| 久久99热这里只频精品6学生| 51国产日韩欧美| 老熟女久久久| 免费人成在线观看视频色| 国产精品.久久久| 国产精品成人在线| 美女xxoo啪啪120秒动态图| 夜夜骑夜夜射夜夜干| 亚洲精品国产av蜜桃| 亚洲国产精品999| 亚洲一码二码三码区别大吗| 久久久国产一区二区| 欧美日韩国产mv在线观看视频| 黄色 视频免费看| 久热这里只有精品99| 肉色欧美久久久久久久蜜桃| 日韩三级伦理在线观看| 亚洲欧美清纯卡通| 国产色爽女视频免费观看| 精品国产露脸久久av麻豆| 不卡视频在线观看欧美| 精品国产一区二区久久| 日韩免费高清中文字幕av| 十八禁高潮呻吟视频| 免费av不卡在线播放| 一区二区av电影网| 亚洲国产最新在线播放| 夜夜爽夜夜爽视频| 亚洲精品自拍成人| 欧美精品av麻豆av| 香蕉精品网在线| 国产一区二区三区综合在线观看 | 最近手机中文字幕大全| 亚洲少妇的诱惑av| 日本黄大片高清| 国产一区亚洲一区在线观看| 久久久久国产精品人妻一区二区| 国产精品人妻久久久久久| 1024视频免费在线观看| 久久ye,这里只有精品| 97精品久久久久久久久久精品| 一区二区日韩欧美中文字幕 | av在线播放精品| 国产精品久久久久久精品古装| 日韩成人伦理影院| 91国产中文字幕| 美女xxoo啪啪120秒动态图| 国产精品国产三级国产专区5o| 午夜免费男女啪啪视频观看| 国产黄色免费在线视频| 久久久久人妻精品一区果冻| 欧美 日韩 精品 国产| 日本黄色日本黄色录像| 婷婷色综合大香蕉| 一二三四中文在线观看免费高清| 乱码一卡2卡4卡精品| 99国产综合亚洲精品| 日韩不卡一区二区三区视频在线| 好男人视频免费观看在线| 久久精品国产a三级三级三级| 久久女婷五月综合色啪小说| 又粗又硬又长又爽又黄的视频| 国内精品宾馆在线| 99热全是精品| 黄色配什么色好看| 日日撸夜夜添| 国产精品久久久久久精品电影小说| 国产又爽黄色视频| 汤姆久久久久久久影院中文字幕| 在线观看国产h片| 国产精品一二三区在线看| 欧美+日韩+精品| 欧美日韩视频精品一区| 亚洲欧美日韩卡通动漫| 自线自在国产av| 亚洲人成77777在线视频| 99热网站在线观看| 在线免费观看不下载黄p国产| 午夜福利,免费看| 波多野结衣一区麻豆| 免费观看无遮挡的男女| 亚洲国产毛片av蜜桃av| 欧美精品一区二区大全| 黄色怎么调成土黄色| 18禁观看日本| 最近最新中文字幕大全免费视频 | 精品少妇黑人巨大在线播放| 亚洲av男天堂| 免费黄色在线免费观看| 国产精品一区二区在线观看99| 久久久久视频综合| 国产1区2区3区精品| 国产精品久久久久成人av| 国产日韩欧美在线精品| 91国产中文字幕| 美女福利国产在线| 美女中出高潮动态图| 男女边摸边吃奶| 黑人欧美特级aaaaaa片| 国产日韩欧美亚洲二区| 午夜日本视频在线| 亚洲精品av麻豆狂野| 精品亚洲乱码少妇综合久久| 人成视频在线观看免费观看| 搡老乐熟女国产| 国产一级毛片在线| 国产成人精品久久久久久| 黑人猛操日本美女一级片| 久久精品久久久久久久性| 久久久久国产网址| 色哟哟·www| 丰满饥渴人妻一区二区三| 黄色配什么色好看| 777米奇影视久久| 三级国产精品片| 亚洲国产日韩一区二区| 午夜激情久久久久久久| 晚上一个人看的免费电影| 欧美最新免费一区二区三区| 国产视频首页在线观看| 久久久国产欧美日韩av| 日本与韩国留学比较| 美女内射精品一级片tv| 97超碰精品成人国产| 欧美日韩精品成人综合77777| 午夜福利影视在线免费观看| 91成人精品电影| 精品国产露脸久久av麻豆| 少妇的逼好多水| 黄网站色视频无遮挡免费观看| 大香蕉久久网| 在线观看三级黄色| 欧美日韩亚洲高清精品| 欧美精品一区二区大全| 亚洲婷婷狠狠爱综合网| 亚洲精品久久成人aⅴ小说| 在线观看免费日韩欧美大片| 中文字幕精品免费在线观看视频 | 国精品久久久久久国模美| 又黄又爽又刺激的免费视频.| 最新的欧美精品一区二区| 亚洲,一卡二卡三卡| 精品亚洲成国产av| 国产熟女午夜一区二区三区| 亚洲国产色片| 欧美+日韩+精品| 亚洲一级一片aⅴ在线观看| 久久精品夜色国产| 九九在线视频观看精品| 天美传媒精品一区二区| 精品视频人人做人人爽| 日本vs欧美在线观看视频| 免费黄色在线免费观看| 久久人妻熟女aⅴ| 丰满迷人的少妇在线观看| 成年av动漫网址| 国产免费一级a男人的天堂| 久久99热这里只频精品6学生| 国产高清三级在线| 人人妻人人澡人人看| 午夜激情av网站| 欧美日本中文国产一区发布| 又粗又硬又长又爽又黄的视频| 国产老妇伦熟女老妇高清| 国产伦理片在线播放av一区| 五月玫瑰六月丁香| 中文字幕免费在线视频6| 欧美xxⅹ黑人| 在线观看免费高清a一片| 精品一区二区三卡| 日本wwww免费看| 97人妻天天添夜夜摸| 捣出白浆h1v1| www.熟女人妻精品国产 | 女人精品久久久久毛片| 女的被弄到高潮叫床怎么办| 1024视频免费在线观看| 在现免费观看毛片| 欧美激情极品国产一区二区三区 | 亚洲成av片中文字幕在线观看 | 男女高潮啪啪啪动态图| 日本91视频免费播放| freevideosex欧美| 亚洲av免费高清在线观看| 伦精品一区二区三区| av又黄又爽大尺度在线免费看| 在线观看www视频免费| 亚洲第一区二区三区不卡| 黄色怎么调成土黄色| 夜夜爽夜夜爽视频| 日日撸夜夜添| 久久久久久人妻| 亚洲一区二区三区欧美精品| av不卡在线播放| 国产成人一区二区在线| 亚洲内射少妇av| 亚洲精品国产av蜜桃| 日韩av在线免费看完整版不卡| videossex国产| 欧美激情 高清一区二区三区| 日本vs欧美在线观看视频| 建设人人有责人人尽责人人享有的| 99国产综合亚洲精品| 黄片播放在线免费| 国产免费又黄又爽又色| av福利片在线| 欧美成人午夜精品| 国产永久视频网站| 久久国内精品自在自线图片| 欧美最新免费一区二区三区| 综合色丁香网| 精品人妻在线不人妻| 亚洲国产精品专区欧美| 女人被躁到高潮嗷嗷叫费观| 国产精品成人在线| 少妇熟女欧美另类| 国产精品 国内视频| 午夜老司机福利剧场| 亚洲久久久国产精品| 国产一区亚洲一区在线观看| 极品少妇高潮喷水抽搐| 亚洲av成人精品一二三区| 人妻系列 视频| 香蕉丝袜av| 亚洲,一卡二卡三卡| 另类亚洲欧美激情| 搡女人真爽免费视频火全软件| 性色avwww在线观看| 日韩欧美一区视频在线观看| 亚洲国产欧美在线一区| 午夜福利视频在线观看免费| 少妇人妻久久综合中文| 韩国精品一区二区三区 | 久久精品国产a三级三级三级| 只有这里有精品99| 久久人人爽人人片av| 国产成人精品久久久久久| 在线观看免费高清a一片| 熟女av电影| 午夜激情久久久久久久| 亚洲第一av免费看| 爱豆传媒免费全集在线观看| 丁香六月天网| 一区二区av电影网| 九九爱精品视频在线观看| 国产精品国产三级专区第一集| 高清欧美精品videossex| 国产精品.久久久| 国产片内射在线| 亚洲精品国产av蜜桃| 日本黄大片高清| 亚洲高清免费不卡视频| 日韩成人伦理影院| 久久影院123| 国产xxxxx性猛交| 美女国产视频在线观看| 女人被躁到高潮嗷嗷叫费观| av免费观看日本| 最近的中文字幕免费完整| 你懂的网址亚洲精品在线观看| 91在线精品国自产拍蜜月| av视频免费观看在线观看| av又黄又爽大尺度在线免费看| 精品亚洲乱码少妇综合久久| 菩萨蛮人人尽说江南好唐韦庄| 妹子高潮喷水视频| av在线观看视频网站免费| 在线 av 中文字幕| 日本-黄色视频高清免费观看| 免费观看无遮挡的男女| 五月开心婷婷网| 日韩中文字幕视频在线看片| 三级国产精品片| 伊人亚洲综合成人网| av电影中文网址| 久久久国产精品麻豆| xxxhd国产人妻xxx| 欧美激情国产日韩精品一区| 久久热在线av| 考比视频在线观看| 夫妻午夜视频| 午夜福利乱码中文字幕| 国产亚洲欧美精品永久| 国产免费又黄又爽又色| 大片电影免费在线观看免费| 国产日韩欧美在线精品| 国产成人精品婷婷| 日韩视频在线欧美| av又黄又爽大尺度在线免费看| av在线app专区| 又大又黄又爽视频免费| 中文字幕制服av| av播播在线观看一区| 亚洲欧美日韩另类电影网站| 一级,二级,三级黄色视频| 国产精品嫩草影院av在线观看| 全区人妻精品视频| av不卡在线播放| 国产av码专区亚洲av| 好男人视频免费观看在线| 国产精品.久久久| av在线播放精品| 久久精品国产亚洲av涩爱| 久久热在线av| 不卡视频在线观看欧美| 国产av国产精品国产| 国产伦理片在线播放av一区| 欧美日韩视频高清一区二区三区二| 宅男免费午夜| 制服人妻中文乱码| 欧美+日韩+精品| 亚洲情色 制服丝袜| 欧美老熟妇乱子伦牲交| 成人18禁高潮啪啪吃奶动态图| 免费看不卡的av| 伦理电影免费视频| 国产精品人妻久久久久久| 黑人猛操日本美女一级片| h视频一区二区三区| 亚洲内射少妇av| 日韩电影二区| 一级a做视频免费观看| 国产免费又黄又爽又色| 国产亚洲午夜精品一区二区久久| 午夜免费观看性视频| 国产精品三级大全| 欧美变态另类bdsm刘玥| 欧美性感艳星| 亚洲国产日韩一区二区| 亚洲成人手机| 大陆偷拍与自拍| 一级毛片我不卡| 熟女电影av网| av不卡在线播放| 亚洲av日韩在线播放| 色婷婷久久久亚洲欧美| a级毛片在线看网站| 欧美xxxx性猛交bbbb| 男女边摸边吃奶| 免费黄色在线免费观看| 国产又爽黄色视频| 免费久久久久久久精品成人欧美视频 | 欧美成人午夜精品| 如日韩欧美国产精品一区二区三区| 黄网站色视频无遮挡免费观看| 国产黄频视频在线观看| 亚洲国产成人一精品久久久| 亚洲国产精品专区欧美| 高清毛片免费看| 少妇人妻久久综合中文| 日本欧美视频一区| 免费久久久久久久精品成人欧美视频 | 亚洲欧美精品自产自拍| 久久久久网色| 日韩一本色道免费dvd| 99热全是精品| 美女主播在线视频| 丝瓜视频免费看黄片| 国产欧美日韩综合在线一区二区| 亚洲图色成人| 日韩大片免费观看网站| 精品一区二区三卡| 国产高清国产精品国产三级| 人人妻人人爽人人添夜夜欢视频| 黄色毛片三级朝国网站| 亚洲人成77777在线视频| av电影中文网址| 亚洲精品美女久久av网站| 丰满乱子伦码专区| 国产深夜福利视频在线观看| 成人黄色视频免费在线看| 国产成人精品久久久久久| 十分钟在线观看高清视频www| 最近2019中文字幕mv第一页| 成人毛片60女人毛片免费| 99久久精品国产国产毛片| 国产精品国产av在线观看| 免费av中文字幕在线| 国产色婷婷99| 成人国产av品久久久| 啦啦啦啦在线视频资源| 亚洲av中文av极速乱| 欧美国产精品va在线观看不卡| 中文字幕人妻丝袜制服| 日日摸夜夜添夜夜爱| 成年人免费黄色播放视频| 搡老乐熟女国产| 国产一区有黄有色的免费视频| 亚洲国产精品专区欧美| 久久精品国产鲁丝片午夜精品| 欧美日韩视频精品一区| 欧美丝袜亚洲另类| 如日韩欧美国产精品一区二区三区| 狂野欧美激情性bbbbbb| 纵有疾风起免费观看全集完整版| 久久午夜福利片| 波野结衣二区三区在线| 国产成人av激情在线播放| 欧美老熟妇乱子伦牲交| 中国国产av一级| 欧美丝袜亚洲另类| 免费人妻精品一区二区三区视频| 日韩 亚洲 欧美在线| 熟女电影av网| 成人无遮挡网站| 国产熟女欧美一区二区| 美女内射精品一级片tv| 久久精品人人爽人人爽视色| 午夜免费鲁丝| 尾随美女入室| 激情五月婷婷亚洲| 一级爰片在线观看| 国产一区二区在线观看av| 妹子高潮喷水视频| 亚洲激情五月婷婷啪啪| 欧美国产精品va在线观看不卡| 国产白丝娇喘喷水9色精品| 久久ye,这里只有精品| 午夜免费观看性视频| 精品午夜福利在线看| 在线观看人妻少妇| 国产极品粉嫩免费观看在线| 久久精品久久久久久久性| 22中文网久久字幕| 日本欧美国产在线视频| 亚洲第一区二区三区不卡| 亚洲国产最新在线播放| 久久久国产精品麻豆| 亚洲,一卡二卡三卡| 在线精品无人区一区二区三| 国产精品一区二区在线不卡| 男女无遮挡免费网站观看| 国产精品国产三级国产av玫瑰| 一边摸一边做爽爽视频免费| 久久久国产精品麻豆| 在线观看美女被高潮喷水网站| 国产亚洲av片在线观看秒播厂| 免费高清在线观看视频在线观看| 高清欧美精品videossex| 精品国产露脸久久av麻豆| 国产免费一区二区三区四区乱码| 久久久久久久亚洲中文字幕| 久久久久久人人人人人| 最新的欧美精品一区二区| 大码成人一级视频| 在线观看www视频免费| 亚洲成av片中文字幕在线观看 | 肉色欧美久久久久久久蜜桃| 男男h啪啪无遮挡| 国产毛片在线视频| 黑人猛操日本美女一级片| 国产一区有黄有色的免费视频| 26uuu在线亚洲综合色| 激情视频va一区二区三区| 2022亚洲国产成人精品| 久久ye,这里只有精品| 不卡视频在线观看欧美| 国产亚洲精品久久久com| 亚洲激情五月婷婷啪啪| 熟妇人妻不卡中文字幕| 亚洲综合色网址| 亚洲欧美成人精品一区二区| 啦啦啦视频在线资源免费观看| 亚洲欧美精品自产自拍| 另类精品久久| 国产成人免费观看mmmm| av视频免费观看在线观看| 高清不卡的av网站| 天天操日日干夜夜撸| 国产精品久久久久久久久免| xxx大片免费视频| 亚洲欧美成人精品一区二区| videossex国产| 一级,二级,三级黄色视频| 在线亚洲精品国产二区图片欧美| 国产成人精品婷婷| 亚洲国产日韩一区二区| 在线 av 中文字幕| 中文字幕人妻熟女乱码| 国产淫语在线视频| 成人18禁高潮啪啪吃奶动态图| 尾随美女入室| 一级a做视频免费观看| 久久这里有精品视频免费| 丰满乱子伦码专区| 人妻一区二区av| 免费观看性生交大片5| 亚洲欧美一区二区三区黑人 | 亚洲成av片中文字幕在线观看 | 天堂俺去俺来也www色官网| 久久久久久伊人网av| 蜜桃在线观看..| av卡一久久| 91午夜精品亚洲一区二区三区| 九九爱精品视频在线观看| 色哟哟·www| 国产精品.久久久| 国产成人aa在线观看| 最黄视频免费看| 国产xxxxx性猛交| 久久精品久久精品一区二区三区| 插逼视频在线观看| 日韩一区二区三区影片| 日韩大片免费观看网站| 亚洲成av片中文字幕在线观看 | 久久精品久久久久久久性| 18禁裸乳无遮挡动漫免费视频| 国产成人一区二区在线| 热99国产精品久久久久久7| 日本色播在线视频| 国产xxxxx性猛交| 国产av国产精品国产| 交换朋友夫妻互换小说| 九九爱精品视频在线观看| av一本久久久久| 国产成人精品在线电影| 久久精品夜色国产| av片东京热男人的天堂| 亚洲精品一二三| 少妇猛男粗大的猛烈进出视频| 久久免费观看电影| 黑人巨大精品欧美一区二区蜜桃 | 精品酒店卫生间| 国产精品99久久99久久久不卡 | 美女主播在线视频| 久久精品aⅴ一区二区三区四区 | 亚洲欧美一区二区三区国产| 日韩制服骚丝袜av| 久久av网站| 婷婷色综合www| 哪个播放器可以免费观看大片| videos熟女内射| 视频中文字幕在线观看| 亚洲欧美一区二区三区黑人 | 国产成人精品在线电影| 三级国产精品片| 亚洲精品色激情综合| 亚洲精品美女久久久久99蜜臀 | 久久精品国产自在天天线| 久热这里只有精品99| 欧美日韩国产mv在线观看视频| 街头女战士在线观看网站| 巨乳人妻的诱惑在线观看| 免费观看性生交大片5| 22中文网久久字幕| 极品少妇高潮喷水抽搐| 天美传媒精品一区二区| 母亲3免费完整高清在线观看 | 一区二区日韩欧美中文字幕 | 国产在线免费精品| 在线免费观看不下载黄p国产| 中国三级夫妇交换| 日韩电影二区| 制服人妻中文乱码| 成人亚洲精品一区在线观看| 在线免费观看不下载黄p国产| 色吧在线观看| 免费大片黄手机在线观看| 黄色视频在线播放观看不卡| 精品国产一区二区久久| 人妻一区二区av| 久久久久精品性色| 久久这里只有精品19| 久久av网站| 久久久久网色| 日本vs欧美在线观看视频| 99国产精品免费福利视频| 97人妻天天添夜夜摸| 成人国产av品久久久| 午夜福利视频在线观看免费| 免费观看性生交大片5| 亚洲精品视频女| 久久久国产欧美日韩av| 中文字幕av电影在线播放| 91午夜精品亚洲一区二区三区| tube8黄色片| 在线观看免费高清a一片| 国产黄频视频在线观看| 久久婷婷青草|