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

    一種自適應(yīng)增益限的反Q濾波

    2015-03-16 11:06:54張固瀾林進(jìn)王熙明賀振華曹俊興張建軍賀錫雷林凱薛雅娟
    地球物理學(xué)報 2015年7期
    關(guān)鍵詞:子波頻帶振幅

    張固瀾, 林進(jìn), 王熙明, 賀振華, 曹俊興,張建軍, 賀錫雷, 林凱, 薛雅娟

    1 成都理工大學(xué)油氣藏地質(zhì)及開發(fā)工程國家重點實驗室, 成都 610059 2 中國石油集團(tuán)東方地球物理公司新興物探開發(fā)處, 河北涿州 072751 3 陜西延長石油(集團(tuán))有限責(zé)任公司油氣勘探公司, 陜西延安 716000 4 成都信息工程大學(xué)通信工程學(xué)院, 成都 610225

    ?

    一種自適應(yīng)增益限的反Q濾波

    張固瀾1,2, 林進(jìn)3, 王熙明2, 賀振華1, 曹俊興1,張建軍2, 賀錫雷1, 林凱1, 薛雅娟4

    1 成都理工大學(xué)油氣藏地質(zhì)及開發(fā)工程國家重點實驗室, 成都 610059 2 中國石油集團(tuán)東方地球物理公司新興物探開發(fā)處, 河北涿州 072751 3 陜西延長石油(集團(tuán))有限責(zé)任公司油氣勘探公司, 陜西延安 716000 4 成都信息工程大學(xué)通信工程學(xué)院, 成都 610225

    地層的Q吸收會造成地震波振幅衰減、相位畸變,分辨率和信噪比明顯降低.反Q濾波可消除由于地層Q吸收造成的振幅衰減和相位畸變,從而提高地震資料的分辨率;但反Q濾波振幅補(bǔ)償?shù)臄?shù)值不穩(wěn)定性問題會嚴(yán)重降低地震資料的信噪比,并產(chǎn)生很多假象.截止頻率法和穩(wěn)定因子法反Q濾波振幅補(bǔ)償方法雖可控制數(shù)值非穩(wěn)定性問題,但振幅補(bǔ)償函數(shù)的增益限為一個時不變的常數(shù),且與地震數(shù)據(jù)動態(tài)范圍無關(guān),其經(jīng)常會壓制深層地震波的高頻成分,反而降低地震資料的分辨率;因此,本文在研究截止頻率法和穩(wěn)定因子法的基礎(chǔ)上,結(jié)合地震數(shù)據(jù)的動態(tài)范圍對地震記錄分辨率的影響,提出了一種自適應(yīng)增益限的反Q濾波振幅補(bǔ)償方法,其增益限和穩(wěn)定因子都是時變的,且都自適應(yīng)于地震數(shù)據(jù)有效頻帶的截止頻率.合成數(shù)據(jù)和實際數(shù)據(jù)試算表明,本文的自適應(yīng)增益限的反Q濾波方法可恢復(fù)地震信號有效頻帶范圍內(nèi)的能量,且能較好地控制數(shù)值非穩(wěn)定性問題,最終獲得高分辨率和高信噪比的地震數(shù)據(jù).

    反Q濾波; 動態(tài)范圍; 自適應(yīng); 增益限; 穩(wěn)定因子

    1 引言

    地震波在傳播過程中要經(jīng)受與地層有關(guān)的非彈性衰減,因此隨傳播時間的增大,地震波振幅衰減、相位畸變(Futterman, 1962; Kjartansson, 1979; Wang and Guo, 2004),分辨率和信噪比明顯降低.造成地震波非彈性衰減主要有兩方面的原因:一方面是諧振散射,其造成低頻陰影現(xiàn)象;另一方面是地層的Q吸收效應(yīng).

    低頻陰影現(xiàn)象的主要表現(xiàn)是在油氣儲集層下方出現(xiàn)低頻的強(qiáng)能量,從而可直接指示油氣儲層.在國外,Castagna等(2003)利用匹配追蹤的時頻分解方法,展示了低頻陰影在油氣預(yù)測中的良好應(yīng)用前景.國內(nèi),賀振華等(2008)利用基于彌散黏滯理論模型來解釋低頻陰影現(xiàn)象,并利用廣義S變換進(jìn)行油氣檢測,取得理想效果.由于地震波低頻信號的分辨率畢竟有限,因此低頻陰影現(xiàn)象無法準(zhǔn)確指示含油氣的位置.

    反Q濾波(Hale, 1981, 1982; Bickel and Natarajan, 1985;Hargreaves and Calvert,1991; Bano,1996; Wang,2002,2003; 姚振興等,2003;Wang Y H,2006; Wang H Y,2008; 劉財?shù)?2013;陳增保等,2014;張固瀾等,2014)可補(bǔ)償?shù)貙観吸收造成的能量損失,提高地震記錄主頻,拓寬有效頻帶,同時可消除子波時變影響,從而最終提高地震資料的分辨率.反Q濾波振幅補(bǔ)償后,地震記錄主要能量逐漸向高頻方向移動,若對反Q濾波后的地震數(shù)據(jù)進(jìn)行時頻分析,則低頻陰影現(xiàn)象對油氣預(yù)測的分辨率就會提高.

    受硬件存儲精度的影響,反Q濾波振幅補(bǔ)償后數(shù)據(jù)會出現(xiàn)數(shù)值穩(wěn)定性問題.為控制反Q濾波的數(shù)值非穩(wěn)定性問題,有關(guān)學(xué)者提出了增益限(反Q濾波振幅補(bǔ)償函數(shù)的最大值)約束的截止頻率法(Bickel and Natarajan, 1985)和穩(wěn)定因子法(James and Knight,2003;Wang, 2006)用于反Q濾波振幅補(bǔ)償.但這些方法中,振幅補(bǔ)償函數(shù)的增益限時不變,并未與地震數(shù)據(jù)的動態(tài)范圍及其有效頻帶結(jié)合;當(dāng)增益限太小,反Q濾波雖可控制數(shù)值非穩(wěn)定性問題,但深層地震波的高頻能量反而被壓制,因此地震記錄分辨率反而降低;當(dāng)增益限太大,反Q濾波雖可提高分辨率,但由于硬件存儲精度造成的誤差會被明顯放大,因此信噪比降低明顯.綜上,固定增益限約束的反Q濾波振幅補(bǔ)償方法很難取得令人滿意的效果.

    本文基于地震數(shù)據(jù)的動態(tài)范圍,結(jié)合固定增益限約束的截止頻率法和穩(wěn)定因子法反Q濾波,提出一種自適應(yīng)的反Q濾波振幅補(bǔ)償方法.該方法中振幅補(bǔ)償函數(shù)的增益限和穩(wěn)定因子都為時變參數(shù),且與各時刻地震波有效頻帶范圍的截止頻率有關(guān);因此,其可以恢復(fù)地震波有效頻帶范圍內(nèi)的能量,并較好壓制高頻噪聲,從而獲得高分辨率和高信噪比的地震數(shù)據(jù).

    2 反Q濾波理論基礎(chǔ)

    2.1 反Q濾波原理

    基于Q與頻率無關(guān)的假設(shè)和一維雙程傳播波動方程:

    (1)

    U(t+Δt,ω)=

    (2)

    地層Q吸收對應(yīng)的振幅衰減函數(shù)為:

    (3)

    反Q濾波延拓方程可表示為:

    U(t+Δt,ω)=

    (4)

    反Q濾波振幅補(bǔ)償函數(shù)為:

    (5)

    2.2 反Q濾波的數(shù)值非穩(wěn)定性問題

    圖1a為合成數(shù)據(jù);圖1b為按方程(2)進(jìn)行Q吸收正演后的結(jié)果;圖1c按方程(4)進(jìn)行反Q濾波的結(jié)果.Q吸收正演和反Q濾波,所用的Q相同,參考角頻率也相同.

    圖1c中,由于相位校正不存在穩(wěn)定性問題,因此反Q濾波后地震子波恢復(fù)成零相位子波;但反Q濾波振幅補(bǔ)償后出現(xiàn)不穩(wěn)定性問題,1 s后高頻噪聲明顯,信噪比和分辨率都嚴(yán)重降低.出現(xiàn)這種強(qiáng)噪聲的原因有:首先,在地震數(shù)據(jù)野外采集過程中,地震波在地層中傳播時振幅衰減很快,當(dāng)超過一定的傳播時間,地震波的振幅已低于環(huán)境噪聲.其次,設(shè)備的背景噪聲.對野外采集的實際地震數(shù)據(jù)而言,背景噪聲是用于地震數(shù)據(jù)采集的檢波器的存儲精度造成的誤差;對本文的無噪聲合成數(shù)據(jù)而言,背景噪聲是用于地震數(shù)據(jù)處理的計算機(jī)的存儲精度造成的誤差(Wang,2002).反Q濾波振幅補(bǔ)償恢復(fù)了地震信號,但也同樣放大了環(huán)境噪聲、設(shè)備背景噪聲.我們把這種放大噪聲的現(xiàn)象稱為反Q濾波的數(shù)值非穩(wěn)定性問題.

    2.3 地震數(shù)據(jù)的動態(tài)范圍

    圖1 反Q濾波前后效果對比Fig.1 The result of synthetic tarce before and after inverse Q-filter

    (6)

    這里,動態(tài)范圍的單位是dB;Up(t,ω)和Ul(t,ω)都是時變的,不同時刻對應(yīng)不同值.基于對動態(tài)范圍的定義,若t時刻地震波的振幅譜小于Ul(t,ω),它則以存儲設(shè)備背景噪聲H(t,ω)來呈現(xiàn),因此:

    (7)

    定義ωp(t)為地震波振幅譜的峰值角頻率,則U(t,ωp(t))=Up(t,ω);定義ωd(t)為地震波有效頻帶的截止角頻率,則:U(t,ωd(t))=Ul(t,ω),且ωd(t)>ωp(t).因此,方程(7)也可表示為:

    (8)

    考慮到地震波的Q吸收效應(yīng),ωp(t)和ωd(t)不僅與品質(zhì)因子Q、傳播時間t、震源子波的振幅譜U(t,ω)有關(guān),還與設(shè)備可準(zhǔn)確存儲的數(shù)據(jù)動態(tài)范圍d有關(guān).

    為直觀說明設(shè)備存儲精度對地震數(shù)據(jù)的動態(tài)范圍的影響,對峰值頻率為50Hz的Ricker子波(圖2a)進(jìn)行傅里葉變換,其振幅譜如圖2b細(xì)線;利用方程(10),可得其理論上的振幅譜,如圖2b粗線.Ricker子波的表達(dá)式為:

    (9)

    Ricker子波振幅譜理論上的表達(dá)式為:

    (10)

    其中,fp為Ricker子波的峰值頻率.

    由圖2b可知,當(dāng)頻率為222 Hz時,Ricker子波的振幅譜為-140 dB;由于計算機(jī)的存儲精度和傅里葉變換的截斷誤差造成的影響,當(dāng)Ricker子波振幅譜小于-140 dB時,其已不能被準(zhǔn)確記錄,以背景噪聲的形式呈現(xiàn).由圖2可看出:用于數(shù)值模擬的計算機(jī)可準(zhǔn)確存儲的數(shù)據(jù)動態(tài)范圍為140 dB,且其所記錄下來的Ricker子波有效頻帶為[0,222]Hz,有效頻帶的截止頻率為222 Hz.

    利用Q=50,傳播時間t=2 s,結(jié)合方程(3)和(10),可得經(jīng)Q吸收振幅衰減后的振幅譜,如圖3b粗線;對圖3b的振幅譜進(jìn)行反傅里葉變換,得到經(jīng)Q吸收振幅衰減后的波形(圖3a);對圖3a中的波形進(jìn)行傅里葉變換,可得到與計算機(jī)的存儲精度有關(guān)的振幅譜,如圖3b細(xì)線.當(dāng)振幅譜小于-140 dB,其不能被完全記錄下來,此時,子波有效頻帶為[0,130] Hz,有效頻帶的截止頻率為130 Hz.考慮到地層的Q吸收,地震波有效頻帶的截止頻率隨時間增大而逐漸降低,有效頻寬隨時間的增大逐漸變窄.

    若嚴(yán)格按照方程(5)對經(jīng)過地層吸收衰減后的地震波進(jìn)行反Q濾波振幅補(bǔ)償,就會放大背景噪聲,出現(xiàn)數(shù)值非穩(wěn)定性問題,從而嚴(yán)重降低地震資料的信噪比.

    3 時不變增益限的反Q濾波

    時不變增益限約束的截止頻率法和穩(wěn)定因子法(Wang,2006)反Q濾波振幅補(bǔ)償方法可控制反Q濾波振幅補(bǔ)償?shù)姆€(wěn)定性和信噪比,但增益限為一時不變的常數(shù)c.

    圖2 Ricker子波及其振幅譜Fig.2 The Ricker wavelet and its amplitude spectrum

    圖3 Ricker子波經(jīng)Q吸收衰減后的波形及其振幅譜Fig.3 The Ricker wavelet and its amplitude spectrum afer the earth Q-filter

    3.1 截止頻率法反Q濾波

    截止頻率法反Q濾波振幅補(bǔ)償函數(shù)表達(dá)式為:

    (11)

    其中,Q(t)為傳播時間t處對應(yīng)的地層的品質(zhì)因子;B(t,ωc(t),Q(t))=c;ωc(t)是振幅補(bǔ)償函數(shù)F(t,ω,Q(t))取最大值時對應(yīng)的角頻率,稱之峰值角頻率.基于方程(5),則有:

    (12)

    因此,ωc(t)隨時間增大逐漸減小,且ωc(t)與品質(zhì)因子Q(t)、傳播時間t和增益限c有關(guān).

    由方程(11),當(dāng)ω≤ωc(t)時,地震波可以完全被恢復(fù)成震源子波;當(dāng)ω>ωc(t)時,振幅補(bǔ)償函數(shù)被增益限c所代替,因此可有效壓制高頻噪聲.但由于振幅補(bǔ)償函數(shù)F(t,ω,Q(t))在ωc(t)處被截斷,其會造成補(bǔ)償后的地震波振幅譜出現(xiàn)截斷現(xiàn)象,從而造成嚴(yán)重的吉普斯效應(yīng),很難取得滿意的效果.

    3.2 穩(wěn)定因子法反Q濾波

    穩(wěn)定因子法(Wang,2006)反Q濾波振幅補(bǔ)償方法將振幅補(bǔ)償函數(shù)B(t,ω,Q(t))轉(zhuǎn)化為分?jǐn)?shù)形式,并在分子和分母上都加上一個很小的正數(shù),從而可在截止頻率法基礎(chǔ)上進(jìn)一步壓制高頻噪聲,并有效避免吉普斯效應(yīng).穩(wěn)定因子法反Q濾波振幅補(bǔ)償函數(shù)表達(dá)式為:

    (13)

    (14)

    設(shè)B(t,ω,Q(t))為自變量,并對方程(14)求微分,得:

    (15)

    若假設(shè)ωs1(t)是振幅補(bǔ)償函數(shù)S1(t,ω,Q(t))取最大值時對應(yīng)的角頻率,則

    S1(t,ωs1(t),Q(t))=

    (16)

    (17)

    (18)

    將方程(18)代入方程(16),得:

    S1(t,ωs1,Q(t))=0.5×B(t,ωs1(t),Q(t)),

    (19)

    將方程(19)代入方程(17),得:

    (20)

    將方程(20)代入方程(18),得:

    (21)

    因此,方程(13)可表示為:

    (22)

    由方程(22),當(dāng)ω>ωs1(t),振幅補(bǔ)償函數(shù)S1(t,ω,Q(t))隨著角頻率的增大而逐漸減小;因此,利用穩(wěn)定因子法反Q濾波,可有效避免吉普斯效應(yīng),取得較好效果.

    基于方程(5)和(20), 則有:

    (23)

    因此,ωs1(t)與品質(zhì)因子Q(t)、傳播時間t和增益限c有關(guān).

    對于截止頻率法和穩(wěn)定因子法,若增益限c太大, 基于方程 (12),則ωc(t)?ωd(t);基于方程 (23),則ωs1(t)?ωd(t).由于反Q濾波振幅補(bǔ)償后,背景噪聲和高頻環(huán)境噪聲被明顯放大,因此信噪比明顯降低.若增益限c太小,深層的地震波高頻成分明顯被壓制,地震波分辨率不僅沒有提高,反而會降低.因此,穩(wěn)定因子法的反Q濾波也很難取得令人滿意的效果.

    4 自適應(yīng)增益限的反Q濾波

    反Q濾波振幅補(bǔ)償應(yīng)充分補(bǔ)償?shù)卣鸩ㄓ行ьl帶范圍內(nèi)的能量,并有效壓制有效頻帶范圍外的高頻噪聲,從而最終獲得令人滿意的效果.基于對地震數(shù)據(jù)的動態(tài)范圍的分析,若要完全恢復(fù)地震波有效頻帶范圍ω≤ωd(t)內(nèi)的振幅,可嚴(yán)格按照方程(5)對地震波的振幅譜進(jìn)行補(bǔ)償,即當(dāng)ω≤ωd(t)時,振幅補(bǔ)償函數(shù)為:

    (24)

    為壓制有效頻帶范圍外的高頻噪聲,當(dāng)ω>ωd(t)時,振幅補(bǔ)償函數(shù)S2(t,ω,Q(t))應(yīng)隨頻率的增大而逐漸減??;因此,需要使振幅補(bǔ)償函數(shù)S2[t,ω,Q(t)]的增益限滿足:

    (25)

    由于ωd(t)是時變的,且與地震數(shù)據(jù)的動態(tài)范圍緊密聯(lián)系,因此增益限c(t)是時變的,我們稱時變增益限c(t)自適應(yīng)于地震波有效頻帶的截止角頻率ωd(t)或地震數(shù)據(jù)的動態(tài)范圍.

    當(dāng)ω>ωd(t)時,由于穩(wěn)定因子法的振幅補(bǔ)償函數(shù)S1[t,ω,Q(t)]隨ω增大而逐漸減小,因此可以基于穩(wěn)定因子法補(bǔ)償函數(shù)來壓制高頻噪聲.若使:

    (26)

    將方程(26)代入方程(20),則穩(wěn)定因子法補(bǔ)償函數(shù)S1[t,ω,Q(t)]增益限為c滿足:

    (27)

    (28)

    因此,方程(22)變?yōu)椋?/p>

    (29)

    為保證振幅補(bǔ)償函數(shù)S2[t,ω,Q(t)]在ω=ωd(t)處連續(xù),從而進(jìn)一步降低吉普斯效應(yīng),可使:

    (30)

    相應(yīng)地,方程(29)可寫成:

    (31)

    S2(t,ω,Q(t))=

    (32)

    5 合成記錄試算

    將峰值頻率為50Hz的Ricker子波作為震源子波.由于野外地震數(shù)據(jù)的動態(tài)范圍一般都小于60dB,因此為模擬真實地震數(shù)據(jù)的動態(tài)范圍,假設(shè)合成數(shù)據(jù)的動態(tài)范圍為60dB;當(dāng)歸一化后的振幅譜小于-60dB時,假設(shè)其為背景噪聲H(t,ω):

    H(t,ω)=0.5×10-3×[1+rand(0,1)],

    (33)

    rand(0,1)表示取0和1之間的隨機(jī)數(shù).

    圖4a是合成記錄;圖4b是利用Q=50對圖4a進(jìn)行Q吸收衰減模擬后的結(jié)果.圖5是震源子波及經(jīng)過Q吸收衰減后各時刻的地震子波的振幅譜.由于:ωp(0.3)=80π,ωp(0.6)=64π,ωp(0.9)=52π,ωp(1.2)=44π,ωp(1.5)=38π,ωp(1.8)=32π;ωd(0.3)=284π,ωd(0.6)=254π,ωd(0.9)=224π,ωd(1.2)=204π,ωd(1.5)=180π,ωd(1.8)=162π.因此,隨時間增大,地震波振幅譜的峰值頻率和有效頻帶的截止頻率都減小,有效頻帶逐漸變窄.

    利用穩(wěn)定因子法反Q濾波振幅補(bǔ)償,結(jié)合不同的增益限 (20, 40 和60 dB),對圖 4b的數(shù)據(jù)進(jìn)行反Q濾波.圖6是增益限為20 dB時,穩(wěn)定因子法反Q濾波補(bǔ)償后的振幅譜(圖6a)和反Q濾波后的結(jié)果(圖6b).

    圖4 (a)合成數(shù)據(jù); (b)經(jīng)Q吸收衰減后的結(jié)果Fig.4 (a) The synthetic trace considering the dynamic range; (b) Result of earth Q-filter considering the dynamic range

    當(dāng)增益限為20 dB時,由于:ωs1(0.3)=318π,ωs1(0.6)=158π,ωs1(0.9)=106π,ωs1(1.2)=80π,ωs1(1.5)=62π,ωs1(1.8)=52π. 當(dāng)傳播時間大于0.6s,ωs1(t)<2.8×ωp(t).按Widess(1982)對于零相位子波分辨率的表達(dá)式,我們利用圖5和圖6a中各時刻地震子波不同頻帶范圍內(nèi)的振幅譜(也就是不同動態(tài)范圍的振幅譜,且動態(tài)范圍越大,有效頻帶越寬)計算不同時刻地震子波的分辨率,如圖7所示.圖7中,曲線1代表經(jīng)地層Q吸收衰減后,不同時刻的子波分辨率,曲線2代表當(dāng)增益限為20 dB時,穩(wěn)定因子法反Q濾波振幅補(bǔ)償后不同時刻子波的分辨率.當(dāng)傳播時間大于0.6 s時,反Q濾波后子波的分辨率不僅沒有提高反而降低.

    圖8是增益限為40 dB時,穩(wěn)定因子法反Q濾波補(bǔ)償后的振幅譜(圖8a)和反Q濾波后的結(jié)果(圖8b);

    ωs1(0.3)=558π,ωs1(0.6)=280π,ωs1(0.9)=186π,ωs1(1.2)=140π,ωs1(1.5)=112π,ωs1(1.8)=60π.

    當(dāng)傳播時間小于等于0.3s,ωs1(t)?ωd(t),高頻噪聲明顯放大,信噪比降低;當(dāng)傳播時間大于0.9s,ωs1(t)<ωd(t),高頻噪聲雖未明顯放大,但主瓣旁邊出現(xiàn)很多旁瓣,且這種現(xiàn)象隨時間的增大越來越嚴(yán)重.

    圖9是增益限為60dB時,穩(wěn)定因子法反Q濾波補(bǔ)償后的振幅譜(圖9a)和反Q濾波后的結(jié)果(圖9b);

    ωs1(0.3)=804π,ωs1(0.6)=402π,ωs1(0.9)=268π,ωs1(1.2)=202π,ωs1(1.5)=160π,ωs1(1.8)=134π. 當(dāng)傳播時間小于等于0.9s時,ωs1(t)?ωd(t),高頻噪聲放大明顯,信噪比降低;且末端假頻現(xiàn)象嚴(yán)重.

    利用本文的自適應(yīng)增益限的反Q濾波方法,對圖4b所示的數(shù)據(jù)進(jìn)行反Q濾波補(bǔ)償.當(dāng)我們設(shè)定數(shù)據(jù)動態(tài)范圍為60 dB時,不同時刻的振幅補(bǔ)償函數(shù)具有不同的增益限:

    c(0.3)=23 dB,c(0.6)=41 dB,c(0.9)=55 dB;

    c(1.2)=67 dB,c(1.5)=73 dB,c(1.8)=79 dB.

    圖10是當(dāng)數(shù)據(jù)的動態(tài)范圍為60 dB時,本文的自適應(yīng)增益限的反Q濾波在各時刻的振幅補(bǔ)償函數(shù)(圖10a)和反Q濾波振幅補(bǔ)償后的振幅譜(圖10b):各時刻的地震波振幅譜基本接近峰值頻率為50 Hz的Ricker子波的振幅譜,因此其分辨率與峰值頻率為50 Hz的Ricker子波分辨率相近.

    圖11a是當(dāng)數(shù)據(jù)動態(tài)范圍為60 dB時,利用自適應(yīng)增益限的方法對圖4b所示的數(shù)據(jù)進(jìn)行反Q濾波補(bǔ)償?shù)玫降慕Y(jié)果;圖11b是其與圖9b之間的差值.可以看出,本文的自適應(yīng)增益限的方法得到的結(jié)果要好于穩(wěn)定因子法:高頻能量得到了更好的恢復(fù),且高頻噪聲和末端假頻都得到了很好的控制.

    對圖4b的數(shù)據(jù)加入隨機(jī)噪聲,然后利用穩(wěn)定因子法和自適應(yīng)增益限的反Q濾波方法進(jìn)行反Q濾波補(bǔ)償;噪聲如圖12a所示,且噪聲振幅最大值與地震信號振幅最大值的比值為 0.1%.

    圖12b是增益限為60 dB時,穩(wěn)定因子法反Q濾波補(bǔ)償后的效果,可看到在0.3 s,0.6 s和0.9 s處,地震波的最大振幅譜都大于1.0,其破壞了振幅間的相對關(guān)系;另外,還可看到很強(qiáng)的高頻噪聲和末端假頻.圖12c是利用本文的自適應(yīng)增益限的反Q濾波方法補(bǔ)償后的結(jié)果,最大振幅都小于1.0,且成像時間1.8 s處,最大振幅也接近了1.0;因此,本文自適應(yīng)增益限的反Q濾波方法對有效信號恢復(fù)更好,高頻噪聲也得到了很好的壓制,無論是信噪比還是分辨率都要高于圖12b,這在圖12d中也可得到驗證.

    6 應(yīng)用實例

    對圖13a所示的地面地震剖面,利用穩(wěn)定因子法和本文的自適應(yīng)增益限的反Q濾波方法進(jìn)行振幅補(bǔ)償.為了顯示反Q濾波振幅補(bǔ)償后的細(xì)節(jié),用300 ms的自動增益進(jìn)行數(shù)據(jù)顯示.增益限為20 dB時,穩(wěn)定因子法反Q濾波振幅補(bǔ)償?shù)男Ч鐖D13b:相對于反Q濾波之前的記錄而言,分辨率不僅沒有提高,反而降低;增益限為40 dB時,穩(wěn)定因子法反Q濾波振幅補(bǔ)償?shù)男Ч鐖D13c:相對于反Q濾波之前的記錄而言,記錄的分辨率有所提高.假定地震數(shù)據(jù)動態(tài)范圍為40 dB時,利用本文的自適應(yīng)增益限的反Q濾波得到的結(jié)果如圖13d:記錄分辨率明顯提高,且高于圖13c對應(yīng)的增益限為40 dB時,穩(wěn)定因子法反Q濾波振幅補(bǔ)償?shù)男Ч?

    圖5 不同時刻的振幅譜Fig.5 The normalized amplitude spectra before and after the Earth Q-filter amplitude attenuation at different times

    圖6 增益限20 dB,穩(wěn)定因子法補(bǔ)償后的(a)振幅譜和(b)補(bǔ)償后的結(jié)果Fig.6 (a) Normalized amplitude spectra after the stable factor method with 20 dB gain-limit; (b) Result of the stable factor method with 20 dB gain-limit

    圖7 增益限20 dB,穩(wěn)定因子法補(bǔ)償前后各時刻子波的分辨率:(a) 0.3 s, (b) 0.6 s, (c) 0.9 s, (d) 1.2 s, (e) 1.5 s, (f) 1.8 sFig.7 Wavelet resolution before and after inverse Q-filter at (a) 0.3 s, (b) 0.6 s, (c) 0.9 s, (d) 1.2 s, (e) 1.5 s, (f) 1.8 s

    圖8 增益限40 dB,穩(wěn)定因子法補(bǔ)償后的(a)振幅譜和(b)補(bǔ)償后的結(jié)果Fig.8 (a) Normalized amplitude spectra after the stable factor method with 40 dB gain-limit; (b) Result of the stable factor method with 40 dB gain-limit

    圖9 增益限60 dB,穩(wěn)定因子法補(bǔ)償后的(a)振幅譜和(b)補(bǔ)償后的結(jié)果Fig.9 (a) Normalized amplitude spectra after the stable factor method with 60 dB gain-limit; (b) Result of the stable factor method with 60 dB gain-limit

    圖10 數(shù)據(jù)動態(tài)范圍60 dB時,自適應(yīng)增益限反Q濾波:(a)振幅補(bǔ)償函數(shù); (b)補(bǔ)償后的振幅譜Fig.10 (a) Amplitude compensation function of the self-adaptive method at different times; (b) Normalized amplitude spectrum at different time after the self-adaptive method

    圖11 (a) 自適應(yīng)增益限反Q濾波的結(jié)果,數(shù)據(jù)動態(tài)范圍為60 dB; (b) 圖11a和圖9b之間的差值Fig.11 (a) Result of the self-adaptive method with 60 dB dynamic range; (b) The difference between Figs.11a and 9b

    圖12 (a)加入的隨機(jī)噪聲; (b)穩(wěn)定因子法的結(jié)果,增益限60 dB; (c)自適應(yīng)增益限反Q濾波的結(jié)果, 數(shù)據(jù)動態(tài)范圍60 dB; (d) 圖12b 和 12c之差Fig.12 (a) Added random noise; (b) Result of the stable factor method with 60 dB gain-limit; (c) Result of the self-adaptive method with 60 dB dynamic range; (d) The difference between Fig.12b and 12c

    圖13 (a)疊后時間偏移剖面; (b)穩(wěn)定因子法的結(jié)果,增益限20 dB; (c)穩(wěn)定因子法的結(jié)果,增益限為40 dB; (d)本文新方法的結(jié)果,假設(shè)數(shù)據(jù)動態(tài)范圍40 dBFig.13 (a) The post-stack time migration profile; (b) Result of the stable factor method with 20 dB gain-limit; (c) Result of the stable factor method with 40 dB gain-limit; (d) Result of the self-adaptive method with 40 dB dynamic range

    7 結(jié)論

    本文提出了一種自適應(yīng)增益限的反Q濾波振幅補(bǔ)償方法.在該方法中,振幅補(bǔ)償函數(shù)的增益限時變,且與地震數(shù)據(jù)有效頻帶中的截止頻率有關(guān),因此有效頻帶范圍內(nèi)的地震信號可以完全恢復(fù)成震源子波;另外,該方法中的穩(wěn)定因子是時變增益限的函數(shù),其可以很好地控制反Q濾波的數(shù)值穩(wěn)定性問題并避免吉普斯效應(yīng).本文所提出的自適應(yīng)增益限的反Q濾波方法可用于實際地震數(shù)據(jù)的反Q濾波處理,從而獲得高信噪比和高分辨率的結(jié)果.

    Bano M. 1996.Q-phase compensation of seismic records in the frequency domain.BulletinoftheSeismologicalSocietyofAmerica, 86(4): 1179-1186.

    Bickel S H, Natarajan R R.1985. Plane-waveQdeconvolution.Geophysics, 50(9): 1426-1439.

    Castagna J P, Sun S J, Siegfried R M. 2003. Instantaneous spectral analysis: Detection of low-frequency shadows associated with hydrocarbons.TheLeadingEdge, 22(2): 120-127

    Chen Z B, Chen X H, Li J Y, et al. 2014. A band-limited and robust inverseQfiltering algorithm.OilGeophysicalProspecting(in Chinese), 49(1): 68-75.Futterman W I. 1962. Dispersive body waves.J.Geophys.Res., 69(13): 5279-5291.

    Hale D.1981.Qand adaptive prediction error filters. ∥ Stanford Exploration Project Report, 28: 209-231.

    Hale D.1982. An inverseQ-filter. ∥ Stanford Exploration Project Report, 26: 231-243.

    Hargreaves N D, Calvert A J.1991. InverseQfiltering by Fourier transform.Geophysics, 56(4): 519-527.

    He Z H, Xiong X J, Bian L E. 2008. Numerical simulation of seismic low-frequency shadows and its application.AppliedGeophysics(in Chinese), 5(4): 301-306.

    James D I, Knight R J. 2003. Removal of wavelet dispersion from ground-penetrating radar data.Geophysics, 68(3): 960-970.Kjartansson E. 1979. ConstantQwave propagation and attenuation.J.Geophys.Res., 84(B9): 4737-4748.Liu C, Feng X, Zhang J. 2013. A stable inverseQfiltering using the iterative filtering method.OilGeophysicalProspecting(in

    Chinese), 48(6): 890-895.

    Wang Y H. 2002. A stable and efficient approach of inverseQfiltering.Geophysics, 67(2): 657-663.

    Wang Y H. 2003. Quantifying the effectiveness of stabilized inverseQ-filtering.Geophysics, 68(1): 337-345.

    Wang Y H, Guo J. 2004. Modified Kolsky model for seismic attenuation and dispersion.JournalofGeophysicsandEngineering, 1(3): 187-196.Wang Y H. 2006. InverseQ-filter for seismic resolution enhancement.Geophysics, 71(3): 51-60.Wang H Y. 2008. Inverse-Qfiltered migration.Geophysics, 73(1): S1-S6.

    Widess M A.1982. Quantifying resolving power of seismic systems.Geophysics, 47(8): 1160-1173.

    Yao Z X, Gao X, Li W X. 2003. The forwardQmethod for compensating attenuation and frequency dispersion used in the seismic profile of depth domain.ChineseJ.Geophys. (in Chinese), 46(2): 229-230.Zhang G L, He Z H, Wang X M, et al. 2014. Seismic wave dispersion effects and inverseQ-filter phase compensation.ChineseJ.Geophys. (in Chinese), 57(5): 1655-1663, doi: 10.6038/cjg20140528.

    附中文參考文獻(xiàn)

    陳增保, 陳小宏, 李景葉等. 2014. 一種帶限穩(wěn)定的反Q濾波算法. 石油地球物理勘探, 49(1): 68-75.

    賀振華, 熊曉軍, 邊立恩. 2008. 地震低頻伴影的數(shù)值模擬與應(yīng)用. 應(yīng)用地球物理, 5(4): 301-306.

    劉財, 馮晅, 張瑾. 2013. 穩(wěn)定的迭代法反Q濾波. 石油地球物理勘探, 48(6): 890-895.

    姚振興, 高星, 李維新. 2003. 用于深度域地震剖面衰減與頻散補(bǔ)償?shù)姆碤濾波方法. 地球物理學(xué)報, 46(2): 229-230.

    張固瀾, 賀振華, 王熙明等. 2014. 地震波頻散效應(yīng)與反Q濾波相位補(bǔ)償. 地球物理學(xué)報, 57(5): 1655-1663, doi: 10.6038/cjg20140528.

    (本文編輯 何燕)

    A self-adaptive approach for inverseQ-filtering

    ZHANG Gu-Lan1,2, LIN Jin3, WANG Xi-Ming2, HE Zhen-Hua1, CAO Jun-Xing1,ZHANG Jian-Jun2, HE Xi-Lei1, LIN Kai1, XUE Ya-Juan4

    1StateKeyLaboratoryofOilandGasReservoirGeologyandExploitation,ChengduUniversityofTechnology,Chengdu610059,China2BGP,CNPC,HebeiZhuozhou072751,China3Oil&GasExplorationCompany,ShaanxiYanchangPetroleumCo.,Ltd.,ShaanxiYan′an716000,China4CollegeofCommunicationEngineering,ChengduUniversityofInformationTechnology,Chengdu610225,China

    The EarthQ-filter, with frequency-dependent amplitude attenuation and velocity dispersion, can distort seismic wavelet and reduce the seismic resolution. While the inverseQ-filter attempts to eliminate such effect and hence improve the seismic resolution. But the numerical instability of inverseQ-filter amplitude compensation can reduce the SNR (signal-to-noise ratio) and limit the spatial resolution. Although the gain-limit constrained cut-off frequency method and stable factor method for the inverseQ-filter can control the numerical instability and the SNR, in which the gain-limit is time-invariant and is not associated with the dynamic range of seismic data, it usually suppresses high frequencies at later time and reduces the seismic resolution. In order to solve this problem, we should fully compensate energy of seismic data in the effective frequency band and suppress noise beyond it.This work proposes a self-adaptive approach for inverseQ-filter amplitude compensation based on the cut-off frequency and stable factor methods. The gain-limit in this method is time-variant and self-adaptive to the cut-off frequency of the effective frequency band of seismic data; and the stabilizing factor changes in inverse proportion to the square of the self-adaptive gain-limit.The self-adaptive approach can restore energy of seismic data in the effective frequency band and control the numerical instability, and finally achieve high resolution and high SNR data. Synthetic and real data examples demonstrate that the self-adaptive inverseQ-filter can compensate for energy loss without boosting high frequency noise, and produce desirable seismic images with high quality.In order to achieve desirable resolution and high SNR of seismic data by the self-adaptive approach for inverseQ-filtering, theQvalue and the cut-off frequency of the effective frequency band at any time should be relatively accurate.

    InverseQ-filter;Dynamic range; Self-adaptive; Gain-limit; Stable factor

    10.6038/cjg20150727.

    國家自然科學(xué)基金項目 (41174114,41374111), 國家自然科學(xué)基金青年基金項目 (41404101,41404102)聯(lián)合資助.

    張固瀾,男,1982年生,工程師.2008年在成都理工大學(xué)獲得地球探測與信息技術(shù)專業(yè)碩士學(xué)位;2008年起就職于中國石油集團(tuán)東方地球物理公司,主要從事地震資料處理方面的研究及工作.現(xiàn)在成都理工大學(xué)攻讀地球探測與信息技術(shù)專業(yè)博士學(xué)位. E-mail: zhanggulan@cnpc.com.cn

    10.6038/cjg20150727

    P631

    2014-05-27,2015-06-01收修定稿

    張固瀾,林進(jìn),王熙明等. 2015. 一種自適應(yīng)增益限的反Q濾波.地球物理學(xué)報,58(7):2525-2535,

    Zhang G L, Lin J, Wang X M, et al. 2015. A self-adaptive approach for inverseQ-filtering.ChineseJ.Geophys. (in Chinese),58(7):2525-2535,doi:10.6038/cjg20150727.

    猜你喜歡
    子波頻帶振幅
    一類非線性動力系統(tǒng)的孤立子波解
    Wi-Fi網(wǎng)絡(luò)中5G和2.4G是什么?有何區(qū)別?
    單音及部分頻帶干擾下DSSS系統(tǒng)性能分析
    雙頻帶隔板極化器
    電子測試(2017年15期)2017-12-18 07:18:51
    十大漲跌幅、換手、振幅、資金流向
    十大漲跌幅、換手、振幅、資金流向
    十大漲跌幅、換手、振幅、資金流向
    滬市十大振幅
    地震反演子波選擇策略研究
    調(diào)諧放大器通頻帶的計算及應(yīng)用
    亚洲国产精品久久男人天堂| 午夜福利18| 18禁国产床啪视频网站| 十八禁人妻一区二区| e午夜精品久久久久久久| 精品久久久久久久久久免费视频| h日本视频在线播放| 午夜免费成人在线视频| 亚洲性夜色夜夜综合| 午夜精品久久久久久毛片777| 日韩欧美国产在线观看| 免费无遮挡裸体视频| 午夜免费成人在线视频| 日韩欧美免费精品| 国产中年淑女户外野战色| 真实男女啪啪啪动态图| 国产精品久久久久久久电影 | 最近在线观看免费完整版| 成人18禁在线播放| 又紧又爽又黄一区二区| 免费av毛片视频| 丝袜美腿在线中文| 国产视频内射| 免费av观看视频| 一级毛片高清免费大全| 精品一区二区三区人妻视频| 日本精品一区二区三区蜜桃| 国产精华一区二区三区| 国产高清三级在线| 国产精品一区二区三区四区免费观看 | 亚洲内射少妇av| 日本一二三区视频观看| 天堂网av新在线| 在线a可以看的网站| 久久精品国产自在天天线| 香蕉丝袜av| 精品无人区乱码1区二区| av专区在线播放| 国产一区二区三区视频了| 国产真人三级小视频在线观看| 极品教师在线免费播放| 久久久国产成人精品二区| 欧美日本视频| 18禁黄网站禁片免费观看直播| 噜噜噜噜噜久久久久久91| 亚洲最大成人手机在线| 高潮久久久久久久久久久不卡| 又紧又爽又黄一区二区| 天堂动漫精品| 99精品久久久久人妻精品| 国产成人a区在线观看| 在线观看免费视频日本深夜| 亚洲av免费在线观看| 日韩欧美精品免费久久 | 亚洲av二区三区四区| 中文字幕人成人乱码亚洲影| 亚洲五月天丁香| 亚洲久久久久久中文字幕| 岛国在线观看网站| 国产高清激情床上av| 欧美zozozo另类| 国产v大片淫在线免费观看| 一本综合久久免费| 村上凉子中文字幕在线| 亚洲狠狠婷婷综合久久图片| 少妇的丰满在线观看| 91九色精品人成在线观看| 午夜福利在线在线| 91在线精品国自产拍蜜月 | 国产精品,欧美在线| 国产成人影院久久av| 色尼玛亚洲综合影院| 激情在线观看视频在线高清| 91麻豆精品激情在线观看国产| 欧美黑人巨大hd| 精品99又大又爽又粗少妇毛片 | 波多野结衣巨乳人妻| 国产精品综合久久久久久久免费| 国产精品自产拍在线观看55亚洲| 99久久精品国产亚洲精品| 亚洲国产欧美网| 国产三级中文精品| 日本a在线网址| 蜜桃亚洲精品一区二区三区| 天天一区二区日本电影三级| 国产高潮美女av| 少妇人妻一区二区三区视频| 热99在线观看视频| 精品国产美女av久久久久小说| 免费电影在线观看免费观看| 偷拍熟女少妇极品色| 最后的刺客免费高清国语| 亚洲精品456在线播放app | 日韩欧美精品免费久久 | 色视频www国产| 91麻豆av在线| 精品国内亚洲2022精品成人| 老司机深夜福利视频在线观看| 不卡一级毛片| 国产av在哪里看| 国产精品久久久人人做人人爽| 国产精品精品国产色婷婷| 中亚洲国语对白在线视频| 国产中年淑女户外野战色| 国产av一区在线观看免费| 中文资源天堂在线| 老汉色∧v一级毛片| 99国产极品粉嫩在线观看| 波多野结衣高清无吗| 亚洲精品一区av在线观看| 亚洲美女黄片视频| 五月伊人婷婷丁香| 黄色女人牲交| 久久婷婷人人爽人人干人人爱| 高清日韩中文字幕在线| 亚洲av成人av| 一个人免费在线观看的高清视频| 九九热线精品视视频播放| 黄片大片在线免费观看| 国产精品野战在线观看| 国产av在哪里看| 51午夜福利影视在线观看| 校园春色视频在线观看| 国内精品美女久久久久久| 国产黄a三级三级三级人| 国产成人aa在线观看| 少妇的逼水好多| 无遮挡黄片免费观看| 免费人成视频x8x8入口观看| 国产在线精品亚洲第一网站| 成人av一区二区三区在线看| 天美传媒精品一区二区| 真实男女啪啪啪动态图| 最好的美女福利视频网| 免费看十八禁软件| 免费av观看视频| 18禁黄网站禁片免费观看直播| 亚洲精品影视一区二区三区av| 国产精品一区二区免费欧美| 18禁裸乳无遮挡免费网站照片| 国产v大片淫在线免费观看| 国产精品98久久久久久宅男小说| 麻豆国产97在线/欧美| 黄色片一级片一级黄色片| 亚洲无线观看免费| 18禁国产床啪视频网站| 亚洲一区高清亚洲精品| 在线天堂最新版资源| 99在线人妻在线中文字幕| xxxwww97欧美| 九九在线视频观看精品| 亚洲精品国产精品久久久不卡| 国内精品久久久久久久电影| 变态另类成人亚洲欧美熟女| 国产高清视频在线观看网站| 欧美中文日本在线观看视频| 亚洲精品在线美女| a级毛片a级免费在线| 又爽又黄无遮挡网站| 亚洲色图av天堂| 国产精品免费一区二区三区在线| 亚洲欧美激情综合另类| 黄色成人免费大全| 九九热线精品视视频播放| 久久草成人影院| 国产精品1区2区在线观看.| 97超视频在线观看视频| 亚洲成a人片在线一区二区| 在线十欧美十亚洲十日本专区| 久久99热这里只有精品18| 欧美成狂野欧美在线观看| 亚洲中文日韩欧美视频| 九九久久精品国产亚洲av麻豆| 久久久色成人| 亚洲va日本ⅴa欧美va伊人久久| 搡女人真爽免费视频火全软件 | 又黄又粗又硬又大视频| 日韩欧美一区二区三区在线观看| 国产视频内射| 国内精品美女久久久久久| or卡值多少钱| 欧美日韩亚洲国产一区二区在线观看| 全区人妻精品视频| 噜噜噜噜噜久久久久久91| 欧美另类亚洲清纯唯美| 99视频精品全部免费 在线| АⅤ资源中文在线天堂| 国产视频内射| 亚洲av二区三区四区| 天堂动漫精品| 亚洲人成网站在线播| 夜夜看夜夜爽夜夜摸| a级毛片a级免费在线| 婷婷丁香在线五月| 在线天堂最新版资源| 久久精品综合一区二区三区| 国产精品永久免费网站| 亚洲精品亚洲一区二区| 久久久久国内视频| 国产熟女xx| 亚洲成人久久性| 一二三四社区在线视频社区8| 可以在线观看的亚洲视频| 亚洲无线观看免费| 少妇丰满av| 最近视频中文字幕2019在线8| 国产单亲对白刺激| 特级一级黄色大片| 国产欧美日韩一区二区三| 99精品久久久久人妻精品| 两性午夜刺激爽爽歪歪视频在线观看| 精华霜和精华液先用哪个| 听说在线观看完整版免费高清| 日日干狠狠操夜夜爽| 国产精品嫩草影院av在线观看 | 久久99热这里只有精品18| 久久久国产成人免费| netflix在线观看网站| 国产高清激情床上av| 熟女人妻精品中文字幕| 久久久久国产精品人妻aⅴ院| 国产高清激情床上av| 中国美女看黄片| 久久久精品大字幕| 日本熟妇午夜| 操出白浆在线播放| 美女cb高潮喷水在线观看| 亚洲av不卡在线观看| 脱女人内裤的视频| 女人高潮潮喷娇喘18禁视频| 黄片大片在线免费观看| 国产毛片a区久久久久| 99热只有精品国产| 精品欧美国产一区二区三| 久久国产精品人妻蜜桃| 国产v大片淫在线免费观看| 午夜亚洲福利在线播放| 19禁男女啪啪无遮挡网站| 亚洲精品久久国产高清桃花| www国产在线视频色| 中文资源天堂在线| 在线a可以看的网站| 国产视频一区二区在线看| 伊人久久大香线蕉亚洲五| www日本黄色视频网| 91麻豆精品激情在线观看国产| 日本 欧美在线| 午夜激情福利司机影院| 最新美女视频免费是黄的| 亚洲成av人片在线播放无| 全区人妻精品视频| 国产精品免费一区二区三区在线| 首页视频小说图片口味搜索| 欧美一区二区精品小视频在线| 亚洲aⅴ乱码一区二区在线播放| 啦啦啦韩国在线观看视频| 欧美在线黄色| 成人av在线播放网站| 亚洲欧美日韩东京热| 亚洲18禁久久av| 亚洲av电影不卡..在线观看| 2021天堂中文幕一二区在线观| 嫁个100分男人电影在线观看| 91九色精品人成在线观看| 99热只有精品国产| 久久伊人香网站| 欧美又色又爽又黄视频| 欧美色视频一区免费| 欧美不卡视频在线免费观看| 噜噜噜噜噜久久久久久91| 国产又黄又爽又无遮挡在线| 欧美bdsm另类| 成年免费大片在线观看| 日本免费一区二区三区高清不卡| 国产亚洲精品久久久com| 无限看片的www在线观看| 日本精品一区二区三区蜜桃| 一本精品99久久精品77| 男女做爰动态图高潮gif福利片| 手机成人av网站| 国产国拍精品亚洲av在线观看 | 一级黄色大片毛片| 国产伦人伦偷精品视频| 麻豆一二三区av精品| 日韩 欧美 亚洲 中文字幕| 国产美女午夜福利| 国产一级毛片七仙女欲春2| 欧美zozozo另类| 午夜日韩欧美国产| 国产精品一区二区三区四区久久| 亚洲国产中文字幕在线视频| 亚洲乱码一区二区免费版| 色精品久久人妻99蜜桃| 又黄又粗又硬又大视频| 舔av片在线| 久久久久久久亚洲中文字幕 | 日韩av在线大香蕉| 国产午夜精品论理片| 日韩欧美三级三区| 中文字幕人妻熟人妻熟丝袜美 | 精品久久久久久久久久免费视频| 一进一出好大好爽视频| 欧美国产日韩亚洲一区| 亚洲片人在线观看| 91麻豆av在线| 在线视频色国产色| 免费在线观看亚洲国产| 国产乱人视频| 国产精品美女特级片免费视频播放器| 国产午夜福利久久久久久| 久久亚洲真实| 色老头精品视频在线观看| 久久香蕉国产精品| 免费人成在线观看视频色| 午夜免费观看网址| av黄色大香蕉| 国产精品av视频在线免费观看| 免费在线观看亚洲国产| 日韩精品青青久久久久久| 久久性视频一级片| 九九在线视频观看精品| 欧美成人一区二区免费高清观看| 免费无遮挡裸体视频| 亚洲av一区综合| 免费人成在线观看视频色| 国产私拍福利视频在线观看| 在线播放国产精品三级| 精华霜和精华液先用哪个| 在线看三级毛片| 麻豆久久精品国产亚洲av| 欧美在线一区亚洲| 国产亚洲精品综合一区在线观看| 搡老岳熟女国产| 在线观看免费午夜福利视频| 欧美一区二区精品小视频在线| 免费观看的影片在线观看| 老司机在亚洲福利影院| 午夜亚洲福利在线播放| 无限看片的www在线观看| 国产精品爽爽va在线观看网站| 天堂网av新在线| 亚洲激情在线av| 色综合欧美亚洲国产小说| 嫩草影视91久久| 久久久色成人| 99久久久亚洲精品蜜臀av| 美女cb高潮喷水在线观看| 黄色成人免费大全| 久久天躁狠狠躁夜夜2o2o| 国产精品乱码一区二三区的特点| 日本免费a在线| 亚洲午夜理论影院| 久久久精品欧美日韩精品| 窝窝影院91人妻| 岛国视频午夜一区免费看| 精品日产1卡2卡| 国产极品精品免费视频能看的| 国产精品99久久久久久久久| 日韩成人在线观看一区二区三区| 十八禁人妻一区二区| 熟妇人妻久久中文字幕3abv| 88av欧美| 操出白浆在线播放| 看黄色毛片网站| 美女免费视频网站| 母亲3免费完整高清在线观看| 级片在线观看| 999久久久精品免费观看国产| 久久这里只有精品中国| 久久久色成人| 脱女人内裤的视频| 久久久久亚洲av毛片大全| 叶爱在线成人免费视频播放| 99久久精品国产亚洲精品| av天堂中文字幕网| 丰满的人妻完整版| 国产美女午夜福利| 国产精品野战在线观看| 国产午夜精品久久久久久一区二区三区 | 国产免费男女视频| 香蕉丝袜av| 久久久久久久久中文| 老汉色av国产亚洲站长工具| 欧美一级a爱片免费观看看| 精品乱码久久久久久99久播| 国产成人av激情在线播放| 亚洲av美国av| 99热精品在线国产| 国产三级黄色录像| 少妇的逼水好多| 90打野战视频偷拍视频| 在线播放无遮挡| 乱人视频在线观看| 久久精品亚洲精品国产色婷小说| 69人妻影院| 老司机深夜福利视频在线观看| 免费大片18禁| 伊人久久精品亚洲午夜| 女人十人毛片免费观看3o分钟| 一个人看视频在线观看www免费 | 久久香蕉精品热| 日韩高清综合在线| 色哟哟哟哟哟哟| av女优亚洲男人天堂| 老汉色∧v一级毛片| 国产欧美日韩一区二区精品| 中文字幕av成人在线电影| bbb黄色大片| 大型黄色视频在线免费观看| 免费在线观看成人毛片| 亚洲美女黄片视频| 99久久综合精品五月天人人| 久久精品夜夜夜夜夜久久蜜豆| 欧美一级毛片孕妇| 一a级毛片在线观看| 婷婷丁香在线五月| 亚洲美女视频黄频| xxx96com| 成年女人永久免费观看视频| 9191精品国产免费久久| 国产高清有码在线观看视频| 亚洲精品一卡2卡三卡4卡5卡| 免费看十八禁软件| 两个人的视频大全免费| 在线看三级毛片| 亚洲国产精品999在线| 蜜桃亚洲精品一区二区三区| 久久婷婷人人爽人人干人人爱| 午夜久久久久精精品| 欧美黑人欧美精品刺激| 免费观看精品视频网站| 精品久久久久久久毛片微露脸| 91九色精品人成在线观看| 一夜夜www| 婷婷亚洲欧美| 久久久久九九精品影院| 精品电影一区二区在线| 免费看日本二区| 亚洲美女视频黄频| 免费观看精品视频网站| 最近视频中文字幕2019在线8| 天天一区二区日本电影三级| 高清毛片免费观看视频网站| 国产爱豆传媒在线观看| 性欧美人与动物交配| avwww免费| 国产综合懂色| 综合色av麻豆| 白带黄色成豆腐渣| 欧美色视频一区免费| 午夜免费成人在线视频| 国产成人啪精品午夜网站| 精品乱码久久久久久99久播| 午夜日韩欧美国产| 丝袜美腿在线中文| 三级国产精品欧美在线观看| 国产视频内射| 日韩欧美一区二区三区在线观看| 哪里可以看免费的av片| а√天堂www在线а√下载| 国产淫片久久久久久久久 | 久久精品91无色码中文字幕| 午夜视频国产福利| 日韩欧美三级三区| 欧美3d第一页| 精品国产超薄肉色丝袜足j| 国产精品 国内视频| 国产激情欧美一区二区| 久久这里只有精品中国| 免费在线观看日本一区| 少妇的逼好多水| 国产成+人综合+亚洲专区| 久久久国产精品麻豆| 99久久精品一区二区三区| 淫秽高清视频在线观看| 中文字幕av在线有码专区| 我要搜黄色片| 国产伦在线观看视频一区| 久久精品91无色码中文字幕| 女同久久另类99精品国产91| 午夜激情欧美在线| 女人高潮潮喷娇喘18禁视频| 在线播放国产精品三级| 美女黄网站色视频| 欧美高清成人免费视频www| 中文字幕熟女人妻在线| a级一级毛片免费在线观看| 国产一区二区三区在线臀色熟女| 午夜日韩欧美国产| 99久久无色码亚洲精品果冻| 在线a可以看的网站| 国产成人aa在线观看| 伊人久久大香线蕉亚洲五| 亚洲欧美精品综合久久99| 精品久久久久久,| 欧美大码av| 美女高潮喷水抽搐中文字幕| 国产又黄又爽又无遮挡在线| 舔av片在线| 日日夜夜操网爽| 日本a在线网址| 欧美日韩瑟瑟在线播放| 午夜福利视频1000在线观看| 精品一区二区三区视频在线观看免费| 很黄的视频免费| 精品久久久久久,| 日韩欧美精品免费久久 | 最好的美女福利视频网| e午夜精品久久久久久久| 免费av毛片视频| 国产高清三级在线| 少妇的逼好多水| 内射极品少妇av片p| 91久久精品国产一区二区成人 | 日韩 欧美 亚洲 中文字幕| 亚洲成a人片在线一区二区| 无遮挡黄片免费观看| 国产aⅴ精品一区二区三区波| 亚洲av美国av| 亚洲人成网站高清观看| 偷拍熟女少妇极品色| 亚洲无线观看免费| 桃红色精品国产亚洲av| 国产高潮美女av| 成熟少妇高潮喷水视频| avwww免费| a级一级毛片免费在线观看| 极品教师在线免费播放| 亚洲专区国产一区二区| 欧美乱妇无乱码| 成人三级黄色视频| 日本在线视频免费播放| 一个人免费在线观看的高清视频| 国产成+人综合+亚洲专区| 久久性视频一级片| 欧美一区二区精品小视频在线| 亚洲欧美精品综合久久99| 亚洲国产精品sss在线观看| 国产又黄又爽又无遮挡在线| 亚洲美女视频黄频| 久久久久久久久中文| 男女之事视频高清在线观看| av在线蜜桃| 中文字幕精品亚洲无线码一区| 精品国产三级普通话版| 亚洲av成人精品一区久久| 叶爱在线成人免费视频播放| 成人国产综合亚洲| 久久久国产成人免费| av中文乱码字幕在线| 亚洲天堂国产精品一区在线| 国产黄a三级三级三级人| 白带黄色成豆腐渣| 国产亚洲av嫩草精品影院| 亚洲精品一区av在线观看| 午夜福利成人在线免费观看| 九色成人免费人妻av| 日韩欧美在线乱码| 一夜夜www| 国产激情欧美一区二区| 亚洲av一区综合| 三级国产精品欧美在线观看| 免费在线观看亚洲国产| 在线免费观看的www视频| 国产精品野战在线观看| 一区二区三区免费毛片| 国产免费一级a男人的天堂| 老司机福利观看| 日韩欧美三级三区| 岛国在线观看网站| 欧美xxxx黑人xx丫x性爽| 日本与韩国留学比较| 99久久精品国产亚洲精品| 亚洲,欧美精品.| av专区在线播放| 一边摸一边抽搐一进一小说| 日本与韩国留学比较| 色尼玛亚洲综合影院| 18禁在线播放成人免费| 又黄又爽又免费观看的视频| 黄色片一级片一级黄色片| 一级a爱片免费观看的视频| 看免费av毛片| 日韩欧美免费精品| 一进一出抽搐gif免费好疼| 日本黄色片子视频| 免费高清视频大片| 欧美乱色亚洲激情| 一级作爱视频免费观看| 最好的美女福利视频网| 免费在线观看日本一区| 国产精品精品国产色婷婷| 欧美日韩黄片免| 国产成人aa在线观看| 国产91精品成人一区二区三区| 伊人久久精品亚洲午夜| 嫁个100分男人电影在线观看| 精华霜和精华液先用哪个| 中文字幕av在线有码专区| 国产亚洲av嫩草精品影院| av国产免费在线观看| 成人精品一区二区免费| 高清毛片免费观看视频网站| eeuss影院久久| netflix在线观看网站| 日韩亚洲欧美综合| 黄色女人牲交| 久久精品国产亚洲av涩爱 | 欧美黄色淫秽网站| 日韩欧美一区二区三区在线观看| 熟妇人妻久久中文字幕3abv| 中国美女看黄片| 亚洲精品成人久久久久久| 亚洲不卡免费看| 免费看美女性在线毛片视频|