李 倩,狄?guī)妥?魏建新

(1.中國(guó)地震局地球物理勘探中心,河南鄭州 450002;2.中國(guó)石油大學(xué)(北京)油氣資源與探測(cè)國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京102249;3.中國(guó)石油大學(xué)(北京)CNPC物探重點(diǎn)實(shí)驗(yàn)室,北京102249)

基于稀疏約束反演譜分解的縫洞儲(chǔ)層疊后數(shù)據(jù)去噪應(yīng)用效果分析

李 倩1,2,3,狄?guī)妥?,3,魏建新2,3

(1.中國(guó)地震局地球物理勘探中心,河南鄭州 450002;2.中國(guó)石油大學(xué)(北京)油氣資源與探測(cè)國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京102249;3.中國(guó)石油大學(xué)(北京)CNPC物探重點(diǎn)實(shí)驗(yàn)室,北京102249)

高質(zhì)量的地震剖面對(duì)碳酸鹽巖縫洞儲(chǔ)層的流體識(shí)別、儲(chǔ)量預(yù)測(cè)具有重要作用,因此對(duì)疊后地震數(shù)據(jù)進(jìn)行去噪處理很有必要。為此,分析研究了基于稀疏約束反演譜分解的去噪方法在碳酸鹽巖縫洞儲(chǔ)層疊后地震數(shù)據(jù)中的應(yīng)用效果。首先介紹了該方法的基本原理,然后分別將該方法用于數(shù)值模擬數(shù)據(jù)、物理模擬數(shù)據(jù)以及實(shí)際地震數(shù)據(jù)的處理,并分析去噪效果。在數(shù)值模擬測(cè)試中,將該方法與小波閾值去噪法的處理結(jié)果進(jìn)行了比較,表明該方法具有較好的去噪性能和保幅性。物理模擬測(cè)試以及實(shí)際資料的處理結(jié)果均表明,該方法可以較好的壓制碳酸鹽巖縫洞儲(chǔ)層疊后地震數(shù)據(jù)中的隨機(jī)噪聲和偏移噪聲,使地震剖面中碳酸鹽巖孔洞的弱反射特征得到突顯,提高了地震數(shù)據(jù)處理成果的質(zhì)量,為后續(xù)的定量解釋奠定了良好的基礎(chǔ)。

去噪;稀疏約束反演譜分解;碳酸鹽巖儲(chǔ)層;疊后地震數(shù)據(jù);串珠狀特征;稀疏特性

高質(zhì)量的地震剖面更有利于碳酸鹽巖縫洞儲(chǔ)層的解釋,因此有必要進(jìn)行疊后地震數(shù)據(jù)的去噪處理。常用的地震數(shù)據(jù)去噪方法有中值濾波法、頻域?yàn)V波法、f-k濾波法、多項(xiàng)式擬合法、傅里葉變換法、小波變換法和波動(dòng)方程法等[1-4]。2010年,BONAR等[5]提出了將時(shí)間域數(shù)據(jù)轉(zhuǎn)換到稀疏時(shí)頻域的一種譜分解方法進(jìn)行噪聲壓制。當(dāng)使用該方法來達(dá)到去噪目的時(shí),它就相當(dāng)于基追蹤去噪方法[6]。

目前對(duì)基于稀疏變換的去噪方法研究很多[7-11],其基本原理是地震信號(hào)在某一轉(zhuǎn)換域被認(rèn)為是稀疏的,這就意味著可以用少量高能量的值來代表幾乎整個(gè)地震信號(hào)的能量,而噪聲在該轉(zhuǎn)換域被認(rèn)為不聚焦,因此利用稀疏成分重構(gòu)信號(hào)就可以有效地去除噪聲能量。ELBOTH等[7]等將時(shí)頻域去噪方法應(yīng)用于海洋地震數(shù)據(jù)處理,取得了良好的效果。RODRIGUEZ等[9]提出了一種基于稀疏約束時(shí)頻變換的非相關(guān)噪聲衰減方法,并將之用于合成和實(shí)際的微地震數(shù)據(jù)中,去噪效果證明了該方法的優(yōu)越性。HAN等[10-11]在時(shí)頻域內(nèi)地震數(shù)據(jù)重構(gòu)方法的基礎(chǔ)上,在指定域內(nèi)設(shè)置一個(gè)門檻值,使得地震同相軸比噪聲更集中且具有更強(qiáng)的能量,使用該時(shí)頻譜對(duì)地震信號(hào)進(jìn)行重構(gòu),最終成功降低了隨機(jī)噪聲。

