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

    基于反演策略的數(shù)據(jù)自適應衰減補償方法

    2014-03-25 06:51:20王本鋒陳小宏李景葉陳增保
    石油物探 2014年3期
    關(guān)鍵詞:子波雙曲位子

    王本鋒,陳小宏,李景葉,陳增保

    (1.中國石油大學(北京)油氣資源與探測國家重點實驗室,北京102249;2.中國石油大學(北京)海洋石油勘探國家工程實驗室,北京102249)

    由于地下介質(zhì)的非均勻性以及黏滯性對地震波產(chǎn)生了吸收衰減作用(即正Q濾波效應),使其振幅減弱、相位畸變,地震記錄的分辨率降低[1]。Varela等[2]利用Futterman衰減頻散關(guān)系分析了衰減與頻散效應,驗證了地下介質(zhì)的吸收衰減效應。衰減補償是提高地震記錄分辨率,進行精細油藏描述的有效途徑之一[3]。反Q濾波是衰減補償?shù)闹匾椒?其相位校正無條件穩(wěn)定且具有較強的抗噪性,但振幅補償因子隨著傳播時間的增大以及頻率的增加呈指數(shù)放大趨勢,穩(wěn)定性及抗噪性較差[4-7]。Hargreaves等[8]利用類似Stolt頻率-波數(shù)域偏移的方法,高效地實現(xiàn)了反Q濾波算法,有效地校正了由速度頻散引起的相位畸變,但為了避免振幅補償?shù)牟环€(wěn)定性,僅對相位進行了校正。趙建勛等[9]通過反Q濾波與串聯(lián)偏移相結(jié)合,將常Q算法推廣至深變Q值模型;通過串聯(lián)頻散補償算法與限幅補償算法相結(jié)合,更好地提高了地震記錄的分辨率。

    針對反Q濾波振幅補償?shù)牟环€(wěn)定性,Wang[4,6]提出了振幅增益控制方法以及穩(wěn)定的反Q濾波方法(可同時補償振幅以及校正相位),設(shè)計了Gabor域算法以提高計算效率。對含噪地震記錄,振幅增益控制方法在提高分辨率的同時,較大地放大了噪聲,導致振幅補償不穩(wěn)定;穩(wěn)定的反Q濾波方法可較好地處理含噪地震記錄,但一定程度上導致振幅欠補償。Yan等[10]將穩(wěn)定的反Q濾波公式推廣至多分量地震勘探中,利用疊前轉(zhuǎn)換波道集提取S波的Q值,利用疊前縱波道集提取P波的Q值,實現(xiàn)了疊前PP波以及PS波的衰減補償,提高了疊前地震記錄的分辨率。Zhao等[11]為克服振幅補償過程中噪聲被放大的缺陷,在時頻域定義信噪比,對信噪比大于1的區(qū)域進行振幅補償,其它區(qū)域不進行補償,從而在保證信噪比的前提下,提高了地震記錄的分辨率。為減弱反Q濾波對Q值的依賴性,Braga等[12]在小波域?qū)崿F(xiàn)反Q濾波方法,在Q值誤差較大時,補償后的地震記錄能在提高分辨率的同時依然保持原有的AVO特性,從而有利于疊前反演。Zhang等[13]為克服反Q濾波振幅補償不穩(wěn)定或欠補償?shù)娜毕?利用最小二乘策略以及貝葉斯理論,給定子波或提取最小相位子波來迭代反演稀疏的反射系數(shù)序列,提取Q值作相位校正,得到高分辨率的地震記錄。Wang[14]利用Futterman吸收衰減模型,基于爆炸反射面模型推導了正Q濾波公式,利用正則化策略以及反演理論反演得到補償后的地震記錄,克服了常規(guī)反Q濾波方法振幅補償不穩(wěn)定或欠補償?shù)娜毕?且具有一定的抗噪性。Margrave等[15]將維納反褶積進行推廣,提出Gabor反褶積策略;Margrave等[16]在Q值未知以及反射系數(shù)白噪的前提下,對衰減地震記錄的Gabor譜進行雙曲平滑,統(tǒng)計衰減函數(shù)以及最小相位子波,最終得到高分辨率的反射系數(shù)序列。Reine等[17]基于窗口變化與否分析了時頻分析方法度量地震波衰減的穩(wěn)健性,提出變窗口的時頻變換方法可降低衰減估計的不確定性以及估計偏差,使其具有較高的精度,但該時頻分析方法均假設(shè)子波是最小相位的,一定程度上偏離實際情況。

    我們借鑒反Q濾波思想僅進行吸收衰減補償,在Gabor域基于能量均分雙曲平滑方法估計衰減函數(shù)譜,基于非穩(wěn)態(tài)褶積理論[18],利用正則化方法反演得到衰減補償后的高分辨率地震記錄。

    2 方法原理

    2.1 非穩(wěn)態(tài)褶積理論

    假設(shè)震源子波為w(t),反射系數(shù)序列為r(t),根據(jù)褶積原理,理想的無噪聲地震記錄s(t)的表述為

    (1)

    由于地下介質(zhì)的非均勻性以及黏滯性,地震子波在傳播過程中被吸收、衰減產(chǎn)生頻散效應,實際上為非穩(wěn)態(tài)子波,因此基于穩(wěn)態(tài)褶積的處理結(jié)果往往偏離實際情況[16]。Margrave[18]將穩(wěn)態(tài)褶積模型推廣到非穩(wěn)態(tài)褶積情況,公式為

    (2)

    根據(jù)Gabor變換的性質(zhì)知

    (3)

    其中,S(f,τ),R(f,τ),α(f,τ),W(f)分別為衰減地震記錄以及反射系數(shù)的Gabor譜、衰減函數(shù)、穩(wěn)態(tài)子波頻譜。在時頻域,可以較容易地估計出非穩(wěn)態(tài)子波的振幅譜,并在最小相位的假設(shè)下,借助希爾伯特變換求得相位譜,從而實現(xiàn)非穩(wěn)態(tài)反褶積。因此將非穩(wěn)態(tài)褶積公式(2)轉(zhuǎn)換到時頻域:

    (4)

    其中,S(f)為衰減記錄的頻譜,αw(f,τ)為衰減子波的時頻譜。

    由于實際子波大多為混合相位子波,基于最小相位子波假設(shè)的非穩(wěn)態(tài)反褶積方法一定程度上偏離實際情況,因此我們忽略子波的影響,僅消除吸收衰減的效應,提高地震記錄的分辨率。非穩(wěn)態(tài)褶積公式(4)演變?yōu)?/p>

    (5)

    其中,α(f,τ)為復值衰減函數(shù),s0(τ)是原始未衰減的地震記錄。

    公式(5)可離散成矩陣方程組的形式,即

    (6)

    其中,S為包含地震數(shù)據(jù)帶限頻率成分的復值向量,Φ為含有衰減以及傅里葉變換信息的復值矩陣,s0為原始未衰減記錄(補償后記錄),為一實值向量。

    將公式(6)轉(zhuǎn)換為實值方程,即

    (7)

    (8)

    其中,μ為正則化因子,D為微分算子,可選為單位算子、一階微分算子或二階微分算子。

    將(8)式分別對s0求導并令之為0,即可獲得關(guān)于s0的正則化方程:

    (9)

    可利用共軛梯度算法、牛頓法等對正則化方程進行求解,最終得到衰減補償后的地震記錄。但是復值衰減函數(shù)α(f,τ)一般是未知的,如何得到精確的衰減函數(shù),成為衰減補償?shù)年P(guān)鍵問題。

    2.2 衰減函數(shù)估計

    衰減函數(shù)的理論表達式為

    (10)

    其中,Q為品質(zhì)因子,H(·)為沿頻率方向的希爾伯特變換,其假設(shè)衰減函數(shù)為最小相位。如果已知Q值,則衰減函數(shù)具有解析表達式。但是精確的Q值估計相對困難且其誤差會影響衰減補償?shù)木?導致錯誤的解釋結(jié)果。因此,避開Q值估計問題,直接估計衰減函數(shù)可提高補償精度。令c=τf,根據(jù)衰減函數(shù)的理論表達式知,在c等值線上衰減值相同,因此在反射系數(shù)白噪的假設(shè)前提下,對衰減地震記錄的Gabor譜進行雙曲平滑可估計衰減函數(shù)的振幅譜;在最小相位假設(shè)下,利用希爾伯特變換可求得相應的相位譜,利用公式(9)反演得到衰減補償后的地震記錄。

    不同的雙曲平滑方法對衰減函數(shù)的估計精度不同。傳統(tǒng)的雙曲平滑方法[16]將時頻域分成N個雙曲型條帶,條帶邊界為ci=(i-1)dc,i=2,…,N+1,c1=0,dc=cmax/N,如圖1所示。

    圖1 傳統(tǒng)雙曲平滑方法條帶分割a N=10; b N=30

    由圖1a可知,傳統(tǒng)的雙曲平滑分割方法對c值較小的區(qū)域欠采樣,對c值較大的區(qū)域過采樣,而地震信號的主要能量集中在c值較小的區(qū)域,使得估計的衰減函數(shù)譜精度較低。增加劃分條帶的數(shù)目,如圖1b所示,一定程度上可提高估計精度,但是光滑方法屬于統(tǒng)計方法,條帶數(shù)目增多,則每一條帶內(nèi)的樣點數(shù)減少,統(tǒng)計效應減弱。

    另外,以ci=exp[(i-1)dc],i=2,…,N+1,c1=0為條帶邊界,也可將時頻域分成N個雙曲型條帶,如圖2所示,其中dc=(lncmax)/N。得到的雙曲型條帶寬度從低c值到高c值逐漸增加,克服了傳統(tǒng)雙曲光滑方法在c值較小區(qū)域的欠采樣問題。c值較大區(qū)域的過采樣問題也得到了一定程度上的緩解,但是位于低c值條帶內(nèi)的樣點數(shù)減少,統(tǒng)計效應減弱,且其只是運用數(shù)學方法機械地將條帶寬度進行調(diào)整,缺乏物理涵義。

    圖2 對數(shù)域雙曲平滑方法條帶分割a N=10; b N=30

    針對雙曲平滑方法遇到的問題,Li等[19]基于能量均分思想,提出變步長的雙曲平滑方法。由理論衰減函數(shù)振幅譜A=exp(-πτf/Q)=exp(-πc/Q)知,c在[0,cmax]變化時,A的取值范圍為(0,1],將能量區(qū)間進行均分,即dA=1/N,得到條帶邊界ci=-{ln[(i-1)dA]}Q/π,c1=0。品質(zhì)因子Q未知,但由于平滑得到的估計值對Q值不敏感,所以Q取經(jīng)驗常數(shù)值即可。能量均分雙曲平滑分割方法具有一定的物理涵義,其按照能量進行區(qū)域分割與疊前反演抽取近、中、遠偏移距道集類似。圖3展示了分割份數(shù)N=10,30時分別對應的c條帶分割結(jié)果。分割份數(shù)較少時,如圖3a所示,在保持統(tǒng)計效應的同時,克服了c值較小的區(qū)域欠采樣以及c值較大的區(qū)域過采樣問題。

    圖3 能量平滑雙曲平滑方法條帶分割a N=10; b N=30

    將時頻域分割成N個雙曲型條帶之后,在每一條帶內(nèi)估計均值Ek,并將該平均值放置在雙曲條帶的中心線mk上:Ek=mean({|S(τ,f)|,(τ,f)∈Ωk}),其中Ωk為第k個雙曲型條帶,mk=0.5×(ck+ck+1)。

    假設(shè)時頻域(τ,f)點位于[mk,mk+1],其中mk為第k個條帶的中心,由雙曲平滑方法得其時頻譜值Ek,則(τ,f)點處的時頻譜值可由插值方法獲得,即

    (11)

    3 數(shù)值模型算例與實際資料處理

    本文方法為單道處理方法,為驗證方法的可行性,首先設(shè)計反射系數(shù)序列,分別與最小相位以及零相位子波褶積合成無衰減的地震記錄,利用非穩(wěn)態(tài)褶積理論合成衰減的地震記錄;然后基于傳統(tǒng)雙曲平滑以及能量均分雙曲平滑方法進行衰減函數(shù)估計,利用反演策略進行衰減補償。通過對實際資料進行處理,驗證了方法的有效性。

    3.1 理論模型分析

    設(shè)計的反射系數(shù)序列如圖4a所示;零相位子波為30Hz的雷克子波,如圖4b所示;最小相位子波為與零相位子波振幅譜相對應的最小相位子波,如圖4c所示。

    基于褶積理論以及非穩(wěn)態(tài)褶積理論(品質(zhì)因子Q=50)的最小相位子波合成記錄、零相位子波合成記錄及其理論衰減函數(shù)譜分別如圖5,圖6和圖7 所示。由圖5b,圖6b可知,隨著傳播時間的增加,子波的能量逐漸減弱,且高頻損失相對嚴重;由圖7可知,隨著時間的增加和頻率的提高,衰減量越來越大。

    圖4 反射系數(shù)序列(a)與零相位子波(b)以及最小相位子波(c)

    對最小相位子波合成的未衰減以及衰減地震記錄(圖5)進行Gabor譜分析,如圖8所示。由圖8a 可知,未衰減地震記錄Gabor譜的深、淺層能量相對均衡,由于帶限子波的影響,能量主要集中在中、低頻區(qū)域;由圖8b可知,衰減地震記錄的Gabor譜能量隨著傳播時間的增加逐漸減小,且高頻成分衰減相對較快,與理論衰減規(guī)律一致。

    圖5 衰減前(a)、后(b)的最小相位子波合成記錄

    圖6 衰減前(a)、后(b)的零相位子波合成記錄

    圖7 理論衰減函數(shù)譜

    圖8 最小相位子波合成的未衰減記錄Gabor譜(a)以及衰減記錄Gabor譜(b)

    基于傳統(tǒng)的雙曲平滑方法,將時頻區(qū)域分割成N=10,30,50,100份得到的衰減函數(shù)譜如圖9a至圖9d所示。可以看出,隨著分割份數(shù)的增加,估計得到的衰減譜與理論衰減譜(圖7)越來越接近;但是條帶數(shù)越多,每個條帶內(nèi)的采樣點越少,統(tǒng)計效應越弱?;谀芰烤值碾p曲平滑方法(N=10)估計得到的衰減譜如圖10所示,該方法在分割份數(shù)較少的情況下與理論衰減模型(圖7)具有較好的一致性,在保持統(tǒng)計性的同時,具有較高的精度。由估計的衰減函數(shù)振幅譜在最小相位假設(shè)下利用希爾伯特變換得到相位譜,基于公式(9)反演得到補償后的地震記錄,如圖11所示。

    圖9 傳統(tǒng)雙曲平滑方法不同分割份數(shù)估計的衰減譜a N=10; b N=30; c N=50; d N=100

    圖10 能量均分雙曲平滑方法估計的衰減譜(N=10)

    分析圖9至圖11可見,傳統(tǒng)的雙曲平滑方法估計的衰減譜的精度隨著分割份數(shù)的增加而提高,基于該衰減譜反演得到的補償記錄的精度也相應提高;基于能量均分的雙曲平滑方法在分割份數(shù)較少、保證統(tǒng)計性的前提下,估計得到的衰減譜(圖10)與理論衰減譜(圖7)具有較高的一致性,利用該衰減譜反演得到的補償后的地震記錄(圖11g)與未衰減地震記錄(圖11a)具有較高的一致性。綜合分析可知,基于能量均分的雙曲平滑方法相對傳統(tǒng)的雙曲平滑方法具有較大的優(yōu)勢,估計得到的衰減譜精度較高,反演得到的補償后地震記錄與原始未衰減記錄具有較好的一致性。

    圖11 理論模型數(shù)據(jù)衰減補償前后的地震記錄a 原始未衰減地震道; b 品質(zhì)因子Q=50時的衰減記錄; c—f 基于傳統(tǒng)雙曲平滑方法在分割份數(shù)分別為10,30,50,100時估計的衰減譜反演得到的補償后地震記錄; g 基于能量均分雙曲平滑方法估計的衰減譜反演得到的補償后地震記錄

    將能量均分雙曲平滑方法以及反演策略應用于零相位子波合成的地震記錄(圖6)中,處理過程與處理最小相位子波合成記錄類似,補償后的地震記錄如圖12所示。由圖12可見,基于能量均分雙曲平滑方法以及反演策略得到的補償?shù)卣鹩涗浥c原始未衰減地震記錄具有較好的一致性,補償?shù)木容^高。理論模型數(shù)據(jù)試算的結(jié)果表明,我們研究并提出的方法對地震子波的相位沒有特殊要求,且不需要精確Q值的信息,衰減函數(shù)從地震記錄的Gabor譜中統(tǒng)計得到,利用反演的策略得到衰減補償后的地震數(shù)據(jù)避免了振幅補償?shù)牟环€(wěn)定性,實現(xiàn)了提高地震記錄分辨率的目的。由于忽略子波的影響,該方法具有數(shù)據(jù)自適應性。

    3.2 實際地震資料處理

    圖13a為某油田實際疊后地震資料,共有251道,每道651個采樣點,時間采樣率為2ms??梢钥闯?由于吸收衰減作用,深層地震反射同相軸能量相對較弱,波形畸變且連續(xù)性較差(黑色橢圓以及箭頭所示位置)。利用本文方法,首先從地震記錄的Gabor譜中估計衰減函數(shù)譜,在衰減函數(shù)最小相位的假設(shè)下,利用希爾伯特變換求得相位譜,利用公式(9)實現(xiàn)地震記錄的衰減補償。補償后的地震剖面如圖13b所示,深層反射能量相對增強,波形得到一定程度的校正,且連續(xù)性變好,剖面分辨率得到提高。實際地震資料處理驗證了本文方法的有效性。

    圖12 零相位子波合成記錄以及補償結(jié)果a 原始地震記錄; b Q=50對應的衰減地震記錄; c 補償后的地震記錄

    圖13 實際地震資料衰減補償前(a)、后(b)的剖面對比

    4 結(jié)論與認識

    我們基于非穩(wěn)態(tài)地震道模型,利用能量均分雙曲平滑方法從衰減地震記錄的Gabor譜中估計衰減函數(shù),基于反演策略得到補償后的地震數(shù)據(jù)。研究得出以下幾點結(jié)論與認識:

    1) 該方法忽略子波的影響,對子波的相位無要求,適用于任意相位子波合成的地震記錄,克服了Gabor反褶積中子波最小相位假設(shè)的缺陷;

    2) 該方法不需要精確的Q值信息,衰減信息從觀測地震記錄的Gabor譜中提取,克服了常規(guī)反Q濾波需要精確Q值的弊端;

    3) 該方法基于反演思想得到補償后的地震數(shù)據(jù),避免了常規(guī)反Q濾波振幅補償?shù)牟环€(wěn)定性;

    4) 相對傳統(tǒng)雙曲平滑方法,基于能量均分雙曲平滑方法估計衰減譜的精度高,據(jù)此反演得到的衰減補償記錄分辨率更高;

    5) 該方法沒有消除子波的影響,建議對補償后的地震記錄進行穩(wěn)態(tài)盲反褶積處理,同時得到子波信息以及高分辨率的反射系數(shù)序列;

    6) 該方法為單道處理方法,沒有考慮橫向約束,有待拓展到多道處理,提高橫向的連續(xù)性;

    7) 基于衰減函數(shù)最小相位假設(shè)和利用希爾伯特變換求取相位譜時會帶入計算誤差,下一步需要研究如何減弱最小相位假設(shè)和提高精度。

    參 考 文 獻

    [1] Futterman W I.Dispersive body waves[J].Journal of Geophysical Research,1962,67(13):5279-5291

    [2] Varela C L,Rosa A L R,Ulrych T J.Modeling of attenuation and dispersion[J].Geophysics,1993,58(8):1167-1173

    [3] Van der Baan M.Bandwidth enhancement:inverse Q filtering or time-varying Wiener deconvolution?[J].Geophysics,2012,77(4):V133-V142

    [4] Wang Y.A stable and efficient approach of inverse Q filtering[J].Geophysics,2002,67(2):657-663

    [5] Wang Y.Quantifying the effectiveness of stabilized inverse Q filtering[J].Geophysics,2003,68(1):337-345

    [6] Wang Y.Inverse Q-filter for seismic resolution enhancement[J].Geophysics,2006,71(3):V51-V60

    [7] Zhang X,Han L,Zhang F,et al.An inverse Q-filter algorithm based on stable wavefield continuation[J].Applied Geophysics,2007,4(4):263-270

    [8] Hargreaves N,Calvert A.Inverse Q filtering by Fourier transform[J].Geophysics,1991,56(4):519-527

    [9] 趙建勛,倪克森.串聯(lián)反Q濾波及其應用[J].石油地球物理勘探,1992,27(6):722-730

    Zhao J X,Ni K S.Cascaded inverse Q filtering and the application[J].Oil Geophysical Prospecting,1992,27(6):722-730

    [10] Yan H,Liu Y.Estimation of Q and inverse Q filtering for prestack reflected PP- and converted PS-waves[J].Applied Geophysics,2009,6(1):59-69

    [11] Zhao Y,Liu Y,Li X,et al.An approach of inverse Q filtering considering time-frequency domain signal-to-noise ratio[J].74th EAGE Conference & Exhibition,2012,P103

    [12] Braga I L S,Moraes F S.High-resolution gathers by inverse Q filtering in the wavelet domain[J].Geophysics,2013,78(2):V53-V61

    [13] Zhang C,Ulrych T J.Seismic absorption compensation:a least squares inverse scheme[J].Geophysics,2007,72(6):R109-R114

    [14] Wang S.Attenuation compensation method based on inversion[J].Applied Geophysics,2011,8(2):150-157

    [15] Margrave G F,Dong L,Gibson P,et al.Gabor deconvolution:extending Wiener’s method to nonstationarity[R].Calgary:CREWES Research Report,2003

    [16] Margrave G F,Lamoureux M P,Henley D C.Gabor deconvolution:estimating reflectivity by nonstationary deconvolution of seismic data[J].Geophysics,2011,76(3):W15-W30

    [17] Reine C,van der Baan M,Clark R.The robustness of seismic attenuation measurements using fixed-and variable-window time-frequency transforms[J].Geophysics,2009,74(2):WA123-WA135

    [18] Margrave G F.Theory of nonstationary linear filtering in the Fourier domain with application to time-variant filtering[J].Geophysics,1998,63(1):244-259

    [19] Li F,Wang S,Chen X.Gabor deconvolution-hyperbolic smoothing with variable -step sampling[J].75thEAGE Conference & Exhibition,2013,Th-06-04

    猜你喜歡
    子波雙曲位子
    中國科學技術(shù)館之“雙曲隧道”
    軍事文摘(2021年22期)2022-01-18 06:22:48
    一類非線性動力系統(tǒng)的孤立子波解
    換位思考
    雙曲型交換四元數(shù)的極表示
    一階雙曲型偏微分方程的模糊邊界控制
    幸存位子
    基于雙曲和代數(shù)多項式的HC-Bézier曲線
    地震反演子波選擇策略研究
    基于倒雙譜的地震子波估計方法
    魔晶重現(xiàn)
    欧美日韩在线观看h| 久久人人爽人人片av| 国产又色又爽无遮挡免| 80岁老熟妇乱子伦牲交| 成年女人在线观看亚洲视频 | 少妇丰满av| 熟妇人妻久久中文字幕3abv| 国产成人午夜福利电影在线观看| 联通29元200g的流量卡| 亚洲美女搞黄在线观看| 最近手机中文字幕大全| 国产亚洲午夜精品一区二区久久 | 看十八女毛片水多多多| 日韩在线高清观看一区二区三区| 欧美极品一区二区三区四区| 色综合站精品国产| 免费av不卡在线播放| 伦精品一区二区三区| 男人和女人高潮做爰伦理| 18禁在线播放成人免费| 亚洲人成网站在线观看播放| av在线老鸭窝| 国产探花在线观看一区二区| 麻豆av噜噜一区二区三区| 日韩亚洲欧美综合| 国模一区二区三区四区视频| 性插视频无遮挡在线免费观看| 夫妻午夜视频| av又黄又爽大尺度在线免费看| 我的女老师完整版在线观看| 日韩av免费高清视频| 美女黄网站色视频| 天天躁夜夜躁狠狠久久av| 熟女电影av网| 亚洲成人av在线免费| 黄色欧美视频在线观看| 成年人午夜在线观看视频 | xxx大片免费视频| 精品一区二区三区人妻视频| 国产一区二区三区综合在线观看 | 欧美xxⅹ黑人| av在线观看视频网站免费| 91aial.com中文字幕在线观看| 99久久九九国产精品国产免费| 国产精品人妻久久久久久| 亚洲国产精品成人综合色| 色播亚洲综合网| 亚洲国产色片| 中文字幕亚洲精品专区| 日韩视频在线欧美| 亚洲欧美成人精品一区二区| 日韩av不卡免费在线播放| 一级片'在线观看视频| 午夜福利在线观看吧| 亚洲精品久久午夜乱码| 日日啪夜夜撸| 噜噜噜噜噜久久久久久91| av在线蜜桃| 成人欧美大片| 国产成人一区二区在线| 亚洲经典国产精华液单| 免费观看精品视频网站| 亚洲人成网站高清观看| 中文字幕亚洲精品专区| 在现免费观看毛片| 韩国av在线不卡| av女优亚洲男人天堂| 国产亚洲最大av| 国产亚洲av嫩草精品影院| 春色校园在线视频观看| 成人无遮挡网站| 日韩av在线大香蕉| 欧美性猛交╳xxx乱大交人| 国产高潮美女av| 日韩av在线免费看完整版不卡| av网站免费在线观看视频 | 亚洲成人中文字幕在线播放| 国产成人91sexporn| 美女大奶头视频| 亚洲av一区综合| 午夜激情欧美在线| 亚洲欧美日韩东京热| 日韩中字成人| 男插女下体视频免费在线播放| 国产 一区 欧美 日韩| 夜夜看夜夜爽夜夜摸| 最近最新中文字幕免费大全7| 能在线免费看毛片的网站| 午夜久久久久精精品| 人人妻人人澡人人爽人人夜夜 | 熟女人妻精品中文字幕| 晚上一个人看的免费电影| 大片免费播放器 马上看| 校园人妻丝袜中文字幕| 亚洲av成人av| 国产精品久久久久久精品电影| 久久99热这里只有精品18| 69人妻影院| av国产久精品久网站免费入址| 免费在线观看成人毛片| 嫩草影院入口| 99久国产av精品国产电影| 亚洲精品乱久久久久久| 91久久精品电影网| 97超视频在线观看视频| 亚洲欧美精品专区久久| 久久久久精品久久久久真实原创| 五月玫瑰六月丁香| 日韩三级伦理在线观看| 国产亚洲午夜精品一区二区久久 | 淫秽高清视频在线观看| 久久久久九九精品影院| 又黄又爽又刺激的免费视频.| 在线天堂最新版资源| 亚洲欧洲国产日韩| 一个人观看的视频www高清免费观看| 亚洲国产成人一精品久久久| 亚洲,欧美,日韩| 国产精品人妻久久久久久| 国产成人午夜福利电影在线观看| 人人妻人人看人人澡| 80岁老熟妇乱子伦牲交| 免费看不卡的av| 91精品一卡2卡3卡4卡| 久久久久免费精品人妻一区二区| 中国美白少妇内射xxxbb| 国产欧美日韩精品一区二区| 免费观看精品视频网站| 亚洲成色77777| 国产亚洲5aaaaa淫片| 日韩欧美三级三区| av在线播放精品| 五月天丁香电影| 啦啦啦啦在线视频资源| 极品少妇高潮喷水抽搐| 免费观看精品视频网站| 亚洲欧美日韩东京热| 网址你懂的国产日韩在线| 国产黄片美女视频| 午夜福利高清视频| 看非洲黑人一级黄片| 欧美三级亚洲精品| 亚洲成人中文字幕在线播放| 国产精品三级大全| 少妇裸体淫交视频免费看高清| 国产黄片视频在线免费观看| 国产精品熟女久久久久浪| 欧美日韩亚洲高清精品| 亚洲国产色片| av在线老鸭窝| 小蜜桃在线观看免费完整版高清| kizo精华| 秋霞伦理黄片| 三级毛片av免费| 99久国产av精品国产电影| 18禁在线无遮挡免费观看视频| 在线免费观看的www视频| 欧美激情在线99| 三级毛片av免费| 尤物成人国产欧美一区二区三区| 日韩成人伦理影院| 熟妇人妻久久中文字幕3abv| 精品国产露脸久久av麻豆 | 爱豆传媒免费全集在线观看| 亚洲乱码一区二区免费版| 国产亚洲5aaaaa淫片| 麻豆久久精品国产亚洲av| 国产乱人视频| 美女大奶头视频| 成人特级av手机在线观看| 亚洲欧美清纯卡通| 夜夜爽夜夜爽视频| 国产免费又黄又爽又色| 只有这里有精品99| 插阴视频在线观看视频| 波多野结衣巨乳人妻| 日韩伦理黄色片| 五月天丁香电影| 99re6热这里在线精品视频| 搡女人真爽免费视频火全软件| 午夜福利视频1000在线观看| 午夜免费男女啪啪视频观看| 美女黄网站色视频| 亚洲av不卡在线观看| 色网站视频免费| 色哟哟·www| 嘟嘟电影网在线观看| 女人久久www免费人成看片| 爱豆传媒免费全集在线观看| 国产精品国产三级国产专区5o| 夫妻性生交免费视频一级片| 乱人视频在线观看| 亚洲人成网站高清观看| 国内揄拍国产精品人妻在线| 亚洲精品成人av观看孕妇| 高清视频免费观看一区二区 | 欧美97在线视频| 久久99蜜桃精品久久| 天天躁日日操中文字幕| 尾随美女入室| 国产亚洲91精品色在线| 亚洲av在线观看美女高潮| av网站免费在线观看视频 | 中国美白少妇内射xxxbb| 亚洲av.av天堂| 建设人人有责人人尽责人人享有的 | 中文精品一卡2卡3卡4更新| 免费在线观看成人毛片| 成人亚洲欧美一区二区av| 一级毛片电影观看| 成人美女网站在线观看视频| 亚洲欧美成人综合另类久久久| 能在线免费观看的黄片| 久久久久国产网址| 国产精品精品国产色婷婷| 国产在视频线在精品| 亚洲精品视频女| 欧美成人一区二区免费高清观看| 欧美+日韩+精品| 丰满人妻一区二区三区视频av| 麻豆乱淫一区二区| 亚洲欧美一区二区三区黑人 | 九九在线视频观看精品| 欧美性猛交╳xxx乱大交人| 十八禁国产超污无遮挡网站| 亚洲欧美日韩无卡精品| 一级黄片播放器| 亚洲成色77777| 狂野欧美激情性xxxx在线观看| 成年av动漫网址| 亚洲精品aⅴ在线观看| 69人妻影院| 成人欧美大片| 亚洲av免费在线观看| 女人被狂操c到高潮| 国产有黄有色有爽视频| 搡老妇女老女人老熟妇| 亚洲人成网站在线观看播放| 国产精品人妻久久久影院| 一级av片app| 国产高潮美女av| 性插视频无遮挡在线免费观看| 午夜视频国产福利| 欧美精品国产亚洲| 蜜桃久久精品国产亚洲av| 91午夜精品亚洲一区二区三区| 国模一区二区三区四区视频| 白带黄色成豆腐渣| 嫩草影院新地址| 欧美潮喷喷水| 久久久久久久久久黄片| 色尼玛亚洲综合影院| 国产精品福利在线免费观看| 人妻少妇偷人精品九色| 日本黄大片高清| 毛片一级片免费看久久久久| 久久久久久久久中文| 亚洲一级一片aⅴ在线观看| 特大巨黑吊av在线直播| 只有这里有精品99| 亚洲精品第二区| 国产av码专区亚洲av| 视频中文字幕在线观看| 亚洲欧美成人精品一区二区| 久久久精品94久久精品| 免费少妇av软件| 久久这里只有精品中国| 国产v大片淫在线免费观看| 欧美精品一区二区大全| 免费观看的影片在线观看| 免费大片黄手机在线观看| 亚洲国产欧美人成| 成年女人在线观看亚洲视频 | 日本一二三区视频观看| h日本视频在线播放| 最近视频中文字幕2019在线8| 丝袜喷水一区| 菩萨蛮人人尽说江南好唐韦庄| 国产中年淑女户外野战色| 亚洲欧美一区二区三区黑人 | 精品人妻偷拍中文字幕| 久久久色成人| 国产av码专区亚洲av| 国产av国产精品国产| 在线播放无遮挡| 高清av免费在线| 亚洲精品一区蜜桃| 亚洲怡红院男人天堂| 美女高潮的动态| 成人毛片60女人毛片免费| 超碰av人人做人人爽久久| 大又大粗又爽又黄少妇毛片口| 亚洲国产精品专区欧美| 国产91av在线免费观看| 26uuu在线亚洲综合色| 午夜激情福利司机影院| 成人毛片60女人毛片免费| 一级毛片我不卡| 国产成人精品福利久久| 欧美激情在线99| 高清视频免费观看一区二区 | 国产白丝娇喘喷水9色精品| 青青草视频在线视频观看| 久久精品熟女亚洲av麻豆精品 | 成年av动漫网址| 亚洲熟女精品中文字幕| 成年免费大片在线观看| 神马国产精品三级电影在线观看| 99热6这里只有精品| 欧美xxⅹ黑人| 在现免费观看毛片| 国产精品一区二区性色av| 插阴视频在线观看视频| 国产精品久久久久久久电影| 丰满少妇做爰视频| 波野结衣二区三区在线| videossex国产| 亚洲色图av天堂| av在线播放精品| 欧美人与善性xxx| 久久久精品欧美日韩精品| 日本免费在线观看一区| 欧美日韩国产mv在线观看视频 | 丰满乱子伦码专区| 久久精品国产亚洲网站| 免费少妇av软件| 精品一区二区免费观看| 亚洲欧美中文字幕日韩二区| 亚洲欧美成人精品一区二区| 日韩制服骚丝袜av| 国产人妻一区二区三区在| 国产精品一区二区三区四区免费观看| 一级片'在线观看视频| 国产91av在线免费观看| 日韩人妻高清精品专区| 国产综合精华液| 亚洲av成人精品一区久久| 五月伊人婷婷丁香| 国产精品久久久久久精品电影小说 | 又粗又硬又长又爽又黄的视频| 99久国产av精品| 黄色一级大片看看| 亚洲熟妇中文字幕五十中出| 晚上一个人看的免费电影| 国产女主播在线喷水免费视频网站 | 亚洲一区高清亚洲精品| 黄色配什么色好看| 国产精品一区二区三区四区久久| 国产精品伦人一区二区| 两个人的视频大全免费| 日日干狠狠操夜夜爽| 亚洲在线观看片| 美女大奶头视频| 精品欧美国产一区二区三| 亚洲在久久综合| 亚洲内射少妇av| 成年av动漫网址| 亚洲精品亚洲一区二区| 成人漫画全彩无遮挡| 日韩视频在线欧美| 亚洲成人精品中文字幕电影| 午夜福利视频1000在线观看| 精品少妇黑人巨大在线播放| 亚洲18禁久久av| 久久久色成人| 免费人成在线观看视频色| 亚洲国产精品sss在线观看| 欧美另类一区| 亚洲国产最新在线播放| 人妻一区二区av| 亚洲一级一片aⅴ在线观看| 国产v大片淫在线免费观看| 国内少妇人妻偷人精品xxx网站| 国模一区二区三区四区视频| 国产伦一二天堂av在线观看| 又粗又硬又长又爽又黄的视频| 久久久久久久国产电影| 亚洲精品自拍成人| 国产亚洲精品久久久com| 不卡视频在线观看欧美| 色综合站精品国产| 中文字幕人妻熟人妻熟丝袜美| 深夜a级毛片| 男人和女人高潮做爰伦理| 亚洲欧美日韩东京热| 男女国产视频网站| 国产 一区 欧美 日韩| 日本午夜av视频| 最近2019中文字幕mv第一页| .国产精品久久| 18禁在线播放成人免费| 大又大粗又爽又黄少妇毛片口| 26uuu在线亚洲综合色| 99热全是精品| 国产熟女欧美一区二区| av在线天堂中文字幕| 亚洲在线观看片| 亚洲国产av新网站| 亚洲av电影不卡..在线观看| 22中文网久久字幕| 午夜福利视频1000在线观看| 91狼人影院| 99热网站在线观看| 麻豆av噜噜一区二区三区| 熟女电影av网| 嫩草影院新地址| 亚洲人成网站在线观看播放| 国产精品国产三级国产av玫瑰| 亚洲精品中文字幕在线视频 | 欧美日韩在线观看h| 能在线免费观看的黄片| 夜夜看夜夜爽夜夜摸| 91久久精品国产一区二区成人| 午夜免费观看性视频| 久久久久精品久久久久真实原创| 久久人人爽人人片av| 精品国产露脸久久av麻豆 | 色尼玛亚洲综合影院| 少妇人妻一区二区三区视频| 亚洲国产精品成人久久小说| 亚洲经典国产精华液单| 日本av手机在线免费观看| 99久久精品一区二区三区| 97超视频在线观看视频| 欧美高清性xxxxhd video| 国产av在哪里看| 乱系列少妇在线播放| 欧美xxⅹ黑人| 亚洲高清免费不卡视频| 精品酒店卫生间| 久久97久久精品| 亚洲在久久综合| 亚洲av电影在线观看一区二区三区 | 少妇人妻精品综合一区二区| 精品熟女少妇av免费看| 在线观看美女被高潮喷水网站| 成人美女网站在线观看视频| 亚洲国产av新网站| 晚上一个人看的免费电影| 中文欧美无线码| 狂野欧美白嫩少妇大欣赏| 日本黄大片高清| 能在线免费看毛片的网站| 中文资源天堂在线| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美成人a在线观看| 国产伦理片在线播放av一区| 一级a做视频免费观看| 国产免费又黄又爽又色| 国产精品国产三级国产av玫瑰| 别揉我奶头 嗯啊视频| 看非洲黑人一级黄片| 国产精品国产三级国产专区5o| 久久热精品热| 男人爽女人下面视频在线观看| 国产毛片a区久久久久| 欧美性感艳星| 亚洲精品一二三| 国内精品宾馆在线| 伦理电影大哥的女人| 亚洲av一区综合| 国产精品麻豆人妻色哟哟久久 | 午夜精品国产一区二区电影 | 黄色日韩在线| 69人妻影院| 美女被艹到高潮喷水动态| 一级毛片 在线播放| 日本三级黄在线观看| 哪个播放器可以免费观看大片| 国产亚洲av嫩草精品影院| 三级男女做爰猛烈吃奶摸视频| 久久韩国三级中文字幕| 天天一区二区日本电影三级| 18禁在线播放成人免费| 日韩在线高清观看一区二区三区| 一夜夜www| 国产精品久久久久久精品电影小说 | 亚洲av电影在线观看一区二区三区 | 欧美日韩亚洲高清精品| 欧美日韩在线观看h| 成人鲁丝片一二三区免费| 边亲边吃奶的免费视频| 久久久久网色| 国产av码专区亚洲av| av在线播放精品| 免费大片黄手机在线观看| 最近中文字幕高清免费大全6| 国产成人一区二区在线| 中国国产av一级| 日韩在线高清观看一区二区三区| 国产国拍精品亚洲av在线观看| 联通29元200g的流量卡| 国产av不卡久久| 又黄又爽又刺激的免费视频.| 国产精品一及| 久久久a久久爽久久v久久| 国产精品国产三级国产av玫瑰| 乱码一卡2卡4卡精品| 免费电影在线观看免费观看| 亚洲精品乱码久久久v下载方式| 精品久久久久久电影网| 色尼玛亚洲综合影院| 亚洲美女视频黄频| 色尼玛亚洲综合影院| 免费av毛片视频| 十八禁国产超污无遮挡网站| 亚洲av免费在线观看| 一级爰片在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日本一二三区视频观看| 国产有黄有色有爽视频| 国产亚洲精品久久久com| 欧美xxⅹ黑人| 精品一区二区三区视频在线| 国产高潮美女av| 精品人妻熟女av久视频| 日韩av在线免费看完整版不卡| 少妇人妻精品综合一区二区| 亚洲精品日韩在线中文字幕| 久久人人爽人人片av| 免费高清在线观看视频在线观看| 国产视频内射| 黄片无遮挡物在线观看| 女的被弄到高潮叫床怎么办| av国产免费在线观看| 亚洲精品国产av蜜桃| 日韩不卡一区二区三区视频在线| 久热久热在线精品观看| 美女大奶头视频| 国产91av在线免费观看| 精品一区二区免费观看| 观看免费一级毛片| 又粗又硬又长又爽又黄的视频| 精品国内亚洲2022精品成人| 午夜激情久久久久久久| 大又大粗又爽又黄少妇毛片口| 精品久久久久久久久av| 午夜免费观看性视频| or卡值多少钱| 亚洲精华国产精华液的使用体验| 亚洲精品国产成人久久av| 久久亚洲国产成人精品v| 日韩成人av中文字幕在线观看| 黄色一级大片看看| videossex国产| 亚洲电影在线观看av| 亚洲国产精品成人久久小说| 亚洲精品自拍成人| 日本黄色片子视频| 在线观看一区二区三区| videos熟女内射| 国产黄a三级三级三级人| 成人欧美大片| 久久久久久国产a免费观看| 久久久精品免费免费高清| 激情 狠狠 欧美| 亚洲最大成人中文| 国模一区二区三区四区视频| 国产午夜福利久久久久久| 51国产日韩欧美| 日本与韩国留学比较| 久久久久久久国产电影| 国产不卡一卡二| 亚洲天堂国产精品一区在线| 一级片'在线观看视频| 久久久欧美国产精品| 久久综合国产亚洲精品| 伦理电影大哥的女人| 国产亚洲av嫩草精品影院| 亚洲精品亚洲一区二区| 99久国产av精品| 最近最新中文字幕免费大全7| 国产综合懂色| 日本猛色少妇xxxxx猛交久久| 中文字幕av在线有码专区| 观看美女的网站| 精品不卡国产一区二区三区| 亚洲国产色片| 人人妻人人澡欧美一区二区| 欧美成人精品欧美一级黄| 国产精品久久久久久精品电影| 成人一区二区视频在线观看| 成人特级av手机在线观看| 搞女人的毛片| 最后的刺客免费高清国语| 国内少妇人妻偷人精品xxx网站| 又爽又黄无遮挡网站| 99热网站在线观看| 成人综合一区亚洲| 一级二级三级毛片免费看| 成人亚洲精品一区在线观看 | 国产老妇女一区| 乱系列少妇在线播放| 深夜a级毛片| 国产综合懂色| 久久久久久久午夜电影| 伊人久久精品亚洲午夜| 丰满少妇做爰视频| 午夜福利高清视频| 极品教师在线视频| 99热网站在线观看| 久久久精品免费免费高清| 嫩草影院新地址| 国产精品久久久久久久电影| 97精品久久久久久久久久精品| 少妇被粗大猛烈的视频| 欧美一级a爱片免费观看看| 99九九线精品视频在线观看视频| .国产精品久久| 亚洲欧美一区二区三区国产| 精品久久久久久久久亚洲| av线在线观看网站| 国产亚洲午夜精品一区二区久久 | 全区人妻精品视频| 熟女人妻精品中文字幕|