張顯文,范廷恩,胡光義,高云峰,王宗俊

(中海油研究總院,北京100027)

地下介質(zhì)對地震波能量的吸收作用是地層的內(nèi)在屬性,地層吸收特性與巖性、孔隙度、孔隙流體等密切相關(guān)。尤其是地層含油氣后,地震波頻率成分會異常衰減,表現(xiàn)在低頻能量相對增加,高頻能量相對減弱。因此,地震頻譜衰減分析可以作為油氣檢測的重要技術(shù)手段[1],如鶯歌海盆地、瓊東南盆地和遼東灣等海區(qū)的天然氣藏預(yù)測技術(shù)都將地層吸收特性作為識別含氣信息的有效標(biāo)志之一[2]。

Mitchell等(1996)[3]提出了地震信號能量吸收分析的方法,分時窗對地震記錄進(jìn)行頻譜分析,通過高頻段和低頻段頻率組分衰減速率的異常變化進(jìn)行油氣預(yù)測。為了減弱反射系數(shù)對吸收分析的影響,Lichman等(2003)[4]提出了瞬時子波吸收分析油氣預(yù)測技術(shù),通過復(fù)賽譜技術(shù)提取時窗內(nèi)的瞬時子波,根據(jù)瞬時子波振幅譜高頻端能量的變化來評估儲層的含油氣性,提高了油氣預(yù)測的精度。董寧等(2008)[5]通過小波變換得到地震信號的振幅譜,根據(jù)地震波能量及對應(yīng)的頻率擬合出能量與頻率的衰減梯度,求取振幅衰減梯度因子進(jìn)行吸收分析。Hu等(2009)[6]通過雙相介質(zhì)正演模擬,探索了含油氣儲層對地震波的吸收衰減作用。Jiang等(2010)[7]采用小波變換技術(shù),利用能量吸收分析方法進(jìn)行了吸收分析油氣預(yù)測。刁瑞等(2011)[8]采用基于廣義S變換的吸收衰減分析技術(shù)對理論模型和實(shí)際數(shù)據(jù)進(jìn)行油氣識別,取得了較好的預(yù)測效果。

常規(guī)吸收分析技術(shù)都是基于常時窗長度計算吸收因子,對于強(qiáng)振幅型油氣藏具有較好的應(yīng)用效果,但對于中等強(qiáng)度振幅型油氣藏存在一定的多解性。由于含油氣地層在地震記錄上通常表現(xiàn)為低頻的地震反射特征,我們將局部頻率算法引入到吸收分析技術(shù)中,提出了基于自適應(yīng)時窗的瞬時子波吸收分析技術(shù)(AWEA)。首先闡明AWEA技術(shù)原理及其實(shí)現(xiàn)步驟;然后用含氣砂體模型測試分析AWEA算法的抗噪性;進(jìn)一步通過某氣藏實(shí)際地震資料的應(yīng)用試驗,對比分析AWEA技術(shù)與常規(guī)吸收分析技術(shù)的油氣預(yù)測效果及實(shí)用性。

1 技術(shù)原理與實(shí)現(xiàn)步驟

AWEA技術(shù)基于疊后純波地震數(shù)據(jù),算法實(shí)現(xiàn)逐道進(jìn)行,對于每一道地震記錄,首先利用局部頻率算法計算每一采樣點(diǎn)的自適應(yīng)時窗,再基于復(fù)賽譜技術(shù)提取自適應(yīng)時窗內(nèi)地震記錄的瞬時子波頻譜,最后通過對頻譜的高頻段和低頻段進(jìn)行指數(shù)擬合求取吸收因子。

1.1 自適應(yīng)時窗的計算

常規(guī)吸收分析技術(shù)在計算吸收因子時,每個采樣點(diǎn)的時窗長度相同,而沒有考慮含油氣儲層與非儲層地震反射的頻率特征差異。由于地層含油氣后具有相對低頻的地震反射特征,其瞬時波長相對較長,因此,吸收分析時窗的選取應(yīng)與地震記錄的頻率特性相關(guān)。

在計算吸收系數(shù)之前,首先計算地震記錄S(t)的瞬時頻率ω(t)[9]:

(1)

