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

    時空域同步擠壓小波變換聯(lián)合面波壓制方法

    2022-08-31 13:07:32林浩然邢磊劉懷山李倩倩張洪茂
    地球物理學(xué)報 2022年9期
    關(guān)鍵詞:面波波數(shù)壓制

    林浩然,邢磊,2*,劉懷山,2,李倩倩,張洪茂

    1 中國海洋大學(xué)海底科學(xué)與探測技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,青島 266100 2 青島海洋科學(xué)與技術(shù)試點(diǎn)國家實(shí)驗(yàn)室海洋礦產(chǎn)資源評價與探測技術(shù)功能實(shí)驗(yàn)室,青島 266071

    0 引言

    面波是地震資料中最為常見的規(guī)則干擾波,其是由體波到達(dá)地表或介質(zhì)分界面時,在一定條件下相互干涉并疊加產(chǎn)生的次聲波.面波具有低頻、低速、強(qiáng)振幅、能量衰減緩慢等特點(diǎn),在炮集記錄中呈近似線性分布,面波的存在嚴(yán)重影響中深層有效信號,顯著降低了地震資料的信噪比(王偉奇等,2022).

    在地震資料處理過程中,研究人員主要依據(jù)面波低頻、低視速度的特點(diǎn),提出了許多壓制面波的方法,如高通濾波、FK濾波、Radon域?yàn)V波、KL變換域?yàn)V波、維納濾波、SVD濾波、徑向道分解濾波、二維時間導(dǎo)數(shù)濾波等.這些方法雖具有一定的壓制效果,但是往往會對有效信號造成一定程度的損傷.為了獲得更好的面波壓制效果,提高地震資料信噪比,Deighan和Watts (1997)提出將一維小波變換應(yīng)用于面波壓制工作中.小波變換具有良好時頻分辨能力,以及信號奇異性捕獲能力(陳永芮等,2021).近年來,在小波變換的基礎(chǔ)上,具有更好逼近高維奇異信號能力的多尺度幾何分析方法,如脊波(Ridgelet)變換(Candès and Donoho,1999;包乾宗等,2007)、曲波(Curvelet)變換(Candès and Donoho,2004)等被引入面波壓制方法研究領(lǐng)域中來.其優(yōu)越之處在于能更好的利用面波與有效波的視速度差異,捕捉地震信號的奇異特性,同時能更好的刻畫圖像邊緣信息.然而這些時頻分析方法都受到海森堡—蓋博測不準(zhǔn)原理的影響,導(dǎo)致其時頻分辨能力受到一定限制.為克服這一問題并提高它們表示非平穩(wěn)信號的能力,基于銳化時頻表示技術(shù)(Kodera et al.,1976)和時頻重分配思想(Auger and Flandrin,1995),Daubechies和Maes(1996)、Daubechies等(2011)提出了同步擠壓算法,并證明其可以應(yīng)用于地震數(shù)據(jù)處理領(lǐng)域.劉晗等(2016)提出將同步擠壓小波變換應(yīng)用于面波壓制中.時頻重分配算法現(xiàn)階段已被推廣至二維(Yang et al.,2015;Lu et al.,2016)和三維空間(Zhu et al.,2021)中,且有多種變體(Clause et al.,2015;Lu and Yang,2018)和改進(jìn)(陳文超等,2009;Oberlin and Meignen,2017;Daubechies et al.,2016).

    隨著處理要求的提高和處理技術(shù)的進(jìn)步,近幾年來提出了一些采用多種手段聯(lián)合壓制面波的方法.包乾宗等(2011)提出了聯(lián)合小波變換和脊波變換的面波壓制方法.馬見青和李慶春(2011)提出了S變換和TT變換聯(lián)合壓制面波.畢云云等(2017)提出了基于離散曲波變換字典和二維局部離散余弦變換字典聯(lián)合壓制面波的方法.徐陽等(2018)提出了廣義S變換與二維離散小波變換聯(lián)合壓制面波方法,該方法對包含面波成分的低頻信號再次進(jìn)行二維離散小波變換,成功提取出了低頻信息里包含的有效信號,減少了在有效壓制面波干擾的同時對有效反射波的損害.路鵬飛等(2020)提出了曲波變換與傅里葉變換聯(lián)合壓制面波方法.這些多手段聯(lián)合壓制面波的方法,在盡可能保護(hù)有效信號的原則下,最大限度地將有效信號和面波分離,進(jìn)一步提升了地震資料的信噪比.

    本文首先在前人研究的基礎(chǔ)上將同步擠壓小波變換推廣至一維空間域,后根據(jù)面波干擾頻率低、視速度低的特點(diǎn),采用高精度時頻分析方法同步擠壓小波時域和空域變換聯(lián)合壓制面波干擾.從時間、空間、頻率、波數(shù)維度對地震數(shù)據(jù)進(jìn)行分析,利用時域同步擠壓小波變換結(jié)合時頻域?yàn)V波器將面波與有效信號分離,并對分離出的面波數(shù)據(jù)體使用空域同步擠壓小波變換進(jìn)行低頻有效信號的重提取,通過理論模型測試驗(yàn)證了所提聯(lián)合方法的有效性,最后通過實(shí)際資料處理與一些常用方法的處理結(jié)果相比較,證明了所提出的聯(lián)合方法處理結(jié)果信噪比高、保幅性好的優(yōu)點(diǎn).

    1 方法原理

    1.1 一維空間域同步擠壓小波變換

    同步擠壓小波變換(Synchrosqueezing Wavelet Transform, SWT)是由Daubechies等(2011)引入的一種針對非平穩(wěn)信號的稀疏表示.SWT是基于小波變換的時頻重分配,其類似于經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition, EMD)可以將信號分解為本征模態(tài)函數(shù),優(yōu)于EMD的是,SWT具有更堅實(shí)的理論基礎(chǔ)(Thakur et al.,2013).

    同步擠壓算法通過沿頻率軸壓縮信號的能量來聚焦時頻表示,這種后處理方法的時頻分辨率極高且具有嚴(yán)格的連續(xù)可逆性.該變換可將時變信號分解為如式(1)所示的固有模態(tài):

    (1)

    式中,Am(t)為瞬時振幅,η(t)代表噪聲或誤差,M為模態(tài)數(shù)量,φm(t)為瞬時相位.每個模態(tài)的瞬時頻率是通過瞬時相位的偏差(Boashash,2003)來計算的,計算方程為:

    (2)

    現(xiàn)進(jìn)行空間域同步擠壓小波變換推導(dǎo),給定一空變信號s(x):

    s(x)=A(x)cos(kx+φ0),(3)

    式中,x為波動距離,φ0為初相,A為振幅,(kx+φ0)為信號相位,k為空域角頻率也稱為角波數(shù):

    k=2πξ,(4)

    式中,ξ為空間頻率.

    若需要獲取空變信號的空間頻率即波數(shù),通常對于空間采樣信號通過傅里葉變換來進(jìn)行空間-波數(shù)域變換,得到此空間采樣信號的波數(shù)譜,其一般定義式為:

    (5)

    給定母小波函數(shù)φ,對空變信號s(x)進(jìn)行連續(xù)小波變換,得到空間-尺度域小波系數(shù)Ws(a,b)(Stéphane,2009):

    (6)

    式中a為尺度序列,b為空間偏移.根據(jù)Plancherel定理,式(6)可以寫成波數(shù)域表達(dá)形式:

    (7)

    (8)

    對Ws(a,b)來說,|Ws(a,b)|表示各空間模態(tài)分量的小波系數(shù)(Li and Liang,2012):

    (9)

    和時間偏移τ一樣,當(dāng)給定一個空間偏移b,信號能量主要沿尺度軸擴(kuò)散.忽略沿空間軸的能量擴(kuò)散,對任意Ws(a,b)不為零的尺度-空間點(diǎn)(a,b),其瞬時空間頻率估計表示為:

    (10)

    當(dāng)把地震信號同一時刻各道的數(shù)據(jù)序列看作是反映一定范圍內(nèi)地層及其地質(zhì)參數(shù),便構(gòu)成了地震信號的空間序列.地震信號的空間序列與時/深剖面是一一對應(yīng)的,對于某一探區(qū),不同時刻的空間序列信號對應(yīng)于相應(yīng)時/深剖面的不同深度地層位置(黨建武和黃建國,2003).

    其實(shí)際情況下的一般離散型式為:

    (11)

    式中,i為時序采樣點(diǎn),q為地震記錄的道號.

    (12)

    離散形式的瞬時頻率估計為:

    (13)

    需要注意的是,式(13)中得到的瞬時空間頻率估計往往并不是瞬時空間頻率本身(只有在輸入信號為理想簡諧信號時才為瞬時頻率),他是我們用來計算實(shí)際空間頻率的校正量.

    基于上述概念推導(dǎo)出空間同步擠壓小波變換(Spatial Domain Synchrosqueezing Wavelet Transform, SSWT)的離散形式為:

    (14)

    這里把a(bǔ)和ω看成連續(xù)函數(shù),空間同步擠壓小波變換的連續(xù)形式為:

    (15)

    式中,A(x)={a;Ws(a,x)≠0}.

    式(14)、(15)將時頻重分配思想推廣至空間-波數(shù)域,且同步擠壓只發(fā)生在空間頻率軸.

    空間域同步擠壓小波變換得到的同步擠壓空間波數(shù)表示,可以使用類似于時域同步擠壓小波變換的重構(gòu)方式,重構(gòu)得到原始信號:

    (16)

    如前文所述,野外采集得到的地震信號是時空維度的,通過相速度我們可以建立起地震時域序列與空域序列的聯(lián)系:

    (17)

    式中,k為角波數(shù),Ω為角頻率,Cp為相速度.由于地震橫波如面波具有很強(qiáng)的頻散現(xiàn)象其傳播速度取相速度,地震體波的相速度和群速度基本相差不大,因此我們可以將式(17)中的相速度近似規(guī)定為地震波傳播速度,且地震波波速越大,其角波數(shù)越小.對于任意t時刻,有:

    (18)

    式中,Ω、φt為常數(shù).

    信號s(x)的連續(xù)小波變換(Continuous wavelet transform, CWT)為:

    (19)

    瞬時空間頻率為:

    (20)

    信號s(x)的SSWT為:

    (21)

    通過式(19)、(21)建立起了CWT、SSWT與地震波波速之間的聯(lián)系,并將其展示在空間波數(shù)譜中.

    由于面波具有視速度低的特性,其在炮集記錄中呈現(xiàn)近乎垂直的圓錐狀,與有效信號差異明顯,因此在對地震數(shù)據(jù)進(jìn)行空間域同步擠壓小波變換時,面波干擾將主要出現(xiàn)在高波數(shù)區(qū),據(jù)此可選合適濾波方式對低頻段面波干擾和有效信號分離,恢復(fù)部分低頻有效信號.

    1.2 時域同步擠壓小波變換與空域同步擠壓小波變換聯(lián)合面波壓制方法

    綜合上述方法原理,針對面波與有效信號的振幅、頻率、視速度差異,采用時域同步擠壓小波變換與空域同步擠壓小波變換聯(lián)合方法壓制面波干擾,具體實(shí)現(xiàn)步驟如下:

    (1)對目標(biāo)單炮記錄進(jìn)行頻譜分析,確定有效信號和面波干擾的頻帶.選取最優(yōu)母小波函數(shù),按道對輸入的單炮記錄進(jìn)行時域同步擠壓小波正變換(Time Domain Synchrosqueezing Wavelet Transform, TSWT),得到每道信號對應(yīng)的時頻譜.

    (2)根據(jù)面波和有效信號在時頻譜上能量顯示差異,面波能量位于低頻段沿時間軸聚集,并呈條帶狀分布,有效信號能量沿頻率軸聚集且沿頻率軸正方向收斂,同時結(jié)合頻譜分析得到的面波與有效信號各自的頻帶,在時頻譜上確定面波分布區(qū).

    (3)選擇合適時頻窗口將面波區(qū)域內(nèi)的數(shù)據(jù)充零.后進(jìn)行時域同步擠壓小波反變換,得到去除面波干擾的單道時域信號,按道依次重復(fù)步驟1~3,即可得到去除面波干擾后的地震記錄.同時,將原始單炮記錄與此記錄匹配相減,得到一次濾波濾除的面波數(shù)據(jù)體.

    (4)對分離出來的低頻數(shù)據(jù)體按時序進(jìn)行SSWT,選擇合適空頻窗口將面波區(qū)域內(nèi)的系數(shù)充零,反變換為各時刻點(diǎn)對應(yīng)的空變信號,即提取出在時頻域?yàn)V波中損傷的部分有效信號.將其疊加到第3步得到的一次濾波后的地震記錄中,完成了對面波干擾的聯(lián)合壓制,有效保護(hù)了低頻段反射波信號.

    本文基于Gauss、Bump、Cmhat、Morlet四種母小波函數(shù)進(jìn)行了信號重構(gòu)(重構(gòu)誤差見圖1).結(jié)果表明Morlet小波的重構(gòu)誤差遠(yuǎn)小于其他小波,且能較好地表征信號時頻特征(圖2).因此本文選用Morlet小波作為母小波函數(shù).

    圖1 四種母小波函數(shù)重構(gòu)誤差比較Fig.1 Comparison of reconstruction errors of four mother wavelet functions

    圖2 四種母小波函數(shù)的TSWT時頻譜Fig.2 The time frequency spectrum of TSWT of four mother wavelet functions

    2 理論模型測試

    為驗(yàn)證所提方法的有效性,首先對正演得到的理論模型進(jìn)行處理,模型參數(shù)見表1.

    表1 地層模型參數(shù)Table 1 Formation model parameters

    圖3 理論模型合成記錄Fig.3 Theoretical model synthesis record

    圖3為采用彈性波波場模擬方式生成的理論模型合成記錄,單震源單邊激發(fā),震源位于地表,二維均勻各項(xiàng)同性介質(zhì),共設(shè)計了2層水平反射層,采樣間隔為1 ms,共有750個采樣點(diǎn),總記錄長度為0.75 s,80道接收,最小偏移距0 m,道間距6.25 m,最大偏移距493.75 m.面波主頻為20 Hz,反射波主頻為50 Hz,整個合成記錄不含隨機(jī)噪聲.

    圖4為理論模型第35道數(shù)據(jù)進(jìn)行TSWT后得到的時頻譜.時頻分析后可得,根據(jù)面波與有效信號能量聚集差異和頻率差異可以將其區(qū)分,但面波與有效信號在部分頻段存在重合.因此,當(dāng)進(jìn)行面波區(qū)充零處理時將損傷部分低頻有效信號(徐陽等,2018).

    圖4 理論模型合成記錄第35道的時頻譜Fig.4 The time frequency spectrum of channel 35 of theoretical model synthesis record

    對理論模型合成記錄進(jìn)行TSWT,確定面波區(qū)域并將其充零,得到不含面波干擾的模型記錄,如圖5a所示.其與原始理論模型合成記錄匹配相減,得到一次濾波后濾出的面波干擾,如圖5b所示.圖5a與圖3對比發(fā)現(xiàn),其反射波的同相軸振幅發(fā)生了改變,分析可得是由于充零過程導(dǎo)致低頻有效信號收到了損傷,這部分有效信號在圖5b中也得到了體現(xiàn).

    圖5 (a) 經(jīng)TSWT壓制面波后的理論模型記錄; (b) TSWT分離出的面波記錄; (c) 經(jīng)SSWT重構(gòu)出的有效信號; (d) 經(jīng)聯(lián)合方法處理所得到的最終記錄Fig.5 (a) Record of theoretical model after suppressing ground roll by TSWT; (b) Ground roll separated by TSWT; (c) Valid signal reconstructed by SSWT; (d) Final record of theoretical model after suppressing ground roll by combination method

    為恢復(fù)被損傷的低頻有效信號,對經(jīng)TSWT后分離出的面波數(shù)據(jù)體再次按時序?qū)臻g序列進(jìn)行SSWT,得到空間波數(shù)譜.如圖6所示為0.32 s時刻的空間波數(shù)譜對比.

    可見對0.32 s的空變信號序列進(jìn)行同步擠壓小波變換后得到的空間波數(shù)譜較連續(xù)小波變換后能量沿波數(shù)軸有了明顯的集中,同步擠壓小波變換后的空頻分辨率遠(yuǎn)高于連續(xù)小波變換.如圖6c所示,偏移距在0~100 m能量較強(qiáng)且基本對應(yīng)同一波數(shù),結(jié)合如圖5的輸入數(shù)據(jù)不難得到偏移距為0~100 m時有效信號同相軸接近水平,說明其視速度相差不大,因此對應(yīng)在空間波數(shù)譜上的波數(shù)接近一致,而面波視速度較低且在同相軸不水平連續(xù),因此在空間波數(shù)譜上位于高波數(shù)區(qū)域且能量軸傾斜發(fā)散.這正是在空間波數(shù)譜上區(qū)分面波和有效信號的依據(jù).對分離的面波數(shù)據(jù)體按時序進(jìn)行SSWT,抽取同樣時刻數(shù)據(jù)得到的空間波數(shù)譜如圖6f所示.在時頻域被損傷的有效信號即模型中的兩套水平層位的反射波部分能量對應(yīng)一波數(shù)軸刻度集中呈長條狀分布,而面波干擾能量團(tuán)在空間波數(shù)譜中高波數(shù)區(qū)呈傾斜“堅果”狀分布,如圖6e、f中紅框內(nèi)所示.圖6充分說明了SSWT的進(jìn)行空頻分析的可行性和其較連續(xù)小波變換更卓越的空頻分辨能力.

    將面波區(qū)充零,并重構(gòu)即可提取低頻段部分被損傷的有效信號.據(jù)圖6所示,面波干擾得到有效壓制同時有效信號的同相軸連續(xù)性和形態(tài)得到極大改善,與理論模型合成記錄接近.同時,本文還使用TSWT對各方法的濾波結(jié)果進(jìn)行了單頻地震數(shù)據(jù)提取,提取出的四組數(shù)據(jù)從低頻至高頻涵蓋了面波與有效信號的優(yōu)勢頻帶,并計算了各目標(biāo)頻段以及整個炮集記錄的的信噪比(具體算法請見附錄).如圖7所示,對比原始記錄、TSWT濾波后記錄和使用聯(lián)合方法所得各有效頻段信噪比結(jié)果可得,經(jīng)聯(lián)合方法處理后得到的地震記錄的分頻記錄和整體記錄的信噪比均高于前者,尤其是在低頻段.以上結(jié)果證明了SSWT有效恢復(fù)在TSWT濾波過程中被損傷的部分反射波,TSWT和SSWT的聯(lián)合方法提升了地震記錄的信噪比.

    圖6 0.32 s時刻信號空間波數(shù)譜對比(a) 面波分離后合成記錄的波數(shù)-歸一化振幅譜; (b) 面波分離后合成記錄經(jīng)空間連續(xù)小波變換后得到的空間波數(shù)譜; (c) 面波分離后合成記錄經(jīng)SSWT后得到的空間波數(shù)譜; (d) 分離出的面波記錄的波數(shù)-歸一化振幅譜; (e) 分離出的面波記錄經(jīng)空間連續(xù)小波變換后得到的空間波數(shù)譜; (f) 分離出的面波記錄經(jīng)SSWT后得到的空間波數(shù)譜.Fig.6 Comparison of spatial wavenumber spectrum of signal at 0.32 s(a) The wavenumber-normalized amplitude spectrum of the synthetic record after ground roll separation; (b) The spatial wavenumber spectrum of the synthetic record after ground roll separation by Spatial Domain Continuous Wavelet Transform (SCWT); (c) The spatial wavenumber spectrum of the synthetic record after ground roll separation by SSWT; (d) The wavenumber-normalized amplitude spectrum of the ground roll; (e) The spatial wavenumber spectrum of the ground roll by SCWT; (f) The spatial wavenumber spectrum of the ground roll by SSWT.

    圖8 聯(lián)合方法面波壓制效果對比圖(a) Alaska地區(qū)第20炮原始記錄; (b) 時域同步擠壓小波變換壓制面波干擾后的記錄; (c) 使用聯(lián)合方法壓制面波干擾后的記錄.Fig.8 Comparison of ground roll suppression effects of combined method(a) The original record of the 20th source in Alaska;(b) The record after ground roll suppression by TSWT;(c) The record after ground roll suppression by combined method.

    圖9 SSWT的輸入與輸出(a) TSWT后濾除的低頻記錄; (b) SSWT重構(gòu)結(jié)果.Fig.9 Input and output of SSWT(a) Low frequency record filtered byTSWT;(b) Reconstruction result of SSWT.

    圖10 時空域同步擠壓小波變換切片(a) TSWT濾波前15 Hz頻率切片; (b) TSWT濾波后15 Hz頻率切片; (c) SSWT濾波前0.001 rad·m-1波數(shù)切片; (d) SSWT濾波后0.001 rad·m-1波數(shù)切片.Fig.10 Slices of SSWT and TSWT(a) 15 Hz frequency slice before TSWT; (b) 15 Hz frequency slice after TSWT; (c) 0.001 rad·m-1 wavenumber slice before SSWT; (d) 0.001 rad·m-1 wavenumber slice after SSWT.

    圖11 原始單炮記錄第30道TSWT時頻譜和分離出的低頻數(shù)據(jù)體第320 ms SSWT空間波數(shù)譜(a)和(c) 面波區(qū)充零前; (b)和(d) 面波區(qū)充零后.Fig.11 The time frequency spectrum of channel 30 of TSWT and spatial wavenumber spectrum at 320 ms of SSWT of the separated low-frequency data(a) and (c) Before zero filling in the ground roll area; (b) and (d) After zero filling in the ground roll area.

    圖12 3組數(shù)據(jù)的頻譜對比Fig.12 Spectrum comparison of three datasets

    3 實(shí)際資料處理

    本文選取美國Alaska地區(qū)某一測線實(shí)際數(shù)據(jù)進(jìn)行處理,震源采用炸藥震源,總道數(shù)5568道(曾祥堃,2015).抽取第20炮作為輸入單炮記錄,每炮共96道,道間距33.53 m,最小偏移距55 m,最大偏移距5525 m,3000個采樣點(diǎn),采樣間隔2 ms,記錄長度5 s,如圖8a所示,面波呈“圓錐狀”存在,近偏移距的淺層反射波和遠(yuǎn)偏移距中、深層反射波與面波干擾重疊,有效信號難以識別,信噪比嚴(yán)重降低.按道依次進(jìn)行TSWT,得到各道的時頻譜,結(jié)合數(shù)據(jù)分析結(jié)果,確定面波分布區(qū)和時頻濾波窗口將其充零.整個記錄的面波分離完成后,將得到的不含面波干擾的地震記錄(圖8b),與原始記錄匹配相減,即可得到經(jīng)TSWT后濾除的低頻數(shù)據(jù)體,如圖9a所示.低頻數(shù)據(jù)體中含有若干條反射波同相軸,結(jié)合圖10a、b的TSWT濾波前后的15Hz頻率切片來看,分析該現(xiàn)象是由于低頻段有效信號與面波存在重疊,且為徹底壓制面波干擾而擴(kuò)大使用時頻濾波窗口等綜合原因?qū)е?

    為減少對有效信號的損害,將時域同步擠壓小波變換分離出的低頻數(shù)據(jù)體按時序進(jìn)行SSWT,得到空間波數(shù)譜,此時空域的地震數(shù)據(jù)將轉(zhuǎn)換至空間、波數(shù)域中.在空間波數(shù)譜中確定面波區(qū)將其充零,如圖11所示.圖10c、d為SSWT濾波前后空頻濾波窗口內(nèi)0.001 rad·m-1波數(shù)切片,濾波前對應(yīng)的重構(gòu)結(jié)果未混有有效信號,濾波后該波數(shù)對應(yīng)的面波干擾基本被完全壓制,證實(shí)了有效信號與面波干擾視速度差異的SSWT方法的優(yōu)越性.

    有效信號提取結(jié)果如圖9b所示,可見明顯同相軸.將重構(gòu)結(jié)果疊加到TSWT壓制面波干擾后的記錄中,得到聯(lián)合面波壓制方法的最終結(jié)果,如圖8c所示.圖8c中反射波能量尤其是面波覆蓋區(qū)和中深層段,顯著增強(qiáng).通過圖12中3組數(shù)據(jù)頻譜對比,可以更直觀地看出使用聯(lián)合方法后,低頻段能量被有效恢復(fù).根據(jù)圖13信噪比對比結(jié)果來看,經(jīng)聯(lián)合方法處理后的低頻段數(shù)據(jù)以及整體地震數(shù)據(jù)的信噪比均高于原始數(shù)據(jù)和單一使用TSWT處理后的結(jié)果.從實(shí)際資料的處理結(jié)果可以看出本文所述面波壓制聯(lián)合方法在有效壓制面波的同時能極大程度地保留低頻段有效信號.

    4 討論

    在這一部分中,主要討論時域和空域同步擠壓小波變換聯(lián)合方法的優(yōu)勢和通過實(shí)際結(jié)果呈現(xiàn)出的該方法尚未解決的缺點(diǎn)以及未來的研究方向.

    4.1 聯(lián)合方法實(shí)際應(yīng)用優(yōu)勢

    時域和空域同步擠壓小波變換聯(lián)合面波壓制方法是高精度時頻分析方法在地震資料處理面波壓制領(lǐng)域的綜合應(yīng)用.如前文分析,理論上使用TSWT一次濾波后可以將面波徹底分離,再對提取出的面波進(jìn)行SSWT,重提取部分低頻有效信號.對于實(shí)際數(shù)據(jù)時頻分析來說,精確定位面波區(qū),徹底分離面波干擾,避免過多損傷非目標(biāo)頻段的有效信號才是主要目的.

    時域同步擠壓算法可以有效解決傳統(tǒng)時頻分析方法能量沿頻率軸發(fā)散的問題.雖然未使用同步擠壓算法的時頻分析方法得到的時頻譜看起來能量成團(tuán)分布,但是相鄰能量出現(xiàn)多處重疊,這導(dǎo)致以往的時頻域?yàn)V波方法為了徹底分離面波,存在擴(kuò)大濾波窗口(系數(shù)充零區(qū))的現(xiàn)象,會損傷非目標(biāo)頻段的有效信號,而TSWT會有效改善這一問題,但仍未徹底解決.因此時頻窗口選擇范圍越精確,信號的保幅性越好,雖然這會犧牲一定的計算效率.

    本文所提出的SSWT重提取在一次濾波中被濾除的有效信號,充分利用了反射波和面波的視速度差異.結(jié)合圖5c和圖9b的重構(gòu)結(jié)果,是聯(lián)合方法中SSWT可以恢復(fù)部分被損傷的低頻有效信號的直接證據(jù),這一優(yōu)勢是單一方法和其他聯(lián)合面波壓制方法不具備的.圖14的四種方法處理結(jié)果對比再次證實(shí)了上述觀點(diǎn).圖14a中可見基于視速度差異的FK濾波對面波壓制效果較為徹底但是會損傷淺層有效信號如直達(dá)波,同時會產(chǎn)生大量空間假頻;圖14b和圖14c對比證明了TSWT的面波壓制能力和保幅能力優(yōu)于CWT但還是丟失了部分中深層有效信號;圖14d使用聯(lián)合方法得到的處理結(jié)果面波壓制力和保幅性兼優(yōu),壓制效果最為理想.

    4.2 聯(lián)合方法的缺點(diǎn)及未來研究方向

    如上述分析在使用同步擠壓小波變換后對時頻、空頻濾波窗口內(nèi)系數(shù)充零時沒有尋找到最優(yōu)處理方式.本研究下一步可以改進(jìn)濾波方式,在閾值迭代去噪、自適應(yīng)時窗、深度學(xué)習(xí)等領(lǐng)域進(jìn)行探索.

    在SSWT重構(gòu)有效信號環(huán)節(jié),主要依據(jù)的是面波和有效信號視速度的差異即在地震剖面中同相軸傾斜程度不同.然而由于面波的頻散效應(yīng)部分相鄰道往往會形成低傾角同相軸,這些假層位同相軸振幅往往較小,但會不可避免的出現(xiàn)在重構(gòu)結(jié)果中,如圖5c所示.如何去除這些假層位低傾角同相軸是未來研究的重點(diǎn)問題.

    在圖13中45 Hz的信噪比對比中聯(lián)合方法處理后的地震數(shù)據(jù)信噪比略低于僅使用TSWT處理.我們抽取了三組45 Hz地震數(shù)據(jù)及其時頻能量切片,發(fā)現(xiàn)不存在明顯的面波干擾,如圖15b、e所示,因此判斷是由于重構(gòu)誤差或是信噪比算法誤差所致.聯(lián)合方法處理得到的整炮地震記錄信噪比是最高的未與理論相背.

    5 結(jié)論

    本文將同步擠壓小波變換推廣至空間域,并對其實(shí)現(xiàn)過程進(jìn)行了公式推導(dǎo)與實(shí)際驗(yàn)證.并通過時域與空域同步擠壓小波變換聯(lián)合方法壓制面波干擾,并對理論和實(shí)際資料進(jìn)行處理,調(diào)整最優(yōu)參數(shù),取得了較好的處理效果.該聯(lián)合方法的優(yōu)點(diǎn)在于:

    (1)使用高精度時頻分析方法時域同步擠壓小波變換,結(jié)合時頻域?yàn)V波與空頻域?yàn)V波,對面波壓制徹底.

    (2)對一次濾波后分離的面波干擾再次進(jìn)行空域同步擠壓小波變換,從空間、波數(shù)域中判別有效信息與面波干擾的地震屬性差異,能有效提取部分低頻反射波,切實(shí)保護(hù)了低頻段有效信息,提高了地震數(shù)據(jù)的信噪比.

    (3)本文方法的實(shí)現(xiàn),證明了在面波壓制中進(jìn)行時空域?yàn)V波,從而恢復(fù)低頻段有效信號的方法的優(yōu)越性.

    圖13 野外單炮記錄各有效頻段信噪比對比Fig.13 Comparison of signal-to-noise ratio of each valid frequency band recorded by actualdata

    圖14 四種面波壓制方法實(shí)際資料處理結(jié)果對比(a) FK濾波; (b) CWT; (c) TSWT; (d) TSWT & SSWT.Fig.14 Comparison of actual data processing results of four ground roll suppression methods(a) FK filtering; (b) CWT; (c) TSWT; (d) TSWT & SSWT.

    圖15 實(shí)際記錄45 Hz切片(a) TSWT 45 Hz 切片; (b) SSWT重構(gòu)結(jié)果45 Hz切片; (c) 聯(lián)合方法最終記錄45 Hz切片; (d) TSWT 45 Hz時頻譜; (e) SSWT重構(gòu)結(jié)果 45 Hz時頻譜; (f) 聯(lián)合方法最終記錄45 Hz時頻譜.Fig.15 45 Hz slices of actual record(a) TSWT; (b) Reconstruction result of SSWT; (c) Final record with combination method; (d) The time frequency spectrum of TSWT at 45 Hz; (e) The time frequency spectrum of reconstruction result of SSWT at 45 Hz; (f) The time frequency spectrum of final record with combination method at 45 Hz.

    附錄

    本文計算信噪比的方法為相關(guān)時間移法及其改進(jìn)算法(牛聰?shù)?2006;王紅玲,2007).

    Qi,i+1(ti,i+1)=max{Qi,i+1(τ)},(A1)

    其中相鄰道間最大的互相關(guān)值Qi,i+1(ti,i+1)作為信號的功率譜:

    (A2)

    (A3)

    N為地震道數(shù),使用各道的自相關(guān)值作為地震記錄的總能量,則噪聲的能量為:

    (A4)

    信噪比為:

    (A5)

    當(dāng)對含有面波干擾的原始地震記錄進(jìn)行信噪比估算時,需要將使用聯(lián)合方法濾波后的最終記錄作為有效信號能量Es的輸入記錄.總能量E則用原始地震記錄通過相關(guān)時移法計算.

    猜你喜歡
    面波波數(shù)壓制
    聲場波數(shù)積分截斷波數(shù)自適應(yīng)選取方法
    一種基于SOM神經(jīng)網(wǎng)絡(luò)中藥材分類識別系統(tǒng)
    電子測試(2022年16期)2022-10-17 09:32:26
    gPhone重力儀的面波頻段響應(yīng)實(shí)測研究
    地震研究(2021年1期)2021-04-13 01:04:56
    自適應(yīng)相減和Curvelet變換組合壓制面波
    一種新型無人機(jī)數(shù)據(jù)鏈抗壓制干擾技術(shù)的研究
    空射誘餌在防空壓制電子戰(zhàn)中的應(yīng)用
    無人機(jī)(2018年1期)2018-07-05 09:51:02
    一種舊物品擠壓成型機(jī)
    科技資訊(2016年12期)2016-05-30 05:07:58
    對GPS接收機(jī)帶限高斯噪聲壓制干擾的干擾帶寬選擇分析
    重磁異常解釋的歸一化局部波數(shù)法
    基于聲場波數(shù)譜特征的深度估計方法
    一区二区三区免费毛片| 成人手机av| 高清在线视频一区二区三区| 日日摸夜夜添夜夜爱| 亚州av有码| 最新的欧美精品一区二区| 国产成人精品一,二区| 久久久精品区二区三区| 毛片一级片免费看久久久久| 精品久久国产蜜桃| 日本色播在线视频| 久久国产精品大桥未久av| 夜夜爽夜夜爽视频| 九九爱精品视频在线观看| 王馨瑶露胸无遮挡在线观看| 日韩强制内射视频| 视频区图区小说| freevideosex欧美| www.av在线官网国产| 欧美变态另类bdsm刘玥| 久久久久久久国产电影| 国产成人aa在线观看| 久久99热这里只频精品6学生| 午夜福利网站1000一区二区三区| 日韩在线高清观看一区二区三区| 亚洲精品国产av蜜桃| 中文字幕精品免费在线观看视频 | 欧美精品一区二区免费开放| 在线天堂最新版资源| 少妇的逼水好多| 91精品国产九色| 曰老女人黄片| 王馨瑶露胸无遮挡在线观看| 九色成人免费人妻av| 3wmmmm亚洲av在线观看| 你懂的网址亚洲精品在线观看| 久热这里只有精品99| 99久久人妻综合| 欧美日韩av久久| 毛片一级片免费看久久久久| xxxhd国产人妻xxx| 久久人人爽人人爽人人片va| 日本wwww免费看| 美女大奶头黄色视频| 欧美激情极品国产一区二区三区 | 亚洲五月色婷婷综合| 亚洲成人手机| 制服人妻中文乱码| 精品久久久久久久久av| 欧美精品一区二区免费开放| 22中文网久久字幕| 亚洲人成77777在线视频| 人妻 亚洲 视频| 久久精品国产鲁丝片午夜精品| 99久久精品国产国产毛片| 亚洲av电影在线观看一区二区三区| 亚洲精品色激情综合| 街头女战士在线观看网站| 一级毛片黄色毛片免费观看视频| .国产精品久久| 色5月婷婷丁香| 国产成人freesex在线| 国产淫语在线视频| 自线自在国产av| 亚洲国产欧美日韩在线播放| 午夜激情av网站| 在线 av 中文字幕| 亚洲av免费高清在线观看| 亚洲综合精品二区| 高清不卡的av网站| 亚洲精品久久午夜乱码| 在线观看国产h片| 一级毛片电影观看| 国产男人的电影天堂91| 亚洲av成人精品一二三区| 欧美亚洲 丝袜 人妻 在线| 三级国产精品欧美在线观看| 日本免费在线观看一区| 久久久久久伊人网av| 亚洲激情五月婷婷啪啪| 亚洲精品乱码久久久久久按摩| 老司机影院毛片| 看十八女毛片水多多多| 伦理电影大哥的女人| 午夜福利,免费看| 国产精品人妻久久久久久| 成人午夜精彩视频在线观看| av黄色大香蕉| 亚洲不卡免费看| 妹子高潮喷水视频| 日韩 亚洲 欧美在线| 亚洲av免费高清在线观看| 中文字幕精品免费在线观看视频 | 国产一区二区在线观看av| 国产欧美日韩综合在线一区二区| 成人国产麻豆网| 精品久久久久久久久亚洲| 交换朋友夫妻互换小说| 亚洲成色77777| 99国产精品免费福利视频| 男人添女人高潮全过程视频| 亚洲图色成人| av女优亚洲男人天堂| 亚洲欧美清纯卡通| 婷婷色综合www| 一本大道久久a久久精品| 国产亚洲精品久久久com| 我要看黄色一级片免费的| 精品久久久精品久久久| 精品久久久精品久久久| 免费观看的影片在线观看| 久久99热6这里只有精品| 国产精品嫩草影院av在线观看| 能在线免费看毛片的网站| 99九九在线精品视频| 好男人视频免费观看在线| 男人爽女人下面视频在线观看| 插阴视频在线观看视频| 我要看黄色一级片免费的| 大片免费播放器 马上看| 精品国产国语对白av| 国产免费视频播放在线视频| 综合色丁香网| 一级二级三级毛片免费看| 国产亚洲最大av| 欧美另类一区| 国产精品三级大全| 欧美性感艳星| 人人妻人人添人人爽欧美一区卜| 亚洲精品国产色婷婷电影| 久久久精品94久久精品| a级毛片免费高清观看在线播放| 国产男人的电影天堂91| 啦啦啦在线观看免费高清www| 色哟哟·www| 黑人猛操日本美女一级片| 人人澡人人妻人| 观看美女的网站| av国产精品久久久久影院| 国产高清国产精品国产三级| 久久久久人妻精品一区果冻| 亚洲,一卡二卡三卡| a 毛片基地| 亚洲av免费高清在线观看| 国产精品久久久久久精品电影小说| 色婷婷久久久亚洲欧美| 丝袜脚勾引网站| 亚洲精品成人av观看孕妇| 国产日韩欧美视频二区| 国产日韩欧美视频二区| 国产亚洲一区二区精品| 亚洲成人一二三区av| 丰满饥渴人妻一区二区三| 亚洲av综合色区一区| 久久精品久久久久久久性| 性色av一级| 特大巨黑吊av在线直播| 午夜精品国产一区二区电影| 九九爱精品视频在线观看| 国产又色又爽无遮挡免| 精品酒店卫生间| 亚洲久久久国产精品| 夜夜看夜夜爽夜夜摸| 精品国产露脸久久av麻豆| 日本免费在线观看一区| 男的添女的下面高潮视频| 亚洲精品自拍成人| 制服诱惑二区| 又大又黄又爽视频免费| 午夜激情久久久久久久| 免费久久久久久久精品成人欧美视频 | 爱豆传媒免费全集在线观看| 日韩中文字幕视频在线看片| 免费av中文字幕在线| av一本久久久久| 午夜老司机福利剧场| 午夜久久久在线观看| av国产精品久久久久影院| 婷婷色综合www| 国产一级毛片在线| 午夜免费男女啪啪视频观看| 人妻 亚洲 视频| 看免费成人av毛片| 欧美xxⅹ黑人| 免费人妻精品一区二区三区视频| 国产精品秋霞免费鲁丝片| 久久久国产一区二区| 一级二级三级毛片免费看| 王馨瑶露胸无遮挡在线观看| 热99国产精品久久久久久7| 韩国av在线不卡| 高清在线视频一区二区三区| 久久久久久久亚洲中文字幕| 国产色爽女视频免费观看| 亚洲国产精品999| 大陆偷拍与自拍| 亚洲av国产av综合av卡| 3wmmmm亚洲av在线观看| 一本色道久久久久久精品综合| 成人影院久久| 一本久久精品| 建设人人有责人人尽责人人享有的| 91成人精品电影| 在线精品无人区一区二区三| 国产成人精品无人区| freevideosex欧美| 少妇被粗大猛烈的视频| 又粗又硬又长又爽又黄的视频| 久久人人爽av亚洲精品天堂| 国产成人aa在线观看| 特大巨黑吊av在线直播| 免费少妇av软件| 伦理电影大哥的女人| 纵有疾风起免费观看全集完整版| 2022亚洲国产成人精品| 69精品国产乱码久久久| 亚洲国产精品专区欧美| 一级a做视频免费观看| 欧美日韩视频高清一区二区三区二| 久久 成人 亚洲| 久久国内精品自在自线图片| 国产 精品1| 亚洲精品456在线播放app| 国产乱人偷精品视频| 免费人成在线观看视频色| 蜜桃在线观看..| 日韩熟女老妇一区二区性免费视频| 久久人妻熟女aⅴ| 国产精品国产三级专区第一集| 观看av在线不卡| 一本色道久久久久久精品综合| 亚洲精品国产色婷婷电影| videossex国产| 秋霞在线观看毛片| 美女cb高潮喷水在线观看| 亚洲国产色片| 午夜福利视频精品| 大香蕉97超碰在线| 国产黄色视频一区二区在线观看| 亚洲国产成人一精品久久久| 久久久国产欧美日韩av| 欧美老熟妇乱子伦牲交| 国产熟女欧美一区二区| 18在线观看网站| 少妇丰满av| 一级片'在线观看视频| 综合色丁香网| 日本欧美国产在线视频| 丰满饥渴人妻一区二区三| 十八禁高潮呻吟视频| 中国美白少妇内射xxxbb| 黑丝袜美女国产一区| 国产精品 国内视频| 草草在线视频免费看| 九色成人免费人妻av| 91久久精品国产一区二区成人| 五月天丁香电影| 亚洲av综合色区一区| 一个人免费看片子| 高清视频免费观看一区二区| 久久人妻熟女aⅴ| 观看av在线不卡| 久久久国产一区二区| 久久 成人 亚洲| 精品一区二区三区视频在线| 女的被弄到高潮叫床怎么办| 国产色婷婷99| 国产成人精品福利久久| 午夜日本视频在线| 日韩大片免费观看网站| 一级毛片电影观看| 精品久久久精品久久久| 热re99久久精品国产66热6| 亚洲国产毛片av蜜桃av| 不卡视频在线观看欧美| 午夜福利网站1000一区二区三区| 日韩伦理黄色片| 成人18禁高潮啪啪吃奶动态图 | 一本久久精品| 热re99久久精品国产66热6| 日韩视频在线欧美| av有码第一页| 最新的欧美精品一区二区| 久久久久视频综合| 亚洲精品av麻豆狂野| 在线免费观看不下载黄p国产| 老女人水多毛片| 亚洲av福利一区| 黄色怎么调成土黄色| 曰老女人黄片| 国产在线一区二区三区精| 久久av网站| 欧美xxⅹ黑人| 搡女人真爽免费视频火全软件| 青春草亚洲视频在线观看| 亚洲av免费高清在线观看| 午夜精品国产一区二区电影| 熟女av电影| 美女大奶头黄色视频| 亚洲av.av天堂| 91精品国产国语对白视频| 国产精品国产三级国产专区5o| 国产色爽女视频免费观看| av.在线天堂| 99热网站在线观看| 欧美日韩国产mv在线观看视频| 在线亚洲精品国产二区图片欧美 | 亚洲国产av新网站| 天美传媒精品一区二区| 久久久久国产精品人妻一区二区| 亚洲av成人精品一区久久| av.在线天堂| www.色视频.com| 日日摸夜夜添夜夜爱| 免费黄色在线免费观看| 99久久精品一区二区三区| 永久免费av网站大全| 自拍欧美九色日韩亚洲蝌蚪91| 丰满少妇做爰视频| 久久精品国产亚洲av天美| 99国产精品免费福利视频| 亚洲精品日韩在线中文字幕| 精品久久久久久久久av| av女优亚洲男人天堂| 高清在线视频一区二区三区| 精品人妻熟女毛片av久久网站| 亚洲av国产av综合av卡| 国产精品女同一区二区软件| 日韩一区二区三区影片| 亚洲人成网站在线观看播放| 两个人的视频大全免费| 国产男女内射视频| 久久国产亚洲av麻豆专区| 亚洲精品一区蜜桃| 人人妻人人澡人人爽人人夜夜| 日韩中字成人| 免费大片18禁| 少妇 在线观看| 秋霞伦理黄片| 熟女电影av网| 黑人高潮一二区| 91精品国产国语对白视频| 韩国高清视频一区二区三区| 午夜免费男女啪啪视频观看| a级毛片在线看网站| 欧美精品国产亚洲| 色婷婷久久久亚洲欧美| 亚洲av成人精品一区久久| 欧美最新免费一区二区三区| 69精品国产乱码久久久| 亚洲精品av麻豆狂野| 九九爱精品视频在线观看| 国产精品一区二区三区四区免费观看| 国产极品天堂在线| 在线观看国产h片| 少妇 在线观看| 日韩一区二区视频免费看| 精品熟女少妇av免费看| 最近最新中文字幕免费大全7| 一区二区三区乱码不卡18| 免费黄频网站在线观看国产| 国产高清不卡午夜福利| 99热网站在线观看| 日本爱情动作片www.在线观看| 国产一级毛片在线| 久久精品国产亚洲av涩爱| 五月伊人婷婷丁香| 一区二区三区四区激情视频| 午夜精品国产一区二区电影| 国产一区二区三区综合在线观看 | 国产成人午夜福利电影在线观看| 国产永久视频网站| 日韩av免费高清视频| 热99国产精品久久久久久7| 午夜影院在线不卡| 人成视频在线观看免费观看| 黑人猛操日本美女一级片| 亚洲精品成人av观看孕妇| 色5月婷婷丁香| av在线观看视频网站免费| 黄片播放在线免费| 狠狠精品人妻久久久久久综合| 日本黄色日本黄色录像| 欧美日韩成人在线一区二区| 在线观看美女被高潮喷水网站| 韩国高清视频一区二区三区| 国产成人freesex在线| 亚洲国产av影院在线观看| 精品少妇内射三级| 少妇精品久久久久久久| 我的老师免费观看完整版| 久久久精品94久久精品| 丝袜在线中文字幕| 肉色欧美久久久久久久蜜桃| 大片电影免费在线观看免费| 熟女人妻精品中文字幕| 欧美老熟妇乱子伦牲交| 美女主播在线视频| 亚洲av免费高清在线观看| 熟女av电影| 一级,二级,三级黄色视频| 久久青草综合色| 性色av一级| 久久久久久久久久成人| freevideosex欧美| 久久97久久精品| 精品亚洲成a人片在线观看| 永久免费av网站大全| 一级片'在线观看视频| 久久精品国产亚洲网站| 久久鲁丝午夜福利片| 一级毛片 在线播放| 99热6这里只有精品| 伊人久久精品亚洲午夜| 精品国产乱码久久久久久小说| 一二三四中文在线观看免费高清| 内地一区二区视频在线| 国产精品熟女久久久久浪| av在线老鸭窝| 亚洲美女黄色视频免费看| av电影中文网址| 久久av网站| 亚州av有码| 蜜桃国产av成人99| 日本黄大片高清| 久久精品久久精品一区二区三区| 黄色欧美视频在线观看| 亚洲色图 男人天堂 中文字幕 | 搡女人真爽免费视频火全软件| 亚洲欧美一区二区三区国产| 久久精品国产亚洲av涩爱| 亚洲欧美成人综合另类久久久| 大香蕉97超碰在线| 亚洲欧美成人精品一区二区| 久久久久精品久久久久真实原创| 在线 av 中文字幕| 人妻一区二区av| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产精品成人久久小说| 伊人久久精品亚洲午夜| 色吧在线观看| 麻豆乱淫一区二区| 波野结衣二区三区在线| 国产男女内射视频| 婷婷成人精品国产| 日本免费在线观看一区| 精品一区在线观看国产| 中文字幕免费在线视频6| 亚洲第一区二区三区不卡| 日韩强制内射视频| 99热这里只有精品一区| 成人毛片60女人毛片免费| 丰满迷人的少妇在线观看| 高清毛片免费看| 美女大奶头黄色视频| 日本爱情动作片www.在线观看| 国产日韩欧美亚洲二区| 亚洲av中文av极速乱| 国产精品久久久久久av不卡| 成年女人在线观看亚洲视频| 亚洲色图 男人天堂 中文字幕 | 欧美日韩在线观看h| 久久午夜福利片| 成人午夜精彩视频在线观看| 日韩精品有码人妻一区| 亚洲国产色片| 国产成人精品福利久久| 一级片'在线观看视频| 国产精品久久久久久久久免| 亚洲欧美精品自产自拍| 免费大片18禁| 天天影视国产精品| 哪个播放器可以免费观看大片| av.在线天堂| 日本与韩国留学比较| 99久国产av精品国产电影| 十八禁网站网址无遮挡| 大码成人一级视频| 少妇人妻精品综合一区二区| 肉色欧美久久久久久久蜜桃| 黄色怎么调成土黄色| 男女无遮挡免费网站观看| 日韩不卡一区二区三区视频在线| 少妇熟女欧美另类| 国产av码专区亚洲av| 欧美性感艳星| 午夜福利视频在线观看免费| 亚洲av中文av极速乱| 大片免费播放器 马上看| 91精品一卡2卡3卡4卡| 精品久久蜜臀av无| 国产一级毛片在线| 视频在线观看一区二区三区| 亚洲综合精品二区| 国产在线免费精品| 亚洲国产精品一区三区| 亚洲精品久久成人aⅴ小说 | 成人国产av品久久久| 国产精品久久久久久久电影| 黄色怎么调成土黄色| 两个人免费观看高清视频| 精品午夜福利在线看| 妹子高潮喷水视频| a级毛片在线看网站| 成年人午夜在线观看视频| 99精国产麻豆久久婷婷| 日韩成人av中文字幕在线观看| 日韩电影二区| 日本vs欧美在线观看视频| 婷婷色av中文字幕| av网站免费在线观看视频| 最后的刺客免费高清国语| 亚洲国产日韩一区二区| 亚洲情色 制服丝袜| 考比视频在线观看| 日本wwww免费看| 亚洲av综合色区一区| 日韩一本色道免费dvd| 一本久久精品| www.av在线官网国产| av一本久久久久| 在线观看一区二区三区激情| 99热国产这里只有精品6| 人妻系列 视频| 国产免费一级a男人的天堂| 精品少妇内射三级| 在线精品无人区一区二区三| 国产精品三级大全| 老女人水多毛片| 国产精品国产三级国产av玫瑰| 中文字幕人妻熟人妻熟丝袜美| 日韩伦理黄色片| 亚洲av不卡在线观看| 免费高清在线观看视频在线观看| 国产片特级美女逼逼视频| 色网站视频免费| 亚洲,一卡二卡三卡| 久久av网站| 97在线视频观看| 欧美一级a爱片免费观看看| 亚洲国产av影院在线观看| 久久久久久人妻| av福利片在线| 国产精品久久久久久久久免| 亚洲av国产av综合av卡| 人妻少妇偷人精品九色| 边亲边吃奶的免费视频| 国产精品熟女久久久久浪| 久久人妻熟女aⅴ| 中文字幕精品免费在线观看视频 | 蜜桃在线观看..| 中文字幕最新亚洲高清| 国产又色又爽无遮挡免| 免费久久久久久久精品成人欧美视频 | 国产精品国产三级专区第一集| 亚洲av免费高清在线观看| 国产有黄有色有爽视频| av国产久精品久网站免费入址| 亚洲国产精品999| 人人澡人人妻人| 亚洲无线观看免费| 久久午夜综合久久蜜桃| 午夜免费观看性视频| 五月玫瑰六月丁香| 国产免费又黄又爽又色| 日韩中文字幕视频在线看片| 丁香六月天网| 国产成人免费观看mmmm| 午夜福利视频在线观看免费| 免费不卡的大黄色大毛片视频在线观看| 国产在线一区二区三区精| 亚洲国产av新网站| videos熟女内射| 一本色道久久久久久精品综合| 国产亚洲精品第一综合不卡 | 免费高清在线观看视频在线观看| 亚洲欧美中文字幕日韩二区| 全区人妻精品视频| 亚洲欧美清纯卡通| 99国产综合亚洲精品| 91精品国产国语对白视频| 大陆偷拍与自拍| 免费观看的影片在线观看| 国产高清不卡午夜福利| 久久国内精品自在自线图片| 精品人妻偷拍中文字幕| 天堂俺去俺来也www色官网| 观看av在线不卡| 久久久久久久久久久丰满| 国产亚洲欧美精品永久| 色5月婷婷丁香| 91精品国产国语对白视频| 国产高清不卡午夜福利| 日韩欧美精品免费久久| 日本黄大片高清| 水蜜桃什么品种好| freevideosex欧美| 久久久精品区二区三区| 久久99蜜桃精品久久| 飞空精品影院首页| 又大又黄又爽视频免费| 亚洲美女搞黄在线观看| 少妇人妻久久综合中文| 涩涩av久久男人的天堂| 中文精品一卡2卡3卡4更新| 成年女人在线观看亚洲视频| 男女高潮啪啪啪动态图| 午夜视频国产福利| 中国国产av一级| 免费大片黄手机在线观看| 亚洲成人手机| 国产精品一区二区在线观看99| 亚洲婷婷狠狠爱综合网| 九九在线视频观看精品| 人妻少妇偷人精品九色|