地震信號(hào)在稀疏約束譜域中是稀疏的,而隨機(jī)噪聲與子波庫無關(guān),不會(huì)在稀疏譜中產(chǎn)生大的相關(guān),因此使用稀疏約束反演譜對(duì)信號(hào)進(jìn)行重構(gòu)可以壓制地震信號(hào)中的隨機(jī)噪聲。本文將基于稀疏約束反演譜分解的去噪方法用于碳酸鹽巖縫洞儲(chǔ)層疊后地震數(shù)據(jù)處理,利用數(shù)學(xué)模型、物理模型以及實(shí)際資料進(jìn)行驗(yàn)證,結(jié)果表明該方法能夠有效壓制隨機(jī)噪聲,減少有效信號(hào)的損失。

1 方法原理

1.1 基于稀疏約束反演譜分解的去噪方法

通常,地震褶積模型能夠被描述為一個(gè)地震記錄由地震子波w(t)和反射系數(shù)序列r(t)的褶積再加上隨機(jī)噪聲n(t)組成:

(1)

式中:“*”表示褶積運(yùn)算。地震褶積模型在地震勘探中有重要意義,卻無法表示地震信號(hào)頻率隨時(shí)間的變化。因此,褶積模型被拓展為一個(gè)多子波褶積模型,即地震記錄由一系列不同頻率的地震子波和與之相對(duì)應(yīng)的頻率依賴的偽反射系數(shù)序列褶積之后相加再加上隨機(jī)噪聲得到,此模型被稱之為非平穩(wěn)地震褶積模型。

(2)

式中:k可以認(rèn)為與子波主頻線性相關(guān);J表示參與計(jì)算的反射系數(shù)向量或者子波的總個(gè)數(shù);wk(t)是主頻與角標(biāo)k相對(duì)應(yīng)的地震子波;rk(t)是頻率依賴的反射系數(shù)序列。

根據(jù)數(shù)學(xué)理論,方程(2)可以改寫為:

(3)

即:

(4)

