• <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
    999久久久国产精品视频| 在线观看免费日韩欧美大片| videosex国产| 国产成人精品久久二区二区91| 超色免费av| 操出白浆在线播放| 免费久久久久久久精品成人欧美视频| 日日爽夜夜爽网站| 在线观看66精品国产| 99久久人妻综合| 成年动漫av网址| 黄色视频不卡| 桃红色精品国产亚洲av| 成年女人毛片免费观看观看9 | 黄色丝袜av网址大全| 99国产精品99久久久久| xxxhd国产人妻xxx| 午夜激情av网站| 大型黄色视频在线免费观看| 嫩草影视91久久| 啦啦啦免费观看视频1| 日韩欧美三级三区| 亚洲av熟女| 又大又爽又粗| 老司机午夜十八禁免费视频| 成人三级做爰电影| 国产精品电影一区二区三区 | 国产亚洲欧美在线一区二区| 丰满的人妻完整版| 亚洲精品久久成人aⅴ小说| 99re6热这里在线精品视频| 亚洲片人在线观看| 另类亚洲欧美激情| 18禁国产床啪视频网站| 精品人妻熟女毛片av久久网站| 色94色欧美一区二区| 中亚洲国语对白在线视频| 亚洲精品自拍成人| 午夜精品国产一区二区电影| 成人国产一区最新在线观看| 男女午夜视频在线观看| 美女午夜性视频免费| 在线免费观看的www视频| 性色av乱码一区二区三区2| 黄片大片在线免费观看| 欧美日韩成人在线一区二区| 午夜老司机福利片| 免费在线观看亚洲国产| 狂野欧美激情性xxxx| 国产男女内射视频| 搡老乐熟女国产| 黑人欧美特级aaaaaa片| 免费在线观看日本一区| 亚洲av片天天在线观看| 欧美 亚洲 国产 日韩一| 一级a爱视频在线免费观看| 亚洲专区国产一区二区| 欧美人与性动交α欧美软件| 50天的宝宝边吃奶边哭怎么回事| 怎么达到女性高潮| 国产成人啪精品午夜网站| 免费在线观看日本一区| 高潮久久久久久久久久久不卡| 黑人欧美特级aaaaaa片| 丁香六月欧美| 国产精品99久久99久久久不卡| 国产深夜福利视频在线观看| 午夜福利乱码中文字幕| 成在线人永久免费视频| 久久性视频一级片| 欧美成人免费av一区二区三区 | 下体分泌物呈黄色| 乱人伦中国视频| 精品乱码久久久久久99久播| 久久精品国产亚洲av高清一级| 交换朋友夫妻互换小说| 最新的欧美精品一区二区| 午夜精品在线福利| 免费在线观看亚洲国产| 精品欧美一区二区三区在线| 在线播放国产精品三级| 91大片在线观看| 国产成人免费无遮挡视频| 午夜91福利影院| 国产亚洲精品一区二区www | 精品国产乱码久久久久久男人| 男女下面插进去视频免费观看| 久久午夜亚洲精品久久| 9色porny在线观看| 黑人巨大精品欧美一区二区蜜桃| 黄色视频,在线免费观看| 少妇的丰满在线观看| 女性被躁到高潮视频| avwww免费| 视频在线观看一区二区三区| 中文字幕av电影在线播放| 乱人伦中国视频| 狠狠婷婷综合久久久久久88av| 大码成人一级视频| 一边摸一边做爽爽视频免费| 久久精品人人爽人人爽视色| 久久久久久亚洲精品国产蜜桃av| 一级黄色大片毛片| 国产精品.久久久| 免费在线观看亚洲国产| 欧美国产精品va在线观看不卡| 国产高清激情床上av| 超色免费av| 一级毛片女人18水好多| 欧美+亚洲+日韩+国产| 亚洲性夜色夜夜综合| aaaaa片日本免费| 一本大道久久a久久精品| 国产不卡一卡二| 中亚洲国语对白在线视频| 亚洲情色 制服丝袜| 飞空精品影院首页| 亚洲精品中文字幕一二三四区| 午夜成年电影在线免费观看| 中文字幕制服av| 成在线人永久免费视频| 制服人妻中文乱码| 精品人妻1区二区| 天堂中文最新版在线下载| 免费久久久久久久精品成人欧美视频| 国产精品亚洲一级av第二区| 老司机影院毛片| 欧美一级毛片孕妇| 一区二区三区激情视频| 老司机影院毛片| 一级片'在线观看视频| 亚洲欧美日韩另类电影网站| videosex国产| www.熟女人妻精品国产| 好男人电影高清在线观看| 俄罗斯特黄特色一大片| 亚洲情色 制服丝袜| 嫩草影视91久久| 日本a在线网址| 国产精品 欧美亚洲| 亚洲精品国产区一区二| 男女之事视频高清在线观看| 国产成人av激情在线播放| 久久人人97超碰香蕉20202| 久久中文字幕一级| 亚洲欧洲精品一区二区精品久久久| 亚洲精品国产区一区二| 99国产精品一区二区蜜桃av | 俄罗斯特黄特色一大片| 亚洲情色 制服丝袜| 丝瓜视频免费看黄片| 性色av乱码一区二区三区2| 久久久久精品国产欧美久久久| 亚洲国产欧美日韩在线播放| 一夜夜www| 国产片内射在线| 每晚都被弄得嗷嗷叫到高潮| 久久久久久久久久久久大奶| 真人做人爱边吃奶动态| 精品一区二区三卡| 欧美精品一区二区免费开放| av一本久久久久| 国产精品亚洲av一区麻豆| 91国产中文字幕| 嫁个100分男人电影在线观看| √禁漫天堂资源中文www| 欧美 亚洲 国产 日韩一| 婷婷精品国产亚洲av在线 | 亚洲欧美激情综合另类| 国产人伦9x9x在线观看| 美女国产高潮福利片在线看| 成人国产一区最新在线观看| 亚洲精品美女久久av网站| 香蕉国产在线看| 夜夜躁狠狠躁天天躁| 国产亚洲av高清不卡| 国产又爽黄色视频| 久久人妻福利社区极品人妻图片| 男人操女人黄网站| 国产精品一区二区免费欧美| 精品福利观看| 无遮挡黄片免费观看| 一个人免费在线观看的高清视频| 欧美色视频一区免费| 香蕉丝袜av| 国产亚洲欧美98| 91精品国产国语对白视频| 12—13女人毛片做爰片一| 啦啦啦 在线观看视频| 免费观看a级毛片全部| 狂野欧美激情性xxxx| 成在线人永久免费视频| 国产色视频综合| 成人影院久久| 久久久精品区二区三区| 老熟女久久久| 人妻 亚洲 视频| 国产一区二区三区在线臀色熟女 | 另类亚洲欧美激情| 免费高清在线观看日韩| 日韩欧美免费精品| 欧美精品人与动牲交sv欧美| 日韩有码中文字幕| 欧美日韩视频精品一区| 一级片免费观看大全| 操美女的视频在线观看| 极品人妻少妇av视频| 中文字幕另类日韩欧美亚洲嫩草| 天天添夜夜摸| 久久精品人人爽人人爽视色| 国产伦人伦偷精品视频| 亚洲人成电影免费在线| 美女 人体艺术 gogo| 很黄的视频免费| 国产精品免费大片| xxxhd国产人妻xxx| 欧美激情久久久久久爽电影 | 搡老熟女国产l中国老女人| 亚洲欧美日韩另类电影网站| 国产在线一区二区三区精| 99国产精品免费福利视频| 久久性视频一级片| 777米奇影视久久| 99久久99久久久精品蜜桃| 丰满饥渴人妻一区二区三| 日韩制服丝袜自拍偷拍| 免费高清在线观看日韩| 久久精品国产99精品国产亚洲性色 | 国产单亲对白刺激| 婷婷丁香在线五月| 美女扒开内裤让男人捅视频| 久久精品亚洲熟妇少妇任你| 色在线成人网| 久久久久久久午夜电影 | 最近最新免费中文字幕在线| 高清av免费在线| 国产精品1区2区在线观看. | 黄片大片在线免费观看| 国产精品一区二区在线不卡| 一边摸一边做爽爽视频免费| 日韩欧美三级三区| 国产熟女午夜一区二区三区| 真人做人爱边吃奶动态| 搡老岳熟女国产| 在线观看午夜福利视频| 亚洲av片天天在线观看| 午夜福利欧美成人| 欧美乱妇无乱码| 搡老熟女国产l中国老女人| 少妇被粗大的猛进出69影院| 久99久视频精品免费| 免费日韩欧美在线观看| 亚洲成人国产一区在线观看| 久久香蕉激情| 国产成人影院久久av| 国产成人精品久久二区二区免费| 波多野结衣av一区二区av| 国产精品 国内视频| 久久热在线av| 自线自在国产av| 夜夜爽天天搞| 精品国内亚洲2022精品成人 | 久久久水蜜桃国产精品网| 99热国产这里只有精品6| 中文字幕高清在线视频| 人人妻,人人澡人人爽秒播| 日本a在线网址| 国内久久婷婷六月综合欲色啪| 成人永久免费在线观看视频| 久久人妻av系列| 在线观看免费视频日本深夜| 日韩欧美一区二区三区在线观看 | 国产成人一区二区三区免费视频网站| 真人做人爱边吃奶动态| 久久久精品区二区三区| 亚洲精品国产精品久久久不卡| 18禁黄网站禁片午夜丰满| 脱女人内裤的视频| 国产精品综合久久久久久久免费 | 美女 人体艺术 gogo| 国产99久久九九免费精品| 亚洲精品中文字幕在线视频| 亚洲五月天丁香| 精品久久久精品久久久| 大陆偷拍与自拍| 黄色怎么调成土黄色| 80岁老熟妇乱子伦牲交| 午夜久久久在线观看| 精品久久久久久,| 人人妻人人澡人人看| 精品国产超薄肉色丝袜足j| 国产精品久久久久久人妻精品电影| 丝袜美足系列| 90打野战视频偷拍视频| 亚洲成人手机| 美女高潮到喷水免费观看| 两个人看的免费小视频| 欧美黄色片欧美黄色片| 欧美成狂野欧美在线观看| 亚洲自偷自拍图片 自拍| 人人妻人人澡人人爽人人夜夜| 国产单亲对白刺激| 在线观看www视频免费| av网站免费在线观看视频| 午夜免费鲁丝| 午夜视频精品福利| 少妇裸体淫交视频免费看高清 | 美女午夜性视频免费| 夜夜爽天天搞| 人人妻人人添人人爽欧美一区卜| 老熟女久久久| 交换朋友夫妻互换小说| 久久午夜综合久久蜜桃| 国产欧美日韩一区二区三区在线| 香蕉丝袜av| www.熟女人妻精品国产| 男人舔女人的私密视频| 免费女性裸体啪啪无遮挡网站| 欧美另类亚洲清纯唯美| 国产亚洲精品久久久久久毛片 | 91精品国产国语对白视频| 美女午夜性视频免费| a在线观看视频网站| 99久久99久久久精品蜜桃| 午夜福利欧美成人| 香蕉国产在线看| 午夜日韩欧美国产| 久久狼人影院| 一级作爱视频免费观看| 又紧又爽又黄一区二区| av超薄肉色丝袜交足视频| 色老头精品视频在线观看| 免费女性裸体啪啪无遮挡网站| 咕卡用的链子| 久久人妻福利社区极品人妻图片| 无限看片的www在线观看| 黑丝袜美女国产一区| 男男h啪啪无遮挡| 大香蕉久久网| 99国产精品99久久久久| 老司机影院毛片| 一区二区三区国产精品乱码| 一个人免费在线观看的高清视频| 亚洲欧美一区二区三区久久| 久久中文字幕一级| av不卡在线播放| 狠狠婷婷综合久久久久久88av| 韩国av一区二区三区四区| 国产免费现黄频在线看| 欧美日韩国产mv在线观看视频| 国产精品一区二区在线观看99| 亚洲色图综合在线观看| 成人免费观看视频高清| 在线观看免费高清a一片| 国产午夜精品久久久久久| 50天的宝宝边吃奶边哭怎么回事| 亚洲熟妇熟女久久| 国产成人精品久久二区二区免费| 国产av精品麻豆| 久久精品国产综合久久久| 免费高清在线观看日韩| 无限看片的www在线观看| 精品国产乱码久久久久久男人| av电影中文网址| 午夜福利在线免费观看网站| 精品久久久精品久久久| а√天堂www在线а√下载 | 欧美色视频一区免费| 美女视频免费永久观看网站| 欧美精品一区二区免费开放| 国产精品综合久久久久久久免费 | 午夜福利在线免费观看网站| 亚洲成人手机| 欧美激情高清一区二区三区| 国产精品久久久久成人av| 国内毛片毛片毛片毛片毛片| 超碰成人久久| 一级片免费观看大全| 国产精品免费视频内射| 青草久久国产| 老熟妇乱子伦视频在线观看| 99国产综合亚洲精品| 国产99久久九九免费精品| 国产一卡二卡三卡精品| 精品午夜福利视频在线观看一区| 少妇粗大呻吟视频| 亚洲精品中文字幕一二三四区| 脱女人内裤的视频| 欧美在线黄色| 三级毛片av免费| 欧美激情久久久久久爽电影 | 老汉色∧v一级毛片| 日韩人妻精品一区2区三区| 黄色毛片三级朝国网站| 国产一卡二卡三卡精品| 亚洲综合色网址| 国产日韩欧美亚洲二区| 黄片大片在线免费观看| 无限看片的www在线观看| 少妇被粗大的猛进出69影院| 久久这里只有精品19| 水蜜桃什么品种好| 亚洲中文av在线| 国产精品久久电影中文字幕 | 99精国产麻豆久久婷婷| 亚洲精品国产色婷婷电影| 精品亚洲成国产av| 男人操女人黄网站| av不卡在线播放| 久久香蕉国产精品| 国产在线精品亚洲第一网站| av线在线观看网站| 亚洲aⅴ乱码一区二区在线播放 | 91麻豆av在线| 亚洲少妇的诱惑av| 99热国产这里只有精品6| 动漫黄色视频在线观看| 母亲3免费完整高清在线观看| 欧美日韩国产mv在线观看视频| 久久久水蜜桃国产精品网| 我的亚洲天堂| 国产精品久久久人人做人人爽| 国产精品秋霞免费鲁丝片| av网站在线播放免费| 伊人久久大香线蕉亚洲五| 欧美成狂野欧美在线观看| 久久精品亚洲熟妇少妇任你| 欧美黑人欧美精品刺激| 免费久久久久久久精品成人欧美视频| 精品人妻1区二区| 国产国语露脸激情在线看| 欧美亚洲日本最大视频资源| 色精品久久人妻99蜜桃| 亚洲久久久国产精品| 成人18禁在线播放| 日韩成人在线观看一区二区三区| 欧美成狂野欧美在线观看| 777久久人妻少妇嫩草av网站| 嫩草影视91久久| 亚洲视频免费观看视频| 美女 人体艺术 gogo| 一区二区三区激情视频| 欧美+亚洲+日韩+国产| 狠狠婷婷综合久久久久久88av| 高潮久久久久久久久久久不卡| av欧美777| 18禁裸乳无遮挡动漫免费视频| 老熟女久久久| 在线观看日韩欧美| 久久国产精品人妻蜜桃| 99国产精品99久久久久| 午夜福利乱码中文字幕| 人人妻人人澡人人看| 亚洲五月色婷婷综合| 免费女性裸体啪啪无遮挡网站| 淫妇啪啪啪对白视频| 777米奇影视久久| 丝瓜视频免费看黄片| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成人精品久久二区二区免费| 欧美黄色片欧美黄色片| 啦啦啦视频在线资源免费观看| 国内毛片毛片毛片毛片毛片| 亚洲欧美精品综合一区二区三区| 欧美人与性动交α欧美精品济南到| 日本五十路高清| 国产男女内射视频| 捣出白浆h1v1| 超碰成人久久| 久久精品亚洲av国产电影网| 日本五十路高清| 法律面前人人平等表现在哪些方面| 国产精品二区激情视频| 女人爽到高潮嗷嗷叫在线视频| 欧美精品人与动牲交sv欧美| 欧美 日韩 精品 国产| 十分钟在线观看高清视频www| 久久精品aⅴ一区二区三区四区| 午夜免费鲁丝| 午夜亚洲福利在线播放| 亚洲成人免费电影在线观看| 免费在线观看亚洲国产| 午夜福利一区二区在线看| 精品国内亚洲2022精品成人 | 大码成人一级视频| 美女高潮喷水抽搐中文字幕| 精品久久久精品久久久| 一边摸一边抽搐一进一小说 | 高清av免费在线| 日本精品一区二区三区蜜桃| 天天躁日日躁夜夜躁夜夜| 精品福利永久在线观看| 久久久久久人人人人人| 亚洲综合色网址| 日本精品一区二区三区蜜桃| 国产成+人综合+亚洲专区| 久久久久国内视频| 日韩制服丝袜自拍偷拍| 国产精品 欧美亚洲| 岛国毛片在线播放| 色婷婷久久久亚洲欧美| 青草久久国产| 久久国产乱子伦精品免费另类| 成人影院久久| 成在线人永久免费视频| 在线观看免费视频日本深夜| 很黄的视频免费| 80岁老熟妇乱子伦牲交| 亚洲成人手机| 午夜福利在线观看吧| 欧美+亚洲+日韩+国产| 亚洲国产中文字幕在线视频| 国产三级黄色录像| 精品福利观看| 三级毛片av免费| 久热这里只有精品99| 国产成人免费无遮挡视频| 高清黄色对白视频在线免费看| 国产一区二区三区视频了| 18禁裸乳无遮挡动漫免费视频| 亚洲,欧美精品.| 欧美乱码精品一区二区三区| 欧美亚洲 丝袜 人妻 在线| 乱人伦中国视频| www.熟女人妻精品国产| 国产一区二区三区在线臀色熟女 | 亚洲第一青青草原| 在线国产一区二区在线| 1024视频免费在线观看| 国产精品影院久久| 黑人猛操日本美女一级片| 久久精品人人爽人人爽视色| 一二三四在线观看免费中文在| 热99re8久久精品国产| 两个人免费观看高清视频| 中文字幕高清在线视频| 亚洲一区二区三区欧美精品| 午夜免费成人在线视频| av天堂久久9| 国产av精品麻豆| 日韩中文字幕欧美一区二区| 十分钟在线观看高清视频www| 亚洲第一青青草原| 国产精品久久久久成人av| 欧美日韩视频精品一区| 精品国产一区二区三区四区第35| av国产精品久久久久影院| 最近最新中文字幕大全免费视频| 国产高清视频在线播放一区| 亚洲国产中文字幕在线视频| 成人亚洲精品一区在线观看| 女同久久另类99精品国产91| 最新在线观看一区二区三区| 亚洲久久久国产精品| www日本在线高清视频| 在线观看日韩欧美| 国产av一区二区精品久久| 色综合婷婷激情| 十八禁人妻一区二区| 无遮挡黄片免费观看| 亚洲精品久久午夜乱码| 久久久久久久精品吃奶| 精品电影一区二区在线| 久久久久精品国产欧美久久久| 99国产综合亚洲精品| 亚洲精品在线美女| 婷婷丁香在线五月| 成人免费观看视频高清| 在线观看免费视频日本深夜| 亚洲精品乱久久久久久| 啦啦啦 在线观看视频| 亚洲精品一卡2卡三卡4卡5卡| aaaaa片日本免费| 人人妻人人爽人人添夜夜欢视频| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美乱妇无乱码| 国产又色又爽无遮挡免费看| 热re99久久国产66热| 欧美日韩成人在线一区二区| 丁香欧美五月| 亚洲精品av麻豆狂野| 亚洲av熟女| 亚洲熟女毛片儿| 视频区图区小说| 亚洲国产毛片av蜜桃av| 久久亚洲真实| 日日爽夜夜爽网站| 成在线人永久免费视频| 国产精品乱码一区二三区的特点 | 中文亚洲av片在线观看爽 | 欧美 日韩 精品 国产| 色婷婷av一区二区三区视频| 午夜福利欧美成人| 丰满迷人的少妇在线观看| 久久天堂一区二区三区四区| 亚洲 欧美一区二区三区| 建设人人有责人人尽责人人享有的| а√天堂www在线а√下载 | 捣出白浆h1v1| 国产成+人综合+亚洲专区| 电影成人av| 妹子高潮喷水视频| 热99re8久久精品国产| 亚洲欧美色中文字幕在线| 免费久久久久久久精品成人欧美视频| 亚洲成人国产一区在线观看| 免费一级毛片在线播放高清视频 | 一进一出抽搐gif免费好疼 | 久久ye,这里只有精品| 成年人免费黄色播放视频| 天天躁日日躁夜夜躁夜夜| 亚洲五月天丁香| 999久久久国产精品视频| 大片电影免费在线观看免费| 亚洲av美国av|