陳德元

（1.國(guó)家現(xiàn)代地質(zhì)勘查工程技術(shù)研究中心，河北 廊坊065000；2.中國(guó)地質(zhì)科學(xué)院地球物理地球化學(xué)勘查研究所，河北 廊坊065000）

1 引言

未加約束條件的反演呈現(xiàn)出能量分散特征，添加約束條件之后的反演結(jié)果難以真實(shí)反映真實(shí)反射系數(shù)，主要是由于Cauchy 約束會(huì)壓制反射系數(shù)，L2 模約束條件的運(yùn)用能夠提升算法的精確性。這一算法在實(shí)踐運(yùn)用中得到了良好的運(yùn)用效果，提升了地震勘探的分辨率。

2 譜反演技術(shù)概述

譜反演技術(shù)運(yùn)用中能夠?qū)?duì)調(diào)諧厚度的薄層進(jìn)行有效的成像處理，是一種采用頻譜地分解技術(shù)。反射系數(shù)序列奇偶分解能夠顯著提升分辨率，提出建立譜反演的目標(biāo)函數(shù)，提出對(duì)薄層反射系數(shù)大小、極性與位置等分析中具有重要運(yùn)用價(jià)值。運(yùn)用中可以不建立反射系數(shù)的數(shù)學(xué)假設(shè)，也可以不建立先驗(yàn)?zāi)Ｐ停贿M(jìn)行任何層位約束與測(cè)井資料強(qiáng)制約束，在低于調(diào)諧厚度的地震薄層分辨上具有顯著優(yōu)勢(shì)[1]。

3 譜反演技術(shù)分析

3.1 建立譜反演目標(biāo)函數(shù)

譜反演運(yùn)用過(guò)程中能夠避免地震子波對(duì)地震記錄的不良影響，運(yùn)用地震子波與地震記錄進(jìn)行數(shù)據(jù)分析，以此得到反射系數(shù)序列，進(jìn)行地震數(shù)據(jù)分析。結(jié)合函數(shù)奇偶的分解理論，可以將反射系數(shù)分解為奇反射系數(shù)對(duì)與偶反射系數(shù)對(duì)，進(jìn)行兩者之間的對(duì)比。

通過(guò)對(duì)奇反射系數(shù)對(duì)的模擬研究結(jié)果分析可見(jiàn)，隨著厚度的降低，振幅呈現(xiàn)出先增大后減小的現(xiàn)象。偶反射系數(shù)對(duì)的研究則表現(xiàn)出相反的特征，隨著厚度的降低，振幅呈現(xiàn)出先減小后增大的現(xiàn)象，在實(shí)際的分析過(guò)程中一般包括2 個(gè)奇偶反射系數(shù)，對(duì)地震響應(yīng)分析能夠提供的數(shù)據(jù)支持。

3.2 譜反演目標(biāo)函數(shù)線性化

譜反演目標(biāo)函數(shù)中具有多項(xiàng)未知元素，包括反射系數(shù)大小、位置等，不同因素之間互相影響，研究中具有很大的不確定性以及多解性，算法的全局尋優(yōu)性能對(duì)精確求解具有重要影響作用。為了便于研究，可以采用線性方程組研究。不同求解反射系數(shù)位置，利用線性方程能夠得到反射系數(shù)的反射系數(shù)與奇反射系數(shù)，通過(guò)對(duì)其重構(gòu)處理之后能夠得出反射系數(shù)的解。求解可以采用最小二乘QR 分解、LSCG 方法與SVD 方法與共軛梯度法等。

4 線性譜反演的局限性

反射系數(shù)頻譜運(yùn)用中使用0～Nyquist 頻率的地震數(shù)據(jù)進(jìn)行反射系數(shù)的反演。研究中若假設(shè)每個(gè)時(shí)間采樣點(diǎn)對(duì)應(yīng)一個(gè)反射系數(shù)值，則實(shí)際計(jì)算量遠(yuǎn)遠(yuǎn)大于真實(shí)反射系數(shù)個(gè)數(shù)，對(duì)計(jì)算提出了更高的要求。

