刁 瑞

(中國(guó)石化勝利油田分公司物探研究院，山東 東營(yíng) 257022)

0 引言

地震波在地下傳播過(guò)程中受到地層吸收衰減的影響，導(dǎo)致地震波高頻能量損失和相位畸變，降低了地震資料的分辨率[1]。對(duì)于薄互層油氣藏和復(fù)雜隱蔽巖性油氣藏，常規(guī)的地震反射波資料難以滿足薄砂體儲(chǔ)層精細(xì)刻畫的要求，特別是中深層油氣藏面臨著難以精細(xì)預(yù)測(cè)薄層砂體發(fā)育情況的難題[2]，需要針對(duì)性地開(kāi)展地震資料提高分辨率處理。提高分辨率是地震資料處理的重要環(huán)節(jié)之一，提高分辨率處理方法包括：預(yù)測(cè)反褶積、分頻帶反褶積、零相位反褶積、子波整形反褶積、反Q濾波、譜模擬反褶積、譜白化、井約束提頻、近地表Q值吸收補(bǔ)償?shù)确椒╗3-8]。反褶積方法種類很多，主要區(qū)別之一在于對(duì)地震子波的假設(shè)和估計(jì)地震子波的方法，通常假設(shè)地震信號(hào)是平穩(wěn)信號(hào)，反射系數(shù)序列是白噪序列，子波是最小相位[9]。各種提高地震分辨率處理方法在實(shí)際地震資料處理中發(fā)揮了舉足輕重的作用，并取得了良好的應(yīng)用效果。

隨著勘探開(kāi)發(fā)程度的不斷提高，成熟油區(qū)薄互層油氣藏、隱蔽巖性油氣藏逐漸成為勘探開(kāi)發(fā)的重點(diǎn)，薄砂體儲(chǔ)層對(duì)提高地震分辨率處理方法提出了更高的要求。針對(duì)不同品質(zhì)的地震數(shù)據(jù)，應(yīng)該采用何種提高分辨處理方法和參數(shù)，以及地震數(shù)據(jù)的分辨率應(yīng)該提高到何種程度等問(wèn)題，制約著提高分辨率處理效果。劉浩杰等[10]從客觀量化角度對(duì)地震資料分辨率與頻譜特征參數(shù)的關(guān)系進(jìn)行了研究。李曙光等[11]從剖面和頻譜特征對(duì)3種提高分辨率方法進(jìn)行了探討。芮擁軍[12]分析了相對(duì)保幅處理的難點(diǎn)，對(duì)振幅處理、噪聲壓制和提高分辨率等處理技術(shù)的保幅性進(jìn)行了分析。錢榮鈞[13]分析了分辨率的分類和定義，并給出了計(jì)算空間分辨率的方法。萬(wàn)歡等[14]通過(guò)頻譜、信噪比和相似系數(shù)等參數(shù)，對(duì)比分析了4種提高分辨率技術(shù)，在疊前地震資料中取得了較好的應(yīng)用效果。

常規(guī)地震資料提高分辨率處理過(guò)程中，由于缺乏系統(tǒng)的提高分辨率效果定量評(píng)價(jià)分析方法和軟件，需要對(duì)提高分辨率處理方法和處理參數(shù)等進(jìn)行大量的重復(fù)試驗(yàn)，在進(jìn)行提高分辨率處理結(jié)果對(duì)比評(píng)價(jià)中，人為因素影響較大，難以真實(shí)客觀地定量評(píng)價(jià)處理方法和參數(shù)的選擇。筆者針對(duì)薄砂體精細(xì)描述的地質(zhì)需求，開(kāi)展了提高地震分辨率效果定量評(píng)價(jià)方法研究，并研發(fā)了軟件系統(tǒng)，從而高效優(yōu)選提高分辨率處理技術(shù)、流程和參數(shù)，全面可靠的監(jiān)控提高分辨率處理流程，實(shí)現(xiàn)地震數(shù)據(jù)提高分辨率處理的效果定量評(píng)價(jià)分析。效果定量評(píng)價(jià)技術(shù)在新疆車排子地區(qū)沙灣組薄層砂體油藏進(jìn)行了應(yīng)用，對(duì)提高分辨率處理進(jìn)行全流程的效果定量評(píng)價(jià)分析，提高分辨率處理后識(shí)別薄層能力顯著增強(qiáng)，為后續(xù)的油藏描述和井位部署奠定了堅(jiān)實(shí)的基礎(chǔ)。

