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

    聲波方程參數(shù)化模式及多參數(shù)全波形反演去耦合化策略

    2018-11-07 01:13:40何兵紅方伍寶胡光輝劉定進(jìn)孫思宇
    石油物探 2018年5期

    何兵紅,方伍寶,胡光輝,劉定進(jìn),孫思宇

    (1.中國石油化工股份有限公司石油物探技術(shù)研究院,江蘇南京211103;2.中國石油大學(xué)(北京),北京102249)

    隨著計(jì)算機(jī)硬件的發(fā)展和高性能計(jì)算能力的提高,基于波動理論的地震全波形反演逐步成為研究熱點(diǎn)[1]。全波形反演充分利用了地震數(shù)據(jù)的運(yùn)動學(xué)特征和動力學(xué)特征,得到較高分辨率的速度模型,從而提高地震成像精度[2-5]。隨著全波形反演研究的不斷深入,多參數(shù)全波形反演輔助儲層預(yù)測成為新的研究方向[6-9]。目前理論研究主要聚焦于彈性介質(zhì)多參數(shù)全波形反演[10-13],但是由于當(dāng)前采集得到的地震數(shù)據(jù)以縱波為主,且在應(yīng)用推廣中彈性波全波形反演受計(jì)算能力和存儲能力的限制而難以有所突破[14-16],因而有必要開展聲波方程的多參數(shù)全波形反演研究。

    多參數(shù)間的耦合是聲波方程多參數(shù)全波形反演面臨的最主要問題之一。JANNANE等[17]討論了L2范數(shù)下目標(biāo)函數(shù)與長波長、中波長和短波長擾動之間的關(guān)系。但是該研究存在一定的局限,主要問題是最大偏移距只有1760m,限制了長波長的反演。ZHOU[18]等在速度-阻抗參數(shù)化模式下先固定阻抗值,利用潛波建立低波數(shù)速度場,再固定速度值,利用反射波進(jìn)行高波數(shù)阻抗反演。董良國等[19]則通過變密度聲波方程討論了不同地震數(shù)據(jù)子集對應(yīng)的目標(biāo)函數(shù)隨物性變化的非線性特征,提出了分步驟、分尺度的多參數(shù)全波形反演方案。張廣智等[20-21]依靠初始模型對速度-密度進(jìn)行順序反演。YANG等[22-23]研究了聲波方程中密度-速度參數(shù)化模式下的參數(shù)耦合現(xiàn)象,并給出了分偏移距的反演策略。石玉梅等[24]采用頻率域單程聲波延拓法進(jìn)行波場模擬并反演得到地層密度和體積模量,提高了氣藏成像的準(zhǔn)確性。

    因多參數(shù)全波形反演中目標(biāo)函數(shù)的非線性源自波動方程的非線性,故需對波場產(chǎn)生機(jī)制進(jìn)行理論分析,進(jìn)而為多參數(shù)全波形反演提供理論指導(dǎo)。本文首先利用格林函數(shù)積分求解波動方程,在線性近似條件下推導(dǎo)了地震波振幅在6種參數(shù)化模式下(模量-密度、速度-密度、阻抗-密度、阻抗-速度、模量-速度、模量-阻抗)的敏感核函數(shù);然后研究了不同參數(shù)化模式下參數(shù)擾動的輻射模式;再從理論上分析了波場隨參數(shù)擾動的角度變化特征,對比了6種參數(shù)化下模式參數(shù)耦合特性,并提出了針對不同參數(shù)化模式的多參數(shù)全波形反演弱耦合化策略;最后進(jìn)行了數(shù)值實(shí)驗(yàn)分析。

    1 聲波方程多參數(shù)耦合模式

    1.1 速度-密度參數(shù)化模式

    地震正演是地震反演的基礎(chǔ)。聲波方程正演通常不考慮介質(zhì)密度的變化,只采用速度來描述地下介質(zhì)的物性。但實(shí)際地層中,密度也是重要的物性參數(shù),同時介質(zhì)的變化通常伴隨著密度的變化。考慮時間域變密度聲波方程為:

    (1)

    式中:v(x)為速度;ρ(x)為密度;P(x,t)為壓力場;f(x,t)為時間域震源;xs為震源位置坐標(biāo),δ(x-xs)表示將震源置于坐標(biāo)點(diǎn)xs處。

    根據(jù)介質(zhì)的擾動理論,總介質(zhì)(m)可以分為背景介質(zhì)(m0)和擾動介質(zhì)(δm),其中,m=m0+δm,對應(yīng)的總波場(P)可以分解為背景波場(P0)和擾動波場(δP),其中,P=P0+δP。在背景介質(zhì)中產(chǎn)生的背景波場滿足波動方程:

    (2)

    式中:ρ0(x),v0(x)分別為背景密度和背景速度。其中,

    (3)

    擾動波場滿足方程:

    (4)

    式中:ζ(x,t)為虛震源;P(x,t)為壓力場,在Born近似條件下,背景波場遠(yuǎn)遠(yuǎn)大于擾動波場,虛震源項(xiàng)可以近似為:

    (5)

    式中:P0(x,ω)為背景波場。擾動波場對應(yīng)的積分解可以用對應(yīng)方程的格林函數(shù)G0(r,t;x)表示:

    (6)

    式中:P0(x,t)為背景波場。

    在頻率域,卷積項(xiàng)對應(yīng)于乘積項(xiàng),各個參數(shù)統(tǒng)一用模型的相對變化量δm/m表示,(6)式可以表示為:

    (7)

    整理(7)式可以得到速度相對擾動量和密度相對擾動量引起的波場擾動:

    (8)

    其中,Kvρ-v(x)和Kvρ-ρ(x)分別為速度項(xiàng)和密度項(xiàng)在速度-密度參數(shù)化模式下的敏感核函數(shù):

    速度敏感核函數(shù)Kvρ-v(x)表示地球介質(zhì)用速度和密度表示時速度的相對變化量產(chǎn)生的擾動波場的強(qiáng)度。同理,密度敏感核函數(shù)Kvρ-ρ(x)表示在速度-密度參數(shù)化模式下密度的相對變化量產(chǎn)生的擾動波場的強(qiáng)度。

    當(dāng)用不同的參數(shù)化模式描述地球介質(zhì)時,各個參數(shù)的敏感核函數(shù)具有不同的特征。聲波方程常用的4種參數(shù)化模式包含速度-密度、模量-密度、阻抗-密度和阻抗-速度。本文又?jǐn)U展了模量-速度、阻抗-模量兩種參數(shù)化模式。我們從最常用的速度-密度聲波方程出發(fā)推導(dǎo)了6種參數(shù)化模式的敏感核函數(shù)。

    1.2 模量-密度參數(shù)化模式

    當(dāng)采用體積模量和密度對地球介質(zhì)進(jìn)行模型參數(shù)化時,在Born近似條件下,總擾動波場是體積模量和密度相對變化量引起的擾動波場的線性疊加。

    (11)

    在弱擾動條件下,存在如下關(guān)系:

    (12)

    在速度-密度和模量-密度兩種參數(shù)化模式下總波場擾動量相等,故有:

    (13)

    因此,可以得到模量-密度參數(shù)化模式下模量和密度的敏感核函數(shù):

    1.3 其它參數(shù)化模式

    依據(jù)模量-密度參數(shù)化模式下敏感核函數(shù)推導(dǎo)流程,可以得到其它4種參數(shù)化模式的敏感核函數(shù)。

    1) 阻抗-密度參數(shù)化模式下阻抗和密度的敏感核函數(shù):

    (16)

    (17)

    2) 阻抗-速度參數(shù)化模式下阻抗和速度的敏感核函數(shù):

    (18)

    (19)

    3) 模量-速度參數(shù)化模式下模量和速度的敏感核函數(shù)為:

    (20)

    (21)

    4) 模量-阻抗參數(shù)化模式下模量和阻抗的敏感核函數(shù):

    (22)

    (23)

    本部分只給出這4種參數(shù)化模式下敏感核函數(shù)的解析表達(dá)式,敏感核函數(shù)的物理意義將在下一部分結(jié)合參數(shù)輻射模式進(jìn)行闡述。

    2 參數(shù)輻射模式及全波形反演策略

    經(jīng)過推導(dǎo)得到了不同參數(shù)化模式下擾動介質(zhì)對擾動場的貢獻(xiàn)量,同時因?yàn)榻橘|(zhì)的擾動存在一定的角度特征,所以根據(jù)不同參數(shù)化模式下的敏感核函數(shù)得到各種擾動介質(zhì)對波場貢獻(xiàn)量及角度特征,并制定了多參數(shù)全波形反演策略。

    在三維介質(zhì)中,存在關(guān)系式:

    (24)

    式中:θ為入射波和散射波之間的夾角。

    因此速度-密度參數(shù)化模式下密度的敏感核函數(shù)可以表示為:

    (25)

    結(jié)合密度的敏感核函數(shù)表達(dá)式(10),得到速度和密度的輻射模式表達(dá)形式:

    (26)

    同理,我們可以分別得到其它參數(shù)化模式下的輻射模式表達(dá)式:

    將垂直向下的方向視作入射方向,利用輻射模式的表達(dá)式得到空間各個方向的散射波能量。在各向同性介質(zhì)中,任意兩個垂直方向的平面內(nèi)擾動參數(shù)的輻射模式具有完全等價性。

    圖1中箭頭表示波的傳播方向,入射波和反射波夾角為θ,反射波與x軸的夾角為φ。為了清楚地表示各參數(shù)的輻射模式,我們將依次展示單個垂直剖面內(nèi)(φ=135°)的各個參數(shù)的輻射情況。

    圖1 地震波入射反射示意

    圖2a中速度輻射為標(biāo)準(zhǔn)的圓形,表明速度擾動引起的波場擾動在各個方向上振幅大小一樣。圖2b 中,密度擾動引起的波場擾動隨著反射角的增大而減小。在小角度范圍內(nèi)相同的速度和密度相對擾動量引起的波場擾動量相當(dāng),因此僅依靠小角度地震數(shù)據(jù)難以區(qū)分速度和密度。因大角度地震數(shù)據(jù)主要由速度擾動產(chǎn)生,故可以通過大角度地震數(shù)據(jù)實(shí)現(xiàn)速度反演。全波形反演策略為先通過大角度地震數(shù)據(jù)反演速度,再利用小角度地震數(shù)據(jù)同時反演速度和密度。

    圖3為速度-密度模式下,波動方程數(shù)值模擬計(jì)算得到的擾動波場。背景速度和密度分別為3500m/s,2300kg/m3,擾動點(diǎn)位于模型的中心,擾動速度值為背景值的80%。炮點(diǎn)置于地表中心位置,位于擾動點(diǎn)正上方(本文中出現(xiàn)的6種參數(shù)化模式下的擾動波場的擾動點(diǎn)位置與圖3完全相同,且采用了同樣的觀測方式)。由圖3可知,在速度和密度相對擾動量相同的前提下,擾動引起的小角度波場擾動量基本一致,而密度擾動引起的大角度波場擾動量逐漸變?nèi)?這與圖2中的分析結(jié)果完全一致。

    圖2 φ=135°時速度-密度輻射情況a 速度;b 密度

    圖3 速度-密度模式下擾動波場a 速度擾動引起的擾動波場;b 密度擾動引起的擾動波場

    圖4a中模量輻射為圓形,表明模量擾動對波場擾動的貢獻(xiàn)量在各個方向上相等。圖4b中密度擾動對波場擾動的貢獻(xiàn)量隨著角度的增加先減小后增加,大角度的地震數(shù)據(jù)并不能解決模量和密度之間的耦合。當(dāng)θ為90°時密度的改變量對擾動場的貢獻(xiàn)為零(圖5),可以先利用中角度地震數(shù)據(jù)先反演模量,然后利用小角度和大角度地震數(shù)據(jù)同時反演模量和密度。

    圖6為阻抗-密度輻射情況,阻抗擾動對擾動場的貢獻(xiàn)量在各個方向相等,密度擾動在θ為小角度時對擾動場的貢獻(xiàn)很小,相同的密度相對擾動量在θ為大角度時和阻抗相對擾動量對擾動波場的貢獻(xiàn)量相當(dāng)(圖7)。對于阻抗-密度模式,可以先利用小角度地震數(shù)據(jù)實(shí)現(xiàn)阻抗反演,再利用大角度地震數(shù)據(jù)進(jìn)行阻抗和密度的同時反演。

    圖8為阻抗-速度輻射情況,可以看出,阻抗擾動對小角度地震數(shù)據(jù)的擾動波場貢獻(xiàn)較大,對大角度地震數(shù)據(jù)的擾動波場貢獻(xiàn)較小。與此相反,速度擾動對大角度地震數(shù)據(jù)的擾動波場貢獻(xiàn)較大,而對小角度地震數(shù)據(jù)的擾動波場貢獻(xiàn)較小。圖9中數(shù)值模擬計(jì)算結(jié)果驗(yàn)證了圖8推導(dǎo)結(jié)果的正確性,可以首先利用大角度地震數(shù)據(jù)進(jìn)行速度反演,再利用小角度地震數(shù)據(jù)進(jìn)行阻抗反演。

    圖10為模量-速度輻射模式,模量擾動主要引起小角度的擾動波場,在該參數(shù)化模式下,速度擾動對大角度和小角度擾動波場貢獻(xiàn)較大,圖11驗(yàn)證了該結(jié)論。圖12為模量-阻抗輻射模式,模量擾動主要引起大角度的擾動波場,在這種參數(shù)化模式下,阻抗擾動對大角度和小角度擾動波場貢獻(xiàn)較大,對中角度的擾動波場貢獻(xiàn)較小。圖13進(jìn)一步驗(yàn)證了該結(jié)論。

    圖4 φ=135°時模量-密度輻射情況a 模量;b 密度

    圖5 模量-密度模式下擾動波場a 模量擾動引起的擾動波場;b 密度擾動引起的擾動波場

    圖6 φ=135°時阻抗-密度輻射情況a 阻抗;b 密度

    圖8 φ=135°時阻抗-速度輻射情況a 阻抗;b 速度

    圖9 阻抗-速度模式下擾動波場a 阻抗擾動引起的擾動波場;b 速度擾動引起的擾動波場

    圖10 φ=135°時模量-速度輻射情況a 模量;b 速度

    圖11 模量-速度模式下擾動波場a 模量擾動引起的擾動波場;b 速度擾動引起的擾動波場

    圖12 φ=135°時模量-阻抗輻射情況a 模量;b 阻抗

    圖13 模量-阻抗模式下擾動波場a 模量擾動引起的擾動波場;b 阻抗擾動引起的擾動波場

    綜合以上6種參數(shù)化模式,在理想狀態(tài)下具備所有角度的擾動波場時,阻抗-速度模式下能夠利用大角度地震數(shù)據(jù)和小角度地震數(shù)據(jù)分別進(jìn)行速度和阻抗反演,該種模式下兩種參數(shù)耦合化程度較弱,是聲波方程多參數(shù)反演中較為理想的參數(shù)化模式。阻抗-密度和速度-密度輻射圖具有類似的特征,可先利用大角度或者小角度地震數(shù)據(jù)進(jìn)行單參數(shù)反演,然后進(jìn)行兩種參數(shù)的同時反演。在模量-密度、模量-速度以及模量-阻抗3種參數(shù)化模式下進(jìn)行多參數(shù)反演則難以區(qū)分不同參數(shù)的作用,多解性問題較嚴(yán)重。

    3 數(shù)值實(shí)驗(yàn)分析

    理論上認(rèn)為大角度地震數(shù)據(jù)是指120°≤θ<180°的地震數(shù)據(jù),當(dāng)θ=180°時,該數(shù)據(jù)已經(jīng)不是反射波而是透射波。在地面地震觀測方式中,大角度總是相對的,在數(shù)值模擬計(jì)算中只能盡可能提高偏移距與目標(biāo)層深度的比值,來增加最大入射角度。

    本文對速度-密度、模量-密度、阻抗-密度和速度-阻抗4種參數(shù)化模式進(jìn)行了數(shù)值測試及分析。圖14 為真實(shí)模型和對應(yīng)的初始模型,初始模型由將平滑算法作用于真實(shí)模型得到。模型設(shè)計(jì)時,保證了兩種參數(shù)化模式下波動方程正演所得到的模擬數(shù)據(jù)一致。圖15為單炮觀測數(shù)據(jù),我們首先將全波場數(shù)據(jù)作為全波形反演的觀測數(shù)據(jù)進(jìn)行速度-密度同時反演。在圖16速度-密度反演中,速度和密度真實(shí)異常得到了較好的反演,但是在速度中出現(xiàn)了密度異常引起的速度異常,并且密度項(xiàng)中速度異常引起的變化與密度本身的異常大小接近,即速度-密度同時反演時出現(xiàn)強(qiáng)烈的耦合現(xiàn)象。

    本文還對模量-密度以及阻抗-密度參數(shù)化模式進(jìn)行了數(shù)值測試。在模量-密度反演中模量的高低波數(shù)信息反演結(jié)果較為理想,但是密度的反演結(jié)果穩(wěn)定性較差(圖17)。

    圖14 真實(shí)模型及初始模型對比a 真實(shí)速度模型;b 初始速度模型;c 真實(shí)密度模型;d 初始密度模型

    圖15 單炮觀測數(shù)據(jù)

    在阻抗與密度的反演中,阻抗以小尺度信息為主,密度項(xiàng)中的部分低波數(shù)信息得到了恢復(fù)(圖18)。但總體上這兩種參數(shù)化反演中參數(shù)的耦合化現(xiàn)象比速度-密度參數(shù)化模式下的反演結(jié)果更為嚴(yán)重。

    相比之下,速度-阻抗反演模式下,速度和阻抗表現(xiàn)出非常弱的耦合。同時,也不難發(fā)現(xiàn)兩組反演中的速度表現(xiàn)出不同的特征。在速度-密度參數(shù)化模式下得到的速度異常具有清晰的邊界,而速度-阻抗參數(shù)化模式下速度異常的中、低波數(shù)較為豐富,較高波數(shù)成分缺失,邊界模糊(圖19),該現(xiàn)象與前面的理論推導(dǎo)吻合。

    圖16 全波形反演速度(a)全波形反演密度(b)

    圖17 真實(shí)模量(a)、初始模量(b)與利用全波場觀測數(shù)據(jù)同時反演模量(c)和密度(d)的結(jié)果

    圖18 真實(shí)阻抗(a)、初始阻抗(b)與利用全波場觀測數(shù)據(jù)同時反演阻抗(c)和密度(d)的結(jié)果

    速度-阻抗參數(shù)化模式下,速度異常主要與大角度地震數(shù)據(jù)有關(guān),而阻抗主要與小角度地震數(shù)據(jù)有關(guān)。為了進(jìn)一步降低速度和阻抗之間的耦合特征,本文采用中角度和大角度地震數(shù)據(jù)作為全波形反演的觀測數(shù)據(jù)(圖20)進(jìn)行速度-阻抗的同時反演。反演結(jié)果如圖21所示,可以看出速度中的高波數(shù)成分進(jìn)一步減少,中低波數(shù)成分更加突出。同時阻抗對速度的干擾降低,即速度和阻抗之間的耦合現(xiàn)象在速度反演中得到了緩解。由于缺失了小角度地震數(shù)據(jù),相對于圖19d中的全波場觀測數(shù)據(jù)的阻抗反演,僅采用中、大角度地震數(shù)據(jù)進(jìn)行阻抗反演不但阻抗異常未能準(zhǔn)確的反演,而且速度對阻抗的影響變大。因而在進(jìn)行分角度反演時,需對反演結(jié)果進(jìn)行合理的取舍。

    圖19 真實(shí)阻抗(a)、初始阻抗(b)與利用全波場觀測數(shù)據(jù)同時反演速度(c)和阻抗(d)的結(jié)果

    圖20 單炮觀測數(shù)據(jù)(中角度和大角度地震數(shù)據(jù))

    圖21 利用中、大角度地震觀測數(shù)據(jù)同時反演速度(a)和阻抗(b)

    4 討論與認(rèn)識

    本文從理論上推導(dǎo)了聲波方程不同參數(shù)化模式的敏感核函數(shù),并分析了各種模式下參數(shù)的角度特征。速度-密度參數(shù)化模式下小角度地震數(shù)據(jù)的參數(shù)耦合嚴(yán)重不易區(qū)分,大角度地震數(shù)據(jù)則能夠較好的區(qū)分;模量-密度參數(shù)化模式下小角度和大角度地震數(shù)據(jù)參數(shù)耦合都很嚴(yán)重,中角度地震數(shù)據(jù)耦合現(xiàn)象有所減弱;阻抗-密度參數(shù)化模式下大角度地震數(shù)據(jù)耦合嚴(yán)重,小角度地震數(shù)據(jù)各參數(shù)容易區(qū)分;速度-阻抗參數(shù)化模式下小角度和大角度地震數(shù)據(jù)參數(shù)耦合化現(xiàn)象都相對較弱。同時還分析了模量-速度和模量阻抗兩種參數(shù)化模式,模量-速度參數(shù)化模式下小角度地震數(shù)據(jù)參數(shù)耦合嚴(yán)重,模量-阻抗參數(shù)化模式下大角度地震數(shù)據(jù)參數(shù)耦合嚴(yán)重。綜上可知在理論上速度-阻抗參數(shù)化模式是最為理想的參數(shù)化模式。

    通過理論模型的數(shù)值模擬計(jì)算分析了常用的速度-密度和速度-阻抗參數(shù)化模式下多參數(shù)反演。本文采用的同時反演策略,順序反演會將所有的波場擾動歸于其中一種參數(shù)擾動,容易導(dǎo)致不穩(wěn)定性。數(shù)值模擬計(jì)算結(jié)果表明:速度-阻抗參數(shù)化模式下速度和阻抗的耦合化現(xiàn)象明顯減弱。相對于速度-密度參數(shù)化模式反演得到的速度,速度-阻抗參數(shù)化模式下反演速度呈現(xiàn)出低波數(shù)特性,這是由于該模式下主要是大角度地震數(shù)據(jù)作用于速度反演。

    理論推導(dǎo)在實(shí)際應(yīng)用中往往受限。例如實(shí)際采集中最大角度有限,特別是陸上資料采集偏移距受限程度較大。此外,地震數(shù)據(jù)的大角度和小角度在理論上也沒有定量的界限,因此需要根據(jù)觀測系統(tǒng)制定合理的反演策略。

    窝窝影院91人妻| 美女国产高潮福利片在线看| www日本在线高清视频| 国产伦一二天堂av在线观看| 精品一区二区三区av网在线观看| 黄色片一级片一级黄色片| 久久性视频一级片| 1024视频免费在线观看| 老鸭窝网址在线观看| 大型av网站在线播放| 日韩av在线大香蕉| 麻豆国产av国片精品| 久久午夜综合久久蜜桃| 亚洲中文字幕一区二区三区有码在线看 | 老司机靠b影院| 色综合亚洲欧美另类图片| 亚洲七黄色美女视频| 日本vs欧美在线观看视频| 国产激情欧美一区二区| 亚洲成人久久性| 精品乱码久久久久久99久播| 如日韩欧美国产精品一区二区三区| 国内久久婷婷六月综合欲色啪| 亚洲va日本ⅴa欧美va伊人久久| 99久久久亚洲精品蜜臀av| 国产成年人精品一区二区| 欧美成人一区二区免费高清观看 | 在线免费观看的www视频| 一个人观看的视频www高清免费观看 | 韩国av一区二区三区四区| 国产99白浆流出| 亚洲视频免费观看视频| 国产色视频综合| 日韩欧美国产在线观看| 国产亚洲精品久久久久久毛片| 国产在线观看jvid| 黄频高清免费视频| 欧美黄色淫秽网站| 日韩精品中文字幕看吧| 欧美在线黄色| 搡老妇女老女人老熟妇| 欧美中文综合在线视频| 激情在线观看视频在线高清| 色综合亚洲欧美另类图片| 一级,二级,三级黄色视频| 成人亚洲精品av一区二区| 熟妇人妻久久中文字幕3abv| 无遮挡黄片免费观看| 嫩草影院精品99| 日韩有码中文字幕| 国产亚洲精品av在线| 日韩欧美国产一区二区入口| 中文字幕另类日韩欧美亚洲嫩草| 少妇的丰满在线观看| 成人永久免费在线观看视频| 欧美黑人精品巨大| 一个人观看的视频www高清免费观看 | 一二三四在线观看免费中文在| 一a级毛片在线观看| 亚洲七黄色美女视频| 91麻豆av在线| 亚洲欧美日韩另类电影网站| 黄片小视频在线播放| 国产成人影院久久av| 高清在线国产一区| 国产精品99久久99久久久不卡| 久久狼人影院| 一本久久中文字幕| 人人澡人人妻人| 精品一区二区三区视频在线观看免费| 亚洲无线在线观看| 成人亚洲精品av一区二区| 亚洲午夜精品一区,二区,三区| 久久国产精品男人的天堂亚洲| 男女之事视频高清在线观看| 91麻豆av在线| 亚洲电影在线观看av| 国产精品免费一区二区三区在线| 91国产中文字幕| 久久人人爽av亚洲精品天堂| 自线自在国产av| 色综合站精品国产| 欧美丝袜亚洲另类 | 男女下面插进去视频免费观看| 成年女人毛片免费观看观看9| 免费在线观看视频国产中文字幕亚洲| 一区福利在线观看| 欧美成狂野欧美在线观看| 亚洲 欧美一区二区三区| 少妇粗大呻吟视频| 18禁美女被吸乳视频| 国产成人欧美在线观看| 国产乱人伦免费视频| 久久香蕉精品热| 悠悠久久av| 亚洲成人久久性| 亚洲欧美日韩无卡精品| 欧美黑人欧美精品刺激| 日韩av在线大香蕉| 美女 人体艺术 gogo| 国产成人啪精品午夜网站| 免费久久久久久久精品成人欧美视频| 亚洲第一电影网av| 97人妻天天添夜夜摸| 亚洲激情在线av| 国产精华一区二区三区| 日韩大尺度精品在线看网址 | 操出白浆在线播放| 精品一区二区三区av网在线观看| 岛国在线观看网站| 国产视频一区二区在线看| 欧美在线黄色| АⅤ资源中文在线天堂| 色综合站精品国产| 露出奶头的视频| 亚洲国产精品成人综合色| 国产成人一区二区三区免费视频网站| 久久伊人香网站| 亚洲av成人一区二区三| 一进一出抽搐gif免费好疼| 久久国产乱子伦精品免费另类| 少妇被粗大的猛进出69影院| 国产精品爽爽va在线观看网站 | 悠悠久久av| 亚洲全国av大片| 久久人人精品亚洲av| 一a级毛片在线观看| 在线永久观看黄色视频| 久久久精品欧美日韩精品| 亚洲人成伊人成综合网2020| 亚洲av五月六月丁香网| 免费av毛片视频| 欧美色视频一区免费| 丝袜美足系列| 国产高清激情床上av| 亚洲五月婷婷丁香| 亚洲 欧美一区二区三区| 在线观看免费视频日本深夜| 久久狼人影院| 久久久精品欧美日韩精品| 真人一进一出gif抽搐免费| 99久久综合精品五月天人人| 久久久久久免费高清国产稀缺| 亚洲精品美女久久av网站| 黄色片一级片一级黄色片| 少妇熟女aⅴ在线视频| 黄色女人牲交| 操美女的视频在线观看| 女同久久另类99精品国产91| 国产av精品麻豆| 在线av久久热| 色哟哟哟哟哟哟| 国产一区二区三区在线臀色熟女| 亚洲国产精品久久男人天堂| 男人舔女人的私密视频| 日日干狠狠操夜夜爽| 伦理电影免费视频| 午夜福利高清视频| 国产亚洲精品久久久久5区| 亚洲国产精品久久男人天堂| 人人妻,人人澡人人爽秒播| 人成视频在线观看免费观看| 国产精品电影一区二区三区| 国产av一区二区精品久久| 国产亚洲av嫩草精品影院| 99国产精品99久久久久| 久久亚洲精品不卡| 国内精品久久久久久久电影| 女同久久另类99精品国产91| 日韩视频一区二区在线观看| 午夜久久久久精精品| 国产爱豆传媒在线观看| 99久久无色码亚洲精品果冻| 联通29元200g的流量卡| 亚洲精品成人久久久久久| www.www免费av| 欧美高清成人免费视频www| 日韩精品有码人妻一区| 无人区码免费观看不卡| 国产熟女欧美一区二区| 偷拍熟女少妇极品色| 色尼玛亚洲综合影院| 两个人的视频大全免费| 岛国在线免费视频观看| 欧美潮喷喷水| 国产高清激情床上av| 国产淫片久久久久久久久| 色5月婷婷丁香| 精品一区二区免费观看| 在线看三级毛片| 亚洲色图av天堂| 色哟哟·www| 成人亚洲精品av一区二区| 亚洲电影在线观看av| 99九九线精品视频在线观看视频| 男女那种视频在线观看| 精品久久久久久久久亚洲 | 亚洲 国产 在线| а√天堂www在线а√下载| 久久久久久久亚洲中文字幕| avwww免费| 国产精品一区www在线观看 | 一级黄片播放器| 欧美最新免费一区二区三区| 午夜亚洲福利在线播放| 窝窝影院91人妻| 成人国产麻豆网| 成人av在线播放网站| 伊人久久精品亚洲午夜| 亚洲av成人av| 亚洲国产精品久久男人天堂| 亚洲国产精品sss在线观看| 亚洲欧美清纯卡通| 亚洲最大成人中文| 亚洲av日韩精品久久久久久密| 别揉我奶头 嗯啊视频| 久久久久久九九精品二区国产| 欧美在线一区亚洲| 在线观看舔阴道视频| 国产精品免费一区二区三区在线| 成熟少妇高潮喷水视频| 国产精品久久久久久亚洲av鲁大| 国产成人影院久久av| 欧美+日韩+精品| 亚洲av免费在线观看| 国产69精品久久久久777片| 精品不卡国产一区二区三区| 日本免费a在线| 精品一区二区三区av网在线观看| 少妇被粗大猛烈的视频| 国产69精品久久久久777片| 自拍偷自拍亚洲精品老妇| 亚洲av免费在线观看| 国产亚洲精品久久久com| 乱码一卡2卡4卡精品| 色综合色国产| 色在线成人网| 狠狠狠狠99中文字幕| 欧美色欧美亚洲另类二区| 在线国产一区二区在线| 一进一出抽搐gif免费好疼| 国产免费av片在线观看野外av| 亚洲黑人精品在线| 桃红色精品国产亚洲av| 成人欧美大片| 成人av一区二区三区在线看| 国产av不卡久久| 久久久国产成人免费| 亚洲在线自拍视频| 国产免费男女视频| 国产 一区精品| 亚洲av一区综合| 日本黄色片子视频| 日韩在线高清观看一区二区三区 | 最近中文字幕高清免费大全6 | 国产毛片a区久久久久| 国产黄a三级三级三级人| 国产精品久久电影中文字幕| 午夜免费成人在线视频| 成人性生交大片免费视频hd| 夜夜看夜夜爽夜夜摸| 岛国在线免费视频观看| 亚洲无线观看免费| 综合色av麻豆| 窝窝影院91人妻| 久久精品国产清高在天天线| 中文字幕免费在线视频6| 欧美+亚洲+日韩+国产| www.www免费av| 91狼人影院| 男女下面进入的视频免费午夜| 日韩强制内射视频| 久久久国产成人免费| 日日夜夜操网爽| 欧美不卡视频在线免费观看| 国产精品av视频在线免费观看| 狂野欧美激情性xxxx在线观看| 午夜福利视频1000在线观看| 国产av一区在线观看免费| 国产成人一区二区在线| 久久热精品热| 精品午夜福利在线看| 亚洲四区av| 久久久久久久亚洲中文字幕| 草草在线视频免费看| 国产视频内射| 日本五十路高清| 精品欧美国产一区二区三| 成年版毛片免费区| 免费人成视频x8x8入口观看| 两性午夜刺激爽爽歪歪视频在线观看| 18禁在线播放成人免费| 人人妻人人看人人澡| 亚洲国产欧洲综合997久久,| 国产在线精品亚洲第一网站| 97超视频在线观看视频| 国产综合懂色| 亚洲综合色惰| 日韩欧美在线二视频| 村上凉子中文字幕在线| 最新在线观看一区二区三区| 男女啪啪激烈高潮av片| 国产精品美女特级片免费视频播放器| av在线蜜桃| 俺也久久电影网| 久久久久国产精品人妻aⅴ院| 国产精品爽爽va在线观看网站| 免费观看精品视频网站| 美女高潮喷水抽搐中文字幕| 天堂影院成人在线观看| 日本 av在线| 国产69精品久久久久777片| 成人国产一区最新在线观看| 日韩亚洲欧美综合| 免费在线观看影片大全网站| 精品一区二区三区视频在线| 不卡一级毛片| 久久午夜亚洲精品久久| 久久精品91蜜桃| 少妇猛男粗大的猛烈进出视频 | 大又大粗又爽又黄少妇毛片口| 久久精品国产清高在天天线| 乱码一卡2卡4卡精品| 搡女人真爽免费视频火全软件 | 色尼玛亚洲综合影院| 精华霜和精华液先用哪个| 麻豆av噜噜一区二区三区| 此物有八面人人有两片| 久久久久性生活片| 日本与韩国留学比较| 国产成人aa在线观看| 日本-黄色视频高清免费观看| 在线观看午夜福利视频| 三级男女做爰猛烈吃奶摸视频| 久久久久久久精品吃奶| 精品福利观看| 男女视频在线观看网站免费| 在线观看66精品国产| 国产精品久久久久久av不卡| 免费av不卡在线播放| 少妇猛男粗大的猛烈进出视频 | 国产精品不卡视频一区二区| av在线天堂中文字幕| 18+在线观看网站| 精品久久久久久久久av| 黄色欧美视频在线观看| 欧美精品国产亚洲| 欧美一级a爱片免费观看看| 88av欧美| 99热这里只有是精品50| 亚洲精品色激情综合| 十八禁国产超污无遮挡网站| 亚洲av成人av| 国产精品野战在线观看| 狂野欧美激情性xxxx在线观看| 成人精品一区二区免费| 色综合站精品国产| 搡老熟女国产l中国老女人| 舔av片在线| 国产精品不卡视频一区二区| 少妇裸体淫交视频免费看高清| netflix在线观看网站| 国产一区二区亚洲精品在线观看| 人人妻人人看人人澡| 久久久久久大精品| 日韩在线高清观看一区二区三区 | 成人综合一区亚洲| 亚洲精品粉嫩美女一区| 国产精品人妻久久久久久| or卡值多少钱| 国产免费男女视频| 免费人成视频x8x8入口观看| 小说图片视频综合网站| 国产亚洲精品久久久com| 国产午夜福利久久久久久| 久久九九热精品免费| 国产一区二区三区视频了| 久久午夜亚洲精品久久| 中文字幕av成人在线电影| 狂野欧美白嫩少妇大欣赏| 22中文网久久字幕| 97超级碰碰碰精品色视频在线观看| 人人妻人人澡欧美一区二区| 久久久国产成人精品二区| 99国产极品粉嫩在线观看| 中文字幕精品亚洲无线码一区| 22中文网久久字幕| 桃色一区二区三区在线观看| av在线观看视频网站免费| 欧美丝袜亚洲另类 | 欧美中文日本在线观看视频| 成人av一区二区三区在线看| 97碰自拍视频| 欧美在线一区亚洲| 啦啦啦韩国在线观看视频| 国内揄拍国产精品人妻在线| 真人一进一出gif抽搐免费| 亚洲经典国产精华液单| a级毛片免费高清观看在线播放| 九九爱精品视频在线观看| 91午夜精品亚洲一区二区三区 | 日韩精品中文字幕看吧| 久久久久久久午夜电影| 又爽又黄a免费视频| 最新中文字幕久久久久| 一级黄片播放器| 亚洲中文字幕一区二区三区有码在线看| 国产精品爽爽va在线观看网站| 99热这里只有是精品50| 色噜噜av男人的天堂激情| 成人av在线播放网站| 特级一级黄色大片| 亚洲内射少妇av| 亚洲久久久久久中文字幕| 午夜福利在线在线| 丝袜美腿在线中文| 91狼人影院| or卡值多少钱| 亚洲成av人片在线播放无| 午夜免费成人在线视频| av视频在线观看入口| 99久久精品一区二区三区| 午夜a级毛片| 午夜激情福利司机影院| av福利片在线观看| 亚洲av日韩精品久久久久久密| 十八禁国产超污无遮挡网站| 美女大奶头视频| 国产v大片淫在线免费观看| 日本 av在线| 久久国产精品人妻蜜桃| 欧美一区二区国产精品久久精品| 熟妇人妻久久中文字幕3abv| 国产精品久久视频播放| 悠悠久久av| 韩国av在线不卡| 亚洲成人中文字幕在线播放| 中文字幕人妻熟人妻熟丝袜美| 亚洲一级一片aⅴ在线观看| 免费人成在线观看视频色| 中文字幕熟女人妻在线| 国产精品一区二区三区四区免费观看 | 久久久成人免费电影| 欧美色视频一区免费| 天天一区二区日本电影三级| 精品午夜福利在线看| 成熟少妇高潮喷水视频| 国产高潮美女av| 一个人观看的视频www高清免费观看| 小蜜桃在线观看免费完整版高清| 久久久成人免费电影| 久久草成人影院| 久久精品人妻少妇| 毛片女人毛片| 国语自产精品视频在线第100页| 我要搜黄色片| 国产精品人妻久久久久久| 日本黄大片高清| 国产免费一级a男人的天堂| 我要搜黄色片| 中文字幕av在线有码专区| 久久精品国产亚洲av香蕉五月| 欧美激情久久久久久爽电影| 黄片wwwwww| 成人鲁丝片一二三区免费| 22中文网久久字幕| 日韩亚洲欧美综合| 九九在线视频观看精品| 色精品久久人妻99蜜桃| 欧美日韩瑟瑟在线播放| 日本爱情动作片www.在线观看 | 日韩大尺度精品在线看网址| 91在线观看av| 日韩精品有码人妻一区| 午夜激情欧美在线| 亚洲成人久久性| 简卡轻食公司| 在线观看av片永久免费下载| 国产精品嫩草影院av在线观看 | 九色国产91popny在线| 香蕉av资源在线| 久久久久精品国产欧美久久久| 婷婷六月久久综合丁香| 国产日本99.免费观看| 国产一级毛片七仙女欲春2| 午夜福利在线观看吧| 国产探花极品一区二区| 男女边吃奶边做爰视频| 久久久精品大字幕| 三级男女做爰猛烈吃奶摸视频| 热99在线观看视频| 91麻豆av在线| 日日啪夜夜撸| 俺也久久电影网| 日韩中文字幕欧美一区二区| 欧美成人a在线观看| 欧美+日韩+精品| 免费大片18禁| 色5月婷婷丁香| 午夜亚洲福利在线播放| 欧美绝顶高潮抽搐喷水| 俄罗斯特黄特色一大片| 色综合婷婷激情| av在线天堂中文字幕| 天天一区二区日本电影三级| 亚洲国产精品成人综合色| 欧美日韩瑟瑟在线播放| 久久午夜福利片| 成人毛片a级毛片在线播放| 中出人妻视频一区二区| 波多野结衣高清作品| 中国美白少妇内射xxxbb| 欧美成人一区二区免费高清观看| 男人舔奶头视频| 久久精品国产亚洲av天美| 久久精品影院6| 日本熟妇午夜| 久久九九热精品免费| 无遮挡黄片免费观看| 又黄又爽又刺激的免费视频.| 国产一区二区三区在线臀色熟女| 麻豆成人av在线观看| 亚洲人成伊人成综合网2020| 久久久久久久久大av| 日韩在线高清观看一区二区三区 | 十八禁国产超污无遮挡网站| 日韩欧美免费精品| 国产极品精品免费视频能看的| 成人鲁丝片一二三区免费| 在线免费观看的www视频| 亚洲精华国产精华液的使用体验 | 亚洲七黄色美女视频| 三级国产精品欧美在线观看| 成熟少妇高潮喷水视频| 男女边吃奶边做爰视频| .国产精品久久| 很黄的视频免费| 九九久久精品国产亚洲av麻豆| 美女被艹到高潮喷水动态| 日韩av在线大香蕉| 91午夜精品亚洲一区二区三区 | 琪琪午夜伦伦电影理论片6080| 日韩精品中文字幕看吧| 国产美女午夜福利| 国产免费av片在线观看野外av| 日韩欧美在线乱码| 在线观看午夜福利视频| 午夜福利18| 久久人妻av系列| 99九九线精品视频在线观看视频| av.在线天堂| 九九爱精品视频在线观看| avwww免费| 国产精品永久免费网站| 亚洲综合色惰| 天堂网av新在线| 天堂√8在线中文| 少妇的逼水好多| 制服丝袜大香蕉在线| 黄色配什么色好看| 久久久久九九精品影院| 人妻夜夜爽99麻豆av| 午夜激情福利司机影院| 一级a爱片免费观看的视频| 在线播放国产精品三级| 精品久久久久久久久av| 国产伦在线观看视频一区| 国产精品久久久久久av不卡| 欧美日韩黄片免| 久99久视频精品免费| 国产真实伦视频高清在线观看 | 高清日韩中文字幕在线| 真实男女啪啪啪动态图| 99在线人妻在线中文字幕| 国产亚洲精品综合一区在线观看| av视频在线观看入口| 免费观看在线日韩| 国产麻豆成人av免费视频| 日本与韩国留学比较| 九九爱精品视频在线观看| 亚洲av成人av| 国产伦在线观看视频一区| 欧美黑人巨大hd| 亚洲一区二区三区色噜噜| aaaaa片日本免费| 搡女人真爽免费视频火全软件 | 国产精品久久电影中文字幕| 88av欧美| 国产午夜精品论理片| 亚洲aⅴ乱码一区二区在线播放| 99视频精品全部免费 在线| 欧美色视频一区免费| 午夜免费成人在线视频| 特大巨黑吊av在线直播| 国产蜜桃级精品一区二区三区| 亚洲av成人精品一区久久| or卡值多少钱| 又黄又爽又免费观看的视频| 亚洲最大成人手机在线| 亚洲精品一区av在线观看| 又爽又黄a免费视频| 校园人妻丝袜中文字幕| 大型黄色视频在线免费观看| 国产av一区在线观看免费| 搡老妇女老女人老熟妇| 成人国产一区最新在线观看| 男女那种视频在线观看| 亚洲黑人精品在线| 欧美日本亚洲视频在线播放| 国产精品电影一区二区三区| 亚洲精品久久国产高清桃花| 婷婷色综合大香蕉| netflix在线观看网站| 国产淫片久久久久久久久| 久久精品国产清高在天天线|