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

    三維似French模型高精度彈性波正演數(shù)值模擬及其波場(chǎng)特征

    2010-09-07 02:24:30劉阿男劉思彤陳可洋
    關(guān)鍵詞:快照波場(chǎng)縱波

    劉阿男,劉思彤,陳可洋,楊 微

    (1.東北石油大學(xué)地球科學(xué)學(xué)院,黑龍江大慶 163318; 2.大慶油田有限責(zé)任公司勘探開發(fā)研究院,黑龍江大慶163712; 3.大慶油田有限責(zé)任公司第六采油廠,黑龍江大慶 163712)

    三維似French模型高精度彈性波正演數(shù)值模擬及其波場(chǎng)特征

    劉阿男1,劉思彤1,陳可洋2,楊 微3

    (1.東北石油大學(xué)地球科學(xué)學(xué)院,黑龍江大慶 163318; 2.大慶油田有限責(zé)任公司勘探開發(fā)研究院,黑龍江大慶163712; 3.大慶油田有限責(zé)任公司第六采油廠,黑龍江大慶 163712)

    基于三維一階雙曲型各向異性介質(zhì)速度-應(yīng)力彈性波動(dòng)方程,給出各向異性介質(zhì)向各向同性介質(zhì)轉(zhuǎn)化的條件,采用高階交錯(cuò)網(wǎng)格有限差分法和PML吸收邊界條件,實(shí)現(xiàn)高精度三維彈性波正演數(shù)值模擬.以均勻介質(zhì)為例,研究彈性波在三維空間的波場(chǎng)特征;以三維似French模型為例,合成不同時(shí)刻的三分量波場(chǎng)快照和共炮點(diǎn)記錄,并對(duì)比分析三維各向同性介質(zhì)和各向異性介質(zhì)在該模型中的波場(chǎng)響應(yīng)特征.數(shù)值結(jié)果驗(yàn)證了方法的正確性和有效性,同時(shí)具有精度高和邊界吸收效果明顯的優(yōu)點(diǎn).

    三維彈性介質(zhì);高精度正演數(shù)值模擬;PML吸收邊界條件;似French模型;高階交錯(cuò)網(wǎng)格有限差分法

    0 引言

    彈性波正演數(shù)值模擬在多分量地震資料處理、觀測(cè)系統(tǒng)采集設(shè)計(jì)、處理參數(shù)試驗(yàn)、解釋方案驗(yàn)證等方面發(fā)揮越來(lái)越重要作用,是地震學(xué)和地震勘探學(xué)研究領(lǐng)域中的前沿課題之一,也是難度較大的課題之一[1].近些年,隨計(jì)算機(jī)計(jì)算能力的增強(qiáng)該技術(shù)得到飛躍發(fā)展.彈性波正演數(shù)值模擬技術(shù)大多是基于二維波場(chǎng)觀測(cè)研究地震勘探問(wèn)題,原因在于二維彈性波場(chǎng)的計(jì)算速度快、算法易實(shí)現(xiàn)、剖面易顯示,而實(shí)際的地球介質(zhì)是三維空間,二維波場(chǎng)觀測(cè)忽略三維構(gòu)造體的側(cè)向效應(yīng)等很多重要信息,無(wú)法準(zhǔn)確、真實(shí)地刻畫三維空間的彈性波場(chǎng)響應(yīng)特征.只有通過(guò)三維彈性波正演數(shù)值模擬,才能全面地認(rèn)識(shí)彈性介質(zhì)與彈性波場(chǎng)之間的耦合關(guān)系及彈性波場(chǎng)的空間變化特征,后者對(duì)于研究和認(rèn)識(shí)地球介質(zhì)的各向異性特征尤為重要,然而三維彈性波場(chǎng)正演問(wèn)題研究困難,其原因在于:三維空間的計(jì)算變量較多;空間網(wǎng)格節(jié)點(diǎn)配置關(guān)系復(fù)雜;邊界吸收問(wèn)題處理較難;剖面立體顯示困難;計(jì)算量和存儲(chǔ)量巨大,三維聲波正演數(shù)值模擬的算法要比彈性波正演數(shù)值模擬的算法更為簡(jiǎn)單、易實(shí)現(xiàn).

    目前,對(duì)三維聲波正演問(wèn)題已得到豐富的研究成果[2-7],如頻散問(wèn)題[8-15](精度問(wèn)題)、邊界吸收問(wèn)題[16]、穩(wěn)定性問(wèn)題[17]等,而對(duì)三維彈性波正演問(wèn)題的研究相對(duì)較少,且大多是基于均勻單層介質(zhì)研究彈性波場(chǎng)的傳播特征(波場(chǎng)響應(yīng)),如IgelH等[18]討論三維各向異性介質(zhì)中彈性波數(shù)值模擬的有限差分實(shí)現(xiàn)方法;張文生等[19]對(duì)三維正交各向異性介質(zhì)彈性波場(chǎng)進(jìn)行有限差分模擬;裴正林[20-21]采用高階交錯(cuò)網(wǎng)格有限差分法對(duì)三維各向異性介質(zhì)中彈性波場(chǎng)進(jìn)行高精度數(shù)值模擬;董良國(guó)等[17]提出三維任意各向異性介質(zhì)彈性波動(dòng)方程的穩(wěn)定性條件;夏凡等[22-23]提出三維均勻介質(zhì)彈性波動(dòng)方程的吸收邊界條件及其高精度數(shù)值計(jì)算方法.

    筆者借助Linux集群計(jì)算機(jī)的高速運(yùn)算能力和大容量存儲(chǔ)能力,采用高階交錯(cuò)網(wǎng)格有限差分法求解三維各向異性介質(zhì)彈性波動(dòng)方程,通過(guò)數(shù)值實(shí)例分析彈性波在三維空間經(jīng)典似French速度模型中的三分量波場(chǎng)響應(yīng)特征,給出三維各向異性介質(zhì)向各向同性介質(zhì)轉(zhuǎn)化的條件,同時(shí)使用任意方位波場(chǎng)切片技術(shù)(包括波場(chǎng)快照切片、模擬記錄切片等),實(shí)現(xiàn)三分量彈性波場(chǎng)全方位可視化顯示,研究三分量彈性波場(chǎng)的傳播規(guī)律,以此驗(yàn)證方法的準(zhǔn)確有效性.

    1 基本理論

    通常三維二階各向異性介質(zhì)(VTI介質(zhì))彈性波動(dòng)方程可表示為

    式中:U、V、W分別為x、y、z方向質(zhì)點(diǎn)的振動(dòng)位移;C為三維彈性張量(矩陣組合形式C也可稱為物性矩陣,一般情況下沿主對(duì)角線方向是對(duì)稱正定的);ρ為彈性介質(zhì)密度;t為時(shí)間.參數(shù)均為空間坐標(biāo)的函數(shù).

    當(dāng)各向異性介質(zhì)對(duì)稱軸與觀測(cè)系統(tǒng)不一致時(shí),需要進(jìn)行觀測(cè)坐標(biāo)旋轉(zhuǎn)或者物性矩陣旋轉(zhuǎn),以達(dá)到合理觀測(cè)目的.此時(shí)需要用到Bond矩陣旋轉(zhuǎn)變換,對(duì)物性矩陣旋轉(zhuǎn)運(yùn)算方程為C′=MCMT,其中M為正交正定矩陣,C′為變換后的物性矩陣.

    根據(jù)Thomsen參數(shù)的定義,準(zhǔn)縱波相速度Vp0=C33/ρ;準(zhǔn)橫波相速度Vs0=C44/ρ;度量準(zhǔn)縱波各向異性強(qiáng)度的參數(shù)ε=(C11-C33)/(2C33);度量準(zhǔn)橫波各向異性強(qiáng)度的參數(shù)γ=(C66-C44)/(2C44);影響垂直于VTI介質(zhì)對(duì)稱軸方向附加準(zhǔn)縱波速度的參數(shù)δ=[(C13+C44)2-(C33-C44)2]/[2C33(C33-C44)].

    為實(shí)現(xiàn)三維各向異性介質(zhì)向各向同性介質(zhì)轉(zhuǎn)化,借用各向同性介質(zhì)中的泊松比ν概念,建立三維VTI介質(zhì)中準(zhǔn)縱波和準(zhǔn)橫波相速度之間的關(guān)系式,即

    同時(shí)可推得

    為便于應(yīng)用高階交錯(cuò)網(wǎng)格有限差分法,將三維二階各向異性彈性波動(dòng)方程(式(1))變換為一階雙曲型速度—應(yīng)力彈性波動(dòng)方程的微分形式,即

    采用高階交錯(cuò)網(wǎng)格有限差分法(時(shí)間2階、空間10階差分精度)對(duì)式(2)進(jìn)行高精度差分離散(形成差分行式),對(duì)離散方程中的各變量進(jìn)行交錯(cuò)網(wǎng)格節(jié)點(diǎn)配置,可得到三維高精度彈性波正演遞推算子.在人工截?cái)噙吔缣?采用統(tǒng)一格式的最佳匹配層吸收邊界條件(吸收層厚度取為20個(gè)節(jié)點(diǎn)),消除或削弱截?cái)噙吔缣幍姆瓷洳?要保證邊界計(jì)算過(guò)程穩(wěn)定,還需滿足差分迭代所需的穩(wěn)定性條件[17].

    2 波場(chǎng)特征

    2.1 均勻介質(zhì)三維三分量

    三維均勻各向同性介質(zhì)模型尺寸為500m×500m×500m,3個(gè)方向空間網(wǎng)格長(zhǎng)度為5m,采用純縱波震源(在3個(gè)方向正應(yīng)力場(chǎng)中加入等量外力擾動(dòng))并置震源于模型中央,其最大頻率為40Hz(波形為雷克子波),介質(zhì)的縱波速度為2000m/s;泊松比為0.25;密度為1g/cm3;時(shí)間步長(zhǎng)取0.5ms.滿足穩(wěn)定性條件,波場(chǎng)快照的記錄時(shí)間為0.125s.波場(chǎng)任意方向切片顯示方式的定義:(x=A,yoz)代表沿x方向深度A處且在yoz平面內(nèi)的波場(chǎng)切片(其中y為圖片的橫坐標(biāo);z為圖片的縱坐標(biāo));(t=B,yoz)代表在時(shí)間B節(jié)點(diǎn)處且在yoz平面內(nèi)的波場(chǎng)切片,其他情況依此類推.

    三維任一分量彈性波場(chǎng)在3個(gè)相互垂直方向上的波場(chǎng)快照切片見(jiàn)圖1.由圖1可知:在純縱波震源激勵(lì)下,x方向質(zhì)點(diǎn)速度場(chǎng)中縱波波形與yoz平面對(duì)稱;y方向質(zhì)點(diǎn)速度場(chǎng)中縱波波形與xoz平面對(duì)稱;z方向質(zhì)點(diǎn)速度場(chǎng)中縱波波形與yox平面對(duì)稱,且對(duì)稱面兩側(cè)極性相反.這驗(yàn)證在均勻各向同性介質(zhì)中,三維純縱波波場(chǎng)具有球?qū)ΨQ的特點(diǎn).純橫波三分量彈性波場(chǎng)在任意方向的質(zhì)點(diǎn)振動(dòng)速度和與其相垂直的平面存在對(duì)稱關(guān)系,且極性相一致(見(jiàn)圖2,在x分量質(zhì)點(diǎn)的振動(dòng)速度處加入外力擾動(dòng)).數(shù)值計(jì)算結(jié)果表明計(jì)算精度較高、頻散現(xiàn)象很弱.該三維三分量彈性波場(chǎng)特征分析方法同樣適合于表征三維各向異性介質(zhì)準(zhǔn)縱波和準(zhǔn)橫波波場(chǎng)的特點(diǎn),受各向異性參數(shù)的影響,準(zhǔn)縱波和準(zhǔn)橫波波形不再是球?qū)ΨQ關(guān)系,而是橢圓形甚至是具有不規(guī)則特征的對(duì)稱關(guān)系.

    2.2 三維似French模型

    圖1 0.125ms時(shí)刻縱波震源情況下三維各向同性介質(zhì)彈性波波場(chǎng)快照

    采用似French模型的速度模型,包含1個(gè)水平層界面、1個(gè)傾斜斷層、1個(gè)隆起構(gòu)造、1個(gè)凹陷構(gòu)造(見(jiàn)圖3).模型尺寸為500m×500m×500m,3個(gè)方向空間網(wǎng)格長(zhǎng)度為5m,震源位于模型地表中央位置,其最大頻率為50Hz,每層密度均為1g/cm3,時(shí)間步長(zhǎng)為0.5ms,滿足穩(wěn)定性條件,接收時(shí)間為1s.對(duì)于三維各向同性介質(zhì)模型,上層介質(zhì)縱波速度為2000m/s,下層介質(zhì)縱波速度為2500m/s,泊松比為0.25;對(duì)于三維各向異性介質(zhì)模型,上層介質(zhì)準(zhǔn)縱波速度為2000m/s,下層介質(zhì)準(zhǔn)縱波速度為2500m/s,泊松比為0.25,ε為0.1,δ為0.2,γ為0.

    圖2 0.125ms時(shí)刻混合源情況下三維各向同性介質(zhì)彈性波波場(chǎng)快照(x分量)

    圖3 三維似French模型(單位:m)

    在0.20,0.25s時(shí),三維各向同性介質(zhì)和各向異性介質(zhì)波場(chǎng)快照中的x、y、z方向的中間位置,且沿平行于yoz、xoz、xoy平面的波場(chǎng)快照切片(只給出三維垂直分量的波場(chǎng)快照)見(jiàn)圖4和圖5.考慮三維似French速度模型與三維各向同性介質(zhì)彈性波場(chǎng)快照,經(jīng)過(guò)對(duì)比容易識(shí)別直達(dá)P波、傾斜界面的反射PP波和反射PS波、水平層界面引起的反射PP波和反射PS波、隆起構(gòu)造的反射PP波和反射PS波、凹陷構(gòu)造引起的回轉(zhuǎn)反射PP波和反射PS波;彈性波傳播至下層介質(zhì)時(shí)形成的透射PP波和透射PS波,以及彈性波在隆起構(gòu)造和凹陷構(gòu)造之間相互作用形成的復(fù)雜波(包括轉(zhuǎn)換波能量).與三維各向異性介質(zhì)彈性波場(chǎng)快照對(duì)比,彈性波場(chǎng)更為復(fù)雜,波場(chǎng)快照中包含由各向異性特性引起的直達(dá)S波,模型中構(gòu)造形成的反射SP波、反射SS波、透射SP波、透射SS波,以及這些波在構(gòu)造之間相互作用形成的復(fù)雜波形,直達(dá)S波的波場(chǎng)能量比其他反射波能量要強(qiáng).由圖4可知,在不同的切片中,凹陷構(gòu)造和隆起構(gòu)造的側(cè)向傳播效應(yīng)較為明顯.三維各向同性介質(zhì)和各向異性介質(zhì)數(shù)值模擬記錄的在不同位置處的波場(chǎng)切片(見(jiàn)圖5),兩者最明顯的特征在于各向異性介質(zhì)中存在直達(dá)S波及其轉(zhuǎn)換波,其波場(chǎng)情況明顯比各向同性介質(zhì)更為復(fù)雜.由此可知,三維彈性波在三維空間中進(jìn)行傳播,僅考慮二維空間的彈性波傳播不夠準(zhǔn)確和全面.另外,地震波到達(dá)模型邊界時(shí)無(wú)邊界反射波形成,三維彈性波場(chǎng)的傳播過(guò)程中無(wú)明顯的頻散現(xiàn)象,說(shuō)明數(shù)值計(jì)算方法具有較高精度.

    圖4 三維彈性波場(chǎng)快照切片顯示

    圖5 三維彈性波模擬記錄切片顯示

    3 結(jié)論

    (1)采用高階交錯(cuò)網(wǎng)格有限差分法求解三維一階雙曲型各向異性介質(zhì)速度—應(yīng)力彈性波動(dòng)方程及統(tǒng)一格式的PML吸收邊界條件方程構(gòu)建三維高精度彈性波正演遞推算子,給出三維VTI介質(zhì)與三維各向同性介質(zhì)之間的轉(zhuǎn)換關(guān)系,實(shí)現(xiàn)了三維均勻介質(zhì)和似French模型的彈性波高精度正演數(shù)值模擬.

    (2)使用任意方位彈性波場(chǎng)切片技術(shù)實(shí)現(xiàn)三維三分量彈性波場(chǎng)全方位可視化,研究三分量彈性波場(chǎng)在三維空間中的波場(chǎng)特征及其傳播規(guī)律,縱波或準(zhǔn)縱波具有平面對(duì)稱、極性相反的特點(diǎn);橫波或準(zhǔn)橫波具有平面對(duì)稱、極性相同的特點(diǎn).

    (3)三維彈性波場(chǎng)正演數(shù)值模擬技術(shù)可以更加全面地表征彈性波在全三維空間中的傳播過(guò)程,計(jì)算方法準(zhǔn)確有效.

    [1] 裴正林,牟永光.地震波傳播的數(shù)值模擬[J].地球物理學(xué)進(jìn)展,2004,19(4):933-941.

    [2] 劉金平,張向君,董建平.波動(dòng)方程顯式短算子三維正演模擬[J].吉林大學(xué)學(xué)報(bào):地球科學(xué)版,2004,34(1):138-141.

    [3] 王潤(rùn)秋,張明,費(fèi)建博,等.精細(xì)積分法三維地震波正演模擬[J].勘探地球物理進(jìn)展,2006,29(6):394-397.

    [4] 何兵壽,魏修成,劉洋.三維波動(dòng)方程的數(shù)值頻散關(guān)系及其疊前和疊后數(shù)值模擬[J].石油大學(xué)學(xué)報(bào):自然科學(xué)版,2001,25(1):67-71.

    [5] 熊曉軍,賀振華,黃德濟(jì).三維波動(dòng)方程正演及模型應(yīng)用研究[J].石油物探,2005,44(6):554-556.

    [6] 熊高君,張學(xué)工,賀振華,等.三維定位原理與三維反射波場(chǎng)模擬[J].礦物巖石,2002,22(3):93-97.

    [7] 熊高君,賀振華,張毅祥.三維混合延拓一步法波動(dòng)方程正演模擬[J].成都理工學(xué)院學(xué)報(bào),1999,26(1):48-51.

    [8] Alford R M,Kelly K R,Boore D M.Accuracy of finite-difference modeling of the acoustic wave equation[J].Geophysics,1974,39 (6):834-842.

    [9] Marfurt KJ.Accuracy of finite difference and finite element modeling of the scalar and elastic equations[J].Geophysics,1984,49(5): 533-549.

    [10] Dablain M A.The application of high-differencing to the scalar wave equation[J].Geophysics,1986,51(1):54-66.

    [11] Fei T,Larner K.Elimination of numerical dispersion in finite-difference modeling and migration by flux-corrected transport[J].Geophysics,1995,60(6):1830-1842.

    [12] 董良國(guó),李培明.地震波傳播數(shù)值模擬中的頻散問(wèn)題[J].天然氣工業(yè),2004,24(6):53-56.

    [13] 吳國(guó)忱,王華忠.波場(chǎng)模擬中的數(shù)值頻散分析與校正策略[J].地球物理學(xué)進(jìn)展,2005,20(1):58-65.

    [14] 寧剛,熊章強(qiáng),陳持.波動(dòng)方程有限差分正演模擬誤差來(lái)源分析[J].物探與化探,2008,32(2):203-206.

    [15] 董良國(guó),馬在田,曹景忠,等.一階彈性波方程交錯(cuò)網(wǎng)格高階差分解法[J].地球物理學(xué)報(bào),2000,43(3):411-419.

    [16] 陳可洋,楊微.優(yōu)化的三維地震波旁軸近似吸收邊界條件[J].勘探地球物理進(jìn)展,2009,32(3):179-181,206.

    [17] 董良國(guó),馬在田,曹景忠.一階彈性波方程交錯(cuò)網(wǎng)格高階差分解法穩(wěn)定性研究[J].地球物理學(xué)報(bào),2000,43(6):856-864.

    [18] Igel H,Morap,Riolletb.Anisotropic wave propagation through finite-difference grids[J].Geophysics,1995,60(4):1203-1216.

    [19] 張文生,宋海斌.三維正交各向異性介質(zhì)三分量高精度有限差分正演模擬[J].石油地球物理勘探,2001,36(4):422-432.

    [20] 裴正林.三維各向異性介質(zhì)中彈性波方程交錯(cuò)網(wǎng)格高階有限差分法數(shù)值模擬[J].石油大學(xué)學(xué)報(bào):自然科學(xué)版,2004,28(5):23-29.

    [21] 裴正林.三維各向同性介質(zhì)彈性波方程交錯(cuò)網(wǎng)格高階有限差分法模擬[J].石油物探,2005,44(4):308-315.

    [22] Xia F,Dong L G,Ma Z T.Absorbing boundary conditions in 3D elastic-wave numerical modeling[J].Chinese J.Geophys,2004,47 (1):132-136.

    [23] Xia F,Dong L G,Ma Z T.The numerical modeling of 3-D elastic wave equation using a high-order,staggered-grid,finite difference scheme[J].Applied Geophysics,2004,1(1):38-41.

    High-precision elastic w ave forw ard numerical simulation and its w ave field characteristic analysis on 3D like French model/2010,34(4):45-50

    LIU A-nan1,LIU Si-tong1,CHEN Ke-yang2,YANG Wei3
    (1.Geoscience College,Northeast Petroleum University,Daqing,Heilongjiang163318,China;2. Ex ploration and Development Institute ofDaqing Oilf ield Company L td.,Daqing,Heilongjiang 163712,China;3.Oil Recovery Plant N o.6,Daqing Oilf ield Corp.L td.,Daqing,Heilongjiang 163712,China)

    3D elastic wave forward numerical simulation is the key means for studying seismic wave’s wave field characteristic and its propagating rules in complex 3D spatial elastic wave medium.Based on 3D one-order hyperbolic anisotropic medium velocity-stress elastic wave equation,the paper gives out the transformation condition from anisotropic medium to isotropic medium,and uses high-order staggered-grid finite-difference method and PML absorbing boundary condition to achieve high-precision 3D elastic wave forward numerical simulation.Take an isotropic medium as an example to study the wave field characteristic in 3D space,and take 3D like French model as another model to obtain three-component snapshots with different time and common-shot record,compare and analyze the wave field response in this model of 3D isotropic medium and anisotropic medium.The numerical result verifies the method’s correctness and effectiveness,and shows the advantages of high-precision and obvious perfect boundary absorption effect.

    3D elastic medium ; high-precision forward numerical simulation ; PML absorbing boundarycondition ; French like model ; high-order staggered-grid finite-difference met hod

    book=4,ebook=364

    TE19

    A

    1000-1891(2010)04-0045-06

    2010-06-22;審稿人:王云專;編輯:任志平,張兆虹

    國(guó)家“973”重大基礎(chǔ)研究項(xiàng)目(2001CB209104)

    劉阿男(1973-),女,碩士,工程師,主要從事油氣成藏方面的研究.

    猜你喜歡
    快照波場(chǎng)縱波
    EMC存儲(chǔ)快照功能分析
    天津科技(2022年5期)2022-05-31 02:18:08
    彈性波波場(chǎng)分離方法對(duì)比及其在逆時(shí)偏移成像中的應(yīng)用
    黃257井區(qū)疊前縱波方位各向異性裂縫分布預(yù)測(cè)
    創(chuàng)建磁盤組備份快照
    交錯(cuò)網(wǎng)格與旋轉(zhuǎn)交錯(cuò)網(wǎng)格對(duì)VTI介質(zhì)波場(chǎng)分離的影響分析
    基于Hilbert變換的全波場(chǎng)分離逆時(shí)偏移成像
    數(shù)據(jù)恢復(fù)的快照策略
    變截面階梯桿中的縱波傳播特性實(shí)驗(yàn)
    旋轉(zhuǎn)交錯(cuò)網(wǎng)格VTI介質(zhì)波場(chǎng)模擬與波場(chǎng)分解
    一張“快照”搞定人體安檢
    91成人精品电影| 高清视频免费观看一区二区| 成人国产一区最新在线观看| 欧洲精品卡2卡3卡4卡5卡区| 免费在线观看完整版高清| 伊人久久大香线蕉亚洲五| 久久久久久久国产电影| 久久久久久人人人人人| 欧美日韩瑟瑟在线播放| 国产免费男女视频| 一本综合久久免费| 国产在视频线精品| 国产精品一区二区免费欧美| 一级毛片高清免费大全| 久久久久久久久免费视频了| 久久国产精品人妻蜜桃| 看免费av毛片| 亚洲 欧美一区二区三区| 亚洲一码二码三码区别大吗| 欧美黑人欧美精品刺激| 黄色女人牲交| 国产在线观看jvid| 女人高潮潮喷娇喘18禁视频| 亚洲精品在线观看二区| 久久久久久免费高清国产稀缺| 老司机亚洲免费影院| 一级毛片高清免费大全| 脱女人内裤的视频| 国产在线精品亚洲第一网站| www.熟女人妻精品国产| 一级片'在线观看视频| 色综合欧美亚洲国产小说| 日本一区二区免费在线视频| 亚洲欧美激情综合另类| 婷婷丁香在线五月| 久久久久久亚洲精品国产蜜桃av| 日本wwww免费看| 精品无人区乱码1区二区| 国产精品久久久久久人妻精品电影| 久久人妻福利社区极品人妻图片| 国产真人三级小视频在线观看| 咕卡用的链子| 欧美黄色淫秽网站| 亚洲aⅴ乱码一区二区在线播放 | 国产精品免费视频内射| 超碰成人久久| 亚洲专区字幕在线| 亚洲一区二区三区不卡视频| 亚洲欧美一区二区三区久久| 法律面前人人平等表现在哪些方面| 国产亚洲欧美98| 亚洲成国产人片在线观看| av片东京热男人的天堂| 亚洲七黄色美女视频| 69精品国产乱码久久久| 丝袜人妻中文字幕| 夜夜爽天天搞| 午夜福利免费观看在线| 亚洲av成人av| 在线永久观看黄色视频| 两个人免费观看高清视频| 两个人免费观看高清视频| 国产精品国产av在线观看| 99国产极品粉嫩在线观看| 午夜精品国产一区二区电影| 欧美精品高潮呻吟av久久| 国产欧美日韩一区二区三| 岛国毛片在线播放| 叶爱在线成人免费视频播放| 男女免费视频国产| 日韩大码丰满熟妇| 欧美日韩视频精品一区| 99精品久久久久人妻精品| 国产av又大| 少妇裸体淫交视频免费看高清 | 在线永久观看黄色视频| 人妻 亚洲 视频| 日本精品一区二区三区蜜桃| 欧美在线黄色| 91成人精品电影| 久久久久久久国产电影| 一级a爱视频在线免费观看| 在线观看舔阴道视频| 丰满迷人的少妇在线观看| 中文字幕人妻熟女乱码| 久久久国产精品麻豆| 国产不卡一卡二| 久热爱精品视频在线9| 欧美精品av麻豆av| 一级作爱视频免费观看| 视频区欧美日本亚洲| 国产成人系列免费观看| 一二三四在线观看免费中文在| 国产日韩一区二区三区精品不卡| 男人的好看免费观看在线视频 | 亚洲欧美日韩高清在线视频| 国产av又大| 老司机影院毛片| 国产不卡一卡二| 热re99久久精品国产66热6| 色综合婷婷激情| 久久久久精品国产欧美久久久| 夫妻午夜视频| 国产极品粉嫩免费观看在线| 在线观看免费日韩欧美大片| 交换朋友夫妻互换小说| 色精品久久人妻99蜜桃| 久久午夜亚洲精品久久| 嫁个100分男人电影在线观看| 这个男人来自地球电影免费观看| 欧美精品人与动牲交sv欧美| 99精国产麻豆久久婷婷| 国产欧美亚洲国产| 久久香蕉国产精品| 99精国产麻豆久久婷婷| 黄色视频不卡| 91麻豆精品激情在线观看国产 | 国产激情欧美一区二区| 成人影院久久| 亚洲欧洲精品一区二区精品久久久| 高清在线国产一区| 国产精品一区二区在线不卡| 高潮久久久久久久久久久不卡| 日日摸夜夜添夜夜添小说| 亚洲av成人一区二区三| 女人爽到高潮嗷嗷叫在线视频| 亚洲色图av天堂| 精品一区二区三区视频在线观看免费 | 18禁裸乳无遮挡免费网站照片 | 国产成人精品久久二区二区免费| 久久香蕉国产精品| 国产男女超爽视频在线观看| 99精品久久久久人妻精品| 国产精品一区二区在线不卡| 欧美日韩瑟瑟在线播放| 12—13女人毛片做爰片一| 乱人伦中国视频| 搡老熟女国产l中国老女人| 成人18禁高潮啪啪吃奶动态图| √禁漫天堂资源中文www| 精品午夜福利视频在线观看一区| 男女高潮啪啪啪动态图| 国产1区2区3区精品| 国产高清国产精品国产三级| 黄片播放在线免费| www.熟女人妻精品国产| 国产精品电影一区二区三区 | 黑人巨大精品欧美一区二区mp4| 在线观看舔阴道视频| 一本综合久久免费| 欧美乱色亚洲激情| 免费在线观看视频国产中文字幕亚洲| 99国产精品99久久久久| 一本大道久久a久久精品| 色综合婷婷激情| 老司机亚洲免费影院| av免费在线观看网站| 成人特级黄色片久久久久久久| 欧美+亚洲+日韩+国产| 男女午夜视频在线观看| 一a级毛片在线观看| 亚洲久久久国产精品| 性色av乱码一区二区三区2| 午夜福利在线观看吧| x7x7x7水蜜桃| 精品国内亚洲2022精品成人 | 成人av一区二区三区在线看| 国产真人三级小视频在线观看| 国产精品免费一区二区三区在线 | 制服诱惑二区| 老司机午夜十八禁免费视频| 两人在一起打扑克的视频| 美女扒开内裤让男人捅视频| 我的亚洲天堂| 国产精品乱码一区二三区的特点 | 免费久久久久久久精品成人欧美视频| 91老司机精品| 欧美性长视频在线观看| 国产精品久久久久久精品古装| 欧美黑人精品巨大| 成年人免费黄色播放视频| 超色免费av| 精品欧美一区二区三区在线| 久久久国产一区二区| 99精品在免费线老司机午夜| 免费在线观看视频国产中文字幕亚洲| 中文亚洲av片在线观看爽 | 精品国产一区二区久久| 欧美另类亚洲清纯唯美| 人人澡人人妻人| 激情视频va一区二区三区| 国产精品二区激情视频| 99精品久久久久人妻精品| 亚洲情色 制服丝袜| 激情视频va一区二区三区| 亚洲,欧美精品.| 桃红色精品国产亚洲av| 成人18禁在线播放| 成年版毛片免费区| 久久久精品区二区三区| 欧美日韩亚洲高清精品| 国产在线一区二区三区精| 国产在线观看jvid| 国产真人三级小视频在线观看| 99re6热这里在线精品视频| 久久久久久久午夜电影 | 亚洲精品粉嫩美女一区| 久久精品人人爽人人爽视色| 男女免费视频国产| 亚洲色图 男人天堂 中文字幕| aaaaa片日本免费| 婷婷丁香在线五月| 伦理电影免费视频| 亚洲三区欧美一区| 国产精品久久久人人做人人爽| 免费观看人在逋| 国产精品美女特级片免费视频播放器 | 午夜久久久在线观看| 精品高清国产在线一区| 黄网站色视频无遮挡免费观看| 午夜免费成人在线视频| 久久中文看片网| 最近最新中文字幕大全免费视频| 久99久视频精品免费| 久久精品国产99精品国产亚洲性色 | 99精品久久久久人妻精品| 久久久久国产一级毛片高清牌| 免费在线观看日本一区| 国产精品一区二区精品视频观看| 嫩草影视91久久| 一进一出好大好爽视频| 日韩欧美免费精品| 很黄的视频免费| 日日摸夜夜添夜夜添小说| 国产精品国产高清国产av | 老司机靠b影院| 国产在线观看jvid| 欧美日本中文国产一区发布| 欧美激情 高清一区二区三区| 黑丝袜美女国产一区| 在线观看免费日韩欧美大片| 精品午夜福利视频在线观看一区| 热re99久久国产66热| 热re99久久精品国产66热6| 精品一区二区三卡| 在线天堂中文资源库| 日韩人妻精品一区2区三区| 久久久精品国产亚洲av高清涩受| 国产人伦9x9x在线观看| 757午夜福利合集在线观看| 嫩草影视91久久| 一本一本久久a久久精品综合妖精| 69av精品久久久久久| 亚洲人成伊人成综合网2020| 久久人妻av系列| 又黄又爽又免费观看的视频| 男女午夜视频在线观看| 国产亚洲精品第一综合不卡| 久久精品熟女亚洲av麻豆精品| av在线播放免费不卡| 久久中文字幕人妻熟女| 99精国产麻豆久久婷婷| 后天国语完整版免费观看| 亚洲熟女精品中文字幕| 欧美精品一区二区免费开放| 99精品久久久久人妻精品| 国产熟女午夜一区二区三区| 成人三级做爰电影| 国产精品久久久久成人av| 大型黄色视频在线免费观看| 欧美老熟妇乱子伦牲交| 日本a在线网址| 久久久久久亚洲精品国产蜜桃av| 亚洲 欧美一区二区三区| 黄网站色视频无遮挡免费观看| 日韩欧美一区二区三区在线观看 | 波多野结衣一区麻豆| 欧美中文综合在线视频| 两个人看的免费小视频| 麻豆成人av在线观看| 精品国产超薄肉色丝袜足j| 午夜影院日韩av| 老汉色av国产亚洲站长工具| 久久午夜亚洲精品久久| 黄色怎么调成土黄色| 不卡一级毛片| 80岁老熟妇乱子伦牲交| 日韩视频一区二区在线观看| 香蕉国产在线看| 国产成人欧美在线观看 | 精品一品国产午夜福利视频| 国产一区二区激情短视频| 亚洲一码二码三码区别大吗| 国产欧美日韩综合在线一区二区| 免费av中文字幕在线| 国产在视频线精品| 国产精品综合久久久久久久免费 | 亚洲黑人精品在线| 久久久久视频综合| 国产成人精品久久二区二区免费| 妹子高潮喷水视频| 亚洲片人在线观看| 久久天躁狠狠躁夜夜2o2o| 国产aⅴ精品一区二区三区波| 曰老女人黄片| 欧美在线黄色| 一级片免费观看大全| 欧美一级毛片孕妇| 欧美日韩瑟瑟在线播放| 精品亚洲成国产av| 中文字幕人妻丝袜制服| 国产麻豆69| netflix在线观看网站| 国产精品影院久久| 黄色视频,在线免费观看| 女警被强在线播放| 国产精品欧美亚洲77777| 精品久久久精品久久久| 老司机靠b影院| 丰满的人妻完整版| 久久中文字幕一级| 精品熟女少妇八av免费久了| 成人18禁在线播放| 757午夜福利合集在线观看| av线在线观看网站| 91精品三级在线观看| 丝袜人妻中文字幕| 欧美日韩亚洲国产一区二区在线观看 | 极品人妻少妇av视频| xxxhd国产人妻xxx| 麻豆成人av在线观看| 精品国产一区二区三区四区第35| 午夜福利免费观看在线| 欧美色视频一区免费| 黄色丝袜av网址大全| 又黄又爽又免费观看的视频| 免费看a级黄色片| 岛国毛片在线播放| 国产精品成人在线| av一本久久久久| 欧美不卡视频在线免费观看 | 女人精品久久久久毛片| 亚洲精品av麻豆狂野| 亚洲精品中文字幕在线视频| 国产欧美日韩一区二区三| 午夜亚洲福利在线播放| 人人妻人人澡人人看| 黄片播放在线免费| 男人舔女人的私密视频| 人人妻人人澡人人爽人人夜夜| 欧美丝袜亚洲另类 | 国产深夜福利视频在线观看| 十八禁网站免费在线| 亚洲熟女精品中文字幕| 久久久久久久久免费视频了| 建设人人有责人人尽责人人享有的| 国产一卡二卡三卡精品| 超色免费av| 高潮久久久久久久久久久不卡| 老熟妇仑乱视频hdxx| 黄频高清免费视频| 黄频高清免费视频| 国产极品粉嫩免费观看在线| 真人做人爱边吃奶动态| 在线观看www视频免费| 一本一本久久a久久精品综合妖精| 欧美日韩视频精品一区| 亚洲avbb在线观看| 女性被躁到高潮视频| 国产一卡二卡三卡精品| 欧美乱妇无乱码| xxx96com| 免费在线观看影片大全网站| 在线天堂中文资源库| 国产欧美日韩一区二区精品| 日日爽夜夜爽网站| 丝瓜视频免费看黄片| 男女下面插进去视频免费观看| 国产精品电影一区二区三区 | 精品无人区乱码1区二区| 一边摸一边抽搐一进一出视频| a级片在线免费高清观看视频| 精品一区二区三区四区五区乱码| 热re99久久国产66热| 免费看a级黄色片| cao死你这个sao货| 久久ye,这里只有精品| 悠悠久久av| 欧美日韩中文字幕国产精品一区二区三区 | 黄色女人牲交| 大香蕉久久成人网| 亚洲在线自拍视频| 亚洲精品中文字幕在线视频| 99re6热这里在线精品视频| 午夜精品在线福利| 欧美最黄视频在线播放免费 | 国产97色在线日韩免费| 女人爽到高潮嗷嗷叫在线视频| 国产亚洲精品第一综合不卡| 丝袜美腿诱惑在线| 国产男女内射视频| 热re99久久精品国产66热6| 国产精品二区激情视频| 久久久精品国产亚洲av高清涩受| 欧美国产精品va在线观看不卡| 80岁老熟妇乱子伦牲交| 午夜老司机福利片| 一级黄色大片毛片| 国产av精品麻豆| 精品午夜福利视频在线观看一区| 亚洲专区国产一区二区| 香蕉丝袜av| 精品国内亚洲2022精品成人 | 在线观看免费视频网站a站| 在线观看午夜福利视频| 国产成人av激情在线播放| 999久久久国产精品视频| 极品教师在线免费播放| 成人三级做爰电影| 亚洲熟女毛片儿| 亚洲精品国产精品久久久不卡| 免费不卡黄色视频| 国产精品免费视频内射| 91精品三级在线观看| 精品无人区乱码1区二区| 欧美性长视频在线观看| avwww免费| 在线观看免费午夜福利视频| 午夜91福利影院| avwww免费| 夜夜夜夜夜久久久久| 国产高清视频在线播放一区| 午夜福利在线免费观看网站| 久久精品亚洲精品国产色婷小说| 国产精品综合久久久久久久免费 | 女警被强在线播放| 91九色精品人成在线观看| 久久狼人影院| 色精品久久人妻99蜜桃| 国产精品永久免费网站| 女人被躁到高潮嗷嗷叫费观| 国内毛片毛片毛片毛片毛片| 欧美精品av麻豆av| av线在线观看网站| 欧美在线一区亚洲| 亚洲在线自拍视频| 精品久久久久久,| 成人手机av| 中文字幕精品免费在线观看视频| 18禁国产床啪视频网站| 亚洲欧美激情综合另类| 侵犯人妻中文字幕一二三四区| 国产免费现黄频在线看| 国产亚洲欧美精品永久| 少妇粗大呻吟视频| 天堂中文最新版在线下载| 亚洲在线自拍视频| 免费少妇av软件| 精品无人区乱码1区二区| 超碰成人久久| 国产国语露脸激情在线看| 欧美精品亚洲一区二区| 国产一区二区三区在线臀色熟女 | 一区二区三区精品91| 亚洲午夜理论影院| tube8黄色片| 母亲3免费完整高清在线观看| 欧美日韩国产mv在线观看视频| 国产精品久久久av美女十八| 亚洲黑人精品在线| 黄色怎么调成土黄色| 欧美+亚洲+日韩+国产| 男人的好看免费观看在线视频 | 久久精品人人爽人人爽视色| 在线播放国产精品三级| 久久精品国产a三级三级三级| 十八禁人妻一区二区| 欧美日韩国产mv在线观看视频| 亚洲熟妇中文字幕五十中出 | 两人在一起打扑克的视频| 一级黄色大片毛片| 国产精品av久久久久免费| 人人澡人人妻人| 99re在线观看精品视频| 国产麻豆69| 欧美日韩黄片免| 在线看a的网站| 免费一级毛片在线播放高清视频 | 最新美女视频免费是黄的| 亚洲国产毛片av蜜桃av| 日韩有码中文字幕| 国产亚洲精品久久久久久毛片 | 国产成人免费无遮挡视频| 国产免费av片在线观看野外av| 视频区图区小说| 一级黄色大片毛片| 亚洲九九香蕉| 91老司机精品| 免费少妇av软件| 国产精品美女特级片免费视频播放器 | 99热网站在线观看| 亚洲精品在线美女| 久久99一区二区三区| 亚洲色图综合在线观看| 国产欧美日韩一区二区三| 久久这里只有精品19| 欧美成人午夜精品| 18禁美女被吸乳视频| xxxhd国产人妻xxx| 在线观看免费视频网站a站| 首页视频小说图片口味搜索| 欧美亚洲日本最大视频资源| cao死你这个sao货| 亚洲中文日韩欧美视频| 午夜福利,免费看| 亚洲第一欧美日韩一区二区三区| 国产在线一区二区三区精| 嫁个100分男人电影在线观看| 下体分泌物呈黄色| 成人免费观看视频高清| 精品一区二区三区av网在线观看| 波多野结衣av一区二区av| 别揉我奶头~嗯~啊~动态视频| 欧美精品亚洲一区二区| 香蕉丝袜av| 午夜免费鲁丝| 午夜福利欧美成人| 国产成人av教育| 动漫黄色视频在线观看| 国产成人精品久久二区二区91| 免费一级毛片在线播放高清视频 | 国产麻豆69| 天堂俺去俺来也www色官网| 电影成人av| 欧美 亚洲 国产 日韩一| 国产三级黄色录像| 午夜久久久在线观看| 国产xxxxx性猛交| 亚洲人成伊人成综合网2020| 在线观看www视频免费| 高清黄色对白视频在线免费看| 99国产精品99久久久久| 99国产精品免费福利视频| 国产成人欧美在线观看 | 一区二区三区精品91| 亚洲色图av天堂| 久久中文字幕人妻熟女| 久久国产亚洲av麻豆专区| 91精品三级在线观看| 国产日韩一区二区三区精品不卡| 色在线成人网| 免费在线观看亚洲国产| 欧美亚洲日本最大视频资源| 色精品久久人妻99蜜桃| 精品福利观看| 男女下面插进去视频免费观看| 韩国av一区二区三区四区| 欧美日韩中文字幕国产精品一区二区三区 | 国产极品粉嫩免费观看在线| 一区二区三区激情视频| 视频在线观看一区二区三区| 国产在视频线精品| 亚洲中文av在线| 精品一区二区三区av网在线观看| 一进一出好大好爽视频| 夫妻午夜视频| 啦啦啦免费观看视频1| 亚洲免费av在线视频| 亚洲伊人色综图| 亚洲色图 男人天堂 中文字幕| 美女国产高潮福利片在线看| а√天堂www在线а√下载 | 99在线人妻在线中文字幕 | tube8黄色片| 91麻豆av在线| 天天躁狠狠躁夜夜躁狠狠躁| 中文字幕人妻熟女乱码| 欧美日韩亚洲高清精品| 久久久国产成人免费| 国产高清视频在线播放一区| 嫁个100分男人电影在线观看| 国产精品1区2区在线观看. | 99久久精品国产亚洲精品| www.自偷自拍.com| 69精品国产乱码久久久| 悠悠久久av| 国产欧美日韩综合在线一区二区| 三上悠亚av全集在线观看| 国产极品粉嫩免费观看在线| 欧美日韩亚洲高清精品| 在线观看日韩欧美| 三上悠亚av全集在线观看| 一进一出抽搐动态| 久久人人97超碰香蕉20202| 美女 人体艺术 gogo| 嫩草影视91久久| 日本vs欧美在线观看视频| av天堂在线播放| 女人精品久久久久毛片| 久久国产乱子伦精品免费另类| av片东京热男人的天堂| 亚洲全国av大片| 精品国产一区二区久久| 极品少妇高潮喷水抽搐| 在线av久久热| 国产又爽黄色视频| 欧美中文综合在线视频| 深夜精品福利| 亚洲午夜理论影院| 国产精品av久久久久免费| 97人妻天天添夜夜摸| 两性夫妻黄色片| 亚洲av成人不卡在线观看播放网| 操美女的视频在线观看| 久久国产乱子伦精品免费另类|