式中:H(t)為地震記錄S(t)的希爾伯特變換;H′(t)和S′(t)分別是H(t)和S(t)的導(dǎo)數(shù);τ為中間變量。對瞬時頻率取倒數(shù)得到瞬時波長,再乘以常數(shù)C得到自適應(yīng)時窗的長度ATW:

(2)

上述常規(guī)求取瞬時頻率的算法容易受到噪聲的影響,會產(chǎn)生無意義的負(fù)值[10]?？紤]到通過正則化反演計算地震記錄瞬時頻率的算法具有很好的穩(wěn)定性和抗噪性,因此,本文采用基于正則化反演實(shí)現(xiàn)的局部頻率算法[11]:

(3)

式中:n代表(1)式中的分子;D為由(1)式中的分母構(gòu)成的對角算子;S為平滑整形算子,通常選取高斯時窗;λ為D的模的平方,控制反演迭代的快速收斂;I為恒等算子。

將(3)式計算求得的ω(t)代入(2)式,得到相應(yīng)的自適應(yīng)時窗。此時,每個采樣點(diǎn)的自適應(yīng)時窗的計算與地震信號的頻率特征建立了聯(lián)系。

1.2 瞬時子波譜分解

通過上述步驟1.1得到每個時間采樣點(diǎn)的自適應(yīng)時窗后,進(jìn)一步可以通過時窗內(nèi)地震信號的頻譜分析來求取吸收因子,但該種算法受反射系數(shù)的影響較大,難以區(qū)分強(qiáng)振幅異常是由地層含氣引起還是由高阻抗引起[2]。Lichman等(2003)[4]提出的瞬時子波吸收分析技術(shù)利用復(fù)賽譜技術(shù)提取時窗內(nèi)地震記錄的瞬時子波,通過對其頻譜分析來計算吸收因子,可以減弱反射系數(shù)對吸收分析的影響,提高吸收分析油氣預(yù)測的精度。

對于每道地震記錄的采樣點(diǎn)τi,利用傅里葉變換計算自適應(yīng)時窗內(nèi)地震記錄的頻譜S(τi,f),作為該樣點(diǎn)的瞬時頻譜,記為IFS(f)。進(jìn)一步,利用復(fù)賽譜技術(shù)計算該采樣點(diǎn)的瞬時子波:

(4)

在時間域,地震記錄為子波和反射系數(shù)的褶積;在復(fù)賽譜域,子波在原點(diǎn)附近,而反射系數(shù)和隨機(jī)噪聲在遠(yuǎn)離原點(diǎn)的邊緣附近。因此,在復(fù)賽譜域可以通過“低通濾波(Low Filter)”的方式分離子波與反射系數(shù)。所以，復(fù)賽譜參數(shù)的設(shè)置與帶通濾波參數(shù)相似,用參數(shù)C1,C2,C3,C4控制,如圖1所示。

圖1 復(fù)賽譜域低通濾波示意圖解

對于最小相位子波,位于復(fù)賽譜域正半軸區(qū)域,取C1=C2=0;對于最大相位子波,位于復(fù)賽譜域負(fù)半軸區(qū)域,取C3=C4=0;通常選取C1=C2=-0.08s,C3=C4=1.0s。一般情況下子波和反射系數(shù)的分離區(qū)較大,因此子波提取結(jié)果對復(fù)賽譜參數(shù)的取值不敏感。由此瞬時子波振幅譜為

(5)

式中：LLowFilter代表低通濾波算子。由于傅里葉變換的算法穩(wěn)定性受時窗內(nèi)樣點(diǎn)數(shù)限制,因此,在譜分解前要先對地震記錄進(jìn)行時域和頻域插值。首先,為了保證地震記錄內(nèi)插前、后頻率組分不變,通過頻率域添零達(dá)到期望的奈奎斯特頻率實(shí)現(xiàn)時域內(nèi)插。類似地,為了使得有效頻帶范圍具有足夠的頻率樣點(diǎn)數(shù),通過時間域零填充實(shí)現(xiàn)頻域內(nèi)插,進(jìn)而實(shí)現(xiàn)穩(wěn)定的瞬時譜分解流程。

1.3 非線性反演衰減分析

