基于F-X預(yù)測(cè)濾波的基追蹤算法在雀莫錯(cuò)地區(qū)低信噪比資料處理中的應(yīng)用

王 若,張 帥,劉 曉,祝有海,李志偉

(中國(guó)地質(zhì)調(diào)查局非常規(guī)油氣地質(zhì)重點(diǎn)實(shí)驗(yàn)室,北京100083)

羌塘盆地位于青藏高原北部,面積約為1.8×105km2。盆地內(nèi)有北羌塘坳陷、南羌塘坳陷和中央隆起帶3個(gè)次級(jí)構(gòu)造單元。發(fā)育了巨厚的中生代海相沉積地層,三疊系、侏羅系地層廣泛分布。近些年的研究表明,上三疊統(tǒng)地層為羌塘盆地油氣生成和天然氣水合物成藏最好的目的層系[1-5]。此次研究的主要目的層為上三疊統(tǒng)甲丕拉組(T3j)、波里拉組(T3b)和巴貢組(T3bg)地層。在晚三疊世沉積時(shí)期,研究區(qū)內(nèi)整體以淺水三角洲-碳酸鹽巖臺(tái)地相沉積為主,沉積巖性主要為碎屑巖和碳酸鹽巖。其中,研究區(qū)南部抱布德一帶沉積環(huán)境主要為臺(tái)地相和三角洲相,沉積水體較淺,水動(dòng)力較強(qiáng),物源供應(yīng)較豐富,地層厚度大;研究區(qū)北部沉積環(huán)境以水體相對(duì)較深的陸棚相為主,物源供應(yīng)相對(duì)減少,地層厚度減薄,總體上呈現(xiàn)出基本對(duì)稱(chēng)的地層結(jié)構(gòu)[6-8]。

雀莫錯(cuò)地區(qū)位于羌塘盆地的北羌塘坳陷帶東部的雀莫錯(cuò)與各拉丹東山之間,區(qū)內(nèi)地形起伏大,地形切割強(qiáng)烈,地貌形態(tài)多樣,既有高山區(qū)的冰川侵蝕作用形成的地貌,也有山麓地帶發(fā)育的冰磧堤以及低山區(qū)的谷地、洼地或湖泊、河谷相間地貌。由于地表?xiàng)l件復(fù)雜,研究區(qū)內(nèi)采集得到的地震資料信噪比較低,采用傳統(tǒng)的處理流程得到的地震剖面,其分辨率與橫向連續(xù)性均較差。而反褶積作為一種反演反射系數(shù)序列的方法,可應(yīng)用于雀莫錯(cuò)地區(qū)地震資料的優(yōu)化處理,從而有效解決上述問(wèn)題。

由于地震道的帶限特性,反褶積問(wèn)題通常是病態(tài)的。由此得到的反射系數(shù)形式通常有無(wú)窮多組,需要施加適當(dāng)?shù)募s束從而獲得合適的結(jié)果[9]。在此前的研究中,柯西準(zhǔn)則[10]和L1范數(shù)[11]已被廣泛應(yīng)用于稀疏反射系數(shù)反演。MALLAT等[12]于1993年首先介紹了基于過(guò)完備字典的信號(hào)分解理論,并提出了匹配追蹤算法(MP)。匹配追蹤算法將信號(hào)擴(kuò)展為一個(gè)冗余字典中的小波(或原子)的線性疊加。通過(guò)尋找最符合條件的小波并計(jì)算它與原始地震道之間的殘差,反復(fù)迭代直至殘差的能量低于預(yù)設(shè)的閾值,從而得到反射系數(shù)序列的最優(yōu)解[13-16]?；诙嗤ǖ赖钠ヅ渥粉櫵惴?MCMP)則進(jìn)一步提高了空間的連續(xù)性,并且獲得了高分辨率的時(shí)頻譜[17]。在地震反射系數(shù)反演的方法中,基追蹤算法是一個(gè)相對(duì)較新的算法[18-20]。在使用該方法進(jìn)行反演的過(guò)程中,由奇、偶偶極子組成字典用于分解反射系數(shù)對(duì),通過(guò)使用對(duì)數(shù)障礙法求得系數(shù)向量,并將獲得的系數(shù)乘以字典從而得到反射系數(shù)序列。將該方法用于求解欠定線性方程組時(shí),獲得的解是唯一并且稀疏的,這一特性說(shuō)明基追蹤算法適用于求解稀疏反射系數(shù)反演問(wèn)題。然而,以上介紹的傳統(tǒng)方法是一個(gè)基于單道的過(guò)程,即每一道的反射系數(shù)獨(dú)立進(jìn)行反演,因此獲得的反射系數(shù)剖面通常缺乏橫向連續(xù)性。

