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

    二維層狀VTI介質(zhì)試射射線追蹤方法

    2016-08-05 09:00:09李永博王大勇李慶春李建國
    物探化探計(jì)算技術(shù) 2016年3期

    李永博, 吳 瓊, 王大勇, 李慶春, 李建國

    (1.中國地質(zhì)科學(xué)院 地球物理地球化學(xué)勘查研究所,廊坊 065000;2.長安大學(xué) 地質(zhì)工程與測繪學(xué)院,西安 710054;3.東方地球物理公司 新興物探開發(fā)處,涿州 072751)

    ?

    二維層狀VTI介質(zhì)試射射線追蹤方法

    李永博1, 吳瓊1, 王大勇1, 李慶春2, 李建國3

    (1.中國地質(zhì)科學(xué)院地球物理地球化學(xué)勘查研究所,廊坊065000;2.長安大學(xué)地質(zhì)工程與測繪學(xué)院,西安710054;3.東方地球物理公司新興物探開發(fā)處,涿州072751)

    摘要:地震波場正演模擬是研究地震波傳播規(guī)律的有效手段,為地震勘探提供基本的理論依據(jù)。為提高地震勘探的精度及勘探能力,研究復(fù)雜介質(zhì)中地震波傳播問題顯得尤為重要。針對各向異性介質(zhì)中的地震波場正演問題,基于各向異性介質(zhì)和試射射線追蹤方法的基本理論,分析了橫向各向同性(VTI)介質(zhì)中存在的群相關(guān)系,研究了VTI介質(zhì)中的Snell定律及角度轉(zhuǎn)換關(guān)系,并引入了透射系數(shù)和反射系數(shù),具體討論了試射射線追蹤算法,通過數(shù)值模擬實(shí)現(xiàn)了地震波場正演,客觀合理地模擬了地震波在VTI介質(zhì)中的傳播。

    關(guān)鍵詞:波場正演; 橫向各向同性介質(zhì); 射線追蹤

    0引言

    地震波場正演模擬,為地震勘探提供了重要的理論基礎(chǔ)。常規(guī)的地震勘探方法幾乎都是以均勻介質(zhì)模型為假設(shè)前提,然而地下介質(zhì)基本上都是各向異性的,這使地震波的傳播更加復(fù)雜多變[1-2]。因此,研究地震波在各向異性介質(zhì)中的傳播規(guī)律就顯得尤為重要[4]。

    常見的地震波場正演方法包含射線追蹤法和波動方程法,其中射線追蹤方法是一種快速有效的地震波場正演模擬方法。20世紀(jì)80年代末以來,隨著Kirchhoff積分疊前深度偏移方法的成功運(yùn)用,作為其基礎(chǔ)算法之一的射線追蹤方法也得到了很大的推進(jìn)和發(fā)展[3]。隨后,國內(nèi)、外學(xué)者不斷將各向同性介質(zhì)中的射線追蹤方法引入到各向異性介質(zhì)中,并加以研究和發(fā)展,使各向異性介質(zhì)中的射線追蹤方法日臻完善[4-8,11-16]。

    VTI(Vertical Transversely Isotropy)介質(zhì)是一種常見的各向異性介質(zhì),是典型的具有垂向?qū)ΨQ軸的橫向各向同性介質(zhì)[8]。在VTI介質(zhì)中,地震波的傳播速度及透反射的能量分配受到介質(zhì)各向異性的影響發(fā)生變化。因此,研究VTI介質(zhì)射線追蹤方法對解決介質(zhì)的各向異性問題具有重要意義。這里采用經(jīng)典的試射射線追蹤方法對VTI介質(zhì)進(jìn)行正演試算,并引入透反射系數(shù)計(jì)算地震波的能量分配。

    1方法原理

    研究各向異性介質(zhì),首先要確定介質(zhì)的各向異性參數(shù)和彈性參數(shù)的對應(yīng)關(guān)系,這樣才有利于正確分析各向異性介質(zhì)中影響地震波傳播的因素,進(jìn)而進(jìn)行正演模擬及數(shù)據(jù)處理等研究。VTI介質(zhì)的各向異性性質(zhì)可以用Thomsen參數(shù)[11]或彈性系數(shù)來表示,但后者的物理意義相對不明確。因此,這里采用Thomsen參數(shù)來描述VTI介質(zhì)的各向異性性質(zhì)。

    1.1群相關(guān)系

    圖1 各向異性介質(zhì)中群相關(guān)系示意圖Fig.1 The group-phase relations in        VTI media

    根據(jù)介質(zhì)弱各向異性的假設(shè),給出弱各向異性介質(zhì)中P波和SV波相速度的表達(dá)式[11]:

    (1)

    式中:vp0為P波的垂直速度;vS0為S波的垂直速度;ε表示P波各向異性程度;δ表示P波和S波的變異系數(shù),ε、δ為Thomsen參數(shù)。

    Berryman[17]給出了由相速度轉(zhuǎn)化為群速度的表達(dá)式,如式(2)所示。

    (2)

    Byun[18]給出了一種計(jì)算相速度的表達(dá)式,見式(3)。

    (3)

    由式(3)推導(dǎo)出相角和群角存在如下關(guān)系:

    (4)

    將式(2)代入式(4)可得到關(guān)系式,如式(5)所示。

    (5)

    趙愛華等[9]推導(dǎo)并驗(yàn)證了弱各向異性介質(zhì)中由群角轉(zhuǎn)換為相角的表達(dá)式:

    (6)

    1.2透射與反射

    由于VTI介質(zhì)的對稱軸是垂向的,故Snell定律可表達(dá)為[6]:

    (7)

    式中:P為射線參數(shù);α為入射角。

    VTI介質(zhì)中透反射的角度關(guān)系如圖2所示[8],圖2中αinc、αtra、αref分別表示波矢量入射角、透射角、反射角,φinc、φtra、φref分別表示射線入射角、透射角、反射角,θinc、θtra、θref分別表示波矢量入射相角、透射相角、反射相角,αdip為界面傾角。

    李建國等[8]歸納推導(dǎo)并歸納了VTI介質(zhì)中透反射的角度關(guān)系,即:

    (8)

    Ruger[19]給出了VTI介質(zhì)透反射系數(shù)的近似公式,各向異性介質(zhì)中的透反射系數(shù)包括各向同性項(xiàng)和各向異性項(xiàng)兩部分。以P波為例,VTI介質(zhì)中PP波反射系數(shù)近似解表示為RPP,其各向同性項(xiàng)和各向異性項(xiàng)分別表示為RiPP、RaPP,則:

    RPP(θ)=RiPP(θ)+RaPP(θ)

    (9)

    式中,各向同性項(xiàng)為:

    各向異性項(xiàng)為:

    2試射射線追蹤關(guān)鍵技術(shù)

    試射法是一種經(jīng)典的傳統(tǒng)的兩點(diǎn)間射線追蹤方法,即給定震源和接收點(diǎn)位置,根據(jù)Snell定律調(diào)整從震源出發(fā)的射線直至收斂到接收點(diǎn)。

    VTI介質(zhì)試射法射線追蹤的算法流程示意圖如圖3所示。其中,VTI介質(zhì)射線追蹤實(shí)現(xiàn)過程中,由入射射線計(jì)算透反射射線是整個算法的核心。

    圖2 VTI介質(zhì)中透反射示意圖Fig.1 The transmission and reflection in VTI media

    以P波入射為例,射線及角度關(guān)系如圖2所示,詳細(xì)討論VTI介質(zhì)中由入射射線計(jì)算反射射線的實(shí)現(xiàn)過程。

    圖3 VTI介質(zhì)試射射線追蹤算法流程示意圖Fig.3 Flow chart of shooting ray-tracing      algorithm in VTI media

    1) 給定P波初始射線入射角φinc,根據(jù)式(6)將φinc轉(zhuǎn)化為波矢量入射角θinc。

    2) 將θinc分別代入式(1)和式(8)計(jì)算P波相速度vP1及入射角αinc,則由式(7)可得反射點(diǎn)的射線參數(shù)p為式(10)。

    (10)

    3) 當(dāng)P波入射發(fā)生反射時,滿足Snell定律為:

    (11)

    對于PP反射波, 反射波矢量的相角θref滿足方程:

    (12)

    求解方程(12),得到θref。

    (13)

    VTI介質(zhì)中,計(jì)算PP透射波、PSV反射波及PSV透射波的思路與計(jì)算PP反射波的思路一致。

    在建模方面,針對VTI介質(zhì)結(jié)構(gòu)特點(diǎn),作者采用層狀建模的方法,并對連續(xù)性不同的界面采用相應(yīng)的擬合方法,對光滑連續(xù)的界面采用多項(xiàng)式擬合,對含間斷點(diǎn)的線性界面(如斷層)采用分段線性擬合。

    3模型試算

    3.1水平層狀模型

    水平層狀VTI介質(zhì)的模型示意圖及參數(shù)如圖4(a)所示。

    由圖4(b)可以看出,VTI介質(zhì)試射射線追蹤法與波動方程有限差分法計(jì)算的PP波旅行時相對誤差不超過1.2 %,兩者的相對誤差隨著炮檢距的增大而減小,這是由于波動方程有限差分法的誤差主要分布在炮點(diǎn)附近且計(jì)算精度受網(wǎng)格大小的影響;由圖4(c)、圖4(d)可知,兩種方法得到的單炮記錄中PP反射波基本一致,綜上可知,VTI介質(zhì)試射射線追蹤是正確、可行、有效的。

    3.2斷層模型

    VTI介質(zhì)斷層模型示意圖及參數(shù)如圖5所示,斷層模型的橫縱向斷距均為100 m。

    圖4 VTI介質(zhì)射線追蹤法與波動方程法對比Fig.4 Comparison between and shooting ray-tracing method and wave equation method in VTI media(a) 模型示意圖及參數(shù)(參數(shù)來自文獻(xiàn)[4]);(b) 試射法與波動方程法PP波旅行時相對誤差;(c) 試射射線追蹤PP波單炮記錄;(d) 波動方程法單炮記錄(數(shù)值模擬來自文獻(xiàn)[10]的程序)

    由圖6可知,斷層的PSV反射波的射線路徑及合成記錄上出現(xiàn)盲區(qū),斷層上下盤的反射清晰,且炮集記錄中出現(xiàn)斷面反射,說明試射法數(shù)值模擬結(jié)果能夠正確反映斷層的形態(tài),由此可見,VTI介質(zhì)試射射線追蹤對二維間斷界面模型是有效的。

    3.3透鏡體模型

    VTI介質(zhì)透鏡體模型示意圖及參數(shù)如圖7所示,透鏡體深度范圍為350 m~450 m。

    由圖8可以看出,透鏡體頂面和底面的PP反射波射線路徑及合成記錄的清晰可見,底面PP反射波的極性受到各向異性的影響發(fā)生變化,且在合成記錄上能夠較容易分辨出透鏡體的形態(tài)。結(jié)果表明,VTI介質(zhì)試射射線追蹤對較為復(fù)雜的起伏模型是可行的。

    圖5 模型示意圖及參數(shù)Fig.5 Model diagram and parameters

    4結(jié)論

    試射法是一種經(jīng)典的射線追蹤方法,通過對VTI介質(zhì)的正演模擬,得到以下結(jié)論:

    1)歸納總結(jié)了VTI介質(zhì)速度及角度的表達(dá)式及群相轉(zhuǎn)化公式,并給出了VTI介質(zhì)試射射線追蹤算法中由入射射線計(jì)算透射和反射射線的具體實(shí)現(xiàn)方法。

    圖6 斷層模型PSV波射線追蹤Fig.6 Result of fault model PSV wave ray tracing(a) PSV波單炮射線路徑;(b) PSV波單炮記錄;(c) PSV波炮集記錄

    圖7 模型示意圖及參數(shù)Fig.7 Model diagram and parameters

    2)在VTI介質(zhì)射線追蹤過程中加入透反射系數(shù)的計(jì)算,改善了大多數(shù)各向異性介質(zhì)射線追蹤方法不考慮能量分配的問題。

    3)通過模型試算可知,VTI介質(zhì)試射射線追蹤方法是可行有效的。對起伏界面和間斷界面模型,VTI射線追蹤都能夠較準(zhǔn)確反映界面的形態(tài)和介質(zhì)結(jié)構(gòu)。同時,程序能夠?qū)崿F(xiàn)初至波、PP波及PSV波的單一波或多波組合的正演模擬。由于建模方法,作者提出的正演模擬對復(fù)雜模型還有一定的局限性,有待進(jìn)一步改進(jìn)。此外,該射線追蹤方法可推廣至其他各項(xiàng)異性介質(zhì)(如HTI介質(zhì))及三維正演模擬。

    圖8 透鏡體模型PP波射線追蹤Fig.8 Result of lens model PP wave ray tracing (a) PP波單炮射線路徑;(b) PP波單炮記錄;(c) PP波炮集記錄

    參考文獻(xiàn):

    [1]滕吉文,張中杰,王愛武,等.彈性介質(zhì)各向異性研究沿革現(xiàn)狀與問題[J].地球物理學(xué)進(jìn)展,1992,7(4):14-28.

    TENG J W,ZHANG Z J,WANG A W,et al.The study of anisotropy in elastic medium:evolution,present situation and questions[J].Progress in Geophysics,1992,7(4):14-28.(In Chinese)

    [2]張中杰.地震各向異性研究進(jìn)展[J].地球物理學(xué)進(jìn)展,2002,17(2):281-293.

    ZHANG Z J.A review of the seismic anisotropy and its applications[J].Progress in Geophysics, 2002,17(2):281-293.(In Chinese)

    [3]高爾根,徐果明,趙燚.一種任意界面的逐段迭代射線追蹤方法[J].石油地球物理勘探,1998,33(1):54-60.

    GAO E G,XU G M,ZHAO Y.Segmentally-iterative ray tracing method for any interface[J].OGP,1998,33(1):54-60.(In Chinese)

    [4]孔選林,李錄明,羅省賢,等.各向異性介質(zhì)中地震波射線正演[J].物探化探計(jì)算技術(shù),2008,30(3):178-184.

    KONG X L,LI L M,LUO S X,et al.Seismic wave ray forward in anisotropic medium[J].Computing Techniques for Geophysical and Geochemical Exploration,2008,30(3):178-184.(In Chinese)

    [5]張文生,何樵登,朱建偉,等.橫向各向同性介質(zhì)中群速度的計(jì)算[J].物探化探計(jì)算技術(shù),1997,19(2):97-102.

    ZHANG W S,HE J D,ZHU J W,et al.Calculation of group velocity in transversely isotropic media[J].Computing Techniques for Geophysical and Geochemical Exploration,1997,19(2):97-102. (In Chinese)

    [6]鄧懷群,劉雯林,趙正茂.橫向各向同性介質(zhì)中縱波和轉(zhuǎn)換橫波的快速射線追蹤方法[J].石油物探,2000,39(4):1-11.

    DENG H Q,LIU W L,ZHAO Z M.Fast ray-tracing method for compressional and converted waves in transversely isotropic media[J].GPP,2000,39(4):1-11.(In Chinese)

    [7]趙愛華,張美根,丁志峰.橫向各向同性介質(zhì)中地震波走時模擬[J].地球物理學(xué)報,2006,49(6):1762-1769.

    ZHAO A H,ZHANG M G,DING Z F.Seismic traveltime computation for transversely isotropic media[J].Chinese J.Geophys,2006,49(6):1762-1769.(In Chinese)

    [8]李建國,李彥鵬,郭曉玲.VTI介質(zhì)試射射線追蹤[J].石油地球物理勘探,2010,45(4):491-496.

    LI J G,LI Y P,GUO X L.VTI medium test-firing ray tracing[J].OGP,2010,45(4):491-496.(In Chinese)

    [9]趙愛華,丁志峰.一種弱各向異性介質(zhì)地震波群速度的近似表示新方法[J].地球物理學(xué)進(jìn)展,2005,20(4):916-919.

    ZHAO A H,DING Z F.New approximate expressions of Seismic group velocities for weakly anisotropic media[J].Progress in Geophysics,2005,20(4):916-919.(In Chinese)

    [10]周學(xué)明.交錯網(wǎng)格高階差分?jǐn)?shù)值模擬及疊前逆時偏移[D].西安:長安大學(xué),2010.

    ZHOU X M.Staggered-grid High-order Finite-difference Numerical Simulation and Prestack Reverse-time Migration[D].Xi’an:Chang’an University,2010.(In Chinese)

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

    [12]GAJEWSKI D.PSENCIK I.Computation of high-frequency seismic wavefields in 3-D laterrally inhomogeneous anisotropic media[J].Geophys.J.R.astro.,1987,91:383-411.

    [13]GRECHKA V.Y.MCMECHAN G.A.3-D two-point ray tracing for heterogeneous,weakly transversely isotropic media[J].Geophysics,1996,61(6):1883-1894.

    [15]ROGISTER Y.SLAWINSKI M.A.Analytic solution of raytracing equations for a linearly inhomogeneous and elliptically anisotropic velocity model[J].Geophysics,2005,70(5):37-41.

    [16]DEHGHAN K.FARRA V.NICOLETIS L.Approximate ray tracing for qP-waves in inhomogeneous layered media with weak structural[J].Geophysics,2007,72(5):35-46.

    [17]BERRYMAN J.G.Long-wave elastic anisotropy in transversely isotropic media[J].Geophysics,l979,44:896-917.

    [18]BYUN B.S..Seismic parameters for transversely isotropic media[J].Geophysics,1984,49(11):1908-1914.

    [19]RUGER A.P-wave reflection coefficients for transversely isotropic models with vertical and horizontal axis of symmetry[J].Geophysics,1997,62:713-722.

    收稿日期:2015-03-13改回日期:2015-06-23

    基金項(xiàng)目:國家自然科學(xué)基金(41374145,41304105)

    作者簡介:李永博(1986-),男,碩士,主要研究方向?yàn)榈卣鹂碧胶碗姶欧碧?,E-mail:liyongbo@igge.cn。

    文章編號:1001-1749(2016)03-0396-07

    中圖分類號:P 631.4

    文獻(xiàn)標(biāo)志碼:A

    DOI:10.3969/j.issn.1001-1749.2016.03.16

    Shooting ray tracing in 2-D VTI media

    LI Yong-bo1, WU Qiong1, WANG Da-yong1, LI Qing-chun2, LI Jian-guo3

    (1. Institute of Geophysical and Geochemical Exploration, CAGS, Langfang065000, China;2. College of Geology Engineering and Geomatics, Chang'an University, Xi’an710054, China;3.Bureau of Geophysical prospecting,Zhuozhou072751,China)

    Abstract:Seismic wave numerical simulation is an effective measure to study seismic wave propagation, and it provides theoretical basis for seismic exploration. In order to improve the precision and ability of seismic exploration, it is particularly important to research seismic wave propagation in complex media. According to the seismic forward modeling in anisotropic media, from the theories of anisotropic media and shooting ray-tracing method, and on the basis of existing studies, the group-phase relations of velocity and angle are analyzed respectively in vertical transversely isotropic (VTI) media. Then the Snell theorem and the angular relation are studied in VTI media, and both the reflection coefficient and the transmission coefficient are brought in. The shooting ray-tracing algorithm is discussed specifically, and the seismic forward modeling is achieved by numerical simulation, which can simulate the seismic wave propagation more objectively and reasonably in anisotropic media.

    Key words:numerical simulation; VTI; ray tracing

    欧美av亚洲av综合av国产av | 在线观看三级黄色| 精品一区二区免费观看| 天美传媒精品一区二区| 在线观看国产h片| 精品久久蜜臀av无| 最近最新中文字幕免费大全7| 亚洲人成电影观看| 亚洲一卡2卡3卡4卡5卡精品中文| 成人影院久久| 少妇 在线观看| 在线观看人妻少妇| 伦理电影大哥的女人| av网站免费在线观看视频| 电影成人av| 亚洲av男天堂| 国产亚洲午夜精品一区二区久久| 国产精品久久久久久精品电影小说| 亚洲伊人色综图| 搡老岳熟女国产| 五月开心婷婷网| 免费人妻精品一区二区三区视频| 国产精品一区二区在线不卡| 秋霞在线观看毛片| 少妇被粗大的猛进出69影院| 亚洲国产欧美一区二区综合| 日韩av免费高清视频| 国产免费一区二区三区四区乱码| 自线自在国产av| 人妻一区二区av| 男女免费视频国产| 在线观看三级黄色| av女优亚洲男人天堂| 91成人精品电影| av有码第一页| 69精品国产乱码久久久| 欧美日韩视频高清一区二区三区二| av女优亚洲男人天堂| 黑人巨大精品欧美一区二区蜜桃| 90打野战视频偷拍视频| 一区二区av电影网| 午夜日韩欧美国产| 最近2019中文字幕mv第一页| 欧美av亚洲av综合av国产av | 免费高清在线观看视频在线观看| 欧美少妇被猛烈插入视频| 午夜福利视频精品| 欧美人与性动交α欧美软件| 韩国av在线不卡| 亚洲av在线观看美女高潮| 国产午夜精品一二区理论片| 久久久久久久久久久免费av| 亚洲精品久久久久久婷婷小说| 久久这里只有精品19| 国产片内射在线| 看免费av毛片| 亚洲精品国产av成人精品| 青春草国产在线视频| 九九爱精品视频在线观看| 婷婷成人精品国产| 亚洲第一av免费看| 高清视频免费观看一区二区| 青草久久国产| 欧美少妇被猛烈插入视频| av国产久精品久网站免费入址| 国产成人免费无遮挡视频| a 毛片基地| 伊人久久国产一区二区| 九草在线视频观看| 一级片免费观看大全| 久久久精品国产亚洲av高清涩受| 综合色丁香网| 在线免费观看不下载黄p国产| 人人妻,人人澡人人爽秒播 | 国产又色又爽无遮挡免| 久久久久久久久久久免费av| 亚洲天堂av无毛| av又黄又爽大尺度在线免费看| 99九九在线精品视频| 秋霞在线观看毛片| 国产午夜精品一二区理论片| 久久精品亚洲av国产电影网| 亚洲精品久久午夜乱码| 亚洲av综合色区一区| 中文字幕亚洲精品专区| 国产毛片在线视频| 国产又色又爽无遮挡免| 亚洲av国产av综合av卡| 国产一区二区 视频在线| 成年人午夜在线观看视频| 国产精品久久久久久精品古装| 男女无遮挡免费网站观看| 看十八女毛片水多多多| 香蕉国产在线看| 国产亚洲欧美精品永久| 99久国产av精品国产电影| 飞空精品影院首页| 宅男免费午夜| 亚洲第一区二区三区不卡| 亚洲欧美成人综合另类久久久| 久久久久精品国产欧美久久久 | 亚洲激情五月婷婷啪啪| 久久久国产一区二区| 久久99热这里只频精品6学生| 好男人视频免费观看在线| 青春草国产在线视频| 国产色婷婷99| 国产成人精品在线电影| 国产精品成人在线| 夫妻午夜视频| 一本久久精品| 18禁动态无遮挡网站| 国产又爽黄色视频| 亚洲av日韩精品久久久久久密 | 五月开心婷婷网| 日韩人妻精品一区2区三区| 最黄视频免费看| 久久久精品94久久精品| 亚洲成人一二三区av| 国产伦理片在线播放av一区| 国产一区二区三区av在线| av一本久久久久| 日韩一区二区三区影片| 欧美变态另类bdsm刘玥| 天天操日日干夜夜撸| 国产国语露脸激情在线看| 伦理电影免费视频| 美女国产高潮福利片在线看| tube8黄色片| 日韩精品免费视频一区二区三区| 国产av码专区亚洲av| 日韩不卡一区二区三区视频在线| 国产成人免费无遮挡视频| 国产欧美日韩一区二区三区在线| 亚洲人成网站在线观看播放| 男女边摸边吃奶| 日韩大片免费观看网站| 亚洲欧美清纯卡通| 亚洲少妇的诱惑av| 久久人妻熟女aⅴ| 99精品久久久久人妻精品| 久久午夜综合久久蜜桃| 少妇 在线观看| 国产精品久久久av美女十八| 色网站视频免费| 制服诱惑二区| 搡老乐熟女国产| 亚洲欧美中文字幕日韩二区| 亚洲精品久久成人aⅴ小说| 51午夜福利影视在线观看| 久久久久精品久久久久真实原创| 91精品伊人久久大香线蕉| 女人被躁到高潮嗷嗷叫费观| 伊人久久国产一区二区| 国产av码专区亚洲av| 精品一区二区免费观看| 我要看黄色一级片免费的| 丝瓜视频免费看黄片| 免费在线观看完整版高清| 久久精品国产亚洲av高清一级| 亚洲av成人精品一二三区| 日本av手机在线免费观看| av女优亚洲男人天堂| 91成人精品电影| 午夜福利在线免费观看网站| 一本色道久久久久久精品综合| av国产久精品久网站免费入址| 如何舔出高潮| 不卡视频在线观看欧美| 69精品国产乱码久久久| 久久韩国三级中文字幕| 狠狠精品人妻久久久久久综合| 国产成人午夜福利电影在线观看| 一级,二级,三级黄色视频| 成年女人毛片免费观看观看9 | 久久久久人妻精品一区果冻| 国产黄色视频一区二区在线观看| h视频一区二区三区| 看免费成人av毛片| 久久天躁狠狠躁夜夜2o2o | 亚洲国产精品一区三区| 欧美日韩av久久| 日日啪夜夜爽| 在线天堂最新版资源| 69精品国产乱码久久久| 日日摸夜夜添夜夜爱| 亚洲熟女毛片儿| 男女无遮挡免费网站观看| 亚洲成人av在线免费| 男人操女人黄网站| 国产成人av激情在线播放| 亚洲欧美中文字幕日韩二区| 久热这里只有精品99| 叶爱在线成人免费视频播放| 桃花免费在线播放| netflix在线观看网站| 飞空精品影院首页| 欧美激情极品国产一区二区三区| 婷婷色综合大香蕉| netflix在线观看网站| 国产爽快片一区二区三区| 国产伦理片在线播放av一区| 最黄视频免费看| 一本久久精品| 亚洲欧美色中文字幕在线| 久久精品人人爽人人爽视色| 免费不卡黄色视频| 亚洲伊人久久精品综合| 最近中文字幕高清免费大全6| 国产日韩欧美视频二区| 亚洲精品一二三| 亚洲精品视频女| 多毛熟女@视频| 午夜激情av网站| 久久狼人影院| 伊人亚洲综合成人网| www.精华液| 女人爽到高潮嗷嗷叫在线视频| videos熟女内射| 波野结衣二区三区在线| 爱豆传媒免费全集在线观看| 97人妻天天添夜夜摸| 日韩电影二区| 久久久亚洲精品成人影院| svipshipincom国产片| 国产精品99久久99久久久不卡 | 国产精品人妻久久久影院| 男男h啪啪无遮挡| 纯流量卡能插随身wifi吗| 日韩,欧美,国产一区二区三区| 青草久久国产| 高清不卡的av网站| 国产在线一区二区三区精| 捣出白浆h1v1| 中文欧美无线码| 亚洲一级一片aⅴ在线观看| 高清av免费在线| 老司机在亚洲福利影院| 80岁老熟妇乱子伦牲交| 街头女战士在线观看网站| 青青草视频在线视频观看| 黄网站色视频无遮挡免费观看| 日韩欧美一区视频在线观看| 亚洲av国产av综合av卡| 精品少妇久久久久久888优播| 老司机靠b影院| 九草在线视频观看| 欧美日韩一级在线毛片| 色播在线永久视频| 国产免费又黄又爽又色| 高清黄色对白视频在线免费看| 十分钟在线观看高清视频www| 中文字幕精品免费在线观看视频| 国产免费又黄又爽又色| 1024视频免费在线观看| 人人妻人人添人人爽欧美一区卜| 久久久久久人妻| 国产av精品麻豆| e午夜精品久久久久久久| 夜夜骑夜夜射夜夜干| 中文字幕制服av| 黄色一级大片看看| 大码成人一级视频| 亚洲精品在线美女| 啦啦啦 在线观看视频| bbb黄色大片| 国语对白做爰xxxⅹ性视频网站| 欧美av亚洲av综合av国产av | 久久韩国三级中文字幕| 97在线人人人人妻| 黄频高清免费视频| www.精华液| 国产 精品1| 午夜免费观看性视频| 亚洲精品日韩在线中文字幕| 啦啦啦中文免费视频观看日本| 亚洲av在线观看美女高潮| 精品国产一区二区三区久久久樱花| 中文字幕色久视频| 欧美成人午夜精品| 国产精品.久久久| 操美女的视频在线观看| 亚洲精品日本国产第一区| 另类亚洲欧美激情| 又大又黄又爽视频免费| 欧美 亚洲 国产 日韩一| 一个人免费看片子| 亚洲国产欧美一区二区综合| av国产精品久久久久影院| 亚洲美女搞黄在线观看| 黑丝袜美女国产一区| 啦啦啦在线观看免费高清www| 日韩 亚洲 欧美在线| 黄片小视频在线播放| 69精品国产乱码久久久| 免费在线观看视频国产中文字幕亚洲 | 最近的中文字幕免费完整| 免费看不卡的av| 国产精品无大码| 亚洲,一卡二卡三卡| 伊人久久国产一区二区| 亚洲一区二区三区欧美精品| 精品国产国语对白av| 亚洲第一av免费看| 丝袜人妻中文字幕| 91精品伊人久久大香线蕉| 女人高潮潮喷娇喘18禁视频| 母亲3免费完整高清在线观看| 亚洲,一卡二卡三卡| 天天影视国产精品| 久久性视频一级片| 人人澡人人妻人| 亚洲专区中文字幕在线 | 国产成人a∨麻豆精品| av在线老鸭窝| 欧美精品一区二区免费开放| av电影中文网址| 亚洲人成网站在线观看播放| 国产欧美日韩综合在线一区二区| 咕卡用的链子| avwww免费| 成人黄色视频免费在线看| 毛片一级片免费看久久久久| 热99久久久久精品小说推荐| 在线观看免费午夜福利视频| 一二三四在线观看免费中文在| 亚洲精品久久成人aⅴ小说| 国产精品一国产av| 老司机亚洲免费影院| 老熟女久久久| 国产精品国产三级专区第一集| 9191精品国产免费久久| 欧美变态另类bdsm刘玥| 国产97色在线日韩免费| 国产一级毛片在线| 最近中文字幕2019免费版| 国产男人的电影天堂91| 国产精品一区二区精品视频观看| 老汉色av国产亚洲站长工具| 成年女人毛片免费观看观看9 | 色综合欧美亚洲国产小说| 中文乱码字字幕精品一区二区三区| 亚洲国产精品999| 日韩制服骚丝袜av| 精品少妇一区二区三区视频日本电影 | 欧美av亚洲av综合av国产av | 亚洲婷婷狠狠爱综合网| 国产精品久久久av美女十八| 婷婷色综合大香蕉| 日本欧美国产在线视频| 高清视频免费观看一区二区| 亚洲免费av在线视频| 老司机在亚洲福利影院| 美女高潮到喷水免费观看| 亚洲精品国产一区二区精华液| 亚洲精品中文字幕在线视频| 欧美xxⅹ黑人| 我的亚洲天堂| 波多野结衣av一区二区av| 亚洲精品自拍成人| www.熟女人妻精品国产| 晚上一个人看的免费电影| 久久女婷五月综合色啪小说| 欧美人与性动交α欧美精品济南到| 啦啦啦在线免费观看视频4| 欧美在线黄色| 国产精品偷伦视频观看了| 嫩草影院入口| 18禁国产床啪视频网站| 国产一区亚洲一区在线观看| 黄网站色视频无遮挡免费观看| 亚洲第一青青草原| 成年女人毛片免费观看观看9 | 宅男免费午夜| 极品人妻少妇av视频| 国产精品免费视频内射| videosex国产| 久久ye,这里只有精品| 天天躁夜夜躁狠狠躁躁| 亚洲av欧美aⅴ国产| 国产一区二区 视频在线| 十八禁高潮呻吟视频| 成人毛片60女人毛片免费| 国产亚洲最大av| 久久久久久久大尺度免费视频| 99精国产麻豆久久婷婷| 亚洲国产中文字幕在线视频| 国产亚洲av高清不卡| 成人亚洲欧美一区二区av| 激情五月婷婷亚洲| 九九爱精品视频在线观看| 美女大奶头黄色视频| 久久97久久精品| 国产又爽黄色视频| 亚洲色图 男人天堂 中文字幕| 热re99久久国产66热| 久久国产精品男人的天堂亚洲| 91国产中文字幕| 欧美乱码精品一区二区三区| kizo精华| 女人精品久久久久毛片| 老司机深夜福利视频在线观看 | 日韩成人av中文字幕在线观看| 国产亚洲av片在线观看秒播厂| 免费人妻精品一区二区三区视频| 亚洲欧美一区二区三区久久| 国产午夜精品一二区理论片| 十分钟在线观看高清视频www| 国产精品国产三级国产专区5o| 亚洲av日韩精品久久久久久密 | 精品一区二区三区四区五区乱码 | 亚洲一区二区三区欧美精品| 婷婷色综合www| 黄色一级大片看看| 欧美国产精品一级二级三级| 黑人巨大精品欧美一区二区蜜桃| 国产一区二区三区av在线| 亚洲一区中文字幕在线| 午夜福利免费观看在线| 日本爱情动作片www.在线观看| 日韩电影二区| 亚洲av男天堂| 狠狠精品人妻久久久久久综合| 这个男人来自地球电影免费观看 | 亚洲人成网站在线观看播放| 亚洲国产av新网站| 超碰97精品在线观看| 国产成人精品无人区| 久久97久久精品| 1024香蕉在线观看| 日韩中文字幕视频在线看片| 亚洲伊人色综图| 美女主播在线视频| 亚洲欧美激情在线| 亚洲国产精品一区二区三区在线| 精品卡一卡二卡四卡免费| 高清视频免费观看一区二区| 在线观看国产h片| 国产精品香港三级国产av潘金莲 | 日本猛色少妇xxxxx猛交久久| 亚洲图色成人| 99精品久久久久人妻精品| 久久久久久久大尺度免费视频| 久久性视频一级片| 久久久久国产精品人妻一区二区| 成人午夜精彩视频在线观看| 久久久国产精品麻豆| bbb黄色大片| 久久久久久久精品精品| 国产精品嫩草影院av在线观看| 少妇精品久久久久久久| 午夜日本视频在线| 日日爽夜夜爽网站| 国产精品久久久人人做人人爽| 香蕉国产在线看| 亚洲欧美一区二区三区久久| 国产淫语在线视频| 日本av免费视频播放| 人妻人人澡人人爽人人| 亚洲视频免费观看视频| 少妇人妻 视频| 水蜜桃什么品种好| 老汉色av国产亚洲站长工具| 一区二区三区激情视频| 亚洲欧美成人精品一区二区| 黑人欧美特级aaaaaa片| 深夜精品福利| 91精品伊人久久大香线蕉| 国产一卡二卡三卡精品 | av国产久精品久网站免费入址| 人人妻人人爽人人添夜夜欢视频| 搡老乐熟女国产| 国产色婷婷99| xxx大片免费视频| 少妇精品久久久久久久| 国产精品亚洲av一区麻豆 | 大陆偷拍与自拍| 亚洲精品av麻豆狂野| 亚洲国产日韩一区二区| 国产日韩一区二区三区精品不卡| 热99久久久久精品小说推荐| 激情五月婷婷亚洲| 青春草视频在线免费观看| 精品久久蜜臀av无| 久久国产精品男人的天堂亚洲| 国产成人免费观看mmmm| 亚洲国产精品国产精品| 免费不卡黄色视频| 麻豆乱淫一区二区| 母亲3免费完整高清在线观看| 美女国产高潮福利片在线看| 交换朋友夫妻互换小说| 免费在线观看视频国产中文字幕亚洲 | 国精品久久久久久国模美| 在线精品无人区一区二区三| 中国国产av一级| 久久狼人影院| 欧美日韩亚洲国产一区二区在线观看 | 国产一区二区三区av在线| 老司机在亚洲福利影院| 中文字幕人妻丝袜制服| 咕卡用的链子| 免费高清在线观看日韩| 亚洲美女视频黄频| 中文字幕制服av| 一级片免费观看大全| 高清不卡的av网站| 中文字幕人妻丝袜一区二区 | 激情视频va一区二区三区| 日韩一本色道免费dvd| 我要看黄色一级片免费的| 亚洲av日韩精品久久久久久密 | 亚洲精品久久久久久婷婷小说| 老司机深夜福利视频在线观看 | 亚洲精品国产av蜜桃| 如日韩欧美国产精品一区二区三区| 高清欧美精品videossex| 美女午夜性视频免费| 日韩大片免费观看网站| 老司机靠b影院| 国产伦人伦偷精品视频| 精品久久久精品久久久| 人人妻人人澡人人爽人人夜夜| 一区二区三区四区激情视频| 欧美av亚洲av综合av国产av | 国语对白做爰xxxⅹ性视频网站| 久久久久久久久久久免费av| 亚洲色图综合在线观看| 人人妻人人爽人人添夜夜欢视频| 七月丁香在线播放| 午夜日韩欧美国产| 亚洲精品一二三| 欧美日韩福利视频一区二区| av女优亚洲男人天堂| 一本一本久久a久久精品综合妖精| 99热国产这里只有精品6| 老司机影院成人| 看十八女毛片水多多多| 啦啦啦中文免费视频观看日本| 亚洲视频免费观看视频| 久久这里只有精品19| 最近手机中文字幕大全| 丝袜美腿诱惑在线| 国产精品三级大全| 国产精品香港三级国产av潘金莲 | 成年女人毛片免费观看观看9 | 国产精品欧美亚洲77777| 18禁动态无遮挡网站| 老鸭窝网址在线观看| 久久精品久久精品一区二区三区| 人人妻人人爽人人添夜夜欢视频| 亚洲伊人久久精品综合| 999久久久国产精品视频| 欧美精品人与动牲交sv欧美| 国产老妇伦熟女老妇高清| 深夜精品福利| 亚洲精品国产色婷婷电影| 日本91视频免费播放| 菩萨蛮人人尽说江南好唐韦庄| 999精品在线视频| 色精品久久人妻99蜜桃| 精品卡一卡二卡四卡免费| 欧美激情极品国产一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 国产日韩欧美亚洲二区| 在线观看免费日韩欧美大片| 国产无遮挡羞羞视频在线观看| 国产黄色免费在线视频| 精品卡一卡二卡四卡免费| 人体艺术视频欧美日本| 久久久久精品国产欧美久久久 | 色精品久久人妻99蜜桃| 亚洲一区中文字幕在线| 久久精品人人爽人人爽视色| 在线观看免费高清a一片| 欧美变态另类bdsm刘玥| 亚洲欧洲精品一区二区精品久久久 | 91成人精品电影| 国产欧美日韩综合在线一区二区| 天堂俺去俺来也www色官网| 青草久久国产| 一边摸一边抽搐一进一出视频| 欧美xxⅹ黑人| 国精品久久久久久国模美| 国产在线视频一区二区| 亚洲av电影在线观看一区二区三区| 高清在线视频一区二区三区| 日韩不卡一区二区三区视频在线| 免费黄色在线免费观看| 精品人妻熟女毛片av久久网站| 国产精品免费大片| 欧美在线黄色| 亚洲美女视频黄频| 飞空精品影院首页| 国产黄色免费在线视频| 久久人人爽人人片av| 免费女性裸体啪啪无遮挡网站| 热re99久久国产66热| 午夜免费男女啪啪视频观看| 中文字幕色久视频| 亚洲三区欧美一区| 国产精品偷伦视频观看了| 久久精品国产综合久久久| videosex国产| 国产视频首页在线观看| 亚洲视频免费观看视频| 亚洲av福利一区| 丁香六月天网| 亚洲欧洲日产国产| 啦啦啦啦在线视频资源| 亚洲欧美精品综合一区二区三区| 久久精品人人爽人人爽视色| 久久午夜综合久久蜜桃| 亚洲第一av免费看| 色综合欧美亚洲国产小说|