許自龍,孟繁舉，唐 勇，宋 林，夏洪瑞

(1.中國石油化工股份有限公司石油物探技術(shù)研究院，江蘇南京211103;2.大慶鉆探工程公司地球物理勘探一公司研究院，黑龍江大慶163357;3.北京金陽普泰石油技術(shù)股份有限公司，北京100107;4.中國石油化工集團(tuán)公司石油工程地球物理有限公司江漢分公司物探研究中心，湖北潛江433100)

疊前反演能夠得到比疊后反演更準(zhǔn)確的彈性參數(shù)，可提高儲(chǔ)層預(yù)測(cè)尤其是流體預(yù)測(cè)精度，因而逐步成為儲(chǔ)層預(yù)測(cè)的關(guān)鍵技術(shù)。在生產(chǎn)實(shí)踐中，疊前反演已逐步成為地震資料解釋性處理的常規(guī)環(huán)節(jié)。疊前反演主要分為疊前AVO反演、疊前彈性反演和疊前地震全波形反演。目前，反演技術(shù)的發(fā)展仍以修正和完善各種反演算法為主，眾多專家學(xué)者提出了許多卓有成效的反演算法并在實(shí)踐中取得了較好的效果[1-7]。當(dāng)然，疊前反演的準(zhǔn)確性不僅取決于反演算法，也取決于參與反演數(shù)據(jù)的品質(zhì)。對(duì)于疊前反演所用地震數(shù)據(jù)品質(zhì)的改善問題，已有一些地球物理工作者進(jìn)行了有益的嘗試并取得了成功的經(jīng)驗(yàn)[8-12]，但實(shí)際生產(chǎn)中目前大都仍然采用常規(guī)處理的流程和技術(shù)。

疊前反演是在CMP道集的單道記錄上進(jìn)行，這種處理失去了疊加技術(shù)這個(gè)最強(qiáng)有力去噪技術(shù)的支撐，因而對(duì)參與疊前反演的數(shù)據(jù)品質(zhì)有更高的要求。在常規(guī)處理流程中，一切以求出可供解釋的偏移結(jié)果為主要目的，較少涉及為地震儲(chǔ)層描述而進(jìn)行保持反射波動(dòng)力學(xué)特征的處理。另外，由于預(yù)處理之后的地震數(shù)據(jù)中存有的剩余時(shí)差和殘余噪聲大都可由偏移來消除，故在常規(guī)處理流程中，對(duì)由采集因素產(chǎn)生的誤差及由處理流程、處理技術(shù)及不當(dāng)參數(shù)而產(chǎn)生的處理誤差有相當(dāng)?shù)娜菰S度。而對(duì)上述誤差，疊前反演算法對(duì)疊前反演數(shù)據(jù)不可能有偏移算法對(duì)偏移數(shù)據(jù)那樣的容許度。這是因?yàn)榀B前反演處理應(yīng)用的數(shù)據(jù)是單道記錄，這些單道記錄沒有經(jīng)過保真處理，也沒有經(jīng)過疊加過程去除上述各種誤差，因而沒有滿足疊前反演嚴(yán)格的應(yīng)用條件，其將給疊前反演結(jié)果帶來非常大的陷阱。

由于偏移速度不可能準(zhǔn)確提供，偏移公式自身也僅僅是地震波傳播規(guī)律的近似，因此經(jīng)偏移后CRP道集中的有效信號(hào)不可能嚴(yán)格排齊。這樣的數(shù)據(jù)不能完全滿足AVO反演的應(yīng)用前提(道集內(nèi)的反射同相軸為水平)，由此而得到的反演結(jié)果會(huì)帶有一些錯(cuò)誤信息。