針對(duì)雀莫錯(cuò)地區(qū)采集到的地震數(shù)據(jù)信噪比低,傳統(tǒng)的處理流程得到的地震剖面連續(xù)性較差的特點(diǎn),本文將F-X預(yù)測(cè)濾波[21]與傳統(tǒng)的基追蹤算法結(jié)合起來(lái),對(duì)雀莫錯(cuò)地區(qū)的二維地震資料進(jìn)行了有針對(duì)性的處理,將預(yù)測(cè)算子與反褶積流程相結(jié)合,最終獲得了分辨率高、地下結(jié)構(gòu)連續(xù)性較好的地震反射剖面。

1 基于F-X預(yù)測(cè)濾波的基追蹤算法

假設(shè)地層介質(zhì)為一系列具有固定速度和密度的水平層狀介質(zhì),地震記錄道可表示為反射系數(shù)和子波的褶積:

(1)

式中:d(t)為實(shí)際觀測(cè)地震記錄;w(t)為地震子波;r(t)為待求反射系數(shù);n(t)為噪聲。子波一般可通過(guò)高階統(tǒng)計(jì)法求得[22-23]。

為了利用基追蹤算法求解反射系數(shù),可將反射系數(shù)剖面看作一系列地層反射系數(shù)對(duì)的疊加,而反射系數(shù)對(duì)可進(jìn)一步表示為偶反射對(duì)(re)和奇反射對(duì)(ro)的加權(quán)疊加,如圖1所示。

圖1 反射系數(shù)對(duì)可表示為偶反射對(duì)(re)和奇反射對(duì)(ro)的加權(quán)疊加

因此,反射系數(shù)序列可表示為:

式中:Δt為采樣率;nΔt為兩個(gè)反射層之間的時(shí)間厚度;an,l和bn,l表示相應(yīng)的反射系數(shù)對(duì)的權(quán)重系數(shù);l用于調(diào)整該奇或偶反射對(duì)在地震道上的位置,可上移或下移,稱(chēng)為地震道時(shí)序?？紤]到所有可能的地層厚度,n的范圍可以從1變化到表示最大地層厚度的常數(shù)nmax。

進(jìn)一步將(2a)式寫(xiě)為矩陣矢量相乘的形式,可以得到:

(3)

式中:r表示反射系數(shù)的向量形式;m是由an,l和bn,l構(gòu)成的系數(shù)的向量形式;G′是由奇、偶偶極子組成的楔形反射矩陣。因?yàn)樾ㄐ文Ｐ褪怯糜诜纸夥瓷湎禂?shù)的過(guò)完備字典,系數(shù)向量m通常比r更加稀疏,意味著m更適于在稀疏約束條件下求解。

根據(jù)(1)式所描述的褶積模型,(3)式等號(hào)兩側(cè)都可以與子波進(jìn)行褶積,從而得到地震道的矩陣矢量相乘形式:

(4)

式中:d表示地震道的向量形式;G是子波與G′褶積得到的楔形地震響應(yīng)矩陣。從公式(3)可知待求的未知量為an,l及bn,l,與公式(4)中的未知量相同,因此可以利用基追蹤算法基于輸入的地震道信號(hào),已知的子波及楔形反射矩陣求得這些未知量,并且代入公式(3),從而獲得最終的反射系數(shù)。

在反演過(guò)程中,系數(shù)向量m的初始值由下式求得:

(5)

基于初始模型,m可以通過(guò)最小化目標(biāo)函數(shù)((6)式)求得:

(6)

式中:λ為控制數(shù)據(jù)殘差及L1范數(shù)模型約束兩者平衡的權(quán)衡參數(shù)。在基追蹤算法中,(6)式可以轉(zhuǎn)換為一個(gè)線性規(guī)劃問(wèn)題,并可以利用原始對(duì)偶對(duì)數(shù)障礙法加以求解[24-25]。得到矢量m后,反射系數(shù)可以通過(guò)公式(3)求得。