對于每個時間采樣點(diǎn),按照步驟1.1和1.2得到自適應(yīng)時窗內(nèi)的瞬時子波頻譜后,便可以通過頻譜分析來計算其吸收因子。按照Lichman等(2003)[4]提出的理論,巖石孔隙內(nèi)的流體主要引起子波低頻段能量變化,而氣體對子波的高頻吸收起主要影響作用。因此,對于得到的瞬時子波振幅譜IWFS(f)與地層吸收的關(guān)系(圖2),滿足如下方程:

(6)

式中:ω為角頻率;ω0為瞬時子波主頻;Qf為低頻吸收因子,對孔隙流體敏感;Qg為高頻吸收因子,

圖2 瞬時子波頻譜與地層吸收關(guān)系

對孔隙內(nèi)存在氣體敏感;c1,c2,c3為常數(shù)。

根據(jù)方程(6),對于每個采樣點(diǎn)瞬時子波振幅譜的吸收因子計算是通過求解(7)式的最小化問題:

(7)

式中:ω0為有效信號的高截頻。利用非線性反演算法求解方程(7)的最小化問題[12],得到每個時間采樣點(diǎn)的吸收因子。

1.4 AWEA技術(shù)實(shí)現(xiàn)步驟

AWEA算法的實(shí)現(xiàn)分為3步:

1) 基于正則化反演算法計算目標(biāo)儲層地震記錄的局部頻率,進(jìn)一步計算自適應(yīng)時窗;

2) 對于每個樣點(diǎn),利用復(fù)賽譜技術(shù)計算時窗內(nèi)的瞬時子波頻譜;

3) 對頻譜高頻段和低頻段進(jìn)行非線性反演,計算高、低頻吸收因子。

2 測試分析與應(yīng)用試驗

由于吸收因子的求取是基于地震波頻率特性,而地震信號的頻譜易受噪聲的影響,因此,我們首先通過含噪砂體模型進(jìn)行算法的抗噪性分析,再通過某氣藏實(shí)際地震資料的應(yīng)用進(jìn)行吸收分析油氣預(yù)測,并與常規(guī)吸收分析技術(shù)進(jìn)行效果對比。

2.1 算法抗噪性分析

為了分析算法的抗噪性,構(gòu)建如圖3所示的含氣砂體模型。模型中間一層的砂體厚度為30m(CDP范圍為1～80),在兩邊的致密砂巖(速度為3200m/s)中間夾著含氣砂巖(CDP范圍為26～55),Q值為25;上、下圍巖的泥巖速度為2400m/s。

圖3 簡單砂體模型(CDP:1～80)

圖4給出了含氣砂體模型的無噪聲合成記錄(圖4a;CDP:26～55;SNR=∞)和AWEA技術(shù)計算的高頻吸收Qg剖面(圖4b),可見Qg剖面對含氣砂巖有較好的吸收異常刻畫。圖5是添加了少許隨機(jī)噪聲后的合成記錄(圖5a;SNR=20)和AWEA技術(shù)計算的Qg剖面,可見少許的隨機(jī)噪聲對含氣砂巖的吸收特征影響很小。隨著隨機(jī)噪聲的進(jìn)一步增加(圖6a;SNR=10),相應(yīng)Qg剖面(圖6b) 上含氣砂巖內(nèi)部的吸收特征受到一定的影響,但是對含氣砂巖整體吸收特征的刻畫影響較小。但當(dāng)隨機(jī)噪聲較強(qiáng)時(圖7a;SNR=5),所獲Qg剖面(圖7b)上含氣砂巖內(nèi)部和砂體邊界的吸收特征都受到了較嚴(yán)重的影響,此時利用吸收特征進(jìn)行油氣預(yù)測需要考慮噪聲的因素。模型測試分析結(jié)

果表明,AWEA算法具有較好的抗噪性,但當(dāng)隨機(jī)噪聲較強(qiáng)時,算法應(yīng)用需要重視前期的地震資料去噪處理。

2.2 應(yīng)用試驗效果分析