1 方法原理

如何實(shí)現(xiàn)提高分辨率不同技術(shù)、流程及參數(shù)的最優(yōu)化選擇，一直是提高分辨率處理所面臨的難題。在提高分辨率處理效果評(píng)價(jià)方面，針對(duì)常規(guī)頻譜和剖面對(duì)比的優(yōu)選難題，提出了對(duì)不同數(shù)據(jù)進(jìn)行多種屬性計(jì)算、多種屬性組合優(yōu)選的效果定量評(píng)價(jià)方法。

首先，選用不同提高分辨率方法和參數(shù)進(jìn)行數(shù)據(jù)處理；然后，通過(guò)數(shù)學(xué)公式計(jì)算測(cè)試數(shù)據(jù)地質(zhì)目的層屬性參數(shù)，包括：振幅、頻率、相位、信噪比、時(shí)頻特征、相干特征、振幅比特征、波形特征、斷層特征等；再次，綜合分析所有屬性信息，對(duì)提高分辨率處理結(jié)果進(jìn)行效果定量評(píng)價(jià)；最后，從不同測(cè)試結(jié)果中優(yōu)選最理想效果，確定最佳提高分辨率方法和處理參數(shù)。通過(guò)技術(shù)(參數(shù))差異性繪圖，結(jié)合屬性提取、合成記錄標(biāo)定等多因素綜合分析，從而實(shí)現(xiàn)對(duì)提高分辨率處理前后效果優(yōu)劣的評(píng)價(jià)，達(dá)到高效勘探與開(kāi)發(fā)的目的。

1.1 頻譜類方法

地震波在地層中傳播，隨著傳播距離的增加，高頻成份相對(duì)于低頻而言更容易衰減，包含不同頻率成份的頻譜則能夠反映地震數(shù)據(jù)的分辨率。振幅譜主要通過(guò)主頻和頻寬兩個(gè)參數(shù)來(lái)描述，是一個(gè)時(shí)間域信號(hào)在頻率域的表達(dá)方式，可以對(duì)時(shí)間域信號(hào)進(jìn)行快速傅里葉變換得到。

通過(guò)傅里葉變換算法可以求得地震數(shù)據(jù)的振幅譜和分貝譜，從而計(jì)算主頻、有效頻寬、倍頻程、低截頻值、高截頻值、不同頻段能量等信息，分析不同提高分辨率方法和參數(shù)的頻譜屬性，從而實(shí)現(xiàn)最優(yōu)參數(shù)的優(yōu)選。振幅譜和分貝譜計(jì)算公式如下：

dB(f)=20log10|H(f)|

(1)

式中，dB(f)是分貝譜，|H(f)|是歸一化后的振幅譜，H(ω)是振幅譜，h(t)是地震記錄。

對(duì)不同提高分辨率方法和參數(shù)的地震數(shù)據(jù)分別進(jìn)行傅里葉變換，在變換得到的頻譜中求取主頻、頻寬等屬性，優(yōu)勢(shì)頻寬以-18 dB為標(biāo)準(zhǔn)，有效頻寬以-24 dB為標(biāo)準(zhǔn)。主頻越高，有效頻帶越寬，倍頻程越大，則地震分辨率越高。隨著“兩寬一高”地震技術(shù)的發(fā)展，越來(lái)越重視低頻段信息的恢復(fù)和保護(hù)，低截頻值越小，則表示低頻段信息越豐富，倍頻程也越大，從而分辨率越高。

1.2 子波類方法

地震子波提取方法主要包括：自相關(guān)方法、同態(tài)譜方法、譜模擬方法、高階統(tǒng)計(jì)量方法、高階譜方法、時(shí)空變方法和非線性理論方法等，地震子波提取方法各有優(yōu)缺點(diǎn)，適用于不同信噪比的地震數(shù)據(jù)。

通過(guò)地震記錄的高階累積量可以近似得到地震子波的高階累積量，根據(jù)地震記錄的功率譜計(jì)算得到子波振幅譜，利用地震記錄的雙譜重構(gòu)地震子波的相位。地震子波頻譜的計(jì)算公式如下[15]：

B(ω)=|B(ω)|ejφb(ω)。

(2)

地震子波的相位譜與雙譜之間的關(guān)系為：

ψ3b(ω1,ω2)=φb(ω1)+φb(ω2)-φb(ω1+ω2),

(3)