噪聲是與有效信號(hào)如影隨行的干擾信息，地震資料處理中的所有處理技術(shù)都涉及如何盡可能地減少噪聲對(duì)有效信號(hào)的影響。遺憾的是，目前仍沒有一種去噪技術(shù)能實(shí)現(xiàn)真正的信噪分離。更為嚴(yán)重地是，目前噪聲消除幾乎全部采用多道識(shí)別多道消除的思路，經(jīng)這種去噪思路去噪處理后的記錄不可避免地存有擬合噪聲。即只要應(yīng)用了多道消除思想，其一定利用相鄰道的信息，處理記錄的動(dòng)力學(xué)特征必然會(huì)受到破壞。這樣的記錄應(yīng)用于疊前反演，結(jié)果的真實(shí)性難以保證。

在地震數(shù)據(jù)采集和預(yù)處理過程中，CMP道集中的記錄存有的振幅、相位誤差，這些誤差同樣影響疊前反演結(jié)果。

綜上所述，常規(guī)處理后的地震數(shù)據(jù)不能直接應(yīng)用于疊前反演，而需應(yīng)用新的技術(shù)對(duì)疊前反演數(shù)據(jù)進(jìn)行疊前反演預(yù)處理，以滿足疊前反演的需要。疊前反演數(shù)據(jù)優(yōu)化處理的內(nèi)容主要包括吸收補(bǔ)償，頻率、相位校正，剩余時(shí)差校正和去噪處理等。我們僅針對(duì)剩余時(shí)差校正及隨機(jī)噪聲消除(高保真去噪)問題展開討論。

1 疊前反演數(shù)據(jù)的優(yōu)化處理

1.1 基于常規(guī)多項(xiàng)式擬合的剩余時(shí)差校正技術(shù)

目前進(jìn)行剩余時(shí)差校正的辦法通常為速度調(diào)整和靜校正方法，這兩種方法均在一定程度上減少剩余時(shí)差，使信號(hào)同相軸相對(duì)排齊，但也都存有一定的缺陷。由于現(xiàn)行的動(dòng)校/偏移算法均無法完全準(zhǔn)確描述地震波的運(yùn)動(dòng)規(guī)律，其速度調(diào)整只是相對(duì)概念，不可能將同一道集內(nèi)的有效信息完全排齊。靜校方法采用的是單道整體移動(dòng)，而疊前反演數(shù)據(jù)存在的剩余時(shí)差是由動(dòng)校、靜校及其它因素而產(chǎn)生，每道中不同時(shí)刻的有效信息存有時(shí)變性的剩余時(shí)差，用靜校方法實(shí)施之后會(huì)產(chǎn)生顧此失彼的應(yīng)用效果。

實(shí)現(xiàn)剩余時(shí)差校正的關(guān)鍵因素有二：①每道所接收信息的正確位置;②如何將有效信息校正到正確位置。我們知道，單道信息自身無法確定信息位置是否正確，只能靠多道來識(shí)別。由疊加原理知，疊加結(jié)果可相對(duì)準(zhǔn)確地消除剩余時(shí)差。因而，對(duì)疊前反演CRP道集中每道記錄而言，將某一CRP道集中每道記錄分別與對(duì)應(yīng)的疊加道分時(shí)窗互相關(guān)，可求出各道記錄中每個(gè)時(shí)窗的剩余時(shí)差，然后采用濾波插值，即可消除各道記錄中每一時(shí)刻的剩余時(shí)差[13]。

上述技術(shù)在建立模型道時(shí)采用線性擬合技術(shù)，而地震有效信息并非呈線性，由此建立的模型道并不能準(zhǔn)確描述地震道信息的正確位置，因而需采用更為精確的算法來建立模型道。常規(guī)多項(xiàng)式擬合技術(shù)[14]是用曲線來擬合地震有效信息以達(dá)到消除噪聲的技術(shù)，借用其用曲線擬合地震有效信息的思想來建立模型道，便可克服文獻(xiàn)[13]中用直線擬合建立模型道存在的不足。另外，文獻(xiàn)[13]中采用分時(shí)窗求時(shí)差，并用拉伸/壓縮(濾波插值技術(shù))消除時(shí)差，這種方法僅考慮了時(shí)窗內(nèi)的整體效果，故不能使信號(hào)同相軸真正排齊;同時(shí)，拉伸/壓縮技術(shù)將改變信號(hào)的波形特征。為克服這些缺陷，我們采用分波形求時(shí)差，然后采用靜態(tài)時(shí)移方法消除時(shí)差。

