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

    基于獨立分量分析基的地震隨機(jī)噪聲壓制

    2016-04-26 08:33:22孫成禹
    石油物探 2016年2期

    孫成禹,邵 婕,藍(lán) 陽,唐 杰

    (中國石油大學(xué)(華東)地球科學(xué)與技術(shù)學(xué)院,山東青島266580)

    ?

    基于獨立分量分析基的地震隨機(jī)噪聲壓制

    孫成禹,邵婕,藍(lán)陽,唐杰

    (中國石油大學(xué)(華東)地球科學(xué)與技術(shù)學(xué)院,山東青島266580)

    摘要:傳統(tǒng)的獨立分量分析(Independent Component Analysis,ICA)去噪方法假設(shè)地震記錄的相鄰道含有相同的隨機(jī)噪聲,僅適用于同相軸較平的地震記錄,去噪效果并不顯著。為了改善ICA方法對高斯隨機(jī)噪聲的壓制效果,首先通過構(gòu)造度量數(shù)據(jù)非高斯性的目標(biāo)函數(shù)求取地震數(shù)據(jù)的ICA基,將數(shù)據(jù)轉(zhuǎn)換至ICA域;然后通過貝葉斯方法構(gòu)造出滿足非高斯分布的閾值函數(shù),進(jìn)行閾值法去噪處理。為了滿足獨立分量分析的假設(shè)條件,將地震數(shù)據(jù)進(jìn)行分塊處理,并假設(shè)每個數(shù)據(jù)塊與整體的數(shù)據(jù)含有相似的數(shù)據(jù)結(jié)構(gòu)。理論模型及實際資料試算結(jié)果表明,該方法可以有效地壓制剖面中的高斯隨機(jī)噪聲,對含復(fù)雜界面的數(shù)據(jù)也十分有效,具有較好的應(yīng)用價值。

    關(guān)鍵詞:獨立分量分析;ICA基;高斯隨機(jī)噪聲壓制;貝葉斯閾值函數(shù)

    壓制隨機(jī)噪聲的方法有很多,如均值濾波、中值濾波、基于小波變換或曲波變換的閾值去噪方法以及基于獨立分量分析(ICA)的去噪方法等。但隨著地震勘探程度的不斷提高及勘探目標(biāo)的日益復(fù)雜,對去噪技術(shù)的要求也越來越高。獨立分量分析技術(shù)是實現(xiàn)盲源分離的有效手段,它考慮了數(shù)據(jù)的高階統(tǒng)計特性,因而能更加全面地表達(dá)數(shù)據(jù)的本質(zhì)特征[1]。劉喜武等[2]假設(shè)地震記錄中的有效信號和隨機(jī)噪聲具有統(tǒng)計意義的相互獨立性質(zhì)且服從非高斯分布,利用某道記錄的多次觀測結(jié)果或鄰近多道記錄,通過獨立分量分析算法直接從地震記錄中分離出有效信號及非高斯隨機(jī)噪聲[3]。此后很多學(xué)者對獨立分量分析的算法進(jìn)行了改進(jìn),如李大衛(wèi)等[4]提出了模擬退火獨立分量分析方法,并將其應(yīng)用于共反射點道集中有效信號和干擾信號的分離,提高了分離信號的質(zhì)量。張銀雪等[5]為了改善常規(guī)固定步長獨立分量分析算法在疊前資料中的去噪效果,提出了基于混沌粒子群優(yōu)化(Particle Swarm Optimization,PSO)的改進(jìn)ICA算法。這些改進(jìn)的獨立分量分析方法能有效地壓制地震數(shù)據(jù)中的高斯隨機(jī)噪聲,但通常假設(shè)輸入的相鄰道記錄中包含相同的隨機(jī)噪聲,這在實際生產(chǎn)中很難得到滿足。為了改善ICA技術(shù)壓制加性高斯隨機(jī)噪聲的效果,呂文彪等[6]先應(yīng)用兩步特征值分解法去除部分加性噪聲的影響,再利用ICA算法成功地分離出原信號。郭科等[7]利用改進(jìn)的ICA預(yù)處理算法去除加性高斯噪聲,采用特征矩陣聯(lián)合對角化(Joint Approximate Diagonalization of Eigenmatrices)盲分離算法分離出有效信號和非高斯分布的隨機(jī)噪聲。王維強(qiáng)等[8]利用經(jīng)驗?zāi)B(tài)分解(Empirical Mode Decomposition,EMD)能將信號自適應(yīng)分解為不同尺度振動模態(tài)的優(yōu)點及獨立分量分析能提取獨立源信號的優(yōu)勢,提出了一種EMD和ICA相結(jié)合的去噪算法,很好地實現(xiàn)了地震有效信號和隨機(jī)噪聲的分離。

    上述方法雖然改善了ICA算法的去噪效果,但通常只適用于同相軸較平的疊后記錄,且運(yùn)算效率比較低。本文通過構(gòu)造度量數(shù)據(jù)非高斯性的目標(biāo)函數(shù)及尋優(yōu)算法確定輸入數(shù)據(jù)的ICA基,將地震記錄轉(zhuǎn)換至ICA域進(jìn)行去噪處理,以進(jìn)一步改善ICA方法對高斯隨機(jī)噪聲的壓制效果。

    1獨立分量分析(ICA)

    1.1基函數(shù)

    在信號處理領(lǐng)域,為了有效地分析數(shù)據(jù)的本征結(jié)構(gòu),人們通常將其表述為一系列基向量的線性組合,提出了各種不同的基函數(shù),如傅里葉變換中的正弦、余弦基,以及小波基、曲波基等。這些基函數(shù)通常有確定的表達(dá)式,可以將數(shù)據(jù)變換到不同的域中進(jìn)行分析。由于不同數(shù)據(jù)所隱含的數(shù)據(jù)結(jié)構(gòu)并不相同,確定的基函數(shù)無法突出信號的某些特征,因此,我們定義包含基函數(shù)所需性質(zhì)的目標(biāo)函數(shù),對輸入數(shù)據(jù)進(jìn)行訓(xùn)練,通過優(yōu)化算法確定適用于任意輸入數(shù)據(jù)的基函數(shù)。

    設(shè)數(shù)據(jù)x(k,l)的大小為M×N,k=1,2,…,M,l=1,2,…,N,可表示為K個基向量pj(k,l)的線性組合:

    (1)

    式中:yj為對應(yīng)基向量pj(k,l)的表示系數(shù),yT=[y1,…,yK],P=[p1(k,l),…,pK(k,l)]。通過(1)式,觀測信號被投影到一組基向量張成的空間內(nèi)。由此看來,確定一組能夠捕捉輸入數(shù)據(jù)大部分結(jié)構(gòu)的基函數(shù)是有效描述信號特征問題的關(guān)鍵,可以通過獨立分量分析的方法來獲得基向量pj(k,l)及其對應(yīng)的表示系數(shù)yj。

    1.2ICA基

    獨立分量分析本質(zhì)上也是一種優(yōu)化算法,它可以確定統(tǒng)計意義上相互獨立的基向量[9]。中心極限定理表明,如果觀測信號可以表示為一組基向量的組合,且具有最小的高斯性,信號便是獨立的。因此,統(tǒng)計獨立和非高斯性是等價的,構(gòu)造ICA目標(biāo)函數(shù)的關(guān)鍵便是度量分離結(jié)果的非高斯性[6]。當(dāng)(1)式中表示系數(shù)yj的非高斯性達(dá)到最大時,所確定的基向量相互獨立。

    度量非高斯性的方法有很多,如高階累積量、熵及互信息等[10]。本文選用負(fù)熵的方法,其定義如下:

    (2)

    其中,yG是均值為0、與y具有相同方差的高斯隨機(jī)變量,H(y)是隨機(jī)變量的熵,其定義式為:

    (3)

    式中,p(y)為y的概率密度函數(shù)。由(2)式可知,當(dāng)y服從高斯分布時,J(y)=0;非高斯變量的熵H(y)為非負(fù)值,非高斯性越強(qiáng),其值越大。

    (4)

    其中,E{·}為隨機(jī)變量的均值;G(·)是非二次函數(shù),其表達(dá)形式有很多[12],本文選用的公式為:

    (5)

    式中:a,b為常數(shù);c為取值較小、控制數(shù)值穩(wěn)定性的常數(shù)(c≈0.1)。

    為了求取J(y)的最大值,HYVRINEN[11]提出了固定點算法,又稱為FastICA方法。該算法首先對原始數(shù)據(jù)x進(jìn)行預(yù)處理,主要是進(jìn)行零均值化及白化,其中白化可通過主成分分析(Principal Component Analysis,PCA)方法實現(xiàn)。用主成分分析法對零均值化處理后的數(shù)據(jù)的協(xié)方差矩陣進(jìn)行特征值分解,可實現(xiàn)數(shù)據(jù)的降維處理。這一過程用線性變換V來表示,則有:

    (6)

    (7)

    通常將V稱作為PCA基(由主成分分析算法求得)[13]。從公式(7)中選取L

    (8)

    通過牛頓迭代法求解方程(8),最終得到FastICA算法的迭代公式為:

    (9)

    (10)

    將求得的ICA基及其逆矩陣作用于相應(yīng)的數(shù)據(jù),便可以實現(xiàn)空間域與ICA域間的相互映射。

    圖1 地震記錄分塊及重排列

    由ICA的約束條件可知,觀測信號的個數(shù)應(yīng)不少于源信號的個數(shù),但對于地震數(shù)據(jù)而言,通常只有一個觀測結(jié)果。傳統(tǒng)的ICA方法通常選取某道記錄鄰近的幾道記錄作為輸入,并假設(shè)其含有相同的隨機(jī)噪聲。但這種假設(shè)并不符合實際情況,去噪效果也不顯著。本文借鑒圖像處理中的分塊處理方式[15-16],將地震記錄進(jìn)行分塊處理,得到相同大小的子剖面。假設(shè)這些子剖面包含的數(shù)據(jù)結(jié)構(gòu)與原始剖面相似,將每一子剖面作為觀測信號進(jìn)行ICA處理,便可以得到反映地震記錄整體性質(zhì)的基函數(shù)。在實際處理中,為了簡化運(yùn)算過程,通常將分塊得到的二維子剖面按一定的順序(圖1中實線箭頭所示)轉(zhuǎn)化為一維數(shù)組。此外,為了節(jié)省運(yùn)算時間,并不對所有的一維子塊進(jìn)行處理,而是隨機(jī)地選取k0塊(k0≥1000)構(gòu)成一個新的輸入矩陣,如圖1所示。將此輸入矩陣作為上述ICA算法的輸入,求出分離矩陣W及混合矩陣W-1,進(jìn)而求出ICA基。圖2為模擬的正演記錄以及通過上述方法得到的ICA基。在計算過程中,地震數(shù)據(jù)被劃分為16×16(樣點數(shù))大小的子塊。為了提高運(yùn)算效率,從劃分的數(shù)據(jù)塊中隨機(jī)地選取1500塊作為ICA算法的輸入。為了形象地顯示基函數(shù),將求得的基函數(shù)(長度n=256)變換為16×16的矩陣。

    圖2 正演記錄(a)和ICA基(b)

    2基于ICA基的閾值法去噪

    野外觀測的含噪聲地震記錄如下:z=s+n

    在實際應(yīng)用中,有效信號合理的概率密度估計

    p(s)及其參數(shù)選取問題是去噪技術(shù)的關(guān)鍵,但對于各接收道而言,數(shù)據(jù)之間存在著復(fù)雜的統(tǒng)計依賴性,難以確定統(tǒng)一的p(s)函數(shù)。為此,借助獨立分量分析的特有性質(zhì),使原始數(shù)據(jù)變換后各分量間的統(tǒng)計依賴性更小(即相互獨立),且使分離后各信號的非高斯性達(dá)到最大。即:Pz=Ps+Pn

    (13)其中,P為上述ICA算法確定的變換算子,它將數(shù)據(jù)由空間域映射到ICA域。對于任意的P,Pn仍服從高斯分布,因此,可以選取具有非高斯分布特點的概率密度函數(shù)作為先驗信息,并通過(12)式中的貝葉斯方法實現(xiàn)對隨機(jī)噪聲的壓制。

    服從非高斯分布的概率密度函數(shù)有多種不同的形式,本文選用的形式如

    (14)其中,d為向量s的標(biāo)準(zhǔn)差,α為控制p(sj)形式的參數(shù),其表達(dá)式如下:

    (15)其中,k=d2p(0)2。當(dāng)α→∞時,公式(14)為拉普拉斯分布。

    將公式(14)代入公式(12)中,并求其最大值,推導(dǎo)后得:

    (16)

    圖3 基于ICA基的去噪處理算法流程

    除了上述閾值函數(shù)外,我們也可以考慮其它形式的概率分布,從而得到不同的閾值算法[19]。圖3為基于ICA基的去噪處理算法流程,其主要包括如下步驟。1) 將大小為M×N的地震數(shù)據(jù)進(jìn)行分塊,即選取大小為I×I的窗函數(shù)對原始記錄從左到右或從上到下做重疊分塊,并按照其存儲位置將其轉(zhuǎn)化成I×I×1的向量(按列順序)。所選用的窗函數(shù)如下:

    (17)

    2) 從分塊處理得到的數(shù)據(jù)塊中隨機(jī)選取k0塊(k0≥1000)作為FastICA算法的輸入,求出ICA域的分離矩陣W及混合矩陣W-1,進(jìn)而求出將數(shù)據(jù)由空間域映射至ICA域的基向量P-1。

    3) 利用求得的ICA基將地震數(shù)據(jù)映射到ICA域,并采用蒙特卡羅方法[20]求取噪聲的方差估算值σ,然后通過閾值方法壓制隨機(jī)噪聲。

    4) 通過P將處理后的數(shù)據(jù)由ICA域映射到空間域,并按照先前的分塊順序?qū)?shù)據(jù)進(jìn)行重新組合。

    5) 由于ICA算法固有的不確定性,處理后數(shù)據(jù)的振幅值相對于原始數(shù)據(jù)整體相差一個倍數(shù),因此需要對處理后的數(shù)據(jù)進(jìn)行重構(gòu),但這并不改變數(shù)據(jù)間振幅的相對關(guān)系。

    3模型測試及實際資料處理

    3.1理論模型測試

    為了驗證上述算法的去噪效果,建立如圖4a所示的速度模型,并通過正演模擬得到圖4b所示的地震剖面。該模型既包含了常見的水平及傾斜界面,又包含了復(fù)雜的彎曲界面、斷層及透鏡體。在原始剖面中加入信噪比為1.08的高斯隨機(jī)噪聲(圖5a),圖5b為采用傳統(tǒng)的獨立分量分析方法去除隨機(jī)噪聲后的剖面;圖5c為采用本文方法去噪后的剖面,去噪前數(shù)據(jù)被劃分為10×10大小的子塊,并隨機(jī)地選取1000塊數(shù)據(jù)作為ICA算法的輸入。

    對比圖5b和圖5c可知,本文方法的去噪效果要明顯好于傳統(tǒng)方法,剖面中的大部分隨機(jī)噪聲都得到了有效壓制,且有效信號的損失較少,即使是那些彎曲程度較大的復(fù)雜同相軸,有效信號的損失也比較少,相比之下,傳統(tǒng)方法在復(fù)雜同相軸處的有效信號損失較為嚴(yán)重。圖5d為含噪聲的地震剖面與本文方法去噪后的地震剖面差值,可以看出差值剖面中殘留的有效信號較少。

    圖4 理論速度模型(a)和正演地震剖面(b)

    圖5 理論模型的去噪效果a 含噪聲的原始剖面; b 傳統(tǒng)方法去噪后的剖面; c 本文方法去噪后的剖面; d 本文方法去除的隨機(jī)噪聲

    對圖5c中虛線處(CDP2400)的地震道進(jìn)行波形分析,結(jié)果如圖6所示。圖6a為上述地震道正演模擬的波形及加噪后的波形對比圖,可以看出,深層的弱信號完全淹沒在噪聲中,無法得到有效識別。圖6b為該地震道正演模擬的波形及去噪后的波形對比圖,可以看出深層的弱信號得到了有效恢復(fù),虛線圈內(nèi)所示透鏡體處的信號幾乎沒有損失,但斷層處的信號略有損失,這主要是在去噪過程中直接將小于閾值的幅值賦零造成的。

    圖6 圖5c中虛線處(CDP2400)地震道去噪效果分析a 有效信號及含噪聲數(shù)據(jù)的波形分析; b 有效信號及去噪后數(shù)據(jù)的波形分析

    3.2實際資料處理

    選取A地區(qū)實際地震資料,利用本文獨立分量分析方法對其進(jìn)行去噪處理。數(shù)據(jù)被劃分為10×10大小的子塊,并隨機(jī)地選取1500塊數(shù)據(jù)作為ICA算法的輸入。去噪結(jié)果如圖7所示,其中,圖7a為含噪聲的原始地震剖面,圖7b為采用傳統(tǒng)獨立分量分析方法去噪后的剖面,圖7c為采用本文方法去噪后的剖面。對比圖7a至圖7c可以看出,兩種方法都能使原始剖面中的隨機(jī)噪聲得到有效壓制,但本文方法的去噪效果要好于傳統(tǒng)方法,整個剖面上的同相軸變得更加連續(xù)、光滑,如圖中紅線框內(nèi)所示。圖7d為本文方法去除的隨機(jī)噪聲,即圖7a和圖7c兩圖的差值剖面,可以看出,本文方法去噪后損失的有效信號比較少。

    為了精確地顯示去噪處理對原始數(shù)據(jù)信噪比的改善情況,計算了兩種方法去噪前后地震剖面的信噪比譜[21],結(jié)果如圖8所示。可以看出,兩種方法去噪后各個頻帶內(nèi)的信噪比都得到了改善,但本文方法去噪后的信噪比要大于傳統(tǒng)方法去噪后的信噪比。計算信噪比譜的公式如下:

    (18)

    圖7 A地區(qū)實際資料去噪效果a 含噪聲的地震剖面; b 原有方法去噪后的剖面; c 本文方法去噪后的剖面; d 本文方法去除的隨機(jī)噪聲

    式中:Ps(f)為有效信號的功率譜,即第i道記錄與第i+1道記錄的互相關(guān);Pz(f)為n道記錄的平均功率譜。

    (19)

    (20)

    其中,Zi(f)為第i(i=1,2,…,n)道原始記錄的振幅譜。

    圖9a為B地區(qū)含噪聲的地震剖面,該剖面同相軸較為傾斜,且含有斷點現(xiàn)象,隨機(jī)噪聲使剖面中的斷點不清晰,無法有效地追蹤整個斷層的延伸情況,如圖9a中紅色箭頭所示。圖9b為采用本文方法去噪后的剖面,同相軸變得光滑、連續(xù),隨機(jī)噪聲得到了有效的壓制,斷點清晰可見,更加有利于小斷層的追蹤。圖9c為采用傳統(tǒng)方法去噪后的剖面,雖然隨機(jī)噪聲得到了一定程度的壓制,但與圖9b 相比,去噪效果較差。圖9d為本文方法去除的隨機(jī)噪聲。

    圖8 A地區(qū)實際資料去噪前后的信噪比譜

    為了進(jìn)一步分析本文方法的去噪效果,截取圖9a 中虛線框內(nèi)的斷層區(qū)域進(jìn)行局部顯示,如圖10所示。由圖10可見,去噪前虛線框內(nèi)藍(lán)色箭頭所示的斷點位置無法進(jìn)行有效識別(圖10a),采用本文方法去噪后不僅可以準(zhǔn)確地確定虛線框內(nèi)斷點的位置,實線框內(nèi)的斷點也更加清晰(圖10b),而原有方法去噪后部分?jǐn)帱c處仍較為模糊,如圖10c中虛線框內(nèi)箭頭所示。求取圖9中兩種方法去噪前后地震剖面的信噪比譜,如圖11所示。由圖11可以看出,對于較為復(fù)雜的地震剖面,原有方法去噪后的結(jié)果信噪比改善不大,而本文方法去噪后各個頻帶內(nèi)的信噪比都得到大大提升。

    圖9 B地區(qū)實際資料去噪效果a 含噪聲的地震剖面; b 本文方法去噪后的剖面; c 原有方法去噪后的剖面; d 本文方法去除的隨機(jī)噪聲

    圖10 B地區(qū)實際資料去噪前后的局部顯示a 含噪聲的地震剖面; b 本文方法去噪后的剖面; c 原有方法去噪后的剖面

    圖11 B地區(qū)實際資料去噪前后的信噪比譜

    4結(jié)束語

    本文通過構(gòu)造度量數(shù)據(jù)非高斯性的目標(biāo)函數(shù),利用尋優(yōu)算法確定適應(yīng)數(shù)據(jù)自身特點的ICA基,將數(shù)據(jù)轉(zhuǎn)換至ICA域進(jìn)行去噪處理,獲得以下結(jié)論與認(rèn)識:

    1) 通過獨立分量分析使分離后有效信號的非高斯性最大,隨機(jī)噪聲仍滿足高斯分布,根據(jù)貝葉斯理論推導(dǎo)出滿足該分布情況下的閾值函數(shù),結(jié)合閾值法去噪技術(shù),可以壓制剖面中的高斯隨機(jī)噪聲。

    2) 與傳統(tǒng)的基于獨立分量分析的去噪算法不同在于,本文方法對數(shù)據(jù)進(jìn)行分塊處理,假設(shè)每一個數(shù)據(jù)塊所含的基函數(shù)與整體數(shù)據(jù)的基函數(shù)相同,僅選取部分?jǐn)?shù)據(jù)塊進(jìn)行訓(xùn)練,運(yùn)算效率得以提高。

    3) 通過構(gòu)造的基函數(shù)將數(shù)據(jù)轉(zhuǎn)換至ICA域進(jìn)行去噪處理,可以有效地避免傳統(tǒng)獨立分量分析方法僅適用于較平同相軸的問題,使其適用于復(fù)雜剖面。

    需要注意的是,本文選用的閾值函數(shù)簡單、實用性強(qiáng),但其直接將小于閾值的幅值賦零,會損害地震數(shù)據(jù)中的部分弱信號。因此,下一步研究目標(biāo)是推導(dǎo)出更加符合地震數(shù)據(jù)特征的概率分布函數(shù),得到更加復(fù)雜的閾值函數(shù),用于含弱信號地震數(shù)據(jù)的隨機(jī)噪聲壓制。

    參考文獻(xiàn)

    [1]COMMON P.Independent component analysis,a new concept?[J].Signal Processing,1994,36(3):287-314

    [2]劉喜武,劉洪,李幼銘.快速獨立分量變換與去噪初探[J].中國科學(xué)院研究生院學(xué)報,2003,20(4):488-492

    LIU X W,LIU H,LI Y M.Independent component transform and its testing application on seismic noise elimation[J].Journal of the Graduate School of the Chinese Academy of Sciences,2003,20(4):488-492

    [3]彭才,朱仕軍,孫建庫,等.基于獨立成分分析的地震數(shù)據(jù)去噪[J].勘探地球物理進(jìn)展,2007,30(1):30-33

    PENG C,ZHU S J,SUN J K,et al.Seismic denoising based on independent component analysis method[J].Progress in Exploration Geophysics,2007,30(1):30-33

    [4]李大衛(wèi),尹成,謝兵.模擬退火獨立分量分析方法及其應(yīng)用[J].石油物探,2007,46(1):24-27

    LI D W,YIN C,XIE B.Simulated annealing independent component analysis method and its application[J].Geophysical Prospecting for Petroleum,2007,46(1):24-27

    [5]張銀雪,田學(xué)民.基于改進(jìn)PSO-ICA的地震信號去噪方法[J].石油地球物理勘探,2012,47(1):56-62

    ZHANG Y X,TIAN X M.Seismic denoising based on improved PSO-ICA method[J].Oil Geophysical Prospecting,2012,47(1):56-62

    [6]呂文彪,尹成,張白林,等.利用獨立分量分析法去除地震噪聲[J].石油地球物理勘探,2007,42(2):132-136

    LV W B,YIN C,ZHANG B L,et al.Seismic denoising using independent component analysis method[J].Oil Geophysical Prospecting,2007,42(2):132-136

    [7]郭科,陳聆,陳輝,等.基于JADE算法的地震資料隨機(jī)噪音盲分離方法研究[J].地學(xué)前緣,2011,18(3):302-309

    GUO K,CHEN L,CHEN H,et al.The method of seismic random noise blind separation based on JADE[J].Earth Science Frontiers,2011,18(3):302-309

    [8]王維強(qiáng),楊國權(quán).基于EMD與ICA的地震信號去噪技術(shù)研究[J].石油物探,2012,51(1):19-29

    WANG W Q,YANG G Q.The study of seismic denoising based on EMD and ICA[J].Geophysical Prospecting for Petroleum,2012,51(1):19-29

    [10]劉杰.基于獨立分量分析的地震盲反褶積方法及應(yīng)用研究[D].山東東營:中國石油大學(xué)(華東),2008

    LIU J.Study on independent component analysis- based seismic blind deconvolution method and its application[D].Dongying:China University of Petroleum (East China),2008

    [13]芮挺,沈春林,TIAN Q,等.ICA與PCA特征提取能力的比較分析[J].模式識別與人工智能,2005,18(1):124-128

    RUI T,SHEN C L,TIAN Q,et al.Comparison and analysis on ICA & PCA's ability in feature extraction[J].Pattern Recognition and Artificial Intelligence,2005,18(1):124-128

    [14]HANSON M A.Invexity and the Kuhn-Tucker theorem[J].Journal of Mathematical Analysis and Application,1999,236(2):594-604

    [15]劉琚,孫健德,徐宏吉.盲信號處理理論與應(yīng)用[M].北京:科學(xué)出版社,2013:203-208

    LIU J,SUN J D,XU H J.Theory and application of blind signal processing[M].Beijing:Science Press,2013:203-208

    [16]NIKOLAOS M,STANTHAKI T.Optimal contrast correction for ICA-based fusion of multimodal images[J].IEEE Sensors Journal,2008,6(1):1-21

    [17]SANCHIS C,HANSSEN A.Sparse code shrinkage for signal enhancement of seismic data[J].Geophysics,2011,76(6):V151-V167

    [18]HOYER P.Independent compent analysis in image denoising[D].Espoo:Helsinki University of Technology,1999

    [19]SENDUR,SELESNICK I W.Bivariate shrinkage functions for wavelet-based denoising exploiting interscale dependency[J].IEEE Transactions on Signal Processing,2002,50(11):2744-2756

    [20]張恒磊,劉天佑.Curvelet 域蒙特卡羅估計的隨機(jī)噪聲衰減[J].西南石油大學(xué)學(xué)報(自然科學(xué)版),2011,33(4):64-68

    ZHANG H L,LIU T Y.Seismic random noise attenuation via Monte Carlo estimator in curvelet domain[J].Journal of Southwest Petroleum University(Science & Technology Edition),2011,33(4):64-68

    [21]張軍華,藏勝濤,周振曉,等.地震資料信噪比定量計算及方法比較[J].石油地球物理勘探,2009,44(4):481-486

    ZHANG J H,ZANG S T,ZHOU Z X,et al.Quantitative computation and comparison of S/N ratio in seismic data[J].Oil Geophysical Prospecting,2009,44(4):481-486

    (編輯:戴春秋)

    Seismic random noise suppression based on independent component analysis basis functions

    SUN Chengyu,SHAO Jie,LAN Yang,TANG Jie

    (SchoolofGeosciences,ChinaUniversityofPetroleum,Qingdao266580,China)

    Abstract:Traditional seismic denoising methods based on independent component analysis (ICA) assume that adjacent traces in seismic record contain the same random noise.The denoising effect based on this assumption is not significant,and it is only suitable for seismic records with flat event.Therefore,we convert the seismic data to ICA domain by constructing an ICA basis according to the objective function measuring the non-Gaussianity and set a threshold to reduce the random noise.This threshold function is obtained by Bayesian method.In order to satisfy the assumption of independent component analysis,we carried out the partition processing on the seismic data and supposed that each data block has the similar structures with the whole data.The processing results of theoretical and real data show that this method can suppress the random noise effectively and it is suitable for seismic data with complex interface.So the denoising method based on independent component analysis has the obvious advantage and good application value for Gaussian random noise suppression.

    Keywords:independent component analysis,independent component analysis (ICA) basis,Gaussian random noise suppression,Bayesian threshold function

    文章編號:1000-1441(2016)02-0196-09

    DOI:10.3969/j.issn.1000-1441.2016.02.005

    中圖分類號:P631

    文獻(xiàn)標(biāo)識碼:A

    基金項目:國家自然科學(xué)基金(41374123)和國家科技重大專項(2011ZX05006-002)項目共同資助。

    作者簡介:孫成禹(1968—),男,教授,主要從事地震勘探理論和方法研究。

    收稿日期:2015-06-03;改回日期:2015-08-10。

    This research is financially supported by the National Science Foundation of China (Grant No.41374123) and the National Science and Technology Major Project of China (Grant No.2011ZX05006-002).

    18+在线观看网站| 国产乱人视频| 亚洲国产日韩一区二区| 秋霞伦理黄片| 久久精品国产鲁丝片午夜精品| 熟女人妻精品中文字幕| 婷婷色综合www| 深爱激情五月婷婷| 一级毛片 在线播放| 日韩精品有码人妻一区| 91aial.com中文字幕在线观看| 边亲边吃奶的免费视频| 美女脱内裤让男人舔精品视频| 一个人看视频在线观看www免费| 成人美女网站在线观看视频| 中文字幕人妻熟人妻熟丝袜美| 极品少妇高潮喷水抽搐| 日本黄色片子视频| 久久久久国产精品人妻一区二区| 免费观看的影片在线观看| 亚洲综合色惰| av国产免费在线观看| 国产黄色免费在线视频| 九九在线视频观看精品| 波野结衣二区三区在线| 色综合色国产| 国产探花极品一区二区| 色5月婷婷丁香| 十八禁网站网址无遮挡 | 少妇的逼好多水| 国产精品国产三级国产av玫瑰| 身体一侧抽搐| 成人免费观看视频高清| 欧美激情极品国产一区二区三区 | 蜜桃在线观看..| 熟女av电影| 久热这里只有精品99| 中文字幕人妻熟人妻熟丝袜美| 精品亚洲成a人片在线观看 | 一区在线观看完整版| 国产精品成人在线| 国产女主播在线喷水免费视频网站| 国产女主播在线喷水免费视频网站| 天天躁夜夜躁狠狠久久av| 1000部很黄的大片| 两个人的视频大全免费| 三级国产精品欧美在线观看| 国产精品久久久久久久电影| .国产精品久久| 波野结衣二区三区在线| 久久国产精品大桥未久av | 国产欧美日韩精品一区二区| 在线观看人妻少妇| 国产乱人偷精品视频| 国产精品蜜桃在线观看| 日本色播在线视频| 熟女电影av网| 一级毛片黄色毛片免费观看视频| 人妻制服诱惑在线中文字幕| 亚洲av二区三区四区| 国产一区二区三区综合在线观看 | 99国产精品免费福利视频| 国产精品99久久久久久久久| 日韩成人伦理影院| 亚洲最大成人中文| h日本视频在线播放| 中文天堂在线官网| 汤姆久久久久久久影院中文字幕| 精品久久久久久久久亚洲| 日本vs欧美在线观看视频 | 身体一侧抽搐| 两个人的视频大全免费| 亚洲av中文av极速乱| 日韩av不卡免费在线播放| 亚洲精品乱码久久久久久按摩| 2022亚洲国产成人精品| 国产午夜精品久久久久久一区二区三区| 久久99精品国语久久久| 久久精品国产自在天天线| 菩萨蛮人人尽说江南好唐韦庄| 免费不卡的大黄色大毛片视频在线观看| 国产亚洲欧美精品永久| 亚洲内射少妇av| 国产精品欧美亚洲77777| 亚洲精品久久午夜乱码| 毛片一级片免费看久久久久| 午夜老司机福利剧场| 国产精品久久久久久精品电影小说 | 人人妻人人添人人爽欧美一区卜 | 婷婷色av中文字幕| 嘟嘟电影网在线观看| 日日啪夜夜爽| 久久久久久久亚洲中文字幕| 国产色婷婷99| 欧美老熟妇乱子伦牲交| 亚洲欧美日韩东京热| 女性生殖器流出的白浆| 在线天堂最新版资源| 日韩欧美 国产精品| av在线观看视频网站免费| 亚洲av成人精品一二三区| 久久久成人免费电影| 有码 亚洲区| 成人影院久久| 日本av手机在线免费观看| 国产高清有码在线观看视频| 精品少妇久久久久久888优播| 99久久综合免费| 乱系列少妇在线播放| av一本久久久久| 观看美女的网站| 嫩草影院入口| 国产精品熟女久久久久浪| 麻豆精品久久久久久蜜桃| 国产亚洲一区二区精品| 久久久久久久国产电影| 精品亚洲成a人片在线观看 | 精品久久久久久电影网| 爱豆传媒免费全集在线观看| 久久ye,这里只有精品| 毛片一级片免费看久久久久| 久久午夜福利片| 婷婷色麻豆天堂久久| 国产精品国产三级专区第一集| 精华霜和精华液先用哪个| 亚洲国产av新网站| 丰满人妻一区二区三区视频av| 97在线视频观看| 青春草亚洲视频在线观看| 午夜免费鲁丝| 久久6这里有精品| 18禁在线无遮挡免费观看视频| 18禁裸乳无遮挡免费网站照片| 一级毛片我不卡| 黄色视频在线播放观看不卡| 亚洲综合精品二区| 99久国产av精品国产电影| 国产成人aa在线观看| 99热这里只有是精品在线观看| 色视频www国产| 亚洲av不卡在线观看| 国产久久久一区二区三区| 国产高清三级在线| 久久久欧美国产精品| 日韩国内少妇激情av| 夜夜看夜夜爽夜夜摸| 夜夜看夜夜爽夜夜摸| 大话2 男鬼变身卡| 日本-黄色视频高清免费观看| 亚洲精品自拍成人| av在线观看视频网站免费| 在线观看一区二区三区| 日韩中字成人| 国产精品国产三级国产av玫瑰| 免费久久久久久久精品成人欧美视频 | 亚洲精品456在线播放app| 久久久色成人| 久久久久性生活片| av线在线观看网站| 日韩电影二区| 99国产精品免费福利视频| 99视频精品全部免费 在线| 男的添女的下面高潮视频| 国产午夜精品久久久久久一区二区三区| 3wmmmm亚洲av在线观看| 中国三级夫妇交换| 国产男人的电影天堂91| 观看美女的网站| 精品久久久久久久末码| 国产精品一区www在线观看| 在线亚洲精品国产二区图片欧美 | 免费观看在线日韩| 午夜福利高清视频| 中文欧美无线码| 爱豆传媒免费全集在线观看| 黑人猛操日本美女一级片| 亚洲av电影在线观看一区二区三区| 久久国产亚洲av麻豆专区| 久久国内精品自在自线图片| 国产精品一二三区在线看| 蜜桃久久精品国产亚洲av| 免费人成在线观看视频色| 色综合色国产| 亚洲电影在线观看av| 国产精品国产三级专区第一集| 蜜臀久久99精品久久宅男| 精品久久久久久久久亚洲| 亚洲国产精品专区欧美| 看非洲黑人一级黄片| 大片电影免费在线观看免费| 观看美女的网站| 插阴视频在线观看视频| 国产乱来视频区| 亚洲成人一二三区av| 国产免费又黄又爽又色| 人人妻人人看人人澡| 国产伦在线观看视频一区| 99精国产麻豆久久婷婷| 日本-黄色视频高清免费观看| 国产亚洲5aaaaa淫片| 99热国产这里只有精品6| 成人漫画全彩无遮挡| 一本久久精品| 男人狂女人下面高潮的视频| 联通29元200g的流量卡| 亚洲色图av天堂| freevideosex欧美| 国产精品免费大片| 能在线免费看毛片的网站| av国产久精品久网站免费入址| 亚洲人与动物交配视频| 日日啪夜夜爽| 免费看光身美女| 久久久亚洲精品成人影院| 色吧在线观看| 久久综合国产亚洲精品| 成年人午夜在线观看视频| 高清毛片免费看| 欧美激情国产日韩精品一区| 边亲边吃奶的免费视频| 精品国产乱码久久久久久小说| 亚洲精品国产色婷婷电影| 99久久精品国产国产毛片| 中文字幕亚洲精品专区| 91狼人影院| 免费大片黄手机在线观看| 亚洲av日韩在线播放| 草草在线视频免费看| 精品一区在线观看国产| 亚洲va在线va天堂va国产| 欧美成人精品欧美一级黄| 成人一区二区视频在线观看| 人体艺术视频欧美日本| 观看av在线不卡| 久久影院123| 免费黄色在线免费观看| av在线老鸭窝| 久久久久久久久久人人人人人人| 国产极品天堂在线| 一级片'在线观看视频| 美女cb高潮喷水在线观看| 狠狠精品人妻久久久久久综合| 亚洲精华国产精华液的使用体验| 国产免费一级a男人的天堂| 高清日韩中文字幕在线| 国产成人免费无遮挡视频| 国产精品蜜桃在线观看| 人人妻人人爽人人添夜夜欢视频 | 蜜桃亚洲精品一区二区三区| 97超碰精品成人国产| 亚洲精品乱码久久久久久按摩| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲综合精品二区| 国产一级毛片在线| 欧美高清成人免费视频www| 日韩欧美 国产精品| 自拍欧美九色日韩亚洲蝌蚪91 | 色视频在线一区二区三区| 少妇熟女欧美另类| 成年女人在线观看亚洲视频| 亚洲欧美精品自产自拍| 亚洲怡红院男人天堂| 国产免费视频播放在线视频| 香蕉精品网在线| 国产在线男女| 精品人妻偷拍中文字幕| 又粗又硬又长又爽又黄的视频| 高清欧美精品videossex| 老师上课跳d突然被开到最大视频| 不卡视频在线观看欧美| 亚洲精品亚洲一区二区| 国产精品欧美亚洲77777| 亚洲va在线va天堂va国产| 精品亚洲成a人片在线观看 | 我的女老师完整版在线观看| 啦啦啦中文免费视频观看日本| 欧美丝袜亚洲另类| 99久久精品一区二区三区| 人人妻人人澡人人爽人人夜夜| 少妇被粗大猛烈的视频| 色视频在线一区二区三区| 超碰97精品在线观看| kizo精华| 亚洲精品自拍成人| 亚洲精品国产色婷婷电影| 国产 精品1| 精品久久久噜噜| 免费观看av网站的网址| 日韩制服骚丝袜av| 日韩成人av中文字幕在线观看| 51国产日韩欧美| 亚洲av中文字字幕乱码综合| 欧美日韩综合久久久久久| 精品久久久久久久久亚洲| 国产高清国产精品国产三级 | 三级经典国产精品| 成人黄色视频免费在线看| 中国国产av一级| 水蜜桃什么品种好| 国产色婷婷99| 国产一区二区三区综合在线观看 | 七月丁香在线播放| 干丝袜人妻中文字幕| 六月丁香七月| 麻豆精品久久久久久蜜桃| 国产又色又爽无遮挡免| 寂寞人妻少妇视频99o| 777米奇影视久久| 国内精品宾馆在线| 又爽又黄a免费视频| 网址你懂的国产日韩在线| 街头女战士在线观看网站| 97在线视频观看| 一本一本综合久久| 高清毛片免费看| 夫妻性生交免费视频一级片| 亚洲av.av天堂| 啦啦啦视频在线资源免费观看| 国产精品蜜桃在线观看| 国产精品.久久久| 黑人高潮一二区| 99久国产av精品国产电影| 亚洲国产精品国产精品| 一级毛片黄色毛片免费观看视频| 久久久久久伊人网av| 国语对白做爰xxxⅹ性视频网站| 最近最新中文字幕大全电影3| 成年美女黄网站色视频大全免费 | 国产91av在线免费观看| 大香蕉久久网| 美女主播在线视频| 在线观看免费高清a一片| 国产精品蜜桃在线观看| 久久久久久久久久人人人人人人| 少妇人妻久久综合中文| 天天躁日日操中文字幕| 国产日韩欧美亚洲二区| 国产伦理片在线播放av一区| 亚洲成人中文字幕在线播放| 啦啦啦视频在线资源免费观看| av线在线观看网站| 99视频精品全部免费 在线| 中国三级夫妇交换| 久久久色成人| 国产伦精品一区二区三区四那| 永久免费av网站大全| 久久久精品94久久精品| 青春草国产在线视频| 黄色配什么色好看| 在线观看一区二区三区| 亚洲婷婷狠狠爱综合网| 国产成人a区在线观看| 三级国产精品欧美在线观看| 九色成人免费人妻av| 少妇 在线观看| 亚洲性久久影院| 在线观看一区二区三区| 亚洲婷婷狠狠爱综合网| 国产精品蜜桃在线观看| 丝瓜视频免费看黄片| 国产精品99久久久久久久久| 国内揄拍国产精品人妻在线| 性色avwww在线观看| 亚洲第一av免费看| 22中文网久久字幕| 成人影院久久| 超碰97精品在线观看| www.av在线官网国产| 国产成人免费观看mmmm| 日韩 亚洲 欧美在线| 岛国毛片在线播放| 99热网站在线观看| a 毛片基地| 91狼人影院| 3wmmmm亚洲av在线观看| 日本欧美国产在线视频| 在线播放无遮挡| 又大又黄又爽视频免费| 一区二区三区免费毛片| 大又大粗又爽又黄少妇毛片口| 亚洲激情五月婷婷啪啪| 人人妻人人添人人爽欧美一区卜 | 一个人免费看片子| 欧美精品一区二区大全| 亚洲av电影在线观看一区二区三区| 亚洲国产成人一精品久久久| 美女脱内裤让男人舔精品视频| 一级片'在线观看视频| 妹子高潮喷水视频| 日本黄色日本黄色录像| 国产一区亚洲一区在线观看| 建设人人有责人人尽责人人享有的 | 精品一区二区三卡| 日本与韩国留学比较| 午夜老司机福利剧场| 久久久精品免费免费高清| 又大又黄又爽视频免费| 高清欧美精品videossex| 青春草视频在线免费观看| 亚洲精品久久久久久婷婷小说| a级毛片免费高清观看在线播放| 免费黄网站久久成人精品| 亚洲第一区二区三区不卡| 成人漫画全彩无遮挡| 女的被弄到高潮叫床怎么办| 成年美女黄网站色视频大全免费 | 国产淫语在线视频| 久久99蜜桃精品久久| 97热精品久久久久久| 小蜜桃在线观看免费完整版高清| 少妇人妻一区二区三区视频| 精品久久久久久久末码| 国产午夜精品一二区理论片| 婷婷色综合大香蕉| 午夜老司机福利剧场| 国产淫语在线视频| 成人黄色视频免费在线看| 国产精品久久久久成人av| 亚洲,欧美,日韩| 久久久久久伊人网av| 少妇熟女欧美另类| 亚洲成人av在线免费| 国产亚洲午夜精品一区二区久久| 超碰97精品在线观看| 国产有黄有色有爽视频| 欧美少妇被猛烈插入视频| 亚洲精品456在线播放app| 波野结衣二区三区在线| 亚洲一级一片aⅴ在线观看| 少妇 在线观看| 老司机影院毛片| 麻豆国产97在线/欧美| 亚洲av综合色区一区| 国产精品久久久久久久电影| 国产亚洲欧美精品永久| 亚洲av综合色区一区| 久久国产精品大桥未久av | 日韩中文字幕视频在线看片 | 男的添女的下面高潮视频| 久久99热这里只频精品6学生| 涩涩av久久男人的天堂| 国产精品成人在线| 热99国产精品久久久久久7| 亚洲最大成人中文| 国产亚洲午夜精品一区二区久久| 国产爽快片一区二区三区| 永久免费av网站大全| 国产久久久一区二区三区| 国产亚洲av片在线观看秒播厂| 边亲边吃奶的免费视频| 亚洲av电影在线观看一区二区三区| 老熟女久久久| 在线观看av片永久免费下载| 高清不卡的av网站| 大又大粗又爽又黄少妇毛片口| 少妇人妻一区二区三区视频| 男人和女人高潮做爰伦理| 久久这里有精品视频免费| 久久精品久久久久久噜噜老黄| 2021少妇久久久久久久久久久| 十八禁网站网址无遮挡 | 欧美激情极品国产一区二区三区 | 久久精品国产亚洲网站| 交换朋友夫妻互换小说| 免费观看的影片在线观看| 国产探花极品一区二区| 欧美性感艳星| 欧美高清成人免费视频www| 有码 亚洲区| 成人美女网站在线观看视频| 色婷婷av一区二区三区视频| 亚洲人成网站在线播| 国产真实伦视频高清在线观看| 国产中年淑女户外野战色| 大话2 男鬼变身卡| 夜夜骑夜夜射夜夜干| 国产高清不卡午夜福利| 91久久精品国产一区二区成人| 色5月婷婷丁香| 亚洲精品色激情综合| 中文字幕亚洲精品专区| 男的添女的下面高潮视频| 少妇猛男粗大的猛烈进出视频| 午夜免费男女啪啪视频观看| 亚洲国产av新网站| 国产视频内射| 两个人的视频大全免费| 六月丁香七月| 日韩在线高清观看一区二区三区| 99久久中文字幕三级久久日本| 99热6这里只有精品| 免费观看av网站的网址| 久久精品久久精品一区二区三区| 精品一区二区三卡| 一区二区av电影网| 亚洲精品第二区| 永久免费av网站大全| 91精品国产国语对白视频| 亚洲欧美成人精品一区二区| 大片电影免费在线观看免费| 免费黄网站久久成人精品| 午夜日本视频在线| 国产老妇伦熟女老妇高清| 欧美激情极品国产一区二区三区 | 亚洲av男天堂| 91精品国产九色| 2022亚洲国产成人精品| 国产老妇伦熟女老妇高清| 大陆偷拍与自拍| 国产成人免费无遮挡视频| 久久毛片免费看一区二区三区| 亚洲av福利一区| 亚洲av不卡在线观看| 国内少妇人妻偷人精品xxx网站| 美女内射精品一级片tv| 欧美+日韩+精品| 一级毛片电影观看| 精品久久久噜噜| 我的女老师完整版在线观看| 久久精品国产鲁丝片午夜精品| 精华霜和精华液先用哪个| 国产乱来视频区| 精品人妻偷拍中文字幕| 亚洲av免费高清在线观看| 夫妻性生交免费视频一级片| 国产黄片视频在线免费观看| 联通29元200g的流量卡| 婷婷色综合大香蕉| 日韩一区二区三区影片| 精品人妻偷拍中文字幕| 成年人午夜在线观看视频| 九九在线视频观看精品| 91久久精品国产一区二区成人| 97超碰精品成人国产| 夫妻午夜视频| 美女中出高潮动态图| 日韩制服骚丝袜av| 午夜精品国产一区二区电影| 日本色播在线视频| 欧美精品亚洲一区二区| 男女免费视频国产| 中文字幕精品免费在线观看视频 | 久久国产精品大桥未久av | 久久青草综合色| 伦精品一区二区三区| 青春草亚洲视频在线观看| 午夜视频国产福利| 在线观看三级黄色| 日韩视频在线欧美| xxx大片免费视频| 99热这里只有是精品50| 色吧在线观看| videos熟女内射| 国产精品秋霞免费鲁丝片| 精品99又大又爽又粗少妇毛片| 亚洲精品,欧美精品| 少妇 在线观看| 亚洲精品视频女| 国产视频内射| 国产一区二区三区av在线| 久热久热在线精品观看| 国产真实伦视频高清在线观看| 国产精品国产三级国产专区5o| av在线观看视频网站免费| 美女高潮的动态| 亚洲精品,欧美精品| 日本黄色片子视频| 日本与韩国留学比较| 国产v大片淫在线免费观看| 国产一区二区在线观看日韩| 国产v大片淫在线免费观看| 欧美精品一区二区免费开放| 嘟嘟电影网在线观看| 男女下面进入的视频免费午夜| 欧美高清性xxxxhd video| 99热网站在线观看| 久久久a久久爽久久v久久| 久久久久视频综合| 在线观看一区二区三区激情| 自拍欧美九色日韩亚洲蝌蚪91 | videos熟女内射| 国产精品人妻久久久影院| 久热这里只有精品99| 激情五月婷婷亚洲| 国产成人一区二区在线| 亚洲,一卡二卡三卡| 成人亚洲欧美一区二区av| 亚洲av免费高清在线观看| 中文天堂在线官网| 99热这里只有是精品在线观看| 亚洲欧美精品专区久久| 亚洲av不卡在线观看| 日本vs欧美在线观看视频 | 欧美一级a爱片免费观看看| 黑人猛操日本美女一级片| 欧美+日韩+精品| 少妇猛男粗大的猛烈进出视频| 欧美日韩一区二区视频在线观看视频在线| 在线 av 中文字幕| 中文字幕av成人在线电影| 一级黄片播放器| 亚洲熟女精品中文字幕| 丝瓜视频免费看黄片| 99久久中文字幕三级久久日本| 日韩电影二区| 亚洲国产精品一区三区| 女性生殖器流出的白浆| 99久国产av精品国产电影| 中文字幕av成人在线电影| 亚洲人与动物交配视频| 欧美精品一区二区免费开放| 成年av动漫网址| 少妇丰满av| 一个人免费看片子| 99热国产这里只有精品6|