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

    基于CRP道集的頻寬一致性子波反褶積技術(shù)

    2015-09-28 02:39:19劉力輝
    巖性油氣藏 2015年6期
    關(guān)鍵詞:頻寬反褶積子波

    明 君,劉力輝,丁 燕,李 聰

    (1.中海石油(中國)有限公司天津分公司,天津300452;2.北京諾克斯達(dá)石油科技有限公司,北京100192)

    技術(shù)方法

    基于CRP道集的頻寬一致性子波反褶積技術(shù)

    明君1,劉力輝2,丁燕2,李聰2

    (1.中海石油(中國)有限公司天津分公司,天津300452;2.北京諾克斯達(dá)石油科技有限公司,北京100192)

    將疊后復(fù)賽譜域子波提取技術(shù)用于疊前CRP道集子波提取和疊前提高分辨率處理時(shí),由于道集內(nèi)遠(yuǎn)近偏移距道頻寬不一致,影響了子波的空間一致性,也影響了疊前子波反褶積提高分辨率的效果。為此,對復(fù)賽譜域子波提取技術(shù)進(jìn)行了改進(jìn),使其提取的子波具有空間一致性,并將改進(jìn)后的復(fù)賽譜域子波提取技術(shù)與道集譜平衡技術(shù)相結(jié)合,形成了一種道集內(nèi)頻寬一致性子波反褶積技術(shù)。該技術(shù)增加了子波反褶積希望輸出的人為定義,使其能適應(yīng)低信噪比資料的疊前提高分辨率處理。

    復(fù)賽譜域子波提??;譜平衡;頻寬一致的子波反褶積;疊前提高分辨率

    0 引言

    CRP道集優(yōu)化是疊前反演和疊前屬性提取的基礎(chǔ)[1],因此,基于CRP道集的優(yōu)化處理越來越受到人們的重視。CRP道集優(yōu)化主要包括提高信噪比、提高分辨率及增強(qiáng)保真度等方面的處理[2],其中,提高分辨率處理是關(guān)鍵。AVO屬性分析的基礎(chǔ)是佐布里茲方程,而該方程是基于單一界面的,沒有考慮層間干涉,所以在薄層地區(qū)CRP道集優(yōu)化中的提高分辨率處理對AVO屬性分析至關(guān)重要[3]。此外,CRP道集優(yōu)化中的提高分辨率處理若能拓寬有效頻帶,對疊后進(jìn)一步提高分辨率也具有重要意義,但疊前道集一般信噪比不高,是疊前提高分辨率處理的一個(gè)難點(diǎn)。

    確定性子波反褶積、時(shí)變譜白化和反Q濾波是常規(guī)疊后提高分辨率處理的三大技術(shù)。對于疊前CRP道集上的提高分辨率處理,雖然時(shí)變譜白化和反Q濾波技術(shù)都具有一定的效果,但時(shí)變譜白化技術(shù)要求資料有一定的信噪比,而反Q濾波技術(shù)也受到多種因素的影響,如穩(wěn)定性、噪聲和Q值誤差等。因此,針對低信噪比資料,在實(shí)際應(yīng)用中尋找一種簡便且穩(wěn)定的疊前提高分辨率處理技術(shù),具有很重要的現(xiàn)實(shí)意義。確定性子波反褶積中的同態(tài)反褶積[4]是一種非線性變換,對子波的相位沒有要求,它通過將地震頻譜轉(zhuǎn)換成復(fù)賽譜,可有效地分離出地震子波和反射系數(shù)。該方法抗噪性相對較強(qiáng),應(yīng)用較廣,但將它用于疊前提高分辨率處理時(shí)往往效果不甚理想。這主要是由于其假設(shè)子波不隨時(shí)間和空間發(fā)生變化[5-7],而實(shí)際疊前CRP道集資料中遠(yuǎn)近偏移距道的子波頻帶寬度、振幅能量及相位均會(huì)隨著偏移距的變化而出現(xiàn)較大的改變[5,8],顯示出明顯的空間不一致性,導(dǎo)致使用傳統(tǒng)的同態(tài)反褶積方法達(dá)不到理想的提高分辨率效果。為此,筆者對復(fù)賽譜域混合相位子波提取的同態(tài)反褶積方法進(jìn)行改進(jìn),并將它與CRP道集的譜平衡技術(shù)相結(jié)合,形成一種新的頻寬一致性子波反褶積技術(shù),用于疊前道集的提高分辨率處理。

    1 方法原理

    1.1復(fù)賽譜域子波提取技術(shù)

    忽略噪聲項(xiàng),地震記錄褶積模型在頻率域表示為

    式中:x(ω),w(ω)和ξ(ω)分別為地震記錄、子波和反射系數(shù)的傅氏變換。

    通過對式(1)取對數(shù),可以將其轉(zhuǎn)化為線性系統(tǒng),即

    對式(2)作反傅氏變換,則有

    式中:x(t),w(t)和ξ(ω)分別為地震記錄、子波和反射系數(shù)的時(shí)間序列的復(fù)賽譜序列。

    由于子波和反射系數(shù)序列“平滑度”的差別,子波的復(fù)賽譜一般在原點(diǎn)附近,而反射系數(shù)序列的復(fù)賽譜遠(yuǎn)離原點(diǎn)[9-11]?;谶@種認(rèn)識(shí),在復(fù)賽譜域設(shè)計(jì)低通濾波器就可以實(shí)現(xiàn)子波與反射系數(shù)的分離,達(dá)到子波提取的目的。相同振幅譜而相位譜不同的地震子波對應(yīng)不同的最大或最小相位分量。在復(fù)賽譜域,遵照一定的規(guī)則,可將子波振幅譜分割為2部分,分別對應(yīng)于子波的最大和最小相位分量,分割比例的不同決定了子波的不同相位,由此可以得到最小相位子波、零相位子波和一系列混合相位子波。根據(jù)地震記錄信噪比譜,設(shè)定地震記錄的希望輸出,對地震記錄進(jìn)行確定性子波反褶積,再依據(jù)反褶積后地震記錄方差模的大小和分辨率情況,最終從一系列相同振幅譜而相位譜不同的地震子波中確定理想的混合相位子波。

    1.2頻寬一致性子波反褶積技術(shù)

    頻寬一致性子波反褶積技術(shù)主要是將改進(jìn)后的復(fù)賽譜域子波提取技術(shù)與道集譜平衡技術(shù)進(jìn)行了結(jié)合,其方法原理為:首先進(jìn)行道集內(nèi)子波提取及道集處理,在CRP道集內(nèi)取標(biāo)準(zhǔn)道,采用Gabor變換的時(shí)頻譜對遠(yuǎn)偏移距道進(jìn)行譜平衡處理,并提取子波;其次,在道集間選取多個(gè)提取的子波形成子波集,采用“優(yōu)化平均”法得到一個(gè)空間一致性子波;最后,根據(jù)井-震匹配和預(yù)測目的,人為定義希望輸出子波的振幅譜,并求取反子波,從而獲得分辨率得到了提高的疊前數(shù)據(jù)。

    1.2.1疊前道集譜平衡技術(shù)

    在實(shí)際應(yīng)用中,采用同態(tài)反褶積方法對有些資料進(jìn)行提高分辨率處理的效果好,而有些效果差。仔細(xì)分析后發(fā)現(xiàn),提高分辨率效果差的資料大多都存在道集內(nèi)遠(yuǎn)近偏移距道頻寬不一致的現(xiàn)象,即遠(yuǎn)偏移距道波形“發(fā)胖”[12],這會(huì)造成統(tǒng)計(jì)子波振幅譜時(shí)空間的不一致性。對有些資料作道集內(nèi)遠(yuǎn)近偏移距道的頻譜分析時(shí)會(huì)發(fā)現(xiàn),遠(yuǎn)偏移距道的有效頻帶比近偏移距道的有效頻帶窄,即遠(yuǎn)近偏移距道的頻寬不一致,而且有些差別還比較大。為此可采用譜平衡技術(shù)來提高復(fù)賽譜域子波提取的精度,從而提高反褶積效果。

    所謂譜平衡技術(shù)[2],就是將信噪比高的近、中偏移距道作為標(biāo)準(zhǔn)道,以標(biāo)準(zhǔn)道的時(shí)頻特征為標(biāo)準(zhǔn),將其他偏移距道的時(shí)頻特征[13]向其靠攏,使得遠(yuǎn)近偏移距道的頻寬保持一致。具體步驟為:①將欲補(bǔ)償?shù)牡兰M(jìn)行時(shí)頻展開,并參考標(biāo)準(zhǔn)道的時(shí)頻譜,在時(shí)頻譜上逐點(diǎn)補(bǔ)償遠(yuǎn)偏移距道所丟失的高頻信息,從而得到遠(yuǎn)近偏移距道頻寬一致的道集;②進(jìn)行濾波及道均衡,使各偏移距道的振幅能量保持一致。

    1.2.2復(fù)賽譜統(tǒng)計(jì)“優(yōu)化平均”

    將復(fù)賽譜域子波提取技術(shù)用于疊前道集時(shí),其重點(diǎn)在于如何求取空間一致的地震子波的振幅譜,使得其提取的子波是一種具有代表性的綜合子波(圖1)。對于一個(gè)三維資料,在特定的目的層段,選用一定的線距、道距和道集內(nèi)不同的偏移距道作為分析樣板,即可求得一個(gè)子波集。由于該子波集中不同地震道的信噪比與頻譜特征均不相同,因此得出的子波空間一致性往往不太理想。為得到一個(gè)能代表整體數(shù)據(jù)特征的“典型”子波,筆者采用一種“剔除平均”的方法,根據(jù)井-震匹配和整體反褶積的效果,人為優(yōu)選出一種子波,這樣提取的子波不僅具有全區(qū)的代表性,而且用其所做的反褶積抗噪性較強(qiáng)。

    對比圖1(a)與圖1(b)可以發(fā)現(xiàn),多道統(tǒng)計(jì)平均后得到的子波振幅譜主頻有所提高,頻帶得到了拓寬。從圖1(c)可以看出,多道“優(yōu)選”提取的子波旁瓣得到了壓制,采用這種優(yōu)選的子波進(jìn)行反褶積處理能夠達(dá)到提高分辨率的目的。

    圖1 多道子波振幅譜與單道子波振幅譜對比Fig.1 Comparison between multi-trace wavelet amplitude spectrum and single-trace wavelet amplitude spectrum

    1.2.3人為定義希望輸出

    常規(guī)確定性子波反褶積和時(shí)變譜白化技術(shù)都具有抗噪能力低的局限,但子波反褶積有一個(gè)明顯的優(yōu)勢是能定義希望輸出子波。最常規(guī)的希望輸出子波為Ricker子波,其通過控制Ricker子波主頻來控制希望輸出頻譜的特征,但Ricker子波的頻帶較窄,不能適應(yīng)靈活、復(fù)雜的解釋性處理要求。為此,筆者采用了一種人為定義希望輸出頻譜的方法,即手工定義希望輸出的子波振幅譜(圖2),使其主頻偏向低頻,頻帶得到拓寬,具有俞氏子波的特征,這樣就可以壓制高頻噪聲的抬升,將反褶積后道集的信噪比控制在合理的范圍內(nèi),從而提高低信噪比資料的分辨率和AVO屬性的穩(wěn)定性。

    對比圖2(a)與圖2(b)可看出,通過人為定義希望輸出子波振幅譜進(jìn)行頻寬一致性子波反褶積處理后的振幅譜主頻有所提高,頻帶得到了拓寬,高頻噪聲未抬升,表明頻寬一致性子波反褶積技術(shù)具有一定的抗噪能力。

    圖2 人為定義希望輸出子波振幅譜Fig.2 Artificial definition of the desired output wavelet amplitude spectrum

    2 模型試算

    為說明遠(yuǎn)近偏移距道頻寬不一致對提高分辨率效果的影響,設(shè)計(jì)了如下實(shí)驗(yàn)?zāi)P停旱貙雍穸纫来稳?00m,32.5m,30m,20m,20m,20m,20m,60m;對應(yīng)的縱波速度依次為3100m/s,2500m/s,2800m/s,3000m/s,2700m/s,3300m/s,3000m/s,3400m/s;橫波速度依次為1280m/s,1265m/s,1255.023m/s,1224.745m/s,1653.406m/s,1697.836m/s,1224.745m/s,1749.286m/s;密度依次為2.3 g/cm3,2.12 g/cm3,2.01 g/cm3,2.12 g/cm3,2.07 g/cm3,2.21 g/cm3,2.10 g/cm3,2.20 g/cm3;道間距為40 m,雷克子波的主頻為50 Hz。

    動(dòng)校正后的道集如圖3(a)所示,可看出遠(yuǎn)偏移距道波形“發(fā)胖”。圖3(b)是對圖3(a)模型進(jìn)行復(fù)賽譜域子波反褶積處理后的效果,可看出遠(yuǎn)偏移距道的分辨率未得到有效提高。圖3(c)是對圖3(a)模型進(jìn)行頻寬一致性反褶積處理后的效果,可以看出遠(yuǎn)、近偏移距道的分辨率均得到了有效提高。

    圖3 模型數(shù)據(jù)試算Fig.3 Numerical experiments

    3 實(shí)際資料應(yīng)用

    渤中X油田位于渤海灣盆地渤中凹陷與黃河口凹陷的分界處,該區(qū)沙三段目的層段為湖底扇沉積,發(fā)育的微相類型為辮狀水道,主要巖性為含礫砂巖、砂礫巖和泥巖互層。工區(qū)內(nèi)有20多口井,為一開發(fā)井區(qū),斷塊復(fù)雜,且橫向沉積相變快,儲(chǔ)層厚度為3~5 m,而地震主頻僅為29 Hz。因此,地震資料分辨率低是該區(qū)儲(chǔ)層預(yù)測的主要難題。

    圖4為渤中X油田的實(shí)際道集資料。從疊前道集上看,近偏移距道有多次波,地震資料信噪比較低,主頻僅約20 Hz,原始地震資料的分辨率遠(yuǎn)達(dá)不到薄層反演的要求。此外,由于動(dòng)校正拉伸使得道集資料中遠(yuǎn)近偏移距道頻寬不一致,遠(yuǎn)偏移距道波形“發(fā)胖”,能量降低,給傳統(tǒng)反褶積方法的應(yīng)用造成嚴(yán)重困難。

    圖4 渤中X油田原始CRP道集及相應(yīng)頻譜Fig.4 The seismic profile of CRP gather and its frequency spectrum of Bozhong X oilfield

    對渤中X油田的實(shí)際道集資料進(jìn)行疊前提高分辨率處理,首先采用均值濾波對原始道集(圖4)作去噪優(yōu)化處理,去噪后的道集如圖5(a)所示;然后對去噪后的道集作復(fù)賽譜域子波反褶積處理[圖5(b)]和頻寬一致性子波反褶積處理[圖5(c)]。由于近道信噪比低,所以將譜平衡的模板道選擇為15~20道。對比圖5(b)和圖5(c)可以看出,復(fù)賽譜域子波反褶積處理后遠(yuǎn)偏移距道的分辨率未得到顯著提高,而頻寬一致性子波反褶積處理后遠(yuǎn)、近偏移距道的分辨率均得到了顯著提高。

    圖5 CRP道集提高分辨率效果的對比Fig.5 Comparison of increased resolution of CRP gathers

    與去噪處理后的疊加剖面[圖6(a)]相比,復(fù)賽譜域子波反褶積處理后的疊加剖面[圖6(b)]變化不大,而頻寬一致性子波反褶積處理后的疊加剖面[圖6(c)]目的層主頻由原來的20 Hz增大至40 Hz,頻帶拓寬為53 Hz,可達(dá)到的時(shí)間分辨率為12.2 ms,為后續(xù)將沙三段振幅屬性用于沉積相等研究提供了依據(jù)[14-16]。

    圖6 疊加剖面提高分辨率效果對比Fig.6 Comparison of increased resolution of poststack profiles

    圖7 振幅切片提高分辨率效果對比Fig.7 Comparison of increased resolution of poststack slices

    將提取的沙三段振幅切片與沉積相(圖7)進(jìn)行對比后可看出:原始振幅切片[圖7(a)]信噪比低,沉積模式不清晰;去噪后振幅切片[圖7(b)]信噪比有所提高,連續(xù)性得到增強(qiáng),但沉積模式仍然模糊不清;復(fù)賽譜域子波反褶積振幅切片[圖7(c)]能夠顯現(xiàn)一定的沉積模式;頻寬一致性子波反褶積振幅切片[圖7(d)]與沉積相的模式[圖7(e)]對應(yīng)最好,扇三角洲前緣泥巖發(fā)育區(qū)、扇根和扇中主體輪廓均較清晰,更好地反映了沉積相的平面展布。因此,頻寬一致性子波反褶積技術(shù)對提高低信噪比資料的分辨率具有更好的效果。

    4 結(jié)論

    (1)實(shí)驗(yàn)?zāi)P妥C明,遠(yuǎn)近偏移距道頻寬不一致,采用常規(guī)復(fù)賽譜域子波提取技術(shù)不能有效提高CRP道集的分辨率。為此,對復(fù)賽譜域子波提取技術(shù)進(jìn)行了改進(jìn),并將它與道集譜平衡技術(shù)相結(jié)合,形成了頻寬一致性子波反褶積技術(shù)。

    (2)實(shí)際資料應(yīng)用證明,利用頻寬一致性子波反褶積技術(shù)提取的子波具有空間一致性,能夠有效壓制高頻噪聲的抬升,對提高CRP道集低信噪比資料的分辨率具有較好的效果。

    (References):

    [1]鮑熙杰,趙海波.CRP道集優(yōu)化處理對疊前反演的影響[J].石油地質(zhì)與工程,2013,27(1):39-41. Bao Xijie,Zhao Haibo.CRP gathers optimization effect on prestack inversion[J].Petroleum Geology and Engineering,2013,27(1):39-41.

    [2]劉力輝,楊曉,丁燕,等.基于巖性預(yù)測的CRP道集優(yōu)化處理[J].石油物探,2013,52(5):34-39. Liu Lihui,Yang Xiao,Ding Yan,et al.Advanced process of CRP gather based on the prediction of lithology[J].Geophysical Prospecting for Petroleum,2013,52(5):34-39.

    [3]趙偉,陳小宏,李景葉.薄互層調(diào)諧效應(yīng)對AVO的影響[J].石油物探,2006,45(6):570-573. Zhao Wei,Chen Xiaohong,Li Jingye.The influence of thin interbed tuning effect of AVO[J].Geophysical Prospecting for Petroleum,2006,45(6):570-573.

    [4]王萬里,李國發(fā),桂金詠.混合相位子波有色反褶積[J].巖性油氣藏,2013,25(3):82-86. Wang Wanli,Li Guofa,Gui Jinyong.Colored deconvolution of mixed phase wavelet[J].Lithologic Reservoirs,2013,25(3):82-86.

    [5]Barnes A E.Another look at NMO stretch[J].Geophysics,1992,57(5):749-751.

    [6]李秉富,魏長江.地震資料動(dòng)校正拉伸畸變分析[J].青島海洋大學(xué)學(xué)報(bào),1994(3):39-44. Li Bingfu,Wei Changjiang.An analysis of NMO correction lengthening[J].Journal of Ocean University of Qingdao,1994(3):39-44.

    [7]趙波,丈政軍.除動(dòng)校正拉伸影響的方法[J].石油地球物理勘探,1995,30(3):417-421. Zhao Bo,Zhang Zhengjun.A method for removing NMO correction lengthening[J].Oil Geophysical Prospecting,1995,30(3):417-421.

    [8]黃緒德.反褶積與地震道反演[M].北京:石油工業(yè)出版社,1992:256-278. Huang Xude.De-convolution and seismic inversion[M].Beijing:Petroleum Industry Press,1992:256-278.

    [9]李國發(fā),彭蘇萍,高日勝,等.復(fù)賽譜域提取混合相位子波的方法[J].天然氣工業(yè),2005,25(1):85-87. Li Gufa,Peng Suping,Gao Risheng,et al.Method of mixed phase wavelet extraction in cepstrum domain[J].Gas Industries,2005,25(1):85-87.

    [10]Ulrych T J.Application of homomorphic deconvolution to seismology[J].Geophysics,1971,36(4):650-660.

    [11]Stoffa P L,Buhl P,Bryan G M.The application of homomorphic deconvolution to shallow-water marine seismology-PartⅠ:Models[J].Geophysics,1974,39(4):401-416.

    [12]Xu Y,Chopra S.Improving AVO fidelity by NMO stretching and offset-dependent tuning corrections[J].The Leading Edge,2007,26(12):1548-1551.

    [13]刁瑞,尚新民,芮擁軍,等.時(shí)頻域譜模擬反褶積方法研究[J].巖性油氣藏,2013,25(1):116-121. Diao Rui,Shang Xinmin,Rui Yongjun,et al.Study on spectral modeling deconvolution in time-frequency domain[J].Lithologic Reservoirs,2013,25(1):116-121.

    [14]陳可洋,吳清嶺,李來林,等.松遼盆地三維地震資料連片處理關(guān)鍵技術(shù)及其應(yīng)用效果分析[J].巖性油氣藏,2012,24(2):87-91. Chen Keyang,Wu Qingling,Li Lailin,et al.Key techniques and case study of Multi-survey joint processing of 3-D seismic data in Songliao Basin[J].Lithologic Reservoirs,2012,24(2):87-91.

    [15]張猛剛,洪忠,竇玉壇,等.時(shí)頻分析在蘇里格地區(qū)含氣性檢測中的應(yīng)用[J].巖性油氣藏,2013,25(5):76-85. Zhang Menggang,Hong Zhong,Dou Yutan,et al.Application of time-frequency analysis technology to the gas detection in Sulige area[J].Lithologic Reservoirs,2013,25(5):76-85.

    [16]陶云光,袁剛,陳新安,等.二維地震連片相對保幅處理及儲(chǔ)層預(yù)測技術(shù)在塔西南地區(qū)的應(yīng)用[J].巖性油氣藏,2007,19(1):96-100. Tao Yunguang,Yuan Gang,Chen Xin'an et al.Application of relative preserved amplitude of 2D Merging seismic data processing and reservoir prediction in Southwest Tarim[J].Lithologic Reservoirs,2007,19(1):96-100.

    (本文編輯:王會(huì)玲)

    Wavelet deconvolution with spatial consistency based on CRP gathers

    Ming Jun1,Liu Lihui2,Ding Yan2,Li Cong2
    (1.Tianjin Branch Company,CNOOC,Tianjin 300452,China;2.Beijing Rockstar Oil&Techology Company Ltd.,Beijing 100192,China)

    When poststack wavelet extraction method of cepstrum domain is used for prestack wave extraction of CRP gathers and prestack raising resolution processing,the far and near bandwidth of the gather is not consistent,which influences the spatial consistency of wavelet and the effect of the resolution of prestack wavelet deconvolution.Therefore,this paper improved the wavelet extraction technology by using multi-cepstrum average,made the extraction wavelets have spatial consistency,and combined this improved wavelet extraction technology of cepstrum domain with spectral balance technology,to form a wavelet extraction technology which can keep the bandwidth consistent between the traces. This technology added artificial definition of the output of wavelet deconvolution,which makes it can adapt to prestack extraction frequency processing of the data with low signal to noise ratio.

    wavelet extraction of cepstrum domain;spectral balance technology;wavelet deconvolution keeping bandwidth consistant;prestack raising resolution processing

    P631.4

    A

    1673-8926(2015)06-0097-07

    2015-07-20;

    2015-09-17

    中國海洋石油總公司“十二五”科技重大專項(xiàng)“渤海典型低孔低滲油藏勘探開發(fā)關(guān)鍵技術(shù)研究與實(shí)踐”(編號:CNOOC-KJ-125-ZDXM-07)資助

    明君(1970-),男,博士,高級工程師,主要從事油田地球物理綜合解釋和方法研究。地址:(300452)天津市濱海新區(qū)塘沽閘北路1號609信箱渤海石油研究院。E-mail:mingjun@cnooc.com.cn。

    猜你喜歡
    頻寬反褶積子波
    一類非線性動(dòng)力系統(tǒng)的孤立子波解
    基于AMESim的電液伺服系統(tǒng)頻寬分析
    “心理頻寬”決定你的幸福
    意林(2016年23期)2017-01-04 20:32:58
    NB—IoT揭露更多LTE物聯(lián)網(wǎng)技術(shù)細(xì)節(jié)
    深水波浪破碎特征影響因素的實(shí)驗(yàn)研究
    反褶積試井技術(shù)在計(jì)算低滲儲(chǔ)層滲透率中的應(yīng)用
    地震反演子波選擇策略研究
    保持信噪比的相位分解反褶積方法研究
    基于反褶積與編碼激勵(lì)的長輸管道損傷檢測
    基于倒雙譜的地震子波估計(jì)方法
    亚洲一区二区三区色噜噜| 窝窝影院91人妻| 久99久视频精品免费| 国产成人av教育| 久99久视频精品免费| 国产精品美女特级片免费视频播放器 | 男人舔女人的私密视频| 久久香蕉精品热| 国产精品98久久久久久宅男小说| 日本黄色视频三级网站网址| 久久久久免费精品人妻一区二区| 欧美日本视频| 中文字幕最新亚洲高清| 日韩欧美国产在线观看| 老司机午夜十八禁免费视频| 很黄的视频免费| 校园春色视频在线观看| 一本大道久久a久久精品| 欧美又色又爽又黄视频| 国产黄a三级三级三级人| 给我免费播放毛片高清在线观看| 国产片内射在线| 国产乱人伦免费视频| 在线观看一区二区三区| 色精品久久人妻99蜜桃| 国产av一区在线观看免费| 国产精品国产高清国产av| 午夜福利18| 天天添夜夜摸| 性欧美人与动物交配| netflix在线观看网站| 91字幕亚洲| 国产一区二区在线av高清观看| 久久久久九九精品影院| 精品国产乱子伦一区二区三区| 一区二区三区激情视频| 国产精品,欧美在线| 又大又爽又粗| 一本精品99久久精品77| 亚洲精品中文字幕在线视频| 欧美日韩中文字幕国产精品一区二区三区| 久久久久性生活片| www日本在线高清视频| 无遮挡黄片免费观看| av在线天堂中文字幕| 精品一区二区三区四区五区乱码| 亚洲午夜理论影院| 天天躁夜夜躁狠狠躁躁| 90打野战视频偷拍视频| 在线看三级毛片| 国产成人精品久久二区二区免费| 午夜福利18| 在线观看美女被高潮喷水网站 | 国产亚洲精品久久久久久毛片| www日本黄色视频网| 51午夜福利影视在线观看| 成年女人毛片免费观看观看9| 国产伦一二天堂av在线观看| 男人舔女人下体高潮全视频| 亚洲精品久久国产高清桃花| 怎么达到女性高潮| 日本五十路高清| 亚洲性夜色夜夜综合| 99久久精品国产亚洲精品| 午夜日韩欧美国产| 美女午夜性视频免费| 99国产精品一区二区蜜桃av| 日韩欧美在线二视频| 国语自产精品视频在线第100页| 国产精品av视频在线免费观看| 国产精品久久久久久人妻精品电影| 色综合欧美亚洲国产小说| 不卡av一区二区三区| 国产精品久久久av美女十八| 久久精品91无色码中文字幕| 精品一区二区三区av网在线观看| 男人舔女人的私密视频| 亚洲精品中文字幕一二三四区| 老司机福利观看| 国内精品一区二区在线观看| 日韩欧美三级三区| 99热6这里只有精品| 日韩欧美在线乱码| 日韩精品免费视频一区二区三区| 给我免费播放毛片高清在线观看| 免费观看精品视频网站| 久久欧美精品欧美久久欧美| 狠狠狠狠99中文字幕| 亚洲国产中文字幕在线视频| 久久久久久免费高清国产稀缺| 天堂动漫精品| 国产亚洲欧美98| 又黄又粗又硬又大视频| 国产一区在线观看成人免费| 免费看十八禁软件| 丁香六月欧美| 午夜影院日韩av| 日本a在线网址| 精品久久久久久成人av| 三级国产精品欧美在线观看 | bbb黄色大片| netflix在线观看网站| 久久国产乱子伦精品免费另类| 国产黄色小视频在线观看| 婷婷丁香在线五月| 精品国产超薄肉色丝袜足j| 国产欧美日韩精品亚洲av| 免费在线观看黄色视频的| 亚洲av片天天在线观看| 亚洲精品在线观看二区| 99久久无色码亚洲精品果冻| 观看免费一级毛片| 久久久久亚洲av毛片大全| 久久人妻福利社区极品人妻图片| 色综合欧美亚洲国产小说| av福利片在线观看| 久久久久久国产a免费观看| 免费在线观看成人毛片| 亚洲色图av天堂| 国产精品,欧美在线| tocl精华| 久久精品亚洲精品国产色婷小说| 国产精品一及| 美女 人体艺术 gogo| 国产熟女xx| 欧美在线黄色| 亚洲欧美日韩高清专用| 亚洲av中文字字幕乱码综合| 亚洲成人中文字幕在线播放| 露出奶头的视频| 毛片女人毛片| 中文字幕av在线有码专区| 三级男女做爰猛烈吃奶摸视频| 两个人免费观看高清视频| 午夜免费激情av| 亚洲精品av麻豆狂野| 国产精品98久久久久久宅男小说| 亚洲精品美女久久久久99蜜臀| 在线观看免费午夜福利视频| 国产精品一区二区三区四区免费观看 | 国语自产精品视频在线第100页| 1024香蕉在线观看| av片东京热男人的天堂| 手机成人av网站| 黄片小视频在线播放| 欧美日韩福利视频一区二区| 亚洲国产精品合色在线| 亚洲成人国产一区在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲自偷自拍图片 自拍| 男女做爰动态图高潮gif福利片| 高清毛片免费观看视频网站| 国产午夜精品论理片| 男人舔女人下体高潮全视频| 黄色片一级片一级黄色片| 免费在线观看视频国产中文字幕亚洲| 亚洲av美国av| 老司机午夜福利在线观看视频| 国内精品一区二区在线观看| 久久香蕉国产精品| 欧美色欧美亚洲另类二区| 曰老女人黄片| 欧美性猛交黑人性爽| 99久久综合精品五月天人人| 一区二区三区激情视频| 九九热线精品视视频播放| 色av中文字幕| 成人国产一区最新在线观看| 日韩高清综合在线| 亚洲国产欧美网| 欧美在线黄色| 女警被强在线播放| 一边摸一边做爽爽视频免费| a级毛片a级免费在线| e午夜精品久久久久久久| 国产av麻豆久久久久久久| 在线a可以看的网站| 亚洲电影在线观看av| 国产精品久久久av美女十八| 亚洲av电影在线进入| 狂野欧美白嫩少妇大欣赏| 午夜日韩欧美国产| 99国产极品粉嫩在线观看| 99国产综合亚洲精品| 97超级碰碰碰精品色视频在线观看| 亚洲国产欧美网| 女人高潮潮喷娇喘18禁视频| 亚洲性夜色夜夜综合| 欧美3d第一页| 无限看片的www在线观看| 美女大奶头视频| 90打野战视频偷拍视频| 麻豆国产av国片精品| 91老司机精品| 变态另类成人亚洲欧美熟女| 成人av一区二区三区在线看| 欧美一区二区精品小视频在线| 久久精品国产亚洲av香蕉五月| 午夜亚洲福利在线播放| 亚洲黑人精品在线| 国产区一区二久久| 国内精品久久久久精免费| 一本综合久久免费| 身体一侧抽搐| 亚洲国产看品久久| 淫秽高清视频在线观看| 久久久国产成人免费| 精品一区二区三区av网在线观看| 成人永久免费在线观看视频| 婷婷精品国产亚洲av在线| 亚洲18禁久久av| 国模一区二区三区四区视频 | 一卡2卡三卡四卡精品乱码亚洲| 窝窝影院91人妻| 日本黄大片高清| АⅤ资源中文在线天堂| 成人18禁在线播放| 美女大奶头视频| 好男人在线观看高清免费视频| 亚洲欧美日韩东京热| 亚洲男人天堂网一区| 女人被狂操c到高潮| 国产99白浆流出| 18禁美女被吸乳视频| 50天的宝宝边吃奶边哭怎么回事| 美女免费视频网站| 国产成人av教育| 在线观看66精品国产| 黄色视频不卡| 久9热在线精品视频| svipshipincom国产片| 99在线人妻在线中文字幕| 亚洲九九香蕉| 99热6这里只有精品| www.熟女人妻精品国产| 欧美乱色亚洲激情| 悠悠久久av| 久久精品国产亚洲av香蕉五月| 啪啪无遮挡十八禁网站| 国产亚洲欧美98| 中文字幕高清在线视频| 国产亚洲精品久久久久久毛片| 久久亚洲精品不卡| 夜夜爽天天搞| 亚洲欧洲精品一区二区精品久久久| 亚洲中文字幕一区二区三区有码在线看 | 国产精品九九99| 三级毛片av免费| 一个人观看的视频www高清免费观看 | 伦理电影免费视频| 精品一区二区三区av网在线观看| 久久国产精品影院| 国产高清视频在线观看网站| 1024手机看黄色片| 人妻丰满熟妇av一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 国产91精品成人一区二区三区| 18禁美女被吸乳视频| 国产视频内射| 欧美久久黑人一区二区| 成人高潮视频无遮挡免费网站| 欧美三级亚洲精品| 国产精品,欧美在线| 精品久久久久久久人妻蜜臀av| 国产午夜福利久久久久久| 亚洲av美国av| 真人一进一出gif抽搐免费| 九色国产91popny在线| 小说图片视频综合网站| 久久亚洲精品不卡| 桃红色精品国产亚洲av| 一本精品99久久精品77| ponron亚洲| 国产精品精品国产色婷婷| 久久这里只有精品19| 国产亚洲精品av在线| 少妇人妻一区二区三区视频| 亚洲国产欧美一区二区综合| 午夜久久久久精精品| 欧美三级亚洲精品| 国产午夜精品久久久久久| 99久久国产精品久久久| 看黄色毛片网站| 人人妻,人人澡人人爽秒播| 最近视频中文字幕2019在线8| 成人国语在线视频| 久久精品亚洲精品国产色婷小说| 啦啦啦免费观看视频1| a级毛片a级免费在线| 免费在线观看亚洲国产| 日韩国内少妇激情av| 国产精品久久久人人做人人爽| 别揉我奶头~嗯~啊~动态视频| 久久久久国产一级毛片高清牌| 丰满人妻熟妇乱又伦精品不卡| 精品少妇一区二区三区视频日本电影| 久久99热这里只有精品18| 这个男人来自地球电影免费观看| 免费一级毛片在线播放高清视频| 中文字幕人成人乱码亚洲影| 一本一本综合久久| 国产精品野战在线观看| 亚洲精品一区av在线观看| 老鸭窝网址在线观看| 精品午夜福利视频在线观看一区| 日韩欧美在线二视频| 国产午夜精品论理片| 黄色 视频免费看| 色综合站精品国产| 欧美成人免费av一区二区三区| 精品少妇一区二区三区视频日本电影| 久久久水蜜桃国产精品网| 国产区一区二久久| 国产亚洲av高清不卡| 欧美日韩中文字幕国产精品一区二区三区| 成人国语在线视频| 黑人巨大精品欧美一区二区mp4| 国产一区二区在线av高清观看| 又爽又黄无遮挡网站| 国产视频一区二区在线看| 19禁男女啪啪无遮挡网站| 可以免费在线观看a视频的电影网站| 亚洲午夜理论影院| 国内揄拍国产精品人妻在线| 欧美乱妇无乱码| 中文资源天堂在线| 色噜噜av男人的天堂激情| 后天国语完整版免费观看| 国产午夜精品论理片| 国产蜜桃级精品一区二区三区| 国产熟女午夜一区二区三区| 欧美乱色亚洲激情| 国产精品免费一区二区三区在线| 国产麻豆成人av免费视频| 在线国产一区二区在线| 午夜a级毛片| 别揉我奶头~嗯~啊~动态视频| 搞女人的毛片| 88av欧美| 精品欧美一区二区三区在线| 日本成人三级电影网站| 一二三四社区在线视频社区8| 亚洲成a人片在线一区二区| 变态另类成人亚洲欧美熟女| 久久久久久久午夜电影| www.精华液| 欧美一区二区精品小视频在线| 欧美成人性av电影在线观看| 久久久久久久久免费视频了| 村上凉子中文字幕在线| 国产精品久久视频播放| 日日爽夜夜爽网站| 男人舔女人下体高潮全视频| 亚洲成av人片免费观看| 一级毛片高清免费大全| 五月伊人婷婷丁香| 亚洲欧美日韩高清专用| 亚洲最大成人中文| 变态另类丝袜制服| 国产97色在线日韩免费| 亚洲人成网站在线播放欧美日韩| 天天添夜夜摸| 中文在线观看免费www的网站 | 性色av乱码一区二区三区2| 欧美久久黑人一区二区| 久久国产乱子伦精品免费另类| 男插女下体视频免费在线播放| 91在线观看av| 香蕉av资源在线| e午夜精品久久久久久久| 夜夜看夜夜爽夜夜摸| 香蕉丝袜av| 99久久精品热视频| 日本五十路高清| 亚洲片人在线观看| bbb黄色大片| 欧美日韩黄片免| 18美女黄网站色大片免费观看| www国产在线视频色| 国产精品久久久久久亚洲av鲁大| 国产欧美日韩一区二区精品| 身体一侧抽搐| 久久九九热精品免费| 久久久久久久久免费视频了| 91麻豆精品激情在线观看国产| 久久草成人影院| 亚洲国产精品sss在线观看| 无人区码免费观看不卡| 免费在线观看影片大全网站| 在线观看午夜福利视频| 国产视频一区二区在线看| 免费看十八禁软件| 一区二区三区激情视频| 久久午夜亚洲精品久久| 99久久久亚洲精品蜜臀av| 久久九九热精品免费| 国产一区二区三区视频了| 久久午夜亚洲精品久久| 国产精品综合久久久久久久免费| 午夜老司机福利片| 日本黄大片高清| 久久九九热精品免费| 老汉色∧v一级毛片| 国产欧美日韩一区二区精品| 国产三级中文精品| 久久久久久亚洲精品国产蜜桃av| 18禁黄网站禁片免费观看直播| 免费看日本二区| 国产成人系列免费观看| 欧美国产日韩亚洲一区| 久久久久久久久久黄片| 成人欧美大片| 成在线人永久免费视频| 国产精品亚洲美女久久久| 性色av乱码一区二区三区2| www日本在线高清视频| 两个人视频免费观看高清| 欧美激情久久久久久爽电影| 亚洲国产欧美一区二区综合| 亚洲人成网站在线播放欧美日韩| 给我免费播放毛片高清在线观看| 成人av一区二区三区在线看| x7x7x7水蜜桃| 一级a爱片免费观看的视频| 欧美成人免费av一区二区三区| 日韩免费av在线播放| 俄罗斯特黄特色一大片| 亚洲av电影不卡..在线观看| 手机成人av网站| 国产又色又爽无遮挡免费看| 俺也久久电影网| 午夜久久久久精精品| 国产精品一区二区精品视频观看| 丰满的人妻完整版| 给我免费播放毛片高清在线观看| 精品久久久久久久久久久久久| 色老头精品视频在线观看| 男人舔奶头视频| 深夜精品福利| 禁无遮挡网站| 亚洲国产精品久久男人天堂| 欧美 亚洲 国产 日韩一| 亚洲美女视频黄频| 精品熟女少妇八av免费久了| 亚洲国产看品久久| 免费看美女性在线毛片视频| 又爽又黄无遮挡网站| 国产视频一区二区在线看| 国产一区二区在线观看日韩 | 国产精品电影一区二区三区| 久久人人精品亚洲av| 丝袜人妻中文字幕| 亚洲成av人片在线播放无| 两个人免费观看高清视频| 好男人在线观看高清免费视频| 欧美乱色亚洲激情| 日本精品一区二区三区蜜桃| 中出人妻视频一区二区| 亚洲中文字幕日韩| 免费观看人在逋| 精品久久久久久久久久免费视频| 日韩成人在线观看一区二区三区| 成人手机av| 亚洲va日本ⅴa欧美va伊人久久| 两个人免费观看高清视频| 99riav亚洲国产免费| 人妻丰满熟妇av一区二区三区| 大型黄色视频在线免费观看| 丝袜美腿诱惑在线| 亚洲精品国产一区二区精华液| 精品国产亚洲在线| 欧美性猛交╳xxx乱大交人| 伦理电影免费视频| 女人高潮潮喷娇喘18禁视频| 一本久久中文字幕| 国产伦在线观看视频一区| 免费av毛片视频| 女警被强在线播放| 中文字幕av在线有码专区| 日韩欧美国产一区二区入口| or卡值多少钱| svipshipincom国产片| 精品国产超薄肉色丝袜足j| 国产午夜精品久久久久久| 久久久精品国产亚洲av高清涩受| 免费一级毛片在线播放高清视频| 免费在线观看影片大全网站| 美女黄网站色视频| 久久香蕉激情| 又黄又爽又免费观看的视频| 久久久水蜜桃国产精品网| 美女扒开内裤让男人捅视频| 国产爱豆传媒在线观看 | 亚洲avbb在线观看| 在线观看www视频免费| 免费高清视频大片| 亚洲人成电影免费在线| 免费无遮挡裸体视频| 一个人免费在线观看的高清视频| 法律面前人人平等表现在哪些方面| 一区二区三区高清视频在线| 久久久久久久精品吃奶| 麻豆一二三区av精品| 一区福利在线观看| 成人av在线播放网站| 亚洲自偷自拍图片 自拍| 亚洲精品美女久久久久99蜜臀| 熟女少妇亚洲综合色aaa.| 女人高潮潮喷娇喘18禁视频| 久久香蕉激情| 熟女电影av网| 最近最新中文字幕大全免费视频| 亚洲av成人一区二区三| 长腿黑丝高跟| 一本精品99久久精品77| 国产成人系列免费观看| 91国产中文字幕| 欧美高清成人免费视频www| 麻豆久久精品国产亚洲av| 亚洲,欧美精品.| 一本一本综合久久| 日本熟妇午夜| 露出奶头的视频| 50天的宝宝边吃奶边哭怎么回事| 国产不卡一卡二| 精品无人区乱码1区二区| 日韩中文字幕欧美一区二区| 狠狠狠狠99中文字幕| 亚洲成av人片在线播放无| 久久久久久久午夜电影| 国产99久久九九免费精品| 久久99热这里只有精品18| 免费看日本二区| 免费看a级黄色片| 18禁黄网站禁片午夜丰满| 麻豆久久精品国产亚洲av| 俺也久久电影网| 亚洲欧美精品综合一区二区三区| 亚洲精品中文字幕一二三四区| 无遮挡黄片免费观看| 久久久久久亚洲精品国产蜜桃av| 国产精品免费一区二区三区在线| 国产亚洲精品av在线| 丝袜人妻中文字幕| 国产三级在线视频| 国产人伦9x9x在线观看| 99国产综合亚洲精品| 女同久久另类99精品国产91| 亚洲欧美激情综合另类| 免费电影在线观看免费观看| 国产伦一二天堂av在线观看| 国产精华一区二区三区| 精品久久久久久,| 亚洲人成网站高清观看| 日本a在线网址| 999精品在线视频| 精品午夜福利视频在线观看一区| 天堂动漫精品| 亚洲最大成人中文| 国产av又大| 97人妻精品一区二区三区麻豆| 黑人欧美特级aaaaaa片| 在线观看美女被高潮喷水网站 | 久久中文字幕一级| 午夜两性在线视频| 国产成人欧美在线观看| 亚洲成人精品中文字幕电影| 麻豆成人午夜福利视频| 91老司机精品| 亚洲 欧美 日韩 在线 免费| 亚洲一码二码三码区别大吗| 91在线观看av| 听说在线观看完整版免费高清| 亚洲中文日韩欧美视频| 十八禁网站免费在线| www.熟女人妻精品国产| 亚洲 国产 在线| 欧美乱妇无乱码| 欧美三级亚洲精品| 久久久久国内视频| 欧美日韩精品网址| 国产成人精品久久二区二区91| 好看av亚洲va欧美ⅴa在| 我要搜黄色片| 亚洲一区二区三区不卡视频| 天天添夜夜摸| 午夜日韩欧美国产| 国产精品自产拍在线观看55亚洲| 看黄色毛片网站| 国产精品自产拍在线观看55亚洲| 午夜日韩欧美国产| or卡值多少钱| 成人永久免费在线观看视频| 亚洲狠狠婷婷综合久久图片| 亚洲国产日韩欧美精品在线观看 | 国产黄a三级三级三级人| 欧美日韩精品网址| 夜夜夜夜夜久久久久| 久久久久久久久久黄片| 成人国产综合亚洲| 亚洲人与动物交配视频| 亚洲激情在线av| 久久精品亚洲精品国产色婷小说| 亚洲色图 男人天堂 中文字幕| 久久精品夜夜夜夜夜久久蜜豆 | 十八禁网站免费在线| 777久久人妻少妇嫩草av网站| 老司机午夜十八禁免费视频| 黄片小视频在线播放| www日本在线高清视频| 精品无人区乱码1区二区| 深夜精品福利| 欧美精品亚洲一区二区| e午夜精品久久久久久久| 久久久水蜜桃国产精品网|