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

    球面波PP反射系數(shù)的頻變特征研究

    2016-11-08 03:04:04李京南王尚旭董春暉袁三一王靜波
    地球物理學(xué)報 2016年10期
    關(guān)鍵詞:平面波反射系數(shù)入射角

    李京南, 王尚旭*, 董春暉, 袁三一, 王靜波

    1 中國石油大學(xué)(北京)油氣資源與探測國家重點實驗室, 北京 102249 2 中國石油大學(xué)(北京)CNPC物探重點實驗室, 北京 102249 3 中國石油化工股份有限公司勘探分公司, 成都 610041

    ?

    球面波PP反射系數(shù)的頻變特征研究

    李京南1,2, 王尚旭1,2*, 董春暉1,2, 袁三一1,2, 王靜波3

    1 中國石油大學(xué)(北京)油氣資源與探測國家重點實驗室, 北京102249 2 中國石油大學(xué)(北京)CNPC物探重點實驗室, 北京102249 3 中國石油化工股份有限公司勘探分公司, 成都610041

    與平面波反射系數(shù)相比,球面波反射系數(shù)可以更精確地描述實際地震波的反射特征.近些年關(guān)于球面波的研究主要聚焦于球面波反射系數(shù)隨入射角的變化規(guī)律,很少對球面波反射系數(shù)隨頻率的變化(頻變)做詳細研究.為了更全面地了解球面波的反射機制,本文研究了兩層彈性介質(zhì)中球面波PP反射系數(shù)(幅值和相位)的頻變規(guī)律.文中基于經(jīng)典的Sommerfeld積分構(gòu)造球面波PP反射系數(shù),通過自適應(yīng)的Gauss-Kronrod求積算法對其進行穩(wěn)定的數(shù)值計算.數(shù)值試驗發(fā)現(xiàn),對于不同的介質(zhì)參數(shù),球面波反射系數(shù)表現(xiàn)出了復(fù)雜的頻變規(guī)律.尤其是當(dāng)平面波反射系數(shù)為零時,對應(yīng)的球面波反射系數(shù)是非零的,且球面波反射系數(shù)的相位隨頻率增加在高頻趨近于90°或-90°,即此時球面反射波相對于入射波會有90°的相位旋轉(zhuǎn).對四類AVO模型的測試表明,球面波反射系數(shù)與平面波反射系數(shù)在臨界角附近和低頻時差異很大.

    球面波; 平面波; 反射系數(shù); 頻變

    1 引言

    近十年來,關(guān)于反射球面波的研究再次引起了地球物理學(xué)家的注意,研究內(nèi)容主要聚焦于大偏移距的球面波AVO特征分析和AVO參數(shù)反演.這些研究的前提需要計算球面波反射系數(shù),其計算方法大致可以分為兩類:一類是通過Sommerfeld積分構(gòu)建反射系數(shù)并對其進行數(shù)值計算(Haase, 2004; Ursenbach et al., 2007; 黃饒等, 2009);另一類是利用反射率法或有限差分法進行全波場地震模擬來求取反射系數(shù)(印興耀等, 2006; Alhussain et al., 2008; Zhu and McMechan, 2012; 楊心超, 2012).研究發(fā)現(xiàn)在臨界角附近平面波反射系數(shù)和球面波反射系數(shù)之間存在很大的差異,此時平面波理論是失效的.鑒于此,Alhussain等(2008),Skopintseva等(2011),Zhu和McMechan(2012)利用球面波反射系數(shù)進行AVO反演,與常規(guī)的基于平面波反射系數(shù)的AVO反演相比,反演結(jié)果得到了顯著的改善.

    上述研究主要是利用大入射角度的球面波反射系數(shù)隨入射角變化的特征,未對球面波反射系數(shù)與頻率的關(guān)系做深入分析.根據(jù)經(jīng)典的幾何地震學(xué),平面波反射系數(shù)是與頻率無關(guān)的,但球面波反射系數(shù)是隨頻率變化(頻變)的.忽略球面波反射系數(shù)的頻變特征可能會導(dǎo)致不精確的、甚至錯誤的地下構(gòu)造解釋,尤其是對于淺層、低頻率和反射界面曲率很大的情況.因此,為了更精確地描述實際地震波的反射特征、充分利用地震數(shù)據(jù)含有的有用信息,對球面波反射系數(shù)頻變特征進行研究具有重要的意義.

    目前并沒有論文和著作闡述彈性介質(zhì)中球面波反射系數(shù)的頻變機制.本文對兩層半無限大彈性介質(zhì)中球面波PP反射系數(shù)的頻變特征進行系統(tǒng)地研究,并與平面波反射系數(shù)進行對比.文中首先基于Sommerfeld積分建立球面波PP反射系數(shù)的積分表達式,采用自適應(yīng)的Gauss-Kronrod求積算法(Shampine, 2008)對其進行穩(wěn)定的數(shù)值計算.然后通過數(shù)值測試揭示球面波PP反射系數(shù)復(fù)雜的頻變特征,并對頻變特征進行分類.最后討論該研究的意義.

    2 方法原理

    對于圖1所示的兩均勻半空間彈性介質(zhì)模型,當(dāng)激發(fā)的P波震源位于上層介質(zhì)時(圖1中的S點),根據(jù)Aki和Richards (1980),在接收位置P處的球面反射P波位函數(shù)可表示為

    (1)

    圖1 本文研究的兩層介質(zhì)模型的觀測系統(tǒng).S和S*為點震源及其對應(yīng)的虛震源,P是檢波器位置.實線表示P波,虛線表示SV波Fig.1 Observation system of the two-layer model investigated in this study. S and S* are point source and the corresponding image source, respectively; P is a recording position. Solid lines represent P waves and dashed lines represent SV waves.

    其中,ρ1和ρ2分別為上層和下層介質(zhì)的密度,β1和β2分別為上層和下層介質(zhì)的SV波速度,α2為下層介質(zhì)的P波速度,i2是P波的透射角度,j1和j2分別是SV波的反射角度和透射角度.

    球面反射P波的位移與位函數(shù)有如下關(guān)系

    (3)其中,ur和uz分別為反射P波沿徑向和垂直方向的位移分量.

    本文考慮反射P波沿射線方向的位移分量,可表示為

    uref=ursinθ+uzcosθ

    ×eiω ξ1(h+z)dp,

    (4)

    其中,J1是一階貝塞爾函數(shù),θ為射線入射角,注意這個角度為真實的入射角度,上述的i1為積分變量中的入射角.

    (5)由于虛源入射波的位函數(shù)可以解析地表示為

    (6)則虛源入射波沿射線方向的位移可以容易地得到

    (7)其中,R=(h+z)/cosθ,即地震波由震源出發(fā)到接收點的傳播距離.

    將(4)和(7)代入(5)式,則有

    (8)

    (8)式的積分路徑中含有奇點p=1/α1,這給數(shù)值計算帶來了困難.參考Ursenbach等(2007)的做法,通過如下的變量代換來消除奇點:p=sini1/α1,ξ1=cosi1/α1,于是有(p/ξ1)dp=-d(cosi1)/α1,cosi1∈[1,0]∪[0,i∞].將上述代換代到(8)中,并令x=cosi1,經(jīng)過一系列推導(dǎo)和整理,最后可以得到球面反射P波的位移反射系數(shù)為

    (9)

    根據(jù)Snell定理,可用sini1表示(2)式中的cosi2、cosj1、cosj2和p,再用cosi1表示sini1,由于x=cosi1,最后便可以獲得(9)式中的RPP(x).

    3 數(shù)值結(jié)果及分析

    3.1球面波反射系數(shù)的頻變規(guī)律

    為深入探討球面波PP反射系數(shù)的頻變規(guī)律,下面針對不同的參數(shù)模型繪制其隨頻率的變化曲線并歸納分析.這里只考慮臨界角前較小的入射角度.

    對各種參數(shù)模式的頻變規(guī)律進行了數(shù)值測試,從測試結(jié)果中選取兩個典型的例子進行說明,即α1<α2和α1>α2兩種情況,如表1所示.其他參數(shù)為h=500 m,z=500 m,入射角以0°為例.圖2和圖3分別為當(dāng)α1<α2和α1>α2時的球面波PP反射系數(shù)頻變曲線.為了便于比較不同參數(shù)時球面波反射系數(shù)的頻變規(guī)律,對反射系數(shù)幅值進行規(guī)零化,即用球面波和平面波反射系數(shù)的幅值與對應(yīng)的平面波反射系數(shù)幅值的差值作為縱軸,如圖2a、3a所示,|Rsw|、|Rpw|分別表示球面波反射系數(shù)和平面波反射系數(shù)的幅值,則不同介質(zhì)參數(shù)的平面波反射系數(shù)都由黑色虛線表示;反射系數(shù)的相位不做處理,如圖2b、3b所示.從圖2和圖3中可以看出,對于不同的介質(zhì)參數(shù),球面波PP反射系數(shù)的頻變規(guī)律是很復(fù)雜的,通過詳細地研究分析,可以歸納總結(jié)為如下.

    表1 測試模型參數(shù)Table 1 Parameters of tested models

    當(dāng)α1<α2時(見圖2).

    1) 當(dāng)ρ1=2.50 g·cm-3時,平面波反射系數(shù)為零,而球面波反射系數(shù)是非零的,且其幅值隨頻率增加而單調(diào)遞減,并在高頻趨近于平面波反射系數(shù)的幅值;球面波反射系數(shù)的相位隨著頻率增加而減小,并在高頻趨近于90°,這說明球面反射P波與入射波有90°的相位差.這一現(xiàn)象尚未在目前的文獻中見諸.

    2) 當(dāng)ρ1=2.42 g·cm-3時,平面波反射系數(shù)為較小的正數(shù)(0.0163),球面波反射系數(shù)的幅值在很低頻率處大于平面波反射系數(shù),但在其他的頻率小于平面波反射系數(shù),且隨頻率增加而在高頻趨近于平面波反射系數(shù)的幅值;球面波反射系數(shù)的相位隨頻率增加從一個較大的相位值(約120°)單調(diào)遞減地趨近于平面波反射系數(shù)的相位(0°).當(dāng)ρ1=1.90 g·cm-3時,平面波反射系數(shù)為較大的正數(shù)(0.1364),球面波反射系數(shù)的幅值小于平面波反射系數(shù)的幅值,且隨頻率增大而增大并在高頻趨近平面波反射系數(shù)的幅值;球面波反射系數(shù)的相位隨頻率增加從一個較小的相位值(約40°)逐漸趨近平面波反射系數(shù)的相位(0°).

    圖2 α1<α2時球面波PP反射系數(shù)幅值(a)和相位(b)的頻變規(guī)律α1=2000 m·s-1, α2=2500 m·s-1, ρ2=2 g·cm-3. Fig.2 Frequency-dependent characteristics of spherical-wave PP reflection coefficient when α1<α2. Magnitude (a) and phase (b)

    圖3 α1>α2時球面波PP反射系數(shù)幅值(a)和相位(b)的頻變規(guī)律α1=2500 m·s-1, α2=2000 m·s-1, ρ2=2.5 g·cm-3.Fig.3 Frequency-dependent characteristics of spherical-wave PP reflection coefficient when α1>α2. Magnitude (a) and phase (b)

    3) 當(dāng)ρ1=2.65 g·cm-3時,平面波反射系數(shù)為負數(shù)(-0.0291),球面波反射系數(shù)幅值總是大于對應(yīng)的平面波反射系數(shù)幅值,且隨頻率增加而單調(diào)遞減,并在高頻趨近平面波反射系數(shù)的幅值;球面波反射系數(shù)的相位與平面波反射系數(shù)的相位在低頻處有較小的差異,且隨頻率增加兩者變得幾乎相等(由于此時平面波反射系數(shù)為負數(shù),其相位為180°).

    當(dāng)α1>α2時(見圖3).

    1) 當(dāng)ρ1=2.00 g·cm-3時,平面波反射系數(shù)為零,但球面波反射系數(shù)是非零的,且其幅值隨頻率增加而單調(diào)遞減,并在高頻趨近于平面波反射系數(shù)的幅值;球面波反射系數(shù)的相位隨著頻率增加而減小,并在高頻趨近于-90°,這說明球面反射P波與入射波相比有90°的相位旋轉(zhuǎn).

    2) 當(dāng)ρ1=2.08 g·cm-3時,平面波反射系數(shù)為絕對值較小的負數(shù)(-0.0196),球面波反射系數(shù)的幅值在低頻處大于平面波反射系數(shù),但在其他的頻率小于平面波反射系數(shù),且隨頻率增加在高頻趨近平面波反射系數(shù)的幅值;球面波反射系數(shù)的相位隨頻率增加從約-50°單調(diào)遞減地趨近平面波反射系數(shù)的相位(-180°).當(dāng)ρ1=2.40 g·cm-3時,平面波反射系數(shù)為絕對值較大的負數(shù)(-0.0909),球面波反射系數(shù)的幅值隨頻率增大而單調(diào)遞增并在高頻趨近平面波反射系數(shù)的幅值;球面波反射系數(shù)的相位隨頻率增加而減小,并在高頻趨近平面波反射系數(shù)的相位.

    3) 當(dāng)ρ1=1.80 g·cm-3時,平面波反射系數(shù)為正數(shù)(0.0526),球面波反射系數(shù)幅值總是大于對應(yīng)的平面波反射系數(shù),且隨頻率增加而單調(diào)遞減,并在高頻趨近平面波反射系數(shù)的幅值;球面波反射系數(shù)的相位與平面波反射系數(shù)的相位在低頻處有較小的差異,且隨頻率增加兩者的差異逐漸減小,并在較高頻率趨于相等.

    從上述研究中可以看出,球面波反射系數(shù)是頻變的,平面波反射系數(shù)是非頻變的.需要特別說明的是,平面波反射系數(shù)是非頻變的只是針對均勻平面波而言.實際上,非均勻平面波是頻變的.Sommerfeld積分將球面波分解為均勻平面波和非均勻平面波的疊加,正是由于非均勻平面波的頻變,導(dǎo)致了最后疊加得到的球面反射波是頻變的,因此球面波反射系數(shù)是頻變的.

    3.2單頻的AVO特征

    這部分以AVO分析為例探討球面波反射系數(shù)頻變特征的應(yīng)用.AVO分析廣泛應(yīng)用于地震勘探中,但通?;谄矫娌ɡ碚摚缁赯oeppritz方程或其線性近似.然而平面波理論預(yù)測的反射系數(shù)在臨界角附近是失效的,這促使人們開始關(guān)注更接近真實情況的球面波的AVO響應(yīng).典型地,Haase(2004)、Alhussain等(2008)、李勝軍等(2012)研究了球面PP波反射系數(shù)隨入射角的變化規(guī)律,但其主要是在時間域求取球面波反射系數(shù),即分析的是帶限球面波反射系數(shù)隨入射角變化的特征.為了研究單頻的球面波AVO特征,這里根據(jù)(9)式計算不同頻率時球面諧波反射系數(shù)隨入射角變化的曲線,進而分析頻率對球面波AVO特征的影響.

    測試了四類AVO模型,其中第一、二類模型根據(jù)Haase(2004)所用的模型做了小的修改,第三、四類模型為Castagna等(1998)所用的模型,具體參數(shù)如表2所示.其他的參數(shù)為h=500 m,z=500 m, 入射角范圍為0~80°.根據(jù)上述方法分別計算四類模型參數(shù)在地震頻段內(nèi)三個頻率5 Hz、30 Hz和80 Hz(分別代表低頻、中頻、高頻)時的球面波反射系數(shù)隨入射角變化的曲線,為便于比較,也給出了根據(jù)(2)式計算的平面波反射系數(shù)隨入射角變化的曲線,如圖4—7所示,其中PRC代表平面波反射系數(shù).對比四個圖可以發(fā)現(xiàn),第一、二類AVO的幅值和相位都有較大的頻變;第三、四類AVO的幅值有很小的頻變,相位的頻變也較小.

    表2 第1—4類AVO模型參數(shù)Table 2 Parameters of Class 1—4 AVO models

    圖4 第一類AVO,PP反射系數(shù)幅值(a)和相位(b)隨入射角的變化Fig.4 Class 1 AVO, PP reflection coefficient magnitude (a) and phase (b) variation with incident angles

    圖5 第二類AVO,PP反射系數(shù)幅值(a)和相位(b)隨入射角的變化Fig.5 Class 2 AVO, PP reflection coefficient magnitude (a) and phase (b) variation with incident angles

    圖6 第三類AVO,PP反射系數(shù)幅值(a)和相位(b)隨入射角的變化Fig.6 Class 3 AVO, PP reflection coefficient magnitude (a) and phase (b) variation with incident angles

    圖7 第四類AVO,PP反射系數(shù)幅值(a)和相位(b)隨入射角的變化Fig.7 Class 4 AVO, PP reflection coefficient magnitude (a) and phase (b) variation with incident angles

    對于第三、四類AVO模型,如圖6和圖7所示,在這兩種模型參數(shù)下不會產(chǎn)生首波,因此AVO曲線都比較光滑,沒有出現(xiàn)在第一、二類AVO情況下臨界角后的振蕩現(xiàn)象.中高頻(30 Hz和80 Hz)球面波反射系數(shù)的幅值和平面波反射系數(shù)的幅值差異很小,幾乎是相等的;低頻時(5 Hz)球面波反射系數(shù)幅值與平面波反射系數(shù)幅值之間有較小的差異.另外,不同頻率的球面波反射系數(shù)的相位與平面波反射系數(shù)的相位之間存在一定的差異,且頻率越小這種差異越大,尤其是對于低頻,在某些入射角兩者差異較大.總的來說,頻率對第三、四類AVO模型的球面波AVO響應(yīng)的影響較小,只有當(dāng)頻率很低時這種影響才相對較大.

    綜上所述,對于四類AVO模型,球面波反射系數(shù)與平面波反射系數(shù)在臨界角附近或低頻時差異很大.因此在臨界角附近或低頻時,需要考慮球面波反射的頻變效應(yīng).

    4 討論

    從上述的數(shù)值測試(圖2和圖3)中可以發(fā)現(xiàn),球面波PP反射系數(shù)的頻變規(guī)律是很復(fù)雜的.對于上述的兩種典型情況(α1<α2和α1>α2),球面波PP反射系數(shù)幅值的頻變規(guī)律大體是一致的,其相位的頻變規(guī)律存在較小的差異(比較圖2b和3b).尤其地,當(dāng)平面波PP反射系數(shù)為零時,球面波PP反射系數(shù)卻是非零的,且球面波反射系數(shù)的相位隨頻率增加在高頻趨近于90°或-90°.值得注意的是,這種現(xiàn)象在基于平面波理論的幾何地震學(xué)中是不會發(fā)生的,這豐富了我們對波傳播規(guī)律的認(rèn)識.按照幾何射線理論或平面波理論,當(dāng)平面波入射到界面上某一點后,根據(jù)Snell定理,其反射只對應(yīng)著一條確定的射線路徑.但是根據(jù)Sommerfeld積分,球面波可以分解為平面波的疊加,具體地說,球面波入射到界面上某一點后的反射,可以看作是圍繞著中心射線的一束射線對應(yīng)的平面波反射的加權(quán)疊加(Alhussain et al., 2008),這束射線中的每一條射線都有不同的入射角度,并且都滿足Snell定理.因此盡管圖2和圖3的測試中射線角度是0°,但球面反射波卻是包含0°在內(nèi)的一定范圍入射角的平面反射波的疊加,對于圖2中藍色虛線或圖3中紅色實線對應(yīng)的參數(shù),非零度的平面波反射系數(shù)是非零的,故最后得到的球面波反射系數(shù)也是非零的.

    球面波反射系數(shù)的頻變規(guī)律可以帶給我們很多啟發(fā)性的思考.首先,這一發(fā)現(xiàn)可以豐富我們對地震波反射機制的認(rèn)識,有助于更好地了解實際地下介質(zhì)中地震波的傳播規(guī)律.通過測試可以看出球面波PP反射系數(shù)的頻變規(guī)律在低頻、淺層尤為顯著,因此其對于近地表球面波傳播及反射規(guī)律的研究更有幫助;另外,近些年來隨著低頻勘探技術(shù)的發(fā)展,現(xiàn)已可采集到低至1.5 Hz的地震數(shù)據(jù)(Wang et al., 2013),在如此低的頻率下球面波的頻變效應(yīng)將會很明顯,更精確的勘探需要考慮球面波反射系數(shù)的頻變.另一方面,球面波反射系數(shù)的頻變特征與地下介質(zhì)參數(shù)是緊密相關(guān)的,利用其可能為反演地下介質(zhì)的彈性參數(shù)提供一種新的思路,可以參考王靜波(2014)利用球面聲波反射系數(shù)進行反演的思路,基于(9)式用全局反演方法來反演介質(zhì)的彈性參數(shù).關(guān)于球面波反射系數(shù)頻變特征的具體應(yīng)用需要進一步地研究.

    5 結(jié)論

    本文根據(jù)Sommerfeld積分公式建立了兩層彈性介質(zhì)中球面波PP反射系數(shù)的積分表達式,基于此研究分析了不同參數(shù)情況下球面波PP反射系數(shù)的頻變特征,并對其總結(jié)分類;探討了四類AVO模型的單頻球面波反射系數(shù)隨入射角的變化規(guī)律.概括地,得到了以下幾方面的結(jié)論和認(rèn)識.

    1) 當(dāng)α1<α2、平面波反射系數(shù)為較大的正數(shù)或當(dāng)α1>α2、平面波反射系數(shù)為絕對值較大的負數(shù)時,球面波反射系數(shù)幅值小于平面波反射系數(shù)幅值,并隨頻率的增加而增大,在高頻趨于平面波反射系數(shù);球面波反射系數(shù)相位大于平面波反射系數(shù)相位,且隨頻率增加而減小,在高頻趨于平面波反射系數(shù)相位.

    2) 當(dāng)α1<α2、平面波反射系數(shù)為較小的正數(shù)或當(dāng)α1>α2、平面波反射系數(shù)為絕對值較小的負數(shù)時,球面波反射系數(shù)幅值在很低頻率時大于平面波反射系數(shù)幅值,在其他頻率時小于平面波反射系數(shù)幅值,并隨頻率的增加在高頻趨于平面波反射系數(shù);球面波反射系數(shù)相位大于平面波反射系數(shù)相位,且隨頻率增加而減小,在高頻趨于平面波反射系數(shù)相位.

    3) 無論α1<α2還是α1>α2,當(dāng)平面波反射系數(shù)為零時,球面波反射系數(shù)不為零:球面波反射系數(shù)幅值大于平面波反射系數(shù)幅值,并隨頻率的增加而減小,在高頻趨于平面波反射系數(shù);球面波反射系數(shù)的相位隨頻率增加在高頻趨近于90°或-90°.

    4) 當(dāng)α1<α2、平面波反射系數(shù)為負數(shù)或當(dāng)α1>α2、平面波反射系數(shù)為正數(shù)時,球面波反射系數(shù)幅值大于平面波反射系數(shù)幅值,并隨頻率的增加而減小,在高頻趨于平面波反射系數(shù);球面波反射系數(shù)相位與平面波反射系數(shù)相位在低頻時有較小差異,且隨頻率增加兩者趨于相等.

    5) 對于第一、二類AVO模型,球面波反射系數(shù)與平面波反射系數(shù)在臨界角附近和低頻時差異都很大;對于第三、四類AVO模型,球面波反射系數(shù)與平面波反射系數(shù)僅在頻率很低時才有相對大的差異.

    致謝感謝休斯頓大學(xué)黃龍博士對論文撰寫提供的幫助.

    Aki K, Richards P G. 1980. Quantitative Seismology: Theory and Methods, Vol. 1. San Francisco, CA: W. H. Freeman & Co.

    Alhussain M, Gurevich B, Urosevic M. 2008. Experimental verification of spherical-wave effect on the AVO response and implications for three-term inversion.Geophysics, 73(2): C7-C12.

    Alulaiw B, Gurevich B. 2013. Analytical wavefront curvature correction to plane-wave reflection coefficients for a weak-contrast interface.GeophysicalProspecting, 61(1): 53-62.

    Brekhovskikh L M. 1980. Waves in Layered Media. New York: Academic Press.

    Brekhovskikh L M, Godin O A. 1999. Acoustics of Layered Media II, Point Sources and Bounded Beams. 2nd ed. Berlin Heidelberg: Springer. Castagna J P, Swan H W, Foster D J. 1998. Framework for AVO gradient and intercept interpretation.Geophysics, 63(3): 948-956.

    Du S T, Yu K Y. 1985. Lecture notes on the theory of seismic wave dynamics (in Chinese). Dongying: Geophysical Exploration Department of East China Petroleum Institute.

    Haase A B. 2004. Spherical wave AVO modeling of converted waves in isotropic media. ∥ 74thSEG Annual Meeting. Denver, Colorado: SEG, 263-266.

    Huang R, Chen X H, Li J Y. 2009. Seismogram synthesis method based on spherical wave and AVO analysis.JournalofChinaUniversityofPetroleum(EditionofNaturalScience) (in Chinese), 33(3): 49-52.

    Li S J, Gao J H, Yong X S, et al. 2012. Spherical-wave AVO analysis of reflected P-wave.ChineseJ.Geophys. (in Chinese), 55(10): 3459-3466, doi: 10.6038/j.issn.0001-5733.2012.10.029. Shampine L F. 2008. Vectorized adaptive quadrature in Matlab.J.Comput.Appl.Math., 211(2): 131-140.

    Skopintseva L, Ayzenberg M, LandrM, et al. 2011. Long-offset AVO inversion of PP reflections from plane interfaces using effective reflection coefficients.Geophysics, 76(6): C65-C79.

    Sommerfeld A. 1909. über die Ausbreitung der Wellen in der drahtlosen Telegraphie.AnnalenderPhysik, 333(4): 665-736.

    Ursenbach C P, Haase A B, Downton J E. 2007. Efficient spherical-wave AVO modeling.TheLeadingEdge, 26(12): 1584-1589.

    Wang J B. 2014. Frequency-dependence of spherical-wave reflection and transmission coefficients: Analysis and Application [Ph. D. thesis] (in Chinese). Beijing: China University of Petroleum (Beijing).

    Wang J F, Ma T, Liu Z G, et al. 2013. A new method to extend the range of low-frequency vibroseis sweeps. ∥ 83rdAnnual Meeting. Houston, Texas: SEG, 254-258.

    Weyl H. 1919. Ausbreitung elektromagnetischer Wellen über einem ebenen Leiter.AnnalenderPhysik, 365(21): 481 -500.

    Xu J X, Wu L. 1997. The theoretical study of the elastic wave: on the first term of spherical wave.GeophysicalProspectingforPetroleum(in Chinese), 36(1): 1-12. Yang X C. 2012. Study of characteristics of spherical-wave reflection coefficients [M. S. thesis] (in Chinese). Beijing: China University of Petroleum (Beijing).

    Yin X Y, Zhao J, Zhang F C, et al. 2006. Method of generating prestack synthetic seismograms of spherical wave.JournalofChinaUniversityofPetroleum(EditionofNaturalScience) (in Chinese), 30(1): 26-32.

    Zhu X F, McMechan G A. 2012. Elastic inversion of near- and postcritical reflections using phase variation with angle.Geophysics, 77(4): R149-R159.

    附中文參考文獻

    杜世通, 俞康胤. 1985. 地震波動力學(xué)理論講義. 東營: 華東石油學(xué)院物探教研室.

    黃饒, 陳小宏, 李景葉. 2009. 基于球面波的地震記錄合成方法及AVO分析. 中國石油大學(xué)學(xué)報(自然科學(xué)版), 33(3): 49-52.

    李勝軍, 高建虎, 雍學(xué)善等. 2012. 球面波的反射P波AVO分析. 地球物理學(xué)報, 55(10): 3459-3466, doi: 10.6038/j.issn.0001-5733.2012.10.029.

    王靜波. 2014. 球面波反射和透射系數(shù)頻變特征分析與應(yīng)用[博士論文]. 北京: 中國石油大學(xué)(北京).

    徐基祥, 吳律. 1997. 彈性波理論研究—關(guān)于球面波首項(Ⅰ). 石油物探, 36(1): 1-12.

    楊心超. 2012. 球面波反射系數(shù)特征研究[碩士論文]. 北京: 中國石油大學(xué)(北京).

    印興耀, 趙劍, 張繁昌等. 2006. 球面波疊前地震記錄的合成方法. 中國石油大學(xué)學(xué)報(自然科學(xué)版), 30(1): 26-32.

    (本文編輯胡素芳)

    Study on frequency-dependent characteristics of spherical-wave PP reflection coefficient

    LI Jing-Nan1, 2, WANG Shang-Xu1, 2*, DONG Chun-Hui1, 2, YUAN San-Yi1, 2, WANG Jing-Bo3

    1StateKeyLaboratoryofPetroleumResourceandProspecting,ChinaUniversityofPetroleum(Beijing),Beijing102249,China2CNPCKeyLabofGeophysicalExploration,ChinaUniversityofPetroleum(Beijing),Beijing102249,China3SinopecExplorationCompany,Chengdu610041,China

    Spherical-wave reflection coefficient (SRC) describes the seismic wave reflection more precisely and realistically than the plane-wave reflection coefficient (PRC). This study mainly focuses on the frequency-dependent characteristics of SRC, which are usually not considered. To better understand the reflection mechanism of spherical waves, we systematically investigate the frequency-dependent characteristics of SRC (P-wave only) in a two-layer elastic medium model with a planar interface. We use the classical Sommerfeld integral to construct the SRC and compute it with adaptive Gauss-Kronrod quadrature. Numerical simulations show that for different media with various properties, the frequency-dependence of SRC exhibits very complicated responses. Significant deviations between SRC and PRC are observed at low frequencies. Especially when the PRC is zero, the corresponding SRC is non-zero; with increasing frequency, the phase of the SRC approaches +90° or -90°, meaning a 90° phase rotation after reflection. Tests on Class 1-4 AVO models suggest that large discrepancy between SRC and PRC exists in the vicinity of the critical angle and at low frequencies.

    Spherical wave; Plane wave; Reflection coefficient; Frequency dependent

    10.6038/cjg20161025.

    國家重點基礎(chǔ)研究發(fā)展計劃(973計劃)項目(2013CB228600)和國家自然科學(xué)基金(41304108)聯(lián)合資助.

    李京南,男,1988年生,博士研究生,主要從事地震波傳播、地震資料處理方面的研究.E-mail:li.jnan@163.com

    王尚旭,教授,博士生導(dǎo)師,主要從事地震波傳播理論、巖石物理、地震波正演模擬及反演研究.E-mail:wangsx@cup.edu.cn

    10.6038/cjg20161025

    P631

    2016-03-17,2016-07-23收修定稿

    李京南, 王尚旭, 董春暉等. 2016. 球面波PP反射系數(shù)的頻變特征研究. 地球物理學(xué)報,59(10):3810-3819,

    Li J N, Wang S X, Dong C H, et al. 2016. Study on frequency-dependent characteristics of spherical-wave PP reflection coefficient.ChineseJ.Geophys. (in Chinese),59(10):3810-3819,doi:10.6038/cjg20161025.

    猜你喜歡
    平面波反射系數(shù)入射角
    一般三棱鏡偏向角與入射角的關(guān)系
    Landau-Lifshitz方程平面波解的全局光滑性
    5G OTA測量寬帶平面波模擬器的高效優(yōu)化方法與應(yīng)用
    多道隨機稀疏反射系數(shù)反演
    石油物探(2020年6期)2020-11-25 02:38:46
    預(yù)制圓柱形鎢破片斜穿甲鋼靶的破孔能力分析*
    用經(jīng)典定理證明各向異性巖石界面異常入射角的存在
    復(fù)合函數(shù)漸變傳輸線研究
    基于GPU并行運算的超聲平面波成像仿真
    電子制作(2016年11期)2016-11-07 08:43:45
    沙質(zhì)沉積物反射系數(shù)的寬帶測量方法
    基于反射系數(shù)的波導(dǎo)結(jié)構(gòu)不連續(xù)位置識別
    中国美女看黄片| 夜夜躁狠狠躁天天躁| 亚洲av五月六月丁香网| 日本与韩国留学比较| 亚洲av美国av| www.www免费av| 亚洲精品456在线播放app | 亚洲av不卡在线观看| 中文字幕久久专区| 999久久久精品免费观看国产| 国产精品香港三级国产av潘金莲| а√天堂www在线а√下载| 在线免费观看不下载黄p国产 | svipshipincom国产片| 一个人看视频在线观看www免费 | 制服人妻中文乱码| 国产视频一区二区在线看| 国产成人系列免费观看| 叶爱在线成人免费视频播放| 99久久精品国产亚洲精品| 可以在线观看毛片的网站| 久久国产精品影院| 免费电影在线观看免费观看| 男人的好看免费观看在线视频| 亚洲人成网站在线播| 欧美+日韩+精品| 国产精品久久久久久亚洲av鲁大| 欧美性感艳星| 亚洲精品456在线播放app | 天堂av国产一区二区熟女人妻| 国产精品久久久久久久电影 | 成人欧美大片| 中文字幕av成人在线电影| 长腿黑丝高跟| 精品欧美国产一区二区三| 国产又黄又爽又无遮挡在线| 99视频精品全部免费 在线| av欧美777| 亚洲欧美一区二区三区黑人| 一卡2卡三卡四卡精品乱码亚洲| 国产精品亚洲美女久久久| 神马国产精品三级电影在线观看| 99国产综合亚洲精品| 狂野欧美白嫩少妇大欣赏| 一级毛片女人18水好多| 夜夜看夜夜爽夜夜摸| 日本 av在线| 欧美zozozo另类| 亚洲第一欧美日韩一区二区三区| 亚洲不卡免费看| 欧美日韩瑟瑟在线播放| 亚洲美女视频黄频| 一本久久中文字幕| 国产成人av激情在线播放| 午夜两性在线视频| 亚洲熟妇中文字幕五十中出| 女警被强在线播放| 国产精品亚洲美女久久久| 国产三级在线视频| 岛国在线免费视频观看| 黄色丝袜av网址大全| h日本视频在线播放| 色噜噜av男人的天堂激情| 91在线精品国自产拍蜜月 | 国产精品久久久久久精品电影| 99久久无色码亚洲精品果冻| 制服人妻中文乱码| 国产乱人视频| a级毛片a级免费在线| 熟女人妻精品中文字幕| 丁香六月欧美| 嫩草影视91久久| a级毛片a级免费在线| 99久久精品国产亚洲精品| 搡老熟女国产l中国老女人| 国产野战对白在线观看| а√天堂www在线а√下载| 午夜福利成人在线免费观看| 欧美不卡视频在线免费观看| 久久中文看片网| 国产成人影院久久av| 亚洲天堂国产精品一区在线| 动漫黄色视频在线观看| 久久久久久国产a免费观看| 国产中年淑女户外野战色| 一区二区三区国产精品乱码| 国产成人系列免费观看| a在线观看视频网站| 制服人妻中文乱码| 伊人久久大香线蕉亚洲五| 欧美日韩中文字幕国产精品一区二区三区| 国产亚洲欧美98| 欧美乱色亚洲激情| 久久99热这里只有精品18| 午夜福利在线观看吧| 精品人妻一区二区三区麻豆 | а√天堂www在线а√下载| 亚洲人成网站在线播| 在线观看66精品国产| 日韩欧美一区二区三区在线观看| 1024手机看黄色片| 中文字幕av在线有码专区| 久久久久亚洲av毛片大全| 乱人视频在线观看| 久久久久精品国产欧美久久久| 在线观看免费视频日本深夜| 波多野结衣高清作品| 国产淫片久久久久久久久 | 亚洲最大成人中文| 免费电影在线观看免费观看| 国产日本99.免费观看| 99国产精品一区二区三区| 在线观看一区二区三区| 午夜精品在线福利| 欧美乱码精品一区二区三区| 亚洲片人在线观看| 国产欧美日韩一区二区精品| 国产伦在线观看视频一区| 人人妻人人澡欧美一区二区| 亚洲五月天丁香| 精品人妻一区二区三区麻豆 | 久久久久国内视频| 精品久久久久久久久久免费视频| 久久久久亚洲av毛片大全| 精品免费久久久久久久清纯| 男女那种视频在线观看| 婷婷精品国产亚洲av在线| 制服丝袜大香蕉在线| 午夜福利成人在线免费观看| 国产高潮美女av| 欧美性猛交黑人性爽| 夜夜夜夜夜久久久久| www日本在线高清视频| 国产精品av视频在线免费观看| 性色avwww在线观看| 淫妇啪啪啪对白视频| 在线观看av片永久免费下载| 日本 欧美在线| 真实男女啪啪啪动态图| 日本五十路高清| 日本a在线网址| 欧美中文综合在线视频| 老熟妇仑乱视频hdxx| 中文字幕熟女人妻在线| 在线播放国产精品三级| 伊人久久大香线蕉亚洲五| 18禁美女被吸乳视频| 在线播放国产精品三级| av中文乱码字幕在线| 免费人成视频x8x8入口观看| 国内精品久久久久精免费| 好男人电影高清在线观看| 欧美一区二区亚洲| 观看免费一级毛片| 日本五十路高清| 天天躁日日操中文字幕| 国产激情偷乱视频一区二区| 国产一区二区在线观看日韩 | 白带黄色成豆腐渣| 有码 亚洲区| 亚洲一区高清亚洲精品| 精品国产超薄肉色丝袜足j| 国产精品三级大全| 色尼玛亚洲综合影院| 日本黄大片高清| 亚洲国产欧美人成| 久久伊人香网站| 不卡一级毛片| 99久久无色码亚洲精品果冻| 欧美高清成人免费视频www| 美女免费视频网站| 在线观看午夜福利视频| 搡老熟女国产l中国老女人| 99久久精品热视频| 非洲黑人性xxxx精品又粗又长| 国产精品乱码一区二三区的特点| 亚洲国产欧洲综合997久久,| 精品99又大又爽又粗少妇毛片 | 白带黄色成豆腐渣| 久久精品影院6| 久久精品国产99精品国产亚洲性色| 国产精品99久久99久久久不卡| 欧美黑人欧美精品刺激| 老司机午夜福利在线观看视频| 久久国产精品影院| 观看美女的网站| av黄色大香蕉| 日韩欧美精品v在线| 中文字幕精品亚洲无线码一区| 少妇高潮的动态图| 嫩草影院精品99| 欧美极品一区二区三区四区| 亚洲专区中文字幕在线| 十八禁人妻一区二区| 少妇的逼水好多| 一a级毛片在线观看| 亚洲成av人片免费观看| 国产精品一区二区三区四区免费观看 | 亚洲精品影视一区二区三区av| 高潮久久久久久久久久久不卡| 国产男靠女视频免费网站| 中国美女看黄片| 色综合婷婷激情| 国产精品久久久人人做人人爽| 无人区码免费观看不卡| 欧美最新免费一区二区三区 | 性色av乱码一区二区三区2| 高清在线国产一区| 久久精品人妻少妇| 欧美成人a在线观看| 可以在线观看毛片的网站| 男人舔女人下体高潮全视频| 18禁黄网站禁片午夜丰满| h日本视频在线播放| 亚洲av日韩精品久久久久久密| 一卡2卡三卡四卡精品乱码亚洲| 每晚都被弄得嗷嗷叫到高潮| 99国产精品一区二区蜜桃av| 欧美高清成人免费视频www| 色综合亚洲欧美另类图片| 老司机午夜十八禁免费视频| 精品国产超薄肉色丝袜足j| 噜噜噜噜噜久久久久久91| 午夜a级毛片| 久久性视频一级片| 国产精品国产高清国产av| 长腿黑丝高跟| 国产综合懂色| 久久精品国产自在天天线| 亚洲av免费高清在线观看| 国产熟女xx| 国产欧美日韩一区二区三| bbb黄色大片| 国产精品乱码一区二三区的特点| av欧美777| 午夜福利视频1000在线观看| 丰满乱子伦码专区| 国内毛片毛片毛片毛片毛片| 狂野欧美激情性xxxx| 一区二区三区国产精品乱码| 日本与韩国留学比较| 日韩欧美精品免费久久 | 国内久久婷婷六月综合欲色啪| 免费无遮挡裸体视频| 中文字幕人妻丝袜一区二区| 亚洲成人久久爱视频| 韩国av一区二区三区四区| av中文乱码字幕在线| 国产三级黄色录像| 午夜激情福利司机影院| 欧美绝顶高潮抽搐喷水| 国内揄拍国产精品人妻在线| 男人的好看免费观看在线视频| 波多野结衣巨乳人妻| 在线观看一区二区三区| 99在线视频只有这里精品首页| 国产伦在线观看视频一区| 69人妻影院| 国产精品久久久久久久久免 | 久久久久免费精品人妻一区二区| 久久精品夜夜夜夜夜久久蜜豆| av视频在线观看入口| 亚洲第一电影网av| 亚洲 欧美 日韩 在线 免费| 久久香蕉国产精品| 又爽又黄无遮挡网站| 变态另类成人亚洲欧美熟女| 国产一区二区激情短视频| 久久草成人影院| 婷婷精品国产亚洲av| 免费电影在线观看免费观看| bbb黄色大片| 老司机在亚洲福利影院| 国产国拍精品亚洲av在线观看 | 人妻丰满熟妇av一区二区三区| 日本免费一区二区三区高清不卡| 制服丝袜大香蕉在线| 国产精品 欧美亚洲| 日本黄色视频三级网站网址| 亚洲无线在线观看| 最新在线观看一区二区三区| 香蕉久久夜色| 在线观看一区二区三区| 亚洲成av人片免费观看| 又黄又粗又硬又大视频| 九九久久精品国产亚洲av麻豆| 成人性生交大片免费视频hd| 久久精品国产综合久久久| 亚洲真实伦在线观看| 俄罗斯特黄特色一大片| 亚洲人成伊人成综合网2020| 国产精品99久久久久久久久| 欧美bdsm另类| 午夜福利在线观看吧| 99热这里只有是精品50| 久久欧美精品欧美久久欧美| 深爱激情五月婷婷| 校园春色视频在线观看| 日韩国内少妇激情av| 男女下面进入的视频免费午夜| 最近最新免费中文字幕在线| 夜夜躁狠狠躁天天躁| 中文字幕av成人在线电影| 香蕉丝袜av| 国产精品一区二区免费欧美| 日韩免费av在线播放| 欧美zozozo另类| 淫妇啪啪啪对白视频| 国产av在哪里看| 夜夜爽天天搞| 日韩欧美国产在线观看| 国内精品美女久久久久久| 噜噜噜噜噜久久久久久91| 久久人妻av系列| 少妇人妻一区二区三区视频| 日本a在线网址| 在线播放国产精品三级| 一级黄色大片毛片| 丰满乱子伦码专区| 一卡2卡三卡四卡精品乱码亚洲| 亚洲中文字幕一区二区三区有码在线看| 国产三级黄色录像| 亚洲人成网站在线播| 免费在线观看影片大全网站| АⅤ资源中文在线天堂| 99热6这里只有精品| 国产精品一区二区三区四区免费观看 | av福利片在线观看| www.999成人在线观看| 国产麻豆成人av免费视频| 在线播放无遮挡| 亚洲国产精品999在线| 日韩av在线大香蕉| 深夜精品福利| 日日干狠狠操夜夜爽| 精品一区二区三区av网在线观看| 精品电影一区二区在线| a级毛片a级免费在线| 特大巨黑吊av在线直播| 亚洲最大成人中文| 亚洲美女黄片视频| 日本 欧美在线| 少妇的逼水好多| 欧美最新免费一区二区三区 | 熟女人妻精品中文字幕| 国内精品久久久久久久电影| 99热6这里只有精品| 狂野欧美白嫩少妇大欣赏| 狂野欧美激情性xxxx| 日本撒尿小便嘘嘘汇集6| 国产色爽女视频免费观看| 小说图片视频综合网站| 色在线成人网| 18禁裸乳无遮挡免费网站照片| 美女被艹到高潮喷水动态| 在线观看一区二区三区| x7x7x7水蜜桃| 亚洲精品一卡2卡三卡4卡5卡| 国产精品,欧美在线| 中文字幕精品亚洲无线码一区| 婷婷精品国产亚洲av| 亚洲av日韩精品久久久久久密| 亚洲人成网站在线播| 免费在线观看影片大全网站| 精品免费久久久久久久清纯| 尤物成人国产欧美一区二区三区| 欧美一区二区精品小视频在线| 91在线精品国自产拍蜜月 | 国产免费av片在线观看野外av| bbb黄色大片| 网址你懂的国产日韩在线| 老汉色av国产亚洲站长工具| 女人被狂操c到高潮| 色老头精品视频在线观看| 亚洲熟妇中文字幕五十中出| 狂野欧美激情性xxxx| 三级国产精品欧美在线观看| 国产精品精品国产色婷婷| 蜜桃久久精品国产亚洲av| 日本黄色视频三级网站网址| 天天添夜夜摸| 国产在视频线在精品| 男人舔奶头视频| 男女视频在线观看网站免费| 看片在线看免费视频| 亚洲国产精品sss在线观看| 久久九九热精品免费| 搞女人的毛片| 最好的美女福利视频网| 亚洲国产精品成人综合色| 免费人成在线观看视频色| 午夜a级毛片| 欧美日韩福利视频一区二区| 亚洲欧美日韩东京热| 国产精品三级大全| 国产精品久久电影中文字幕| 人妻夜夜爽99麻豆av| 老鸭窝网址在线观看| 麻豆一二三区av精品| 欧美av亚洲av综合av国产av| 精品欧美国产一区二区三| 日韩大尺度精品在线看网址| 亚洲av二区三区四区| 日本精品一区二区三区蜜桃| 久久精品国产99精品国产亚洲性色| 12—13女人毛片做爰片一| 免费人成视频x8x8入口观看| 午夜日韩欧美国产| 91麻豆精品激情在线观看国产| www国产在线视频色| 欧美一区二区精品小视频在线| 中亚洲国语对白在线视频| 人人妻人人看人人澡| 欧美午夜高清在线| 久久久久性生活片| 色哟哟哟哟哟哟| 亚洲七黄色美女视频| 国内精品久久久久久久电影| 怎么达到女性高潮| 欧美最新免费一区二区三区 | 99精品久久久久人妻精品| 国产av麻豆久久久久久久| 久久这里只有精品中国| 欧美色欧美亚洲另类二区| 动漫黄色视频在线观看| 久久6这里有精品| 国产精品野战在线观看| 久久这里只有精品中国| 18禁黄网站禁片午夜丰满| 老熟妇乱子伦视频在线观看| 久久久久久久久久黄片| x7x7x7水蜜桃| 日本熟妇午夜| 草草在线视频免费看| 欧美日韩福利视频一区二区| 18+在线观看网站| 最好的美女福利视频网| 夜夜爽天天搞| 免费在线观看日本一区| 999久久久精品免费观看国产| 成人三级黄色视频| 免费观看的影片在线观看| 蜜桃久久精品国产亚洲av| 婷婷亚洲欧美| 国产一区二区三区在线臀色熟女| 黄色女人牲交| 真人一进一出gif抽搐免费| 一区二区三区激情视频| 叶爱在线成人免费视频播放| 亚洲在线观看片| 高清日韩中文字幕在线| 又爽又黄无遮挡网站| 日韩人妻高清精品专区| 亚洲五月婷婷丁香| 国产乱人伦免费视频| or卡值多少钱| 国产激情欧美一区二区| 久久久久久久久久黄片| 欧美国产日韩亚洲一区| 美女大奶头视频| 51国产日韩欧美| 欧美一区二区精品小视频在线| 禁无遮挡网站| 老熟妇乱子伦视频在线观看| 久久久精品欧美日韩精品| 午夜精品一区二区三区免费看| 国产高清激情床上av| 国产v大片淫在线免费观看| 久久久久亚洲av毛片大全| 久99久视频精品免费| 国产极品精品免费视频能看的| 欧美xxxx黑人xx丫x性爽| 国产探花极品一区二区| 99国产综合亚洲精品| 欧美高清成人免费视频www| 亚洲性夜色夜夜综合| 在线播放无遮挡| 国产精品日韩av在线免费观看| 亚洲av电影不卡..在线观看| 亚洲国产精品999在线| 日韩有码中文字幕| 成人特级黄色片久久久久久久| 欧美成人a在线观看| 又黄又爽又免费观看的视频| 久久久国产精品麻豆| 国产精品电影一区二区三区| 一进一出好大好爽视频| 亚洲avbb在线观看| 亚洲片人在线观看| 亚洲国产中文字幕在线视频| 村上凉子中文字幕在线| 国产精品美女特级片免费视频播放器| 国产精品国产高清国产av| 亚洲熟妇熟女久久| 国产久久久一区二区三区| 国产三级在线视频| 国产美女午夜福利| 人妻夜夜爽99麻豆av| 国产精品嫩草影院av在线观看 | 亚洲无线在线观看| 脱女人内裤的视频| 国产精品,欧美在线| 久久人妻av系列| 宅男免费午夜| av黄色大香蕉| 97超视频在线观看视频| 中文在线观看免费www的网站| 精品国产美女av久久久久小说| 人妻久久中文字幕网| 日韩 欧美 亚洲 中文字幕| 免费看日本二区| 一级a爱片免费观看的视频| 国产成人欧美在线观看| 村上凉子中文字幕在线| 国产97色在线日韩免费| 精品人妻1区二区| 国内精品久久久久久久电影| 成人国产一区最新在线观看| 国产欧美日韩一区二区精品| 人妻夜夜爽99麻豆av| 欧美乱妇无乱码| 欧美区成人在线视频| 亚洲国产精品成人综合色| 欧美日本视频| 国产精品自产拍在线观看55亚洲| 午夜免费男女啪啪视频观看 | 18禁美女被吸乳视频| 国产精品一区二区免费欧美| 精华霜和精华液先用哪个| 无人区码免费观看不卡| 我的老师免费观看完整版| 国产精品嫩草影院av在线观看 | 婷婷六月久久综合丁香| 日韩中文字幕欧美一区二区| 国模一区二区三区四区视频| 国产老妇女一区| 国产成人影院久久av| 亚洲成人免费电影在线观看| 色精品久久人妻99蜜桃| 亚洲av二区三区四区| 在线观看午夜福利视频| 九九在线视频观看精品| 亚洲最大成人手机在线| 在线天堂最新版资源| 欧美乱妇无乱码| 久久久久久九九精品二区国产| 国产亚洲欧美98| 亚洲成人精品中文字幕电影| 欧美黄色淫秽网站| 在线观看av片永久免费下载| 少妇人妻精品综合一区二区 | 亚洲五月婷婷丁香| 精品一区二区三区视频在线观看免费| 男人和女人高潮做爰伦理| 我的老师免费观看完整版| 最近在线观看免费完整版| 伊人久久大香线蕉亚洲五| 国产97色在线日韩免费| 日韩 欧美 亚洲 中文字幕| 日韩欧美在线二视频| 久久99热这里只有精品18| 国产aⅴ精品一区二区三区波| 99在线人妻在线中文字幕| 真人做人爱边吃奶动态| 最近视频中文字幕2019在线8| 一个人看视频在线观看www免费 | 免费av观看视频| 日本a在线网址| 国产毛片a区久久久久| 国产精品久久视频播放| 国产极品精品免费视频能看的| 嫁个100分男人电影在线观看| 不卡一级毛片| 精品无人区乱码1区二区| 亚洲国产精品sss在线观看| 欧美极品一区二区三区四区| 久久婷婷人人爽人人干人人爱| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美绝顶高潮抽搐喷水| 欧美zozozo另类| 99视频精品全部免费 在线| 免费看十八禁软件| 三级男女做爰猛烈吃奶摸视频| 色精品久久人妻99蜜桃| 精品一区二区三区视频在线观看免费| 真人一进一出gif抽搐免费| 在线观看av片永久免费下载| 色视频www国产| 亚洲精品成人久久久久久| 亚洲久久久久久中文字幕| 久久婷婷人人爽人人干人人爱| 国产综合懂色| 亚洲国产欧洲综合997久久,| 麻豆一二三区av精品| 欧美一区二区亚洲| 婷婷亚洲欧美| 国产成人影院久久av| av在线蜜桃| 最近在线观看免费完整版| 午夜福利高清视频| 日本黄色片子视频| 国产精品av视频在线免费观看| 中文在线观看免费www的网站| 欧美黑人欧美精品刺激| 国内精品久久久久精免费| 18禁在线播放成人免费| 亚洲色图av天堂| 男人的好看免费观看在线视频| 12—13女人毛片做爰片一| 亚洲国产色片| 亚洲专区中文字幕在线| 女警被强在线播放| 国产精品久久久久久亚洲av鲁大| 国产单亲对白刺激| 香蕉av资源在线|