以下以二維資料為例(三維資料則為二維擬合的擴(kuò)展)進(jìn)行分析。

1.1.1 常規(guī)多項(xiàng)式擬合建立模型道

為求取準(zhǔn)確的疊加道以校正剩余時(shí)差，首先將疊前反演CRP道集疊加。因有些疊加道也存有誤差，應(yīng)對(duì)疊加道進(jìn)行校正，才能得到更為準(zhǔn)確的模型道。

根據(jù)文獻(xiàn)[14]，對(duì)k疊加道上的i時(shí)刻tk,i(0)，應(yīng)用常規(guī)多項(xiàng)式擬合公式，通過計(jì)算相鄰疊加道上參與擬合的時(shí)刻tk,i(x)，可擬合出經(jīng)過校正的疊加模型道。這種方法既消除了疊加道自身的誤差，又增強(qiáng)了相鄰道之間的連續(xù)性。常規(guī)多項(xiàng)式擬合公式為

(1)

式中：tk,i(x)為道相對(duì)號(hào)為x關(guān)于第k道時(shí)刻i的擬合時(shí)間值(以樣點(diǎn)數(shù)表示);x為相鄰道關(guān)于第k道的相對(duì)位置,x=-X,-X+1,…,-1,0,1,…,X-1,X(其中X為擬合記錄道的一半);c1和c2為一次多項(xiàng)式和二次多項(xiàng)式系數(shù);tk,i(0)為第k道時(shí)刻i擬合時(shí)，0表示擬合時(shí)的中心道，即第k道在所有擬合道的相對(duì)道號(hào);k=1,2,…,N(其中N為某CRP道集中的總道數(shù));i=1,2,…,L(其中L為單道樣點(diǎn)數(shù))。

通過掃描方式，可求得(1)式中的c1，c2，進(jìn)而求出相鄰各疊加道參與擬合的時(shí)間tk,i(x)，繼而可擬合出模型道。

文獻(xiàn)[14]基于去噪處理而提出，對(duì)參與擬合的數(shù)據(jù)采用中值約束下的矢量分解技術(shù)能夠很好地達(dá)到這一目的。實(shí)施過程中，在需向中值矢量的垂直向量作投影時(shí)，不可避免地保留部分噪聲。我們僅用常規(guī)多項(xiàng)式擬合的思想建立模型道，且僅以模型道為標(biāo)準(zhǔn)，確定對(duì)應(yīng)道集內(nèi)所有記錄道上有效信息的正確位置。對(duì)此，疊加結(jié)果較中值約束下的矢量分解結(jié)果更具有代表意義。

在這里采用疊加方式生成疊加模型道：

(2)

式中：S(tk,i)為疊加模型道;s(tk,i(x))為在疊加道x時(shí)刻tk,i(x)的樣點(diǎn)值;其它同(1)式。

1.1.2 分波形求取剩余時(shí)差及靜態(tài)校正

在文獻(xiàn)[13]中，對(duì)波形位置的校正采用分時(shí)窗的辦法實(shí)施。即對(duì)模型道和道集中記錄進(jìn)行分時(shí)窗處理，應(yīng)用互相關(guān)原理求取兩者的時(shí)差，再應(yīng)用拉伸/壓縮的辦法消除時(shí)差。這樣處理將導(dǎo)致以下缺點(diǎn)：

1) 消除時(shí)差是以時(shí)窗為整體，并不能嚴(yán)格校平道集內(nèi)每道記錄的有效信息。

2) 應(yīng)用拉伸/壓縮處理，將改變記錄原有頻率和相位特征，給反演帶來不利因素。

為此，將分時(shí)窗處理改為分波形處理。為此定義：每道記錄中兩次過零點(diǎn)的部分稱之為一個(gè)波形。時(shí)差校正以每個(gè)波形為單位進(jìn)行。于是有分波形求時(shí)差公式為

