陳 平， 譚輝煌， 秦 童， 閆 濤

( 中海石油(中國(guó))有限公司 天津分公司，天津 300459 )

0 引言

地震資料的分辨率是巖性、構(gòu)造—巖性油氣藏勘探中儲(chǔ)層描述的關(guān)鍵，從測(cè)井?dāng)?shù)據(jù)中獲得更高分辨率的譜藍(lán)化方法不依賴模型且實(shí)現(xiàn)便捷，在高分辨率儲(chǔ)層預(yù)測(cè)、有色波阻抗反演等研究應(yīng)用中取得較好效果。WALDEN A T等[1]發(fā)現(xiàn)井上反射系數(shù)藍(lán)譜特征，提出利用一階自回歸滑動(dòng)平均模型(ARMA)擬合反射系數(shù)的藍(lán)色項(xiàng)；趙波等[2]將藍(lán)色項(xiàng)算子化應(yīng)用在譜模擬反褶積中，彌補(bǔ)常規(guī)反褶積無(wú)法得到有色反射系數(shù)的缺陷；BLACHE-FRASER G[3]將井上藍(lán)色項(xiàng)算子直接應(yīng)用到地震或反演數(shù)據(jù)，以改善砂組單元和氣水界面的地震響應(yīng)，形成譜藍(lán)化方法(Spectrum Blueing, SB)；KAZEMEINI S H等[4]將譜藍(lán)化方法推廣應(yīng)用到疊前。近些年，彭更新等[5]、楊瑞召等[6]、紀(jì)甜甜等[7]、陳文雄等[8]、李賢兵等[9]分別利用譜藍(lán)化方法提高薄互層或隔層的識(shí)別能力，以滿足儲(chǔ)層描述的需求，取得一定的效果。人們對(duì)反射系數(shù)的統(tǒng)計(jì)特性討論較多(頻譜、直方圖、自相關(guān)等)，但缺少相應(yīng)的正演研究，在譜藍(lán)化方法應(yīng)用前還需厘清藍(lán)譜特征和沉積儲(chǔ)層的內(nèi)在聯(lián)系，為譜藍(lán)化提高地震資料的薄儲(chǔ)層響應(yīng)能力奠定應(yīng)用基礎(chǔ)。

譜藍(lán)化方法本質(zhì)上應(yīng)對(duì)白化的反射系數(shù)譜進(jìn)行有色改造[2]。在未盡量消除子波影響或原始地震資料頻帶窄時(shí)，將井上反射系數(shù)的藍(lán)色項(xiàng)直接作用于原始地震數(shù)據(jù)，可能造成地震資料的低頻缺失及旁瓣振蕩，出現(xiàn)假的分辨率。同時(shí)，譜反演是在頻率域求解反射系數(shù)、壓縮子波的線性反演方法。當(dāng)反演目標(biāo)函數(shù)為平方差時(shí)，是線性規(guī)劃問(wèn)題[10]，當(dāng)目標(biāo)函數(shù)是目標(biāo)范式時(shí)，是壓縮感知問(wèn)題[11]。二者求解的算法不同，都要估算子波，只能得到接近白譜的拓頻結(jié)果，無(wú)法反映實(shí)際地層結(jié)構(gòu)下薄層的真實(shí)響應(yīng)。

筆者分析渤海萊州灣地區(qū)已鉆井反射系數(shù)特征，結(jié)合譜藍(lán)化和譜反演的特點(diǎn)，提出基于譜反演的譜藍(lán)化拓頻方法及技術(shù)流程；建立不同地層結(jié)構(gòu)反射系數(shù)模型，探討沉積儲(chǔ)層對(duì)反射系數(shù)譜的影響，采用譜反演的譜藍(lán)化拓頻方法能夠提高薄層響應(yīng)，在同期復(fù)合砂體儲(chǔ)層描述中具有應(yīng)用優(yōu)勢(shì)。

1 基本原理

1.1 反射系數(shù)藍(lán)譜特征分析

渤海萊州灣地區(qū)新近系為典型的砂泥巖不等厚沉積[12-13]，統(tǒng)計(jì)兩口實(shí)鉆井反射系數(shù)(見圖1)。概率擬合曲線偏離高斯分布，具有近似拉普拉斯的分布特征(見圖2(a、c))，即稀疏性較強(qiáng)，反射系數(shù)集中于零值附近。頻譜在地震頻帶內(nèi)(0～200 Hz)低頻弱、高頻強(qiáng)，隨頻率增加呈對(duì)數(shù)形態(tài)上升，形成有色的藍(lán)譜(見圖2(b、d))。

