層析靜校正在FZ地區三維地震資料處理中的應用

李　娜

(大慶油田有限責任公司　勘探開發研究院，大慶　163712)

?

層析靜校正在FZ地區三維地震資料處理中的應用

李娜

(大慶油田有限責任公司勘探開發研究院，大慶163712)

FZ區塊近地表條件復雜，導致單炮記錄復雜多變，存在嚴重的靜校正問題。結合以往的處理經驗可知，折射波、曲時面、微測井及高程靜校正效果都不能滿足目標區的靜校正要求。因此，需要采用層析靜校正方法反演近地表結構，獲得較準確的靜校正量，以保證構造的成像精度。這里簡要說明層析靜校正方法的基本原理，介紹FZ區塊近地表條件和單炮記錄的特點，針對FZ區塊資料的特殊性，嘗試不同替換速度對靜校正效果的影響并進行對比分析，從而證明初至波拾取質量及替換速度的選取都起到了非常關鍵的作用。

層析反演； 靜校正； 復雜地表； 替換速度； 初至波

0　引言

FZ地區地表高程變化劇烈，巖性復雜多變，近地表速度縱橫向變化大。地下陡傾角地層及淺層高速反射層導致大部分資料出現早于直達波到達的頂點偏移的雙曲線，部分地區低降速帶由多層組成，導致多次折射發育，并且該地區原始資料信噪比低，冰帽、淺層高速折射和異常有效波交織在一起，初至波能量不均衡，沒有穩定的可進行連續追蹤的折射界面，致使基于層狀介質模型的折射波靜校正技術在應用上受到很大限制。因此，采用初至波層析反演靜校正方法解決FZ地區復雜的靜校正問題。首先此方法可以將初至的直達波、繞射波和折射波綜合考慮，不需要進行區分，避開了折射波靜校正必須首先識別折射波的限制；其次利用非線性層析反演技術，詳細反演每一觀測位置上的速度信息，求取準確的近地表速度模型[1]，從而計算靜校正量，能夠很好地解決地表起伏問題和速度劇烈變化問題對靜校正的影響。

1　層析反演靜校正方法原理

初至波走時層析反演靜校正的過程主要包括：①建立精細的近地表深度-速度模型，為靜校正計算奠定良好基礎；②對模型進行正演模擬計算，采用射線追蹤方法計算初至波旅行時；③用計算得到的初至波旅行時和實際拾取的初至波旅行時進行比較計算，得到模型的修正量，經過反復迭代最終得到比較精確的近地表模型；④利用得到的近地表模型計算靜校正量[2]。

層析反演是將復雜近地表地質模型微元化，并假設單個微元內介質速度穩定不變，介質之間速度不同，用網絡法進行射線正演，獲得表層速度模型，當微元趨于無限小時，認為它能夠真實地描述表層結構模型。

假設一個初始速度模型，第j條射線的旅行時可以表示如下：

(1)

式中：tj為第j條射線的初至波旅行時；si為第i個網格的慢度；dij為第j條射線在第i個網格中的射線路徑長度。

將式(1)離散后, 可得到代數方程組的矩陣形式，如式(2)所示。

T=D S

(2)

式中：T為炮點到檢波點的旅行時矩陣；D為與射線路徑有關的距離矩陣;S為地下介質的慢度矩陣；其中，D矩陣未知，初始速度模型S(0)預先根據微測井、小折射、大折射以及其他地質和地球物理資料構建，然后利用射線追蹤方法正演求得初至波的射線路徑，并且通過式(1)計算得到初始模型下的理論旅行時T(0)，通過比較地震波初至旅行時T和正演得到的旅行時T(0)，求出旅行時殘差矩陣ΔT為式(3)。

ΔT=T-T0

(3)

由線性方程組式(3)可得反演方程：

ΔT=DΔS

(4)

式中，ΔS為速度單元的慢度修正量，射線路徑距離矩陣D已經由初始速度模型求出，因此可以求解出ΔS。

通過得到的慢度修正量ΔS對初始速度模型S(0)進行修正，得到新的速度模型S，再利用正演求得新的射線路徑D和旅行時殘差ΔT，從而求出新的速度模型修正量。以上過程不斷迭代，直至實際旅行時與理論旅行時的殘差ΔT滿足精度要求為止。由此得到與真實速度模型相近似的近地表模型，進而可以求出較為精確的靜校正量[3]。

