散射成像在小尺度縫洞體識別中的應用
——以川中古隆起燈影組為例

姜曉宇 宋 濤 甘利燈 戴曉峰 丁 騫 周曉越

(中國石油勘探開發研究院，北京 100083)

0 引言

2011年高石梯1井取得突破，掀起了川中古隆起燈影組氣藏勘探開發的熱潮。受沉積及巖溶作用共同控制的燈影組縫洞型儲層具有非均質性強、縫洞體尺度小的特點。儲層原生基質孔隙欠發育，次生溶蝕縫洞發育是氣井高產的主要因素[1]。受埋藏深、短時期風化巖溶作用影響，燈影組主要發育小尺度巖溶縫洞單元，地震剖面響應特征不明顯，是燈影組氣藏縫洞單元地震刻畫的主要困難。

目前，國內、外學者在對塔里木盆地奧陶系的研究中積累了大量利用地震資料預測碳酸鹽巖巖溶縫洞體的經驗[2-6]。大多以應用常規疊后地震資料的方法為主，對大尺度縫洞體取得了較好的預測效果，對于小尺度縫洞，則受到限制。當地下存在溶洞、裂縫、地層尖滅點等不連續體時，將會產生散射波(繞射波)[7]。散射成像利用散射波信息，能夠突出不連續地質體，增強小尺度縫洞體的地震響應，對于洞穴、裂縫等不連續地質體的準確成像具有重要意義。然而，常規的地震處理技術把反射波作為主要的研究對象，而能量較弱的散射波通常被視為噪聲被去除。即使對散射波進行偏移歸位，也常被能量較強的反射波掩蓋，導致目標繞射體的成像分辨率較低。為此要充分利用地震資料中的散射波信息，需要應用專門針對散射波的處理流程，或者將散射波從全波場中分離出來實現單獨成像。

全方位局部角度域成像技術是在地下局部角度域(Local Angle Domian，LAD)實現Kirchhoff偏移。LAD成像的概念[8]最早提出時，由于炮檢距域成像方法的流行并未受到重視。之后，Audebert[9]完善了LAD系統，給出了LAD的4個角度參數[9]；Xu等[10]提出了針對復雜介質的共散射角成像方法；Brandsberg-Dahl等[11]總結了散射角、方位角域地震速度分析方法，并用于地震成像和速度建模。Koren等[12]在全方位角成像技術的理論研究和應用推廣方面首次將共反射角偏移引入LAD；受散射成像方法的影響，Koren等[13]在LAD全方位成像理論框架下闡述了傾角道集和反射角道集的產生方法及其用途；Koren等[14]系統闡述了全方位LAD成像技術的原理，并展示了其在縫洞成像和儲層描述的顯著優勢。中國學者對LAD成像方法研究也做了很多工作。Cheng等[15]基于對地震波局部方向特征的分析，介紹了方位保真局部角度域Kirchhoff疊前時間偏移成像方法以及在各向同性、方位各向異性分析中的應用；段鵬飛等[16]研究了基于射線的LAD偏移算法，用TI介質對比分析了兩種不同的射線追蹤方法，并進一步發展了LAD高斯束疊前深度偏移方法。

為了增強燈影組氣藏小尺度縫洞體地震反射特征、預測小尺度縫洞體的分布，本文基于全方位LAD深度偏移處理方法，對全方位共傾角道集進行散射加權處理，得到散射成像數據體，能夠突出小尺度縫洞體的地震反射特征。根據小尺度縫洞體在散射成像數據剖面上表現的強能量特征，提取散射成像數據的能量體屬性，刻畫縫洞體儲層的分布特征。

1 全方位LAD成像方法

全方位LAD成像方法將地表記錄地震數據映射到地下成像點的LAD，從成像點向地面進行射線追蹤成像，所有的射線都參與成像，避免了散射能量的損失，同時記錄了地下介質成像點的真正角度域信息(方位、方向、能量、傾角等)，產生兩個互補的全方位、三維的角度域道集：共反射角道集和共傾角道集。

全方位LAD處理方法有如下優點：①該方法能克服共炮檢距道集中存在的假象，提高地震成像精度； ②由于成像的優勢，該方法得到的道集包含了地層全部真實的方位角信息，并且能夠得到真實的振幅異常，繼而反映出地下速度和巖性變化[17-18]； ③全方位共傾角道集能夠根據地下不同反射波場區分出反映連續界面反射特征的鏡像能量和反映不連續反射特征的散射能量，通過加權疊加得到散射成像數據體，突出洞穴等非連續地質體的幾何特征。

對于地下連續界面上的反射點，在傾角道集上表現為“碗狀”結構，把連續界面上反射點所產生的能量叫做鏡像能量；對于地下獨立繞射點，在傾角道集上顯示為近似的一條直線，而對于斷層所對應的棱鏡波，會在傾角道集固定的方位角上表現為能量的規則分布特征，它們所對應的能量統稱為散射能量，反映的是地下溶洞、特殊巖性體、斷裂等地質異常信息。將這些信息分離出來進行分析和應用，對于研究區小尺度縫洞體的刻畫有重要意義。