一個(gè)含有平層和傾斜結(jié)構(gòu)的合成地震數(shù)據(jù)(圖2a)被作為輸入來(lái)驗(yàn)證傳統(tǒng)基追蹤方法的有效性。圖2b為采用傳統(tǒng)基追蹤算法得到的反射系數(shù)剖面,可見(jiàn),其結(jié)果非常干凈,并具有連續(xù)的結(jié)構(gòu)。

圖2 合成地震數(shù)據(jù)(a)及采用傳統(tǒng)基追蹤算法得到的反射系數(shù)剖面(b)

以上合成地震數(shù)據(jù)的實(shí)例表明了基追蹤算法的有效性。然而,當(dāng)輸入地震道的信噪比較低時(shí),反演得到的反射系數(shù)剖面橫向連續(xù)性將會(huì)變差。圖3a為添加了隨機(jī)噪聲后的合成地震數(shù)據(jù)(信噪比為2.0),圖3b為利用基追蹤算法得到的反射系數(shù)結(jié)果?？梢园l(fā)現(xiàn),直接使用基追蹤算法求解時(shí),可以去除絕大多數(shù)的噪聲,但一些有效的地震信號(hào)也會(huì)被破壞,從而影響了剖面的橫向連續(xù)性。為了進(jìn)一步提高在低信噪比背景下得到的反射系數(shù)的連續(xù)性,將F-X預(yù)測(cè)濾波與傳統(tǒng)基追蹤算法相結(jié)合用于低信噪比地震數(shù)據(jù)的反射系數(shù)求取。

F-X預(yù)測(cè)濾波方法基于以下假設(shè),即地震信號(hào)圖像可以由頻域中一系列線性相干的反射疊加組成。在此假設(shè)條件下,地震信號(hào)是嚴(yán)格的正弦形式,并在頻域具有可預(yù)測(cè)性。因此,對(duì)某一固定頻率f,某一道地震信號(hào)可以通過(guò)其相鄰的地震道預(yù)測(cè)獲得:

(7)

圖3 添加噪聲的合成地震數(shù)據(jù)(信噪比為2)(a)與采用傳統(tǒng)基追蹤算法得到的反射系數(shù)剖面(b)

地震信號(hào)。

由于輸入的地震信號(hào)和反演得到的反射系數(shù)有相似的空間連續(xù)性,根據(jù)其可預(yù)測(cè)性,反射系數(shù)可表示為:

(8)

式中:Ri和R′分別表示頻域中平滑之后及之前的反射系數(shù)。

基于F-X預(yù)測(cè)濾波的基追蹤算法具體步驟描述如下[26-27]。

步驟1:使用傳統(tǒng)的基追蹤算法計(jì)算反射系數(shù)r′。

步驟2:對(duì)求得的r′應(yīng)用F-X預(yù)測(cè)濾波得到濾波后的反射系數(shù)r,從而有效去除反射系數(shù)剖面上的噪聲,增強(qiáng)剖面的橫向連續(xù)性。需要指出的是,F-X預(yù)測(cè)濾波在頻域中進(jìn)行,而基追蹤算法在時(shí)間域中實(shí)現(xiàn),因此,在實(shí)際操作中,必須先利用傅里葉變換將反射信號(hào)轉(zhuǎn)換到頻域,待預(yù)測(cè)濾波完成后再變換回時(shí)間域。

步驟3:根據(jù)(1)式將濾波后得到的反射系數(shù)r與已知的子波褶積,得到合成地震記錄d。

步驟4:利用(5)式重新計(jì)算系數(shù)向量m的初始值。

步驟5:重復(fù)基追蹤算法流程,根據(jù)步驟4中得到的最新初始值m,迭代最小化目標(biāo)函數(shù)直至獲得滿(mǎn)意的結(jié)果。

根據(jù)地震道的可預(yù)測(cè)性,采用基于F-X預(yù)測(cè)濾波的基追蹤算法可以有效提高最終反射系數(shù)剖面的橫向連續(xù)性。為了驗(yàn)證算法的效果,將添加噪聲的合成地震記錄(圖3a)作為輸入地震信號(hào)來(lái)進(jìn)行效果分析,結(jié)果如圖4所示,可以看出,相比于傳統(tǒng)的基追蹤算法,采用基于F-X預(yù)測(cè)濾波的基追蹤算法能夠獲得更干凈和更連續(xù)的反射系數(shù)剖面。