式中：ω為頻率；φb(ω)為地震子波相位譜；ψ3b(ω1,ω2)為根據(jù)地震信號(hào)三階累積量加時(shí)窗計(jì)算得到的雙譜。

在地質(zhì)目標(biāo)層位劃取時(shí)窗，提取地震子波后，一方面，可以計(jì)算地震子波的主頻、頻寬、相位等，另一方面，對(duì)提高分辨率處理前、后的子波進(jìn)行互相關(guān)，可以計(jì)算得到不同方法或參數(shù)所對(duì)應(yīng)的子波一致性定量數(shù)值，從而直觀評(píng)價(jià)提高分辨率的效果。

在子波頻譜評(píng)價(jià)方面，通過(guò)變換得到振幅譜和相位譜后，可以根據(jù)振幅譜計(jì)算得到主頻、頻寬、低截頻值、高截頻值、倍頻程等，與頻譜類方法的評(píng)價(jià)標(biāo)準(zhǔn)一致。另外，根據(jù)相位譜評(píng)價(jià)提高分辨率處理對(duì)相位的改變程度。

在子波一致性評(píng)價(jià)方面，設(shè)定子波一致性評(píng)價(jià)標(biāo)準(zhǔn)，低于該標(biāo)準(zhǔn),則表明提高分辨率處理破壞了子波一致性;高于該標(biāo)準(zhǔn)，則表明提高分辨率處理在壓縮地震子波的同時(shí)，保持了較好的子波一致性。

1.3 信噪比類方法

信噪比分析方法包括：能量疊加法、頻譜估計(jì)法、功率譜估算法、相關(guān)法和特征值法等[16]。信噪比低的地震資料優(yōu)選頻譜或功率譜法。信噪比是衡量地震資料品質(zhì)的重要參數(shù)，通過(guò)分析提高分辨率處理前后數(shù)據(jù)的信噪比，可以評(píng)價(jià)不同提高分辨率方法和參數(shù)的處理效果，從而根據(jù)信噪比分析結(jié)果優(yōu)選最佳參數(shù)。

通過(guò)信噪比分析可以避免過(guò)多高頻噪聲的不利影響，從而確保擴(kuò)展頻帶寬度的可信度。由于不同層位或同一層位的地震波能量差異巨大，通常情況下在某一時(shí)窗內(nèi)進(jìn)行信噪比分析，能量疊加法信噪比分析公式如下：

式中:M為分析時(shí)窗的采樣點(diǎn)數(shù),N為地震道數(shù),dij是地震記錄振幅。

選取多種信噪比計(jì)算方法，進(jìn)行定量化的信噪比分析，可以對(duì)提高分辨率結(jié)果進(jìn)行效果評(píng)價(jià)。在多種方法達(dá)到相同提高分辨率效果時(shí)，通過(guò)信噪比高低進(jìn)行評(píng)價(jià)，從而指導(dǎo)方法或參數(shù)的優(yōu)選。

1.4 時(shí)頻分析方法

時(shí)頻分析方法解決了一維傅里葉變換中時(shí)間域和頻率域分離的不足，可以同時(shí)分析不同時(shí)刻包含的各頻率成份，以及不同頻率成份隨傳播時(shí)間的變化特征[17]。

修正S變換的表達(dá)式如下：

exp[-fp(t-τ)2/2λ]exp(-i2πft)dt,

(5)

式中:GSTN(τ,f)是二維時(shí)頻譜，h(t)是地震記錄，λ、p是修正S變換調(diào)節(jié)參數(shù)，可以獲得不同分辨率的時(shí)頻分析結(jié)果。

對(duì)比提高分辨率處理前、后的時(shí)頻譜特征，一方面，可以分析不同時(shí)刻的頻譜特征，避免了傅里葉變換時(shí)窗大小的限制，能夠更加精細(xì)地分析頻率變化特征，另一方面，可以分析不同頻率成份隨時(shí)間的變化情況，即不同頻段能量的時(shí)間變化曲線，計(jì)算公式如下：

(6)

式中:t為時(shí)間，單位ms；f為頻率，f1和f2為頻段的起始頻率和截止頻率，單位Hz；GSTN(t,f)為地震記錄的二維時(shí)頻譜。提高分辨率處理后，高頻段能量曲線的數(shù)值明顯增大，表明時(shí)頻譜中優(yōu)勢(shì)頻帶和有效頻帶寬度增加，而高頻段能量曲線的變化趨勢(shì)保持不變，表明地震記錄能量得到很好的恢復(fù)和補(bǔ)償，并且時(shí)頻譜的淺中深層能量關(guān)系匹配較好，保持了原始數(shù)據(jù)的頻譜特征和時(shí)頻特性。

