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

    正交各向異性介質(zhì)反射系數(shù)精確解

    2020-10-17 07:44:08張雪瑩孫鵬遠(yuǎn)馬學(xué)軍李夢(mèng)琦
    石油地球物理勘探 2020年5期
    關(guān)鍵詞:反射系數(shù)偏振矢量

    張雪瑩 孫鵬遠(yuǎn) 馬學(xué)軍 蘆 俊 李夢(mèng)琦

    (①中國地質(zhì)大學(xué)(北京)能源學(xué)院,北京100083;②東方地球物理公司物探技術(shù)研究中心,河北涿州072751;③中國石化西北油田分公司勘探開發(fā)研究院,烏魯木齊830011;④中國地質(zhì)大學(xué)(北京)地球物理與信息技術(shù)學(xué)院,北京100083)

    0 引言

    中國陸相沉積巖大多具有薄互層特征,當(dāng)受到構(gòu)造運(yùn)動(dòng)的影響時(shí),會(huì)發(fā)育垂向或近似垂向的高角度裂縫,呈現(xiàn)出典型的正交各向異性特征。作為一種重要的含油氣儲(chǔ)層,裂縫型儲(chǔ)層中發(fā)育的裂縫既可以作為油氣儲(chǔ)集空間,又可以作為滲流通道[1]。研究正交各向異性的AVO 響應(yīng)對(duì)裂縫型儲(chǔ)層的精細(xì)刻畫有重要意義。

    裂縫儲(chǔ)層的各向異性是指含裂縫地層的彈性性質(zhì)在不同方向上具有差異,且與誘導(dǎo)裂縫的構(gòu)造應(yīng)力有關(guān)[2]。國外學(xué)者在各向異性領(lǐng)域的研究起步較早,Christoffel[3]提出用彈性張量描述介質(zhì)的各向異性,因此各向異性波動(dòng)方程稱為Christoffel方程。Crampin[4]首次發(fā)現(xiàn)各向異性介質(zhì)中的地震橫波分裂現(xiàn)象,并分析了地球上常見的各向異性介質(zhì)的種類與特征。在此基礎(chǔ)上,Hudson[5-6]推導(dǎo)了由于裂縫存在而產(chǎn)生的一階和二階擾動(dòng)量,給出了含定向裂縫介質(zhì)的彈性模量的計(jì)算方法,但只適用于高頻情況。Schoenberg等[7-8]在位移不連續(xù)和旋轉(zhuǎn)不變性的假設(shè)下,忽略裂縫的形狀與結(jié)構(gòu),提出了線性滑移理論,并對(duì)多組裂縫尺度進(jìn)行了數(shù)值模擬。Thomsen[9]提出了弱各向異性理論,給出了表征弱各向異性VTI介質(zhì)的參數(shù),并據(jù)此推導(dǎo)了相速度和群速度的近似表達(dá)式。當(dāng)介質(zhì)中存在與地層走向垂直的裂隙時(shí),使用正交各向異性模型描述更合理。Tsvankin[10]基于Thomsen的TI介質(zhì)參數(shù),分別在三個(gè)對(duì)稱面內(nèi)求解Christoffel方程,給出了正交各向異性介質(zhì)的參數(shù)化方法,并推導(dǎo)了弱正交各向異性介質(zhì)中qP波相速度近似表達(dá)式?;赥svankin提出的正交各向異性參數(shù),Bakulin等[11-13]對(duì)兩種正交各向異性模型進(jìn)行了研究,分別在VTI背景下嵌入一組垂直裂縫和兩組相互垂直的裂縫,運(yùn)用Schoenberg的柔度建模理論構(gòu)建剛度系數(shù)矩陣。Rokhlin[14]基于Christoffel方程和邊界條件提出了一種計(jì)算一般各向異性介質(zhì)反射系數(shù)和透射系數(shù)的方法;Rüger[15-17]推導(dǎo)了弱各向異性介質(zhì)中反射系數(shù)近似公式,并運(yùn)用縱波方位AVO 檢測(cè)裂縫。

    國內(nèi)學(xué)者在各向異性領(lǐng)域的研究起步較晚,但是在正交各向異性方面的研究也有一些進(jìn)展,傅旦丹等[18]研究了偽譜法在正交各向異性介質(zhì)中的應(yīng)用,分析了偽譜法的穩(wěn)定性,并將Backus方法與Hudson理論結(jié)合,推出了周期性薄互層(PTL)與廣泛擴(kuò)容各向異性(EDA)組合的等效正交各向異性彈性常數(shù)的計(jì)算公式。張文生等[19]利用有限差分方法模擬正交各向異性介質(zhì)中地震波場(chǎng);魏建新[20]應(yīng)用物理模擬方法研究了正交各向異性介質(zhì)的波場(chǎng)特征;李娜[21]基于Tsvankin的各向異性參數(shù)研究了正交各向異性介質(zhì)的相速度和群速度。許茜茹等[22]通過耦合S波射線理論和基于迭代的各向異性相速度與偏振矢量的高階近似解,得到了適用于正交各向異性介質(zhì)以qP波入射產(chǎn)生的二階反射、透射系數(shù)的計(jì)算公式。張繁昌等[23]針對(duì)兩組裂縫垂直的正交各向異性介質(zhì),基于Schoenberg的柔度建模理論構(gòu)建剛度矩陣,推導(dǎo)了反射系數(shù)公式,并運(yùn)用Fourier級(jí)數(shù)預(yù)測(cè)了裂縫方位。唐杰等[24]以兩組垂直裂縫的正交各向異性介質(zhì)為研究對(duì)象,研究了干燥和完全水飽和的地震波響應(yīng)特征。郭愷等[25]研究了正交各向異性介質(zhì)的多參數(shù)建模方法;單俊臻等[26]推導(dǎo)了一種新的PP 波HTI介質(zhì)反射系數(shù)一階擾動(dòng)公式,突破了傳統(tǒng)反射系數(shù)近似公式僅適用于小炮間距的局限。秦海旭等[27]將兩組裂縫垂直斜交的各向異性地層等效為單斜介質(zhì),推導(dǎo)了qP波方程并進(jìn)行了正演模擬。

    在上述研究中,基于波動(dòng)方程的正演模擬方法得到的正交各向異性的地震波場(chǎng)非常復(fù)雜,多種模式的波相互混雜,再加上難以完全消除的數(shù)值頻散與邊界效應(yīng),無法用于定量分析正交各向異性介質(zhì)的AVO 響應(yīng)特征。另外,弱正交各向異性假設(shè)下的反射系數(shù)近似公式的精度沒有得到有效驗(yàn)證,無法適用于裂縫發(fā)育程度高或者VTI背景各向異性程度較強(qiáng)的地層。

    針對(duì)上述問題,本文針對(duì)VTI背景介質(zhì)中發(fā)育的一組直立裂縫誘導(dǎo)的正交各向異性,采用Tsvankin提出的各向異性參數(shù)構(gòu)建剛度系數(shù)矩陣,再根據(jù)Christoffel方程求解出慢度矢量和偏振矢量,并結(jié)合彈性波傳播理論的邊界位移連續(xù)條件和應(yīng)力連續(xù)條件,導(dǎo)出了計(jì)算精確反射系數(shù)和透射系數(shù)的方法。使用十二個(gè)各向異性模型分析背景VTI介質(zhì)各向異性強(qiáng)度和裂縫發(fā)育程度對(duì)正交各向異性介質(zhì)反射系數(shù)的影響。為了驗(yàn)證本文計(jì)算結(jié)果的正確性,與Rüger的反射系數(shù)近似公式的計(jì)算結(jié)果進(jìn)行了對(duì)比。

    1 正交各向異性介質(zhì)彈性剛度矩陣

    本文的正交各向異性介質(zhì)是指在VTI背景介質(zhì)中發(fā)育有一組垂直裂縫,如圖1的下層所示。

    圖1中x1軸垂直于裂縫面,對(duì)稱面x2Ox3平行于裂縫面,x1Ox3垂直于裂縫面,x1Ox2平行于水平界面,θ 為入射波傳播方向與x3軸的夾角(入射角),φ 為入射平面與x1軸的夾角(方位角)。文中用上標(biāo)“U、L”分別代表上層介質(zhì)和下層介質(zhì),下標(biāo)“1、2、3”對(duì)應(yīng)軸x1、x2、x3軸向。相速度傳播方向?yàn)椴ㄇ胺ㄏ?,慢度為相速度的倒?shù)。

    正交各向異性介質(zhì)的彈性剛度系數(shù)矩陣為

    Tsvankin[9]定義的正交各向異性介質(zhì)中三個(gè)不同的正交對(duì)稱面上的各向異性參數(shù)為

    式中:VP0為沿著垂直方向(x3方向)傳播的P波速度;VS0為沿著垂直方向傳播且沿著x2方向偏振的S 波 速 度;ε(1,2)、γ(1,2)、δ(1,2,3)表 示 不 同 平 面 內(nèi) 的Thomsen型各向異性參數(shù),其中上標(biāo)“1、2、3”分別代表x2Ox3、x1Ox3和x1Ox2平面[11]。

    可根據(jù)以上Tsvankin的各向異性參數(shù)構(gòu)建彈性剛度系數(shù)矩陣,即

    2 正交各向異性介質(zhì)彈性波反射/透射系數(shù)精確解

    Christoffel方程是由波動(dòng)方程導(dǎo)出,用于研究地震波的慢度和偏振矢量等傳播特征[28]?;贑hristoffel方程及其變形,結(jié)合Snell定律求解出反射波和透射波的慢度矢量和偏振矢量,并根據(jù)邊界的位移連續(xù)條件和應(yīng)力連續(xù)條件,求出反射系數(shù)和透射系數(shù)精確解。

    2.1 入射P波的慢度矢量計(jì)算

    一般均勻完全彈性各向異性介質(zhì)中地震波時(shí)間域Christoffel方程[3]為

    式中:E 為單位矩陣;ρ為介質(zhì)的密度;P 為波的偏振矢量。

    根據(jù)廣義虎克定律,可得

    式中s1、s2、s3為慢度矢量的三個(gè)分量。

    當(dāng)界面上、下介質(zhì)為正交各向異性或是更高對(duì)稱性的介質(zhì)時(shí),假設(shè)入射波以圖1方式入射,則入射波波前的法向單位矢量為

    式(4)兩邊同乘V2/ρ可得[29]

    其展開形式為

    此時(shí)求解式(8)就等價(jià)于相速度V(特征值)和偏振矢量P(特征向量)的求解問題。為滿足非零解條件,Christoffel矩陣的行列式值為零,即變成求解V2的一元三次方程

    解式(11)可得入射波的相速度VP、VS1、VS2(VP>VS1>VS2),則入射P波的慢度為

    2.2 反射波和透射波的慢度矢量計(jì)算

    若用si表示不同波的慢度矢量(i=0、1、…、6,分別代表入射P 波、反射P、S1、S2波、透射P、S1、S2波),則

    根據(jù)Snell定律,對(duì)于所有的入射、反射和透射波,與界面平行的兩個(gè)慢度矢量的分量是相等的,即

    式(4)的未知解是慢度矢量和偏振矢量,對(duì)于反射波和透射波慢度矢量的分量si3,為滿足非零解條件,則需求解方程

    把上層介質(zhì)和下層介質(zhì)的參數(shù)以及已知的si1和si2代入式(16),可以解出上、下介質(zhì)中的六個(gè)si3。

    2.3 反射波和透射波的偏振矢量計(jì)算

    若用Pi表示偏振矢量(i=0、1、…、6,分別代表入射P波、反射P、S1、S2波、透射P、S1、S2波),即

    式(16)已求得慢度矢量,再結(jié)合偏振矢量的歸一化,對(duì)于上層介質(zhì)有

    對(duì)于下層介質(zhì)有

    運(yùn)用Zhou等[30]的方法求解式(18)和式(19),可得到反射波和透射波的偏振矢量。

    對(duì)于求得的7個(gè)偏振矢量,需要通過其與相速度的關(guān)系進(jìn)行偏振矢量符號(hào)的校正。根據(jù)Chen[31]的校正方法,對(duì)于入射、反射和透射的P 波,若偏振矢量P 與慢度矢量s 的點(diǎn)乘大于0,則偏振矢量的符號(hào)為正;對(duì)于反射和透射的S波,入射面的法向量N 為慢度矢量s 與入射波波前法向矢量n 的叉乘,若偏振矢量P 與入射面法向量N 的點(diǎn)乘為正,則偏振矢量的符號(hào)為正。

    2.4 反射系數(shù)和透射系數(shù)的求解

    根據(jù)彈性波傳播理論,當(dāng)入射P 波入射到界面時(shí),應(yīng)滿足邊界的位移連續(xù)條件和應(yīng)力連續(xù)條件。根據(jù)Rokhlin[14]給出的邊界條件可以導(dǎo)出反射系數(shù)和透射系數(shù)計(jì)算公式,即

    式中

    3 模型分析

    為了分析不同類型各向異性參數(shù)對(duì)反射系數(shù)的影響,先將正交各向異性分別退化到用單組裂縫弱度參數(shù)描述的HTI介質(zhì)和用Thomsen參數(shù)描述的VTI介質(zhì)。然后,分兩種情況進(jìn)行數(shù)值分析:①固定裂縫弱度參數(shù)不變,變化背景介質(zhì)的VTI各向異性參數(shù);②固定背景VTI各向異性參數(shù)不變,變化裂縫弱度參數(shù)。

    為了驗(yàn)證本文中正交各向異性反射系數(shù)計(jì)算方法的有效性,將計(jì)算結(jié)果與Rüger在正交各向異性介質(zhì)中x1Ox3與x2Ox3面的PP波近似反射系數(shù)進(jìn)行對(duì)比;再將退化的HTI和VTI介質(zhì)的反射系數(shù)與Rüger近似公式所得到的反射系數(shù)進(jìn)行對(duì)比。

    3.1 正交各向異性參數(shù)與VTI各向異性參數(shù)和裂縫弱度參數(shù)的關(guān)系

    Bakulin等[12]給出了正交各向異性參數(shù)與背景介質(zhì)的VTI各向異性參數(shù)以及裂縫弱度參數(shù)的關(guān)系為式中:下標(biāo)“b”表示VTI背景介質(zhì);g=V2S0/V2P0;ZN、ZV和ZH是無量綱參數(shù),分別表示裂縫的法向弱度、垂直切向弱度以及水平切向弱度[32]。

    3.2 P波從低阻抗介質(zhì)入射到高阻抗介質(zhì)

    假設(shè)模型的上層為低阻抗各向同性介質(zhì),下層為高阻抗的正交各向異性介質(zhì)。根據(jù)式(3),剛度系數(shù)矩陣可以用P 波垂直速度VP0、S 波垂直速度VS0、密度ρ以及各向異性參數(shù)表示。

    模型上層介質(zhì)參數(shù)為:VP0=3150m/s,VS0=1615m/s,ρ=2322kg/m3,ε(1)=ε(2)=γ(1)=γ(2)=δ(1)=δ(2)=δ(3)=0。

    模型下層介質(zhì)參數(shù)為:VP0=3310m/s,VS0=1697m/s,ρ=2351kg/m3,其各向異性參數(shù)如表1或表2所示。

    3.2.1 裂縫弱度參數(shù)固定不變,變化背景介質(zhì)的VTI各向異性參數(shù)

    首先,固定三個(gè)裂縫弱度參數(shù)不變,改變背景介質(zhì)的VTI各向異性參數(shù),分析其對(duì)反射系數(shù)的影響。表1 參數(shù)由式(23)計(jì)算得到,裂縫弱度參數(shù)ZN、ZV和ZH分別為0.215、0.120、0.090;VTI各向異性參數(shù)由弱變強(qiáng):模型A1的背景為各向同性介質(zhì),εb=0、γb=0、δb=0,此時(shí)正交各向異性介質(zhì)退化為HTI介質(zhì);模型A2 的背景為弱各向異性VTI介質(zhì),εb=0.1、γb=0.12、δb=0.07;模型A3的背景為強(qiáng)各向異性VTI介質(zhì),εb=0.3、γb=0.25、δb=0.15。

    固定θ=20°,φ 從0°變化到90°時(shí)模型A1、A2、A3的反射系數(shù)RPP、RPS1、RPS2曲線如圖2所示。可以看出:

    (1)PP波反射系數(shù)的方位各向異性較弱;從模型A1到模型A3,VTI各向異性的增強(qiáng)導(dǎo)致RPP增大,模型A2比模型A1增大了40%,模型A3比模型A1增大了86%。

    (2)反射系數(shù)RPS1方位各向異性很強(qiáng),且隨方位角增大而減??;從模型A1到模型A3,VTI各向異性的增強(qiáng)導(dǎo)致RPS1減小,最大時(shí)模型A2比模型A1減小了33%,模型A3比模型A1減小了83%;

    (3)反射系數(shù)RPS2也具有很強(qiáng)的方位各向異性,曲線變化趨勢(shì)與RPS1沿著方位呈大體鏡像關(guān)系,但是量值上有差異;從模型A1到模型A3,VTI各向異性的增強(qiáng)導(dǎo)致RPS2減小,最大時(shí)模型A2 比模型A1減小了30%,模型A3比模型A1減小了70%。

    表1 裂縫弱度參數(shù)不變、背景介質(zhì)的VTI各向異性參數(shù)變化時(shí)的正交各向異性參數(shù)

    固定φ=60°,模型A1、A2、A3的反射系數(shù)RPP、RPS1、RPS2隨θ的變化曲線如圖3所示。由圖可知:

    (1)小入射角時(shí)(θ<10°),反射系數(shù)RPP受背景介質(zhì)VTI各向異性參數(shù)的影響不大。但是,隨著θ的增大,模型A1 的RPP隨θ 增大而減小,且在θ≈35°時(shí)發(fā)生極性轉(zhuǎn)變;從模型A1到模型A3,VTI各向異性增強(qiáng)導(dǎo)致大入射角的RPP增大。

    (2)隨著θ的增大,反射系數(shù)RPS1在θ<30°時(shí)都呈增大的趨勢(shì),而在30°~45°范圍內(nèi)模型A2、A3開始減小,模型A3在θ≈35°出現(xiàn)了極性轉(zhuǎn)變;從模型A1到模型A3,VTI各向異性參數(shù)的增大導(dǎo)致RPS1減小,在θ=30°時(shí)模型A2比模型A1減小了50%,模型A3比模型A1減小了75%。

    (3)隨著θ的增大,反射系數(shù)RPS2與RPS1變化趨勢(shì)一致;從模型A1到模型A3,VTI各向異性系數(shù)的增大導(dǎo)致RPS2減小,在θ=30°時(shí)模型A2 比模型A1減小了37.5%,模型A3比模型A1減小了75%。

    φ 在0°~90°內(nèi)變化、θ在臨界角范圍以內(nèi),P波由各向同性介質(zhì)入射到正交各向異性模型A1、A2和A3時(shí)的反射系數(shù)RPP、RPS1和RPS2的等值線圖分別如圖4、圖5、圖6所示,可以看出:

    (1)隨著VTI背景各向異性程度的增強(qiáng),RPP等值線曲率變化不大,只在反射系數(shù)的量值上有整體增大,且靠近臨界角時(shí)增加幅度最大。

    (2)RPS1和RPS2隨方位變化的程度遠(yuǎn)大于RPP,都有隨著背景介質(zhì)VTI各向異性強(qiáng)度增強(qiáng)而減小的趨勢(shì),且等值線曲率變化也較大。

    由于背景介質(zhì)中插入的裂縫弱度值不大,裂縫誘導(dǎo)的HTI各向異性屬于弱各向異性。針對(duì)此類正交各向異性介質(zhì)的AVO 分析,單純的PP波對(duì)裂縫的發(fā)育以及背景介質(zhì)的VTI各向異性程度的敏感度遠(yuǎn)小于PS1波和PS2波。

    圖3 φ=60°時(shí)模型A1、A2和A3的R PP(a)、R PS1(b)、R PS2(c)隨θ的變化曲線

    圖4 入射P波由各向同性到各向異性模型A1(a)、A2(b)、A3(c)反射系數(shù)R PP等值線圖

    圖5 入射P波由各向同性到各向異性模型A1(a)、A2(b)、A3(c)反射系數(shù)R PS1等值線圖

    圖6 入射P波由各向同性到各向異性模型A1(a)、A2(b)、A3(c)反射系數(shù)R PS2等值線圖

    3.2.2 背景介質(zhì)VTI各向異性系數(shù)固定不變,變化裂縫弱度參數(shù)

    由式(23)可知,當(dāng)固定背景VTI介質(zhì)的三個(gè)各向異性參數(shù)不變(εb=0.1,γb=0.12,δb=0.07,與表1中模型A2的背景介質(zhì)VTI各向異性參數(shù)取值相同)時(shí),改變?nèi)齻€(gè)裂縫弱度值會(huì)影響HTI各向異性的強(qiáng)弱,弱度越大則HTI各向異性越強(qiáng)。表2為下層高阻抗介質(zhì)的背景VTI各向異性參數(shù)固定不變、變化裂縫弱度參數(shù)計(jì)算的正交各向異性參數(shù)值。模型B1的裂縫弱度參數(shù)ZN=ZV=ZH=0,此時(shí)正交各向異性介質(zhì)退化為VTI各向異性介質(zhì);模型B2的裂縫弱度參數(shù)ZN=0.215、ZV=0.120、ZH=0.090,呈弱HTI各向異性,等同于表1中的模型A2;模型B3的裂縫弱度參數(shù)ZN=0.215、ZV=0.120、ZH=0.090,裂縫發(fā)育程度較高,HTI各向異性也較強(qiáng)。

    固定θ=20°,φ 在0°~90°變化時(shí)模型B1、B2、B3的反射系數(shù)RPP、RPS1、RPS2曲線如圖7所示。由圖7可知:

    (1)裂縫的發(fā)育會(huì)導(dǎo)致PP 波的反射系數(shù)變小,模型B2比模型B1減小了32%,模型B3比模型B1減小了64%;即使在裂縫弱度較強(qiáng)的情況下,PP 波的反射系數(shù)的方位各向異性也是比較微弱的。

    表2 背景VTI各向異性參數(shù)不變、裂縫弱度參數(shù)變化時(shí)的正交各向異性參數(shù)

    圖7 θ=20°時(shí)模型B1、B2和B3的反射系數(shù)R PP(a)、R PS1(b)、R PS2(c)隨φ 的變化曲線

    (2)反射系數(shù)RPS1受φ 的影響較大,在φ=0°時(shí),即平行裂縫面時(shí),RPS1最大;在φ=90°時(shí),即垂直裂縫面時(shí),RPS1為0;從模型B1到模型B3,裂縫弱度的增強(qiáng)導(dǎo)致RPS1增大,最大時(shí)模型B2比模型B1增加了300%,模型B3 比模型B1 增加了600%。當(dāng)裂縫不發(fā)育時(shí)(模型B1),RPS1不存在方位各向異性,等同于PSV 波的反射系數(shù)。

    (3)反射系數(shù)RPS2在φ=0°時(shí)最小,φ=90°時(shí)最大;從模型B1 到模型B3,裂縫弱度的增強(qiáng)導(dǎo)致RPS2增加,最大時(shí)模型B3比模型B2增加了100%。

    固定φ=60°時(shí),模型B1、B2、B3 的反射系數(shù)RPP、RPS1、RPS2隨θ 的 變 化 曲 線 如 圖8 所 示,由 圖可知:

    (1)θ 在0°~45°(臨界角以內(nèi))范圍內(nèi),反射系數(shù)RPP在θ<5°時(shí)無明顯差別,隨著θ增大,RPP表現(xiàn)不同,裂縫弱度越大,RPP減小越快;最大時(shí)模型B2比模型B1減小了75%,模型B3比模型B1減小了

    100%。

    (2)隨著θ增大,RPS1在θ<25°時(shí)都呈增大的趨勢(shì),而在25°~45°范圍內(nèi)模型B1、B2 開始減小,模型B1在θ≈30°出現(xiàn)了極性轉(zhuǎn)變;從模型B1到模型B3,裂縫弱度的增強(qiáng)導(dǎo)致RPS1增大,θ=25°時(shí)模型B2比模型B1增大了100%,模型B3比模型B1增大了300%。

    (3)隨著θ 增大,當(dāng)裂縫發(fā)育時(shí),RPS2都呈負(fù)極性方向增大的趨勢(shì),在θ=0°時(shí)最小,在θ=45°時(shí)最大;從模型B1 到模型B3,裂縫弱度的增強(qiáng)導(dǎo)致RPS2增大,最大時(shí)模型B3比模型B2增大了100%。

    φ 在0°~90°內(nèi)變化、θ在臨界角范圍以內(nèi),P波由各向同性介質(zhì)入射到正交各向異性模型B1、B2和B3時(shí)的反射系數(shù)RPP、RPS1和RPS2的等值線圖分別如圖9、圖10、圖11所示,可以看出:隨著裂縫弱度的增強(qiáng),RPP在正極性方向減小,RPS1和RPS2在負(fù)極性方向增大。總體上,即使裂縫誘導(dǎo)的HTI各向異性較強(qiáng),RPP的方位各向異性也比較弱,其對(duì)于裂縫發(fā)育的敏感度遠(yuǎn)小于RPS1和RPS2。

    圖8 φ=60°時(shí)模型B1、B2和B3的反射系數(shù)R PP(a)、R PS1(b)、R PS2(c)隨θ的變化曲線

    圖9 入射P波由各向同性到各向異性模型B1(a)、B2(b)和B3(c)的反射系數(shù)R PP等值線圖

    圖10 入射P波由各向同性到各向異性模型B1(a)、B2(b)和B3(c)的反射系數(shù)R PS1等值線圖

    圖11 入射P波由各向同性到各向異性模型B1(a)、B2(b)和B3(c)的反射系數(shù)R PS2等值線圖

    3.3 P波從高阻抗介質(zhì)入射到低阻抗介質(zhì)

    為了進(jìn)一步探討界面上、下阻抗對(duì)反射系數(shù)的影響,將A、B系列模型的上、下層介質(zhì)的縱橫波速度和密度交換,下層的各向異性參數(shù)保持不變,可得模型C1、C2、C3(表1,裂縫弱度不變,VTI各向異性強(qiáng)度變化)和模型D1、D2、D3(表2,裂縫弱度變化,VTI各向異性參數(shù)不變)。

    固定θ=20°,模型C1、C2、C3的反射系數(shù)RPP、RPS1、RPS2的隨φ 的變化曲線如圖12所示;固定φ=60°,模型C1、C2、C3的反射系數(shù)RPP、RPS1、RPS2隨θ的變化曲線如圖13所示。對(duì)比圖12、圖13與圖2、圖3可知,與下層介質(zhì)為高阻抗介質(zhì)的情況相反,當(dāng)下層介質(zhì)為低阻抗介質(zhì)時(shí),背景介質(zhì)的VTI各向異性強(qiáng)度越強(qiáng),則PP 波的反射系數(shù)在負(fù)極性方向越小,而PS1波與PS2波在正極性方向越大。

    固定θ=20°,模型D1、D2、D3的反射系數(shù)RPP、RPS1、RPS2隨φ 的變化曲線如圖14 所示;固定φ=60°,模型D1、D2、D3的反射系數(shù)RPP、RPS1、RPS2隨θ的變化曲線如圖15所示。對(duì)比圖14、圖15與圖7、圖8可知,裂縫的發(fā)育會(huì)導(dǎo)致PP 波的反射系數(shù)在負(fù)極性方向變大,而PS1波和PS2波的反射系數(shù)則在正極性方向變小,與下層介質(zhì)為高阻抗介質(zhì)的情況相反。

    無論是從低阻抗介質(zhì)入射,還是高阻抗介質(zhì)入射,PP波反射系數(shù)的方位各向異性都遠(yuǎn)小于PS1波和PS2波。

    3.4 與Rüger近似反射系數(shù)公式的對(duì)比

    為了驗(yàn)證本文反射系數(shù)計(jì)算結(jié)果的正確性,與Rüger近似公式計(jì)算的正交各向異性反射系數(shù)、HTI介質(zhì)反射系數(shù)和VTI介質(zhì)反射系數(shù)進(jìn)行對(duì)比。

    圖12 θ=20°時(shí)模型C1、C2和C3的反射系數(shù)R PP(a)、R PS1(b)、R PS2(c)隨φ 的變化曲線

    圖13 φ=60°時(shí)模型C1、C2和C3的反射系數(shù)R PP(a)、R PS1(b)、R PS2(c)隨θ的變化曲線

    圖14 θ=20°時(shí)模型D1、D2和D3的反射系數(shù)R PP(a)、R PS1(b)、R PS2(c)隨φ 的變化曲線

    圖15 φ=60°時(shí)模型D1、D2和D3的反射系數(shù)R PP(a)、R PS1(b)、R PS2(c)隨θ的變化曲線

    3.4.1 正交各向異性介質(zhì)

    對(duì)于P波入射,Rüger[17]僅推導(dǎo)了正交各向異性介質(zhì)RPP在x1Ox3和x2Ox3面的近似反射系數(shù)公式。將模型A2的參數(shù)代入到其公式計(jì)算反射系數(shù)曲線,并與本文的反射系數(shù)精確解的計(jì)算結(jié)果進(jìn)行對(duì)比(圖16)。在x1Ox3面,兩條曲線變化趨勢(shì)一致,當(dāng)θ在0°~20°范圍內(nèi)時(shí)近乎重合(圖16a);在x2Ox3面,兩條曲線變化趨勢(shì)亦一致,當(dāng)θ 在0°~25°范圍內(nèi)時(shí)近乎重合(圖16b)。

    3.4.2 退化到HTI介質(zhì)

    依據(jù)式(23)可以將正交各向異性介質(zhì)退化到HTI介質(zhì)。將模型A1代入Rüger[17]的HTI介質(zhì)反射系數(shù)近似公式計(jì)算反射系數(shù)曲線,并與本文的反射系數(shù)精確解對(duì)比(圖17)。由圖17 可以看出RPP、RPS1和RPS2的Rüger近似解與本文精確解變化趨勢(shì)一致;RPS2在0°~20°范圍內(nèi)近乎重合。

    3.4.3 退化到VTI介質(zhì)

    依據(jù)式(23)可以將正交各向異性介質(zhì)退化到VTI介質(zhì)。將模型B1 參數(shù)代入Rüger的VTI介質(zhì)反射系數(shù)近似公式,并與本文的反射系數(shù)精確解對(duì)比(圖18)。從圖18可以看出,RPP和RPS的近似解與精確解變化趨勢(shì)一致,RPP在0°~25°范圍內(nèi)完全重合,RPS在0°~15°內(nèi)近乎重合。

    圖16 正交各向異性介質(zhì)Rüger近似公式與本文方法計(jì)算的反射系數(shù)曲線對(duì)比

    圖17 HTI介質(zhì)Rüger近似公式與本文方法計(jì)算的反射系數(shù)曲線對(duì)比

    圖18 VTI介質(zhì)Rüger近似公式與本文方法計(jì)算的反射系數(shù)曲線對(duì)比

    由以上分析可知,本文反射系數(shù)精確解與Rüger近似公式在中、小角度范圍內(nèi)近乎重合,整體趨勢(shì)一致。因此,本文反射系數(shù)精確解計(jì)算方法正確。

    4 結(jié)論

    針對(duì)互層背景中發(fā)育一組裂縫的正交各向異性介質(zhì)模型,本文基于Christoffel方程及彈性界面的位移連續(xù)條件和應(yīng)力連續(xù)條件,推導(dǎo)了反射系數(shù)和透射系數(shù)精確求解方法,建立了正交各向異性參數(shù)與背景介質(zhì)VTI各向異性參數(shù)系數(shù)和裂縫參數(shù)之間的關(guān)聯(lián)。通過與Rüger反射系數(shù)近似公式計(jì)算結(jié)果的對(duì)比,驗(yàn)證了本文計(jì)算結(jié)果的正確性。理論模型計(jì)算結(jié)果表明:

    (1)當(dāng)P 波從低阻抗介質(zhì)入射至高阻抗介質(zhì)時(shí),背景介質(zhì)VTI各向異性參數(shù)的變大會(huì)導(dǎo)致PP波反射系數(shù)增大、PS1 波和PS2 波的反射系數(shù)減??;而裂縫弱度的增強(qiáng)會(huì)導(dǎo)致PP 波反射系數(shù)的減小及PS1波和PS2波的反射系數(shù)的增大。

    (2)當(dāng)從高阻抗介質(zhì)入射至低阻抗介質(zhì)時(shí),背景介質(zhì)VTI各向異性系數(shù)的變大會(huì)導(dǎo)致PP波反射系數(shù)減小、PS1波和PS2波的反射系數(shù)增大;而裂縫弱度的增強(qiáng)會(huì)導(dǎo)致PP 波反射系數(shù)的增大、PS1 波和PS2波的反射系數(shù)的減小。

    (3)PP 波反射系數(shù)的方位各向異性要遠(yuǎn)小于PS1波與PS2波。

    猜你喜歡
    反射系數(shù)偏振矢量
    矢量三角形法的應(yīng)用
    多道隨機(jī)稀疏反射系數(shù)反演
    石油物探(2020年6期)2020-11-25 02:38:46
    偏振糾纏雙光子態(tài)的糾纏特性分析
    電子制作(2019年12期)2019-07-16 08:45:20
    基于LabVIEW的偏振調(diào)制激光測(cè)距測(cè)量控制系統(tǒng)
    偏振旋轉(zhuǎn)效應(yīng)在全光緩存器中的應(yīng)用與實(shí)現(xiàn)
    球面波PP反射系數(shù)的頻變特征研究
    基于矢量最優(yōu)估計(jì)的穩(wěn)健測(cè)向方法
    三角形法則在動(dòng)態(tài)平衡問題中的應(yīng)用
    沙質(zhì)沉積物反射系數(shù)的寬帶測(cè)量方法
    基于反射系數(shù)的波導(dǎo)結(jié)構(gòu)不連續(xù)位置識(shí)別
    天堂俺去俺来也www色官网| 涩涩av久久男人的天堂| 人人妻人人澡人人爽人人夜夜| 狠狠婷婷综合久久久久久88av| 99香蕉大伊视频| 成人漫画全彩无遮挡| 国产免费视频播放在线视频| 人妻系列 视频| 丁香六月天网| 欧美xxⅹ黑人| 成人国产av品久久久| 精品国产露脸久久av麻豆| 国产熟女午夜一区二区三区| 久久国产亚洲av麻豆专区| 亚洲成人一二三区av| 亚洲国产看品久久| 亚洲精品aⅴ在线观看| 十分钟在线观看高清视频www| 国产成人免费无遮挡视频| 婷婷色av中文字幕| 欧美激情极品国产一区二区三区 | 18禁观看日本| 国产成人精品在线电影| 亚洲av男天堂| 一本大道久久a久久精品| 伊人久久国产一区二区| 国产免费现黄频在线看| 日韩大片免费观看网站| 欧美国产精品va在线观看不卡| 涩涩av久久男人的天堂| 满18在线观看网站| 高清av免费在线| 99久国产av精品国产电影| 丁香六月天网| 精品福利永久在线观看| 成人影院久久| 青春草国产在线视频| 国产又爽黄色视频| 国产成人精品久久久久久| 一本—道久久a久久精品蜜桃钙片| 精品卡一卡二卡四卡免费| 男人舔女人的私密视频| 精品酒店卫生间| 婷婷色综合大香蕉| 美女福利国产在线| 久久久久久久亚洲中文字幕| 人成视频在线观看免费观看| 青春草亚洲视频在线观看| 99热6这里只有精品| 欧美变态另类bdsm刘玥| 国产极品天堂在线| 久久久久久久精品精品| 丝瓜视频免费看黄片| 国产在线免费精品| 亚洲一级一片aⅴ在线观看| 欧美日韩国产mv在线观看视频| 丝瓜视频免费看黄片| 国产成人精品在线电影| 大片电影免费在线观看免费| 秋霞在线观看毛片| 欧美日本中文国产一区发布| 欧美日韩视频精品一区| 亚洲婷婷狠狠爱综合网| 免费少妇av软件| 男女免费视频国产| 少妇的逼好多水| 人妻少妇偷人精品九色| 久久亚洲国产成人精品v| 最新的欧美精品一区二区| av天堂久久9| 国产乱来视频区| 尾随美女入室| 精品一区二区三区四区五区乱码 | 免费看光身美女| 97精品久久久久久久久久精品| 十八禁高潮呻吟视频| 久久精品国产亚洲av涩爱| 狠狠婷婷综合久久久久久88av| 免费观看a级毛片全部| 精品久久蜜臀av无| 天美传媒精品一区二区| 99久国产av精品国产电影| 妹子高潮喷水视频| 丝袜脚勾引网站| 黄色 视频免费看| 韩国av在线不卡| 日本欧美视频一区| 欧美日本中文国产一区发布| a级毛片在线看网站| 日韩制服骚丝袜av| 国产精品女同一区二区软件| av国产久精品久网站免费入址| 韩国av在线不卡| 亚洲国产毛片av蜜桃av| 亚洲在久久综合| 久久久亚洲精品成人影院| xxx大片免费视频| 亚洲国产精品成人久久小说| 日韩成人伦理影院| 日本午夜av视频| 久久 成人 亚洲| 国产在线一区二区三区精| 九九在线视频观看精品| 精品一区在线观看国产| 久热久热在线精品观看| 久久人妻熟女aⅴ| av又黄又爽大尺度在线免费看| a级毛片黄视频| 午夜视频国产福利| 制服诱惑二区| 亚洲综合精品二区| 国产午夜精品一二区理论片| 婷婷色综合大香蕉| 国产精品嫩草影院av在线观看| 亚洲三级黄色毛片| 久久99蜜桃精品久久| 亚洲色图综合在线观看| 国产精品久久久av美女十八| 天天躁夜夜躁狠狠久久av| 中文字幕最新亚洲高清| 国产在线一区二区三区精| 亚洲第一区二区三区不卡| 国产国拍精品亚洲av在线观看| 在线免费观看不下载黄p国产| 久久亚洲国产成人精品v| 五月伊人婷婷丁香| 国产老妇伦熟女老妇高清| av视频免费观看在线观看| 久久99蜜桃精品久久| 久久国产亚洲av麻豆专区| 新久久久久国产一级毛片| 满18在线观看网站| 91在线精品国自产拍蜜月| 亚洲精品一区蜜桃| 欧美精品一区二区大全| av在线播放精品| 99久久人妻综合| 免费女性裸体啪啪无遮挡网站| 大片电影免费在线观看免费| 三上悠亚av全集在线观看| 欧美最新免费一区二区三区| 丰满乱子伦码专区| 亚洲图色成人| 亚洲激情五月婷婷啪啪| 亚洲婷婷狠狠爱综合网| 国产精品不卡视频一区二区| 亚洲综合色惰| 日韩一区二区三区影片| 国产男女超爽视频在线观看| 亚洲精品国产色婷婷电影| 亚洲成人一二三区av| 久久国产精品大桥未久av| 久久久国产一区二区| 麻豆乱淫一区二区| 精品国产露脸久久av麻豆| 亚洲欧洲日产国产| 自线自在国产av| 男女免费视频国产| 久久 成人 亚洲| 天天影视国产精品| 香蕉国产在线看| 精品国产国语对白av| 国产xxxxx性猛交| 日本与韩国留学比较| 国产精品久久久久久av不卡| 久久狼人影院| 精品国产一区二区久久| 国产乱来视频区| 国产精品一区二区在线观看99| 一级片'在线观看视频| 韩国高清视频一区二区三区| 欧美xxⅹ黑人| 国产综合精华液| 老熟女久久久| 日产精品乱码卡一卡2卡三| 高清不卡的av网站| 亚洲精品久久成人aⅴ小说| 日本av免费视频播放| 咕卡用的链子| 亚洲成av片中文字幕在线观看 | 97在线视频观看| 欧美亚洲日本最大视频资源| 满18在线观看网站| 美女福利国产在线| 免费看不卡的av| 国产极品粉嫩免费观看在线| 最近2019中文字幕mv第一页| 精品国产一区二区三区久久久樱花| 免费久久久久久久精品成人欧美视频 | 一区二区av电影网| 97人妻天天添夜夜摸| 99国产精品免费福利视频| av片东京热男人的天堂| 亚洲成人一二三区av| 一级毛片黄色毛片免费观看视频| 校园人妻丝袜中文字幕| 国产一区有黄有色的免费视频| 国产精品国产三级专区第一集| av国产久精品久网站免费入址| 另类亚洲欧美激情| 十分钟在线观看高清视频www| 日韩一区二区三区影片| 国产免费一级a男人的天堂| 最近中文字幕高清免费大全6| 天天影视国产精品| av免费观看日本| 久久久久网色| 国产免费现黄频在线看| 在线观看美女被高潮喷水网站| 免费黄频网站在线观看国产| 国产一区二区在线观看日韩| 亚洲精品第二区| 老司机影院毛片| 制服人妻中文乱码| 亚洲少妇的诱惑av| 69精品国产乱码久久久| 亚洲精品第二区| 国产高清不卡午夜福利| 久久久久久久精品精品| 亚洲综合精品二区| 99国产精品免费福利视频| 亚洲av福利一区| 香蕉精品网在线| 男女下面插进去视频免费观看 | 成年美女黄网站色视频大全免费| 视频中文字幕在线观看| 少妇被粗大的猛进出69影院 | 亚洲精品中文字幕在线视频| 成人毛片a级毛片在线播放| 久久国产亚洲av麻豆专区| 精品亚洲成国产av| 国产欧美日韩综合在线一区二区| 国产日韩欧美视频二区| 国产福利在线免费观看视频| 咕卡用的链子| 岛国毛片在线播放| 久久人人爽av亚洲精品天堂| h视频一区二区三区| 亚洲av成人精品一二三区| 久久人人爽人人片av| 你懂的网址亚洲精品在线观看| 久久综合国产亚洲精品| 午夜视频国产福利| 又粗又硬又长又爽又黄的视频| 亚洲国产成人一精品久久久| 91aial.com中文字幕在线观看| 在线 av 中文字幕| 亚洲成人av在线免费| 精品第一国产精品| 成人午夜精彩视频在线观看| 亚洲国产欧美在线一区| 精品亚洲成a人片在线观看| 中文欧美无线码| 亚洲精品日本国产第一区| 欧美激情极品国产一区二区三区 | 久久久国产精品麻豆| 久久ye,这里只有精品| 91成人精品电影| 日本午夜av视频| 亚洲欧美一区二区三区国产| 男女下面插进去视频免费观看 | 国产精品 国内视频| 99九九在线精品视频| 欧美激情国产日韩精品一区| 日韩av免费高清视频| 国产在线一区二区三区精| 国产精品无大码| 精品一区二区三区视频在线| 午夜福利在线观看免费完整高清在| 日本av手机在线免费观看| 久久狼人影院| 亚洲av福利一区| 久久精品国产亚洲av天美| 欧美亚洲 丝袜 人妻 在线| 一本大道久久a久久精品| 国产成人精品福利久久| 免费大片18禁| 热99久久久久精品小说推荐| 亚洲内射少妇av| 日日啪夜夜爽| 亚洲欧洲国产日韩| 亚洲国产成人一精品久久久| 免费观看a级毛片全部| 久久人人爽人人片av| 高清黄色对白视频在线免费看| 丝袜美足系列| 亚洲精品美女久久av网站| av一本久久久久| 丝袜喷水一区| 国产成人免费观看mmmm| 亚洲精品一区蜜桃| 国产一区二区在线观看av| 丝袜在线中文字幕| 亚洲熟女精品中文字幕| 999精品在线视频| 日韩人妻精品一区2区三区| 亚洲四区av| 久久久久久久久久人人人人人人| 大片免费播放器 马上看| 国产毛片在线视频| 国产精品一区二区在线不卡| 亚洲精华国产精华液的使用体验| 欧美少妇被猛烈插入视频| 亚洲国产毛片av蜜桃av| 少妇的丰满在线观看| 天天躁夜夜躁狠狠躁躁| 最新中文字幕久久久久| 欧美 日韩 精品 国产| 大香蕉久久成人网| 国产日韩欧美视频二区| 欧美亚洲 丝袜 人妻 在线| 水蜜桃什么品种好| 中文天堂在线官网| 亚洲欧美中文字幕日韩二区| 久久99一区二区三区| 最后的刺客免费高清国语| 大香蕉久久网| 亚洲精品久久午夜乱码| 99国产精品免费福利视频| 成人影院久久| 欧美激情 高清一区二区三区| 精品久久国产蜜桃| 国产xxxxx性猛交| 一本—道久久a久久精品蜜桃钙片| av线在线观看网站| 少妇人妻 视频| 丝袜在线中文字幕| av福利片在线| 欧美成人午夜免费资源| 精品久久久久久电影网| 高清在线视频一区二区三区| 人人妻人人澡人人看| 久热久热在线精品观看| 亚洲天堂av无毛| 三级国产精品片| 国产精品 国内视频| av女优亚洲男人天堂| 寂寞人妻少妇视频99o| 又黄又粗又硬又大视频| 寂寞人妻少妇视频99o| 中文字幕精品免费在线观看视频 | 亚洲欧美中文字幕日韩二区| 国产精品 国内视频| 中文乱码字字幕精品一区二区三区| 免费女性裸体啪啪无遮挡网站| 最后的刺客免费高清国语| 亚洲人成77777在线视频| 卡戴珊不雅视频在线播放| www.av在线官网国产| 校园人妻丝袜中文字幕| 国产成人精品福利久久| 大片免费播放器 马上看| 看十八女毛片水多多多| 好男人视频免费观看在线| 性高湖久久久久久久久免费观看| 美女内射精品一级片tv| 色网站视频免费| av国产久精品久网站免费入址| 亚洲成色77777| 久久 成人 亚洲| 国产在线免费精品| 久久午夜福利片| 亚洲成人av在线免费| 精品国产国语对白av| 国产男女超爽视频在线观看| 中文字幕精品免费在线观看视频 | 国产成人aa在线观看| 国产永久视频网站| 少妇熟女欧美另类| 成人手机av| 欧美丝袜亚洲另类| 26uuu在线亚洲综合色| 亚洲精品国产av蜜桃| 看非洲黑人一级黄片| 狠狠精品人妻久久久久久综合| 看非洲黑人一级黄片| 午夜精品国产一区二区电影| 亚洲伊人久久精品综合| 亚洲国产毛片av蜜桃av| 色哟哟·www| 人妻 亚洲 视频| 欧美少妇被猛烈插入视频| 午夜免费男女啪啪视频观看| 观看美女的网站| 捣出白浆h1v1| 男男h啪啪无遮挡| 久久女婷五月综合色啪小说| 国产在视频线精品| 大话2 男鬼变身卡| 黄色视频在线播放观看不卡| 蜜臀久久99精品久久宅男| 三级国产精品片| 欧美精品一区二区免费开放| 女人被躁到高潮嗷嗷叫费观| 精品一区二区免费观看| 久久国产精品男人的天堂亚洲 | 国产精品国产av在线观看| √禁漫天堂资源中文www| 国产又爽黄色视频| 亚洲精品色激情综合| 高清在线视频一区二区三区| 亚洲成国产人片在线观看| 亚洲图色成人| 亚洲欧美成人综合另类久久久| 亚洲精品乱久久久久久| 精品一品国产午夜福利视频| 丰满迷人的少妇在线观看| 久久精品国产亚洲av涩爱| 中文字幕av电影在线播放| 岛国毛片在线播放| 美国免费a级毛片| 大话2 男鬼变身卡| 国产 精品1| 免费高清在线观看视频在线观看| 欧美 日韩 精品 国产| 亚洲欧洲国产日韩| 欧美精品高潮呻吟av久久| 国产精品人妻久久久影院| 美女中出高潮动态图| 丰满乱子伦码专区| 捣出白浆h1v1| 最后的刺客免费高清国语| 亚洲欧美一区二区三区黑人 | 久久精品aⅴ一区二区三区四区 | 亚洲国产最新在线播放| 午夜福利在线观看免费完整高清在| a级毛片黄视频| 这个男人来自地球电影免费观看 | 国产成人aa在线观看| 蜜臀久久99精品久久宅男| 香蕉精品网在线| 精品一区二区三区四区五区乱码 | 免费在线观看黄色视频的| 亚洲av综合色区一区| 亚洲激情五月婷婷啪啪| 亚洲欧美清纯卡通| tube8黄色片| 国产男女超爽视频在线观看| 一二三四中文在线观看免费高清| 亚洲国产日韩一区二区| 成人漫画全彩无遮挡| 亚洲精品久久成人aⅴ小说| 成人国产麻豆网| 九九爱精品视频在线观看| 草草在线视频免费看| 自线自在国产av| 亚洲天堂av无毛| 韩国av在线不卡| 美女大奶头黄色视频| 亚洲国产最新在线播放| 尾随美女入室| 777米奇影视久久| 波多野结衣一区麻豆| 一本久久精品| 肉色欧美久久久久久久蜜桃| 午夜91福利影院| 日韩精品免费视频一区二区三区 | 亚洲美女视频黄频| 九九爱精品视频在线观看| 少妇被粗大猛烈的视频| 欧美亚洲日本最大视频资源| 久久久亚洲精品成人影院| 婷婷色综合www| 亚洲成人av在线免费| 春色校园在线视频观看| 制服丝袜香蕉在线| 我要看黄色一级片免费的| 日本免费在线观看一区| 女性被躁到高潮视频| 免费女性裸体啪啪无遮挡网站| www日本在线高清视频| av有码第一页| 黄片播放在线免费| 久久久久久久亚洲中文字幕| 日日撸夜夜添| 亚洲第一av免费看| 99热网站在线观看| 欧美97在线视频| 插逼视频在线观看| 久久午夜福利片| videosex国产| 免费观看无遮挡的男女| 99re6热这里在线精品视频| 狂野欧美激情性xxxx在线观看| 一区二区日韩欧美中文字幕 | 九色亚洲精品在线播放| 天天操日日干夜夜撸| 最后的刺客免费高清国语| 久久久久久久久久成人| 欧美精品av麻豆av| 国产色婷婷99| 免费av中文字幕在线| 午夜视频国产福利| 中文精品一卡2卡3卡4更新| 狠狠精品人妻久久久久久综合| 国产精品一国产av| 欧美日韩av久久| 亚洲一码二码三码区别大吗| 性色avwww在线观看| 亚洲成人av在线免费| 两个人看的免费小视频| 日韩成人伦理影院| 制服人妻中文乱码| 国产成人午夜福利电影在线观看| 91久久精品国产一区二区三区| 人妻 亚洲 视频| 成人毛片a级毛片在线播放| av国产久精品久网站免费入址| 国产免费现黄频在线看| 久久 成人 亚洲| 成人免费观看视频高清| 国精品久久久久久国模美| 亚洲国产日韩一区二区| 免费av中文字幕在线| 久久国产精品男人的天堂亚洲 | 美国免费a级毛片| 中文字幕另类日韩欧美亚洲嫩草| 夜夜骑夜夜射夜夜干| 免费女性裸体啪啪无遮挡网站| 亚洲欧美清纯卡通| 国产成人精品久久久久久| 国产精品久久久久久精品电影小说| 国产av码专区亚洲av| 2021少妇久久久久久久久久久| 午夜福利网站1000一区二区三区| 欧美日韩av久久| 欧美亚洲 丝袜 人妻 在线| 丝袜在线中文字幕| 久久这里只有精品19| 国产日韩欧美亚洲二区| 日韩三级伦理在线观看| 久久久精品免费免费高清| 久久久国产一区二区| 秋霞伦理黄片| 国产精品三级大全| 久久久久网色| 国产一区二区三区综合在线观看 | 黄片无遮挡物在线观看| 中文乱码字字幕精品一区二区三区| 2022亚洲国产成人精品| 蜜桃国产av成人99| 亚洲美女视频黄频| 新久久久久国产一级毛片| xxx大片免费视频| 久久这里有精品视频免费| 青春草亚洲视频在线观看| 精品少妇黑人巨大在线播放| 亚洲成人手机| 国产成人精品无人区| 97精品久久久久久久久久精品| 纵有疾风起免费观看全集完整版| 97精品久久久久久久久久精品| 欧美日韩视频精品一区| 精品亚洲乱码少妇综合久久| 黑人高潮一二区| 成年动漫av网址| 国产日韩一区二区三区精品不卡| 亚洲在久久综合| 黄色怎么调成土黄色| videosex国产| 精品一区二区三区视频在线| 亚洲精品久久久久久婷婷小说| 黄色 视频免费看| 亚洲伊人久久精品综合| 青春草视频在线免费观看| 天堂8中文在线网| 亚洲中文av在线| 午夜福利影视在线免费观看| 视频在线观看一区二区三区| 晚上一个人看的免费电影| 美女国产视频在线观看| 欧美97在线视频| 制服人妻中文乱码| 亚洲成色77777| a级毛片黄视频| 91aial.com中文字幕在线观看| 日日啪夜夜爽| av一本久久久久| 成年人免费黄色播放视频| 亚洲av免费高清在线观看| 亚洲精品久久久久久婷婷小说| 三级国产精品片| 国产国语露脸激情在线看| 午夜福利视频在线观看免费| 国产男女内射视频| 国产午夜精品一二区理论片| 免费大片18禁| 亚洲人成77777在线视频| 亚洲欧美中文字幕日韩二区| 免费观看a级毛片全部| 久久久a久久爽久久v久久| 国产视频首页在线观看| 久久久久国产网址| 婷婷色综合www| 超色免费av| 国产在线一区二区三区精| 国产国拍精品亚洲av在线观看| 人人妻人人爽人人添夜夜欢视频| 在线观看免费高清a一片| 亚洲国产看品久久| 丰满迷人的少妇在线观看| 大香蕉久久网| 亚洲精品日本国产第一区| av网站免费在线观看视频| 蜜桃国产av成人99| 桃花免费在线播放| 考比视频在线观看| 在线观看免费高清a一片| av网站免费在线观看视频| 女人被躁到高潮嗷嗷叫费观| 国产成人精品婷婷| 国产成人一区二区在线| 18禁在线无遮挡免费观看视频| 美女大奶头黄色视频| 天天躁夜夜躁狠狠躁躁|