圖4 采用基于F-X預(yù)測(cè)濾波的基追蹤算法得到的合成地震數(shù)據(jù)(添加了噪聲)反射系數(shù)剖面

2 實(shí)際地震數(shù)據(jù)處理及效果分析

以雀莫錯(cuò)地區(qū)QMC13線為例進(jìn)行實(shí)際地震數(shù)據(jù)處理及效果分析。如圖5所示,由于在雀莫錯(cuò)地區(qū)采集時(shí)使用的是可控震源,導(dǎo)致單炮記錄上初至難以拾取,靜校正問(wèn)題突出。同時(shí),由于全區(qū)不同位置上的單炮記錄的信噪比及目的層品質(zhì)均較低,主要存在面波、聲波、線性干擾、異常振幅和高低頻噪聲等(圖6),原始疊加能量縱向衰減快,橫向變化大,平面上同一測(cè)線炮間能量不均衡,導(dǎo)致產(chǎn)狀無(wú)法識(shí)別。

圖5 雀莫錯(cuò)地區(qū)單炮記錄(一)

圖6 雀莫錯(cuò)地區(qū)單炮記錄(二)

在雀莫錯(cuò)地區(qū)二維地震資料的處理流程中,針對(duì)初至難以拾取的問(wèn)題,選取高程靜校正方案解決靜校正問(wèn)題[28-31],同時(shí)測(cè)試了采用多種保真、保幅的疊前去噪方法壓制面波及線性干擾等干擾波的效果。通過(guò)測(cè)試疊前反褶積參數(shù),有效提高了地震資料的分辨率,拓寬了信號(hào)的頻帶。采用速度分析和剩余靜校正迭代處理技術(shù),解決了速度建模以及短波長(zhǎng)問(wèn)題。在做好噪聲壓制和精細(xì)速度分析以提高成像精度的基礎(chǔ)上,采用多次速度調(diào)整偏移,同時(shí)結(jié)合構(gòu)造及斷層的合理性分析,最終建立精細(xì)偏移速度場(chǎng),確保偏移歸位成像效果。圖7為利用上述處理流程得到的QMC13線偏移剖面,可以看出,剖面的分辨率及橫向連續(xù)性不甚理想。

鑒于此,本文采用基于F-X預(yù)測(cè)濾波的基追蹤算法來(lái)實(shí)現(xiàn)地震反射系數(shù)反演,以期提高地震資料的分辨率及連續(xù)性。圖8a和圖8b分別為使用傳統(tǒng)的基追蹤算法及優(yōu)化后的算法得到的反射系數(shù)剖面。雖然直接觀察整體剖面無(wú)法明確分辨反射系數(shù)剖面的改善效果,但對(duì)比圖7和圖8可以看出,在進(jìn)行反射系數(shù)反演后的剖面上,噪聲被有效地去除,整體構(gòu)造形態(tài)和地層展布規(guī)律更清晰、更易于被刻畫(huà)出來(lái)。

圖7 雀莫錯(cuò)地區(qū)處理后的地震剖面

為了進(jìn)一步展示F-X預(yù)測(cè)濾波基追蹤算法的效果,選取圖8a和圖8b剖面中時(shí)間0.15～0.45s、道號(hào)1400～1500的局部剖面進(jìn)行放大展示(圖9)。對(duì)比紅圈中的波組特征可以看出,在使用改進(jìn)算法后,波組更加連續(xù),一些細(xì)小的層位也得到了恢復(fù)。

圖10為反射系數(shù)及原始地震數(shù)據(jù)的全剖面振幅譜,其中黑色曲線表示原始輸入信號(hào)的振幅譜,紅色曲線表示利用基于F-X預(yù)測(cè)濾波的基追蹤算法求得的反射系數(shù)的振幅譜。從圖10可以看出,反演后,振幅譜的高頻部分被很好地補(bǔ)償,頻帶得到了大幅度拓寬,說(shuō)明采用F-X預(yù)測(cè)濾波基追蹤算法可以求解得到連續(xù)的反射系數(shù)剖面,并且能有效地提高剖面的分辨率。

