劉道理 李 坤 楊登鋒 魏旭旺

1.中海石油（中國）有限公司深圳分公司研究院 2.中國石油大學（華東）地球科學與技術學院

0 引言

基于振幅和頻率信息的疊前地震反演是深層儲層流體特征識別的有效途徑[1]。自20世紀70年代以來，國內(nèi)外專家學者依據(jù)巖石物理理論構建儲層流體敏感參數(shù)，疊前地震反演技術在流體檢測及應用方面得到了長足的發(fā)展。目前，疊前地震反演技術大多是基于地震道集的AVO振幅變化信息，該方法在地震解釋和儲層描述中發(fā)揮了重要的作用，然而隨著地震勘探程度的日益提高，僅僅依靠振幅信息難以滿足高精度地震勘探的要求?，F(xiàn)階段，針對復雜的“多重孔—裂隙結構”與“波動誘導流體流動”的地震巖石物理衰減模型構建及分析方法得到了國內(nèi)外地球物理學家的廣泛關注，旨在用于解釋地震波振幅衰減及速度頻散機理，且利用巖石彈性參數(shù)的頻散特征進行儲層描述已成為研究熱點。前人的研究成果表明，孔隙內(nèi)賦存多重流體介質(zhì)時會導致地震振幅衰減和速度頻散[1]。因此，在地質(zhì)背景更加復雜、儲層含流體多樣化的情況下，聯(lián)合地震振幅以及頻率等相關信息，指導儲層預測和流體表征是一項很有意義的研究課題。Taner等[2]在進行復地震道分析時，通過提取瞬時頻率、瞬時振幅屬性發(fā)現(xiàn)了油氣藏下方頻率向低頻方向移動的現(xiàn)象，將其稱作低頻陰影，為利用地震振幅衰減和頻散信息的儲層油氣識別技術奠定了基礎。Castagana等[3]率先將高分辨率匹配追蹤瞬時譜分解技術用于含油氣儲層相關的地震波低頻陰影檢測，并展示了匹配追蹤譜分解技術、低頻陰影識別技術在地震資料分析解釋中良好應用前景，但未明確給出低頻陰影的形成機制。

速度頻散是指地震波在地下介質(zhì)中傳播時相速度隨頻率變化的現(xiàn)象，振幅衰減與速度頻散現(xiàn)象相伴生，國內(nèi)外學者從巖石物理機制上明確了孔隙流體對地震波衰減以及速度頻散現(xiàn)象的影響，為基于頻散屬性的流體識別方法研究奠定了理論基礎[4-9]。Wilson等[10]依據(jù)Chapman裂隙—孔隙微觀結構衰減巖石物理模型，聯(lián)合譜分解技術與Smith—Gidlow近似公式，提出了基于縱波速度頻散屬性的頻變AVO反演方法；Wu Xiaoyang等[11]結合重排平滑偽維格納分布時頻譜分解技術，實現(xiàn)了均勻介質(zhì)情況下的頻變AVO反演應用；張世鑫等[12-13]在疊后地震數(shù)據(jù)反演與Shuey線性AVO近似的基礎上，通過求取縱波速度隨頻率的偏導數(shù)，反演得到頻變縱波速度，并應用于砂巖儲層含氣性識別領域；程冰潔等[14]在均勻介質(zhì)AVO反射特征方程的基礎上，將頻變AVO流體識別技術應用到實際地震資料處理中，在含氣儲層的識別中取得了理想效果；Zhang等[15]基于Biot—Gassmann孔隙彈性理論，在Russell線性AVO近似公式的基礎上，依據(jù)貝葉斯估計框架提出了頻變流體因子地震AVO反演方法。

筆者在前人研究的基礎上，綜合利用巖石物理分析、譜分解等方法，探討了頻變AVO反演方法在我國近海某盆地P探區(qū)的深層儲層流體檢測中的應用效果。針對工區(qū)內(nèi)深層儲層含氣性難以辨識的問題，考慮到僅利用地震振幅信息進行流體檢測存在多解性，研究和發(fā)展了基于頻變屬性的高精度流體識別方法。實際應用結果表明，敏感的頻變流體因子屬性能夠有效區(qū)分儲層中含流體的類型，消除常規(guī)瞬時譜流體檢測方法造成的流體識別假象，進一步降低了勘探風險。

1 基于FAVO反演的頻變流體因子反演原理

1.1 頻變彈性參數(shù)流體敏感性分析

頻變流體因子相比常規(guī)流體識別因子，充分考慮了巖石彈性參數(shù)隨頻率的變化特征，且該頻變特征與儲層含油氣性質(zhì)相關，具有區(qū)分孔隙含流體類型能力的特征參數(shù)。Chapman等[8]在橢球型孔隙及扁平狀裂隙介質(zhì)的基礎上，考慮了巖石孔隙內(nèi)由于壓力不均衡引起的局部流體流動現(xiàn)象，建立了裂隙—孔隙微結構衰減地震巖石物理模型，并推導了利用巖石基質(zhì)、流體黏滯系數(shù)及角頻率等參數(shù)定量表征的飽和巖石等效體積模量和剪切模量，具體表達式為：