其中,Wk和rk(k∈[1J]分別是中心頻率為fi的復(fù)雷克子波wk(t)形成的褶積矩陣以及與其有關(guān)的反射系數(shù)序列。矩陣D表示褶積矩陣庫,由一系列不同頻率的復(fù)雷克子波組成,也叫子波庫。m則被認(rèn)為是包含了所有偽反射系數(shù)序列的矩陣。地震記錄由一系列中心頻率唯一確定的子波庫在特定時(shí)間被分解為不同的頻率成分,得到頻率依賴的反射系數(shù)矩陣m,這個(gè)反射系數(shù)矩陣就可以被認(rèn)為是地震記錄的時(shí)間—頻率譜[5,11-14]。

如果不考慮隨機(jī)噪聲,頻率依賴的反射系數(shù)矩陣可以直接通過對(duì)D求逆得到,即:

(5)

然而,由于D可能不可逆,用這種方法求解不太可行。另外,由于矩陣m的元素個(gè)數(shù)遠(yuǎn)大于地震記錄的元素個(gè)數(shù),因此該方程式是欠定的。為了解這個(gè)欠定方程,得到稀疏時(shí)頻譜,本文采用了對(duì)反問題的解加額外約束項(xiàng)的方式將欠定方程轉(zhuǎn)化為一個(gè)適定方程[13-17]。本文采用L1范數(shù)作為約束項(xiàng)來產(chǎn)生稀疏解,最終的目標(biāo)函數(shù)為:

(6)

方程右邊第一項(xiàng)是不適定項(xiàng),起到了匹配分解結(jié)果與地震記錄的作用,而右邊第二項(xiàng)是約束項(xiàng),確定了反射系數(shù)最終分布的形狀。λ(λ>0)是正則化算子,主要起調(diào)節(jié)偽反射系數(shù)序列稀疏度以及控制解的穩(wěn)定性的作用。

本文采用快速迭代軟閾值算法[16]求得L1范數(shù)約束反問題的解m之后,再將m代入方程(4)中,即可得去噪后的地震數(shù)據(jù):

(7)

2 方法測(cè)試與應(yīng)用效果分析

為了測(cè)試和分析基于稀疏約束反演譜分解的去噪方法在碳酸鹽巖儲(chǔ)層中的應(yīng)用效果,首先將該方法用于數(shù)值模擬數(shù)據(jù)中,并與小波閾值去噪方法進(jìn)行比較,分析該方法的優(yōu)越性。

2.1 數(shù)值模擬數(shù)據(jù)測(cè)試

設(shè)計(jì)一個(gè)用于數(shù)值模擬的碳酸鹽巖孔洞儲(chǔ)層的地質(zhì)模型,如圖1a所示。圖中的孔洞寬為100m,高度依次是20,30,40,50,60,80,100,120m,圍巖的速度和密度分別是6000m/s和2.73g/cm3,孔洞的速度和密度分別為4000m/s和2.56g/cm3。經(jīng)過波動(dòng)方程正演模擬和數(shù)據(jù)處理,得到了疊后地震數(shù)據(jù)。圖1b 為圖1a黑框部分的疊后地震數(shù)據(jù)。模型各層和孔洞的詳細(xì)參數(shù)見表1。

圖1 碳酸鹽巖孔洞的地質(zhì)模型和疊后地震剖面a 地質(zhì)模型; b 對(duì)應(yīng)圖1a黑框內(nèi)數(shù)據(jù)的疊后地震剖面

層名P波速度/(m·s-1)S波速度/(m·s-1)密度/(g·cm-3)層1450025002.50層2600030002.77層3400023002.50孔洞400020002.56

首先從圖1b的疊后地震數(shù)據(jù)中提取第5個(gè)洞中心的單道數(shù)據(jù)(圖2)進(jìn)行處理,圖3,圖4和圖5分別對(duì)比了不同信噪比情況下稀疏約束反演譜分解去噪與小波閾值去噪結(jié)果。由圖2可見,孔洞處的地震反射特征十分明顯。然后在單道數(shù)據(jù)上分別加上高斯白噪聲使數(shù)據(jù)信噪比為2.0,1.0和0.5,如圖3a,圖4a以及圖5a所示,較強(qiáng)的噪聲對(duì)地震信號(hào)產(chǎn)生了影響,容易對(duì)之后的解釋工作造成不良影響。對(duì)加噪后的單道地震數(shù)據(jù)進(jìn)行稀疏約束反演譜分解,由稀疏約束反演譜分解時(shí)頻圖(圖3b,圖4b和圖5b)可以看出,地震信號(hào)的稀疏分布具有非常好的局部性特征,地震信號(hào)的有效部分具有較強(qiáng)的聚焦性,而隨機(jī)噪聲由于與子波庫不相關(guān),在時(shí)頻圖上很少能夠得到呈現(xiàn)。最后利用公式(7)對(duì)稀疏約束反演譜分解結(jié)果進(jìn)行信號(hào)重構(gòu),得到重構(gòu)后的地震信號(hào)(見圖3c,圖4c和圖5c)。由于該方法只利用少量非零系數(shù)進(jìn)行信號(hào)恢復(fù),因此產(chǎn)生的估計(jì)信號(hào)相對(duì)干凈。為了對(duì)比,利用小波閾值去噪法分別對(duì)圖3a,圖4a和圖5a的含噪信號(hào)去噪,得到結(jié)果如圖3d,圖4d和圖5d所示。采用的小波是db6小波函數(shù),使用啟發(fā)式閾值原則(即heursure準(zhǔn)則)產(chǎn)生閾值,并進(jìn)行軟閾值處理。對(duì)信噪比為2.0的含噪信號(hào)進(jìn)行5層分解;對(duì)信噪比為1.0的含噪信號(hào)進(jìn)行6層分解;對(duì)信噪比為0.5的含噪信號(hào)進(jìn)行4層分解。當(dāng)初始數(shù)據(jù)的信噪比較高時(shí),從圖3和圖4中可以看出基于稀疏約束反演譜分解的去噪方法去噪后孔洞處的反射振幅保持較好,反射振幅得到較好的恢復(fù),僅在部分點(diǎn)有微小的振蕩,而且比較光滑;而小波閾值去噪法去噪后孔洞處的反射振幅保持稍差,振蕩點(diǎn)也較多。當(dāng)初始數(shù)據(jù)的信噪比較低時(shí),從圖5可以看出,雖然基于稀疏約束反演譜分解的去噪方法去噪后振蕩點(diǎn)較多,但是對(duì)于目標(biāo)區(qū)域也就是孔洞處的反射振幅恢復(fù)較好,信號(hào)形變較小;而小波閾值去噪法去噪后雖然振蕩點(diǎn)稍微少一些,但是孔洞處振幅恢復(fù)較差,信號(hào)形變也較大?？偟膩碚f,稀疏約束反演譜分解方法較精確地恢復(fù)了單道信號(hào)的大部分頻率成分。盡管局部信號(hào)的能量稍有降低,孔洞處的反射仍然得到了較好的恢復(fù),證明了該方法對(duì)于降低噪聲的有效性。

圖2 從圖1b的疊后地震數(shù)據(jù)中提取的第5個(gè)洞中心的單道數(shù)據(jù)

圖3 信噪比為2的單道數(shù)據(jù)的稀疏約束反演譜分解去噪與小波閾值去噪結(jié)果對(duì)比a 信噪比為2的單道含噪數(shù)據(jù); b 稀疏約束反演譜分解的時(shí)頻圖; c 從稀疏譜域重構(gòu)的地震信號(hào); d 小波閾值去噪結(jié)果

圖4 信噪比為1的單道數(shù)據(jù)的稀疏約束反演譜分解去噪與小波閾值去噪結(jié)果對(duì)比a 信噪比為1單道含噪數(shù)據(jù); b 稀疏約束反演譜分解的時(shí)頻圖; c 從稀疏譜域重構(gòu)的地震信號(hào); d 小波閾值去噪結(jié)果

信號(hào)的信噪比以及均方根誤差通常被作為信號(hào)去噪性能的評(píng)價(jià)標(biāo)準(zhǔn)。信噪比越高,原始信號(hào)和估計(jì)信號(hào)的均方根誤差越小,說明估計(jì)信號(hào)就越接近原始信號(hào),去噪效果就越好。因此我們給出了基于稀疏約束反演譜分解去噪方法和小波閾值去噪法對(duì)不同含噪信號(hào)去噪后的信號(hào)信噪比(SNR)和均方根誤差(RMSE)(如表2所示)。對(duì)比2種方法的信噪比和均方根誤差,發(fā)現(xiàn)基于稀疏約束反演譜分解去噪方法的信噪比較小波閾值去噪法高,而均方根誤差比小波閾值去噪法小,說明基于稀疏約束反演譜分解去噪方法去噪性能較好。

圖5 信噪比為0.5的單道數(shù)據(jù)的稀疏約束反演譜分解去噪與小波閾值去噪結(jié)果對(duì)比a 信噪比為0.5的單道含噪數(shù)據(jù); b 稀疏約束反演譜分解的時(shí)頻圖; c 從稀疏譜域重構(gòu)的地震信號(hào); d 小波閾值去噪結(jié)果

信號(hào)類型基于稀疏約束反演譜分解去噪小波閾值去噪SNRRMSESNRRMSE信噪比為2.0的原始信號(hào)15.98100.034413.02900.0483信噪比為1.0的原始信號(hào)14.57310.040510.03250.0683信噪比為0.5的原始信號(hào)9.29110.07438.19640.0843

為了測(cè)試該去噪方法在整個(gè)地震剖面上的應(yīng)用效果,在圖1b的地震剖面中提取無噪地震數(shù)據(jù),如圖6 所示,然后加上高斯白噪聲使其信噪比為2.0,1.0和0.5,分別利用基于稀疏約束反演譜分解去噪方法和小波閾值去噪法進(jìn)行處理,結(jié)果如圖7,圖8和圖9所示。小波閾值去噪法選擇的小波、閾值、閾值處理方法以及分解層數(shù)與處理單道數(shù)據(jù)時(shí)選擇的參數(shù)相同。表3中給出了對(duì)整個(gè)地震數(shù)據(jù)分別利用2種去噪方法去噪后結(jié)果的信噪比和均方根誤差。結(jié)合圖7,圖8,圖9和表3進(jìn)行綜合分析,發(fā)現(xiàn)基于稀疏約束反演譜分解去噪方法壓制噪聲的效果整體上比小波閾值去噪法效果好。在初始數(shù)據(jù)的信噪比較高時(shí),基于稀疏約束反演譜分解的去噪方法和小波閾值去噪法的去噪效果都較好,2種方法去噪后數(shù)據(jù)的信噪比和均方根誤差十分接近;隨著初始數(shù)據(jù)的信噪比降低,基于稀疏約束反演譜分解的去噪方法表現(xiàn)相對(duì)穩(wěn)定,去噪后數(shù)據(jù)的信噪比和均方根誤差值變化較小,而小波閾值去噪法去噪后數(shù)據(jù)的信噪比和均方根誤差變化較大,去噪效果逐漸變差。在計(jì)算過程中,發(fā)現(xiàn)基于稀疏約束反演譜分解的去噪方法雖然壓制噪聲效果較好,但存在計(jì)算量大以及在求解稀疏反演譜的過程中正則化算子難以確定的缺點(diǎn);而小波閾值去噪法計(jì)算速度快,但對(duì)低信噪比數(shù)據(jù)的去噪效果稍差,且閾值選擇對(duì)去噪效果至關(guān)重要,要根據(jù)具體情況來選擇合適的閾值。總的來說,不同的去噪方法都有其優(yōu)點(diǎn)以及適用條件,沒有哪一種方法具有普遍適用性。

圖6 從圖1b的疊后地震剖面中提取的無噪地震數(shù)據(jù)

圖7 信噪比為2.0的數(shù)值模擬數(shù)據(jù)的稀疏約束反演譜分解去噪與小波閾值去噪結(jié)果對(duì)比a 信噪比為2.0地震數(shù)據(jù); b 小波閾值去噪結(jié)果; c 基于稀疏約束反演譜分解去噪結(jié)果

圖8 信噪比為1.0的數(shù)值模擬數(shù)據(jù)的稀疏約束反演譜分解去噪與小波閾值去噪結(jié)果對(duì)比a 信噪比為1.0的地震數(shù)據(jù); b 小波閾值去噪結(jié)果; c 基于稀疏約束反演譜分解去噪結(jié)果

圖9 信噪比為0.5的數(shù)值模擬數(shù)據(jù)的稀疏約束反演譜分解去噪與小波閾值去噪結(jié)果對(duì)比a 信噪比為0.5的地震數(shù)據(jù); b 小波閾值去噪結(jié)果; c 基于稀疏約束反演譜分解去噪結(jié)果

信號(hào)類型基于稀疏約束反演譜分解去噪小波閾值去噪SNRRMSESNRRMSE信噪比為2.0的原始信號(hào)13.93120.011913.23210.0129信噪比為1.0的原始信號(hào)13.82730.012010.54070.0175信噪比為0.5的原始信號(hào)13.10890.01306.57240.0277

3 物理模型數(shù)據(jù)和實(shí)際數(shù)據(jù)應(yīng)用效果分析

3.1 物理模型數(shù)據(jù)測(cè)試

將基于稀疏約束反演譜分解的去噪方法應(yīng)用于物理模型數(shù)據(jù),并分析其應(yīng)用效果。該模型是塔里木油田哈拉哈塘地區(qū)碳酸鹽巖孔洞儲(chǔ)層地震物理模型[18],是根據(jù)相似比原理設(shè)計(jì)的與實(shí)際地層比較接近的復(fù)雜孔洞模型,其相似比分別是vm∶v=1∶2,lm∶l=1∶20000,fm∶f=10000∶1(其中,vm是模型材料的速度,v是實(shí)際地層速度,lm是模型的尺度,l是實(shí)際地層的尺度,fm是物理模擬實(shí)驗(yàn)中數(shù)據(jù)采集時(shí)所使用的震源的主頻,f是實(shí)際地震數(shù)據(jù)的主頻)。模型重點(diǎn)考慮了哈拉哈塘地區(qū)主要含油氣層系中奧陶統(tǒng)一間房組至鷹山組一段上部地層,主要具有以下4個(gè)特征:①根據(jù)該地區(qū)的實(shí)際地層情況,目的層上部地層被簡(jiǎn)化為水平層狀地層;②目的層之上3個(gè)相鄰的上覆層整體從南向北逐漸變薄,形成尖滅;③縫洞儲(chǔ)集體主要分布在深度為6300～7000m的地層中;④200多個(gè)不同類型的縫洞儲(chǔ)集體被布置在水平方向上的14個(gè)縫洞區(qū)。

圖10為物理模型Inline線的垂向剖面,圖兩側(cè)標(biāo)示了各層的深度,圖中各層中的數(shù)字則表示該層的層速度和密度,并在圖中標(biāo)出了相應(yīng)的地質(zhì)層位(斷層和河道未在圖中標(biāo)出)。

圖11a是利用復(fù)雜三維地震物理模型得到的地震剖面,圖11b是去噪后的地震剖面。由模型制作引起孔洞反射之間的規(guī)則干擾以及在剖面最右側(cè)的偏移噪聲(如圖11中箭頭所示)在很大程度上被消除,使得孔洞反射更加清晰。

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

將基于稀疏約束反演譜分解去噪方法應(yīng)用于實(shí)際碳酸鹽巖縫洞儲(chǔ)層地震剖面,并分析其應(yīng)用效果。圖12是目的層在3850ms到4300ms間的一個(gè)實(shí)際地震剖面,去噪后(如圖12b所示),在疊加剖面上的隨機(jī)噪聲和剖面最右側(cè)的偏移噪聲(如箭頭所示)得到了很大程度的衰減,目的層段內(nèi)的構(gòu)造變的更加清晰,層與層之間的信息更加豐富。剖面中間位置來自于某些孔洞的弱振幅響應(yīng)和噪聲混雜在一起,去噪之后,“串珠狀”特征清晰可見。圖13為圖12中間部分的放大顯示,在圖13a中看到的是一片雜亂反射,無法判斷是孔洞的弱反射還是干擾,而在圖13b中可以清晰的看到一個(gè)個(gè)來自于孔洞的串珠狀的弱反射,對(duì)后續(xù)的解釋以及儲(chǔ)層預(yù)測(cè)都能提供更好的依據(jù)。

圖10 物理模型Inline線的垂向剖面

圖11 物理模型數(shù)據(jù)的稀疏譜分解去噪a 去噪前的地震剖面; b 去噪后的地震剖面

圖12 實(shí)際數(shù)據(jù)的稀疏譜分解去噪a 去噪前的地震剖面; b 去噪后的地震剖面

圖13 圖12中間部分的放大顯示a 去噪前的地震剖面; b 去噪后的地震剖面

4 結(jié)束語

本文在時(shí)頻變換域中對(duì)地震信號(hào)進(jìn)行稀疏分解,得到有效信號(hào)的稀疏約束反演譜,再通過稀疏約束反演譜對(duì)信號(hào)進(jìn)行重構(gòu),得到衰減了隨機(jī)噪聲的地震記錄。將基于稀疏約束反演譜分解的去噪方法用于碳酸鹽巖縫洞儲(chǔ)層疊后地震數(shù)據(jù)正演,通過數(shù)值模擬數(shù)據(jù)與小波閾值去噪法相比較,證明該方法具有較好的去噪性能。而物理模型數(shù)據(jù)及實(shí)際地震數(shù)據(jù)的去噪結(jié)果表明,該方法能夠有效地去除地震剖面上的隨機(jī)噪聲和偏移噪聲,使得孔洞反射更加清晰,而地震有效信號(hào)的能量損失較小,具有一定的保幅性。

致謝:感謝中國(guó)石油天然氣股份有限公司塔里木油田分公司研究院對(duì)縫洞地震物理模型實(shí)驗(yàn)的支持;感謝中國(guó)石油天然氣股份有限公司石油勘探開發(fā)研究院西北分院對(duì)縫洞地震物理模型數(shù)據(jù)處理的幫助！

[1] 劉婷婷,陳陽康.f-x域經(jīng)驗(yàn)?zāi)Ｊ椒纸馀c多道奇異譜分析相結(jié)合去除隨機(jī)噪聲[J].石油物探,2016,55(1):67-75 LIU T T,CHEN Y K.Random noise attenuation based on EMD and MSSA in f-x domain[J].Geophysical Prospecting for Petroleum,2016,55(1):67-75

