成 林

(西安石油大學(xué) 地球科學(xué)與工程學(xué)院，西安 710065)

地面地震高斯射線束正演

成 林

(西安石油大學(xué) 地球科學(xué)與工程學(xué)院，西安 710065)

高斯射線束正演方法是將射線理論和波動方程理論相結(jié)合的一種地震正演方法，充分利用了射線理論計算速度快、波動方程理論計算精度高的優(yōu)點.通過高斯射線束表達式和編程實現(xiàn)地震記錄合成的基礎(chǔ)上，給出了水平地質(zhì)模型的計算實例，并對比分析了高斯射線束法、普通射線法、波動方程法正演的特點.通過理論分析和實際地質(zhì)模型可以得出：(1)高斯射線束正演方法效率高，計算速度快，在復(fù)雜模型中甚至超過普通射線法的計算速度；(2)高斯束法由于在斷點處存在繞射能量，計算中不存在奇異區(qū)(陰影區(qū)、焦散區(qū)等)；(3)高斯束法合成記錄不僅反映地震波的運動學(xué)特點，而且反映波的高頻動力學(xué)特點，計算結(jié)果可與波動方程法計算結(jié)果一致.

高斯射線束；射線追蹤；波動方程；地震正演；合成地震記錄

在石油勘探中，地震勘探一直是很重要的方法之一.地震勘探有反演和正演兩種方法.地震反演是利用地震反射波、繞射波等波場信息來模擬地層的地質(zhì)構(gòu)造的方法.但是地震反演具有不穩(wěn)定性和結(jié)果多解性問題.為了解決這些問題，一般用已有的地質(zhì)信息和理論知識結(jié)合來建立正確的地質(zhì)模型，通過迭代法，逐漸逼近真實解.所以正演技術(shù)在地震勘探中越來越受到關(guān)注和重視[1-5].

正演方法在地震勘探中的發(fā)展和應(yīng)用，確定了波動理論和射線理論兩種地震技術(shù).波動理論和射線理論可以將橫向不均勻介質(zhì)中的高頻波記錄合成，但是對焦散區(qū)、臨界區(qū)等奇異區(qū)域計算出的振幅趨于無窮，因而失效.因此在1982年Cerveny等人又提出了高斯射線束正演方法，這種方法結(jié)合了波的動力學(xué)特征和運動學(xué)特征，不僅適用于均勻介質(zhì)，而且在復(fù)雜的非均質(zhì)介質(zhì)中同樣有著較好的效果.該方法具有速度快、精度高的特點.

1 高斯射線束正演方法理論

1.1 二維高斯射線束表達式

根據(jù)Cerveny等人在二十世紀(jì)八十年代初期提出高斯射線束正演方法，可以得知二維高斯射線束表達式為公式(1)，集中于射線附近的高頻漸進時間調(diào)和解，該解在頻率域具有如式(2)中的表達形式.[6]

(1)

在公式(1)中：μ1表示位移分量；τ11表示應(yīng)力分量(其中包含拉梅常數(shù)λ和μ)；h1，h2和h3分別為坐標(biāo)c1c2和c3對應(yīng)的標(biāo)量因子，且h1=1.這里只研究c2=0的平面.

(2)

圖1 二維高斯射線坐標(biāo)系

公式(2)中各項的意義：

U(s,n,ω,t)表示了縱波位移，(s,n)為沿某一中心射線的射線坐標(biāo)系坐標(biāo)，s為到射線原點的距離，n為到射線的距離，ω為圓頻率即ω=2·π·f，t為時間參量；

K(s)=v(s)Re[p(s)/q(s)]

表示高斯射線的相前曲率，主要反映波前面的幾何形狀；其中p(s)和q(s)為地震波沿中心射線變化的兩復(fù)函數(shù)，并且他們滿足如下微分方程：

(3)

L(s)表示高斯射線束有效寬度的一半；

A(s)是沿中心射線地震波的位移幅值.對于層狀介質(zhì)，假定射線從原點s0出發(fā)，經(jīng)N個任意界面反射回到地面接收點R，則R處接收到的中心射線振幅可以表示為：

(4)