2 定量評(píng)價(jià)軟件系統(tǒng)研發(fā)

目前常用的編程語(yǔ)言主要有C/C++、Fortran、Matlab等語(yǔ)言[18-19]，在此采用C++和Qt編程語(yǔ)言，進(jìn)行提高地震分辨率效果定量評(píng)價(jià)算法開(kāi)發(fā)與軟件研發(fā)，所研發(fā)的軟件具有可移植、跨平臺(tái)的特點(diǎn)。該軟件共包括：頻率分析、波形分析、子波一致性分析、能量分析、信噪比分析、振幅譜分析、時(shí)頻譜分析、子波分析、沿層屬性分析等模塊功能。

對(duì)原始地震記錄進(jìn)行提高分辨率處理，得到不同方法和參數(shù)的一系列結(jié)果，可以對(duì)單個(gè)數(shù)據(jù)進(jìn)行頻譜分析、時(shí)頻分析和子波分析等，也可以將一系列提高分辨率結(jié)果進(jìn)行編號(hào)，計(jì)算屬性特征的變化趨勢(shì)，從而根據(jù)柱狀圖進(jìn)行效果定量評(píng)價(jià)。

對(duì)單個(gè)提高分辨率數(shù)據(jù)進(jìn)行評(píng)價(jià)分析，圖1是地震數(shù)據(jù)剖面顯示和二維時(shí)頻譜分析結(jié)果，在地震剖面上劃取時(shí)窗可以快速實(shí)現(xiàn)頻譜分析、時(shí)頻分析和子波分析等功能。

圖1 地震數(shù)據(jù)顯示及時(shí)頻譜分析Fig.1 Seismic section and time-frequency spectrum

將一系列提高分辨率結(jié)果進(jìn)行編號(hào)，分析不同方法或參數(shù)的屬性變化特征。圖2是頻譜類方法定量評(píng)價(jià)結(jié)果，頻譜類屬性包括：主頻、頻寬、低截頻值、高截頻值、低頻段能量、高頻段能量和倍頻程等。對(duì)不同提高分辨率方法和處理參數(shù)的測(cè)試結(jié)果進(jìn)行編號(hào)，圖中共計(jì)6個(gè)處理參數(shù)。主頻屬性分析結(jié)果顯示，不同處理參數(shù)的主頻基本相同，編號(hào)4～6的主頻略高，約為20 Hz；從頻寬屬性、高截頻值屬性、高頻段能量屬性的分析結(jié)果可以看出，編號(hào)4和5的處理參數(shù)最佳，頻寬約為44 Hz，高截頻值為50 Hz；低截頻值屬性顯示編號(hào)1～4效果最佳，約為6 Hz，但低頻段能量屬性顯示編號(hào)1～5的效果整體比較理想。綜合考慮主頻、頻寬、低截頻值、高截頻值、低頻段能量和高頻段能量，最佳處理參數(shù)是編號(hào)4，主頻為20 Hz，頻寬為44 Hz，低截頻值為6 Hz，高截頻值為50 Hz，約3個(gè)倍頻程。

圖2 頻譜類方法定量評(píng)價(jià)結(jié)果Fig.2 The quantitative evaluation results of spectrum

3 實(shí)際資料應(yīng)用效果

提高地震分辨率效果定量評(píng)價(jià)技術(shù)及軟件系統(tǒng)在新疆車排子工區(qū)處理中進(jìn)行了應(yīng)用，實(shí)現(xiàn)整個(gè)提高分辨率處理過(guò)程的全流程監(jiān)控。新疆車排子工區(qū)在沙灣組發(fā)育薄砂體油藏，老成果資料存在縱橫向分辨能力低、反射軸雜亂、層位追蹤困難、空白反射較多、斷點(diǎn)與斷面不清楚的難題，造成低序級(jí)斷層識(shí)別難度大，儲(chǔ)層精細(xì)描述落實(shí)存在一定風(fēng)險(xiǎn)，迫切需要高分辨率處理地震資料的支撐。