[2] 劉洋,王典,劉財(cái),等.局部相關(guān)加權(quán)中值濾波技術(shù)及其在疊后隨機(jī)噪聲衰減中的應(yīng)用[J].地球物理學(xué)報(bào),2011,54(2):358-367 LIU Y,WANG D,LIU C,et al.Weighted median filter based on local correlation and its application to poststack random noise attenuation[J].Chinese Journal of Geophysics,2011,54(2):358-367

[3] 王姣,李振春,王德營(yíng).基于CEEMD的地震數(shù)據(jù)小波閾值去噪方法研究[J].石油物探,2014,53(2):164-172 WANG J,LI Z C,WANG D Y.A method for wavelet threshold denoising of seismic data based on CEEMD[J].Geophysical Prospecting for Petroleum,2014,53(2):164-172

[4] 胡天躍.地震資料疊前去噪技術(shù)的現(xiàn)狀與未來[J].地球物理學(xué)進(jìn)展,2002,17(2):218-223 HU T Y.The current situation and future of seismic data prestack noise attenuation techniques[J].Progress in Geophysics,2002,17(2):218-223

[5] BONAR D C,SACCHI M D.Complex spectral decomposition via inversion strategies[J].Expanded Abstracts of 80thAnnual Internat SEG Mtg,2010:1408-1412