鏡像成像是垂直法線的鏡像反射能量的疊加體，是一個將真實反映地層傾角的共反射點道集進行疊加的過程，與常規克希霍夫偏移法得到的疊后地震資料相比，全方位LAD的偏移成像方法得到的鏡像成像數據體能夠顯著提高資料的信噪比(圖1)。

圖1 常規克希霍夫(a)與全方位LAD(b)偏移剖面對比

散射成像體是在壓制鏡像能量后的疊加剖面，主要體現非連續地質特征的成像。一般情況下，為了加強儲層界面的成像質量，會適當增大鏡像疊加權重，壓制散射能量，本文為了反映小尺度縫洞體等不連續現象，經過適當的減小鏡像權重進行散射能量加權疊加，得到散射成像體。如圖2所示，gs9井鉆井存在鉆具放空、鉆井液大量漏失現象，顯示該段縫洞體發育。常規疊后成像突出了層界面的成像連續性，而層內小尺度縫洞體沒有明顯的地震響應。散射成像去除了頂面強反射的影響，突出了內部縫洞等不連續體的地震反射特征，在gs9井處，可見類似串珠狀的地震反射特征(圖2b橢圓圈所示)，為小尺度縫洞體的識別及預測奠定了基礎。

圖2 常規疊后(a)與散射成像(b)剖面對比

2 散射成像數據的應用及效果

2.1 井旁地震道頻譜分析

通過井旁地震道的頻譜成像，可幫助建立儲層特征與振幅譜和相位譜的定量關系，使對頻譜成像處理結果的解釋更具物理意義和地質意義。本文首先通過井點模擬方法確定儲層頻譜成像技術能夠識別研究區的縫洞體儲層。如圖3所示，在縫洞體發育的鉆井液漏失段，散射成像頻譜分解結果表現為在15～45Hz頻帶內的強能量譜特征，因此為了突出縫洞體的能量特征，提取此頻帶內的總能量屬性對縫洞體進行描述。

圖3 井旁地震道的頻譜成像分析

總能量體屬性用地震數據有效頻帶范圍內振幅譜的積分來表征，相當于低截頻到高截頻之間的振幅譜的總面積(圖4)。能量體屬性的計算公式為

圖4 能量體屬性示意圖

(1)

式中：A為振幅；f為頻率；fH為高截頻；fL為低截頻。

根據上述公式計算井旁道的能量體，可以看到，能量體屬性高值與鉆井液漏失段位置十分吻合(圖3)，表征了研究區小尺度縫洞體的發育位置。

2.2 應用效果分析

利用上述能量體計算公式，計算整個工區內散射成像數據的能量體屬性。從過gs9井地震和能量體剖面(圖2b、圖5)可以看出，散射數據提取的能量體屬性清晰地刻畫了散射成像剖面上縫洞體的串珠狀反射特征。分別對常規疊后數據和散射成像數據提取能量體屬性，常規疊后數據計算的能量體屬性反映了常規疊后地震數據剖面的儲層界面的連續、強反射特征，無論是常規疊后地震數據還是由此計算出來的能量體屬性都與不能反映縫洞體的特征(圖2a、圖6)。散射成像計算的能量體屬性與井上漏失段位置相符，能夠反映縫洞體的分布特征(圖5、圖7)。

圖5 過gs9井散射成像能量體剖面

圖6 過gs001-H33井常規疊后剖面(a)及其能量體剖面(b)

圖7 過gs001-H33井散射成像剖面(a)及其能量體剖面(b)

平面上，散射成像能量體屬性預測的縫洞體成分散的點狀分布，符合巖溶作用控制的縫洞體分布地質特征，而常規疊后數據計算的能量體屬性在工區邊界呈塊狀分布，沒有明確的地質意義(圖8)。

圖8 散射成像(a)與常規疊后(b)目的層能量體切片對比

3 結論

本文針對川中燈影組儲層小尺度縫洞體預測難題，從處理解釋一體化研究著手，通過全方位局部角度域處理方法得到全方位共傾角道集，對其進行散射加權處理，得到散射成像數據體，突出了小尺度縫洞體的地震反射特征，根據小尺度洞穴體在散射成像數據剖面上表現的強能量特征，提取散射成像數據的能量體屬性，刻畫了縫洞體儲層的分布特征。得到以下結論：

(1)相比于常規Kirchhoff偏移成像方法，全方位局部角度域偏移方法成像效果更好；

(2)散射成像數據體能夠提高地震數據對不連續地質體的分辨能力，反映小尺度縫洞體等地質現象；

(3)以川中燈影組儲層為例，運用散射成像數據體提取的屬性特征比常規的疊后地震數據提取的屬性特征更能反映縫洞體等不連續地質體的分布特征。