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

    基于一種三維低冗余曲波變換和壓縮感知理論的地震數(shù)據(jù)重建

    2017-11-22 01:22:19曹靜杰王尚旭李文斌
    關(guān)鍵詞:曲波冗余度范數(shù)

    曹靜杰, 王尚旭, 李文斌

    (1.中國石油大學(xué)地球物理與信息工程學(xué)院,北京 102249; 2.河北地質(zhì)大學(xué)勘查技術(shù)與工程學(xué)院,河北石家莊 050031;3.河北地質(zhì)大學(xué)信息工程學(xué)院,河北石家莊 050031)

    基于一種三維低冗余曲波變換和壓縮感知理論的地震數(shù)據(jù)重建

    曹靜杰1,2, 王尚旭1, 李文斌3

    (1.中國石油大學(xué)地球物理與信息工程學(xué)院,北京 102249; 2.河北地質(zhì)大學(xué)勘查技術(shù)與工程學(xué)院,河北石家莊 050031;3.河北地質(zhì)大學(xué)信息工程學(xué)院,河北石家莊 050031)

    基于稀疏變換的重建方法是地震數(shù)據(jù)重建的研究熱點,其中稀疏變換的性質(zhì)對重建的效率和質(zhì)量起關(guān)鍵作用。曲波變換對波場數(shù)據(jù)有非常稀疏的表示和可靠的數(shù)值效果,然而三維曲波變換的冗余度在24~32之間,須消耗很多的內(nèi)存和計算時間。提出基于一種低冗余度曲波變換的地震數(shù)據(jù)重建方法,介紹低冗余曲波變換并且分析其優(yōu)點,給出一個解基于分析的1范數(shù)模型的快速迭代閾值方法。數(shù)據(jù)實驗結(jié)果表明:該低冗余變換將三維曲波變換的冗余度降低了60%,能夠有效提高數(shù)據(jù)重建的計算效率;基于低冗余度曲波變換的數(shù)據(jù)重建的計算效率是基于原始曲波變換重建效率的4倍,對于10%的采樣比例仍然能夠得到較好的重建和去噪效果。

    曲波變換; 地震重建; 稀疏優(yōu)化; 1范數(shù)

    由于受到地表障礙、壞道、禁區(qū)和成本等因素的制約,野外地震數(shù)據(jù)在空間上經(jīng)常不滿足采樣定理,不完備的地震數(shù)據(jù)會影響波動方程偏移[1]、表面多次波消除[2]和AVO分析[3]等數(shù)據(jù)處理的效果,因此須重建地震數(shù)據(jù)以獲得完備的數(shù)據(jù)?;谙∈枳儞Q的重建方法計算效率高,數(shù)值效果穩(wěn)定,因此是重建方法重要的研究方向。該方法利用信號在變換域的稀疏性為約束建立反演模型,并且采用稀疏反演方法求解。稀疏變換的性質(zhì)是決定這類方法質(zhì)量的首要因素。地震信號在變換域中的表達越稀疏,則越有利于提高重建質(zhì)量。在地震信號處理領(lǐng)域經(jīng)典的稀疏變換有傅里葉變換[4]和Radon類變換[5-6]。最近曲波變換[7]等多尺度變換被用于地震數(shù)據(jù)重建且展示了很好的數(shù)值效果。多尺度變換比單一尺度的變換能夠展示更多的信息,因此是信號分析的有力工具。曲波變換打破了傳統(tǒng)傅里葉變換和小波變換只能表示單一尺度或者不能稀疏表示多尺度信號的缺點。對于地震數(shù)據(jù)有非常稀疏的表達,在地震數(shù)據(jù)重建[8-10]、去噪[11-12]、多次波消除[13-14]、反褶積[15]、偏移[16]和反射數(shù)據(jù)的增強[17]等問題中都有應(yīng)用并且展示了很好的數(shù)值效果。曲波變換對曲線或曲面形狀信號進行多方向和各項異性的表示,為邊緣不連續(xù)的目標(biāo)提供了幾乎最優(yōu)的表達[18-21]。雖然曲波變換有多尺度性、多方向性以及無須對數(shù)據(jù)劃分時間窗口等優(yōu)點,但是曲波變換的冗余度(冗余度指的是變換域系數(shù)的維數(shù)與原始信號系數(shù)維數(shù)之比)很大。在二維情況下的原始曲波變換的冗余度約為8,三維數(shù)據(jù)變換的冗余度為 24~32,需要很多的內(nèi)存和計算時間,影響大規(guī)模數(shù)據(jù)計算的效率。筆者將一種低冗余三維曲波變換用于三維地震數(shù)據(jù)重建,這種新的低冗余曲波變換[22]比原始曲波變換的冗余度有很大的下降,因此稱之為低冗余度曲波變換,在求解時采用一種解基于分析的1范數(shù)正則化模型的迭代軟閾值算法,對模擬和真實地震數(shù)據(jù)的試驗證明方法的有效性。

    1 三維低冗余度曲波變換

    曲波[23]變換有兩種快速實現(xiàn)方式,分別是基于USFFT(unequally-space fast transform)的快速算法和基于wrapping的快速算法。由于基于wrapping的方法更加簡單和容易執(zhí)行,本文中采用基于wrapping的曲波變換的實現(xiàn)過程?;趙rapping的三維曲波[23]變換包括一個低通逼近子帶和一些包含曲波系數(shù)的曲波基子帶,這些曲波系數(shù)由尺度、位置和方向索引[22]。該變換主要由兩部分組成:多尺度分解和角度分解。首先,假設(shè)輸入的三維數(shù)據(jù)為N=(Nx,Ny,Nz),其中,Nx、Ny、Nz分別表示x、y、z三個方向上的維數(shù),在傅里葉域通過三維Meyer小波將其分解,得到尺度為N,N/2,…,N/2J的以原點為中心的立方體,其中J為尺度的個數(shù)。第二步,將這些立方體和角度窗口相乘,分割成各向異性的楔形體,這些楔形體符合拋物線尺度準(zhǔn)則,將楔形體包裝成立方體后采用三維傅里葉逆變換得到曲波系數(shù)。三維低冗余度曲波變換的具體實現(xiàn)方式如下:①輸入N=(Nx,Ny,Nz)的三維數(shù)據(jù)體,定義劃分的尺度個數(shù)J,在最小尺度的每個面上劃分的角度個數(shù)Na。 ②多尺度分解。在傅里葉域應(yīng)用三維Meyer小波變換,得到尺度為N,N/2,…,N/2J的立方體。③角度分解。在每一個尺度的每個方向上,在傅里葉域?qū)⑿〔⒎襟w和角度窗相乘得到楔形體;將這些楔形體包裝成中心在原點的平行六面體;對平行六面體應(yīng)用三維傅里葉逆變換得到曲波系數(shù)。

    在原始曲波變換中,冗余度主要來源于角度分解。此外,Meyer小波的應(yīng)用方法產(chǎn)生了額外冗余性,而低冗余曲波變換中Meyer小波沒有產(chǎn)生額外冗余性[22]。 新的曲波變換和原始曲波變換在實現(xiàn)上存在許多不同之處,但是其能夠降低額外冗余度的主要原因是Meyer小波變換的應(yīng)用方法不同。

    (1)

    其中ν從0到1逐漸變化并且滿足ν((x)+ν(1-x)=1。Meyer 尺度函數(shù)定義為

    (2)

    最小尺度的小波在傅里葉域的支撐為[-2/3,-1/6]∪[1/6,2/3], 因此超出了Shannon 帶。原始曲波變換隱含的假設(shè)周期邊界條件,計算一個周期信號的小波變換等價于將信號在周期小波基上分解。最小尺度小波超出Shannon 帶的部分被其以ξ=1/2為中心沿著豎軸做鏡像代替,數(shù)據(jù)的支撐是原始數(shù)據(jù)的 4/3倍,對于三維數(shù)據(jù)來說,產(chǎn)生的額外的冗余度為 (4/3)3。在低冗余度曲波變換中,尺度函數(shù)和小波函數(shù)的支撐首先被縮小為原來的3/4,然后為了保持單位劃分的一致性,對最小尺度的小波進行修改來壓制其下降的尾部,小波在[-1/2,-1/4]∪[1/4,1/2]之間變?yōu)橐粋€固定值。這樣對于三維Meyer小波來說,在傅里葉域內(nèi)沒有增加額外的冗余度。

    對于原始曲波變換,三維數(shù)據(jù)的冗余度為 24~32,而低冗余度曲波變換的冗余度約為10[22]。除了在最小尺度上保持方向選擇性和低冗余性,該低冗余曲波變換保持了等距性,并且能夠快速精確重建,即該低冗余度曲波變換是一個Parseval 緊框架,C*C=I,其中C為曲波分解算子,C*為其伴隨。C*也是逆算子,因此該變換具有快速的重構(gòu)算法。關(guān)于低冗余度曲波變換的詳細實現(xiàn)見參考文獻[22]。

    2 基于壓縮感知理論的地震數(shù)據(jù)重建

    地震采集過程是地震信號的采樣過程,因此地震重建問題可以根據(jù)壓縮感知理論求解[7]。壓縮感知理論證明了可以根據(jù)少量的測量數(shù)據(jù)高精度地重建原始數(shù)據(jù),從而降低地震勘探的成本。地震采集過程可以表示為如下數(shù)學(xué)模型:

    Rx+ε=d.

    (3)

    式中,R為采樣矩陣;x為完整的地震數(shù)據(jù);ε為可加性噪聲;d為采樣的數(shù)據(jù)。由于數(shù)據(jù)d是不完備且不滿足采樣定理的,理論上存在無窮多x滿足式(3)。根據(jù)壓縮感知理論,假設(shè)地震數(shù)據(jù)在某個變換域中是稀疏的,即s=Ψx是稀疏的或者可壓縮的,其中Ψ是一個稀疏變換算子,令Ψ為上述的低冗余曲波變換,由于其是緊框架,其共軛轉(zhuǎn)置也是其逆,即Ψ*=Ψ-1,則公式(3)變?yōu)?/p>

    RΨ*s+ε=d.

    (4)

    其中Ψ*為Ψ的共軛轉(zhuǎn)置?;趕的稀疏性假設(shè),可以通過解以下優(yōu)化問題求出s:

    (5)

    其中σ≥0是對噪聲能量的估計,該問題稱為基追蹤問題[24]。 當(dāng)σ=0時,上述模型針對高信噪比數(shù)據(jù)或者不含噪聲數(shù)據(jù)的重建;當(dāng)σ>0時,上述模型可以求解含隨機噪聲數(shù)據(jù)的重建,即可以實現(xiàn)同時重建和去噪,因此上述模型既可以求解不含噪聲數(shù)據(jù)的重建,也可以求解含隨機噪聲數(shù)據(jù)的重建。當(dāng)s被解出后,對于低冗余曲波變換來說,x可以由x=Ψ*s求出。由于問題(5)是一個約束優(yōu)化問題,對于該問題的求解方法很少[10],一般求解其等價的無約束優(yōu)化問題

    (6)

    其中λ是一個正則參數(shù),用來平衡殘差的2范數(shù)和解的1范數(shù),當(dāng)λ選擇合適時,問題(5)和(6)的解相等。為了使問題(6)的解與數(shù)據(jù)模型相匹配,對于不含噪聲數(shù)據(jù)的重建,只須選擇λ很小,接近于0即可;對于含隨機噪聲的重建,若σ已知,則可以通過偏差原則,廣義偏差原則[25]等方法求解與σ匹配的正則參數(shù)λ,若σ已知,則可以采用類似于廣義交叉檢驗原則(GCV)或者L曲線準(zhǔn)則的方法求解合適的λ。由于廣義交叉檢驗原則(GCV)或者L曲線準(zhǔn)則不能保證一定獲得可靠的λ[26],因此λ常常通過人工經(jīng)驗選取來獲得。

    地震數(shù)據(jù)重建屬于大規(guī)模計算問題,因此需要快速穩(wěn)定的計算方法。對于上述1范數(shù)正則化的問題(6)的求解,不能采用經(jīng)典的方法如共軛梯度法或者牛頓法[27-28],傳統(tǒng)的內(nèi)點法等由于用到矩陣求逆,因此計算速度無法與經(jīng)典方法相比,須研究解1范數(shù)正則化模型的快速算法。閾值類方法[29]是針對1范數(shù)正則化模型求解的一類重要方法,該方法根據(jù)1范數(shù)可分的性質(zhì)得到,很多解問題(6)的方法,如不動點法、逼近點法等都含有軟閾值運算。最近研究發(fā)現(xiàn)采用如下的模型的數(shù)值效果比基于模型(6)的效果要好:

    (7)

    該模型稱為基于分析的模型,以時間空間域的數(shù)據(jù)x為變量,問題(6)被稱為基于合成的模型,以變換域中的系數(shù)s為變量。雖然兩者之間可以相互推導(dǎo),但是基于分析的模型(7)的數(shù)值效果要稍好于前者。給出一個解問題(7)的簡單易行的閾值算法,該算法以x為變量進行迭代求解。算法如下:

    (1)輸入采樣的數(shù)據(jù)d,采樣矩陣R,稀疏變換Ψ,最大迭代步數(shù)N,控制停機的參數(shù)ω和最小閾值τmin,令k=0。

    (5)輸出最終解xfinal=xk。

    如果將算法中軟閾值運算變?yōu)橛查撝祷蛘咂渌倪M的閾值運算,上述算法也能得到稀疏解。本文中步長κk固定為1,同時還有更加先進的步長選擇方法如BB法、回溯法等方法,步長的選擇對算法有較大的影響。閾值的下降可以是線性的,也可以是指數(shù)下降,本文中采用指數(shù)下降τk+1=τ0ec(k-1)/(N-1),其中c=ln(τ0/τmin),實際上閾值的下降方式對解的改善不大,最小閾值τmin的大小決定著解的質(zhì)量,因為τmin與噪聲的能量相關(guān)。當(dāng)噪聲很小時,選擇τmin很小,當(dāng)噪聲較大時τmin應(yīng)適當(dāng)調(diào)節(jié)。由于噪聲的能量實際上也是很難估計的,本文中采用嘗試法選擇τmin。本文中算法與固定閾值的閾值類算法不同,實際上是帶有連續(xù)化策略的閾值法,其中令閾值逐漸下降就是優(yōu)化算法中的連續(xù)化策略,這比固定閾值參數(shù)的算法需要更少的迭代次數(shù)和更好的數(shù)值效果,因此本算法實際上由兩層組成,外層控制閾值的下降,內(nèi)層只用1次迭代,已經(jīng)證明本算法比固定最小閾值的方法數(shù)值表現(xiàn)要改進很多[26]。

    3 數(shù)值模擬

    圖1 三維合成數(shù)據(jù)重建實驗Fig.1 3D synthetic data reconstruction experiment

    參數(shù)采樣率50%CandesWoiselle采樣率20%CandesWoiselle采樣率10%CandesWoiselle采樣率5%CandesWoiselle計算時間t/s10661767561907591761082271信噪比/dB203365202013135222171638876601036314116576232

    用一個海洋地震數(shù)據(jù)驗證基于低冗余曲波變換的真實數(shù)據(jù)的重建質(zhì)量。該海洋數(shù)據(jù)為512×128×128,即時間軸上的采樣個數(shù)為512,時間采樣率2 ms,在inline和crossline 方向上的采樣個數(shù)為128,相鄰地震道間隔為25 m。圖2給出海洋數(shù)據(jù)的一個切片,圖3是隨機采樣的數(shù)據(jù),采樣的道數(shù)為總道數(shù)的50%。圖4是采用低冗余曲波變換和迭代軟閾值法的重建結(jié)果,其信噪比為 16.559 6 dB。為了更加清晰的顯示重建的效果,將其中一道原始數(shù)據(jù)和其相應(yīng)的重建道進行對比,如圖5所示,重建數(shù)據(jù)和原始數(shù)據(jù)吻合得很好。

    圖2 真實海洋數(shù)據(jù)Fig.2 Field marine data

    圖3 隨機采樣道數(shù)為50%的海洋數(shù)據(jù)Fig.3 Randomly sampled marine data with 50% traces sampled

    圖4 基于低冗余曲波變換的海洋數(shù)據(jù)重建結(jié)果Fig.4 Reconstructed marine data based on low-redundancy curvelet transform

    圖5 一道原始海洋數(shù)據(jù)與重建數(shù)據(jù)對比Fig.5 Comparison of one original trace and reconstructed trace of marine data

    為了驗證不同采樣比例下該低冗余曲波變換的重建效果,對一個含噪聲的三維海洋數(shù)據(jù)分別采樣20%和10%的地震道進行重建,該試驗同時也驗證了低冗余曲波變換對含噪聲數(shù)據(jù)重建的可靠性。由于計算機內(nèi)存的限制,截取一小部分數(shù)據(jù)進行計算,數(shù)據(jù)為200×200×200,即在時間上有200個采樣點,空間兩個方向上分別有200個采樣,時間采樣率為4 ms,道間距為12.5 m。 圖6(a)為原始數(shù)據(jù)的一個切片顯示,圖6(b)為隨機采樣 20%地震道的顯示,圖6(c)為圖6(b)采用低冗余曲波變換得到的重建結(jié)果。由圖可見,不僅未采樣的地震道被重建出來,而且噪聲的能量得到了壓制。 圖6(d)為原始數(shù)據(jù)和重建結(jié)果的誤差,視覺上有效信號很少。為了測試更低采樣率的數(shù)據(jù)能否得到較好的重建,將采樣率降至10%,即隨機選擇10%的地震道進行重建,如圖7(b)所示。圖7(c)為基于低冗余曲波變換的重建結(jié)果,圖7(d)為原始數(shù)據(jù)和重建數(shù)據(jù)的差。由此可見,圖7(c)的結(jié)果是可以接受的,并且圖7(d)中有效信號很少。為了在頻率中比較,圖8(a)為原始數(shù)據(jù)圖7(a)的一個空間切片的頻譜,圖8(b)為重建數(shù)據(jù)(圖7(c))相同空間切片的頻譜,圖8(c)為圖7(d)相同空間切片的頻譜,可以看出基于低冗余曲波變換的重建可以濾掉高波數(shù)分量和一些高頻分量。本試驗說明:①低冗余曲波變換對于非常低的采樣比例的數(shù)據(jù)也能夠得到穩(wěn)定的結(jié)果,當(dāng)采樣比例為5%時,數(shù)值效果會變差。②通過選擇合適的參數(shù),采用的模型和算法能夠?qū)崿F(xiàn)含噪聲數(shù)據(jù)的重建,即問題(6)和(7)的模型可以解決含隨機噪聲的重建問題。③相比于二維變換,三維變換需要的采樣比例更少。之前的二維數(shù)據(jù)重建研究表明,當(dāng)以曲波變換為稀疏變換,并且采用閾值類算法如凸集投影法[26]求解時,二維數(shù)據(jù)的采樣比例在50%左右才能得到較好的結(jié)果,一旦低于這個比例,則數(shù)值效果會變差,三維變換由于將兩個空間維數(shù)的信息整合起來,具有更稀疏的表達,因此能夠適合更低的采樣比例。同樣,五維數(shù)據(jù)重建也適合更低的采樣比例。

    圖6 三維真實數(shù)據(jù)隨機采樣20%的重建試驗Fig.6 3D field data reconstruction with 20% traces sampled randomly

    圖7 三維真實數(shù)據(jù)隨機采樣10%的重建Fig.7 3D field data reconstruction with 10% traces sampled randomly

    圖8 真實數(shù)據(jù)和重建數(shù)據(jù)的頻譜比較Fig.8 Spectrum comparation of original data and reconstructed data

    4 結(jié)束語

    將一種低冗余度曲波變換用于加速地震數(shù)據(jù)重建,給出了一種快速的迭代閾值法實現(xiàn)地震數(shù)據(jù)重建和去噪,提高了基于曲波變換的地震數(shù)據(jù)重建的計算效率和重建質(zhì)量。數(shù)據(jù)試驗表明,基于低冗余度曲波變換的重建時間約為基于原始曲波變換的1/4,且能夠達到可靠的重建結(jié)果,可以適合更低采樣比例的數(shù)據(jù)重建?;诙喑叨茸儞Q的重建存在一定的局限性,由于曲波變換是局部的、多尺度變換,如果沒有對小尺度構(gòu)造充分的采樣,則基于這種多尺度變換的重建不能得到滿意的效果。

    [1] LIU B, SACCHI M. Minimum weighted norm interpolation of seismic records [J]. Geophysics,2004,69(6):1560-1568.

    [2] NAGHIZADEH M. Parametric reconstruction of multidimensional seismic records [D]. Edmonton:University of Alberta,2009.

    [3] SACCHI M, LIU B. Minimum weighted norm wavefield reconstruction for AVA imaging [J]. Geophysical Prospecting, 2005,53(6):787-801.

    [4] ZWARTJES P, GISOLF A. Fourier reconstruction of marine-streamer data in four spatial coordinates [J]. Geophysics,2006,71(6):V171-V186.

    [5] SACCHI M, VERSCHUUR D, ZWARTJES P. Data reconstruction by generalized deconvolution [C/OL]//2004 SEG Annual Meeting. Denver,Colorado, October 10-15,2004:SEG Technical Program Expanded Abstracts, 2004:1989-1992 [2009-03-05].http://library.seg.org/doi/abs/10.1190/1.1843303.

    [6] 薛亞茹,楊靜,錢步仁. 基于高階稀疏Radon變換的預(yù)測多次波自適應(yīng)相減方法[J]. 中國石油大學(xué)學(xué)報(自然科學(xué)版),2015,39(1):43-49.

    XUE Yaru, YANG Jing, QIAN Buren. Multiples prediction and adaptive subtraction with high-order sparse Radon transform [J]. Journal of China University of Petroleum (Edition of Natural Science),2015,39(1):43-49.

    [7] HERRMANN F, MOGHADDAM P, STOLK C. Sparsity-and continuity-promoting seismic image recovery with curvelet frames[J]. Applied and Computational Harmonic Analysis,2008,24:50-173.

    [8] 曹靜杰,王彥飛,楊長春. 地震數(shù)據(jù)壓縮重構(gòu)的正則化與零范數(shù)稀疏最優(yōu)化方法[J]. 地球物理學(xué)報,2012,55(2):596-607.

    CAO Jingjie, WANG Yanfei, YANG Changchun. Seismic data restoration based on compressive sensing using the regularization and zero-norm sparse optimization[J]. Chinese Journal of Geophysics (in Chinese),2012,55(2):596-607.

    [9] YANG P, GAO J, CHEN W. Curvelet-based POCS interpolation of nonuniformly sampled seismic records[J]. Journal of Applied Geophysics,2012,79(1):90-99.

    [10] 曹靜杰,王本鋒. 基于一種改進凸集投影方法的地震數(shù)據(jù)同時插值和去噪[J]. 地球物理學(xué)報,2015,58(8):2935-2947.

    CAO Jingjie, WANG Benfeng. An improved projection onto convex sets method for simultaneous interpolation and denoising[J]. Chinese Journal of Geophysics (in Chinese),2015,58(8):2935-2947.

    [11] HENNENFENT G, HERRMANN F. Seismic denoising with non-uniformly sampled curvelets [J]. Computing in Science & Engineering,2006,8(3):16-25.

    [12] NEELAMANI R, BAUMSTEIN A, GILLARD D, et al. Coherent and random noise attenuation using the curvelet transform[J]. The Leading Edge,2008,27(2):240-248.

    [13] HERRMANN F, WANG D, HENNENFENT G, et al. Curvelet-based seismic data processing:a multiscale and nonlinear approach[J]. Geophysics, 2008,73(1): A1-A5.

    [14] LIN T, HERRMANN F. Robust estimation of primaries by sparse inversion via one-norm minimization[J]. Geophysics,2013,78(3):R133-R150.

    [15] KUMAR V, HERRMANN F. Deconvolution with curvelet-domain sparsity[C/OL]//2008 SEG Annual Meeting. Las Vegas,Nevada, November 9-14,2008:SEG Technical Program Expanded Abstracts, c2008:1996-2000 [2011-05-07].http://library.seg.org/doi/abs/10.1190/1.3059287.

    [16] CHAURIS H, NGUYEN T. Seismic demigration/migration in the curvelet domain[J]. Geophysics, 2008,73(2):S35-S46.

    [17] KUMAR V, OUEITY J,CLOWES R M, et al. Enhancing crustal reflection data through curvelet denoising[J]. Tectonophysics,2011,508(1):106-116.

    [18] CANDES E, DONOHO D. Curvelets: a surprisingly effective nonadaptive representation for objects with edges[R].Stanford Univ Ca Dept of Statistics, 2000.

    [19] STARKC J, CANDES E, DONOHO D. The curvelet transform for image denoising[J]. IEEE Transaction on Image Processing,2002,11(6):670-684.

    [20] CANDES E, DONOHO D. New tight frames of curvelets and optimal representation of objects with piecewise C2 singularities[J]. Communications on Pure and Applied Mathematics,2004,57(2):219-266.

    [21] CANDES E, DEMANET L, DONOHO, et al. Fast discrete curvelet transforms[J]. Multiscale Modeling & Simulation,2006,5(3):861-899.

    [22] WOISELLE A, STARCK J, FADILI J. 3D Data denoising and inpainting with the low-redundancy fast curvelet transform[J]. Journal of Mathematical Imaging and Vision,2011,39(2):121-139.

    [23] YING L, DEMANET L, CANDES E. curvelet transform[C/OL].//Proc. SPIE wavelets XI, San Diego. 2005, vol. 591413 [2008-9-4]. https://www.spiedigitallibrary.org/conference-proceedings-of-spie/5914/1/3D-discrete-curvelet-transform/10.1117/12.616205.short?SSO=1.

    [24] CHEN S, DONOHO D, SAUNDERS M. Atomic decomposition by basis pursuit[J]. SIAM Journal on Scientific Computing,1998,20(1):33-61.

    [25] 肖庭延,于慎根,王彥飛. 反問題的數(shù)值解法[M]. 北京:科學(xué)出版社,2003.

    [26] YAGOLA A, LEONOV A, TITARENKO V. Data errors and an error estimation for ill-posed problems[J]. Inverse Problems in Science and Engineering,2002,10(2):117-129.

    [27] CAO J, WANG Y, WANG B. Accelerating seismic interpolation with a gradient projection method based on tight frame property of curvelet[J]. Exploration Geophysics,2015,46(3):253-260.

    [28] 王彥飛. 反問題的計算方法及其應(yīng)用[M]. 北京:高等教育出版社,2007.

    [29] CAO J, WANG Y, ZHAO J, et al. A review on restoration of seismic wavefields based on regularization and compressive sensing [J]. Inverse Problems in Science and Engineering,2011,19(5):679-704.

    (編輯 修榮榮)

    Seismicreconstructionwith3Dlow-redundancycurvelettransformandcompressedsensingtheory

    CAO Jingjie1,2, WANG Shangxu1, LI Wenbin3

    (1.CollegeofGeophysicsandInformationEngineeringinChinaUniversityofPetroleum,Beijing102249,China; 2.CollegeofExplorationTechnologyandEngineering,HebeiGEOUniversity,Shijiazhuang050031,China; 3.CollegeofInformationEngineering,HebeiGEOUniversity,Shijiazhuang050031,China)

    Sparse-transform-based seismic data reconstruction is a hot topic in seismic reconstruction, where properties of the sparse transform may influence the results of reconstruction greatly. Curvelet transform is a multi-scale, multi-directional, and local transform which has nearly the sparsest expression for seismic data. However, this transform is a highly redundant transform with redundancy about 24-32 for three dimensional data. To improve the efficiency of curvelet based reconstruction, this paper proposed a low-redundancy curvelet-transform based seismic reconstruction. The new transform was introduced first and its merits for seismic signal processing were analyzed, followed by an iterative thresholding method for analysis-based L1-norm regularized models. Numerical experiments illustrate that the low-redundancy transform can reduce 60% redundancy of the original 3D curvelet transform, thus improves greatly the computational efficiency. The reconstruction computational efficiency based on the low redundancy transform is 4 times of the original curvelet based reconstruction, for example, even for 10% sampling ratio, this low-redundancy curvelet can get acceptable results.

    curvelet transform; seismic reconstruction; sparse optimization; one-norm

    2016-05-23

    國家自然科學(xué)基金項目(41674114);河北省自然科學(xué)基金項目(D2017403027);河北省高校百名優(yōu)秀創(chuàng)新人才支持計劃Ⅲ(SLRC2017024);中國博士后科學(xué)基金資助項目(2016M600171,2017T100137)

    曹靜杰(1982-),男,副教授,博士,碩士生導(dǎo)師,研究方向為地震信號處理和地球物理反演算法。E-mail: cao18601861@163.com。

    1673-5005(2017)05-0061-08

    10.3969/j.issn.1673-5005.2017.05.007

    P 631.4

    A

    曹靜杰,王尚旭,李文斌.基于一種三維低冗余曲波變換和壓縮感知理論的地震數(shù)據(jù)重建[J].中國石油大學(xué)學(xué)報(自然科學(xué)版),2017,41(5):61-68.

    CAO Jingjie, WANG Shangxu, LI Wenbin. Seismic reconstruction with 3D low-redundancy curvelet transform and compressed sensing theory [J]. Journal of China University of Petroleum(Edition of Natural Science), 2017,41(5):61-68.

    猜你喜歡
    曲波冗余度范數(shù)
    一種航天測控冗余跟蹤弧段處理方法
    上海航天(2024年1期)2024-03-08 02:52:28
    林海雪原(五)
    林海雪原(三)
    林海雪原(四)
    曲波變換三維地震數(shù)據(jù)去噪技術(shù)
    上海某基坑工程考慮冗余度的支撐體系設(shè)計
    山西建筑(2017年29期)2017-11-15 02:04:38
    基于加權(quán)核范數(shù)與范數(shù)的魯棒主成分分析
    橋梁設(shè)計的冗余度分析
    矩陣酉不變范數(shù)H?lder不等式及其應(yīng)用
    橋梁設(shè)計的冗余度
    欧美人与性动交α欧美精品济南到| 欧美黑人精品巨大| 欧美日韩精品网址| 在线十欧美十亚洲十日本专区| 脱女人内裤的视频| svipshipincom国产片| 色播亚洲综合网| 欧美三级亚洲精品| 欧美日韩一级在线毛片| 丰满的人妻完整版| 国产成人欧美| 久久欧美精品欧美久久欧美| 免费看a级黄色片| 午夜福利一区二区在线看| 久久久久久免费高清国产稀缺| 国产国语露脸激情在线看| av天堂在线播放| 国产亚洲欧美98| 久久久久久久久久黄片| 久久青草综合色| 亚洲国产看品久久| 国产免费男女视频| 在线观看66精品国产| 老汉色∧v一级毛片| 人人妻,人人澡人人爽秒播| 女人高潮潮喷娇喘18禁视频| 色哟哟哟哟哟哟| 男人舔奶头视频| 麻豆久久精品国产亚洲av| 精品电影一区二区在线| 成人特级黄色片久久久久久久| 91九色精品人成在线观看| 极品教师在线免费播放| 18禁观看日本| 亚洲男人的天堂狠狠| 日韩av在线大香蕉| 久久九九热精品免费| 亚洲免费av在线视频| 久久人妻av系列| 国产亚洲精品一区二区www| 国产亚洲av高清不卡| 一二三四在线观看免费中文在| 国产成人欧美在线观看| 视频区欧美日本亚洲| 麻豆成人av在线观看| 国产亚洲av嫩草精品影院| 久久久久久九九精品二区国产 | 一级a爱片免费观看的视频| 国产真人三级小视频在线观看| 他把我摸到了高潮在线观看| 美女高潮到喷水免费观看| 亚洲成人精品中文字幕电影| 国产精品久久久av美女十八| 亚洲成人精品中文字幕电影| 伊人久久大香线蕉亚洲五| 手机成人av网站| 国产99白浆流出| 国产v大片淫在线免费观看| 中文字幕另类日韩欧美亚洲嫩草| 成人av一区二区三区在线看| 天天躁狠狠躁夜夜躁狠狠躁| 岛国视频午夜一区免费看| 少妇裸体淫交视频免费看高清 | 久久国产乱子伦精品免费另类| 精品高清国产在线一区| bbb黄色大片| 免费看美女性在线毛片视频| 欧美不卡视频在线免费观看 | 51午夜福利影视在线观看| 欧美又色又爽又黄视频| 老司机深夜福利视频在线观看| 日本成人三级电影网站| 韩国精品一区二区三区| 一级a爱片免费观看的视频| 好男人电影高清在线观看| 老司机靠b影院| 男人舔女人的私密视频| 97碰自拍视频| 亚洲三区欧美一区| 亚洲国产日韩欧美精品在线观看 | 国产高清有码在线观看视频 | 国产成年人精品一区二区| 91av网站免费观看| tocl精华| 最近在线观看免费完整版| 国产午夜精品久久久久久| 18禁观看日本| 在线观看舔阴道视频| 国产欧美日韩一区二区精品| 色在线成人网| 亚洲中文日韩欧美视频| 免费在线观看成人毛片| 性色av乱码一区二区三区2| 男女视频在线观看网站免费 | 欧美成狂野欧美在线观看| 欧美zozozo另类| 99国产精品一区二区蜜桃av| www.www免费av| 欧美成狂野欧美在线观看| 在线观看午夜福利视频| 成年人黄色毛片网站| 岛国视频午夜一区免费看| 亚洲精品在线观看二区| 美国免费a级毛片| 成人一区二区视频在线观看| 免费无遮挡裸体视频| 亚洲真实伦在线观看| 十八禁网站免费在线| 大型av网站在线播放| 免费看十八禁软件| 色综合站精品国产| 午夜免费成人在线视频| 日本熟妇午夜| 黄色丝袜av网址大全| 欧美成狂野欧美在线观看| 国产一区二区三区在线臀色熟女| 色在线成人网| 日韩欧美三级三区| 一二三四社区在线视频社区8| 1024视频免费在线观看| 99久久综合精品五月天人人| 精品国产国语对白av| 黑人巨大精品欧美一区二区mp4| 动漫黄色视频在线观看| 婷婷精品国产亚洲av在线| 久久狼人影院| 婷婷六月久久综合丁香| 国产精品影院久久| 身体一侧抽搐| 99国产精品99久久久久| 欧美中文综合在线视频| 亚洲色图 男人天堂 中文字幕| 国产精品av久久久久免费| 亚洲男人天堂网一区| 欧美av亚洲av综合av国产av| 日本熟妇午夜| 国产高清视频在线播放一区| 国产精品,欧美在线| 欧美zozozo另类| 亚洲成国产人片在线观看| 精品久久久久久久人妻蜜臀av| 欧美三级亚洲精品| 免费无遮挡裸体视频| 91九色精品人成在线观看| 国产免费男女视频| 狠狠狠狠99中文字幕| 丰满的人妻完整版| 免费在线观看完整版高清| 亚洲一卡2卡3卡4卡5卡精品中文| 人妻久久中文字幕网| www.精华液| 精品国产乱子伦一区二区三区| 人人妻人人看人人澡| 中文字幕av电影在线播放| 伦理电影免费视频| 久久久久久亚洲精品国产蜜桃av| 大香蕉久久成人网| 一个人免费在线观看的高清视频| 久热这里只有精品99| 日韩欧美免费精品| 免费在线观看亚洲国产| 国产激情偷乱视频一区二区| 老司机靠b影院| 中文字幕人妻丝袜一区二区| 国产成人精品无人区| 亚洲一卡2卡3卡4卡5卡精品中文| 91av网站免费观看| 91麻豆精品激情在线观看国产| 亚洲专区中文字幕在线| 国产视频一区二区在线看| 一区二区三区国产精品乱码| 97超级碰碰碰精品色视频在线观看| 制服诱惑二区| 午夜亚洲福利在线播放| 午夜免费观看网址| 高清毛片免费观看视频网站| 精品久久久久久久久久久久久 | 人人妻人人澡欧美一区二区| 99re在线观看精品视频| 老汉色av国产亚洲站长工具| 免费搜索国产男女视频| 国产熟女午夜一区二区三区| 成人亚洲精品av一区二区| 久久草成人影院| 男男h啪啪无遮挡| 丰满人妻熟妇乱又伦精品不卡| 亚洲成人久久爱视频| 俺也久久电影网| 免费看十八禁软件| 午夜亚洲福利在线播放| 黄片小视频在线播放| 美女免费视频网站| 欧美日韩黄片免| 国产1区2区3区精品| 亚洲avbb在线观看| 高清毛片免费观看视频网站| 在线观看一区二区三区| 身体一侧抽搐| 国产99久久九九免费精品| 麻豆成人av在线观看| 天堂动漫精品| 免费在线观看亚洲国产| 国产成人精品久久二区二区免费| 999久久久国产精品视频| 中文字幕人妻丝袜一区二区| 淫秽高清视频在线观看| 青草久久国产| 亚洲欧美激情综合另类| 一本精品99久久精品77| www.熟女人妻精品国产| 国产精品一区二区精品视频观看| 成年版毛片免费区| 久久久久国产精品人妻aⅴ院| 香蕉久久夜色| 丁香六月欧美| 国产1区2区3区精品| 国产99久久九九免费精品| 制服人妻中文乱码| 人妻久久中文字幕网| 国产视频内射| 午夜老司机福利片| 欧美色视频一区免费| 国产精品久久久久久人妻精品电影| 听说在线观看完整版免费高清| 别揉我奶头~嗯~啊~动态视频| 日日摸夜夜添夜夜添小说| 嫩草影视91久久| 欧美性猛交╳xxx乱大交人| 国产v大片淫在线免费观看| 国产精品1区2区在线观看.| 男女做爰动态图高潮gif福利片| 成人免费观看视频高清| 一进一出抽搐动态| 国产片内射在线| netflix在线观看网站| 精品高清国产在线一区| 亚洲午夜理论影院| 老熟妇乱子伦视频在线观看| 国产精品99久久99久久久不卡| 亚洲一区中文字幕在线| 日韩精品免费视频一区二区三区| 中亚洲国语对白在线视频| 亚洲熟女毛片儿| 国产精品久久久人人做人人爽| 桃红色精品国产亚洲av| 精品国产乱子伦一区二区三区| 午夜福利高清视频| 久久中文字幕人妻熟女| 精品午夜福利视频在线观看一区| 亚洲av五月六月丁香网| 宅男免费午夜| 麻豆国产av国片精品| 一本综合久久免费| 日韩成人在线观看一区二区三区| 国产欧美日韩精品亚洲av| 最近在线观看免费完整版| 91九色精品人成在线观看| 岛国视频午夜一区免费看| 久久草成人影院| 男人舔女人的私密视频| 国产欧美日韩一区二区精品| 久久精品91蜜桃| 两个人视频免费观看高清| 真人做人爱边吃奶动态| 亚洲熟妇中文字幕五十中出| 欧美大码av| 一级a爱片免费观看的视频| 精品熟女少妇八av免费久了| 妹子高潮喷水视频| 国产成人影院久久av| 日本熟妇午夜| 制服丝袜大香蕉在线| 丁香六月欧美| 哪里可以看免费的av片| 国产aⅴ精品一区二区三区波| 久久久国产欧美日韩av| 国产日本99.免费观看| 国产午夜福利久久久久久| av视频在线观看入口| 精品久久久久久久人妻蜜臀av| 国产一卡二卡三卡精品| 一卡2卡三卡四卡精品乱码亚洲| av天堂在线播放| 国产熟女午夜一区二区三区| 亚洲成人免费电影在线观看| 精品久久久久久久末码| 欧美日韩精品网址| 91av网站免费观看| 国产亚洲精品综合一区在线观看 | a级毛片在线看网站| 在线观看舔阴道视频| 久久久水蜜桃国产精品网| 亚洲欧洲精品一区二区精品久久久| 波多野结衣av一区二区av| 精品熟女少妇八av免费久了| 成人av一区二区三区在线看| 后天国语完整版免费观看| 国产精品一区二区精品视频观看| 69av精品久久久久久| 黑人操中国人逼视频| 亚洲国产欧洲综合997久久, | 国产一区二区激情短视频| 天堂影院成人在线观看| 久久天堂一区二区三区四区| 丰满人妻熟妇乱又伦精品不卡| 亚洲第一av免费看| 99久久综合精品五月天人人| 97人妻精品一区二区三区麻豆 | 老汉色∧v一级毛片| 亚洲五月天丁香| 在线十欧美十亚洲十日本专区| 久久久国产精品麻豆| 人妻丰满熟妇av一区二区三区| 丰满的人妻完整版| 国产激情久久老熟女| 丁香欧美五月| 国产97色在线日韩免费| 日本五十路高清| 国产极品粉嫩免费观看在线| 亚洲专区中文字幕在线| 99在线视频只有这里精品首页| 国产亚洲精品久久久久久毛片| 国产精品久久久久久人妻精品电影| 国产伦一二天堂av在线观看| 俄罗斯特黄特色一大片| 好看av亚洲va欧美ⅴa在| 久久久久久国产a免费观看| 亚洲欧美日韩高清在线视频| 午夜免费观看网址| 最近最新中文字幕大全电影3 | 首页视频小说图片口味搜索| 久久久国产欧美日韩av| 国产单亲对白刺激| 国产午夜福利久久久久久| 亚洲熟妇熟女久久| 妹子高潮喷水视频| 久久99热这里只有精品18| 亚洲五月婷婷丁香| svipshipincom国产片| 亚洲一区中文字幕在线| 久久精品aⅴ一区二区三区四区| 18美女黄网站色大片免费观看| 亚洲人成网站高清观看| 天天躁夜夜躁狠狠躁躁| 精品免费久久久久久久清纯| 精品久久久久久成人av| 国产视频一区二区在线看| 亚洲精品美女久久久久99蜜臀| 免费在线观看黄色视频的| 国产成+人综合+亚洲专区| 99re在线观看精品视频| 两性夫妻黄色片| 香蕉久久夜色| 亚洲激情在线av| 国产麻豆成人av免费视频| 国产成+人综合+亚洲专区| 宅男免费午夜| 99热6这里只有精品| 麻豆成人av在线观看| 99国产精品一区二区蜜桃av| 国产精品一区二区精品视频观看| 午夜福利在线观看吧| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产精品成人综合色| 亚洲一区中文字幕在线| 精品久久久久久久末码| 黄网站色视频无遮挡免费观看| 亚洲精品久久国产高清桃花| 亚洲电影在线观看av| 午夜免费激情av| 熟女电影av网| 制服诱惑二区| 999久久久国产精品视频| 日韩有码中文字幕| 午夜福利高清视频| 国产日本99.免费观看| 精品国产亚洲在线| 国内毛片毛片毛片毛片毛片| 99国产精品一区二区蜜桃av| 午夜亚洲福利在线播放| 成人欧美大片| 一区福利在线观看| 99久久国产精品久久久| 亚洲成人免费电影在线观看| 制服丝袜大香蕉在线| 麻豆一二三区av精品| 中文字幕人妻丝袜一区二区| svipshipincom国产片| 日韩成人在线观看一区二区三区| 国产成人一区二区三区免费视频网站| 熟女电影av网| 亚洲av第一区精品v没综合| 精品国产乱子伦一区二区三区| 亚洲美女黄片视频| 高清毛片免费观看视频网站| 欧美 亚洲 国产 日韩一| 男女床上黄色一级片免费看| 男人操女人黄网站| 无限看片的www在线观看| 日韩欧美三级三区| 精品日产1卡2卡| 丝袜在线中文字幕| 99国产精品一区二区蜜桃av| 精品第一国产精品| 精品国产一区二区三区四区第35| 制服诱惑二区| 激情在线观看视频在线高清| 欧美日韩瑟瑟在线播放| 亚洲,欧美精品.| 亚洲欧美激情综合另类| 亚洲国产精品999在线| 老司机福利观看| 久久精品aⅴ一区二区三区四区| 国产片内射在线| 88av欧美| 亚洲av中文字字幕乱码综合 | 国内精品久久久久精免费| 亚洲精品美女久久久久99蜜臀| 可以在线观看毛片的网站| 亚洲中文av在线| 在线观看舔阴道视频| 性色av乱码一区二区三区2| 激情在线观看视频在线高清| 桃红色精品国产亚洲av| 动漫黄色视频在线观看| 亚洲精华国产精华精| 女性被躁到高潮视频| 一级片免费观看大全| 两性午夜刺激爽爽歪歪视频在线观看 | 又黄又爽又免费观看的视频| 久久久久国内视频| 色哟哟哟哟哟哟| 黄色片一级片一级黄色片| 国产野战对白在线观看| √禁漫天堂资源中文www| 国产欧美日韩一区二区三| 午夜免费鲁丝| 亚洲性夜色夜夜综合| 亚洲va日本ⅴa欧美va伊人久久| 国产精品久久久久久人妻精品电影| 一级黄色大片毛片| 男人操女人黄网站| 免费看a级黄色片| 亚洲av电影不卡..在线观看| 宅男免费午夜| 精品高清国产在线一区| 两个人视频免费观看高清| 侵犯人妻中文字幕一二三四区| 制服诱惑二区| 悠悠久久av| 精品一区二区三区四区五区乱码| 一级a爱片免费观看的视频| 日本在线视频免费播放| 天堂√8在线中文| 91在线观看av| 欧美黑人巨大hd| 日韩中文字幕欧美一区二区| 国产国语露脸激情在线看| 精品久久久久久,| 日本三级黄在线观看| 日本精品一区二区三区蜜桃| 欧美日本视频| 精品久久久久久久末码| 老熟妇仑乱视频hdxx| 亚洲成a人片在线一区二区| 国产不卡一卡二| 久久热在线av| 久久久精品欧美日韩精品| 久久久久国产精品人妻aⅴ院| a在线观看视频网站| 久久久久久人人人人人| 成年女人毛片免费观看观看9| 久久精品aⅴ一区二区三区四区| 男女那种视频在线观看| 亚洲人成77777在线视频| 欧美激情久久久久久爽电影| 欧美成狂野欧美在线观看| av天堂在线播放| 国产精品久久电影中文字幕| 天堂动漫精品| 中文字幕精品免费在线观看视频| 这个男人来自地球电影免费观看| 嫩草影视91久久| 欧美日韩亚洲国产一区二区在线观看| 亚洲成av片中文字幕在线观看| 9191精品国产免费久久| 亚洲五月天丁香| 悠悠久久av| 女人被狂操c到高潮| 少妇被粗大的猛进出69影院| 国产精品电影一区二区三区| 久久精品国产亚洲av高清一级| 亚洲成人免费电影在线观看| 亚洲欧美激情综合另类| 国产蜜桃级精品一区二区三区| 看免费av毛片| 性色av乱码一区二区三区2| 亚洲人成网站在线播放欧美日韩| 一进一出抽搐动态| 国产又色又爽无遮挡免费看| 免费在线观看日本一区| 日本a在线网址| 老司机深夜福利视频在线观看| 69av精品久久久久久| 国产一区二区三区视频了| 亚洲免费av在线视频| 色综合亚洲欧美另类图片| 免费一级毛片在线播放高清视频| 久久久久久久精品吃奶| 一区二区三区高清视频在线| 色综合欧美亚洲国产小说| 亚洲成av片中文字幕在线观看| 美国免费a级毛片| 又大又爽又粗| 国产av在哪里看| 国产成人欧美在线观看| av天堂在线播放| 国产精品二区激情视频| 桃红色精品国产亚洲av| 女同久久另类99精品国产91| 99国产精品一区二区蜜桃av| 国产亚洲精品av在线| a级毛片a级免费在线| 精品国产乱子伦一区二区三区| 变态另类丝袜制服| 久久青草综合色| 变态另类丝袜制服| 丁香六月欧美| 亚洲熟妇中文字幕五十中出| 黄片播放在线免费| 国产精品久久久久久精品电影 | 免费电影在线观看免费观看| 中文资源天堂在线| 午夜免费激情av| 国产三级在线视频| 最近最新中文字幕大全电影3 | 窝窝影院91人妻| 久久久久久久精品吃奶| 精品电影一区二区在线| 亚洲国产精品久久男人天堂| 国产三级黄色录像| 欧美精品亚洲一区二区| 久久久精品国产亚洲av高清涩受| 色综合亚洲欧美另类图片| 午夜日韩欧美国产| 99riav亚洲国产免费| 人人妻人人澡欧美一区二区| 观看免费一级毛片| 日本免费a在线| 亚洲 欧美 日韩 在线 免费| 女人被狂操c到高潮| 亚洲久久久国产精品| 少妇熟女aⅴ在线视频| 一级作爱视频免费观看| 欧美黑人精品巨大| 热re99久久国产66热| 免费看十八禁软件| 亚洲美女黄片视频| 亚洲欧美日韩高清在线视频| 国产成人精品无人区| 老熟妇仑乱视频hdxx| 法律面前人人平等表现在哪些方面| 美女高潮到喷水免费观看| 黄色成人免费大全| 两人在一起打扑克的视频| 亚洲男人天堂网一区| 天天一区二区日本电影三级| 一级作爱视频免费观看| 久热这里只有精品99| 91在线观看av| 成人手机av| 国产精品久久久av美女十八| 亚洲 欧美一区二区三区| 老鸭窝网址在线观看| 韩国精品一区二区三区| 91大片在线观看| 丝袜美腿诱惑在线| 欧美亚洲日本最大视频资源| 亚洲av电影不卡..在线观看| 亚洲无线在线观看| 国产视频一区二区在线看| 女人被狂操c到高潮| 亚洲精品美女久久av网站| 国语自产精品视频在线第100页| 久久这里只有精品19| 黄片小视频在线播放| 欧美精品啪啪一区二区三区| 99久久国产精品久久久| 久久中文看片网| 熟妇人妻久久中文字幕3abv| 99精品久久久久人妻精品| 国产av一区二区精品久久| 美国免费a级毛片| 麻豆成人午夜福利视频| 久久精品夜夜夜夜夜久久蜜豆 | 欧美日韩精品网址| 免费在线观看成人毛片| 99精品欧美一区二区三区四区| 亚洲国产精品合色在线| 人人妻人人澡人人看| 宅男免费午夜| 欧美日韩瑟瑟在线播放| 午夜激情福利司机影院| 亚洲精品色激情综合| 白带黄色成豆腐渣| e午夜精品久久久久久久| 国产av不卡久久| 国产精品免费视频内射| 国产精品亚洲av一区麻豆| 午夜亚洲福利在线播放| 亚洲专区中文字幕在线| 最近最新中文字幕大全电影3 | 在线观看舔阴道视频| 欧美色视频一区免费| 日韩欧美一区二区三区在线观看|