山前帶聯合層析反演地震成像

王延光 尚新民 趙勝天 滕厚華 閆 妍 朱憲懷*

(①中國石化勝利油田分公司，山東東營 257022；②Forland Geophysical Services，美國休斯頓 77079)

0 引言

山前帶構造具有良好的圈閉條件，常常是油氣資源富集的重要場所。近年來，中國正在加強包括塔里木盆地、四川盆地和準噶爾盆地在內的山前帶地震勘探。但由于地表地形和地下地質構造的復雜性，山前高陡構造成像一直是地震勘探的難題。山前帶地震成像的主要難題是：①近地表速度估算和靜校正；②去噪；③速度建模；④從真地表出發的各向異性疊前深度偏移。

像中國的準噶爾盆地南緣、塔里木盆地庫車、四川盆地鎮巴和南美洲的安第斯山脈，地表條件極其復雜，實施地震勘探十分困難[1-3]。在這些地區進行地震數據采集，檢波器線距往往比期望布置的稀疏，炮線也很難按施工設計實施。在堅硬的山地巖石上埋置檢波器也是個問題，有時只能埋置在遠離設計檢波點的位置。由于地表采樣間隔過大，地震道低速噪聲假頻嚴重，很難用常規去噪方法壓制[3]。另外，山前帶的近地表問題已經不再是簡單的靜校正問題，而是需要從真地表出發用疊前深度偏移加以解決。

本文通過對相關方法的理論分析，結合準噶爾盆地南緣山前帶米泉三維地震成像的綜合應用研究，展示了CWT抗假頻去噪和聯合層析速度反演在復雜山前帶地震成像中的優勢和應用前景。并期望對其他盆地復雜地表地區的油氣勘探也具有一定的指導意義。

1 工區背景和地震數據特征

準噶爾盆地南緣米泉區塊位于烏魯木齊市東北部(圖1)博格達山西北緣。工區整體地形呈西北低、東南高的走勢，海拔高程為500～3500m。地表起伏大，巖性變化也大，包括四類不同的表層地震地質條件[4]：①北部農田戈壁區，地形相對平緩，地表主要分布黃土和粒徑為1～5cm的礫巖；②中北部低山區，海拔相對較高，地表為黃土和粒徑為2～10cm的礫巖；③中南部山區，海拔較高，地表覆蓋有1～25m厚的黃土和礫巖，礫巖的直徑在4～20cm，個別達到1m左右；④南部高山區，不僅海拔高，而且高差變化大，山勢陡峭，表層為1～5m厚的黃土和礫巖，部分出露二疊系、石炭系的泥灰巖、灰巖，地層傾角達50°～70°，個別地方甚至近于直立，表層風化較嚴重。

圖1 工區構造位置圖

圖2是工區典型的炮記錄，可以清楚地看到近地表劇烈變化引起的靜校正和噪聲干擾。除了存在幾套不同速度的面波干擾外，近炮檢距道的鳴震(M)異常振幅嚴重，低速線性噪聲假頻(A)嚴重，給信號處理帶來了很大的困難。

圖2 工區典型地震單炮記錄

米泉山前帶地下構造屬于典型的雙復雜構造[4-5]，淺表層和深層都很復雜：構造上處于盆地結合部位，存在多處沖斷，變形嚴重。區內發育3～4套推覆構造，變形強烈，地震成像難度很大[6]。

針對米泉地區近地表和地下構造的特征以及地震波場的復雜性，下面重點分析處理中應用的兩項關鍵技術：①CWT抗假頻去噪；②聯合層析速度反演和各向異性疊前深度偏移。

2 復數域小波變換(CWT)去噪

在信號處理中，將數據從一個域變換到另一個域，可使噪聲和信號更容易分開，從而達到有效壓制噪聲的目的。本文采用復數域小波變換(Complex Wavelet Transform，CWT)，將時間—空間函數變換成復數域小波基函數，先在復數域壓制識別出的噪聲，然后再變換回時間—空間域。該方法的最大優點是抗假頻、保真，在壓制噪聲的同時能保留有效信號，尤其是低頻信號。CWT去噪理論在Yu等[7]的文章中有詳細介紹，這里就不贅述。不同之處是，為了同時壓制近炮檢距高能量鳴震和全炮檢距頻散線性噪聲，本文綜合采用了一維和二維復數域小波變換。一維復數域小波變換旨在保持有效信號的同時消除強振幅野值。在實際應用過程中，增加了對實際資料進行尺度和方位的人機連作實時分析，這樣更有利于參數的選擇。

圖3是二維復數小波變換在地震處理中的應用流程。它既可以應用于疊前炮集去噪，也可以對疊前時間偏移的共反射點(CRP)道集或疊前深度偏移的共成像點(CIG)道集去噪，甚至可以對疊加數據體或共炮檢距剖面去噪，可以根據需要重復應用。但應用于炮集數據時，應該注意保留具有一定斜率的斷面波和繞射波，以免影響偏移歸位。