(3)

設(shè)在模型道上滿足式(3)的波形為S(tk,ii),ii∈(1,I1)，若以下(4)式成立

(4)

為避免采用拉伸/壓縮方法產(chǎn)生波形形變，本文采用靜態(tài)時(shí)移的方法消除時(shí)差。

1.2 基于橫向滑動(dòng)的時(shí)空變小波閾值保真去噪

目前，地震資料去噪技術(shù)大都采用多道識(shí)別、多道消除的技術(shù)思路。它利用相鄰道的信息來識(shí)別噪聲,并應(yīng)用相鄰道信息修正欲處理道信息。應(yīng)用這類技術(shù)后，增強(qiáng)相鄰道相似性，但去噪結(jié)果包含了相鄰道信息而損失各道自身個(gè)性，信號(hào)的保真性降低。多項(xiàng)式擬合、中值濾波、RNA(隨機(jī)噪聲衰減)、F-X反褶積等均屬于這類方法。在常規(guī)處理中，疊前去噪后的數(shù)據(jù)還要經(jīng)過疊加/偏移處理，故可不需特別關(guān)注道集內(nèi)各道自身個(gè)性。但疊前反演正是應(yīng)用這些個(gè)性進(jìn)行反演，因此，疊前反演前數(shù)據(jù)的去噪處理亟需單道保真去噪技術(shù)。

地震資料處理中的單道去噪技術(shù)包括帶通濾波、時(shí)頻峰值濾波及小波閾值等去噪技術(shù)。由于地震噪聲和有效信號(hào)在同頻帶上共存，帶通濾波的去噪能力有限;時(shí)頻峰值濾波因應(yīng)用魏格納一維利分布(Winger-Ville Distribution,WVD)產(chǎn)生交叉項(xiàng)的影響,去噪效果不佳;而小波閾值去噪受控于閾值的準(zhǔn)確性。鑒于目前沒有新的單道信噪分離算法，多道識(shí)別、單道消除技術(shù)成為目前實(shí)現(xiàn)保真去噪的一種選擇。

作為單道去噪技術(shù)的小波閾值去噪，其技術(shù)思路是對(duì)小波變換后小波系數(shù)(大的小波系數(shù)通常為信號(hào)，小的小波系數(shù)通常為噪聲)計(jì)算閾值進(jìn)行噪聲消除。實(shí)踐表明，僅采用單道技術(shù)很難求準(zhǔn)所需閾值?；诖耍谖墨I(xiàn)[15]中采用多道計(jì)算閾值、單道消除噪聲的方法進(jìn)行了相關(guān)疊加前的可控震源記錄噪聲消除。

在實(shí)踐中，常用的小波變換閾值去噪公式為：

(5)

方差為該小波系數(shù)與標(biāo)準(zhǔn)值的差的平方。實(shí)踐中，由于標(biāo)準(zhǔn)值難以求取，常取最小尺度下的小波系數(shù)的均值作為方差應(yīng)用。而地震資料噪聲具有時(shí)空變性質(zhì)，這種以噪聲均值為方差的方法十分粗糙，易產(chǎn)生過扼殺或過保留現(xiàn)象。文獻(xiàn)[15]提出的時(shí)空變閾值去噪技術(shù)能夠在相關(guān)疊加前的可控震源數(shù)據(jù)上應(yīng)用并獲得成功，源于它利用相關(guān)疊加前可控震源資料具有水平排齊、多道同幅同相的特點(diǎn)。文獻(xiàn)[15]將相關(guān)疊加前的整個(gè)可控震源道集內(nèi)所有記錄疊加，并以疊加結(jié)果作為標(biāo)準(zhǔn)值，進(jìn)而求取方差，使閾值求取不準(zhǔn)的問題得到改善。

