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

    交錯網(wǎng)格有限差分正演模擬的聯(lián)合吸收邊界

    2018-09-20 11:50:46胡建林宋維琪張建坤邢文軍徐文會
    石油地球物理勘探 2018年5期
    關(guān)鍵詞:波場邊界條件聲波

    胡建林 宋維琪 張建坤 邢文軍 徐文會

    (①中國石油大學(華東)地球科學與技術(shù)學院,山東青島 266555; ②中國石油冀東油田公司勘探開發(fā)研究院,河北唐山 063004)

    1 引言

    復(fù)雜地下介質(zhì)中,地震波的傳播過程繁冗,難以得到解析解,因此,一般是通過正演模擬探究地下地震波的傳播。在地震波正演模擬中,利用波動方程的正演模擬比用運動學的射線追蹤法可獲得更豐富的動力學信息,因此地震波場的數(shù)值模擬是地震波場傳播研究中的重要手段之一[1-8]。

    對聲波方程進行網(wǎng)格離散時,采用交錯網(wǎng)格的差分格式,可較好地壓制數(shù)值頻散,得到更高精度的正演模擬結(jié)果[2]。交錯網(wǎng)格是指把速度及壓力P分別定義于兩套不同網(wǎng)格上的網(wǎng)格系統(tǒng),該方法能更好地處理一階波動方程中速度和壓力的耦合關(guān)系。

    在數(shù)值模擬時,邊界問題就是找到一種合適的邊界條件,使人為產(chǎn)生的計算區(qū)域邊界在實際效果上是透明的,即不產(chǎn)生邊界反射[3]。吸收邊界條件正是一種產(chǎn)生反射非常小甚至不產(chǎn)生反射的人工邊界條件。Smith[9]將Dirichlet邊界條件與Neumann邊界條件結(jié)合起來,由于兩種邊界條件產(chǎn)生的反射符號相反,兩者相加抵消反射,但其計算量較大; Cerjan等[10]、Burns[11]提出基于擴大計算區(qū)域的吸收邊界條件,該方法在擴大后的邊界區(qū)使用阻尼機制減少反射波,以較大的計算量換取邊界的吸收; Kosloff等[12]提出黏性阻尼方法以減少邊界反射,但需進行波數(shù)域的變換。Clayton等[13]在網(wǎng)格邊界附近用單向波動方程模擬聲波和彈性波的傳播,該方法對于入射角小于45°的入射波很有效; Liao等[14]提出邊界外推方法,用來處理聲波和彈性波的邊界吸收問題; 基于頻率域內(nèi)聲波和彈性波的吸收邊界條件,Randall[15]提出分別處理脹縮波和旋轉(zhuǎn)波分量的方法。Keys[16]和Higdon[17]在Clayton等[13]研究的基礎(chǔ)上提出一種改進的吸收邊界條件,成功消除了任意角入射波的邊界反射,并具有良好穩(wěn)定性。Bécache等[18]研究了時間和空間上的m階Higdon邊界體條件,發(fā)現(xiàn)當階數(shù)m較大時,Higdon邊界將變得異常復(fù)雜。

    Berenger[19]在使用時域有限差分法時,為截斷二維網(wǎng)格空間,提出PML吸收邊界概念; Chew等[20]將PML引入彈性波模擬,采用擴展坐標系的PML方法; Rodent等[21]進一步提出CPML(convolutional complex frequency-shifted PML)邊界條件,該邊界在非線性、各向異性介質(zhì)時效果更好; Gedney等[22]提出ADE-PML(the PML with Auxiliary Differential Equations)邊界條件,該方法可將差分擴展到更高時間階數(shù); Liu等[23,24]提出一種混合吸收邊界,采用一種過渡帶方式從內(nèi)部區(qū)域逐漸過渡到邊界,進一步削弱了邊界反射; Gao等[25]回顧了大部分典型的聲波方程邊界條件,并通過數(shù)值實驗對比了各種方法。

    至于將PML邊界與Higdon邊界的聯(lián)合使用,尚未找到相關(guān)文獻。本文將PML邊界與Higdon邊界構(gòu)成聯(lián)合吸收邊界同時使用,比單獨使用PML邊界減小了邊界厚度,也比單獨使用Higdon邊界具有更好的吸收效果。

    2 三維聲波方程交錯網(wǎng)格有限差分正演模擬

    交錯網(wǎng)格有限差分法與常規(guī)網(wǎng)格有限差分法相比,計算效率更高,且能進一步提高數(shù)值模擬的運算精度,有效壓制數(shù)值頻散。三維聲波方程交錯網(wǎng)格有限差分正演模擬主要由傳播方程、震源、邊界條件三部分組成。

    2.1 傳播方程

    在交錯網(wǎng)格有限差分正演中,使用一階壓力—速度方程。在非均勻各向同性介質(zhì)中,三維聲波方程中的一階壓力—速度方程為

    (1)

    (2)

    式中:P為壓力;vx、vy、vz代表x、y、z方向的速度分量;V為縱波速度(聲波速度);ρ為介質(zhì)密度。

    將式(1)、式(2)有限差分離散,即可得到一階聲波方程的時域迭代計算。

    2.2 震源

    震源子波函數(shù)有雷克子波、高斯子波、正弦衰減子波等。本文選用雷克子波

    f(t)=A[1-2(πf0t)2]exp[-(πf0t)2]

    (3)

    式中f0為子波主頻。

    加入震源后,在震源處,式(1)變?yōu)?/p>

    (4)

    式中f(t)為震源子波函數(shù),其他參數(shù)同式(1)。

    2.3 邊界條件

    在數(shù)值模擬中,只能在有限的區(qū)域范圍內(nèi)進行模擬,如果不加任何處理,即在邊界處壓力值和速度值默認設(shè)置為0,那么地震波到達人工邊界處時,將會產(chǎn)生非常強的反射,嚴重干擾內(nèi)部波場。為了減少人工邊界產(chǎn)生的反射波,需要再添加吸收邊界。

    人工吸收邊界主要有兩種吸收方式:匹配吸收邊界條件(最外層單層匹配吸收)和衰減吸收邊界層(有一定厚度的多層衰減)。文中聯(lián)合吸收邊界中的Higdon邊界屬于前者,而PML吸收邊界則是后者。這兩種邊界的具體討論見下文。

    除了上面的三個交錯網(wǎng)格有限差分的重要組成部分,有限差分離散時還需對離散格式的穩(wěn)定性和頻散進行分析。由于這部分并不是本文討論的重點,所以略過。

    3 PML吸收邊界

    PML吸收邊界通過在垂直邊界方向引入阻尼因子,選取適當?shù)淖枘釁?shù),使向邊界外傳播的波在匹配層中不斷衰減,從而使反射波最大程度上得到削弱,實現(xiàn)消除邊界反射的目的。該方法對一階壓力—速度方程時間域進行分裂,得到含有吸收衰減因子的吸收邊界方程組

    (5)

    P=Px+Py+Pz

    (6)

    (7)

    式中:dx、dy、dz分別為x、y、z方向的衰減因子,計算如下[26]

    (8)

    (9)

    式中:L為PML邊界總厚度,lx、ly、lz為PML邊界內(nèi)的點到最外層點的距離;d0為吸收系數(shù)的初始值;R=1.0×10-6為理論反射系數(shù)[27];β為控制PML吸收邊界的吸收效果參數(shù),常常設(shè)置為匹配層內(nèi)的最大縱波速度。

    PML吸收邊界為有一定厚度的多層邊界衰減,以沿-x軸方向傳播的vx和Px分量為例,吸收邊界的模式如圖1所示。

    圖1 PML邊界層沿-x軸方向傳播的vx和Px分量吸收示意圖線的粗細代表能量大小

    由圖1可見,要使PML吸收層達到很好的吸收效果,可將PML層加厚。這種簡單增加層厚度就能實現(xiàn)好的吸收效果是PML邊界的優(yōu)點,同時也會出現(xiàn)計算量問題。由于厚度的增加,使三維模型明顯胖了一圈,導(dǎo)致計算量更為龐大。

    4 Higdon吸收邊界條件

    Higdon[17,28]提出了具有下面形式的吸收邊界條件(以x軸方向上的邊界為例)

    (10)

    式中:m為吸收邊界條件的階數(shù);c是入射波傳播速度;αj是吸收入射角,第j個算子能完全吸收入射角為±αj的波。式中的反射系數(shù)為

    (11)

    式中θ為入射角。

    m=1時,為一階Higdon吸收邊界,式(11)變?yōu)?/p>

    (12)

    m=2時,為二階Higdon吸收邊界,就有

    (13)

    在差分網(wǎng)格上離散后,吸收邊界差分算子為

    (14)

    式中: 系數(shù)βj=cosαj,為正數(shù); Δt、Δx分別為時間步長和空間步長;I為單位算子,滿足

    (15)

    Higdon吸收邊界為最外層邊界匹配吸收,以沿-x軸方向傳播的vx和Px分量為例,吸收邊界的模式如圖2。

    Higdon邊界僅有單層厚度,不會明顯增加計算量,且能消除任意角度入射波的邊界反射,具有良好的穩(wěn)定性。但要想得到更好的吸收效果,必須使用高階的Higdon吸收邊界形式(m>2),而高階時離散化后過于復(fù)雜,這是Higdon邊界的一大缺陷。

    圖2 Higdon邊界沿-x軸方向傳播的vx和Px分量吸收示意圖

    5 聯(lián)合吸收邊界

    5.1 聯(lián)合吸收邊界方法

    由上文可知,PML吸收效果是通過一定厚度邊界中的垂直分量衰減來吸收向邊界外傳播的波,而Higdon邊界是根據(jù)偏微分算子消除一定角度的入射波。根據(jù)PML吸收邊界和Higdon邊界的吸收原理,結(jié)合兩種吸收邊界的優(yōu)劣,研究了新的聯(lián)合吸收邊界,使向邊界傳播的波先經(jīng)過PML邊界進行衰減,然后到達最外層Higdon邊界時,再對殘余的波進行匹配吸收,最大限度地削弱邊界反射。以x軸方向的左邊垂直分量的邊界為例,吸收邊界的模式如圖3所示。

    圖3 聯(lián)合吸收邊界上沿-x軸方向傳播的vx和Px分量吸收示意圖線的粗細代表能量大小

    由于在PML層最外層聯(lián)合使用了Higdon吸收邊界,因此在相同吸收效果下,可減小PML層厚度,從而減少計算量。

    5.2 聯(lián)合吸收邊界方程推導(dǎo)

    針對研究的聯(lián)合吸收邊界方法,對新的吸收邊界方程進行推導(dǎo)。為了最大程度簡化推導(dǎo),先考慮在PML吸收邊界內(nèi)沿x軸方向傳播的平面波,此時PML邊界內(nèi)的一階聲波方程(式(5)~式(7))簡化為

    (16)

    設(shè)vx0、Px0分別為分量vx、Px的幅值,則PML介質(zhì)中傳播的波場可表示為[29]

    (17)

    將式(17)代入式(16),可解得

    (18)

    將式(18)代入式(17)可得PML介質(zhì)中傳播的波場

    (19)

    將式(19)寫為更一般的平面波形式,可得PML邊界中沿著坐標軸傳播的平面波方程,以速度V沿x軸向x>0空間運動的平面波表示為

    (20)

    由上式可推導(dǎo)得到x=0邊界上的Higdon形式的吸收邊界條件算子

    (21)

    將上面的算子應(yīng)用于沿x軸傳播的平面波時,反射波場的值該為零。因此在x=0邊界上,PML吸收層內(nèi),m階的Higdon形式的吸收邊界條件可寫為

    (22)

    同理,可對波場區(qū)域的其他邊界上的PML層內(nèi)的Higdon形式吸收邊界進行類似的推導(dǎo),得到對應(yīng)的Higdon邊界條件。

    5.3 聯(lián)合吸收邊界算法及其計算步驟

    將推導(dǎo)得到的聯(lián)合吸收邊界方程,加入到有限差分離散方程中,通過差分算法的研究,形成了聯(lián)合吸收邊界的交錯網(wǎng)格有限差分算法,算法的波場更新計算如下:

    一、根據(jù)三維聲波方程中的一階壓力—速度方程(式(5))對邊界內(nèi)部與x,y,z有關(guān)的壓力分量單獨進行波場更新;

    二、對PML邊界內(nèi)的壓力分量采用含衰減因子的PML吸收邊界方程組(式(5)、式(8))進行邊界內(nèi)壓力分量波場衰減更新;

    三、對最外層Higdon邊界進行壓力分量波場更新;

    四、使用P=Px+Py+Pz(式(6))對各方向的壓力分量進行求和;

    五、根據(jù)三維聲波方程中的一階壓力—速度方程(式(7))對邊界內(nèi)部與x,y,z有關(guān)的速度分量單獨進行波場更新;

    六、對PML邊界內(nèi)的壓力分量采用含衰減因子的PML吸收邊界方程組(式(7)、式(8))進行邊界內(nèi)速度分量波場衰減更新;

    七、對最外層的Higdon邊界采用Higdon邊界吸收算子(式(22))進行速度分量波場更新;

    根據(jù)上述的波場更新步驟,反復(fù)進行迭代,即可得到邊界吸收良好的聯(lián)合吸收邊界正演模擬。

    6 數(shù)值分析

    設(shè)計均勻各向同性的層狀介質(zhì)模型,介質(zhì)層數(shù)為7層,模型尺寸為500m×600m×700m,模型網(wǎng)格尺寸為5m×5m×5m,網(wǎng)格數(shù)為Nx×Ny×Nz=101×121×141。震源使用40Hz的雷克子波,設(shè)置位置為(250m,300m,0m),接收線過震源,且與y方向平行,同樣間隔5m,共101個檢波器。層狀模型和炮檢分布如圖4,各層參數(shù)分布如表1所示。

    使用高階交錯網(wǎng)格有限差分算法,對層狀介質(zhì)進行正演計算,當PML厚度足夠時,能夠較好地吸收邊界的反射。當PML邊界厚度為15個點,觀察270ms時壓力P的波場快照和480ms長度的炮集記錄,可見有明顯的邊界反射(圖5b、圖6b);在PML邊界外,再用一階Higdon邊界進行吸收,可見邊界反射得到明顯減弱(圖5c、圖6c);在PML邊界外,再用二階Higdon邊界進行吸收,可見邊界反射得到進一步減弱(圖5d、圖6d),此時吸收效果與30個點厚度的PML邊界(圖5a、圖6a)相近。通過結(jié)合Higdon邊界,彌補因邊界厚度不足導(dǎo)致的吸收不足,體現(xiàn)出聯(lián)合吸收邊界的優(yōu)越性。

    從表2正演模擬耗時統(tǒng)計可發(fā)現(xiàn),與單純的PML邊界相比,達到相同的吸收效果時,聯(lián)合吸收邊界大大減小了運算量,PML邊界(15點厚度)聯(lián)合二階Higdon邊界與PML邊界(30點厚度)相比,同樣吸收干凈,但耗時約為后者的62%,大大減少了計算量。

    圖4 層狀模型和炮檢分布示意圖

    層號層網(wǎng)格深度m速度m/s密度g/cm310~15025002.252150~25030002.263250~35035002.284350~45030002.275450~55036002.296550~65038002.307650~70036002.30

    圖5 270ms時接收線剖面上P分量波場快照

    (a)PML邊界(30個點邊界厚度); (b)PML邊界(15個點邊界厚度); (c)PML邊界(15個點邊界厚度)+一階Higdon邊界; (d)PML邊界(15個點邊界厚度)+二階Higdon邊界

    圖6 P分量炮集記錄

    (a)PML邊界(30個點邊界厚度); (b)PML邊界(15個點邊界厚度); (c)PML邊界(15個點邊界厚度)+一階Higdon邊界; (d)PML邊界(15個點邊界厚度)+二階Higdon邊界

    7 結(jié)論

    本文將PML邊界與Higdon邊界結(jié)合起來,得到聯(lián)合吸收邊界,對各向同性多層介質(zhì)的聲波波場模擬,并且通過對比分析,可發(fā)現(xiàn)聯(lián)合吸收邊界具有如下優(yōu)勢:

    (1)聯(lián)合吸收邊界相比單純的PML吸收邊界,可使用較小的PML邊界厚度,使模型大小明顯減小,較明顯地減少了計算量;

    (2)與Higdon邊界相比,聯(lián)合吸收邊界結(jié)合了PML的吸收特性,能夠更有效地吸收邊界反射,獲得更好的模擬結(jié)果。

    猜你喜歡
    波場邊界條件聲波
    一類帶有Stieltjes積分邊界條件的分數(shù)階微分方程邊值問題正解
    帶有積分邊界條件的奇異攝動邊值問題的漸近解
    彈性波波場分離方法對比及其在逆時偏移成像中的應(yīng)用
    愛的聲波 將愛留在她身邊
    中國寶玉石(2018年3期)2018-07-09 03:13:58
    聲波殺手
    自適應(yīng)BPSK在井下鉆柱聲波傳輸中的應(yīng)用
    交錯網(wǎng)格與旋轉(zhuǎn)交錯網(wǎng)格對VTI介質(zhì)波場分離的影響分析
    地震學報(2016年1期)2016-11-28 05:38:36
    基于Hilbert變換的全波場分離逆時偏移成像
    “聲波驅(qū)蚊”靠譜嗎
    旋轉(zhuǎn)交錯網(wǎng)格VTI介質(zhì)波場模擬與波場分解
    老女人水多毛片| 看黄色毛片网站| 乱码一卡2卡4卡精品| 国产精品麻豆人妻色哟哟久久 | kizo精华| 欧美高清成人免费视频www| 久久精品国产亚洲av涩爱 | 超碰av人人做人人爽久久| 舔av片在线| 一区二区三区免费毛片| 夜夜夜夜夜久久久久| 日韩欧美国产在线观看| 一个人看视频在线观看www免费| 九九在线视频观看精品| 欧美日韩在线观看h| 99久久精品热视频| 日本爱情动作片www.在线观看| 久久韩国三级中文字幕| 成人二区视频| 99精品在免费线老司机午夜| 日本熟妇午夜| 夫妻性生交免费视频一级片| 国产一区二区三区av在线 | 成人欧美大片| 亚洲av二区三区四区| 一进一出抽搐gif免费好疼| 亚洲真实伦在线观看| 男人舔女人下体高潮全视频| 偷拍熟女少妇极品色| videossex国产| 亚洲欧美中文字幕日韩二区| 久久久久久久久久久丰满| 看片在线看免费视频| 九草在线视频观看| 草草在线视频免费看| 最近中文字幕高清免费大全6| 麻豆乱淫一区二区| 国产色婷婷99| 日韩人妻高清精品专区| 最近视频中文字幕2019在线8| 亚洲欧美中文字幕日韩二区| 免费av观看视频| 日本免费一区二区三区高清不卡| 久久久精品94久久精品| 麻豆精品久久久久久蜜桃| 麻豆一二三区av精品| 边亲边吃奶的免费视频| 亚洲精品国产成人久久av| 亚洲欧美中文字幕日韩二区| 中文字幕av在线有码专区| 亚洲中文字幕一区二区三区有码在线看| 在线观看一区二区三区| 毛片一级片免费看久久久久| 国产在视频线在精品| 黄色日韩在线| 国产真实伦视频高清在线观看| 波多野结衣高清无吗| 日本黄色视频三级网站网址| 亚洲av.av天堂| 色综合亚洲欧美另类图片| 欧美日韩精品成人综合77777| 亚洲欧美精品专区久久| 天堂中文最新版在线下载 | 十八禁国产超污无遮挡网站| 欧美日韩乱码在线| 久久久成人免费电影| 久久国内精品自在自线图片| 国产午夜福利久久久久久| 欧美+日韩+精品| av在线老鸭窝| 免费不卡的大黄色大毛片视频在线观看 | 亚洲第一电影网av| 美女高潮的动态| 国产精品久久久久久精品电影小说 | 国产高潮美女av| 热99在线观看视频| 狠狠狠狠99中文字幕| 久99久视频精品免费| 欧美另类亚洲清纯唯美| 尾随美女入室| 国产亚洲5aaaaa淫片| 午夜福利在线在线| 国产精品一区二区三区四区免费观看| 日韩欧美 国产精品| 久久久久性生活片| 一区二区三区四区激情视频 | 精品免费久久久久久久清纯| 国产女主播在线喷水免费视频网站 | 特大巨黑吊av在线直播| 欧美色欧美亚洲另类二区| 国产精品麻豆人妻色哟哟久久 | 18禁在线无遮挡免费观看视频| www.av在线官网国产| 久久精品影院6| 日韩av不卡免费在线播放| 亚洲自拍偷在线| 我的女老师完整版在线观看| 欧美精品国产亚洲| 亚洲国产欧洲综合997久久,| 国产精品精品国产色婷婷| 亚洲高清免费不卡视频| 国产成人a∨麻豆精品| 美女内射精品一级片tv| 天天躁日日操中文字幕| 免费人成视频x8x8入口观看| 亚洲av熟女| 在线观看午夜福利视频| av卡一久久| 18禁裸乳无遮挡免费网站照片| 中文字幕免费在线视频6| 欧美日韩国产亚洲二区| 精品熟女少妇av免费看| 又粗又硬又长又爽又黄的视频 | 中文字幕精品亚洲无线码一区| 久久久精品大字幕| av天堂在线播放| 久久久久久国产a免费观看| 久久久久性生活片| 变态另类成人亚洲欧美熟女| 亚洲欧洲日产国产| 好男人在线观看高清免费视频| 国产美女午夜福利| 午夜精品国产一区二区电影 | 欧美潮喷喷水| 欧美精品一区二区大全| 免费不卡的大黄色大毛片视频在线观看 | 日韩在线高清观看一区二区三区| 91aial.com中文字幕在线观看| 国产精品福利在线免费观看| 色5月婷婷丁香| 床上黄色一级片| 国产精品99久久久久久久久| 国产精品久久电影中文字幕| 三级毛片av免费| 啦啦啦啦在线视频资源| 国产在线精品亚洲第一网站| 成人永久免费在线观看视频| 亚洲精品日韩av片在线观看| 中文精品一卡2卡3卡4更新| 丝袜美腿在线中文| 人妻久久中文字幕网| av在线亚洲专区| 一区二区三区高清视频在线| 成人特级av手机在线观看| 国产精品伦人一区二区| 寂寞人妻少妇视频99o| 久久久精品94久久精品| 日本三级黄在线观看| 日本-黄色视频高清免费观看| 97在线视频观看| 男人狂女人下面高潮的视频| 国产极品精品免费视频能看的| 亚洲最大成人手机在线| 99热6这里只有精品| 欧美一区二区亚洲| www.色视频.com| 亚洲自拍偷在线| 亚洲中文字幕一区二区三区有码在线看| 伦精品一区二区三区| 亚洲国产色片| 国产精品久久久久久av不卡| 黄色视频,在线免费观看| 国产私拍福利视频在线观看| www日本黄色视频网| 亚洲精品乱码久久久久久按摩| 在线观看免费视频日本深夜| 国产三级中文精品| 黄色日韩在线| 亚洲人成网站在线观看播放| 日本一本二区三区精品| 亚洲天堂国产精品一区在线| 女人被狂操c到高潮| 久久亚洲精品不卡| 亚洲三级黄色毛片| 国产久久久一区二区三区| 国产精品久久久久久久久免| 日韩一区二区三区影片| 最近2019中文字幕mv第一页| 18禁裸乳无遮挡免费网站照片| 综合色丁香网| or卡值多少钱| 综合色丁香网| 免费大片18禁| 国产黄a三级三级三级人| 日韩欧美一区二区三区在线观看| 亚洲成人中文字幕在线播放| 国产精品福利在线免费观看| 亚洲真实伦在线观看| 日韩一本色道免费dvd| 男的添女的下面高潮视频| 国产精品久久久久久久电影| 小说图片视频综合网站| 国产白丝娇喘喷水9色精品| 久久久久国产网址| 久久鲁丝午夜福利片| 国产免费男女视频| 久久久久免费精品人妻一区二区| 成人午夜精彩视频在线观看| av专区在线播放| 联通29元200g的流量卡| av在线老鸭窝| 中文字幕久久专区| 99热精品在线国产| 亚洲三级黄色毛片| 免费一级毛片在线播放高清视频| 久久久久久久亚洲中文字幕| 免费人成在线观看视频色| 青春草视频在线免费观看| 一个人看视频在线观看www免费| 十八禁国产超污无遮挡网站| 国产一级毛片在线| 中文字幕制服av| 国产久久久一区二区三区| 99久国产av精品| 国产v大片淫在线免费观看| 免费看日本二区| 嘟嘟电影网在线观看| 99久国产av精品| 日本三级黄在线观看| 欧洲精品卡2卡3卡4卡5卡区| 欧美高清性xxxxhd video| 精品久久久久久久久久久久久| 十八禁国产超污无遮挡网站| 联通29元200g的流量卡| 精品久久久久久久人妻蜜臀av| 国产伦精品一区二区三区四那| 亚洲三级黄色毛片| 小蜜桃在线观看免费完整版高清| 小蜜桃在线观看免费完整版高清| 91精品一卡2卡3卡4卡| 成人一区二区视频在线观看| 麻豆一二三区av精品| 尾随美女入室| 少妇高潮的动态图| 小说图片视频综合网站| 亚洲精品日韩在线中文字幕 | 久久久a久久爽久久v久久| 99热这里只有是精品50| 99精品在免费线老司机午夜| 在线播放无遮挡| 亚洲最大成人中文| 久久鲁丝午夜福利片| 国产精品一区二区在线观看99 | 亚洲中文字幕一区二区三区有码在线看| 日韩欧美 国产精品| 在线观看av片永久免费下载| 天堂中文最新版在线下载 | 久久久久网色| 99热全是精品| 欧美另类亚洲清纯唯美| 国产又黄又爽又无遮挡在线| 99热精品在线国产| 日韩国内少妇激情av| 国产午夜福利久久久久久| 国产爱豆传媒在线观看| 男女啪啪激烈高潮av片| 亚洲国产精品sss在线观看| 精品国产三级普通话版| 亚洲av成人av| 日本爱情动作片www.在线观看| 亚洲国产精品成人综合色| 欧美成人一区二区免费高清观看| 91精品一卡2卡3卡4卡| 99热全是精品| 一区二区三区四区激情视频 | 国产 一区 欧美 日韩| 99久久精品一区二区三区| 国产在线男女| 国产高清不卡午夜福利| 日产精品乱码卡一卡2卡三| 久久久久久久久大av| 国产精品爽爽va在线观看网站| 日本欧美国产在线视频| 久久精品国产99精品国产亚洲性色| 国产精品永久免费网站| 中国美白少妇内射xxxbb| 你懂的网址亚洲精品在线观看 | 国产熟女欧美一区二区| 亚洲精品亚洲一区二区| 99热网站在线观看| 一区二区三区四区激情视频 | 九九在线视频观看精品| 哪个播放器可以免费观看大片| 男女边吃奶边做爰视频| 日韩欧美精品免费久久| 亚洲精品乱码久久久v下载方式| 日韩在线高清观看一区二区三区| 亚洲中文字幕日韩| 欧美bdsm另类| 在线观看美女被高潮喷水网站| 精品人妻熟女av久视频| 又黄又爽又刺激的免费视频.| 精品少妇黑人巨大在线播放 | 色综合色国产| 自拍偷自拍亚洲精品老妇| 久久精品影院6| 美女内射精品一级片tv| 免费看av在线观看网站| 精品人妻一区二区三区麻豆| 九草在线视频观看| 人人妻人人澡欧美一区二区| 国国产精品蜜臀av免费| 最后的刺客免费高清国语| 一级av片app| 午夜视频国产福利| 国产午夜福利久久久久久| 精品久久久久久久久亚洲| 亚洲精品色激情综合| 最近2019中文字幕mv第一页| 日韩一区二区视频免费看| 欧美最新免费一区二区三区| 国产真实伦视频高清在线观看| 国产一区二区在线av高清观看| 国产成人影院久久av| 精品日产1卡2卡| 人人妻人人看人人澡| 国产精品久久久久久久电影| 国产伦精品一区二区三区视频9| 亚洲最大成人av| 中文字幕精品亚洲无线码一区| 91精品一卡2卡3卡4卡| 最好的美女福利视频网| 禁无遮挡网站| 久久人妻av系列| 精品久久久久久久人妻蜜臀av| 大香蕉久久网| 亚洲av中文av极速乱| 少妇的逼好多水| 国产精品三级大全| 性色avwww在线观看| 成人av在线播放网站| 99热6这里只有精品| 中文字幕久久专区| 九草在线视频观看| av免费观看日本| 亚洲国产欧美在线一区| 男女视频在线观看网站免费| 免费看日本二区| 亚洲国产日韩欧美精品在线观看| 亚洲精品久久久久久婷婷小说 | 麻豆成人av视频| 亚洲精品乱码久久久久久按摩| 久久久精品94久久精品| 国产亚洲av嫩草精品影院| 99热这里只有精品一区| 成人高潮视频无遮挡免费网站| 亚洲成人中文字幕在线播放| 搡老妇女老女人老熟妇| 99在线视频只有这里精品首页| 免费看日本二区| 天天一区二区日本电影三级| 亚洲欧美日韩无卡精品| 国产成年人精品一区二区| 一边摸一边抽搐一进一小说| 欧美成人免费av一区二区三区| 少妇熟女aⅴ在线视频| 村上凉子中文字幕在线| 成年女人永久免费观看视频| 最近2019中文字幕mv第一页| 十八禁国产超污无遮挡网站| 国语自产精品视频在线第100页| 2021天堂中文幕一二区在线观| 偷拍熟女少妇极品色| 久久精品久久久久久久性| 国产伦一二天堂av在线观看| 91av网一区二区| 狠狠狠狠99中文字幕| 成年av动漫网址| 亚洲色图av天堂| 国产蜜桃级精品一区二区三区| 久久午夜亚洲精品久久| 亚洲人成网站高清观看| 精品一区二区三区人妻视频| 哪里可以看免费的av片| 国产成人aa在线观看| 日韩成人伦理影院| 在线天堂最新版资源| 中文亚洲av片在线观看爽| 日韩一区二区三区影片| 国产美女午夜福利| 级片在线观看| 永久网站在线| 淫秽高清视频在线观看| 夜夜爽天天搞| 亚洲美女搞黄在线观看| 亚洲经典国产精华液单| 久久精品久久久久久噜噜老黄 | 99热网站在线观看| 中国美女看黄片| 久久久久久久亚洲中文字幕| 22中文网久久字幕| 成年女人永久免费观看视频| 久久99热6这里只有精品| 熟妇人妻久久中文字幕3abv| 男人舔女人下体高潮全视频| 一级毛片我不卡| av在线老鸭窝| 此物有八面人人有两片| 一个人免费在线观看电影| 直男gayav资源| av.在线天堂| 国产一区二区激情短视频| 非洲黑人性xxxx精品又粗又长| 麻豆成人av视频| 91av网一区二区| 六月丁香七月| 精品久久久噜噜| 国内少妇人妻偷人精品xxx网站| 色综合站精品国产| 成人二区视频| 婷婷亚洲欧美| 亚洲av免费在线观看| 最近的中文字幕免费完整| 亚洲精华国产精华液的使用体验 | 人妻系列 视频| 人妻制服诱惑在线中文字幕| 国产国拍精品亚洲av在线观看| 午夜久久久久精精品| 亚洲欧美精品专区久久| 两个人视频免费观看高清| 国产在线男女| 午夜福利高清视频| 嫩草影院入口| 又爽又黄无遮挡网站| 日韩国内少妇激情av| 菩萨蛮人人尽说江南好唐韦庄 | 村上凉子中文字幕在线| 97热精品久久久久久| 熟妇人妻久久中文字幕3abv| 黑人高潮一二区| 欧美另类亚洲清纯唯美| 夜夜看夜夜爽夜夜摸| 国产高清视频在线观看网站| 久久精品国产99精品国产亚洲性色| 亚洲国产日韩欧美精品在线观看| 精品午夜福利在线看| 免费在线观看成人毛片| 欧美日韩一区二区视频在线观看视频在线 | 99热全是精品| 你懂的网址亚洲精品在线观看 | 国产av不卡久久| 久久精品国产亚洲av天美| 最好的美女福利视频网| 国产精华一区二区三区| 美女被艹到高潮喷水动态| 国产真实伦视频高清在线观看| 国内精品宾馆在线| 国产av麻豆久久久久久久| 菩萨蛮人人尽说江南好唐韦庄 | 亚洲av二区三区四区| 我的女老师完整版在线观看| 成熟少妇高潮喷水视频| 此物有八面人人有两片| 亚洲欧美中文字幕日韩二区| 18+在线观看网站| 高清日韩中文字幕在线| 一卡2卡三卡四卡精品乱码亚洲| 2022亚洲国产成人精品| 国产精品久久久久久av不卡| 18+在线观看网站| 久久久久久国产a免费观看| 国内精品久久久久精免费| 亚洲一级一片aⅴ在线观看| 级片在线观看| 婷婷六月久久综合丁香| 美女 人体艺术 gogo| 在线观看av片永久免费下载| 国产片特级美女逼逼视频| 九九在线视频观看精品| 丝袜美腿在线中文| 精品少妇黑人巨大在线播放 | 丰满人妻一区二区三区视频av| 国产视频内射| 大型黄色视频在线免费观看| 性色avwww在线观看| 欧美激情国产日韩精品一区| 国产精品日韩av在线免费观看| 99热这里只有精品一区| 国产不卡一卡二| 麻豆成人午夜福利视频| 亚洲国产高清在线一区二区三| 国产成人a∨麻豆精品| 欧美人与善性xxx| 日日撸夜夜添| 亚洲图色成人| av国产免费在线观看| 国产黄色视频一区二区在线观看 | 大香蕉久久网| 少妇裸体淫交视频免费看高清| 青春草视频在线免费观看| 亚洲熟妇中文字幕五十中出| 亚洲第一电影网av| 91av网一区二区| 女人被狂操c到高潮| 中文欧美无线码| 少妇熟女欧美另类| 12—13女人毛片做爰片一| 一夜夜www| 欧美日韩国产亚洲二区| 国产一区亚洲一区在线观看| 中文在线观看免费www的网站| 日本五十路高清| 欧美xxxx黑人xx丫x性爽| 国产探花极品一区二区| 中国国产av一级| 国产成人精品一,二区 | 精品99又大又爽又粗少妇毛片| 亚洲真实伦在线观看| 午夜激情福利司机影院| 中文亚洲av片在线观看爽| 18禁在线无遮挡免费观看视频| 99久久精品一区二区三区| 丰满人妻一区二区三区视频av| 国产综合懂色| 丰满人妻一区二区三区视频av| 国产在视频线在精品| 99热网站在线观看| 一本精品99久久精品77| 国产高清激情床上av| 九九爱精品视频在线观看| 久久精品国产自在天天线| 国产av一区在线观看免费| 国内揄拍国产精品人妻在线| 久久久久性生活片| 人人妻人人澡欧美一区二区| 午夜激情福利司机影院| 日韩人妻高清精品专区| 春色校园在线视频观看| 日韩,欧美,国产一区二区三区 | 亚洲国产欧洲综合997久久,| 久久精品国产鲁丝片午夜精品| 国产蜜桃级精品一区二区三区| 国产亚洲av片在线观看秒播厂 | 中文字幕人妻熟人妻熟丝袜美| 日韩制服骚丝袜av| 午夜久久久久精精品| 国产精品一区二区三区四区久久| 亚洲av成人精品一区久久| 爱豆传媒免费全集在线观看| 又黄又爽又刺激的免费视频.| www.色视频.com| 国产在线精品亚洲第一网站| 国产精品,欧美在线| 国产高清三级在线| 99在线人妻在线中文字幕| 可以在线观看毛片的网站| 九九热线精品视视频播放| 久久久久久久久久成人| 91久久精品国产一区二区成人| 成年版毛片免费区| 99riav亚洲国产免费| 国产熟女欧美一区二区| 老司机影院成人| 国产三级中文精品| 国产单亲对白刺激| 搡女人真爽免费视频火全软件| 日本三级黄在线观看| 高清日韩中文字幕在线| 精品久久久噜噜| 五月玫瑰六月丁香| 久久久久久久久中文| 精品免费久久久久久久清纯| 国产毛片a区久久久久| 能在线免费看毛片的网站| 成人毛片a级毛片在线播放| 国产精品嫩草影院av在线观看| 我的老师免费观看完整版| 国产麻豆成人av免费视频| 国产精品福利在线免费观看| 免费黄网站久久成人精品| 免费电影在线观看免费观看| 国产精品野战在线观看| 丰满乱子伦码专区| 亚洲精华国产精华液的使用体验 | 亚洲久久久久久中文字幕| 亚洲欧美日韩无卡精品| 国产黄色视频一区二区在线观看 | 人人妻人人澡欧美一区二区| 国产v大片淫在线免费观看| 欧美bdsm另类| 久久人妻av系列| 国产精品三级大全| 欧美高清成人免费视频www| 国产精品,欧美在线| 成人午夜精彩视频在线观看| 国产一区二区亚洲精品在线观看| 特级一级黄色大片| 久久韩国三级中文字幕| 日本-黄色视频高清免费观看| 能在线免费观看的黄片| 美女大奶头视频| 亚洲第一电影网av| 成年女人永久免费观看视频| 精品少妇黑人巨大在线播放 | 欧美色视频一区免费| 欧美又色又爽又黄视频| 特大巨黑吊av在线直播| 亚州av有码| 中国美白少妇内射xxxbb| 卡戴珊不雅视频在线播放| 岛国毛片在线播放| 卡戴珊不雅视频在线播放| 人体艺术视频欧美日本| 最近2019中文字幕mv第一页| 亚洲精品久久久久久婷婷小说 | 美女xxoo啪啪120秒动态图| 免费av观看视频| 亚洲av熟女| 欧美激情国产日韩精品一区| 乱人视频在线观看| 天天一区二区日本电影三级| 好男人视频免费观看在线| 色哟哟哟哟哟哟| 啦啦啦啦在线视频资源| 国产高清激情床上av|