• <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è)計的冗余度
    热99re8久久精品国产| 亚洲av电影在线进入| 狠狠婷婷综合久久久久久88av| 热99久久久久精品小说推荐| 99久久人妻综合| 999久久久精品免费观看国产| 久久久精品区二区三区| 久久青草综合色| 亚洲人成电影观看| 99国产精品免费福利视频| 五月开心婷婷网| 亚洲视频免费观看视频| 日韩成人在线观看一区二区三区| 久久精品国产a三级三级三级| 另类精品久久| 蜜桃国产av成人99| 成人av一区二区三区在线看| 国产高清国产精品国产三级| 成人国产av品久久久| 黄片大片在线免费观看| 宅男免费午夜| 国产亚洲一区二区精品| 美国免费a级毛片| www.精华液| 久久久久久久精品吃奶| 国产日韩欧美亚洲二区| netflix在线观看网站| 2018国产大陆天天弄谢| 人人妻人人添人人爽欧美一区卜| 午夜激情久久久久久久| 法律面前人人平等表现在哪些方面| 久久 成人 亚洲| 一级片'在线观看视频| 99国产综合亚洲精品| 国产真人三级小视频在线观看| 国产日韩欧美视频二区| 亚洲va日本ⅴa欧美va伊人久久| 亚洲精品久久午夜乱码| 久久精品亚洲av国产电影网| 99re6热这里在线精品视频| 精品国产亚洲在线| 高清av免费在线| 亚洲综合色网址| 欧美激情极品国产一区二区三区| 乱人伦中国视频| 欧美性长视频在线观看| 免费一级毛片在线播放高清视频 | 亚洲欧美日韩高清在线视频 | 国产真人三级小视频在线观看| 日韩成人在线观看一区二区三区| 午夜福利一区二区在线看| 精品福利观看| 女人精品久久久久毛片| 精品免费久久久久久久清纯 | 免费观看av网站的网址| 性色av乱码一区二区三区2| 国产精品.久久久| 人妻一区二区av| 精品乱码久久久久久99久播| 久久av网站| 黑丝袜美女国产一区| 叶爱在线成人免费视频播放| 亚洲精华国产精华精| 一级片免费观看大全| 黄色视频,在线免费观看| 国产精品久久久久久精品电影小说| 亚洲精品成人av观看孕妇| 在线av久久热| 国产免费现黄频在线看| 天堂中文最新版在线下载| 国产在线免费精品| 久久久久久久精品吃奶| 欧美日韩成人在线一区二区| bbb黄色大片| 人人澡人人妻人| 我要看黄色一级片免费的| 十八禁高潮呻吟视频| 成年人免费黄色播放视频| 少妇裸体淫交视频免费看高清 | 人人妻,人人澡人人爽秒播| 在线 av 中文字幕| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲专区国产一区二区| 国产精品久久久人人做人人爽| 视频区欧美日本亚洲| 99久久国产精品久久久| 我要看黄色一级片免费的| 国产在线观看jvid| 18在线观看网站| 国产精品一区二区在线不卡| 成人av一区二区三区在线看| 久久中文字幕人妻熟女| 日韩中文字幕视频在线看片| 麻豆成人av在线观看| 国产91精品成人一区二区三区 | 久久精品国产a三级三级三级| 黑丝袜美女国产一区| 一本色道久久久久久精品综合| 日韩免费高清中文字幕av| 一个人免费看片子| 1024视频免费在线观看| 国产99久久九九免费精品| 飞空精品影院首页| 免费女性裸体啪啪无遮挡网站| 国产成人一区二区三区免费视频网站| 男人操女人黄网站| 伊人久久大香线蕉亚洲五| 亚洲美女黄片视频| 丁香六月天网| 精品人妻1区二区| 巨乳人妻的诱惑在线观看| 精品一品国产午夜福利视频| 国产伦理片在线播放av一区| 精品乱码久久久久久99久播| 免费久久久久久久精品成人欧美视频| 男人操女人黄网站| 精品一区二区三区四区五区乱码| 99热网站在线观看| 最新在线观看一区二区三区| 大型av网站在线播放| av网站免费在线观看视频| 国产一区二区在线观看av| 国产高清videossex| 欧美久久黑人一区二区| 国产一区二区三区在线臀色熟女 | 最近最新免费中文字幕在线| 无人区码免费观看不卡 | 精品少妇黑人巨大在线播放| 一级片免费观看大全| 亚洲成国产人片在线观看| av电影中文网址| 亚洲美女黄片视频| 黄色视频,在线免费观看| 中文字幕人妻熟女乱码| 久久热在线av| 亚洲精品一卡2卡三卡4卡5卡| 天天影视国产精品| 国产一区二区激情短视频| 久久亚洲真实| 久久久精品国产亚洲av高清涩受| 久久国产精品人妻蜜桃| 岛国在线观看网站| 建设人人有责人人尽责人人享有的| 人人妻,人人澡人人爽秒播| 丁香欧美五月| 国产精品国产av在线观看| 国产亚洲欧美在线一区二区| 视频区图区小说| 无限看片的www在线观看| 欧美黑人精品巨大| 国产极品粉嫩免费观看在线| 岛国在线观看网站| 岛国毛片在线播放| 欧美中文综合在线视频| 三上悠亚av全集在线观看| 国产精品av久久久久免费| 大片电影免费在线观看免费| 一本综合久久免费| 亚洲自偷自拍图片 自拍| 亚洲国产欧美一区二区综合| 男女下面插进去视频免费观看| 亚洲七黄色美女视频| 久久亚洲真实| 久久午夜亚洲精品久久| 国产麻豆69| 多毛熟女@视频| 久久99热这里只频精品6学生| 亚洲精品一卡2卡三卡4卡5卡| 精品人妻1区二区| 亚洲精品粉嫩美女一区| 在线av久久热| 老汉色∧v一级毛片| 久久精品人人爽人人爽视色| 男女高潮啪啪啪动态图| 午夜免费成人在线视频| 国产精品偷伦视频观看了| 热re99久久国产66热| 巨乳人妻的诱惑在线观看| 桃红色精品国产亚洲av| netflix在线观看网站| 久久亚洲精品不卡| 别揉我奶头~嗯~啊~动态视频| 丰满少妇做爰视频| 亚洲男人天堂网一区| 亚洲伊人久久精品综合| 一级a爱视频在线免费观看| 一二三四社区在线视频社区8| 老司机午夜十八禁免费视频| 国产精品免费视频内射| 国产男女超爽视频在线观看| svipshipincom国产片| 12—13女人毛片做爰片一| 免费观看人在逋| 制服人妻中文乱码| 国产精品久久久av美女十八| 中文字幕人妻熟女乱码| 侵犯人妻中文字幕一二三四区| 黄色视频不卡| 正在播放国产对白刺激| 久久久国产欧美日韩av| 男人操女人黄网站| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品粉嫩美女一区| 亚洲欧洲精品一区二区精品久久久| 18禁观看日本| 黄片大片在线免费观看| 美女国产高潮福利片在线看| 人人妻人人澡人人爽人人夜夜| 欧美成人免费av一区二区三区 | 黄色视频在线播放观看不卡| 亚洲人成77777在线视频| 男女高潮啪啪啪动态图| 国产在线观看jvid| 青草久久国产| 精品第一国产精品| 亚洲国产精品一区二区三区在线| 亚洲欧美一区二区三区黑人| a在线观看视频网站| 国产精品秋霞免费鲁丝片| av一本久久久久| 国产成人精品久久二区二区免费| 天天躁夜夜躁狠狠躁躁| 国产午夜精品久久久久久| 美女主播在线视频| 亚洲免费av在线视频| 叶爱在线成人免费视频播放| 桃红色精品国产亚洲av| 操出白浆在线播放| 国产成人系列免费观看| 亚洲欧美色中文字幕在线| 久久久国产精品麻豆| 制服诱惑二区| 亚洲综合色网址| 一本综合久久免费| 三上悠亚av全集在线观看| 性高湖久久久久久久久免费观看| 电影成人av| 色综合欧美亚洲国产小说| 欧美精品啪啪一区二区三区| 亚洲伊人色综图| 91成年电影在线观看| 三级毛片av免费| 亚洲专区国产一区二区| 亚洲精品成人av观看孕妇| 性高湖久久久久久久久免费观看| 国产成人av激情在线播放| 日韩 欧美 亚洲 中文字幕| 中文字幕人妻熟女乱码| 国产日韩欧美视频二区| 久久99一区二区三区| 亚洲 欧美一区二区三区| 亚洲精品久久成人aⅴ小说| 欧美人与性动交α欧美精品济南到| 不卡av一区二区三区| 啦啦啦 在线观看视频| www.精华液| 亚洲精品在线观看二区| 成人黄色视频免费在线看| 国产精品久久久久久精品电影小说| 久久久国产精品麻豆| 高清在线国产一区| 日韩 欧美 亚洲 中文字幕| bbb黄色大片| 亚洲五月色婷婷综合| 欧美日韩亚洲综合一区二区三区_| 久久人人爽av亚洲精品天堂| 另类亚洲欧美激情| 色综合欧美亚洲国产小说| 嫩草影视91久久| 国产精品久久久久久人妻精品电影 | 亚洲五月色婷婷综合| 黑丝袜美女国产一区| 最近最新中文字幕大全电影3 | 欧美成人午夜精品| 国产真人三级小视频在线观看| 亚洲av美国av| 国产成人精品无人区| 国产国语露脸激情在线看| 日日爽夜夜爽网站| 成人国产一区最新在线观看| 国产在线视频一区二区| 在线观看免费午夜福利视频| 精品国产乱子伦一区二区三区| 精品一区二区三区视频在线观看免费 | 老汉色av国产亚洲站长工具| 老司机午夜十八禁免费视频| 97人妻天天添夜夜摸| 成人亚洲精品一区在线观看| av网站在线播放免费| 亚洲精品在线美女| 免费观看av网站的网址| 免费久久久久久久精品成人欧美视频| 精品少妇一区二区三区视频日本电影| 岛国毛片在线播放| 夜夜夜夜夜久久久久| 成人永久免费在线观看视频 | 国产成人精品久久二区二区91| 大片电影免费在线观看免费| 激情在线观看视频在线高清 | 91精品三级在线观看| 国产精品秋霞免费鲁丝片| 麻豆成人av在线观看| 天堂中文最新版在线下载| 国产成人精品久久二区二区免费| 在线天堂中文资源库| 搡老乐熟女国产| 99re在线观看精品视频| 一本—道久久a久久精品蜜桃钙片| 午夜视频精品福利| 久久午夜亚洲精品久久| 另类精品久久| 一区二区三区精品91| 国产精品久久久久久精品电影小说| 丝袜美腿诱惑在线| 欧美乱码精品一区二区三区| 国产男女超爽视频在线观看| 精品少妇内射三级| 国产免费现黄频在线看| 精品国产乱子伦一区二区三区| 1024香蕉在线观看| 国产成人啪精品午夜网站| 国产精品久久久av美女十八| 欧美黄色淫秽网站| 亚洲国产毛片av蜜桃av| 国产成人免费观看mmmm| 91字幕亚洲| 欧美精品一区二区免费开放| 丰满饥渴人妻一区二区三| 亚洲国产欧美一区二区综合| 欧美日韩av久久| 久久人妻福利社区极品人妻图片| 一本综合久久免费| 肉色欧美久久久久久久蜜桃| 久久精品国产综合久久久| 亚洲熟妇熟女久久| 亚洲精品久久午夜乱码| 国产欧美亚洲国产| 啦啦啦中文免费视频观看日本| 国产精品免费一区二区三区在线 | 午夜免费鲁丝| 亚洲成av片中文字幕在线观看| www.精华液| 国产又色又爽无遮挡免费看| 99久久精品国产亚洲精品| 日韩免费av在线播放| tocl精华| 国产精品久久久久久精品古装| 又大又爽又粗| 亚洲欧美激情在线| 成年人黄色毛片网站| 少妇 在线观看| cao死你这个sao货| 美女国产高潮福利片在线看| 亚洲七黄色美女视频| 成人永久免费在线观看视频 | av免费在线观看网站| 19禁男女啪啪无遮挡网站| 精品欧美一区二区三区在线| 亚洲综合色网址| 亚洲精品国产色婷婷电影| 亚洲午夜精品一区,二区,三区| 一个人免费看片子| 日韩欧美免费精品| 免费不卡黄色视频| 91九色精品人成在线观看| tube8黄色片| 久久天躁狠狠躁夜夜2o2o| 制服人妻中文乱码| 欧美人与性动交α欧美软件| 亚洲一区中文字幕在线| 精品人妻在线不人妻| 女警被强在线播放| 色婷婷久久久亚洲欧美| 欧美黄色片欧美黄色片| 久久九九热精品免费| 99九九在线精品视频| 欧美午夜高清在线| 桃红色精品国产亚洲av| 久久国产亚洲av麻豆专区| 国产精品香港三级国产av潘金莲| 久久精品亚洲av国产电影网| 亚洲 欧美一区二区三区| 精品一区二区三卡| 亚洲国产欧美在线一区| 国产精品麻豆人妻色哟哟久久| 国产麻豆69| 午夜免费成人在线视频| 成人国语在线视频| 国产不卡一卡二| 亚洲国产中文字幕在线视频| 午夜福利乱码中文字幕| 成人18禁高潮啪啪吃奶动态图| 国产欧美日韩一区二区三| 正在播放国产对白刺激| 亚洲欧美色中文字幕在线| 欧美精品高潮呻吟av久久| 超碰97精品在线观看| 99精品欧美一区二区三区四区| 欧美日韩成人在线一区二区| 热re99久久国产66热| 午夜福利欧美成人| 免费人妻精品一区二区三区视频| 天堂8中文在线网| 欧美+亚洲+日韩+国产| 国产黄频视频在线观看| 亚洲精品中文字幕在线视频| 午夜精品久久久久久毛片777| 老司机靠b影院| 99热国产这里只有精品6| 成人av一区二区三区在线看| 亚洲伊人久久精品综合| 午夜福利视频在线观看免费| 国产男女内射视频| 宅男免费午夜| 看免费av毛片| 亚洲性夜色夜夜综合| 久久精品国产综合久久久| 中文字幕另类日韩欧美亚洲嫩草| av不卡在线播放| 欧美精品一区二区大全| 麻豆乱淫一区二区| 如日韩欧美国产精品一区二区三区| 午夜两性在线视频| 波多野结衣一区麻豆| 久久人人97超碰香蕉20202| 国产一区二区在线观看av| 久久久国产欧美日韩av| 一级片'在线观看视频| 欧美变态另类bdsm刘玥| xxxhd国产人妻xxx| 精品久久久久久久毛片微露脸| 国产成+人综合+亚洲专区| 一夜夜www| 下体分泌物呈黄色| 99riav亚洲国产免费| 999精品在线视频| 99久久国产精品久久久| 波多野结衣av一区二区av| 亚洲avbb在线观看| 欧美黑人欧美精品刺激| 亚洲午夜理论影院| 精品福利永久在线观看| 妹子高潮喷水视频| 老司机午夜十八禁免费视频| 大型黄色视频在线免费观看| 午夜老司机福利片| 少妇精品久久久久久久| 欧美成人免费av一区二区三区 | 国产在线观看jvid| 国产精品欧美亚洲77777| 国产91精品成人一区二区三区 | 精品国产乱码久久久久久男人| 我要看黄色一级片免费的| 欧美在线一区亚洲| 不卡一级毛片| 国产一区二区在线观看av| 天天影视国产精品| 正在播放国产对白刺激| 国产精品九九99| 久久精品亚洲精品国产色婷小说| 久久亚洲真实| 丝袜在线中文字幕| 久久久精品国产亚洲av高清涩受| 精品国产一区二区三区四区第35| 久久影院123| 黄色视频不卡| 亚洲专区字幕在线| 午夜福利在线观看吧| 欧美日韩黄片免| 日韩一卡2卡3卡4卡2021年| 欧美精品高潮呻吟av久久| 久久久精品94久久精品| 亚洲熟妇熟女久久| av不卡在线播放| 一级,二级,三级黄色视频| 黄色 视频免费看| 亚洲国产av新网站| 美女午夜性视频免费| 最近最新中文字幕大全电影3 | 天堂中文最新版在线下载| 欧美精品亚洲一区二区| 亚洲天堂av无毛| 精品亚洲成国产av| 国产精品 欧美亚洲| 99re6热这里在线精品视频| 精品少妇内射三级| 亚洲av欧美aⅴ国产| 中文亚洲av片在线观看爽 | 99国产精品99久久久久| 国产免费福利视频在线观看| 黄网站色视频无遮挡免费观看| 国产国语露脸激情在线看| 啦啦啦免费观看视频1| 国产高清国产精品国产三级| 日韩熟女老妇一区二区性免费视频| 国产在视频线精品| 久久精品成人免费网站| 超碰97精品在线观看| e午夜精品久久久久久久| 欧美黑人精品巨大| 99精品在免费线老司机午夜| 一夜夜www| 国产亚洲av高清不卡| 三级毛片av免费| 亚洲第一av免费看| 国产精品一区二区免费欧美| 亚洲性夜色夜夜综合| 亚洲欧美精品综合一区二区三区| 69精品国产乱码久久久| 国产真人三级小视频在线观看| 人人澡人人妻人| 国产精品自产拍在线观看55亚洲 | 另类亚洲欧美激情| 国产免费视频播放在线视频| 国产伦人伦偷精品视频| 精品国产亚洲在线| 色婷婷av一区二区三区视频| 久久天堂一区二区三区四区| 大片免费播放器 马上看| 国产福利在线免费观看视频| 成人av一区二区三区在线看| 亚洲第一欧美日韩一区二区三区 | 国产男女内射视频| 久久狼人影院| 精品国内亚洲2022精品成人 | 大香蕉久久成人网| 狠狠精品人妻久久久久久综合| 高清av免费在线| 天堂动漫精品| 老司机亚洲免费影院| 后天国语完整版免费观看| 日韩成人在线观看一区二区三区| 日韩中文字幕视频在线看片| 午夜免费成人在线视频| 女人精品久久久久毛片| 99久久人妻综合| avwww免费| 欧美激情高清一区二区三区| 国产成+人综合+亚洲专区| 亚洲av电影在线进入| 9191精品国产免费久久| 国产免费现黄频在线看| 色播在线永久视频| 美女午夜性视频免费| 精品国产乱码久久久久久男人| 黑人操中国人逼视频| 日韩欧美三级三区| 欧美大码av| 久久毛片免费看一区二区三区| 久久人人爽av亚洲精品天堂| 亚洲,欧美精品.| 一级毛片女人18水好多| 亚洲精品美女久久av网站| 久久99热这里只频精品6学生| 日本黄色日本黄色录像| 国产激情久久老熟女| 亚洲成人国产一区在线观看| 国产一区二区 视频在线| 日本黄色日本黄色录像| 丁香六月欧美| 欧美激情高清一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 一区二区av电影网| 亚洲精品久久成人aⅴ小说| 午夜福利视频精品| 老鸭窝网址在线观看| 久久精品国产a三级三级三级| 丰满人妻熟妇乱又伦精品不卡| 黑丝袜美女国产一区| 精品少妇一区二区三区视频日本电影| 桃花免费在线播放| 十八禁网站免费在线| 欧美性长视频在线观看| 精品国产乱码久久久久久小说| 午夜福利一区二区在线看| 天天躁日日躁夜夜躁夜夜| 一区二区三区精品91| 免费少妇av软件| 啦啦啦中文免费视频观看日本| 久久精品国产亚洲av高清一级| 国产欧美日韩一区二区三区在线| 男女高潮啪啪啪动态图| 免费高清在线观看日韩| 另类精品久久| 美女扒开内裤让男人捅视频| 亚洲av日韩精品久久久久久密| 国产视频一区二区在线看| 可以免费在线观看a视频的电影网站| 91精品三级在线观看| 精品福利观看| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品国产色婷婷电影| 成在线人永久免费视频| 老汉色∧v一级毛片| 国产成人啪精品午夜网站| 国产精品自产拍在线观看55亚洲 | 国产男靠女视频免费网站| 两个人看的免费小视频| 大型av网站在线播放| 好男人电影高清在线观看| 丰满迷人的少妇在线观看| 亚洲av第一区精品v没综合| 亚洲全国av大片| 亚洲欧洲精品一区二区精品久久久| 亚洲九九香蕉| 欧美日韩亚洲高清精品| 国产人伦9x9x在线观看| 免费观看人在逋| 午夜精品国产一区二区电影| 搡老乐熟女国产| 天堂中文最新版在线下载| 好男人电影高清在线观看| 国产精品一区二区免费欧美| 国内毛片毛片毛片毛片毛片| 女人被躁到高潮嗷嗷叫费观|