但是，用于疊前反演的CMP道集通常有幾十甚至數(shù)百道，在整個(gè)道集內(nèi)各道記錄的振幅、頻寬、相位均存有很大變化，CMP道集難以具有水平排齊、多道同幅同相的共性。因此，對(duì)于參與疊前反演的CMP道集而言，直接應(yīng)用文獻(xiàn)[15]的方法不可能求準(zhǔn)時(shí)空變小波閾值。然而，在按偏移距選排的整個(gè)CRP道集內(nèi)，相鄰道記錄之間近似持有上述共性，可采用局部相鄰道疊加的方法來取代整個(gè)道集疊加，進(jìn)而求取時(shí)空變閾值。具體做法為：以某一道為中心，取其左右各N道組成一子道集，在子道集上實(shí)現(xiàn)時(shí)空變閾值去噪;然后橫向滑動(dòng)一道，再進(jìn)行相同的計(jì)算，直到整個(gè)CMP道集。由此，將文獻(xiàn)[15]的方法改進(jìn)為以下方法。

設(shè)σj,i(t)為在第j層尺度下，第i道時(shí)刻t的小波系數(shù)wj,i(t)的方差，則

(6)

式中：rj(t)是某道集在第j尺度下時(shí)刻t的數(shù)學(xué)期望(標(biāo)準(zhǔn)值);wj,i(t)為道集中的第j層尺度下第i道時(shí)刻t的小波系數(shù)?？扇∠噜徣舾傻赖寞B加結(jié)果，特別地還可利用加權(quán)疊加結(jié)果作為數(shù)學(xué)期望rj(t)，其權(quán)用方差的倒數(shù)取代，則某道集第j層尺度下時(shí)刻t的數(shù)學(xué)期望(標(biāo)準(zhǔn)值)可表示為：

(7)

式中，2N+1為第i道相鄰兩邊的求和道數(shù)。顯然，若方差求取準(zhǔn)確，就可求得準(zhǔn)確可靠的加權(quán)疊加結(jié)果。反之，若疊加結(jié)果最佳，σj,i(t)也準(zhǔn)確。

對(duì)于上式，釆用迭代方式來不斷修正σj,i(t)，即先給出一組初始σj,i(t)，代入上式計(jì)算疊加和rj(t)，則剩余噪聲方差為：

(8)

從記錄道起點(diǎn)開始，逐點(diǎn)可計(jì)算出每道每點(diǎn)的方差。這樣便完成了方差的時(shí)空變計(jì)算。將方差σj,i(t)代入(4)式即可完成隨機(jī)噪聲的消除。由上可見：準(zhǔn)確求出方差σj,i(t)時(shí)，應(yīng)用相鄰道小波系數(shù)，但去噪僅應(yīng)用σj,i(t)。這就實(shí)現(xiàn)了具有保真性的多道識(shí)別、單道消除思想。

將求取的σj,i(t)代入(5)式便實(shí)現(xiàn)了小波閾值去噪，然后向下進(jìn)行下一樣點(diǎn)的去噪處理，一直到該道最后一個(gè)樣點(diǎn)。滑動(dòng)一道進(jìn)行同樣處理，直到該道集內(nèi)最后一道。對(duì)所有各層小波系數(shù)進(jìn)行上述處理，然后小波重建便完成了橫向滑動(dòng)的時(shí)空變小波閾值保真去噪。

2 模型試算及實(shí)際資料處理

2.1 理論模型數(shù)據(jù)試算

圖1為應(yīng)用本文提出的疊前反演數(shù)據(jù)優(yōu)化處理技術(shù)對(duì)理論模型數(shù)據(jù)試算的結(jié)果。圖1a是信噪比為1.8的疊前道集，其中第8道和第9道下移了10ms，且信號(hào)能量為相鄰道的0.6倍，以模擬具有時(shí)延、能量變化的含噪道集。圖1b為應(yīng)用本文基于常規(guī)多項(xiàng)式擬合的剩余時(shí)差校正技術(shù)消除剩余時(shí)差后的結(jié)果，可見處理后時(shí)差得以消除。圖1c 是對(duì)圖1b進(jìn)行F-X反褶積去噪處理結(jié)果，圖1b中的隨機(jī)噪聲得到了有效衰減，但第8道和第9道的能量經(jīng)處理后變得和相鄰道相近，這是多道識(shí)別、多道消除產(chǎn)生的結(jié)果，這個(gè)結(jié)果是通過多道信