圖3 二維復數域小波變換在地震處理中的應用流程

為了方便理解，用一個合成地震炮集記錄描述二維小波變換的過程。圖4展示的是一個二維合成地震炮集數據(圖4a)，通過小波變換分解為六個傾角數據(圖4b從左向右分別為-75°, -45°, -15°, 15°, 45°, 75°)。對每個傾角的地震響應(例如-15°,15°)又通過復數域小波變換分解為1～7七個尺度數據(圖5)，尺度0是分解前的數據。尺度越大，視頻率(時間方向的頻率或空間頻率(波數))越低。其他傾角可以做同樣的分解。具體應用中可以選擇不同的傾角和尺度，目標是保留信號、壓制噪聲。在三維情況下，傾角是一個傾斜面，包括方位和斷面傾角。與其他方法(例如FK變換、拉冬變換等)相比，二維小波復數域變換具有以下優點[7-9]：①保真；②可逆；③局部化；④抗假頻；⑤便于自動化。

圖4 二維復數小波變換把二維地震炮集數據(a)分解為 -75°，-45°，-15°，15°，45°，75°六個傾角數據(b)

圖5 二維復數小波變換把二維傾角炮集數據再分解為不同尺度的數據

地震信號處理的難點在于，壓制噪聲的同時既要保留有效信號又能抗假頻。二維復數域小波變換具有這個能力。圖6和圖7是去噪前、后一個實際炮記錄的對比。從炮記錄和相應的振幅譜可以看出，復數域小波變換濾波不僅可以有效地壓制含假頻的低速噪聲，而且有效信號得以保留，尤其是低頻信號(圖7)，而低頻信號對油氣藏分析十分重要。

圖6 實際原始炮點記錄(左)和相應的振幅譜(右)

圖7 經過復數域小波變換濾波后的炮記錄(左)和相應的振幅譜(右)

圖8是復數域小波變換去噪前、后數據疊加剖面對比。可見去噪后剖面的信噪比得到明顯改善，尤其是低頻信號得到加強。這進一步證明了二維復數域小波變換在山前帶壓制噪聲的有效性。

圖8 應用復數小波變換去噪前后的疊加剖面對比

3 聯合層析速度建模和各向異性疊前深度偏移

復雜山前帶地震成像的另一項關鍵技術是速度建模，包括近地表速度建模和淺、中、深層的聯合速度建模[1-2]。

回轉波層析反演(Turning-ray tomography),又稱潛水波層析反演(Diving-wave tomography，DWT)，最早由Zhu等[10]、Stefani[11]、Zhang等[12]提出，是一種行之有效的近地表速度建模方法。但該方法的可靠反演深度一般為最大排列長度的1/8 ～1/4(取決于地下速度的梯度大小和是否存在速度反轉[2,10])，這顯然不能滿足目的層的勘探要求。另一方面，反射波層析反演可以在深度偏移后通過反射波的剩余時差來更新較深的速度模型，但對淺層卻無能為力，因為在淺層缺少同相軸或者道數太少(圖9左)。針對這種情況，Zhu等[13]提出了早期的聯合層析反演思路，即綜合考慮折射波(回轉波或潛水波)層析和反射波層析各自的優點，互相補充(圖9)。首先通過回轉波初至層析反演得到比較可靠的淺層速度模型，與已知的深層速度模型(例如以往的深度偏移速度模型或者從疊前時間偏移的均方根速度轉換成層速度模型)進行融合，作為反射波層析反演的初始速度模型，然后進行反射波層析反演。在反射波層析反演的前幾次迭代過程中，淺層速度約束不變；最后幾次迭代時，放棄淺層速度約束，僅根據反射層析經過的射線路徑進行更新，這是因為初至走時反演得到的淺層速度可能不夠精確，它一方面受到初至拾取誤差的影響，另一方面受到各向異性介質的影響。通過各向異性(TTI)反射波層析反演，可以適當減小淺層速度模型的誤差[1-2]。工業界多年的實踐證明，上述的兩步法是行之有效的[1-2](圖10)。

圖9 早期的聯合層析反演思路 [13-14]

聯合層析速度建模的起始速度是把淺層層析速度體與深層平滑速度體合并(圖10)。首先，兩者合并時需要設計一個過渡帶，使得合并后沒有痕跡。其次，隨著深度的增加，射線角度變窄，使得回轉波層析反演結果邊界效應隨著深度的增加而加大。因此，淺層層析反演速度往往需要進行編輯、平滑、擴展等一系列處理，然后再進行融合(圖11)。同時，由于淺層各向異性等因素，層析反演得到的淺層速度往往偏高，因此還需對層析成像的淺層速度模型進行修改，包括掃描更新和VSP標定，最終以疊前深度偏移后的成像點道集是否拉平為原則進一步優化調整。