[6] CHEN S B,Donoho D L,SAUNDERS M A.Atomic decomposition by basis pursuit[J].SIAM Journal on Scientific Computing,1998,20(1):33-61

[7] ELBOTH T,PRESTERUD I V,HERMANSEN D.Time-frequency seismic data de-noising[J].Geophysical Prospecting,2010,58(3):441-453

[8] SEJDIC E,DJUROVIC I,JIANG J.Time-frequency feature representation using energy concentration:an overview of recent advance[J].Digital Signal Processing,2009,19(1):153-183

[9] RODRIGUEZ I V,BONAR D,SACCHI M.Microseismic data denoising using a 3C group sparsity constrained time-frequency transform[J].Geophysics,2012,77(2):V21-V29

[10] HAN L,SACCHI M D,HAN L G.Spectral decomposition and de-noising via time-frequency and space-wavenumber reassignment[J].Geophysical Prospecting,2014,62(2):244-257

[11] 韓利.高分辨率全譜分解方法研究[D].長(zhǎng)春:吉林大學(xué),2013 HAN L.Research on the methods of high-resolution full spectrum decomposition[D].Changchun:Jilin University,2013

[12] PORTNIAGUINE O,CASTAGNA J.Inverse spectral decomposition[J].Expanded Abstracts of 74thAnnual Internat SEG Mtg,2004:1786-1789