息擬合出各道記錄的振幅，破壞了道集內(nèi)各道記錄的相對(duì)振幅關(guān)系，將給反演結(jié)果帶來致命的錯(cuò)誤。圖1d是應(yīng)用橫向滑動(dòng)時(shí)空變小波閾值保真去噪技術(shù)處理后的結(jié)果，可見去噪結(jié)果較圖1c更為徹底，圖中幾乎沒有噪聲存在，更為可喜的是第8道和第9道保持著原來的能量，與相鄰道能量關(guān)系與去噪前保持一致。

2.2 實(shí)際地震資料處理

圖2給出了實(shí)際地震道集分時(shí)窗和分波形時(shí)差校正處理的效果對(duì)比。圖2a為未進(jìn)行時(shí)差校正的實(shí)際地震資料CMP道集，圖2b為分時(shí)窗校正處理結(jié)果，而圖2c為分波形校正處理結(jié)果。對(duì)比可見，分波形校正較分時(shí)窗校正的時(shí)差校正效果明顯更好。

圖1 理論模型數(shù)據(jù)應(yīng)用本文提出技術(shù)優(yōu)化處理的結(jié)果a 具有時(shí)延和能量變化的含噪疊前道集; b 應(yīng)用本文提出的優(yōu)化處理技術(shù)消除剩余時(shí)差后的結(jié)果; c 對(duì)圖1b進(jìn)行F-X反褶積去噪處理結(jié)果; d 對(duì)圖1b應(yīng)用本文提出的優(yōu)化處理技術(shù)去噪結(jié)果

圖2 實(shí)際地震資料CMP道集(a)的分時(shí)窗(b)和分波形(c)時(shí)差校正處理效果對(duì)比

圖3為對(duì)實(shí)際地震資料含噪道集進(jìn)行不同方法去噪處理的效果對(duì)比。圖3a為含噪道集;圖3b為F-X反褶積去噪結(jié)果;圖3c為文獻(xiàn)[15]的時(shí)空變小波閾值去噪結(jié)果;圖3d為本文橫向滑動(dòng)時(shí)空變小波閾值去噪結(jié)果。F-X反褶積是基于多道預(yù)測(cè)、多道消除思路的去噪技術(shù)，處理時(shí)用相鄰道的信息擬合處理道信息，嚴(yán)重地?fù)p害了相鄰道之間的相對(duì)振幅關(guān)系及波形特征(見圖3b)。而應(yīng)用文獻(xiàn)[15]提出的時(shí)空變小波閾值去噪后(見圖3c)，各道記錄在同一時(shí)間上的相似性增強(qiáng)，破壞了各道記錄的相對(duì)能量關(guān)系及波形特征。由圖3d可見，應(yīng)用本文提出的橫向滑動(dòng)時(shí)空變小波閾值去噪結(jié)果背景干凈，隨機(jī)噪聲得到了有效消除，且很好地保持了有效信息的波形形態(tài)和相對(duì)能量關(guān)系。

圖4為應(yīng)用本文提出的技術(shù)優(yōu)化處理前、后近道疊加剖面(圖中黑色)和遠(yuǎn)道疊加剖面(圖中紅色)疊合的效果對(duì)比。圖4a為優(yōu)化處理前的道集疊加剖面;圖4b為優(yōu)化處理后的道集疊加剖面;圖4c 和圖4d分別為對(duì)應(yīng)圖4a和圖4b中橢圓框的局部放大顯示。由圖4可見，經(jīng)疊前優(yōu)化處理后，近、遠(yuǎn)道疊加結(jié)果更趨重合，這是由于優(yōu)化處理消除了剩余時(shí)差和噪聲，使近道和遠(yuǎn)道記錄的波形趨于一致，進(jìn)而使用于疊前反演分析的反射波振幅來自于同一時(shí)間。

