• <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ì)波場模擬與波場分解
    亚洲无线观看免费| 91在线观看av| 小蜜桃在线观看免费完整版高清| 欧美日韩综合久久久久久 | 麻豆成人午夜福利视频| 最近最新中文字幕大全电影3| 岛国在线观看网站| 男女午夜视频在线观看| 久久久久久久午夜电影| 两性夫妻黄色片| 丰满人妻熟妇乱又伦精品不卡| 不卡av一区二区三区| 好男人电影高清在线观看| 亚洲激情在线av| 亚洲欧美日韩卡通动漫| 欧美成人免费av一区二区三区| netflix在线观看网站| 午夜视频精品福利| 亚洲国产欧美人成| 99久久国产精品久久久| 亚洲av片天天在线观看| 久久欧美精品欧美久久欧美| 九九热线精品视视频播放| 嫁个100分男人电影在线观看| 久久草成人影院| 中文资源天堂在线| 久久欧美精品欧美久久欧美| 黄色视频,在线免费观看| 一进一出抽搐动态| 中文字幕高清在线视频| 久久久久亚洲av毛片大全| 国产午夜精品久久久久久| 美女cb高潮喷水在线观看 | 哪里可以看免费的av片| 亚洲人成电影免费在线| 午夜日韩欧美国产| 亚洲精品在线观看二区| 久久中文字幕人妻熟女| 欧美丝袜亚洲另类 | 欧美一区二区国产精品久久精品| 成人性生交大片免费视频hd| 美女免费视频网站| 一二三四社区在线视频社区8| 国产成人影院久久av| 精品久久久久久成人av| 又紧又爽又黄一区二区| 男人舔女人下体高潮全视频| 久久久久久久精品吃奶| 亚洲第一电影网av| 成人特级黄色片久久久久久久| 97人妻精品一区二区三区麻豆| 香蕉丝袜av| 亚洲人成电影免费在线| 中文字幕av在线有码专区| 午夜福利高清视频| 嫩草影院精品99| 国产成人av教育| 国产伦精品一区二区三区四那| 熟女电影av网| 一进一出好大好爽视频| 禁无遮挡网站| 观看免费一级毛片| 亚洲一区二区三区色噜噜| 久久中文看片网| 麻豆av在线久日| 欧洲精品卡2卡3卡4卡5卡区| 亚洲熟女毛片儿| 精品国产美女av久久久久小说| 久久天躁狠狠躁夜夜2o2o| 亚洲avbb在线观看| 亚洲熟女毛片儿| 国产91精品成人一区二区三区| 欧美成人一区二区免费高清观看 | 国产欧美日韩精品亚洲av| 国产亚洲av嫩草精品影院| 最近最新中文字幕大全电影3| 中文亚洲av片在线观看爽| 欧美午夜高清在线| 久9热在线精品视频| 国产精品电影一区二区三区| 国产三级在线视频| 精品免费久久久久久久清纯| 久久久久国产一级毛片高清牌| 国产精品乱码一区二三区的特点| 午夜福利欧美成人| 欧美乱妇无乱码| 精品久久久久久久久久免费视频| 久久久久久人人人人人| 免费大片18禁| 人妻久久中文字幕网| 欧美激情久久久久久爽电影| 国产精品久久久人人做人人爽| 国产一区二区三区在线臀色熟女| 好男人在线观看高清免费视频| 真实男女啪啪啪动态图| 老司机深夜福利视频在线观看| 最好的美女福利视频网| 国产精品 国内视频| 黄色丝袜av网址大全| 小说图片视频综合网站| 天堂√8在线中文| 别揉我奶头~嗯~啊~动态视频| 国产伦精品一区二区三区四那| 神马国产精品三级电影在线观看| 亚洲av熟女| svipshipincom国产片| 免费av不卡在线播放| 1024香蕉在线观看| 99在线视频只有这里精品首页| 可以在线观看毛片的网站| 日韩欧美国产一区二区入口| cao死你这个sao货| 啦啦啦观看免费观看视频高清| 12—13女人毛片做爰片一| 色综合站精品国产| 一卡2卡三卡四卡精品乱码亚洲| 国语自产精品视频在线第100页| 99久久精品热视频| 男插女下体视频免费在线播放| 日韩大尺度精品在线看网址| 日韩欧美国产一区二区入口| 欧美黑人巨大hd| 亚洲天堂国产精品一区在线| 精品电影一区二区在线| 亚洲七黄色美女视频| svipshipincom国产片| 精品一区二区三区视频在线观看免费| 国产精品av久久久久免费| 国产单亲对白刺激| 身体一侧抽搐| 露出奶头的视频| 性欧美人与动物交配| 国产乱人伦免费视频| 国产私拍福利视频在线观看| 久久久久性生活片| 国产亚洲av高清不卡| 久久久水蜜桃国产精品网| x7x7x7水蜜桃| 成人高潮视频无遮挡免费网站| 一进一出好大好爽视频| 亚洲成a人片在线一区二区| 天天一区二区日本电影三级| 久久这里只有精品中国| 禁无遮挡网站| 人人妻人人看人人澡| 一级a爱片免费观看的视频| 母亲3免费完整高清在线观看| 最近最新中文字幕大全电影3| 日韩免费av在线播放| 性色av乱码一区二区三区2| 国产精品 欧美亚洲| 久久精品91蜜桃| 久久久国产精品麻豆| www国产在线视频色| 精品不卡国产一区二区三区| 最近视频中文字幕2019在线8| 特大巨黑吊av在线直播| 日韩欧美 国产精品| 极品教师在线免费播放| av黄色大香蕉| 亚洲av成人不卡在线观看播放网| 两个人视频免费观看高清| 久久久久精品国产欧美久久久| 色哟哟哟哟哟哟| 舔av片在线| 蜜桃久久精品国产亚洲av| 精品一区二区三区视频在线观看免费| 日日干狠狠操夜夜爽| 国产高清videossex| 12—13女人毛片做爰片一| 欧美日韩乱码在线| 人妻丰满熟妇av一区二区三区| 欧美又色又爽又黄视频| 男插女下体视频免费在线播放| 中文字幕av在线有码专区| 一a级毛片在线观看| 久久性视频一级片| 日韩精品中文字幕看吧| 国产一区二区在线av高清观看| 成年免费大片在线观看| 国产一区二区三区视频了| 免费大片18禁| 国产精品一区二区免费欧美| www.熟女人妻精品国产| 此物有八面人人有两片| 日韩av在线大香蕉| 亚洲av片天天在线观看| 岛国视频午夜一区免费看| 久久这里只有精品19| 亚洲av五月六月丁香网| 亚洲av电影在线进入| 国产精品久久电影中文字幕| 99国产极品粉嫩在线观看| 欧美在线一区亚洲| 亚洲av成人一区二区三| 88av欧美| 亚洲精品456在线播放app | 免费在线观看成人毛片| 成人特级黄色片久久久久久久| 亚洲性夜色夜夜综合| 99久久成人亚洲精品观看| 日本免费一区二区三区高清不卡| 在线观看午夜福利视频| 黑人操中国人逼视频| 久久精品国产亚洲av香蕉五月| 亚洲人与动物交配视频| 亚洲熟女毛片儿| 国内毛片毛片毛片毛片毛片| 精品一区二区三区av网在线观看| 国产又色又爽无遮挡免费看| 在线看三级毛片| 亚洲人成网站在线播放欧美日韩| 欧美成人性av电影在线观看| 天堂av国产一区二区熟女人妻| 国产不卡一卡二| 别揉我奶头~嗯~啊~动态视频| 国产精品一区二区三区四区免费观看 | 久久99热这里只有精品18| 久久久国产成人免费| 日韩国内少妇激情av| 波多野结衣高清作品| 欧美一级毛片孕妇| 男人舔女人的私密视频| 这个男人来自地球电影免费观看| 一个人看视频在线观看www免费 | 久久久久性生活片| 女同久久另类99精品国产91| 国产成人精品无人区| 免费观看精品视频网站| 欧美午夜高清在线| 亚洲avbb在线观看| 日韩欧美精品v在线| 国产精华一区二区三区| 网址你懂的国产日韩在线| 成人特级av手机在线观看| 黄色丝袜av网址大全| 久久人人精品亚洲av| 青草久久国产| 国内精品一区二区在线观看| 亚洲狠狠婷婷综合久久图片| 99精品在免费线老司机午夜| 啦啦啦免费观看视频1| 男女之事视频高清在线观看| 国产伦精品一区二区三区四那| 99国产极品粉嫩在线观看| 精品久久久久久久毛片微露脸| 制服人妻中文乱码| 久久人人精品亚洲av| 两性夫妻黄色片| 色老头精品视频在线观看| 国产黄a三级三级三级人| 国产精品 欧美亚洲| 国产av在哪里看| 99riav亚洲国产免费| 无限看片的www在线观看| 中文字幕人妻丝袜一区二区| av黄色大香蕉| 五月伊人婷婷丁香| 免费在线观看成人毛片| 大型黄色视频在线免费观看| 精品欧美国产一区二区三| 亚洲色图av天堂| 露出奶头的视频| 亚洲av成人精品一区久久| 国产一区二区在线观看日韩 | 久久久久国产一级毛片高清牌| 国产爱豆传媒在线观看| 久久午夜亚洲精品久久| 日韩中文字幕欧美一区二区| 色播亚洲综合网| 中文资源天堂在线| 一个人免费在线观看的高清视频| 成人一区二区视频在线观看| 美女黄网站色视频| 久久伊人香网站| 91字幕亚洲| 免费人成视频x8x8入口观看| 日韩国内少妇激情av| 亚洲中文av在线| av黄色大香蕉| 啦啦啦免费观看视频1| 别揉我奶头~嗯~啊~动态视频| 一个人看视频在线观看www免费 | 两个人视频免费观看高清| 精品国产乱子伦一区二区三区| 美女大奶头视频| 99在线视频只有这里精品首页| 国内精品一区二区在线观看| 听说在线观看完整版免费高清| 中文字幕最新亚洲高清| 美女午夜性视频免费| 日韩有码中文字幕| 在线a可以看的网站| av在线天堂中文字幕| 久久午夜综合久久蜜桃| 午夜a级毛片| 久久这里只有精品19| 在线观看免费视频日本深夜| 黄色视频,在线免费观看| 两个人看的免费小视频| 久久久成人免费电影| 国产97色在线日韩免费| 亚洲成人久久性| 日本 欧美在线| 成人精品一区二区免费| 国产精品久久久人人做人人爽| 国产视频内射| 欧美三级亚洲精品| 久久天躁狠狠躁夜夜2o2o| 在线观看舔阴道视频| 色吧在线观看| 亚洲一区高清亚洲精品| 日韩免费av在线播放| 成人无遮挡网站| 国产熟女xx| 亚洲熟妇中文字幕五十中出| 国产综合懂色| 一区福利在线观看| 老司机午夜十八禁免费视频| 身体一侧抽搐| 18美女黄网站色大片免费观看| 国产精品亚洲美女久久久| 亚洲av成人av| 久久久久精品国产欧美久久久| 久久久国产欧美日韩av| 在线观看一区二区三区| 香蕉国产在线看| 国产伦精品一区二区三区四那| 网址你懂的国产日韩在线| 男女那种视频在线观看| 搞女人的毛片| 色在线成人网| 色播亚洲综合网| 亚洲av电影不卡..在线观看| 精品福利观看| 黄色女人牲交| 久99久视频精品免费| 中文在线观看免费www的网站| 久久久久免费精品人妻一区二区| 亚洲国产精品sss在线观看| 老鸭窝网址在线观看| 人妻丰满熟妇av一区二区三区| 男女午夜视频在线观看| 毛片女人毛片| 在线a可以看的网站| 在线观看一区二区三区| 丝袜人妻中文字幕| 在线国产一区二区在线| 亚洲成人中文字幕在线播放| 男女视频在线观看网站免费| 真实男女啪啪啪动态图| 男人舔女人的私密视频| 亚洲国产欧美人成| e午夜精品久久久久久久| 亚洲欧美精品综合一区二区三区| 国产午夜精品久久久久久| 精品久久蜜臀av无| 亚洲精品在线美女| 日韩欧美免费精品| 99热这里只有是精品50| 久久久国产欧美日韩av| 禁无遮挡网站| 一个人看视频在线观看www免费 | 91av网一区二区| 国产欧美日韩一区二区精品| 免费看光身美女| 十八禁网站免费在线| 欧美3d第一页| 欧美av亚洲av综合av国产av| 色综合欧美亚洲国产小说| av国产免费在线观看| 久久久久国产一级毛片高清牌| 美女大奶头视频| 久久久久国产一级毛片高清牌| 丁香欧美五月| 成熟少妇高潮喷水视频| 久久精品夜夜夜夜夜久久蜜豆| 人妻久久中文字幕网| 亚洲国产看品久久| 成熟少妇高潮喷水视频| 在线永久观看黄色视频| 亚洲av成人精品一区久久| 深夜精品福利| 窝窝影院91人妻| 国产亚洲欧美在线一区二区| 国产亚洲精品av在线| 俄罗斯特黄特色一大片| 两性午夜刺激爽爽歪歪视频在线观看| 日韩av在线大香蕉| 熟女电影av网| 麻豆国产av国片精品| ponron亚洲| 一级作爱视频免费观看| 神马国产精品三级电影在线观看| 又爽又黄无遮挡网站| 亚洲国产精品999在线| 国产v大片淫在线免费观看| 99re在线观看精品视频| 亚洲精品乱码久久久v下载方式 | 国产亚洲av嫩草精品影院| 小说图片视频综合网站| 岛国在线免费视频观看| 国产熟女xx| 国产精品av久久久久免费| 国产熟女xx| 男女之事视频高清在线观看| 国产精品99久久99久久久不卡| 国产精品自产拍在线观看55亚洲| 亚洲在线自拍视频| 中文字幕人妻丝袜一区二区| 美女被艹到高潮喷水动态| 国产精品久久久人人做人人爽| 美女高潮的动态| av国产免费在线观看| 精品久久久久久成人av| 久久中文看片网| 夜夜躁狠狠躁天天躁| 国产一区二区激情短视频| 天堂av国产一区二区熟女人妻| 一个人免费在线观看的高清视频| 日本熟妇午夜| 欧美成人一区二区免费高清观看 | 日本撒尿小便嘘嘘汇集6| 午夜影院日韩av| 婷婷丁香在线五月| 午夜免费观看网址| 可以在线观看毛片的网站| 亚洲av电影在线进入| 国产精品亚洲美女久久久| 成人午夜高清在线视频| 免费在线观看视频国产中文字幕亚洲| a级毛片在线看网站| 亚洲中文字幕日韩| 一进一出抽搐动态| 丁香欧美五月| 天堂影院成人在线观看| 国产午夜精品久久久久久| 亚洲成av人片在线播放无| 性色av乱码一区二区三区2| 精品福利观看| 精品欧美国产一区二区三| 久久精品人妻少妇| 首页视频小说图片口味搜索| 色噜噜av男人的天堂激情| 天堂影院成人在线观看| 亚洲精品国产精品久久久不卡| netflix在线观看网站| 亚洲精品中文字幕一二三四区| 免费av不卡在线播放| 亚洲成av人片免费观看| 久久精品人妻少妇| 色综合欧美亚洲国产小说| 国产一区在线观看成人免费| 国产精品永久免费网站| 久久精品影院6| 黑人操中国人逼视频| 两个人视频免费观看高清| 国产精品98久久久久久宅男小说| 给我免费播放毛片高清在线观看| 好男人电影高清在线观看| 成人鲁丝片一二三区免费| 亚洲国产看品久久| 精品一区二区三区av网在线观看| 久久久精品欧美日韩精品| 欧美精品啪啪一区二区三区| 黄片大片在线免费观看| 每晚都被弄得嗷嗷叫到高潮| 最新在线观看一区二区三区| 国产精品av视频在线免费观看| 国产成人av激情在线播放| 亚洲国产欧美人成| 国内精品一区二区在线观看| 精品久久久久久久末码| ponron亚洲| 国产精品亚洲av一区麻豆| 国产主播在线观看一区二区| 男女那种视频在线观看| 黄色视频,在线免费观看| 亚洲av五月六月丁香网| 最近在线观看免费完整版| 人人妻人人澡欧美一区二区| 丰满人妻熟妇乱又伦精品不卡| 天堂√8在线中文| 国产亚洲精品一区二区www| a级毛片在线看网站| 女同久久另类99精品国产91| 可以在线观看的亚洲视频| 美女大奶头视频| 亚洲天堂国产精品一区在线| 久久久久久大精品| 国产亚洲精品av在线| 欧美性猛交╳xxx乱大交人| 国产激情偷乱视频一区二区| 亚洲熟妇中文字幕五十中出| www.熟女人妻精品国产| 欧美av亚洲av综合av国产av| 男人和女人高潮做爰伦理| 国产午夜福利久久久久久| 免费看十八禁软件| 老司机在亚洲福利影院| 亚洲aⅴ乱码一区二区在线播放| 亚洲专区字幕在线| 成在线人永久免费视频| 一级毛片高清免费大全| 波多野结衣高清作品| 久久午夜综合久久蜜桃| 午夜免费观看网址| 国产高清videossex| 99久久精品热视频| 亚洲国产精品久久男人天堂| 五月伊人婷婷丁香| 91麻豆av在线| 美女免费视频网站| 欧美一级毛片孕妇| 日韩成人在线观看一区二区三区| 色视频www国产| 人人妻人人澡欧美一区二区| 久久久久性生活片| 日日夜夜操网爽| 国产亚洲欧美在线一区二区| 国产精品一区二区免费欧美| 国产高清视频在线观看网站| 看黄色毛片网站| 天天躁日日操中文字幕| 午夜视频精品福利| 欧美另类亚洲清纯唯美| 国产一区二区在线av高清观看| 国内精品一区二区在线观看| 99在线视频只有这里精品首页| 亚洲国产精品sss在线观看| 日韩欧美在线二视频| 久久久久久国产a免费观看| 激情在线观看视频在线高清| 久99久视频精品免费| 在线观看一区二区三区| 国产伦精品一区二区三区四那| 欧美在线一区亚洲| 国内少妇人妻偷人精品xxx网站 | 亚洲av免费在线观看| 国产精品亚洲美女久久久| 日本一二三区视频观看| 欧美av亚洲av综合av国产av| 午夜精品一区二区三区免费看| 搡老妇女老女人老熟妇| 1000部很黄的大片| 色av中文字幕| 国产三级黄色录像| 亚洲欧美精品综合一区二区三区| 99热6这里只有精品| 免费av毛片视频| 国产亚洲精品av在线| 久久久久九九精品影院| 免费看a级黄色片| 中文在线观看免费www的网站| 极品教师在线免费播放| 成年女人看的毛片在线观看| 99久国产av精品| 国产成人福利小说| 欧美日韩精品网址| 美女 人体艺术 gogo| 国产伦精品一区二区三区视频9 | 国产黄片美女视频| x7x7x7水蜜桃| 97超视频在线观看视频| 亚洲最大成人中文| 欧美日韩福利视频一区二区| 久久精品国产亚洲av香蕉五月| 亚洲国产欧美网| 中国美女看黄片| 2021天堂中文幕一二区在线观| 久久人妻av系列| 18禁美女被吸乳视频| 身体一侧抽搐| 精品久久久久久成人av| a级毛片a级免费在线| 亚洲成人免费电影在线观看| 免费看日本二区| 麻豆国产97在线/欧美| 婷婷精品国产亚洲av| 97超级碰碰碰精品色视频在线观看| 97超视频在线观看视频| 少妇熟女aⅴ在线视频| 女警被强在线播放| 国产免费av片在线观看野外av| 午夜日韩欧美国产| 国产精华一区二区三区| 美女午夜性视频免费| 亚洲中文日韩欧美视频| 黄片大片在线免费观看| 巨乳人妻的诱惑在线观看| 麻豆一二三区av精品| 久久精品91无色码中文字幕| av天堂在线播放| 久久热在线av| 国产成人精品久久二区二区91| 成人国产一区最新在线观看| 国产 一区 欧美 日韩| 18禁美女被吸乳视频| 国产乱人视频| 观看免费一级毛片| 国产成人精品久久二区二区91| www.自偷自拍.com| 国产精品一区二区免费欧美| 日韩av在线大香蕉| 午夜免费成人在线视频| 怎么达到女性高潮| 久久天躁狠狠躁夜夜2o2o| 无遮挡黄片免费观看| 美女被艹到高潮喷水动态| 又大又爽又粗| 久久国产乱子伦精品免费另类| 成人欧美大片| 十八禁人妻一区二区| 国产不卡一卡二| 91九色精品人成在线观看|