圖8給出了對某氣藏實(shí)際地震資料應(yīng)用常規(guī)吸收分析技術(shù)和AWEA技術(shù)進(jìn)行吸收分析油氣預(yù)測的結(jié)果。其中圖8a為過P-2井和P-3井的地震連井剖面,包括Ⅰ,Ⅱ,Ⅲ,Ⅳ4個層位,各井鉆遇層位的流體情況如圖中標(biāo)記所示。圖8b是用常規(guī)吸收分析技術(shù)計算的高頻吸收Qg剖面。分析圖8b可以看出:兩口井鉆遇的層位Ⅲ處,高頻吸收Qg剖面上表現(xiàn)出了弱吸收異常值,這也與它們是干層的實(shí)際情況吻合;對于P-3井,層位Ⅰ和Ⅱ處表現(xiàn)為強(qiáng)吸收特征,這與實(shí)際這兩個層位含氣相一致,而Ⅳ層位處吸收異常值較弱,這與該層位實(shí)際為氣層所應(yīng)具有的強(qiáng)吸收特征不符;對于P-2井,層位Ⅰ和Ⅱ處鉆遇水層,但由于地震反射能量較強(qiáng),吸收分析給出了錯誤的強(qiáng)吸收特征,而含氣的Ⅳ層位由于地震反射能量較弱,氣層所應(yīng)該具有的強(qiáng)吸收特征沒有表征出來。因此,常規(guī)吸收分析技術(shù)的油氣預(yù)測效果受地震反射能量強(qiáng)弱的影響,存在較強(qiáng)的多解性。

圖4 砂體模型合成記錄(CDP:26～55;SNR=∞)(a)和AWEA高頻吸收Qg剖面(b)

圖5 砂體模型加噪合成記錄(CDP:26～55;SNR=20)(a)和AWEA高頻吸收Qg剖面(b)

圖6 砂體模型加噪合成記錄(CDP:26～55;SNR=10)(a)和AWEA高頻吸收Qg剖面(b)

圖7 砂體模型加噪合成記錄(CDP:26～55;SNR=5)(a)和AWEA高頻吸收Qg剖面(b)

圖8 某氣藏實(shí)際地震資料吸收分析技術(shù)油氣預(yù)測應(yīng)用實(shí)例a 過P-2井和P-3井的地震連井剖面; b 常規(guī)吸收分析技術(shù)計算的高頻吸收Qg剖面; c AWEA技術(shù)計算的高頻吸收Qg剖面

分析應(yīng)用本文提出的AWEA技術(shù)計算的高頻吸收Qg剖面(圖8c):強(qiáng)吸收異常值有效地預(yù)測了氣層,這與實(shí)際井鉆遇的氣層相吻合;對于常規(guī)吸收分析技術(shù)存在多解性的層位,如P-2井鉆遇的Ⅰ,Ⅱ,Ⅳ層位,圖8c中亦給出了較好的吸收異常的刻畫。由此可見,AWEA技術(shù)引入基于局部波長計算的自適應(yīng)時窗,減弱了非油氣強(qiáng)振幅響應(yīng)對吸收特征分析結(jié)果的影響,降低了吸收分析油氣預(yù)測的多解性。

為了進(jìn)一步說明自適應(yīng)時窗技術(shù)對不同頻率地震響應(yīng)的反應(yīng)能力,圖9展示了獲取圖8c高頻吸收Qg剖面時自適應(yīng)時窗的計算過程。圖9a為基于地震記錄計算的局部頻率剖面,圖9b是進(jìn)而求取的局部波長剖面,設(shè)置了0.75倍半局部波長為自適應(yīng)時窗長度。由于氣層具有低頻的地震反射特征,因此時窗長度較長;而非氣儲層時窗長度相對較短。因此,通過自適應(yīng)時窗的計算對吸收分析間接地起到了增強(qiáng)氣層吸收特征的作用。

圖9 自適應(yīng)時窗計算a 局部頻率剖面; b 局部波長剖面

3 結(jié)束語