圖3 實(shí)際地震資料含噪道集的不同方法去噪處理效果對(duì)比a 含噪道集; b F-X反褶積去噪結(jié)果; c 時(shí)空變小波閾值去噪結(jié)果; d 橫向滑動(dòng)時(shí)空變小波閾值去噪結(jié)果

圖4 優(yōu)化處理前(a)、后(b)近道疊加剖面(黑色)和遠(yuǎn)道疊加剖面(紅色)疊合及其對(duì)應(yīng)的局部放大顯示(c和d)

圖5為應(yīng)用本文技術(shù)優(yōu)化處理前、后道集的AVO響應(yīng)效果對(duì)比。圖5a至圖5d分別為優(yōu)化處理前的CRP道集、正演道集、3°～30°入射角的疊加剖面(圖中插入了測(cè)井曲線)及AVO響應(yīng);圖5e 至圖5h則分別為對(duì)應(yīng)以上各圖的優(yōu)化處理后結(jié)果。對(duì)比優(yōu)化處理前(圖5a)、后(圖5e)的CRP道集，可見圖5e中有效信號(hào)的同相軸校平，層間噪聲得以消除。結(jié)合地層解釋橫線對(duì)比優(yōu)化處理前(圖5b)、后(圖5f)的正演道集，圖5f中的每一解釋層位均與疊加剖面中插入的測(cè)井曲線一致，而圖5b則層位對(duì)比困難。圖5d和圖5h中粉紅線是優(yōu)化前、后正演道集反射振幅隨角度變化趨勢(shì)；而紅線是優(yōu)化前后實(shí)際道集的反射振幅隨角度變化趨勢(shì)。從優(yōu)化處理前(圖5d)、后(圖5h)AVO響應(yīng)的對(duì)比可以看出，經(jīng)優(yōu)化處理后，兩種趨勢(shì)曲線在縱向上更靠近，且沒有改變二者的變化趨勢(shì)。顯然，圖5h的AVO現(xiàn)象更為明顯。

圖5 應(yīng)用本文技術(shù)優(yōu)化處理前、后道集的AVO響應(yīng)效果對(duì)比a 優(yōu)化處理前的CRP道集; b 優(yōu)化處理前的正演道集; c 優(yōu)化處理前的3°～30°入射角疊加剖面; d 優(yōu)化處理前的AVO振幅曲線; e 優(yōu)化處理后的CRP道集; f 優(yōu)化處理后的正演道集; g 優(yōu)化處理后的3°～30°入射角疊加剖面; h 優(yōu)化處理后的AVO振幅曲線

3 結(jié)束語

通過上述理論論述和模型數(shù)據(jù)測(cè)試及實(shí)際地震資料處理驗(yàn)證，可以得出以下結(jié)論。

1) 常規(guī)地震資料處理流程及部分處理技術(shù)已不能滿足疊前反演的需要，疊前反演前必須進(jìn)行數(shù)據(jù)優(yōu)化處理。

2) 有效信息排齊及盡可能地保真去噪是疊前反演數(shù)據(jù)的最基本需求。

3) 利用常規(guī)多項(xiàng)式擬合技術(shù)校正疊加結(jié)果，保證了疊加結(jié)果位置的準(zhǔn)確性;由分波形時(shí)差校正取代分時(shí)窗時(shí)差校正，提高了剩余時(shí)差校正的精度。

4) 多道識(shí)別、單道消除的橫向滑動(dòng)時(shí)空變小波閾值去噪技術(shù)能夠保證去噪過程中的保真性。

5) 利用本文提出的時(shí)差校正及保真去噪技術(shù)實(shí)現(xiàn)了疊前反演數(shù)據(jù)優(yōu)化，在理論模型數(shù)據(jù)及實(shí)際地震資料處理中的成功應(yīng)用證明了這兩項(xiàng)技術(shù)的正確性及適應(yīng)性。

參 考 文 獻(xiàn)

[1] Veeken P C H,Silva M D.Seismic inversion methods and some of their constrains[J].First Break,2004,22(1):47-70