圖1 渤海萊州灣實(shí)鉆井1 000～2 000 m反射系數(shù)Fig.1 Reflection coefficients of drilled wells from 1 000 to 2 000 m in Laizhou Bay, Bohai Sea

圖2 渤海萊州灣實(shí)鉆井反射系數(shù)振幅、頻率統(tǒng)計(jì)特征Fig.2 Characteristics of amplitude and frequency of drilled wells in Laizhou Bay, Bohai Sea

人們對(duì)反射系數(shù)譜形態(tài)特征研究較多，對(duì)藍(lán)譜與沉積儲(chǔ)層相關(guān)關(guān)系分析較少。WALDEN A T等分析反射系數(shù)自相關(guān)譜，認(rèn)為非對(duì)稱旁瓣由反射系數(shù)高頻產(chǎn)生，自相關(guān)譜的正負(fù)峰對(duì)頻譜起到導(dǎo)數(shù)算子作用，從而引起斜坡效應(yīng)[1]，但沒(méi)有探討沉積儲(chǔ)層與藍(lán)譜的內(nèi)在聯(lián)系。反射系數(shù)與工業(yè)上的有色噪聲具有相似性，即譜密度有明顯的峰值且自相關(guān)非單位脈沖[14-15]，可認(rèn)為反射系數(shù)與有色噪聲類似具有時(shí)間序列上的相關(guān)關(guān)系。KL-3、KL-6井反射系數(shù)及高斯白色反射系數(shù)的自相關(guān)譜見圖3。由圖3可見，高斯白色反射系數(shù)表現(xiàn)為尖脈沖，表明任意時(shí)刻不具有相關(guān)關(guān)系，井上反射系數(shù)在-0.4～1.0之間具有多個(gè)極值，說(shuō)明井上反射系數(shù)具有較強(qiáng)的相關(guān)關(guān)系。渤海萊洲灣地區(qū)新近系阻抗區(qū)分性較好，砂巖普遍呈低速、低密特征，與背景泥巖有較好的阻抗區(qū)分性[16-17]，形成明顯的反射界面，儲(chǔ)層的頂界面對(duì)應(yīng)負(fù)反射系數(shù)，儲(chǔ)層的底界面對(duì)應(yīng)正反射系數(shù)，二者振幅強(qiáng)度一致。一定地質(zhì)時(shí)期內(nèi)，儲(chǔ)層或泥巖沉積較為穩(wěn)定，使反射系數(shù)產(chǎn)生時(shí)間序列的相關(guān)關(guān)系，在自相關(guān)關(guān)系上出現(xiàn)峰—谷特征。

圖3 渤海萊洲灣實(shí)鉆井反射系數(shù)的自相關(guān)譜Fig.3 Autocorrelation spectrum of reflection coefficient of drilled wells in Laizhou Bay, Bohai Sea

井上反射系數(shù)服從拉普拉斯分布，位置存在時(shí)間上的關(guān)聯(lián)，同時(shí)藍(lán)譜與地層厚度可能存在聯(lián)系。設(shè)計(jì)4種地層模型(見表1、圖4)，反射系數(shù)服從拉普拉斯分布。當(dāng)?shù)貙咏Y(jié)構(gòu)無(wú)規(guī)律時(shí)，地層模型為服從拉普拉斯分布的隨機(jī)序列(見圖4(a))。當(dāng)?shù)貙訛樯澳嘟换r(shí)，模型2、3、4分別對(duì)應(yīng)地層厚度服從完全隨機(jī)、高斯隨機(jī)和拉普拉斯隨機(jī)(見圖4(b-d))。由圖4可見，模型1的頻譜近似白譜，模型2的頻譜出現(xiàn)斜坡效應(yīng)，但與井上藍(lán)譜不一致(見圖4(b-c))。模型4的頻譜最接近井上反射系數(shù)的頻譜形態(tài)(見圖4(d))。反射系數(shù)模型說(shuō)明，時(shí)間序列相關(guān)關(guān)系只能起到斜坡效應(yīng)，壓制反射系數(shù)一部分低頻，而地層厚度變化主導(dǎo)頻譜的形態(tài)，地層以薄層為主時(shí)出現(xiàn)井上一致的藍(lán)譜形態(tài)。