圖8 雀莫錯(cuò)地區(qū)二維地震數(shù)據(jù)反演得到的反射系數(shù)剖面a 傳統(tǒng)基追蹤算法; b 基于F-X預(yù)測(cè)濾波的基追蹤算法

圖9 雀莫錯(cuò)地區(qū)反演得到的反射系數(shù)局部放大剖面a 傳統(tǒng)基追蹤算法; b 基于F-X預(yù)測(cè)濾波的基追蹤算法

圖10 基于F-X預(yù)測(cè)濾波的基追蹤算法得到的反射系數(shù)振幅譜(紅線)及原始地震數(shù)據(jù)振幅譜(黑線)

利用基于F-X預(yù)測(cè)濾波約束的基追蹤算法得到的反射系數(shù)剖面對(duì)研究區(qū)目的層進(jìn)行構(gòu)造解釋。研究區(qū)整體上呈周低中高的構(gòu)造格局,背斜主要由中央古隆起帶由南向北的擠壓作用形成。該區(qū)現(xiàn)今的構(gòu)造形態(tài)主要以褶皺為主,斷層斷距較小,背斜形態(tài)較完整,上三疊統(tǒng)地層沉積厚度較大且在區(qū)域上分布穩(wěn)定,為天然氣水合物提供了有利的成藏條件。圖11和圖12 分別為采用傳統(tǒng)基追蹤算法和基于F-X預(yù)測(cè)濾波的基追蹤算法得到的反射系數(shù)剖面的構(gòu)造解釋結(jié)果。對(duì)比圖11和圖12可以看出，使用基于F-X預(yù)測(cè)濾波的基追蹤算法處理后的地震剖面分辨率得到提高，波組特征明顯、連續(xù)性增強(qiáng)，能夠更加清晰地刻畫(huà)出研究區(qū)構(gòu)造形態(tài)及地層展布規(guī)律。

3 結(jié)論

傳統(tǒng)基追蹤算法在反演復(fù)雜構(gòu)造區(qū)二維地震資料數(shù)據(jù)時(shí),存在地震剖面的分辨率和連續(xù)性欠佳的問(wèn)題。本文采用基于F-X預(yù)測(cè)濾波的基追蹤算法,在原有基追蹤算法的迭代中加入了F-X預(yù)測(cè)算子,對(duì)每次迭代的結(jié)果進(jìn)行預(yù)測(cè)濾波,在增強(qiáng)結(jié)果的連續(xù)性后將其作為下一次迭代的初值,從而考慮了不同地震道之間的相關(guān)性,在提高分辨率的同時(shí)有效提高了反射系數(shù)剖面的連續(xù)性。

在實(shí)現(xiàn)該算法的過(guò)程中,一方面利用基追蹤算法控制反射系數(shù)的稀疏程度,另一方面利用F-X預(yù)測(cè)濾波器控制反射系數(shù)的橫向連續(xù)性。對(duì)合成記錄進(jìn)行測(cè)試,發(fā)現(xiàn)本方法能夠在保證反射系數(shù)剖面橫向連續(xù)性及分辨率的基礎(chǔ)上,有效地壓制噪聲。

由于羌塘盆地特殊的地理位置及環(huán)境條件約束,地震資料普遍存在品質(zhì)較差(如低信噪比)等問(wèn)題,地震資料的精細(xì)處理成為西藏羌塘地區(qū)油氣及天然氣水合物勘探突破的關(guān)鍵點(diǎn)。雀莫錯(cuò)地區(qū)的實(shí)例應(yīng)用結(jié)果表明,使用該方法處理后的地震反射系數(shù)剖面分辨率得到提高,波組特征明顯、連續(xù)性增強(qiáng),有利于更加清晰地刻畫(huà)出研究區(qū)構(gòu)造形態(tài)及地層展布特征。說(shuō)明采用本方法可以獲得高分辨率及高連續(xù)性的二維地震資料,在類(lèi)似的復(fù)雜環(huán)境及構(gòu)造地區(qū)二維地震資料處理過(guò)程中,具有較好的使用效果與應(yīng)用前景。

參 考 文 獻(xiàn)

[1] 潘佳秋,宋春輝,鮑晶,等.羌塘盆地侏羅系元素地球化學(xué)特征與成鹽層位分析[J].地質(zhì)學(xué)報(bào),2015,89(11):2152-2160