式中Keff表示飽和巖石的等效體積模量，Pa；K表示巖石基質(zhì)的等效體積模量，Pa；ε表示裂隙密度參數(shù)；λ表示巖石基質(zhì)的拉梅參數(shù)，Pa；μ表示巖石基質(zhì)的剪切模量，Pa；r表示孔隙形狀的縱橫比；φ表示巖石的孔隙度；μeff表示飽和巖石等效剪切模量，Pa；i表示虛部；ω表示角頻率，rad/s；τ表示時間松弛因子；Kc、A(ω)和B(ω)表示與上述參數(shù)相關的系數(shù)；k表示為巖石滲透率，mD；η表示流體黏滯系數(shù)，Pa·s；ν表示巖石基質(zhì)的泊松比；ζ表示礦物顆粒尺寸，mm；a表示裂隙半徑，mm；Kf表示流體體積模量，Pa。

根據(jù)Chapman裂隙—孔隙微結構衰減理論，針對飽和流體的巖石的Gassmann流體項、拉梅參數(shù)、體積模量、泊松比、縱波速度、剪切模量、橫波速度等7種彈性參數(shù)的頻變特征及流體敏感性進行了分析，對含氣砂礫巖儲層彈性參數(shù)的頻變特征進行數(shù)值模擬（圖1），分析結果表明，各彈性參數(shù)的歸一化頻變程度大小為：

式中f表示Gassmann流體項，Pa；λsat表示拉梅參數(shù)，Pa；μsat表示剪切模量，Pa；σ表示飽和流體的巖石的泊松比；νp表示縱波速度，m/s；νs表示橫波速度，m/s。

1.2 地震信號時頻分解方法

地震信號在地下介質(zhì)中傳播過程時，振幅和頻率的統(tǒng)計特征是隨時間發(fā)生變化的。時頻譜分解是揭示地震數(shù)據(jù)中所包含的時間和頻率信息的重要手段，如圖2所示，短時傅里葉變換在整個變換過程中時窗長度不變，因此具有單一的時頻分辨率的局限性；連續(xù)小波變換的基函數(shù)是一個自適應的窗口，窗口長度隨頻率增高而縮小，低頻段頻率分辨率較高，高頻段時間分辨率較高，具有多重分辨率的特性。實際應用中綜合考慮時頻分辨率以及計算效率的影響，筆者選用連續(xù)小波變換（Continuous Wavelet Transform，簡稱CWT）應用到頻變流體因子疊前反演中。針對非平穩(wěn)信號f(t)，則其連續(xù)小波變換被定義為：

圖1 頻變彈性參數(shù)的流體敏感性分析圖（氣砂模型）

圖2 傅立葉變換、短時傅里葉變換和連續(xù)小波變換的不同尺度基函數(shù)示意圖

1.3 頻變AVO貝葉斯反演算法

借鑒Wilson等[10]提出的縱波速度頻變屬性反演方法，Zhang 等[15]根據(jù)Russell近似公式建立了頻變流體因子反射系數(shù)方程，并發(fā)展了基于貝葉斯理論的頻變AVO反演算法?；贐iot—Gassmann理論，Russell等[16-17]對飽含流體孔隙介質(zhì)的AVO理論進行了研究，并推導了反射系數(shù)近似公式為：

式中ρ為巖石密度，kg/m3；Δf表示Gassmann流體項的絕對變化量，Pa；Δμsat表示剪切模量的絕對變化量，Pa；Δρ表示密度的絕對變化量，kg/m3；γdry表示干巖石骨架的縱橫波速度比；γsat表示飽和巖石的縱橫波速度比。

由于Gassmann流體項、剪切模量的頻散效應，且地層密度不發(fā)生頻散，將式（5）在參考頻率ω0處進行一階泰勒展開，并令，將上式改寫為：

式中Df表示Gassmann流體項反射率的頻散程度；Dμ表示剪切模量的頻散程度。

考慮到地震記錄是地震子波與反射系數(shù)的褶積，引入非平穩(wěn)頻變子波的概念，即

式中頻率分量ωi和角度分量θ為頻變AVO反演的輸入信息，為了更容易討論反演問題，將式（6）和式（7）改寫成矩陣方程的形式，即

式中N表示輸入的部分角度疊加地震數(shù)據(jù)的數(shù)量，L表示選取的頻率個數(shù)；G表示頻變AVO反演的核矩陣算子；d表示頻率ωi與ω0在第j個入射角時對應的振幅差異向量；表示頻率ωi與ω0在第j個入射角時對應系數(shù)矩陣；為頻率ωi且入射角為θj時對應的子波矩陣，R為待反演頻變屬性組成的列向量。

為了提高頻變AVO地震反演算法的穩(wěn)定性，筆者假設地震AVO數(shù)據(jù)中觀測噪聲服從正態(tài)分布pGauss(d∣R)，待反演模型參數(shù)服從柯西概率密度分布pCauchy(R)先驗信息，通過貝葉斯公式將模型參數(shù)的柯西先驗信息和高斯似然函數(shù)聯(lián)系起來[18-19]，即