表1 渤海萊州灣實(shí)鉆井反射系數(shù)模型

圖4 不同地層模型的反射系數(shù)序列和頻譜Fig.4 Reflection coefficients and amplitude spectrum of different models

1.2 譜藍(lán)化拓頻方法和技術(shù)流程

將地震反射系數(shù)看成白色和有色的反射系數(shù)項(xiàng)的褶積[2]，不含噪聲地震道的褶積模型表示為

x=r*b=rw*rnw*b,

(1)

式中：r為實(shí)際反射系數(shù)；rw為白色反射系數(shù)；rnw為有色反射系數(shù)；x為地震道；b為子波。由于地震子波頻帶較窄，通常原始地震頻譜的有色特征并不明顯。常采用ARMA擬合井上有色反射系數(shù)項(xiàng)，應(yīng)用于地震資料的譜藍(lán)化處理，模擬有色反射系數(shù)z的變換為

F(z)=(1-αz)/(1-βz),

(2)

式中：α和β分別為自回歸、滑動(dòng)平均經(jīng)驗(yàn)參數(shù)。

F(z)的反z變換為有色反射系數(shù)的時(shí)間序列，將零相位處理后為譜藍(lán)化算子，能增強(qiáng)井震反射系數(shù)的相關(guān)關(guān)系，提高地震資料對(duì)薄互層或隔層的識(shí)別能力。譜藍(lán)化算子頻率域響應(yīng)近似為對(duì)數(shù)形態(tài)，本質(zhì)上是一個(gè)壓低頻、抬高頻的過(guò)程，當(dāng)原始地震資料頻帶較窄或本身欠缺低頻時(shí)，可以改造地震子波。藍(lán)化后的子波能提高視覺分辨率，子波旁瓣振蕩更嚴(yán)重，相對(duì)阻抗可能出現(xiàn)假的分辨率。利用10～85 Hz帶通子波及藍(lán)化子波(見圖5(a-b))設(shè)計(jì)頂?shù)追瓷湎禂?shù)相反的二維模型(見圖5(c-d))，及其合成記錄對(duì)應(yīng)的-90°相移剖面(見圖5(e-f))。由圖5(b)可見，子波峰谷比更小、旁瓣振蕩更為劇烈；由圖5(f)可見，出現(xiàn)假的同相軸，容易落入解釋誤區(qū)。因此，地震譜藍(lán)化不能簡(jiǎn)單地直接應(yīng)用藍(lán)化算子，應(yīng)先消除子波的帶限影響，拓寬原始地震的頻帶。

譜反演(Spectrum Inversion, SI)不同于常規(guī)反褶積技術(shù)[18-21]，在稀疏解的約束下，通過(guò)迭代反演的方式避免確定性反褶積算子的求取，提高拓頻的穩(wěn)定性和可靠性。同時(shí)，恢復(fù)的低頻部分相對(duì)常規(guī)反褶積方法有較大優(yōu)勢(shì)[22]，能拓寬相對(duì)頻寬和絕對(duì)頻寬，可減少子波的帶限影響。時(shí)間域反射系數(shù)反演目標(biāo)函數(shù)為

L(r,λ)=|x-r*b|2+λ|r|1,

(3)

式中：λ為正則化系數(shù)；|r|1為反射系數(shù)的L1范數(shù)。

將式(3)通過(guò)傅里葉變換轉(zhuǎn)換到頻率域，得到譜反演的一般式[21-22]：

L(r,λ)=|X-BFr|2+λ|r|1,

(4)

式中：X為地震道的傅里葉變換；F為傅里葉變換矩陣；B為子波的傅里葉變換矩陣形式。

約束反演問(wèn)題可以通過(guò)快速閾值迭代法(FIFTA)或梯度投影法(SPG)[23]求解。譜反演需要估算地震子波，在不考慮相位影響時(shí)，子波振幅譜可近似看成地震振幅譜的包絡(luò)[24-25]，類似地震振幅譜的平滑濾波，忽略白色反射系數(shù)的影響，估算子波的振幅譜West表示為

West≈B·SF(Rnw),