PAN J Q,SONG C H,BAO J,et al.Geochemical characteristics and salt-forming analysis of Jurassic strata in the Qiangtang Basin[J].Acta Geologica Sinica,2015,89(11):2152-2160

[2] 段其發(fā),王健雄,姚華舟,等.青海雀莫錯(cuò)地區(qū)侏羅系雙殼類(lèi)及沉積環(huán)境[J].科學(xué)技術(shù)與工程,2004,4(7):546-550

DUAN Q F,WANG J X,YAO H Z,et al.The Jurassic rocks,bivalves and depositional environments in Quemocuo area,Qinghai Province,west China[J].Science Technology and Engineering,2004,4(7):546-550

[3] 白生海.青海西南部海相侏羅系地層新認(rèn)識(shí)[J].地質(zhì)評(píng)論,1989,35(6):529-536

BAI S H.New recognition of the marine Jurassic strata in southwestern Qinghai[J].Geological Review,1989,35(6):529-536

[4] 陳文西,王劍.晚三疊世—中侏羅世羌塘盆地的形成與演化[J].中國(guó)地質(zhì),2009,36(3):682-693

CHEN W X,WANG J.The formation and evolution of the Qiangtang Basin during the Late Triassic-Middle Jurassic period in northern Tibet[J].Geology in China,2009,36(3):682-693

[5] 李亞林,王成善,伊海生,等.西藏北部新生代大型逆沖推覆構(gòu)造與唐古拉山的隆起[J].地質(zhì)學(xué)報(bào),2006,80(8):1118-1130

LI Y L,WANG C S,YI H S,et al.Cenozoic thrust system and uplifting of the Tanggula mountain,northern Tibet[J].Acta Geologica Sinica,2006,80(8):1118-1130

[6] 張潤(rùn)合,斯春松,陳明,等.西藏羌塘盆地北部拗陷侏羅系層序地層劃分對(duì)比及地質(zhì)意義[J].沉積與特提斯地質(zhì),2015,35(2):1-7

ZHANG R H,SI C S,CHEN M,et al.Sequence stratigraphic division,correlation and geological significance of the Jurassic strata in the northern Qiangtang depression,northern Xizang[J].Sedimentary Geology and Tethyan Geology,2015,35(2):1-7

[7] 張玉修,張開(kāi)均,李勇,等.西藏羌塘盆地東部中-上侏羅統(tǒng)沉積特征及沉積相劃分[J].大地構(gòu)造與成礦學(xué),2007,31(1):52-62

ZHANG Y X,ZHANG K J,LI Y,et al.Characteristics and sedimentary facies of the middle-upper Jurassic clastic rocks in Qiangtang Basin,Tibet[J].Geotectonicaet Metallogenia,2007,31(1):52-62

[8] 方立敏,魯兵,劉池洋,等.羌塘盆地中部隆起的演化及其在油氣勘探中的意義[J].地質(zhì)評(píng)論,2002,48(3):279-283

FANG L M,LU B,LIU C Y,et al.Evolution of central dome in the Qiangtang Basin and its importance in oil-gas exploration[J].Geological Review,2002,48(3):279-283

[9] 馮飛,王征,劉成明,等.基于Shearlet變換稀疏約束地震數(shù)據(jù)重建[J].石油物探,2016,55(5):682-691

FENG F,WANG Z,LIU C M,et al.Seismic data reconstruction based on sparse constraint in the Shearlet domain[J].Geophysical Prospecting for Petroleum,2016,55(5):682-691

[10] SACCHI M D,ULRYCH T J.High-resolution velocity gathers and offset space reconstruction[J].Geophysics,1995,60(4):1169-1177

[11] DEBEYE H W J,VAN RIEL P.Lp-norm deconvolution[J].Geophysical Prospecting,1990,38(4):381-404

[12] MALLAT S,ZHANG Z.Matching pursuit with time-frequency dictionaries[J].IEEE Transactions on Signal Processing,1993,41(12):3397-3415

[13] NGUYEN T.High resolution seismic reflectivity inversion[D].Houston:University of Houston,2008

[14] NGUYEN T,CASTAGNA J.High resolution seismic reflectivity inversion[J].Journal of Seismic Exploration,2010,19(4):303-320