[13] HAN L,HAN L G,ZHAO L.Inverse spectral decomposition with the SPGL1 algorithm[J].Journal of Geophysics and Engineering,2012,9(4):423-427

[14] 劉炳楊,李勝軍,高建虎,等.最小平方約束反演譜分析方法的應(yīng)用效果分析[J].石油物探,2014,53(5):562-569 LIU B Y,LI S J,GAO J H,et al.The application effects of constrained least-square spectrum analysis method[J].Geophysical Prospecting for Petroleum,2014,53(5):562-569

[15] PURYEAR C L,PORTNIAGUINE O N,COBOS C M,et al.Constrained least-squares spectral analysis:application to seismic data[J].Geophysics,2012,77(5):V143-V167

[16] BECK A,TEBOULLE M.A fast iterative shrinkage-thresholding algorithm for linear inverse problems[J].SIAM Journal on Imaging Sciences,2009,2(1):183-202

[17] GARDNER T J,MAGNASCO M O.Sparse time-frequency representations[J].Proceedings of the National Academy of Sciences of the United States of America,2006,103(16):6094-6099

[18] 李倩,狄?guī)妥?魏建新.碳酸鹽巖儲(chǔ)層孔洞體積的地震物理模擬估算[J].石油地球物理勘探,2014,49(6):1147-1156 LI Q,DI B R,WEI J X.Using seismic physical simulation estimate the pore volume of carbonate reservoir[J].Oil Geophysical Prospecting,2014,49(6):1147-1156