(5)

式中：SF為平滑濾波器符號(hào)；Rnw為有色反射系數(shù)的傅里葉變換。

由于估算子波本身是有色的，譜反演無(wú)法得到與井上規(guī)律一致的有色反射系數(shù)譜[2]，因而對(duì)薄層和夾層的分辨能力受到限制，需要結(jié)合井上反射系數(shù)頻譜的有色項(xiàng)進(jìn)一步改進(jìn)。

圖5 帶通子波譜藍(lán)化濾波結(jié)果Fig.5 Blueing results of band-pass wavelet

基于譜反演的譜藍(lán)化拓頻方法是先對(duì)原始地震資料進(jìn)行譜反演，以拓寬頻帶、壓縮子波的帶限影響；再利用井上反射系數(shù)擬合有色項(xiàng)，將藍(lán)譜算子應(yīng)用于譜反演后的拓頻資料。相對(duì)常規(guī)譜藍(lán)化處理，該方法能避免藍(lán)化子波的旁瓣振蕩及造成解釋假象；相對(duì)譜模擬反褶積，該方法不需要對(duì)子波做最小相位假設(shè)，拓頻的穩(wěn)定性和可靠性更好。

在實(shí)際資料處理中，譜反演受原始資料頻帶、信噪比等影響，在中高頻部分(100 Hz以上)出現(xiàn)陷波，需要對(duì)高頻做適當(dāng)截?cái)唷Ｍ瑫r(shí)，譜反演放大隨機(jī)噪聲，高頻截?cái)嗪笠m當(dāng)去噪，提高拓頻資料的信噪比，可采用主成分分析和曲波變換等自適應(yīng)的方法，減少對(duì)構(gòu)造細(xì)節(jié)的損傷。基于該方法的實(shí)際處理技術(shù)流程見圖6。

圖6 基于譜反演的譜藍(lán)化拓頻方法技術(shù)流程Fig.6 Technical process of spectrum blueing method based on spectrum inversion

2 模型測(cè)試

生成200個(gè)服從拉普拉斯分布的隨機(jī)序列，采樣率為2 ms，選取絕對(duì)值大于0的隨機(jī)序列作為儲(chǔ)層底面的反射系數(shù)，儲(chǔ)層的頂面反射系數(shù)與之符號(hào)相反，保證整個(gè)模型的反射系數(shù)振幅統(tǒng)計(jì)特性服從拉普拉斯振幅分布。隨機(jī)生成服從拉普拉斯概率分布的地層厚度，厚度變化范圍為1～10個(gè)樣點(diǎn)，能得到具有藍(lán)譜特征、與井上規(guī)律性一致的反射系數(shù)序列(見圖7(a))。模擬反射系數(shù)頻譜具有和井上接近的藍(lán)譜形態(tài)(見圖8(a))。

圖7 正演模型反射系數(shù)及合成地震記錄Fig.7 Reflection coefficient and seismic records of forward model

圖8 正演模型頻譜Fig.8 Seismic spectrum of forward model

圖9 正演模型譜反演懲罰參數(shù)迭代過(guò)程Fig.9 Iterative process of penalty parameter in spectrum inversion of forward model

使用主頻為35 Hz雷克子波與反射系數(shù)褶積后得到合成地震記錄(見圖7(b))，對(duì)應(yīng)的頻譜見圖8(b)。利用譜反演壓縮子波影響，采用梯度投影算法，sigma設(shè)置為10-1，懲罰參數(shù)λ經(jīng)過(guò)361次迭代后收斂，λ的迭代過(guò)程見圖9，反射系數(shù)的L1范數(shù)為0.756 7。高頻截?cái)嗪蟮淖V反演結(jié)果見圖7(c)，截?cái)囝l率為120 Hz，對(duì)應(yīng)的頻譜見圖8(c)，相對(duì)于原始地震記錄分辨率得到提高。

與模擬反射系數(shù)的頻譜對(duì)比，譜反演的地震頻譜中，低頻段雖具有明顯斜坡，但中高頻段的藍(lán)譜特征表現(xiàn)并不明顯，譜反演不能得到具有明顯藍(lán)譜特征的反射系數(shù)(見圖8(b-c))。繼續(xù)對(duì)譜反演結(jié)果進(jìn)行藍(lán)譜改造，藍(lán)譜改造后的合成記錄和相應(yīng)頻譜見圖7(d)和圖8(d)，頻譜的形態(tài)與模擬的反射系數(shù)譜基本保持一致，比原始合成記錄或譜反演結(jié)果分辨率更高，薄層和弱反射系數(shù)的響應(yīng)特征更清楚。

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