圖10 實現聯合層析反演起始速度融合示意圖。可靠深度是指淺層層析反演的可靠深度，一般與排列長度有關

圖11 聯合層析速度模型起始速度模型與PSDM疊加融合顯示

為了加快速度更新和提高反演的穩定性，聯合層析反演可以在各向同性疊前深度偏移基礎上，先近炮檢距進行各向同性反演。 在近炮檢距反射同相軸拉平后，通過測井資料與解釋層位進行井震標定，然后加大炮檢距進行各向異性參數選擇和各向異性疊前深度偏移，提高地震成像的精度。

需考慮的一個實際問題是，當檢波點距比較大時，折射波和反射波層析之間的“過渡帶”地震道數太少(圖12)，不利于聯合層析反演。經過五維插值(圖12a)，覆蓋次數得到增加，從而增加了淺、中層之間“過渡帶”的近炮檢距道數，聯合層析反演的可靠性得到改善。

圖12 用五維插值增加覆蓋次數，增加淺、中層之間“過渡帶”(橢圓框內)的近炮檢距道數

4 應用效果

米泉三維地震采集的最大排列長度為7000m，部分炮檢距達10000m(圖13a)。通過152口微測井(圖13b)約束，初至層析反演的速度如圖13c所示。根據射線密度的底界(圖13d), 可靠深度(從地表往下)達1500m左右，給聯合層析反演提供了很好的淺層初始速度模型。用微測井數據約束層析反演和方位加權插值得到一個比較精細的三維表層(0～50m)速度體[15]，進一步提高初至波層析反演的精度。但是，如果微測井數據使用不當，經過約束反演后反而會產生所謂的“牛眼”現象。為此先把微測井速度通過空間內插得到一個三維表層速度體，這樣可以提高約束反演的穩定性。如果微測井的深度大部分只有10～20m，則用它進行約束反演的意義不大，因為淺層層析反演網格的垂向增量一般為10m。

圖13 淺層回轉波約束層析反演

圖14是某一縱測線的TTI各向異性參數模型，其中速度是垂直于地層的層速度。各向異性參數δ和ε分別控制地震成像在垂向和橫向上的校正量。在沒有橫波信息的情況下，δ通過鉆井資料、層位劃分和井震對比來估算；ε通過遠炮檢距地震數據估算。傾角和方位角參數模型由反射波層析反演自動計算得到，對各向異性疊前深度偏移和提高傾斜層的成像精度起到非常重要的作用。

圖14 聯合建模更新后的各向異性參數

圖15是分別用初始速度模型和經聯合層析反演更新得到的速度模型進行疊前深度偏移的共成像點道集。雖然工區三維地表條件和地下構造復雜，噪聲嚴重，但通過CWT去噪、聯合層析反演和TTI各向異性疊前深度偏移，共成像點道集基本得到拉平(圖15右)。評價CWT去噪、聯合層析反演和TTI各向異性疊前深度偏移的效果還要看最終疊加剖面。圖16是重新處理(左)與以往處理的疊前深度偏移疊加剖面(右)的對比。重新處理包括CWT去噪、聯合層析速度建模和從真地表出發的TTI各向異性疊前深度偏移。顯然，通過重新處理，低頻信號得到加強，構造形態和斷裂分布更加清晰，有利于后續的解釋工作。

圖15 采用初始速度模型(左)與聯合層析速度更新后模型(右)的疊前深度偏移道集

圖16 疊前深度偏移(轉換為時間)剖面對比

5 認識及展望

通過CWT去噪和聯合層析速度建模在米泉山前帶地震成像中的成功應用，得到以下認識：

(1)三維復數域小波變換可以壓制層間多次波。利用復數域小波變換抗假頻和高保真的優點，并考慮到多次波的傳播路徑在二維和三維空間上的不同，再結合地質構造在不同炮檢距剖面上應該一致的假設條件，可以有效地壓制層間多次波。

(2)目前陸地的全波形反演(FWI)還不成熟，主要問題在于表層0～30m內巖性極其不均勻和地震波散射，引起地震子波極其不穩定和嚴重噪聲干擾。需要通過適當的預處理，改進子波的穩定性和壓制噪聲，用FWI提高淺層速度的精度。

(3)目前聯合層析反演一般分成兩步或多步，在連接上往往需要處理員有一定的經驗才能滿足生產的要求。期望在不久的將來，折射和反射聯合層析反演能夠一步實現。

在本課題的研究中，焦俊如和楊斌做了大量的CWT去噪方法研究和試驗，George McMechan和Sam Gray為聯合層析反演提供了很好的建議，Allen Bertagne為速度建模提供了構造解釋，Ellen Chen參加了時間域重新處理，許琨和馬青提供了軟件開發和技術支持，對此表示感謝。