圖2 波在界面上的反射與透射示意圖

圖3 垂直射線剖面中的振幅曲線

圖4 檢波點R附近的一條射線Ω

式(4)中s為射線路徑上的任意點，Ri為界面反射或透射系數(shù)，ρ為入射線一側(cè)的介質(zhì)密度，v為入射線一側(cè)的地震波速，而ρ′和v′則分別表示生成射線一側(cè)的密度與速度，αi和βi分別表示以地層界面切線方向為水平坐標(biāo)軸的局部坐標(biāo)系中，入射角和生成線與局部坐標(biāo)系X正向的夾角，如圖2.

式中，ε=ε1+iε2.通過改變ε1和ε2，可以將射線的半寬度進行調(diào)節(jié)，確保有較高的計算精度.由方程(2)可見，從射線中心開始波場是以指數(shù)規(guī)律逐漸向外衰減的.如圖3所示，在垂直射線的剖面內(nèi)振幅為鐘形曲線，這就是高斯射線束名稱的由來.

1.2 高斯射線束合成地震記錄

根據(jù)高斯射線束理論計算地震波場時，可分為3步：(1)做射線追蹤；(2)作動力學(xué)追蹤;(3)求得波場值.上述介紹如何在頻率域中計算一條中心射線的高斯射線束.高斯射線束合成地震記錄方法有：頻譜法、褶積法和波包法，其中波包法效率最高.[7-8]

例如未得到某介質(zhì)中R點的波場值，需利用高斯射線束正演將R點附近的波場貢獻疊加，然后轉(zhuǎn)換到時間域.設(shè)射線出射點為O，檢波點為R，并位于射線束的有效半寬度范圍內(nèi)，則R點的振幅A(R),旅行時τ(R)，可由O點的振幅A(O)及旅行時τ(O)近似求得.如圖4所示，Q點為R點的波前與射線Ω的交點，只要知道sQ_sO(即Q點到O點的射線弧長)便近似求得A(R)和τ(R)的值.

射線Ω對檢波點R處波場的貢獻為

(5)

而檢波點R處的總波場就是對所有射線求和的結(jié)果，即

(6)

為了求得合成記錄，我們引入震源函數(shù)：

f(t)=exp[(2πfmt/γ)]2cos(2πfm+ν)

(7)

現(xiàn)在用u(R,t)表示時間域的實際波場，可由方程(7)的傅里葉變換求得：

(8)

(9)

(10)

g為波包，由于f(t)是高頻函數(shù)，所以g為高頻的波包.

f(t)≈exp[(2πfmt/γ)]2cos(2πfm+ν)

(11)

式中fm為子波主頻，γ為控制子波包絡(luò)寬度的參數(shù)，ν為控制子波相位的參數(shù)，此子波也稱Puzyrev子波或高斯包絡(luò)子波.根據(jù)Cerveny 1983的文章，可知波包g的近似解析表達式：

(12)

1.3 二維均勻介質(zhì)線性界面高斯射線束合成記錄p、q值的求取

經(jīng)過運動學(xué)追蹤，就有了中心射線，在此基礎(chǔ)上可以進行高斯射線束的動力學(xué)射線追蹤，即求得p、q的函數(shù)值.函數(shù)p(s)、q(s)在高斯射線束中起著非常重要的作用，它們決定了高斯射線束能量的分布狀態(tài)，也表征沿射線傳播方向的高頻地震波場動力學(xué)特征，也正是因為這一點，才使得高斯射線束正演方法優(yōu)于一般的射線正演方法.主要還是求檢波點處p、q的值.下面進行p、q值的求取在二維均勻介質(zhì)同一個層內(nèi)中，速度是不隨位置發(fā)生變化，則

(13)

故式(13)可以化為：

(14)

解此一階常微分方程得：

(15)

式中C1和C2為任意常數(shù)，選取兩組相互獨立的初始條件

可求得：

(16)

式中v(1)為炮點處速度值；另外，由于在線性邊界上，界面曲率為0，結(jié)合邊界條件式(16)可知檢波點處p、q值為：

(17)