將不同頻率范圍內(nèi)反演的結(jié)果運(yùn)用30Hz 雷克子波進(jìn)行地震記錄合成，結(jié)果顯示其具有較為復(fù)雜的反射系數(shù)，具有薄互層。具有4ms 采樣間隔，在0～Nyquist125Hz 頻率內(nèi)進(jìn)行反演，對(duì)數(shù)據(jù)分析可知其具有較為良好的反演效果，能夠較好地反演薄層，分析5～95Hz 頻率反演內(nèi)反演效果可見(jiàn)，得到了較為分散的能量反演結(jié)果，在真實(shí)反射系數(shù)位置周?chē)霈F(xiàn)不同程度的抖動(dòng)現(xiàn)象，反演頻帶寬度與反演效果呈正比關(guān)系，反演頻帶寬度越高，則能夠得到越好的反演效果，隨著能量的逐漸分散，其中的分辨率也逐漸降低。

在f∈［0，125Hz］、125Hz Nyquist 頻率、4ms 間隔采樣時(shí)，能夠得到較為良好的反演效果。改變反演條件，設(shè)置f∈［0，250Hz］、50Hz Nyquist、2ms 采樣間隔，研究其相應(yīng)的反演效果。結(jié)果顯示具有一定的誤差，反演實(shí)際結(jié)果較為有限。將反射系數(shù)與真實(shí)反射系數(shù)綜合進(jìn)行傅立葉變換，分析其中的頻譜，結(jié)果可見(jiàn)，反射系數(shù)中高頻區(qū)域與真實(shí)數(shù)據(jù)之間存在較大誤差。

分析這些不同數(shù)據(jù)值，可見(jiàn)高頻成分的幅值較接近于0，這與地震子波具有帶限特征之間有著緊密聯(lián)系。S（f）/W（f）求解過(guò)程中具有很大的不確定性，構(gòu)建求解反射系數(shù)方程組進(jìn)行分析。2ms 采樣間隔時(shí)，由于存在一個(gè)隨機(jī)復(fù)數(shù)，使得難以充分反演所有頻率成分。

在具體的數(shù)據(jù)分析過(guò)程中，由于全部的地震資料均帶限，因此，在反演反射系數(shù)過(guò)程中難以利用所有的頻率成分，在分析過(guò)程中應(yīng)當(dāng)將高頻成分、低頻成分去除，可見(jiàn)線性譜反演過(guò)程中難以利用帶限的地震資料得出精確解。為此在研究中可以假定反射系數(shù)稀疏，即非地層介質(zhì)波阻抗分界面反射系數(shù)是0，而地層介質(zhì)波阻抗分界面反射系數(shù)不是0，在這一前提之下進(jìn)行數(shù)值研究，線性譜反演方法運(yùn)用中將每個(gè)采樣點(diǎn)看作一個(gè)反射界面，默認(rèn)反射系數(shù)非稀疏，因此，研究過(guò)程中需要對(duì)反射系數(shù)進(jìn)行稀疏約束。

5 線性譜反演約束方法

地震資料由于帶限，在反演過(guò)程中難以充分利用全部的頻率，線性譜反演假定反射系數(shù)研究中為了得到精確解，要求具有所有頻率成分，這是由于其在每個(gè)時(shí)間采樣點(diǎn)都有存在值。可以約束線性反演條件以解決這一問(wèn)題，進(jìn)行Cauchy 約束并反演，并篩選若干反射系數(shù)，要求應(yīng)當(dāng)高于或者等于真實(shí)反射系數(shù)個(gè)數(shù)，均為已知并具有稀疏特征，利用L2 模約束反演反射系數(shù)。

采用Cauchy 約束的稀疏約束，具有顯著的抗噪與反射系數(shù)位置檢測(cè)效果，運(yùn)用非線性隱式格式，迭代求解過(guò)程中運(yùn)用最小二乘QR 分解法進(jìn)行求解。所求值與真實(shí)值之間具有一定差異，數(shù)值偏小。為此運(yùn)用中篩選若干反射系數(shù)，以此表示真實(shí)反射系數(shù)位置，并反演這些反射系數(shù)，為了得到真實(shí)解，提升算法的穩(wěn)定性，進(jìn)行L2 模約束。

6 分析理論模型試驗(yàn)