[15] 尹繼堯,鐘磊,張吉輝,等.基于連續(xù)小波變換目標(biāo)處理技術(shù)在儲(chǔ)層預(yù)測(cè)中的應(yīng)用[J].石油物探,2016,55(3):433-440

YIN J Y,ZHONG L,ZHANG J H,et al.Target processing by continuous wavelet transform coefficients applied to reservoir prediction[J].Geophysical Prospecting for Petroleum,2016,55(3):433-440

[16] 楊學(xué)亭,劉財(cái),劉洋,等.基于連續(xù)小波變換的時(shí)頻域地震波能量衰減補(bǔ)償[J].石油物探,2014,53(5):523-529

YANG X T,LIU C,LIU Y,et al.The attenuation compensation of seismic wave energy in time-frequency domain based on the continuous wavelet transform[J].Geophysical Prospecting for Petroleum,2014,53(5):523-529

[17] WANG Y.Multichannel matching pursuit for seismic trace decomposition[J].Geophysics,2010,75(4):61-66

[18] ZHANG R.Seismic reflection inversion by basis pursuit[D].Houston:University of Houston,2010

[19] ZHANG R,CASTAGNA J.Sparse-layer inversion with basis pursuit decomposition[J].Geophysics,2011,76(6):R147-R158

[20] ZHANG R,SEN M,SRINIVASAN S.A prestack basis pursuit seismic inversion[J].Geophysics,2013,78(1):R1-R11

[21] WANG Y.Random noise attenuation using forward-backward linear prediction[J].Journal of Seismic Exploration,1999,8(2):133-142

[22] LU X,WANG Y.Mixed-phase wavelet estimation by iterative linear inversion of high-order statistics[J].Journal of Geophysics and Engineering,2007,4(2):184-193

[23] LU X.Seismic ray impedance inversion[D].London:Imperial College London,2010

[24] CHEN S S,DONOHO D L,SAUNDERS M A.Atomic decomposition by basis pursuit[J].Society for Industrial and Applied Mathematics Review,2001,43(1):129-159

[25] DONOHO D L.For most large underdetermined systems of equations,the minimal L1-norm near-solution approximates the sparsest near-solution[J].Communications on Pure and Applied Mathematics,2004,59(7):907-934

[26] WANG R,WANG Y.F-X prediction filtering combined basis pursuit for seismic reflectivity inversion[J].Expanded Abstracts of 75thEAGE Annual Meeting,2013:Tu-P12-14

[27] WANG R.Seismic reflectivity and impedance inversion in multichannel fashion[D].London:Imperial College London,2015

[28] 湯朝陽(yáng),姚華舟,牛志軍,等.羌塘北部拗陷東段晚三疊世地層沉積特征對(duì)比[J].地質(zhì)與資源,2006,15(2):81-88

TANG C Y,YAO H Z,NIU Z J,et al.Sedimentary features and stratigraphical correlation of late Triassic strata in the east of northern Qiangtang Basin depression[J].Geology and Resources,2006,15(2):81-88

[29] 劉池洋,楊興科,任戰(zhàn)利,等.羌塘盆地雀莫錯(cuò)沉降-堆積中心成因:熱力衰減塌陷沉降[J].石油與天然氣地質(zhì),2005,26(2):147-154

LIU C Y,YANG X K,REN Z L,et al.Genesis of subsidence-accumulation centre in Quemocuo area of Qiangtang Basin:caving-in subsidence with thermal attenuation[J].Oil & Gas Geology,2005,26(2):147-154

[30] 譚富文,王劍,王小龍,等.羌塘盆地雁石坪地區(qū)中-晚侏羅世碳、氧同位素特征與沉積環(huán)境分析[J].地球?qū)W報(bào),2004,25(2):119-126

TAN F W,WANG J,WANG X L,et al.Analysis of carbon and oxygen isotope composition and sedimentary environment of the Yanshipingarea of the Qiangtang Basin in middle-late Jurassic[J].Acta Geoscientica Sinica,2004,25(2):119-126

[31] 吳波,徐天吉,唐建明,等.三種反射剩余靜校正方法對(duì)比研究與應(yīng)用[J].石油物探,2012,51(2):172-177

WU B,XU T J,TANG J M,et al.The comparison and application of three reflection residual static correction methods[J].Geophysical Prospecting for Petroleum,2012,51(2):172-177