• <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一区麻豆| 真人一进一出gif抽搐免费| 日本 欧美在线| 日韩视频一区二区在线观看| 国产成年人精品一区二区| 午夜影院日韩av| 亚洲专区国产一区二区| 国产精品野战在线观看| 婷婷丁香在线五月| 麻豆成人午夜福利视频| 日韩欧美精品免费久久| 精品久久国产蜜桃| 亚洲自偷自拍三级| 国产91精品成人一区二区三区| 日韩一本色道免费dvd| 欧美不卡视频在线免费观看| 99久久精品一区二区三区| 啦啦啦观看免费观看视频高清| 啦啦啦啦在线视频资源| 亚洲四区av| 一边摸一边抽搐一进一小说| 一区二区三区四区激情视频 | 亚洲第一电影网av| 亚洲aⅴ乱码一区二区在线播放| 日韩欧美 国产精品| 久久精品国产亚洲av涩爱 | 国产精品久久久久久亚洲av鲁大| 国产精品乱码一区二三区的特点| 一个人看视频在线观看www免费| 精品久久久久久久久久免费视频| 性欧美人与动物交配| 精品午夜福利视频在线观看一区| 欧美色欧美亚洲另类二区| 99热这里只有是精品在线观看| 国产视频内射| 欧美丝袜亚洲另类 | 给我免费播放毛片高清在线观看| 久久6这里有精品| or卡值多少钱| 男女之事视频高清在线观看| 成人特级av手机在线观看| 亚洲熟妇中文字幕五十中出| 亚洲av成人av| 少妇的逼好多水| 成年女人永久免费观看视频| 欧美一区二区亚洲| 亚洲欧美日韩卡通动漫| 中亚洲国语对白在线视频| 啪啪无遮挡十八禁网站| 国产成人影院久久av| 国产人妻一区二区三区在| 可以在线观看的亚洲视频| 在线免费观看不下载黄p国产 | 亚洲一区二区三区色噜噜| 免费av不卡在线播放| 免费观看在线日韩| 亚洲国产高清在线一区二区三| 日本成人三级电影网站| av福利片在线观看| 国产高清不卡午夜福利| 全区人妻精品视频| 免费在线观看影片大全网站| 国产女主播在线喷水免费视频网站 | 亚洲av.av天堂| 欧美日韩精品成人综合77777| 美女免费视频网站| av天堂在线播放| 2021天堂中文幕一二区在线观| 五月伊人婷婷丁香| 亚洲人成网站在线播放欧美日韩| 亚洲真实伦在线观看| 校园人妻丝袜中文字幕| 一级av片app| 两个人的视频大全免费| 精品福利观看| 亚洲av五月六月丁香网| 国产精品98久久久久久宅男小说| 国产老妇女一区| 一级黄色大片毛片| 午夜视频国产福利| 亚洲性久久影院| 少妇高潮的动态图| 国产精品久久久久久久电影| 国产精品1区2区在线观看.| 又紧又爽又黄一区二区| 99久久九九国产精品国产免费| 免费av不卡在线播放| 国产高潮美女av| 国产精品福利在线免费观看| 美女大奶头视频| 久久99热这里只有精品18| 国产91精品成人一区二区三区| 少妇的逼好多水| 成年女人永久免费观看视频| 日本a在线网址| 丝袜美腿在线中文| 99在线视频只有这里精品首页| 欧美成人一区二区免费高清观看| 国产欧美日韩一区二区精品| 狂野欧美白嫩少妇大欣赏| 男人狂女人下面高潮的视频| 欧美成人免费av一区二区三区| 亚洲av电影不卡..在线观看| 老师上课跳d突然被开到最大视频| 欧美最黄视频在线播放免费| 在线观看舔阴道视频| 亚洲经典国产精华液单| 国产精品一区二区性色av| 国产乱人视频| 18禁黄网站禁片免费观看直播| 欧美bdsm另类| 欧美bdsm另类| 亚洲一级一片aⅴ在线观看| 亚洲熟妇熟女久久| 变态另类丝袜制服| 欧美性猛交黑人性爽| 男女之事视频高清在线观看| 在线观看午夜福利视频| 午夜爱爱视频在线播放| 性欧美人与动物交配| 日韩欧美国产在线观看| 干丝袜人妻中文字幕| 别揉我奶头 嗯啊视频| 国产成人a区在线观看| 十八禁国产超污无遮挡网站| 中文在线观看免费www的网站| 精品久久久噜噜| 18+在线观看网站| 日韩,欧美,国产一区二区三区 | 观看免费一级毛片| 精品国内亚洲2022精品成人| 在线播放无遮挡| 亚洲国产精品合色在线| 亚洲专区国产一区二区| 国产精品综合久久久久久久免费| 亚洲自拍偷在线| 国产私拍福利视频在线观看| 中文字幕熟女人妻在线| 国产高清有码在线观看视频| 真人一进一出gif抽搐免费| 中文字幕精品亚洲无线码一区| 欧美极品一区二区三区四区| 日本一本二区三区精品| 在线观看舔阴道视频| 国产黄片美女视频| 亚洲成人久久性| 精品不卡国产一区二区三区| 97碰自拍视频| 他把我摸到了高潮在线观看| 91久久精品国产一区二区三区| 毛片一级片免费看久久久久 | 天堂av国产一区二区熟女人妻| 少妇熟女aⅴ在线视频| 国产乱人视频| 日本 av在线| 有码 亚洲区| 亚洲av一区综合| 午夜福利视频1000在线观看| 亚洲一区高清亚洲精品| 精品久久久久久久久久免费视频| 搞女人的毛片| 永久网站在线| 国产不卡一卡二| 亚洲内射少妇av| 免费看美女性在线毛片视频| 国产人妻一区二区三区在| 亚洲avbb在线观看| 看片在线看免费视频| 国产精品一区二区三区四区免费观看 | xxxwww97欧美| 中文在线观看免费www的网站| 成人永久免费在线观看视频| 夜夜看夜夜爽夜夜摸| 国产精品福利在线免费观看| 久久久久久久久大av| 18禁在线播放成人免费| 国产精华一区二区三区| 成年女人毛片免费观看观看9| 国产综合懂色| 欧美日本视频| 熟女电影av网| 欧美日韩精品成人综合77777| 亚洲熟妇熟女久久| av视频在线观看入口| 熟女人妻精品中文字幕| 精品日产1卡2卡| 美女被艹到高潮喷水动态| 国产精品国产高清国产av| 国产爱豆传媒在线观看| 欧美bdsm另类| 国国产精品蜜臀av免费| 日本一二三区视频观看| 一卡2卡三卡四卡精品乱码亚洲| 国产精品人妻久久久久久| 三级男女做爰猛烈吃奶摸视频| 亚洲美女黄片视频| 亚洲国产欧美人成| 久久久久九九精品影院| 久久草成人影院| 国产伦精品一区二区三区四那| 国产毛片a区久久久久| 国产一区二区三区视频了| 欧美bdsm另类| 舔av片在线| 给我免费播放毛片高清在线观看| 99热这里只有精品一区| 在线播放无遮挡| av在线蜜桃| 最新中文字幕久久久久| 天堂动漫精品| 一个人观看的视频www高清免费观看| 午夜精品久久久久久毛片777| 一本精品99久久精品77| 少妇熟女aⅴ在线视频| 此物有八面人人有两片| 99在线人妻在线中文字幕| 男女视频在线观看网站免费| 麻豆av噜噜一区二区三区| 可以在线观看毛片的网站| av天堂中文字幕网| 国产精品无大码| 亚洲国产欧美人成| 日韩亚洲欧美综合| 日韩欧美在线二视频| 国产精品亚洲美女久久久| 色哟哟哟哟哟哟| 免费av毛片视频| 啪啪无遮挡十八禁网站| 女人被狂操c到高潮| 色精品久久人妻99蜜桃| 欧美高清性xxxxhd video| 久久精品国产鲁丝片午夜精品 | 国产精品,欧美在线| 亚洲国产精品成人综合色| 成人二区视频| 国内精品久久久久精免费| 亚洲国产高清在线一区二区三| 在线观看一区二区三区| 乱人视频在线观看| 国产淫片久久久久久久久| 成人亚洲精品av一区二区| 少妇熟女aⅴ在线视频| 精品人妻1区二区| 国国产精品蜜臀av免费| 男人和女人高潮做爰伦理| 日本三级黄在线观看| 国产伦在线观看视频一区| 免费在线观看影片大全网站| 99热网站在线观看| 毛片一级片免费看久久久久 | 特级一级黄色大片| 人人妻,人人澡人人爽秒播| 十八禁网站免费在线| 午夜福利视频1000在线观看| 国产精品永久免费网站| 欧美人与善性xxx| 国产av在哪里看| 日韩强制内射视频| 99热6这里只有精品| 窝窝影院91人妻| 国产午夜福利久久久久久| 99热只有精品国产| a在线观看视频网站| 99热这里只有精品一区| 亚洲欧美日韩东京热| 国产男靠女视频免费网站| a级毛片免费高清观看在线播放| 99热6这里只有精品| 联通29元200g的流量卡| 国产一区二区三区视频了| 我的女老师完整版在线观看| 99久久无色码亚洲精品果冻| 91精品国产九色| 内射极品少妇av片p| 成人三级黄色视频| 69av精品久久久久久| 久久久久精品国产欧美久久久| 三级男女做爰猛烈吃奶摸视频| 久久精品久久久久久噜噜老黄 | 嫁个100分男人电影在线观看| 级片在线观看| 又紧又爽又黄一区二区| 99热只有精品国产| 给我免费播放毛片高清在线观看| 国产 一区 欧美 日韩| 老熟妇乱子伦视频在线观看| 婷婷色综合大香蕉| 极品教师在线免费播放| 国产视频内射| 不卡视频在线观看欧美| 成年女人毛片免费观看观看9| 窝窝影院91人妻| 男女啪啪激烈高潮av片| 欧美性猛交黑人性爽| 亚洲成人久久爱视频| 真人做人爱边吃奶动态| 日韩强制内射视频| 97热精品久久久久久| 中文字幕人妻熟人妻熟丝袜美| 九九在线视频观看精品| 久久久久久伊人网av| 18禁裸乳无遮挡免费网站照片| 亚洲第一电影网av| 成人特级黄色片久久久久久久| 日日撸夜夜添| 日本欧美国产在线视频| 亚洲欧美日韩高清在线视频| 欧美三级亚洲精品| 久久人人爽人人爽人人片va| 乱系列少妇在线播放| 成人永久免费在线观看视频| 最新中文字幕久久久久| 国产精华一区二区三区| 啦啦啦观看免费观看视频高清| 午夜福利高清视频| 九九热线精品视视频播放| 一级黄片播放器| 国产精品99久久久久久久久| 欧美日韩精品成人综合77777| 欧美成人一区二区免费高清观看| 国产精品一区二区免费欧美| 国产高潮美女av| 最近视频中文字幕2019在线8| 91精品国产九色| 亚洲成人免费电影在线观看| 床上黄色一级片| 国产69精品久久久久777片| 欧美日韩瑟瑟在线播放| 国产精品久久久久久精品电影| 久久亚洲真实| 国产在视频线在精品| 国内少妇人妻偷人精品xxx网站| 我要看日韩黄色一级片| 亚洲一区高清亚洲精品| 免费无遮挡裸体视频| 日韩欧美三级三区| 成人av在线播放网站| 国产伦在线观看视频一区| 人人妻,人人澡人人爽秒播| 亚洲七黄色美女视频| 国产色爽女视频免费观看| 欧美一区二区精品小视频在线| or卡值多少钱| 一卡2卡三卡四卡精品乱码亚洲| 国产熟女欧美一区二区| 少妇的逼好多水| 国产视频内射| 欧美激情久久久久久爽电影| 国产精品免费一区二区三区在线| 亚洲欧美日韩高清专用| 尾随美女入室| 国产亚洲精品久久久久久毛片| 日韩中字成人| 嫩草影院入口| 人妻丰满熟妇av一区二区三区| 搡老岳熟女国产| 日本a在线网址| 成人鲁丝片一二三区免费| 欧美精品啪啪一区二区三区| 一进一出好大好爽视频| 亚洲国产色片| 国产探花极品一区二区| 国产高清激情床上av| 久久午夜亚洲精品久久| 在线观看舔阴道视频| 精品国产三级普通话版| 真人一进一出gif抽搐免费| 国产精品99久久久久久久久| 婷婷六月久久综合丁香| 简卡轻食公司| 欧美绝顶高潮抽搐喷水| 香蕉av资源在线| 国产成人av教育| av福利片在线观看| a级毛片免费高清观看在线播放| 精品午夜福利视频在线观看一区| 熟女电影av网| 我要看日韩黄色一级片| 级片在线观看| 男人舔奶头视频| 亚洲精华国产精华精| 日日摸夜夜添夜夜添小说| 最新中文字幕久久久久| 麻豆国产av国片精品| 国产免费男女视频| 我的女老师完整版在线观看| 国产中年淑女户外野战色| 三级男女做爰猛烈吃奶摸视频| 国产精品三级大全| 国产精品久久久久久久久免| 在线天堂最新版资源| 亚洲狠狠婷婷综合久久图片| 欧美成人免费av一区二区三区| 亚洲人成网站在线播放欧美日韩| 狂野欧美激情性xxxx在线观看| 国产精品久久电影中文字幕| 日韩欧美 国产精品| 一个人免费在线观看电影| 久久人妻av系列| 国产午夜精品久久久久久一区二区三区 | 日韩欧美国产在线观看| 神马国产精品三级电影在线观看| 一级毛片久久久久久久久女| av在线老鸭窝| 最好的美女福利视频网| 欧美zozozo另类| 老女人水多毛片| 少妇人妻一区二区三区视频| 毛片女人毛片| 久久精品国产亚洲av涩爱 | 亚洲在线观看片| 亚洲国产精品sss在线观看| 春色校园在线视频观看| 亚洲,欧美,日韩| 深夜精品福利| 女人被狂操c到高潮| 国产一级毛片七仙女欲春2| 我要看日韩黄色一级片| 国产av一区在线观看免费| 看十八女毛片水多多多| 婷婷六月久久综合丁香| 国产精品嫩草影院av在线观看 | 亚洲av五月六月丁香网| 国产男人的电影天堂91| 91精品国产九色| 色哟哟哟哟哟哟| 精品免费久久久久久久清纯| 国内精品久久久久精免费| 婷婷色综合大香蕉| 乱系列少妇在线播放| 婷婷亚洲欧美| 亚洲无线在线观看| 色尼玛亚洲综合影院| 中文字幕av成人在线电影| 久久天躁狠狠躁夜夜2o2o| 亚洲午夜理论影院| 免费av毛片视频| 两人在一起打扑克的视频| 色哟哟哟哟哟哟| 一区二区三区四区激情视频 | 国产一区二区亚洲精品在线观看| 日本在线视频免费播放| 成人av一区二区三区在线看| 欧美成人性av电影在线观看| 欧美日韩乱码在线| 久久久久久久久久久丰满 | 天堂动漫精品| 精品午夜福利视频在线观看一区| 在线观看av片永久免费下载| 毛片一级片免费看久久久久 | 亚洲精品在线观看二区| 桃红色精品国产亚洲av| 亚洲最大成人中文| 亚洲熟妇熟女久久| 国内久久婷婷六月综合欲色啪| 2021天堂中文幕一二区在线观| 男人的好看免费观看在线视频| 亚洲av中文字字幕乱码综合| 国产私拍福利视频在线观看| 成人二区视频| 国产一区二区三区视频了| 99国产精品一区二区蜜桃av| a级一级毛片免费在线观看| 国产伦精品一区二区三区视频9| 国产亚洲av嫩草精品影院| 特级一级黄色大片| 麻豆国产97在线/欧美| 中文字幕人妻熟人妻熟丝袜美| 欧美日韩精品成人综合77777| 在线免费十八禁| 国产成年人精品一区二区| bbb黄色大片| 欧美日韩黄片免| 99久久无色码亚洲精品果冻| 国产伦一二天堂av在线观看| 国产精品日韩av在线免费观看| 亚洲第一电影网av| 少妇人妻一区二区三区视频| 少妇人妻精品综合一区二区 | 精品午夜福利在线看| 欧美国产日韩亚洲一区| 美女大奶头视频| 国产精品美女特级片免费视频播放器| 夜夜看夜夜爽夜夜摸| 久久亚洲真实| 国产伦一二天堂av在线观看| 成年女人永久免费观看视频| 淫妇啪啪啪对白视频| 久久这里只有精品中国| 亚洲成a人片在线一区二区| 国产免费一级a男人的天堂| 男人的好看免费观看在线视频| 国语自产精品视频在线第100页| 91久久精品国产一区二区三区| 色尼玛亚洲综合影院| 国产女主播在线喷水免费视频网站 | 中文字幕av在线有码专区| 又紧又爽又黄一区二区| 国产伦精品一区二区三区视频9| 两人在一起打扑克的视频| 麻豆久久精品国产亚洲av| 免费人成视频x8x8入口观看| 一个人看视频在线观看www免费| 中文字幕熟女人妻在线| 国产精品野战在线观看| 久久久久久久久久黄片| 九九爱精品视频在线观看| 国产精品久久久久久亚洲av鲁大| 久9热在线精品视频| 欧美最新免费一区二区三区| 两个人视频免费观看高清| 91狼人影院| 亚洲内射少妇av| 黄色欧美视频在线观看| 99热只有精品国产| 精品99又大又爽又粗少妇毛片 | 精品一区二区三区视频在线| 国产黄片美女视频| 欧美一区二区亚洲| 亚洲中文日韩欧美视频| 亚洲内射少妇av| 亚洲三级黄色毛片| 国产成年人精品一区二区| 级片在线观看| 日韩精品有码人妻一区| 久久草成人影院| 亚洲精品成人久久久久久| 国产中年淑女户外野战色| 欧美三级亚洲精品| 性插视频无遮挡在线免费观看| 一个人看视频在线观看www免费| 深爱激情五月婷婷| 日日摸夜夜添夜夜添av毛片 | 夜夜夜夜夜久久久久| 国产精品久久久久久亚洲av鲁大| 成人av在线播放网站| 国产高清视频在线播放一区| 亚洲最大成人手机在线| 国产精品久久久久久久电影| 国产精品久久久久久av不卡| 少妇高潮的动态图| 色播亚洲综合网| 国产白丝娇喘喷水9色精品| а√天堂www在线а√下载| 日本熟妇午夜| 欧美日韩乱码在线| 国产精品一区二区三区四区免费观看 | 亚洲av中文av极速乱 | 亚洲电影在线观看av| 欧美成人a在线观看| 亚洲色图av天堂| а√天堂www在线а√下载| 欧美日韩瑟瑟在线播放| 丝袜美腿在线中文| 国产精品1区2区在线观看.| 天天躁日日操中文字幕| 国产精品国产高清国产av| 制服丝袜大香蕉在线| 成人永久免费在线观看视频| 国产高清有码在线观看视频| 少妇熟女aⅴ在线视频| 91久久精品国产一区二区成人| av在线亚洲专区| 俄罗斯特黄特色一大片| 国产免费男女视频| 一个人免费在线观看电影| 特级一级黄色大片| 超碰av人人做人人爽久久| 一个人看视频在线观看www免费| 一个人免费在线观看电影| 少妇熟女aⅴ在线视频| 免费人成视频x8x8入口观看| .国产精品久久| 99riav亚洲国产免费| 亚洲第一区二区三区不卡| av女优亚洲男人天堂| 日韩欧美国产一区二区入口| xxxwww97欧美| 亚洲成人中文字幕在线播放| 久久精品国产亚洲av天美| 亚洲综合色惰| 黄色欧美视频在线观看| 露出奶头的视频| 内射极品少妇av片p| 天美传媒精品一区二区| 国产白丝娇喘喷水9色精品| 嫁个100分男人电影在线观看| 免费在线观看影片大全网站| 国产中年淑女户外野战色| 少妇人妻一区二区三区视频| 内地一区二区视频在线| 99久久中文字幕三级久久日本| 亚洲三级黄色毛片| 久久久久久大精品| 搡老熟女国产l中国老女人| 乱人视频在线观看| 国产精品国产高清国产av| 高清毛片免费观看视频网站| 我的老师免费观看完整版| 亚洲欧美激情综合另类| 亚洲精品成人久久久久久| 精品人妻一区二区三区麻豆 |