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

    基于多路徑Radon變換的地震數(shù)據(jù)噪聲壓制和波型分離方法研究

    2022-01-25 07:06:02唐歡歡毛偉建
    地球物理學(xué)報(bào) 2022年1期
    關(guān)鍵詞:同相軸多路徑拋物

    唐歡歡,毛偉建*

    1 中國科學(xué)院精密測(cè)量科學(xué)與技術(shù)創(chuàng)新研究院計(jì)算與勘探地球物理研究中心,武漢 430077 2 大地測(cè)量與地球動(dòng)力學(xué)國家重點(diǎn)實(shí)驗(yàn)室,武漢 430077

    0 引言

    抗噪性能良好的Radon變換被廣泛應(yīng)用于圖像分析和信號(hào)重構(gòu)等方面,它是一種積分變換,在地震勘探數(shù)據(jù)處理領(lǐng)域,根據(jù)積分路徑的不同,Radon變換可分為線型、雙曲型以及拋物型(Thorson,1978;Hampson,1986;Yilmaz,1989).線性Radon變換(Linear Radon,L_Radon)又稱為傾斜疊加,它是將時(shí)空域的同相軸振幅沿著偏移距方向以一定的斜率進(jìn)行疊加運(yùn)算,因此線性同相軸在線性Radon變換域的系數(shù)表現(xiàn)為一個(gè)強(qiáng)能量點(diǎn),對(duì)應(yīng)的縱橫坐標(biāo)分別表示同相軸的時(shí)間截距及速度信息;而隨機(jī)噪聲經(jīng)過線性疊加后在變換域分布無規(guī)律且能量極弱,該特征有利于隨機(jī)噪聲壓制(Turner,1990;鞏向博等,2009;Sabbione et al.,2015)、面波壓制(李衛(wèi)忠等,1998;Wang,1999;張智等,2017)以及波場(chǎng)分離(Moon et al.,1986;王維紅等,2006;李志娜等,2014;Lin and Sacchi,2020).但在油氣勘探地震數(shù)據(jù)中,有效信號(hào)為反射波,其時(shí)距曲線為雙曲型,經(jīng)過線性Radon變換后,其變換系數(shù)分布并不集中,失去了變換的稀疏性特征且物理意義不明顯.Hampson(1987)將雙曲Radon變換(Hyperbolic Radon,H_Radon)引入到地震數(shù)據(jù)重建中,但雙曲Radon變換在頻率域不解耦,只能在時(shí)間域進(jìn)行計(jì)算導(dǎo)致效率較低(Spitzer et al.,2001;石穎和王維紅,2012).雙曲型同相軸經(jīng)過部分動(dòng)校正后可近似為拋物型,而拋物Radon變換(Parabolic Radon,P_Radon)在頻率域是解耦的,且在頻率域具有共軛對(duì)稱性使其計(jì)算效率較高,因此拋物Radon變換被廣泛地應(yīng)用于多次波壓制(Foster and Mosher,1992;熊登等,2009;鞏向博等,2014;劉菲菲等,2017)及插值重建(Sacchi and Ulrych,1995).Sacchi和Ulrych (1995)提出的利用貝葉斯參數(shù)反演求解拋物Radon變換系數(shù)方法提高了拋物Radon變換的分辨率,有利于提高多次波壓制(Zhang and Wang,2006;Song et al.,2015;Sun et al.,2019)和數(shù)據(jù)重建效果(劉喜武等,2004;Wang et al.,2019).Johansen(1995)利用二階正交多項(xiàng)式對(duì)疊后水平同相軸地震數(shù)據(jù)進(jìn)行擬合,正交多項(xiàng)式系數(shù)可以表征地震數(shù)據(jù)振幅隨偏移距變化情況,該方法缺少同相軸走時(shí)特征的相關(guān)參數(shù),因此該方法針對(duì)的是水平同相軸數(shù)據(jù),根據(jù)此特征,薛亞茹等(2014)將正交多項(xiàng)式引入到Radon變換中,得到具有保幅特征的2D高階Radon變換;唐歡歡和毛偉健(2014)、唐歡歡等(2018,2020)將其推廣到3D,提出基于3D高階Radon變換的地震數(shù)據(jù)保幅重建方法.在改善積分路徑方面,牛濱華(2001)提出的多項(xiàng)式Radon變換,將積分路徑從線型或拋物型變?yōu)?階多項(xiàng)式曲線,提高了淺層雙曲型同相軸數(shù)據(jù)在Radon域的稀疏程度,該變換的優(yōu)勢(shì)在于不增加計(jì)算量的情況下可以對(duì)復(fù)雜形態(tài)同相軸數(shù)據(jù)進(jìn)行稀疏表示,但不能同時(shí)對(duì)線性、拋物型及復(fù)雜形態(tài)同相軸數(shù)據(jù)進(jìn)行稀疏表示;Trad等(2001)提出了將雙曲Radon變換算子和線性Radon變換算子相結(jié)合的混合Radon變換,該方法缺少這兩種算子的交叉項(xiàng)無法對(duì)復(fù)雜同相軸進(jìn)行精確擬合,且在時(shí)間域進(jìn)行計(jì)算效率較低;魯娥和李慶春(2013)將其應(yīng)用到地震數(shù)據(jù)去噪中來,但該變換與積分路徑相關(guān)的變換參數(shù)只有線性項(xiàng),導(dǎo)致在變換域出現(xiàn)混疊現(xiàn)象.

    以上國內(nèi)外學(xué)者分別從積分路徑的形態(tài)、計(jì)算效率、分辨率以及保幅性方面對(duì)Radon變換進(jìn)行改進(jìn),但現(xiàn)有的Radon變換應(yīng)用于地震數(shù)據(jù)處理中仍然存在一個(gè)問題,即當(dāng)線性同相軸、雙曲型同相軸以及復(fù)雜形態(tài)同相軸同時(shí)存在時(shí),不同形態(tài)同相軸在Radon域系數(shù)不能同時(shí)達(dá)到最佳稀疏狀態(tài),只有同相軸形態(tài)和積分路徑一致的數(shù)據(jù)在Radon域系數(shù)才能集中在一點(diǎn),而其他的系數(shù)則是發(fā)散狀態(tài),在空間上有一定的延展性,這將降低隨機(jī)噪聲壓制效果以及波型分離精度.

    針對(duì)該問題,本文提出了多路徑Radon變換(Multi-path Radon,M_Radon),這里我們把積分路徑僅為線性、雙曲型或拋物型的Radon變換稱為單路徑Radon變換,多路徑Radon變換是指對(duì)每個(gè)同相軸都有線型和拋物型的積分路徑,且具有兩種不同類型的變換參數(shù)來與積分路徑相對(duì)應(yīng),然后利用最小二乘稀疏反演方法對(duì)不同形態(tài)同相軸自適應(yīng)匹配最佳積分路徑,得到在Radon域同時(shí)聚焦的高分辨率系數(shù).該變換有利于提高隨機(jī)噪聲的壓制、面波壓制以及波型分離效果.

    1 方法原理

    在地震數(shù)據(jù)處理中,經(jīng)典拋物Radon正反變換如公式(1a)、(1b)所示(Hampson,1986;Beylkin,1987),是廣義Radon變換的一種特殊形式,公式中d為地震數(shù)據(jù),x為檢波點(diǎn)坐標(biāo),q為拋物型積分路徑的曲率,即速度平方的倒數(shù),m為數(shù)據(jù)在Radon域系數(shù);將積分公式中x2變?yōu)閤,如公式(2a)、(2b)所示,其中p為斜率,即速度的倒數(shù),公式為線性Radon正反變換,也稱為傾斜疊加.公式(1a)、(1b)、(2a)、(2b)分別為

    (1a)

    (1b)

    (2a)

    (2b)

    通過線性Radon或拋物Radon變換,理論上可以將t-x域的線性同相軸或拋物型同相軸變換為Radon域的一個(gè)點(diǎn),如圖1所示,圖1a為一合成的線性同相軸,64個(gè)接收道,道間距為20 m,4 ms時(shí)間采樣率,512個(gè)時(shí)間采樣點(diǎn),合成該數(shù)據(jù)的速度為v=4000 m·s-1,經(jīng)過線性Radon變換后其Radon域系數(shù)如圖1b所示;圖1c為一拋物型同相軸,速度為v=1800 m·s-1,經(jīng)過拋物Radon變換后其Radon域系數(shù)如圖1d所示.從圖中可以看出,t-x域的地震數(shù)據(jù)經(jīng)過合適類型的Radon變換后,其Radon域系數(shù)都各自集中在一個(gè)點(diǎn)附近,從圖1b、d上很容易找到這兩點(diǎn)對(duì)應(yīng)的變換參數(shù)分別為p=0.00025 s·m-1,q=3.08×10-7s2·m-2;再利用公式p=1/v和q=1/v2,可以反向計(jì)算出t-x域地震數(shù)據(jù)線性同相軸和拋物型同相軸的速度分別為v=4000 m·s-1和v=1800 m·s-1,與合成模型數(shù)據(jù)的速度一致,即Radon域系數(shù)分布位置具有明確的物理意義.因此,可以利用Radon變換的稀疏性和明確的物理意義特征進(jìn)行隨機(jī)噪聲壓制、波型分離、多次波壓制以及數(shù)據(jù)重建.

    圖1 單軸數(shù)據(jù)Radon變換效果(a)t-x域線性同相軸;(b)圖(a)數(shù)據(jù)線性Radon變換結(jié)果;(c)t-x域拋物型同相軸;(d)圖(c)數(shù)據(jù)拋物Radon變換結(jié)果.Fig.1 Radon transform result of single-axis data(a)Linear in-phase axis in t-x domain;(b)Linear Radon transform of data in (a);(c)Parabolic in-phase axis in t-x domain;(d)Radon transform of data in (c).

    但是,在實(shí)際地震數(shù)據(jù)中,直達(dá)波、折射波和面波的時(shí)距曲線為線性同相軸,反射波的時(shí)距曲線為雙曲型(經(jīng)部分NMO后可近似為拋物型),在使用線性Radon或拋物Radon變換時(shí),不能使所有波型數(shù)據(jù)在變換域同時(shí)達(dá)到最佳稀疏狀態(tài),如圖2合成數(shù)據(jù)所示,為t-x域地震數(shù)據(jù),包含線型和拋物型同相軸;其線性Radon變換域系數(shù)如圖3a所示,拋物Radon變換域系數(shù)如圖3d所示;可以看出,在使用線性Radon變換時(shí),僅有線性同相軸的Radon域系數(shù)是聚焦、稀疏的,其他拋物型同相軸Radon域系數(shù)并不稀疏;同樣地,使用拋物Radon變換時(shí),稀疏性僅對(duì)拋物型同相軸有效;這是因?yàn)楫?dāng)?shù)卣饠?shù)據(jù)同相軸形態(tài)和積分路徑形態(tài)不一致時(shí),需要一系列的變換參數(shù)和積分路徑組成的基函數(shù)來對(duì)地震數(shù)據(jù)擬合,從而在Radon域剖面上一系列連續(xù)的p值或q值位置都有系數(shù)分布,導(dǎo)致系數(shù)不是稀疏的,且不能從Radon域剖面得到所有同相軸的速度信息.圖3b、e分別是從Radon域系數(shù)反變換回時(shí)空域的數(shù)據(jù);圖3c、f分別為經(jīng)過不同變換后的數(shù)據(jù)與原始數(shù)據(jù)的差剖面,可以看出拋物Radon變換對(duì)線性同相軸有一定的正反變換誤差,特別是在同相軸的起始位置;用公式(3)來計(jì)算這兩種變換的相對(duì)均方誤差:

    圖2 線性和拋物型同相軸同時(shí)存在的地震數(shù)據(jù)Fig.2 Seismic data with linear and parabolic events

    圖3 圖2中地震數(shù)據(jù)Radon變換效果(a)線性Radon正變換系數(shù);(b)線性Radon反變換數(shù)據(jù);(c)線性Radon變換誤差剖面;(d)拋物Radon正變換系數(shù);(e)拋物Radon反變換數(shù)據(jù);(f)拋物Radon變換誤差剖面.Fig.3 Radon transform of seismic data in Fig.2(a)Coefficients of linear Radon forward transform;(b)Data of inverse linear Radon transform;(c)Error of linear Radon transform;(d)Coefficients of parabolic Radon inverse transform;(e)Data of inverse parabolic Radon transform;(f)Error of parabolic Radon inverse transform.

    (3)

    本文提出的多路徑Radon變換能解決該問題,其正反變換公式為

    (4a)

    (4b)

    式中的積分路徑既包含線性項(xiàng)又包含拋物型項(xiàng),線性Radon變換參數(shù)p和拋物Radon變換參數(shù)q同時(shí)存在;因此,對(duì)于含有線型的直達(dá)波、折射、面波及雙曲型反射波(可近似為拋物型)的地震數(shù)據(jù),多路徑Radon變換可以看作線性Radon基函數(shù)和拋物型Radon基函數(shù)的組合.選擇p=0基函數(shù),公式(4)可以對(duì)拋物型同相軸進(jìn)行稀疏表達(dá);選擇q=0的基函數(shù),公式(4)可以對(duì)線性同相軸進(jìn)行稀疏表達(dá);而p和q都不為0的基函數(shù),表示公式(4)的積分路徑為多項(xiàng)式,可以對(duì)時(shí)距曲線不是標(biāo)準(zhǔn)的線型或拋物型同相軸進(jìn)行擬合.

    和傳統(tǒng)Radon變換一樣,多路徑Radon變換也將在頻率域通過反演方法來求解Radon域系數(shù)m,對(duì)公式(4a)、(4b)作Fourier正變換,即得到頻率域多路徑Radon變換公式為

    (5a)

    (5b)

    根據(jù)指數(shù)函數(shù)的性質(zhì),公式(5)可以進(jìn)一步分解為

    (6a)

    (6b)

    (7)

    (8)

    (9)

    (10)

    可以通過加權(quán)迭代方法(Sacchi,1997)進(jìn)行近似求解.

    (11)

    (12)

    (13)

    2 模型數(shù)據(jù)測(cè)試

    2.1 多路徑Radon變換的稀疏表示能力分析

    用多路徑Radon變換對(duì)圖2中的數(shù)據(jù)作正變換,線性積分路徑參數(shù)p的最小值p0=0 s·m-1,采樣間隔dp=2×10-5s·m-1,采樣個(gè)數(shù)np=12;拋物型積分路徑參數(shù)q的最小值q0=0 s2·m-2,采樣間隔dq=2×10-8s2·m-2,采樣個(gè)數(shù)nq=21,3次迭代求解多路徑Radon域系數(shù);系數(shù)分布如圖4所示,圖中有12列,分別表示積分路徑參數(shù)取值為(p0,[q0,q1,…,q20]),(p1,[q0,q1,…,q20]),…,(p11,[q0,q1,…,q20])時(shí)的多路徑Radon域系數(shù),縱坐標(biāo)表示時(shí)間,橫坐標(biāo)表示拋物型積分路徑參數(shù).從圖4中可以看出,與傳統(tǒng)Radon變換不同(圖3a、d),多路徑Radon域系數(shù)可以同時(shí)聚焦,且獨(dú)立分布在兩個(gè)不同區(qū)域上,即圖4中紅色方框標(biāo)識(shí)出的第1列及第11列區(qū)域.第一個(gè)區(qū)域是線性積分路徑參數(shù)p=0×dp區(qū),此時(shí)p=0,表示積分路徑為拋物型,為q剖面,該區(qū)域的Radon系數(shù)有4個(gè),分別對(duì)應(yīng)了4個(gè)拋物型積分參數(shù),能夠與圖2中時(shí)空域四個(gè)拋物型同相軸的起始時(shí)間、速度一一對(duì)應(yīng);第二個(gè)區(qū)域是線性積分路徑參數(shù)p=10×dp區(qū),該區(qū)域的Radon系數(shù)有1個(gè),對(duì)應(yīng)的拋物型積分路徑參數(shù)q=0×dq,表示積分路徑為線性,為p剖面,該位置的Radon系數(shù)與圖2中時(shí)空域線性同相軸起始時(shí)間、速度一一對(duì)應(yīng),且速度可由線性積分路徑參數(shù)p=10×dp計(jì)算得到;因此,通過多路徑Radon不僅可以使時(shí)空域不同形態(tài)的同相軸在Radon域系數(shù)同時(shí)達(dá)到稀疏狀態(tài),而且能夠很容易地識(shí)別這些系數(shù)對(duì)應(yīng)的時(shí)空域同相軸形態(tài)、時(shí)間截距以及速度.這些優(yōu)勢(shì)有助于提高利用Radon變換方法進(jìn)行隨機(jī)噪聲壓制、面波壓制以及波型分離等技術(shù)的處理效果.將多路徑Radon域系數(shù)通過反變換,我們即可得到時(shí)空域地震數(shù)據(jù),如圖5a所示,變換誤差剖面如圖5b所示,利用公式(3)計(jì)算得到相對(duì)均方誤差為1.36×10-4,誤差極小,這是該變換能夠用于地震數(shù)據(jù)處理的基礎(chǔ).

    圖4 多路徑Radon正變換系數(shù)Fig.4 Coefficients in multi-path Radon forward transform

    圖5 多路徑Radon反變換結(jié)果(a)反變換結(jié)果;(b)誤差剖面.Fig.5 Result of multi-path Radon inverse transform(a)Result of inverse transform;(b)Error.

    圖6 三種變換的稀疏表示能力Fig.6 Sparse representation capability of three Radon transforms

    2.2 多路徑Radon變換在模型數(shù)據(jù)隨機(jī)噪聲壓制及波型分離中的應(yīng)用

    多路徑Radon變換的強(qiáng)稀疏表示能力使其在隨機(jī)噪聲壓制和波型分離問題中具有天然優(yōu)勢(shì),先以模型數(shù)據(jù)為例進(jìn)行展示.圖7a為干凈的原始數(shù)據(jù),包含待分離的線性和拋物型兩個(gè)同相軸,同相軸在遠(yuǎn)偏移距形態(tài)接近且有重疊.地震數(shù)據(jù)共有64道,道間距為20 m,每道有300個(gè)時(shí)間采樣點(diǎn),采樣間隔為4 ms.將采取常用的傅里葉變換(FK)、傳統(tǒng)的線性Radon、拋物Radon以及本文提出的多路徑Radon變換方法對(duì)其進(jìn)行分離.圖7b、c分別為該數(shù)據(jù)在線性Radon變換域和拋物Radon變換域的系數(shù)分布圖,可以看出在線性Radon變換中僅有線性同相軸的系數(shù)分布是稀疏、聚焦的,如圖7b中白色箭頭所示,而拋物性同相軸的系數(shù)分布在空間上具有延展性,使二者的系數(shù)有部分重疊,很難將其完全分離開來;類似的情況也出現(xiàn)在拋物Radon變換系數(shù)中,如圖7c所示;圖7d、e、f分別為線性同相軸的FK譜、拋物同相軸的FK譜以及二者混合一起后的FK譜,很明顯這兩個(gè)軸的FK譜能量差異極大,重疊在一起后的FK譜中這兩個(gè)軸的能量也有部分重疊,同樣也不能將其完全分離開;圖7g為多路徑Radon域系數(shù),可以看出線性同相軸和拋物同相軸的系數(shù)分布在不同區(qū)域且同時(shí)稀疏、聚焦,很容易將其分離.

    圖7 時(shí)空域數(shù)據(jù)在不同變換域系數(shù)(a)原始數(shù)據(jù);(b)線性Radon變換域系數(shù);(c)拋物Radon變換域系數(shù);(d)線性軸FK譜;(e)拋物軸FK譜;(f)原始數(shù)據(jù)FK譜;(g)多路徑Radon變換域系數(shù).Fig.7 Coefficients of temporal-spatial data in different transform domains(a)Original data;(b)Coefficients in linear Radon transform domain;(c)Coefficients in parabolic Radon transform domain;(d)FK spectrum of linear event;(e)FK spectrum of parabolic event;(f)FK spectrum of original data;(g)Coefficients of multi-path Radon transform domain.

    表1 四種變換波型分離效果比較Table 1 Comparison of data separation by 4 transforms

    為了測(cè)試本文方法抗噪能力,在圖7a原始干凈數(shù)據(jù)上添加了一定的隨機(jī)噪聲,數(shù)據(jù)的信噪比為-16.09 dB,如圖9a所示.根據(jù)上文三種變換的稀疏表示能力曲線來確定噪聲壓制的閾值,首先將圖7變換域中前30%的強(qiáng)能量系數(shù)保留,其他系數(shù)進(jìn)行壓制,再進(jìn)行下一步的分離.我們比較了本文提出的多路徑Radon變換與傳統(tǒng)線性Radon、拋物Radon變換方法的去噪效果.圖9b、c分別為數(shù)據(jù)在線性Radon、拋物Radon變換域系數(shù),可以看出,由于隨機(jī)噪聲的存在,使得線性同相軸和拋物同相軸的變換域系數(shù)更加難以區(qū)分,進(jìn)而影響下一步的分離;圖9d為數(shù)據(jù)在多路徑Radon變換域系數(shù),可以看出,即使存在能量較強(qiáng)的隨機(jī)噪聲,由于不同類型同相軸在該變換域系數(shù)都同時(shí)稀疏、聚焦,從而仍然能夠較容易將其區(qū)分并分離;線性Radon變換方法的噪聲壓制及波型分離結(jié)果如圖8e、f所示,拋物Radon變換方法的噪聲壓制及波型分離結(jié)果如圖9g、h所示,本文提出的多路徑Radon變換方法的噪聲壓制及波型分離結(jié)果如圖9i、j所示,從圖中可以很明顯看出,線性Radon變換僅對(duì)線性同相軸的噪聲壓制效果較好,拋物Radon變換僅對(duì)拋物同相軸的噪聲壓制有效,而多路徑Radon變換方法對(duì)線性同相軸和拋物同相軸的分離及噪聲壓制同時(shí)有效,同相軸分離完全、無振幅損失且噪聲壓制干凈.三種方法噪聲壓制及分離后數(shù)據(jù)的信噪比如表2所示,本文提出的多路徑Radon變換方法在去噪及分離后,信噪比達(dá)到28 dB左右,與原始-16.09 dB含噪數(shù)據(jù)相比,信噪比提高了44 dB,其效果明顯優(yōu)于線性Radon和拋物Radon變換.

    圖8 波型分離效果(a)(b)傅里葉變換波型分離后的數(shù)據(jù);(c)(d)線性Radon變換波型分離后的數(shù)據(jù);(e)(f)拋物Radon變換波型分離后的數(shù)據(jù);(g)(h)多路徑Radon變換波型分離后的數(shù)據(jù).Fig.8 Waveform separations by different transforms(a)and (b)Fourier transform;(c)and (d)Linear Radon transform;(e)and (f)Parabolic Radon transform;(g)and (h)Multi-path Radon transform.

    圖9 去噪及分離效果(a)含噪數(shù)據(jù);(b)線性Radon變換系數(shù);(c)拋物Radon變換系數(shù);(d)多路徑Radon變換系數(shù);(e)(f)線性Radon變換去噪及分離結(jié)果;(g)(h)拋物Radon變換去噪及分離結(jié)果;(i)(j)多路徑Radon變換去噪及分離結(jié)果.Fig.9 Denoised and separated data by multi-path Radon transform(a)Noised data;(b)Coefficients of linear Radon transform;(c)Coefficients of parabolic Radon transform;(d)Coefficients of multi-path Radon transform;(e)and (f)Denoised and separated data by linear Radon transform;(g)and (h)Denoised and separated data by parabolic Radon transform;(i)and (j)Denoised and separated data by multi-path Radon transform.

    表2 含噪數(shù)據(jù)噪聲壓制及波型分離效果比較Table 2 Comparison of denoised and separated data by three Radon transforms

    2.3 多路徑Radon變換在實(shí)際數(shù)據(jù)噪聲壓制及波型分離中的應(yīng)用

    現(xiàn)將該方法應(yīng)用于塔里木某區(qū)域采集的實(shí)際數(shù)據(jù)面波壓制測(cè)試中,如圖10a所示,共71道,道間距為30 m,每道625個(gè)時(shí)間采樣,采樣間隔為4 ms,強(qiáng)烈的面波表現(xiàn)為線性,將近偏移距的有效信號(hào)湮沒.根據(jù)上文模型數(shù)據(jù)的測(cè)試結(jié)果,可以發(fā)現(xiàn)在線性和拋物型同相軸同時(shí)存在的地震數(shù)據(jù)中,線性Radon變換的處理效果優(yōu)于拋物Radon變換,因此在實(shí)際數(shù)據(jù)面波壓制測(cè)試中,將僅提供線性Radon和本文提出的多路徑Radon變換的處理效果對(duì)比.二者面波壓制效果分別如圖10b、c所示,從圖中可以看出這兩種方法都能夠?qū)⒕€性面波壓制,但多路徑Radon變換方法能夠得到分辨率更高的近偏移距反射波同相軸,如圖中A白色橢圓區(qū)域所示;同時(shí),對(duì)于深層較弱信號(hào),如2.4 s深度B白色橢圓區(qū)域所示,本文方法恢復(fù)的同相軸更為連續(xù);從去除的面波數(shù)據(jù)圖10e、f中紅色箭頭處,可以進(jìn)一步發(fā)現(xiàn)線性Radon變換方法去除了更多有效信號(hào).

    圖10 面波壓制效果(a)原始數(shù)據(jù);(b)線性Radon變換面波壓制結(jié)果;(c)多路徑Radon變換面波壓制結(jié)果;(d)留白區(qū)域;(e)(f)分別為(b)(c)去除的面波.Fig.10 Suppression of surface waves by different transforms(a)Original data;(b)Linear Radon transform;(c)Multi-path Radon transform;(d)Whitespace;(e)and (f)Surface waves removed by (b)and (c),respectively.

    接著將多路徑Radon變換方法應(yīng)用于VSP數(shù)據(jù)的上下行波分離中,原始VSP數(shù)據(jù)如圖11所示,該數(shù)據(jù)共201道,道間距為20 m,每道1905個(gè)時(shí)間采樣,采樣間隔為4 ms,該數(shù)據(jù)較為復(fù)雜,每個(gè)同相軸都不是標(biāo)準(zhǔn)的線性或雙曲型.圖12為三種Radon變換方法的上下行波分離效果:圖12a、b分別為線性Radon變換方法分離的下行波和上行波;圖12c、d分別為拋物Radon變換方法分離的下行波和上行波;圖12e、f分別為多路徑Radon變換方法分離的下行波和上行波.可以看出,多路徑Radon變換方法的分離效果最佳,能夠?qū)崿F(xiàn)上下行波的有效分離;線性Radon變換方法次之,其振幅存在一定的損失,且在上下行波的起始交叉點(diǎn)位置有部分假頻,如圖12a、b中紅色箭頭位置所示;拋物Radon變換方法的效果最差,在近偏移距處傾角較大的同相軸沒有得到有效恢復(fù),且存在拋物型腳印,如圖12c、d中紅色虛框處所示;圖13為近偏移距0~1000 m數(shù)據(jù)分離效果的放大圖,可以更明顯地看出多路徑Radon變換方法上下行波分離的更為徹底且無信號(hào)損失,該實(shí)例進(jìn)一步顯示了本文方法的有效性.

    圖11 原始VSP數(shù)據(jù)Fig.11 Original VSP data

    圖13 VSP數(shù)據(jù)上下行波分離近偏移距放大效果(a)(b)線性Radon變換近偏移距分離效果;(c)(d)拋物Radon變換近偏移距分離效果;(e)(f)多路徑 Radon變換近偏移距分離效果.Fig.13 Enlargement of separation of near-offset upward and downward waves in VSP data (a)and (b)Linear Radon transform;(c)and (d)Parabolic Radon transform;(e)and (f)Multi-path Radon transform.

    3 結(jié)論

    本文提出的基于最小二乘稀疏反演多路徑Radon變換方法,能夠使不同形態(tài)同相軸的Radon域系數(shù)自動(dòng)分離,且保持系數(shù)的稀疏性.該方法解決了在地震數(shù)據(jù)中同時(shí)存在線性和拋物型同相軸時(shí),傳統(tǒng)Radon變換域系數(shù)不能同時(shí)準(zhǔn)確地進(jìn)行稀疏表達(dá)的問題.本文從變換的稀疏程度、計(jì)算效率以及處理效果方面做了詳細(xì)的分析,與傳統(tǒng)線性Radon、拋物Radon相比,在不增加額外計(jì)算量的情況下,多路徑Radon變換能夠更好地識(shí)別不同系數(shù)所對(duì)應(yīng)的同相軸形態(tài)、時(shí)間截距以及速度.這些優(yōu)勢(shì)使Radon變換在隨機(jī)噪聲壓制、面波壓制以及波型分離中獲得更佳的處理效果.

    需要注意的是,在地質(zhì)結(jié)構(gòu)特別復(fù)雜區(qū)域,地震剖面中的同相軸不滿足線性及雙曲型走時(shí)特征時(shí),需要采樣間隔較小的p參數(shù)和q參數(shù)對(duì)其進(jìn)行擬合,即p和q的采樣個(gè)數(shù)增加,使多路徑Radon變換算子規(guī)模成倍加大,從而導(dǎo)致該變換的計(jì)算效率降低,且變換域系數(shù)不再集中于某些點(diǎn)而是在一定范圍內(nèi)進(jìn)行展布,失去稀疏特征.因此,綜合考慮處理效果和計(jì)算效率,多路徑Radon變換在同相軸有明顯線性和雙曲型形態(tài)的地震數(shù)據(jù)中更具優(yōu)勢(shì).

    致謝感謝評(píng)審專家提出的中肯意見和建議,使論文得以完善和豐富.

    猜你喜歡
    同相軸多路徑拋物
    高空拋物罪的實(shí)踐擴(kuò)張與目的限縮
    法律方法(2022年2期)2022-10-20 06:45:28
    多路徑效應(yīng)對(duì)GPS多普勒測(cè)速的影響
    虛同相軸方法及其在陸上地震層間多次波壓制中的應(yīng)用
    關(guān)于拋物-拋物Keller-Segel類模型的全局解和漸近性
    基于5.8G射頻的多路徑識(shí)別技術(shù)應(yīng)用探討
    不要高空拋物!
    高空莫拋物
    一種改進(jìn)的相關(guān)法自動(dòng)拾取同相軸
    一種反射同相軸自動(dòng)拾取算法
    基于5.8GHz多路徑精確識(shí)別方案研究
    免费在线观看视频国产中文字幕亚洲| 人成视频在线观看免费观看| 免费看日本二区| 啦啦啦观看免费观看视频高清| 欧美精品亚洲一区二区| 一级a爱片免费观看的视频| 琪琪午夜伦伦电影理论片6080| 成人手机av| 天堂影院成人在线观看| 久久久久免费精品人妻一区二区 | 国产成人系列免费观看| 99re在线观看精品视频| 91av网站免费观看| 欧美成人免费av一区二区三区| 91在线观看av| 欧美又色又爽又黄视频| 欧美乱妇无乱码| 欧洲精品卡2卡3卡4卡5卡区| av免费在线观看网站| 久久国产精品男人的天堂亚洲| 国产精品av久久久久免费| 在线观看免费午夜福利视频| 99国产精品99久久久久| 身体一侧抽搐| 国产午夜精品久久久久久| 午夜福利视频1000在线观看| 亚洲人成网站高清观看| 12—13女人毛片做爰片一| 88av欧美| 国产精品美女特级片免费视频播放器 | 满18在线观看网站| e午夜精品久久久久久久| 真人做人爱边吃奶动态| 午夜免费成人在线视频| 欧美激情高清一区二区三区| 精品日产1卡2卡| 无限看片的www在线观看| 精品人妻1区二区| 久久久久久久午夜电影| 欧美色欧美亚洲另类二区| 中文字幕高清在线视频| 亚洲全国av大片| 久久久水蜜桃国产精品网| 制服诱惑二区| 中文字幕高清在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 免费人成视频x8x8入口观看| 欧洲精品卡2卡3卡4卡5卡区| 精品国产亚洲在线| 国产三级黄色录像| e午夜精品久久久久久久| 韩国av一区二区三区四区| 男人舔奶头视频| 久久久久久久午夜电影| 成年免费大片在线观看| 欧美乱码精品一区二区三区| 色在线成人网| 亚洲aⅴ乱码一区二区在线播放 | 亚洲熟女毛片儿| 久久天堂一区二区三区四区| 国产欧美日韩一区二区精品| 可以在线观看毛片的网站| 一区二区三区高清视频在线| 亚洲无线在线观看| 首页视频小说图片口味搜索| 国产精品国产高清国产av| 精华霜和精华液先用哪个| 国产伦人伦偷精品视频| 成熟少妇高潮喷水视频| 国产精品98久久久久久宅男小说| 亚洲av美国av| 欧美成狂野欧美在线观看| 欧美久久黑人一区二区| 欧美又色又爽又黄视频| 国产成人精品无人区| 亚洲avbb在线观看| 999久久久精品免费观看国产| 免费在线观看完整版高清| 久久久久九九精品影院| 精品国内亚洲2022精品成人| 婷婷精品国产亚洲av在线| 黄色视频,在线免费观看| 国产亚洲欧美98| 中文资源天堂在线| 欧美黑人精品巨大| 麻豆国产av国片精品| 欧美最黄视频在线播放免费| 国产野战对白在线观看| 日本 欧美在线| 最近最新中文字幕大全电影3 | 一级片免费观看大全| 国内久久婷婷六月综合欲色啪| 国产精品98久久久久久宅男小说| 99在线人妻在线中文字幕| 在线观看免费午夜福利视频| 搡老熟女国产l中国老女人| 91成人精品电影| 一夜夜www| 色av中文字幕| 久久香蕉精品热| 国产精品二区激情视频| 97超级碰碰碰精品色视频在线观看| 女人爽到高潮嗷嗷叫在线视频| 无遮挡黄片免费观看| 亚洲欧美一区二区三区黑人| 91大片在线观看| 亚洲一区二区三区色噜噜| 亚洲精品美女久久av网站| 午夜老司机福利片| 曰老女人黄片| 国产成年人精品一区二区| 亚洲成av人片免费观看| 在线观看舔阴道视频| 久久人人精品亚洲av| 满18在线观看网站| 国产人伦9x9x在线观看| 12—13女人毛片做爰片一| 一个人免费在线观看的高清视频| 中出人妻视频一区二区| 国产精品亚洲一级av第二区| 中文字幕人妻丝袜一区二区| 国产国语露脸激情在线看| 搡老熟女国产l中国老女人| 国产真人三级小视频在线观看| 最近最新免费中文字幕在线| 真人做人爱边吃奶动态| 国产主播在线观看一区二区| 精品熟女少妇八av免费久了| 两性午夜刺激爽爽歪歪视频在线观看 | 国产一卡二卡三卡精品| 国产日本99.免费观看| 色婷婷久久久亚洲欧美| 老汉色av国产亚洲站长工具| 成人亚洲精品av一区二区| 亚洲专区国产一区二区| 久久久国产精品麻豆| 首页视频小说图片口味搜索| 国产99白浆流出| 午夜福利免费观看在线| 久久伊人香网站| 一区二区三区激情视频| 免费看美女性在线毛片视频| 日本一本二区三区精品| 亚洲欧美日韩无卡精品| 国产av又大| 精品卡一卡二卡四卡免费| 亚洲第一av免费看| 久9热在线精品视频| 91在线观看av| 亚洲激情在线av| av福利片在线| 国产精品98久久久久久宅男小说| a级毛片a级免费在线| 最新在线观看一区二区三区| 中文资源天堂在线| 国产一卡二卡三卡精品| 国产成人av教育| 日韩国内少妇激情av| 色综合亚洲欧美另类图片| 99在线视频只有这里精品首页| 欧美乱码精品一区二区三区| 国产成年人精品一区二区| 老司机午夜十八禁免费视频| 一区二区三区激情视频| 一夜夜www| 久久热在线av| 国产在线观看jvid| 精品久久久久久久末码| 窝窝影院91人妻| 欧美av亚洲av综合av国产av| 99久久无色码亚洲精品果冻| 久久国产亚洲av麻豆专区| 丁香欧美五月| av福利片在线| 亚洲专区字幕在线| 精品久久久久久,| 1024视频免费在线观看| 搡老熟女国产l中国老女人| 久久婷婷人人爽人人干人人爱| 国产av又大| 一本大道久久a久久精品| 国产精品久久久久久精品电影 | 97碰自拍视频| 草草在线视频免费看| 亚洲全国av大片| 国产精品国产高清国产av| 黄网站色视频无遮挡免费观看| 在线视频色国产色| 亚洲精品中文字幕一二三四区| 国产久久久一区二区三区| 一进一出抽搐gif免费好疼| 男女之事视频高清在线观看| 国产精品电影一区二区三区| 国产成人系列免费观看| 日韩有码中文字幕| 国产亚洲av高清不卡| 亚洲无线在线观看| 夜夜躁狠狠躁天天躁| 日韩欧美 国产精品| 777久久人妻少妇嫩草av网站| 日韩 欧美 亚洲 中文字幕| 午夜福利视频1000在线观看| 99国产精品一区二区三区| 黄片小视频在线播放| 久久久久久亚洲精品国产蜜桃av| 99久久久亚洲精品蜜臀av| 免费观看人在逋| 亚洲免费av在线视频| 亚洲专区国产一区二区| 18禁美女被吸乳视频| 精品一区二区三区av网在线观看| 亚洲美女黄片视频| 国产熟女午夜一区二区三区| 国产精品免费视频内射| 香蕉丝袜av| 亚洲欧洲精品一区二区精品久久久| 成人国产综合亚洲| 亚洲精品美女久久久久99蜜臀| 亚洲成人久久爱视频| 国产99久久九九免费精品| 成人三级黄色视频| a在线观看视频网站| 一级毛片高清免费大全| 亚洲第一青青草原| 中文字幕精品亚洲无线码一区 | 亚洲自拍偷在线| 午夜亚洲福利在线播放| 久久久久精品国产欧美久久久| 嫩草影视91久久| 69av精品久久久久久| 怎么达到女性高潮| 久久久久久九九精品二区国产 | 午夜精品在线福利| 国产免费男女视频| 露出奶头的视频| √禁漫天堂资源中文www| 悠悠久久av| 免费人成视频x8x8入口观看| 精品卡一卡二卡四卡免费| 免费电影在线观看免费观看| 亚洲av电影不卡..在线观看| 男人舔奶头视频| 男人的好看免费观看在线视频 | 每晚都被弄得嗷嗷叫到高潮| 亚洲七黄色美女视频| 成人欧美大片| 深夜精品福利| 色综合欧美亚洲国产小说| 日韩免费av在线播放| 一边摸一边抽搐一进一小说| 欧美激情久久久久久爽电影| 国产精品av久久久久免费| 夜夜躁狠狠躁天天躁| av欧美777| 精品一区二区三区四区五区乱码| 老司机靠b影院| 精品久久久久久久末码| 激情在线观看视频在线高清| 免费女性裸体啪啪无遮挡网站| 脱女人内裤的视频| 在线免费观看的www视频| 精品久久蜜臀av无| 亚洲美女黄片视频| xxxwww97欧美| 自线自在国产av| a在线观看视频网站| www日本在线高清视频| 色综合站精品国产| 怎么达到女性高潮| 亚洲久久久国产精品| 国产久久久一区二区三区| 欧美日本亚洲视频在线播放| 热99re8久久精品国产| 国产视频一区二区在线看| 欧美性猛交╳xxx乱大交人| 国产成人av激情在线播放| 国产又色又爽无遮挡免费看| 中出人妻视频一区二区| 成人三级黄色视频| 国产单亲对白刺激| 亚洲自拍偷在线| 可以在线观看的亚洲视频| 日本成人三级电影网站| av超薄肉色丝袜交足视频| 无限看片的www在线观看| 国产三级黄色录像| 中文字幕av电影在线播放| 精品国产美女av久久久久小说| 少妇的丰满在线观看| 91字幕亚洲| 午夜福利视频1000在线观看| 亚洲中文字幕日韩| 18禁裸乳无遮挡免费网站照片 | 国产精品一区二区三区四区久久 | 欧美日本亚洲视频在线播放| 午夜激情福利司机影院| 天天躁夜夜躁狠狠躁躁| 1024手机看黄色片| 亚洲一区二区三区色噜噜| 一边摸一边抽搐一进一小说| 亚洲午夜精品一区,二区,三区| av在线天堂中文字幕| 脱女人内裤的视频| 日韩三级视频一区二区三区| 看免费av毛片| 99精品欧美一区二区三区四区| 亚洲一区中文字幕在线| 成人一区二区视频在线观看| 久久国产乱子伦精品免费另类| 三级毛片av免费| 麻豆久久精品国产亚洲av| 亚洲精品久久国产高清桃花| 亚洲专区国产一区二区| 最近最新免费中文字幕在线| 久热爱精品视频在线9| 狠狠狠狠99中文字幕| 香蕉久久夜色| 国产精品av久久久久免费| 一区二区三区国产精品乱码| 国产三级黄色录像| 亚洲中文日韩欧美视频| 久久久久久大精品| 黄色毛片三级朝国网站| 在线天堂中文资源库| 国产在线精品亚洲第一网站| 天堂动漫精品| 欧美一级a爱片免费观看看 | 在线观看日韩欧美| 国产一区二区三区视频了| 18禁裸乳无遮挡免费网站照片 | 老汉色av国产亚洲站长工具| 日本一本二区三区精品| 88av欧美| 看免费av毛片| 美女高潮到喷水免费观看| 国内精品久久久久精免费| 精品久久久久久,| 制服诱惑二区| 久久久久久久久中文| 熟女电影av网| 中文字幕人妻丝袜一区二区| 一个人观看的视频www高清免费观看 | 国产成人精品久久二区二区免费| 国产精品香港三级国产av潘金莲| 麻豆成人午夜福利视频| 国产精品乱码一区二三区的特点| 天堂动漫精品| tocl精华| 日韩高清综合在线| 亚洲人成网站高清观看| 女人高潮潮喷娇喘18禁视频| 亚洲一码二码三码区别大吗| 精品久久久久久成人av| 欧美日韩黄片免| 午夜亚洲福利在线播放| 最新美女视频免费是黄的| 国产精品影院久久| 欧美成人一区二区免费高清观看 | 免费高清在线观看日韩| 在线播放国产精品三级| 成人国产一区最新在线观看| 久久人妻福利社区极品人妻图片| 九色国产91popny在线| 久久国产乱子伦精品免费另类| 大香蕉久久成人网| 久久精品亚洲精品国产色婷小说| 国产精品一区二区三区四区久久 | 亚洲欧美精品综合一区二区三区| 成人亚洲精品一区在线观看| 91麻豆av在线| 久久国产精品男人的天堂亚洲| 婷婷丁香在线五月| 久久亚洲精品不卡| 国产又黄又爽又无遮挡在线| 亚洲国产欧洲综合997久久, | 黄色视频,在线免费观看| 国产亚洲精品av在线| 天天添夜夜摸| 欧美精品啪啪一区二区三区| 成人午夜高清在线视频 | 午夜精品久久久久久毛片777| 999精品在线视频| 后天国语完整版免费观看| 亚洲国产欧美日韩在线播放| 中文字幕精品亚洲无线码一区 | 在线观看日韩欧美| 亚洲精品在线观看二区| 成人国语在线视频| 国产精品香港三级国产av潘金莲| 天堂影院成人在线观看| 久久精品影院6| 日日夜夜操网爽| 亚洲精品国产精品久久久不卡| 免费一级毛片在线播放高清视频| 日韩视频一区二区在线观看| 超碰成人久久| 久久久久久久精品吃奶| 一边摸一边做爽爽视频免费| 麻豆成人午夜福利视频| 日本五十路高清| 黄色片一级片一级黄色片| 人人妻,人人澡人人爽秒播| 日韩一卡2卡3卡4卡2021年| 人人妻人人澡欧美一区二区| 特大巨黑吊av在线直播 | 国产精品电影一区二区三区| 嫩草影院精品99| 一本综合久久免费| 国产成人精品无人区| av免费在线观看网站| 俺也久久电影网| 身体一侧抽搐| 很黄的视频免费| 在线视频色国产色| 亚洲第一电影网av| 国产精品二区激情视频| 亚洲中文日韩欧美视频| 精品免费久久久久久久清纯| 国产精品免费视频内射| 好男人在线观看高清免费视频 | 欧美日韩亚洲国产一区二区在线观看| 一进一出好大好爽视频| 国产精品电影一区二区三区| 久久人妻福利社区极品人妻图片| 国产黄片美女视频| 午夜福利视频1000在线观看| 伊人久久大香线蕉亚洲五| 99久久精品国产亚洲精品| 一级a爱片免费观看的视频| 午夜久久久在线观看| 91麻豆精品激情在线观看国产| 国产aⅴ精品一区二区三区波| 长腿黑丝高跟| 久久精品国产综合久久久| 亚洲精品中文字幕一二三四区| 国语自产精品视频在线第100页| 久久久久九九精品影院| 91字幕亚洲| 国产成人欧美在线观看| 国产免费av片在线观看野外av| 日本一本二区三区精品| 亚洲av电影在线进入| 亚洲va日本ⅴa欧美va伊人久久| av视频在线观看入口| 一级a爱视频在线免费观看| 桃红色精品国产亚洲av| 身体一侧抽搐| 久久九九热精品免费| 精品少妇一区二区三区视频日本电影| 精品乱码久久久久久99久播| 法律面前人人平等表现在哪些方面| 一区福利在线观看| 我的亚洲天堂| 极品教师在线免费播放| 黄色成人免费大全| 一级作爱视频免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 午夜久久久久精精品| 国产黄a三级三级三级人| 91在线观看av| 日本 欧美在线| 国产成人欧美| 首页视频小说图片口味搜索| 色综合亚洲欧美另类图片| 欧美最黄视频在线播放免费| 免费女性裸体啪啪无遮挡网站| 黄色成人免费大全| 操出白浆在线播放| 日韩 欧美 亚洲 中文字幕| 国产又黄又爽又无遮挡在线| 露出奶头的视频| 又黄又爽又免费观看的视频| 欧美成狂野欧美在线观看| 人人妻,人人澡人人爽秒播| 校园春色视频在线观看| 亚洲av成人av| 搞女人的毛片| 男人的好看免费观看在线视频 | 欧美日韩亚洲综合一区二区三区_| 国内毛片毛片毛片毛片毛片| 欧美精品亚洲一区二区| 在线观看免费午夜福利视频| 欧美在线黄色| 精品一区二区三区av网在线观看| 欧美日韩亚洲综合一区二区三区_| 国产精品一区二区精品视频观看| 国产欧美日韩一区二区精品| 亚洲国产精品sss在线观看| 国内久久婷婷六月综合欲色啪| 久久久水蜜桃国产精品网| 好看av亚洲va欧美ⅴa在| 国产精品久久久久久人妻精品电影| 成人特级黄色片久久久久久久| 久久久精品国产亚洲av高清涩受| 国产精品香港三级国产av潘金莲| 久久青草综合色| 免费看美女性在线毛片视频| 亚洲人成网站在线播放欧美日韩| 很黄的视频免费| 婷婷六月久久综合丁香| 国产精品电影一区二区三区| 满18在线观看网站| 黄色 视频免费看| 免费看日本二区| 日韩欧美三级三区| 女人爽到高潮嗷嗷叫在线视频| 麻豆av在线久日| 自线自在国产av| 国产精品久久电影中文字幕| 变态另类丝袜制服| 美女国产高潮福利片在线看| 久久香蕉激情| 精品国产一区二区三区四区第35| 亚洲在线自拍视频| 一夜夜www| 欧美日韩精品网址| 国产国语露脸激情在线看| 午夜a级毛片| 日日爽夜夜爽网站| 人成视频在线观看免费观看| 午夜福利高清视频| 亚洲avbb在线观看| 久久久久久久久久黄片| 最近最新中文字幕大全电影3 | 在线观看www视频免费| 亚洲av第一区精品v没综合| 中文字幕精品免费在线观看视频| 99国产综合亚洲精品| 香蕉久久夜色| 精品国产乱码久久久久久男人| 无限看片的www在线观看| 搡老妇女老女人老熟妇| 欧美性长视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | 久久精品91蜜桃| 国产av又大| 男人舔女人的私密视频| 亚洲人成电影免费在线| 日韩免费av在线播放| 精品一区二区三区四区五区乱码| 波多野结衣高清无吗| xxxwww97欧美| 搡老岳熟女国产| 99热这里只有精品一区 | 在线av久久热| 亚洲国产精品999在线| 欧美另类亚洲清纯唯美| 99在线人妻在线中文字幕| 激情在线观看视频在线高清| 99久久国产精品久久久| 香蕉av资源在线| 久久久久免费精品人妻一区二区 | 欧美精品亚洲一区二区| 很黄的视频免费| 国产精品免费一区二区三区在线| 中文字幕精品免费在线观看视频| 欧美色视频一区免费| 美女午夜性视频免费| 免费在线观看日本一区| 一个人观看的视频www高清免费观看 | 成人国语在线视频| 亚洲熟妇中文字幕五十中出| 757午夜福利合集在线观看| 久99久视频精品免费| 国内揄拍国产精品人妻在线 | 亚洲精品美女久久av网站| 欧美不卡视频在线免费观看 | 黑人操中国人逼视频| 久久久久九九精品影院| 美女大奶头视频| 中文字幕av电影在线播放| 亚洲精华国产精华精| 国产男靠女视频免费网站| 91老司机精品| 国产精品久久久久久人妻精品电影| 日韩精品青青久久久久久| 亚洲天堂国产精品一区在线| 男女之事视频高清在线观看| 久久久久久九九精品二区国产 | 亚洲 欧美 日韩 在线 免费| 精品久久蜜臀av无| 女性生殖器流出的白浆| 久久国产精品影院| 亚洲一区高清亚洲精品| 啦啦啦 在线观看视频| 一级a爱视频在线免费观看| 成年免费大片在线观看| 长腿黑丝高跟| 一二三四在线观看免费中文在| 亚洲成a人片在线一区二区| 十八禁网站免费在线| 国语自产精品视频在线第100页| 97超级碰碰碰精品色视频在线观看| 亚洲久久久国产精品| 国产乱人伦免费视频| 桃色一区二区三区在线观看| 欧美黑人欧美精品刺激| 欧美日韩精品网址| 欧美丝袜亚洲另类 | 久久精品国产清高在天天线| 国产伦在线观看视频一区| 中文字幕人妻熟女乱码| 黄色女人牲交| 日韩国内少妇激情av| 国产成人av教育| 欧美黑人欧美精品刺激| 国产乱人伦免费视频| 国产97色在线日韩免费| 99在线人妻在线中文字幕| 麻豆成人av在线观看| 亚洲欧美精品综合久久99| 一个人免费在线观看的高清视频| 国产一级毛片七仙女欲春2 | 国产亚洲欧美在线一区二区| 一边摸一边做爽爽视频免费|