顯然，初至波層析反演靜校正是正、反演結合的過程，根據實際地震記錄的初至時間反演出近地表速度模型結構，它與初至折射靜校正以及傳統的射線追蹤方法最大的區別是層析法成像采用高密集度的速度單元劃分，因此可以描述更為復雜的速度場。

2　應用實例

2.1關鍵技術

2.1.1初至時間拾取

為解決本區塊復雜的靜校正問題，采用波動方程正演模擬方法，來認識復雜波場的產生機理，從而用來指導丘陵山地部位的初至時間拾取工作。依據實際資料情況，設定模型范圍為10 km×1 km，網格為5 m×5 m，檢波點間距為25 m，檢波點個數400個，炮間距為100 m，總炮數61炮。圖1中，記錄中模擬的A炮和B炮與山地丘陵區的實際單炮記錄相似，通過研究發現，“異常有效波”是陡傾角地層反射波，能夠反映地表高程變化，對于區塊兩端地表高程較高或近地表速度較高的區域的初至時間(圖1中紅色橢圓形標注區域)都需要拾取。

圖1　正演模擬記錄與實際單炮記錄對比Fig.1　The comparison between forward modeling record and single shot records(a) 正演模型；(b) 模擬A炮；(c)模擬B炮 ；(d) 實際A炮；(e)實際B炮

2.1.2替換速度的選擇。

在靜校正處理中，當低降速帶底界水平時，替換速度的大小不會改變反射相位的空間形態，因為替換的量的任意一點是等量的；而當低降速帶底界非水平時，所替換的介質速度與高速層速度不一致，相當于對高速層附加了一層非等厚的另類介質，導致構造形態發生畸變[7]。只有當替換速度等于高速層速度，才能消除低降速帶底界的形態對靜校正量的影響[4]。

研究區塊東部地表被高速老地層覆蓋，利用以往替換速度1 800 m/s不能滿足該區塊的靜校正要求，為了得到更加準確的靜校正量，需要試驗不同基準面和替換速度。模擬實際資料的近地表模型，假設左側地表與高速層頂界之間存在高速體，計算得到的靜校正量變化趨勢如圖2所示，DP表示基準面高程，VR表示替換速度，其中，A到E的曲線分別代表DP=100 m、VR=1 800 m/s；DP=100 m、VR=4 000 m/s；DP=100 m、VR=1 000 m/s；DP=50 m、VR=4 000 m/s；DP=300 m、VR=4 000 m/s的靜校正量曲線。通過分析可知，由于整體量太大，D和E的參數不合理，比較其他三個參數，由于在左側高差大的部位，A和C的平行性不好，且數值偏大，所以最終選擇基準面100 m，替換速度為4 000 m/s。

2.2應用效果分析

不同替換速度的疊加剖面靜校正效果對比：選取研究區塊具有代表性，并能整體控制該區塊低降速帶厚度變化的inline線進行質量監控，從疊加剖面靜校正應用效果來看，圖3(b)比圖3(a)中淺層信息得到明顯改善，基本消除近地表高程及低降速帶帶來的影響，右下部份位置目的層同相軸連續性增強，成像更清晰。

圖2　起伏地表靜校正量示意圖Fig.2　The statics quantity of relief surface

為了保持本次處理目的層的構造特征，解決由于基準面和替換速度不同而產生的整體時差和淺層成像的問題，采用層析反演靜校正與微測井靜校正結合的方式，分別將兩種方法求得的炮點和檢波點的靜校正量，分離成低頻分量和高頻分量，應用層析反演靜校正的高頻分量和微測井靜校正量的低頻分量，在確保靜校正效果的同時避免假構造的產生。從圖4中疊加剖面效果來看，研究區塊內東部地區(剖面右側)出現新的同相軸，說明在淺層老地層的覆蓋下，有新的沉積地層出現。

3　結論

通過對FZ地區地震資料的處理和研究，選定了合適的靜校正處理方法，解決了山前帶及其過度帶靜校正問題，對于這種近地表條件復雜，折射層不穩定，且信噪比較低的資料，初至波層析靜校正方法可以得到相對準確的近地表速度模型，從而計算出更加精確的靜校正量。其中，初至波拾取質量以及替換速度的選取起到了非常關鍵的作用。依據正演模擬結果，單炮初至波時間拾取時，在地表高程較高或近地表速度較高的區域，要拾取有效的遠偏移距異常初至，在多次折射區域，以拾取初至較強層為主。研究區低降速帶厚度和速度橫向變化劇烈，部分地表被高速老地層覆蓋，折射波速度比較高，通過選取不同的替換速度進行對比試驗，選取更合適的替換速度，取得了較好的靜校正處理效果。