設(shè)置模型試驗(yàn)，采用2ms 的時(shí)間采樣間隔，為模型記錄合成地震數(shù)據(jù)并設(shè)置反射系數(shù)。設(shè)置5～150Hz 的頻率范圍，對(duì)反演設(shè)置Cauchy 約束，得到相應(yīng)的反演結(jié)果。將其與未設(shè)置約束條件的反演結(jié)果進(jìn)行對(duì)比，結(jié)果顯示，未加約束條件的反演呈現(xiàn)出能量分散特征，表現(xiàn)為真實(shí)反射系數(shù)周?chē)霈F(xiàn)不同程度的反演，而在添加約束條件之后則能夠促進(jìn)能量的集中，在反射界面中能夠聚集更多的能量。同時(shí)應(yīng)當(dāng)看到，添加約束條件之后的反演結(jié)果難以真實(shí)反映真實(shí)反射系數(shù)，這種誤差的出現(xiàn)主要是由于Cauchy 約束會(huì)壓制反射系數(shù)。

分析Cauchy 約束反演與L2 模約束共同反演之后得到的80 個(gè)反射系數(shù)結(jié)果，可見(jiàn)有效縮短了其中的誤差，已經(jīng)近似于真實(shí)反射系數(shù)。研究10～100Hz 頻率范圍內(nèi)的反演結(jié)果，可見(jiàn)在反演過(guò)程中通過(guò)減少低頻與高頻部分，能夠在很大程度上降低誤差，得出的反演結(jié)果較為理想。但是其誤差數(shù)值依然比5～130Hz 頻率范圍內(nèi)的誤差值要高，通過(guò)這2 種不同頻率范圍反演結(jié)果對(duì)比可知，反演頻率與反演效果精度呈正比，參與反演頻率成分越多，則能夠得到越位精確的反演效果。

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

在實(shí)際疊后地震資料中運(yùn)用譜反演方法，選擇某工區(qū)的實(shí)際地震剖面，具有8000ms 地震道長(zhǎng)度與1ms 采樣率。通過(guò)統(tǒng)計(jì)法得出雷克子波，地震波在不同的時(shí)間段內(nèi)具有高頻衰減現(xiàn)象（見(jiàn)圖1 和圖2）。因此，不同時(shí)間段對(duì)應(yīng)相應(yīng)的子波。通過(guò)數(shù)據(jù)分析可見(jiàn)反射系數(shù)剖面與真實(shí)地層之間一一對(duì)應(yīng)。在對(duì)地震剖面分析中采用肉眼難以觀測(cè)，反射系數(shù)剖面對(duì)人體肉眼觀測(cè)進(jìn)行了有效互補(bǔ)，能夠得出關(guān)于地震剖面的有效斷層細(xì)節(jié)與地層細(xì)節(jié)，以此促進(jìn)地質(zhì)分析中的精細(xì)化解釋。

圖1 不同頻率范圍的反演結(jié)果

圖2 原始地震剖面與反射系數(shù)不同

譜反演分析方式運(yùn)用中在補(bǔ)償高頻成分的同時(shí)，還能夠保留中低頻成分。將反演子波與真實(shí)反射系數(shù)褶積之后得出合成地震記錄，運(yùn)用藍(lán)色標(biāo)記地震道，地震記錄合成過(guò)程中需要運(yùn)用反演子波褶積與反射系數(shù)，運(yùn)用紅色標(biāo)記進(jìn)行記錄，以黑色標(biāo)記表示實(shí)際地震道。兩者之間具有很大的共性，基本一致，可見(jiàn)反演方式能夠?qū)?shí)際地層有限帶寬信息進(jìn)行較為真實(shí)的反映與記錄。將其與真實(shí)反射系數(shù)在褶積之后得到合成記錄，采用藍(lán)色標(biāo)記表示，運(yùn)用紅色標(biāo)記子波與反射系數(shù)褶積之后得到的合成地震記錄，可見(jiàn)譜反演方式的運(yùn)用有效補(bǔ)償了較弱的高頻信息，能夠與測(cè)井?dāng)?shù)據(jù)之間達(dá)到良好的對(duì)應(yīng)效果。

8 結(jié)語(yǔ)

譜反演技術(shù)運(yùn)用過(guò)程中不需要先驗(yàn)信息約束與初始模型，屬于高分辨率反射系數(shù)成像技術(shù)，本文研究方式能夠有效分辨薄層，同時(shí)顯著提升地震頻帶提高分辨率，在對(duì)地震相關(guān)數(shù)據(jù)分析中運(yùn)用效果顯著，研究中主要考慮地震子波的時(shí)變特征，在對(duì)子波的空變特征的分析方面不足，在下一步的反演過(guò)程中可以綜合考慮這一方面的因素，達(dá)到更好的譜反演效果。