[2] Whitcombe D N.Elastic impedance normalization[J].Geophysics,2002,58(1):60-62

[3] Connolly P.Elastic impedance[J].The Leading Edge,1999,18(4):438-452

[4] 張永剛.地震波阻抗反演技術(shù)的現(xiàn)狀和發(fā)展[J].石油物探,2002,41(4):385-390

Zhang Y G.The present and future of wave impedance inversion technique[J].Geophysical Prospecting for Petroleum,2002,41(4):385-390

[5] 管路平.地震疊前反演與直接烴類指示的探討[J].石油物探,2008,47(3):228-234

Guan L P.Pre-stack inversion and direct hydrocarbon indicator[J].Geophysical Prospecting for Petroleum,2008,47(3):228-234

[6] 仝兆岐.儲(chǔ)層地震技術(shù)新進(jìn)展[M].東營:石油大學(xué)出版社,2004:98-135

Tong Z Q.The New Progress of The reservoir Seismic Technique[M].DongYing:Petroleum university Press,2004:98-135

[7] 苑書金.疊前地震反演技術(shù)的進(jìn)展及其在巖性油氣藏勘探中的應(yīng)用[J].地球物理學(xué)進(jìn)展,2007,22(3):879-886

Yuan S J.Progress of pre-stack inversion and application in exploration of the lithological reservoir[J].Process in Geophysics,2007,22(3):879-886

[8] 楊迪琨,胡祥云.含噪聲數(shù)據(jù)反演的概率描述[J].地球物理學(xué)報(bào),2008,51(3):901-907

Yang D K,Hu X Y.Inversion of noisy data by probabilistic methodology[J].Chinese Journal of Geophysics(in Chinese),2008,51(3):901-907

[9] 熊定鈺,趙海珍,陳海云，等.保持地震記錄疊前AVO屬性的噪聲衰減方法[J].石油地球物理勘探，2010,45(6):856-860

Xiong D Y，Zhao H Z，Chen H Y,et al.Noise attention method based on the preserved preserved prestack AVO attribute[J].Oil Geophysical Prospecting,2010,45(6):856-860

[10] 鮑熙杰,趙海波.CRP道集優(yōu)化處理對(duì)疊前反演的影響[J].石油地質(zhì)與工程，2013,27(1):39-41

Bao X J,Zhao H B.effects analysis of CRP gather optimal processing on prectack inversion[J].Exploration and Development Institute，2013,27(1):39-41

[11] Singleton S．The effects of seismic data conditioning on prestack simultaneous impedance inversion[J].The Leading Edge,2009,28(1):772-781

[12] 劉力輝,楊曉,丁燕，等.基于巖性預(yù)測(cè)的CRP道集優(yōu)化處理[J].石油物探，2013,52(5):482-488

Liu L H,Yang X,Ding Y,et al.CRP gather optimization processing based on lithological prediction[J].Geophysical Prospecting for Petroleum,2013,52(5):482-488

[13] 夏洪瑞,陳德剛，周開明.剩余動(dòng)校正量的拾取與消除[J].石油地球物理勘探,1997,32(6):872-877

Xia H R,Chen D G,Zhou K M.Pickup and removal of residual normal moveout[J].Oil Geophysical Prospecting,1997,32(6):872-877

[14] 夏洪瑞,董江偉,鄒少峰,等.常規(guī)二次多項(xiàng)式擬合地震數(shù)據(jù)[J].石油物探，2006,45(5):492-496

Xia H R,Dong J W,Zou S F,et al.Normal quadratic polynomial fitting seismic data[J].Geophysical Prospecting for Petroleum,2006,45(5):492-496

[15] 夏洪瑞,葛川慶,彭濤.小波時(shí)空變閾值去噪方法在可控震源資料處理中的應(yīng)用[J].石油地球物理勘探，2010,45(01):23-27

Xia H R,Ge C Q,Peng T.Application of wavelet time-space-varying threshold denoising method in vibroseis seismic data processing[J].Oil Geophysical Prospecting,2010,45(1):23-27