常規(guī)吸收分析技術(shù)基于常時窗計算吸收因子,吸收分析油氣預(yù)測結(jié)果易受地震反射強(qiáng)弱的影響,盡管瞬時子波吸收分析技術(shù)能夠減弱反射系數(shù)的影響,但在實(shí)際應(yīng)用中這些技術(shù)對非強(qiáng)振幅型油氣藏的油氣預(yù)測結(jié)果仍然存在較強(qiáng)的多解性。本文提出的AWEA油氣預(yù)測技術(shù)通過地震記錄的頻率特性計算吸收分析自適應(yīng)時窗,使得吸收因子的計算充分利用了油氣的低頻地震反射特征,間接地起到了增強(qiáng)氣層吸收特征的作用,尤其是針對非強(qiáng)振幅型油氣藏,油氣預(yù)測的多解性問題將會得到較好的改善。

自適應(yīng)時窗技術(shù)同樣適合于其它時頻分析技術(shù),如短時傅里葉變換、小波變換、S變換等技術(shù)。當(dāng)然,該技術(shù)同時也適合于其它吸收分析油氣預(yù)測技術(shù),如基于地震記錄的吸收分析油氣預(yù)測技術(shù),但瞬時子波吸收分析油氣預(yù)測技術(shù)具有相對更好的油氣預(yù)測效果。

致謝:感謝中科院地質(zhì)與地球物理所姚振興和中海油研究總院李維新、李緒宣、趙偉對論文研究內(nèi)容給予的指導(dǎo)和幫助。

參 考 文 獻(xiàn)

[1] 郭棟,韓文功,王興謀.地層吸收特征參數(shù)分析預(yù)測深層氣藏方法[J].油氣地球物理,2008,6(3):20-24

Guo D,Han W G,Wang X M.Predicting deep gas pool by formation absorption characteristic parameter[J].Petroleum Geophysics,2008,6(3):20-24

[2] 張益明,李緒宣,符力耘.瞬時子波吸收分析技術(shù)的應(yīng)用——以中國近海天然氣檢測為例[J].石油物探,2009,48(4):368-376

Zhang Y M,Li X X,Fu L Y.Application of instantaneous wavelet absorption analysis technology:case study on offshore gas detection in China[J].Geophysical Prospecting for Petroleum,2009,48(4):368-376

[3] Mitchell J T,Derzhi N,Lichman E,et al.Energy absorption analysis:a case study[J].Expanded Abstracts of 66thAnnual Internat SEG Mtg,1996,1785-1788

[4] Lichman E,Goloshubin G.Unified approach to gas and fluid detection on instantaneous seismic wavelets[J].Expanded Abstracts of 73rdAnnual Internat SEG Mtg,2003,1699-1702

[5] 董寧,楊立強(qiáng).基于小波變換的吸收衰減技術(shù)在塔河油田儲層預(yù)測中的應(yīng)用研究[J].地球物理學(xué)進(jìn)展,2008,23(2):533-538

Dong N,Yang L Q.Application of absorption and attenuation based on wavelet transform for prediction of reservoir in Tahe Oilfield[J].Progress in Geophysics,2008,23(2):533-538

[6] Hu X P,Liang Y,Tao C F,et al.Application of hydrocarbon detection technology based on double phase medium theory[J].Expanded Abstracts of 79thAnnual Internat SEG Mtg,2009,593-596

[7] Jiang L,Yang S G,Zhao M J,et al.Frequency-dependent seismic attributes and their application for hydrocarbon detection[J].Expanded Abstracts of 80thAnnual Internat SEG Mtg,2010,1566-1569

[8] 刁瑞,李振春,韓文功,等.基于廣義S變換的吸收衰減分析技術(shù)在油氣識別中的應(yīng)用[J].石油物探,2011,50(3):260-265

Diao R,Li Z C,Han W G,et al.Application of absorption and attenuation analysis technique based on generalized S transform for hydrocarbon identification[J].Geophysical Prospecting for Petroleum,2011,50(3): 260-265

[9] Taner M T,Koehler F,Sheriff R E.Complex seismic trace analysis[J].Geophysics,1979,44(6):1041-1063

[10] White R E.Properties of instantaneous seismic attributes[J].The Leading Edge,1991,56(10):26-32

[11] Sergey F.Local seismic attributes[J].Geophysics,2007,73(3):A29-A33

[12] Byrd R H,Mary E,Jorge N.An interior point algorithm for large-scale nonlinear programming[J].SIAM Journal on Optimization,1999,9(4):877-900

(編輯：陳 杰)