求取式（9）的貝葉斯最大后驗概率密度解（MAP），可以得到目標泛函式為：

式中JG(R)表示服從高斯分布地震數(shù)據(jù)的懲罰項，JCauchy(R)表示由柯西概率密度分布引入的正則化項，σn2表示地震數(shù)據(jù)噪聲的方差，M表示采樣點個數(shù)，σm2表示待反演模型參數(shù)R的方差。

通過目標函數(shù)最優(yōu)化得到最終反演方程式為：

針對上述弱非線性求解問題，筆者依據(jù)反復重加權最小二乘法算法求解式（11），即可實現(xiàn)頻變流體因子的提取。頻變流體因子提取需要聯(lián)合譜分解技術和頻變AVO反演理論，與基于疊前振幅信息的AVO反演方法相比，頻變流體因子反演方法對疊前地震資料品質(zhì)要求更高，需要地震處理人員在地震資料處理過程中最大限度的保存實際波場的AVO振幅與頻率異常信息。

2 實例驗證

為了驗證基于振幅和頻率信息的頻變AVO反演在深層流體識別中的可行性，針對我國近海某盆地P探區(qū)的深層氣藏（有探井A、B、C等3口井，目的層埋深超過4 000 m）的實際資料進行處理。3口井目的層的儲層巖性以中—細砂巖、礫狀砂巖為主，儲層厚度橫向變化大（2～10 m）、孔隙度低（小于10%）、滲透率低（小于10 mD），儲層頂部被一套低速度、低密度的膏鹽層封隔（膏鹽層厚度累計厚度超過200 m），構成了物性變化與膏鹽層封堵并存的隱蔽性深層氣藏。

圖3為過A井、C井的小角度部分疊加地震數(shù)據(jù)剖面（5°～10°）與測井解釋結果，圖中，紅色表示氣層解釋，綠色表示水層，黑色橢圓指示區(qū)域為砂礫巖扇體發(fā)育位置，圈內(nèi)上部的砂礫巖扇體為氣層發(fā)育位置，黑色箭頭指示位置為膏鹽層形成的地震強反射同相軸。從圖3中可見，氣層發(fā)育位置的地震反射振幅較弱，受到膏鹽層地震強反射的影響較嚴重。

圖3 某盆地P探區(qū)過A、C井的小角度部分疊加地震剖面圖

圖4 基于疊前地震反演的Gassmann流體項預測結果圖

圖4為Gassmann流體項反演結果，從圖4中可見，基于振幅的地震反演往往難以實現(xiàn)隱蔽氣藏的有效識別。因此，利用筆者提出的頻變AVO反演方法開展該探區(qū)的深層氣藏的識別，識別結果如圖5所示，該圖為本方法反演得到的頻變Gassmann流體因子剖面圖。

圖5 本方法反演得到的頻變Gassmann流體因子剖面與測井解釋結果圖

據(jù)圖4、5對比分析可知，對比不同參數(shù)對砂礫巖氣藏的識別效果，Gassmann流體因子和頻變流體因子均能夠有效地識別A井含氣儲層，C井為干層，與測井解釋一致。但由于鹽膏層強地震反射軸的存在，相比于Gassmann流體因子反演剖面來看，頻變Gassmann流體因子反演結果受到膏鹽層強地震反射振幅的影響相對較低，具有更高辨識深層氣藏的能力。分析可得，基于疊前地震振幅和頻率信息的頻變流體因子反演在一定程度上降低了膏鹽層對深層氣藏流體識別的影響，提高了疊前地震流體識別的精度。

圖6為頂部砂礫巖扇體的頻變Gassmann流體因子平面切片，圖6中可以看出頻變流體因子清晰地刻畫出A、B井氣層發(fā)育區(qū)域（B井為驗證井），表明基于振幅和頻率信息反演的頻變流體因子可以作為深層氣藏的指示因子，這為深層天然氣層的識別提供了新思路與新方法。

圖6 P探區(qū)頻變流體因子平面展布圖

3 結論

1）基于Chapman裂隙—孔隙微結構衰減理論模型的彈性參數(shù)頻散特征分析對構建流體敏感參數(shù)的頻變特征參數(shù)具有理論指導意義，為頻變AVO反演與深層含氣性識別奠定了理論基礎。

2）連續(xù)小波變換時頻分析方法可以實現(xiàn)部分角度疊加地震數(shù)據(jù)的瞬時譜分解，為頻變AVO反演提供了單頻瞬時譜數(shù)據(jù)，奠定了該方法的數(shù)據(jù)基礎。

3）采用基于柯西約束的貝葉斯頻變AVO反演算法可以實現(xiàn)基于地震資料的敏感頻變特征參數(shù)提取。實際資料處理表明，該方法充分利用了地震資料中蘊含的疊前地震振幅和頻率信息，一定程度上降低了深層儲層地震流體檢測的多解性，提高了深層儲層含氣性識別的精度。