圖4　微測井模型靜校正與層析反演靜校正應用效果對比Fig.4　The application effect comparison of micro logging model static correction and tomography inversion static correction(a)微測井模型靜校正；(b)微測井低頻和層析反演高頻組合靜校正

[1]侯波,楊曉倩,況晏.層析靜校正技術在南方頁巖氣處理中的應用[J].山東理工大學學報：自然科學版,2014,28(2):30-33.

HOU B, YANG X Q, KUANG Y. The application of the tomographic inversion static correction[J]. Journal of Shandong University of Technology:Natural Science Edition, 2014,28(2):30-33.(In Chinese)

[2]蘇世龍,蔡希玲,曾慶芹,等.層析靜校正方法及其在東部勘探中的應用[J].油氣地球物理,2010,8(1):11-16.

SU S L, CAI X L, ZENG Q Q, et al. Tomographic static correction method and its application in the eastern exploration[J]. Petroleum Geophysics, 2010,8(1):11-16.(In Chinese)

[3]楊為民,張云崗,劉原英.層析靜校正技術在山地復雜地區三維地震勘探中的應用[J].中國煤田地質,2007,19(2):63-65.

YANG W M, ZHANG Y G, LIU Y Y. Application of tomographic statics in complicated mountainous area 3D seismic prospecting[J]. Coal Geology of China, 2007,19(2):63-65. (In Chinese)

[4]李錄明,羅省賢. 復雜三維表層模型層析反演與靜校正[J]. 石油地球物理勘探,2003,38(6):636-641.

LI L M, LUO S X. Tomographic inversion for 3-D complex near-surface model and static correction[J]. Oil Geophysical Prospecting, 2003,38(6):636-641. (In Chinese)

[5]秦自耕,李子軒, 庹先國,等. 層析建模技術在庫車坳陷東秋山地區靜校正應用探索[J].地球物理學進展,2015,30(2):769-774.

QIN Z G, LI Z X, TUO X G, et al. Application of chromatography modeling techniques in the Kuqa depression statics exploration area east Akiyama[J]. Progress in Geophysics, 2015,30(2):769-774. (In Chinese)

[6]馮澤元,李培民,唐海忠,等. 利用層析反演技術解決山地復雜區靜校正問題[J].石油物探,2005,4(3):454-457.

FENG Z Y, LI P M, TANG H Z, et al. Solving the static correction problem in mountain complex block using tomographic inversion[J]. Geophysical Prospecting for Petroleum, 2005,4(3):454-457. (In Chinese)

[7]劉伊克,常旭,王輝,等.三維復雜地形近地表速度估算及地震層析靜校正[J].地球物理學報,2001,44(2):272-278.

LIU Y K, CHANG X, WANG H, et al. Estimation of near-surface velocity and seismic tomographic static corrections[J]. Chinese Journal of Geophysics, 2001,44(2):272-278. (In Chinese)

[8]鄭鴻明,呂煥通,婁兵,等.地震勘探近地表異常校正[M].北京:石油工業出版社,2009.

ZHENG H M, LV H T, LOU B, et al. Near surface anomaly correction for seismic exploration[M]. Beijing: Petroleum Industry Press, 2009. (In Chinese)

The application of tomographic static correction for three-dimensional seismic data processing in FZ region

LI Na

(Exploration and Development Research Institute of Daqing Oilfield Co. Ltd., Daqing163712, China)

The complex near surface conditions of FZ, results in single shot record changeable and having serious static correction serious problems. It can be known that refraction static correction, curved surface static correction, micro logging static correction and elevation static correction can't meet the processing demands of static correction in target area, combining with the experience of processing in the past. Therefore, it is necessary to use the tomographic inversion static correction method to acquire the near surface structure mode that obtain more accurate statics, and ensure the accuracy of structure image. This paper briefly introduces the basic principle of tomographic static correction method, introducing the near surface conditions and the characteristics of single shot record of FZ. We research the influence on the static correction effect of different replacement velocity and compare them, aiming at the particularity of FZ.Thus prove that the first-break picking quality and the replace velocity have played a very key role.

tomographic inversion; statics correction; complex surface; replacement velocity; first break

2015-12-18改回日期：2016-02-26

李娜(1989-)，女，本科，從事地震資料處理工作，E-mail：lnlnln2013@petrochina.com.cn。

1001-1749(2016)04-0525-05

P 631.4

A

10.3969/j.issn.1001-1749.2016.04.14