(編輯:朱文杰)

Applicationofdenoisingmethodbasedonsparseconstrainedinversespectraldecompositioninpoststackseismicdataofcave-fracturedreservoirs

LI Qian1,2,3,DI Bangrang2,3,WEI Jianxin2,3

(1.GeophysicalExplorationCenter,ChinaEarthquakeAdministration,Zhengzhou450002,China;2.StateKeyLaboratoryofPetroleumResourceandProspecting,ChinaUniversityofPetroleum,Beijing102249,China;3.CNPCKeyLabofGeophysicsExploration,ChinaUniversityofPetroleum,Beijing102249,China)

High-quality seismic sections are very important for fluid identification and reserve prediction of carbonate reservoirs.Therefore,it is necessary to suppress the noise for poststack seismic data from carbonate reservoirs.The application effect of a denoising method based on sparse constrained inverse spectral decomposition in poststack seismic data from carbonate reservoirs was analyzed.First,the basic principle of the above proposed method was introduced.Then,the proposed method was applied to numerical simulation data,physical modeling data,and field seismic data,and the application effect was analyzed.In the numerical simulation test,the proposed method was compared with the wavelet threshold denoising method,and the results proved that the proposed method has better denoising performance and amplitude preservation than the latter.The application effect of physical modeling data and field seismic data shows that this method can effectively suppress the random noise and migration noise in the stacked seismic section of a carbonate reservoir,highlight the weak reflections from some caves of a carbonate reservoir,improve the quality of the seismic data,and lay a good foundation for the subsequent quantitative interpretation.

denoisng,sparse constrained inverse spectral decomposition,carbonate reservoir,post-stack seismic data,string beads characteristic,sparse property

P631

：A

1000-1441(2017)05-0684-10DOI:10.3969/j.issn.1000-1441.2017.05.009

李倩,狄?guī)妥?魏建新.基于稀疏約束反演譜分解的縫洞儲(chǔ)層疊后數(shù)據(jù)去噪應(yīng)用效果分析[J].石油物探,2017,56(5):693

LI Qian,DI Bangrang,WEI Jianxin,et al.Application of denoising method based on sparse constrained inverse spectral decomposition in poststack seismic data of cavefractured reservoirs

[J].Geophysical Prospecting for Petroleum,2017,56(5):693

2016-03-28;改回日期：2016-10-09。

李倩(1985—),女,工程師,博士,主要從事地震資料處理與解釋方面的研究工作。

國(guó)家科技重大專項(xiàng)(2011ZX05007-006)項(xiàng)目資助。

This research is financially supported by the National Science and Technology Major Project of China (Grant No.2011ZX05007-006).