式中v(i)為第i段射線傳播的速度，τ(i)為相應(yīng)的時間，M為射線上的總段數(shù).

由式(17)可知：

(18)

故有

(19)

2 高斯射線束算法的實現(xiàn)及實例分析

實際計算表明高斯射線束法數(shù)值模擬是有效的，根本不需要在兩點間作射線追蹤，從而使得運動學(xué)追蹤部分的計算速度比普通射線法要快.利用高斯射線束法合成地震記錄可以明顯看出地震波的動力學(xué)特征，使得計算精度比普通射線法要高，而且與波動方程效果相差不大，但計算速度比波動方程法快得多.[9]下面我們通過水平地質(zhì)模型實例來具體分析高斯射線束正演方法的特點.

該模型為一個兩層的水平層狀模型，表1為這兩層介質(zhì)速度、密度、厚度參數(shù)和觀測系統(tǒng)參數(shù)，圖5為此模型的運動學(xué)追蹤圖，圖6為利用高斯射線束法合成的炮集記錄，圖7為利用試射法合成的炮集記錄，圖8為利用波動方程法合成的炮集記錄.在地震勘探中，水平層狀模型是比較理想化且簡單的一種模型，且在構(gòu)造簡單的地區(qū)的地下地質(zhì)構(gòu)造與這個模型基本類似，因此用這個模型可簡單地試出程序的有效性和在檢驗高斯射線束的優(yōu)點時，可以明確表達出.

表1 水平層狀模型觀測參數(shù)及地層參數(shù)

從該模型的運動學(xué)追蹤圖5可以看出：射線從-90°到90°打出，±90°打出的射線沿地面?zhèn)鞑?，形成直達波(圖8地表附近的水平紅線)；其他角度的射線向下傳播，遇到第一個界面發(fā)生了反射和透射，反射波回到地表被檢波器接收到，形成第一層的反射波；透射波繼續(xù)向下傳播，遇到第二個界面再次發(fā)生反射和透射，反射波傳播到地表就形成了第二層的反射波，透射波繼續(xù)向下傳播，由于沒遇到反射界面，此部分能量就損失掉了.從地表接收情況來看，水平層狀地層射線追蹤不存在射線覆蓋不到的陰影區(qū)域，所以對于此種簡單模型用普通的射線方法也可以獲得良好的正演效果，如圖7，可以進行構(gòu)造方面的正演研究.

圖5 水平層狀模型運動學(xué)追蹤圖

圖6 水平層狀高斯射線束法炮集記錄

圖7 水平層狀試射法炮集記錄

圖8 波動方程法炮集記錄

對比圖6、圖7和圖8，可以看出炮集記錄有以下特點：

(1)直達波的時距曲線(3個炮集記錄中同相軸①)為過炮點的直線，反射波的時距曲線為雙曲線，并且在同一地震道上直達波比反射波到達時間早；對于同一條反射波同相軸，零炮檢距處的地震道上旅行時最小，隨炮檢距的增大，旅行時逐漸增大；對比兩條反射波同相軸可知，反射界面越深，反射波時距曲線的曲率越小.這是由于深度越大，反射界面以上地層的均方根速度越大，相應(yīng)的正常時差就越小，因此同相軸就越平緩，曲率就越小.

(2)高斯射線束法合成的地震記錄體現(xiàn)了一定的地震波動力學(xué)特征，和波動方程法合成記錄相似，如圖6和7，各波能量有一定差異；對比兩個反射界面的振幅可知，反射同相軸的能量隨界面深度的增大而減小，這主要是因為界面越深，傳播的射線路徑越長，衰減的能量越多，振幅就越??；對于同一同相軸，零炮檢距附近能量較強，隨炮檢距增大反射能量逐漸減弱，這也是由于炮檢距越大，反射路徑越長，能量衰減越厲害造成的.由此可以看出高斯射線束正演方法不僅考慮了不同深度上能量的變化，還考慮了不同入射角射線能量的差異，所以此方法可以在一定程度上體現(xiàn)地震波的動力學(xué)特征.而用試射法得到的合成記錄，如圖7，各波同相軸的振幅一樣，不能體現(xiàn)地震波傳播的動力學(xué)特征，這是普通射線法的一大缺點.