渤海萊州灣新近系發(fā)育淺水三角洲和曲流河沉積，砂體縱向表現(xiàn)為多期疊置、橫向儲(chǔ)層物性變化較快，受地震分辨率的制約，砂體的疊置關(guān)系并不明確，砂體描述中經(jīng)常出現(xiàn)似斷非斷、似連非連，地震多解性較強(qiáng)；已鉆井揭示的砂體厚度小于10 m，標(biāo)定時(shí)間厚度相對(duì)實(shí)鉆厚度偏大，給勘探和儲(chǔ)量評(píng)價(jià)帶來(lái)較大的影響。根據(jù)頂界面振幅屬性(見圖10)，KL構(gòu)造區(qū)V2油組呈多期河道疊置的特征，早期河道被晚期河道切割，河道的疊置關(guān)系復(fù)雜。

圖10 KL構(gòu)造區(qū)V2油組頂界面振幅屬性Fig.10 Amplitude attribute of top interface of V2 Oil Formation in KL Structrue Area

基于譜反演的譜藍(lán)化拓頻方法進(jìn)行高分辨率處理，V2油組河道砂體疊置關(guān)系得到較好解決。原始地震90°相移資料上，疊置砂體難以區(qū)分(見圖11(a))，儲(chǔ)層對(duì)應(yīng)負(fù)值，泥巖對(duì)應(yīng)正值；譜反演拓頻處理后，河道砂體的接觸關(guān)系改善，但似連非連(見圖11(b))；應(yīng)用該方法后砂體期次和疊置關(guān)系清楚(見圖11(c))，較好地規(guī)避似斷非斷、似連非連的解釋風(fēng)險(xiǎn)，有助于厘清相鄰期次河道砂體的發(fā)育規(guī)律和展布特征。

圖11 V2油組處理前后90°相移剖面Fig.11 90 ° phase shift profiles before and after V2 Oil Formation treatment

萊州灣KL構(gòu)造區(qū)在V3油組有利位置鉆探KL-3、KL-4和KL-5井，鉆遇油層厚度分別為10.3、5.2和11.0 m。KL地區(qū)老資料標(biāo)定后的時(shí)間厚度和實(shí)鉆厚度存在較大的差異(見圖12(a))。新資料井震相關(guān)關(guān)系得到改善，3口井平均相關(guān)因數(shù)從原來(lái)的0.83提高到0.92(歸一化后)(見表2)。KL-3井儲(chǔ)層厚度與振幅響應(yīng)更符合實(shí)際(見圖12(b))，為KL構(gòu)造區(qū)未來(lái)整體評(píng)價(jià)及儲(chǔ)量計(jì)算提供較好的基礎(chǔ)資料。

表2 V3油組井震平均相關(guān)因數(shù)

圖12 V3油組處理前后90°相移剖面Fig.12 90 ° phase shift profiles before and after V3 Oil Formation treatment

4 結(jié)論

(1)反射系數(shù)的藍(lán)譜特征主要受時(shí)間序列的關(guān)聯(lián)性及薄層為主的地層結(jié)構(gòu)共同影響，與萊州灣KL地區(qū)的地質(zhì)沉積認(rèn)識(shí)保持一致，是譜藍(lán)化應(yīng)用的前提和基礎(chǔ)。

(2)相比于單獨(dú)的譜反演或譜藍(lán)化方法，基于譜反演的譜藍(lán)化拓頻方能引入井上先驗(yàn)信息，突出薄層和弱反射系數(shù)的地震響應(yīng)，有效規(guī)避帶限資料形成的解釋缺陷。

(3)譜反演放大高頻隨機(jī)噪聲，應(yīng)用藍(lán)化算子前應(yīng)做截?cái)嗪瓦m當(dāng)去噪；基于譜反演的譜藍(lán)化拓頻方法依賴于褶積的穩(wěn)態(tài)模型，進(jìn)行非穩(wěn)態(tài)的多道聯(lián)合譜反演是下一步研究方向。