圖3是提高分辨率處理前、后地震剖面的效果對(duì)比，圖3a是提高分辨率處理前的原始數(shù)據(jù)，圖3b是采用地表一致性預(yù)測(cè)反褶積(步長(zhǎng)12 ms)和反Q濾波(參數(shù)取20)的多步提頻處理結(jié)果，圖3c是僅采用反Q濾波(參數(shù)取50)的處理結(jié)果。

采用提高分辨率效果定量評(píng)價(jià)軟件對(duì)處理效果進(jìn)行分析，如圖4所示，對(duì)原始數(shù)據(jù)、多步提頻數(shù)據(jù)、反Q濾波數(shù)據(jù)分別編號(hào)1、2、3，其中圖4a為頻譜類定量評(píng)價(jià)結(jié)果，圖4b為信噪比、子波一致性、波形相似性定量評(píng)價(jià)結(jié)果。從圖4a中可以看出：編號(hào)2數(shù)據(jù)的主頻為27 Hz、頻寬為42 Hz、低截頻值為7 Hz、高截頻值為49 Hz，均優(yōu)于編號(hào)3數(shù)據(jù)；編號(hào)2數(shù)據(jù)的低頻段能量和高頻段能量較強(qiáng)，表明低頻段和高頻段信息豐富。從圖4b中可以看出：提頻處理后信噪比略有降低，但在可接受范圍；子波一致性和波形相似性方面，在中頻段兩種方法基本接近，在高頻段編號(hào)2數(shù)據(jù)的子波一致性和波形相似性更優(yōu)。分析結(jié)果顯示多步提頻優(yōu)于反Q濾波，從地震剖面中也可以得到相同的結(jié)論，采用多步提頻的處理結(jié)果較好地保持了河道砂體的特征(紅色橢圓所示)，并且具有較高的保幅性。

在提高地震分辨率處理過(guò)程中，進(jìn)行了全流程的效果定量評(píng)價(jià)，確保優(yōu)選最佳提高分辨率方法和處理參數(shù)，從而獲取最優(yōu)的提高分辨率結(jié)果數(shù)據(jù)。圖5是提高分辨率處理前、后的效果對(duì)比，從提高分辨率前、后的地震剖面中可以看出，提高分辨率處理成果數(shù)據(jù)解決了老成果分辨率低、空白反射多、層位難以追蹤的問(wèn)題，利用提高分辨率處理數(shù)據(jù)可以滿足砂組細(xì)分的需求。

a—原始數(shù)據(jù);b—多步提頻處理結(jié)果;c—反Q濾波處理結(jié)果a—original processing data;b—multistep high resolution data;c—inverse Q filtering data圖3 提高分辨率效果對(duì)比Fig.3 Improved resolution effect contrast

a—頻譜定量評(píng)價(jià)結(jié)果;b—信噪比、子波一致性、波形相似性頻譜定量評(píng)價(jià)結(jié)果a—the quantitative evaluation results of spectrum;b—the quantitative evaluation results of S/N,wavelet consistency,waveform similarity圖4 定量評(píng)價(jià)結(jié)果Fig.4 The quantitative evaluation results

圖5 提高分辨率前(a)后(b)效果對(duì)比Fig.5 The effect contrast of high resolution processing before(a) and after(b)

4 結(jié)論

1)針對(duì)提高地震分辨率效果監(jiān)控的迫切需求，創(chuàng)新研究了提高地震分辨率效果評(píng)價(jià)技術(shù)，從9個(gè)方面的屬性特征對(duì)不同處理方法和參數(shù)進(jìn)行定量評(píng)價(jià)，實(shí)現(xiàn)了提高分辨率評(píng)價(jià)由定性向定量的轉(zhuǎn)變。

2)利用效果定量評(píng)價(jià)技術(shù)對(duì)提高分辨率處理進(jìn)行全流程的監(jiān)控，最終建立針對(duì)不同地質(zhì)目標(biāo)的最佳處理流程，能夠大幅提高薄層砂體的識(shí)別能力。

3)隨著高密度地震技術(shù)的不斷發(fā)展，地震數(shù)據(jù)量呈指數(shù)型增加，提高分辨率效果定量評(píng)價(jià)技術(shù)既能夠提高處理效率，還可以應(yīng)用于海量數(shù)據(jù)的質(zhì)量監(jiān)控和其他處理環(huán)節(jié)的效果評(píng)價(jià)中，具有廣闊的應(yīng)用前景。

致謝：感謝中國(guó)石化勝利油田物探研究院“提高薄互層地震分辨率處理技術(shù)推廣應(yīng)用”項(xiàng)目組的大力支持！