3 結(jié)論

通過計算機，利用C++語言編寫的高斯射線束程序，對水平層狀地質(zhì)模型進行運動學(xué)追蹤和地震記錄的合成，并把其結(jié)果與試射法、波動方程法計算合成記錄相對比得出下列兩個結(jié)論

(1)高斯束法結(jié)合了普通射線正演和波動力學(xué)正演兩種方法的優(yōu)點，同時考慮了波場的運動學(xué)和動力學(xué)特征，適用于復(fù)雜的非均勻介質(zhì)模型，還能考慮介質(zhì)的吸收作用，并且不需要兩點射線追蹤的計算，具有速度快、精度較高的特點，對焦散區(qū)、臨界區(qū)及暗區(qū)等奇異區(qū)域都具有較好的計算效果.

(2)在應(yīng)用方面，可利用高斯射線束地震記錄合成技術(shù)，可對波場進行快速的理論研究，進行野外數(shù)據(jù)采集工作的設(shè)計，而且還可作為野外數(shù)據(jù)采集現(xiàn)場質(zhì)量控制的依據(jù)，也有助于地震資料的處理和解釋.

[1] 李輝峰，徐峰.地震勘探新技術(shù)[M].北京：石油工業(yè)出版社,2009.

[2] 鄧飛，劉超穎.三維射線束快速追蹤及高斯射線束正演[J].石油地球物理勘探，2009，44(2)：158-165.

[3] 江玉樂，雷宛.地球物理信號與信息處理教程[M].北京：地質(zhì)出版社，2006：80-89.

[4] 劉學(xué)才,周熙襄,沙椿.三維高斯射線束法合成三分量VSP記錄[J].石油地球物理勘探,1995,30(5) :669～680.

[5] 趙改善.快速射線追蹤算法[M].北京：石油地球物理勘探雜志，1991.

[6] 馮德益.地震波理論應(yīng)用[M].北京：地震出版社，1988.

[7] 鄧飛，劉超穎，趙波，等.高斯射線束正演與偏移[J].石油地球物理勘探，2009,44(3)：265-269,281.

[8] 郭全仕,鄔達理,唐金良,等.井間地震反射波成像技術(shù)探討[J].石油物探,2005,44(5):439-444.

[9] 劉學(xué)才.三維高斯射線束地震模擬[J].石油地球物理勘探，1995，30(5)：669-680.

[責(zé)任編輯 仲 圓]

Gaussian Beam Forward Modeling on the Surface Seismic

CHENG Lin

(School of Earth Sciences and Engineering, Xi’an Shiyou University, Xi’an 710065, China)

The forward modeling by Gaussian beam method is essentially a seismic forward modeling method that combines wave equation and ray theory, which makes full use of the advantages of the fast calculation of ray theory and the high accuracy of wave equation theory. In this paper, based on Gaussian beam function and programming to achieve the synthesis of seismic records, the calculation example of horizontal geological model is given, and the characteristics of Gaussian beam method, ordinary ray method and wave equation method are compared and analyzed. Through the theoretical analysis and the practical geological model test, it can be concluded that as follows. Firstly, forward modeling by Gaussian beam method is of high efficiency, fast calculation speed, and its calculation speed is even faster than that of the ordinary ray method in the complex model. Secondly, due to the presence of diffraction energy at the break-point, the singular region (region of a shadow area, etc.) does not exist in the calculation of Gaussian beam method. Thirdly and lastly, the synthetic record of Gauss beam method not only reflects the kinematic characteristics of seismic waves, but also reflects the high frequency dynamic characteristics of the wave. The calculation results are consistent with the calculation results of the wave equation method.

Gaussian beam; ray tracing; wave equation; seismic forward modeling; synthetic seismic record

1008-5564(2016)06-0086-06

2016-04-25

中國地質(zhì)調(diào)查局地質(zhì)項目(1212011085499)

成 林(1991—)，男，陜西武功人，西安石油大學(xué)地球科學(xué)與工程學(xué)院碩士研究生，主要從事應(yīng)用地球物理研究

P315.2

A