羅騰騰,徐基祥,秦 臻,孫夕平

(1.中國(guó)石油勘探開(kāi)發(fā)研究院,北京100083;2.三峽大學(xué),湖北宜昌443002)

常規(guī)的偏移算法主要記錄由連續(xù)反射層或主要不連續(xù)點(diǎn)(大斷層)定義的“高能事件”,它是地震解釋人員通常使用的地震信息,常被稱作“鏡像”能量。隨著油氣勘探開(kāi)發(fā)研究的深入,勘探目標(biāo)逐步轉(zhuǎn)向斷層、裂縫、河道、粗糙巖丘邊緣等地下小尺度不連續(xù)地質(zhì)體[1],這些小尺度地質(zhì)體的地震響應(yīng)通常表現(xiàn)為能量較弱的繞射波。

20世紀(jì)50年代,KREY[2]首先利用繞射波信息來(lái)研究地下地質(zhì)異常體如斷層、斷點(diǎn)等。KLEM-MUSATOV等[3]指出,繞射波振幅通常比反射波振幅弱一到兩個(gè)數(shù)量級(jí),且在常規(guī)處理流程中多以反射波作為有效信號(hào),而那些高分辨率、低能量的繞射波信號(hào)常被破壞或壓制。為了充分利用地震記錄中的繞射波信息來(lái)刻畫(huà)地下非均質(zhì)地質(zhì)體,需要在全波場(chǎng)記錄中將繞射波場(chǎng)分離并單獨(dú)成像[4]?；诜瓷洳ê屠@射波在運(yùn)動(dòng)學(xué)和動(dòng)力學(xué)方面的特征差異[5-6],國(guó)內(nèi)外學(xué)者提出了許多不同的繞射波分離成像方法。按照繞射波分離成像的方法原理以及處理手段的差異可將現(xiàn)有方法分為:疊加、濾波、平面波解構(gòu)濾波、聚焦-反聚焦、Radon變換-反Radon變換等。

1) 疊加方法:LANDA等[7]、KANASEWICH等[8]、TSINGAS等[9]提出在共偏移距道集和共斷層點(diǎn)道集上使用繞射波走時(shí)曲線來(lái)描述繞射波相干的總和。LANDA等[10]利用在共繞射點(diǎn)剖面上疊加來(lái)自繞射點(diǎn)處的地震信號(hào)識(shí)別局部非均質(zhì)地質(zhì)體。

2) 濾波方法:KOZLOV等[11]進(jìn)一步修改Kirchhoff偏移算子中的錐形權(quán)函數(shù)使其更好地壓制反射能量,對(duì)地震響應(yīng)中的弱繞射波信息進(jìn)行準(zhǔn)確成像;BANSAL等[12]討論了多種分離反射和繞射波場(chǎng)的方法,主要包括共偏移距傾角濾波、正常時(shí)差校正傾角濾波、正常時(shí)差-傾角時(shí)差校正傾角濾波、特征向量濾波、Radon濾波等,并由合成記錄闡述了不同分離方法的優(yōu)缺點(diǎn);MOSER等[13]在Kirchhoff偏移積分公式中通過(guò)修改加權(quán)函數(shù)來(lái)實(shí)現(xiàn)反穩(wěn)相繞射波提取;KOREN等[14]提出一種繞射波成像新方法,即利用方向角分解設(shè)計(jì)加權(quán)疊加濾波器實(shí)現(xiàn)各向同性/各向異性介質(zhì)中繞射波和反射波的分離,同時(shí)提高連續(xù)構(gòu)造和小尺度地質(zhì)體的成像分辨率;李曉峰等[15]在傳統(tǒng)的Kirchhoff偏移算子中引入反穩(wěn)相濾波器來(lái)壓制反射波能量,凸顯繞射波能量,并在成像過(guò)程中引入校正極性算子實(shí)現(xiàn)繞射波準(zhǔn)確成像;劉培君等[16]借助相關(guān)分析以及高斯束偏移的射線追蹤構(gòu)建反穩(wěn)相濾波算子,采用反穩(wěn)相偏移得到繞射波成像剖面。

3) 平面波解構(gòu)濾波方法:CLAERBOUT[17]最先提出平面波解構(gòu)濾波器的概念,并于1994年[18]進(jìn)一步完善了該濾波器方法;FOMEL[19]改進(jìn)了平面波解構(gòu)濾波器,使所需唯一參數(shù)為局部平面波場(chǎng)的斜率;TANER等[20]模擬了平面波剖面和平面波分解濾波器,并依據(jù)反射波和繞射波在平面波地震記錄上的同相軸連續(xù)性差異在疊前道集中進(jìn)行繞射波分離成像;黃建平等[21]對(duì)平面波分解(PWD)解構(gòu)濾波器使用方法及其在疊前、疊后道集上分離成像繞射波進(jìn)行了系統(tǒng)總結(jié);孔雪等[22]和劉玉金等[23]證實(shí)平面波剖面本質(zhì)上具有區(qū)分繞射波和反射波特性,并利用平面波解構(gòu)濾波技術(shù)壓制反射波來(lái)獲得繞射波;朱生旺等[24]利用局部?jī)A角濾波技術(shù)和預(yù)測(cè)反演技術(shù)修正繞射波場(chǎng)分離時(shí)低傾角信息估計(jì)不足對(duì)平面波分解濾波器的影響,從而提高繞射波成像結(jié)果的橫向分辨率;DECKER等[25]在疊前偏移傾角道集上對(duì)每個(gè)固定傾角的成像剖面應(yīng)用PWD解構(gòu)濾波技術(shù)并疊加繞射波場(chǎng)獲得單獨(dú)的繞射波成像剖面。

4) 聚焦-反聚焦方法:KHAIDUKOV等[26]在疊前共炮點(diǎn)道集將反射波聚焦于鏡像虛震源點(diǎn),而繞射波仍保持發(fā)散的狀態(tài),利用聚焦切除-反聚焦方法成功提取斷層、斷點(diǎn)等小尺度地質(zhì)體,并對(duì)其成像;MOSER等[13]提出兩種繞射波分離成像的方法,一種方法是在疊前深度域利用聚焦濾除反射波實(shí)現(xiàn)繞射波分離與成像,另一種方法是修改偏移算法的核函數(shù),采用反穩(wěn)相偏移壓制反射能量;BERKOVITCH等[27]利用繞射多次聚焦疊加的方法將選取的繞射波最優(yōu)同相軸疊加從而得到主要包含繞射波的疊加剖面;RESHEF等[28]利用角道集上的繞射能量進(jìn)行高分辨率速度分析。

5) Radon變換-反Radon變換方法:KLOKOV等[29]根據(jù)反射波同相軸在傾角域CIGs中表現(xiàn)為開(kāi)口向上的“笑臉”狀,利用混合域Radon變換識(shí)別反射波頂點(diǎn)并消除反射能量,增強(qiáng)繞射信號(hào);KLOKOV等[30]推導(dǎo)了偏移傾角域反射波和繞射波的解析表達(dá)式,并利用混合域Radon變換實(shí)現(xiàn)繞射波場(chǎng)和反射波場(chǎng)分離。該方法計(jì)算成本低,但是利用常規(guī)分辨率Radon變換分離的繞射波場(chǎng)中會(huì)存在大量殘余反射波,國(guó)內(nèi)外學(xué)者致力于提高Radon變換的分辨率來(lái)改善波場(chǎng)分離效果。

Radon變換(Radon Transform,RT)首先由奧地利著名數(shù)學(xué)家J.RADON在1917年提出。20世紀(jì)70年代,美國(guó)斯坦福大學(xué)以CLAERBOUT為代表的地球物理小組首次將Radon變換應(yīng)用于地震勘探中,并對(duì)線性Radon變換(τ-p變換)進(jìn)行了重點(diǎn)研究。THORSON等[31]提出的隨機(jī)反演方法提高了時(shí)間域Radon變換的分辨率;HAMPSON[32]提出頻率域最小平方拋物線Radon變換,該方法因其高效的計(jì)算效率迅速成為當(dāng)時(shí)工業(yè)界的標(biāo)準(zhǔn);SCALES等[33]將迭代重加權(quán)最小二乘算法與預(yù)條件共軛梯度法相結(jié)合,有效求解了線性方程組中的大型稀疏矩陣;SACCHI等[34]提出高精度Radon變換,即以有效信號(hào)的稀疏性作為約束條件利用迭代高分辨率算法求解;CARY[35]指出時(shí)間域Radon變換結(jié)果具有更強(qiáng)的稀疏性而頻率域Radon變換能更準(zhǔn)確地刻畫(huà)有效信號(hào)的頻率特征;HERRMANN[36]提出去假頻高分辨率拋物線Radon變換(DHR),該方法利用低頻數(shù)據(jù)的無(wú)假頻模型解構(gòu)造加權(quán)算子來(lái)壓制高頻數(shù)據(jù)Radon空間的假頻能量,獲得了較高的信噪比和分辨率。

在國(guó)內(nèi),同樣有許多學(xué)者致力于Radon變換研究。吳律[37]全面分析討論了τ-p變換的原理和應(yīng)用;牛濱華等[38]首次提出多項(xiàng)式Radon變換。許多國(guó)內(nèi)學(xué)者[39-44]在提高Radon變換分辨率方面做了大量系統(tǒng)深入的研究工作。

為了解決Radon變換能量團(tuán)泄漏,繞射波場(chǎng)分離不徹底的問(wèn)題,本文在時(shí)間-頻率混合域,引入預(yù)條件算子并構(gòu)建了新的時(shí)變稀疏模型權(quán),在傾角域CIGs上發(fā)展了一種基于預(yù)條件共軛梯度(PCG)算法的混合域高分辨率Radon變換繞射波分離成像方法。文中對(duì)比分析了本文方法與常規(guī)分辨率Radon變換、頻率域高分辨率Radon變換在繞射波分離中的應(yīng)用效果。

1 方法原理

1.1 傾角域共成像點(diǎn)道集波場(chǎng)特征分析

常速介質(zhì)疊后偏移反射面的深度z(x)可以表示為(圖1)[30]:

z(x)=z0+xtanα0

(1)

式中:z0表示傾斜反射界面與z軸交點(diǎn)的深度;α0為反射界面的傾角。在地表y處的自激自收地震響應(yīng)可表示為:

(2)

式中:v表示介質(zhì)速度,t為旅行時(shí)間。

圖1 零偏移距反射示意

由模型坐標(biāo){x,z}和數(shù)據(jù)坐標(biāo){y,t}之間的映射關(guān)系可以得到偏移公式:

(3)

(4)

式中:vm是偏移速度;α為偏移傾角。將公式(2)代入公式(3)和公式(4),并消去y得到傾斜反射界面在傾角域的成像表達(dá)式:

(5)

對(duì)(5)式求偏導(dǎo)可以得到:

(6)

從公式(6)可以看出,當(dāng)偏移速度準(zhǔn)確時(shí),即vm=v且α=α0時(shí),(6)式的偏導(dǎo)數(shù)值為0,即α0為(5)式的一個(gè)極小值點(diǎn)。因此,對(duì)于準(zhǔn)確的偏移速度,在傾角域CIGs上,反射響應(yīng)是一個(gè)具有穩(wěn)相頂點(diǎn)的“笑臉”狀曲線,且穩(wěn)相頂點(diǎn)的位置與地層傾角相對(duì)應(yīng)。

類似地,可以得到地下繞射點(diǎn)(x0,z0)處的傾角域響應(yīng)表達(dá)式:

(7)

zα(α)=z0

(8)

由(8)式可以看出繞射響應(yīng)為一水平直線。

下面用3個(gè)理論模型來(lái)直觀闡明傾角域道集曲線的特點(diǎn),并進(jìn)一步分析反射波和繞射波的差異。圖2a 至圖2c顯示了3個(gè)理論模型,均包含一個(gè)反射界面和一個(gè)繞射體,繞射體位于地下1km處。圖2a中水平反射界面位于地下0.5km處,圖2b和圖2c中反射界面分別向左傾斜和向右傾斜。我們從理論模型0.6km位置處(繞射點(diǎn)正上方)提取傾角域CIGs(圖2d至圖2f),圖2d至圖2f中虛線表示反射界面(水平和30°傾斜反射界面);實(shí)線表示繞射點(diǎn)處的傾角域響應(yīng)曲線。

從圖2d至圖2f可以看出,對(duì)于反射界面來(lái)說(shuō),不論界面是否傾斜,在傾角域道集中響應(yīng)曲線總是保持“微笑”的形狀(如圖中虛線所示),且曲線的穩(wěn)相頂點(diǎn)剛好指示界面的傾角;而在繞射點(diǎn)位置處的繞射響應(yīng)為水平直線(如圖中實(shí)線所示)。

綜上可知,在傾角域CIGs中,反射與繞射響應(yīng)曲線存在明顯的幾何形態(tài)差異,其中反射響應(yīng)為具有穩(wěn)相頂點(diǎn)的凹形曲線,繞射響應(yīng)為水平直線。我們可以利用兩者之間形狀差異將繞射波同相軸單獨(dú)分離。下面以此原理我們采用拋物線Radon變換來(lái)分離傾角域CIGs中反射波和繞射波。

圖2 理論模型及偏移速度正確時(shí)的傾角域共成像點(diǎn)道集示意a,b,c 包含一個(gè)反射界面和一個(gè)繞射點(diǎn)的理論模型; d a 在0.6km處的傾角域CIGs; e b 在0.6km處的傾角域CIGs; f c 在0.6km處的傾角域CIGs

1.2 Radon變換原理

二維Radon變換的定義式為:

(9)

式中:d(t,x)為時(shí)間域地震數(shù)據(jù);m(τ,q)為Radon域數(shù)據(jù);t,τ分別為時(shí)空域的時(shí)間和Radon域的截距時(shí)間;q表示曲率參數(shù);x為炮檢距。將(9)式進(jìn)行離散傅里葉變換,轉(zhuǎn)換到頻率域,可得:

(10)

式中:D(ω,x)和M(ω,q)分別為d(t,x)和m(τ,q)的傅里葉變換結(jié)果;ω為頻率;n為Radon域的曲率個(gè)數(shù)。由(10)式可知,對(duì)于某一頻率ω,Radon變換可以表示為一個(gè)線性方程組的求解問(wèn)題:

D=LM

(11)

式中:L是Radon變換算子,由地震數(shù)據(jù)道的相對(duì)空間位置x0,x1,…,xm-1和Radon變換參數(shù)q0,q1,…,qn-1來(lái)確定,可表示為:

(12)

(13)

(14)

常規(guī)Radon變換采用最小二乘反演計(jì)算,定義目標(biāo)函數(shù):

(15)

對(duì)上述目標(biāo)函數(shù)求極小值,可以得到M的常規(guī)最小二乘解,即:

M=(LTL+βI)-1LTD

(16)

1.3 基于預(yù)條件共軛梯度(PCG)迭代的混合域高分辨率Radon變換

為了提高Radon變換的分辨率,在模型中引入稀疏約束條件,將目標(biāo)函數(shù)改寫(xiě)為:

(17)

式中:R(M)表示稀疏約束項(xiàng)?？紤]到模型空間的噪聲對(duì)Radon變換分辨率的影響比實(shí)際地震數(shù)據(jù)中的非線性噪聲影響更為嚴(yán)重,因此僅引入模型加權(quán)矩陣Wm來(lái)提高模型的擬合度,如分辨率和平滑度。求解方程

(18)

(19)

(20)

令A(yù)=F-1LTLF,b=F-1LTFD,則方程(20)可以寫(xiě)為:

(21)

可得殘差量表達(dá)式為:

(22)

常規(guī)Radon變換采用共軛梯度(CG)算法求解,高分辨率Radon變換采用PCG反演方法,本文基于PCG算法的高分辨率Radon變換處理步驟如下。

1) 由方程(11)計(jì)算得到頻率域常規(guī)分辨率解,經(jīng)傅里葉變換由頻率域變換到時(shí)間域,作為PCG高分辨率迭代的初始解m0。

3)內(nèi)循環(huán)迭代。用PCG法[45]迭代求解方程(20),其中關(guān)鍵求解步驟如下。

② 求解迭代步長(zhǎng)。αj:=(rj,zj)/(Apj,pj),其中,j=0,1,…,nmax為內(nèi)部迭代次數(shù),nmax表示設(shè)置的最大迭代次數(shù)。

③ 計(jì)算新的迭代解。mj+1:=mj+αjpj。

④ 計(jì)算更新殘差。rj+1:=rj-αjApj。

⑥ 求解共軛方向的搜索步長(zhǎng)。βj:=(rj+1,zj+1)/(rj,zj)。

⑦ 計(jì)算下一次迭代的共軛向量。pj+1:=zj+1+βjpj。

當(dāng)多次迭代后達(dá)到事先設(shè)置好的誤差極限時(shí),即可得到(τ,q)的解。

4) 繼續(xù)進(jìn)行外部迭代,求取滿意的高分辨率解。利用步驟3)中的(τ,q)解,重新計(jì)算稀疏約束矩陣,重復(fù)步驟2)和步驟3),繼續(xù)迭代。外部迭代的主要作用是更新稀疏權(quán),增加(τ,q)解的稀疏性。

5) 進(jìn)行Radon反變換,完成反射波切除,得到繞射波分離結(jié)果。

2 數(shù)據(jù)測(cè)試

應(yīng)用的技術(shù)流程主要包括:對(duì)傾角域CIGs進(jìn)行Radon變換,在Radon域?qū)㈤_(kāi)口向上的反射波聚焦于曲率小于0的位置,而繞射波的曲率在零值附近。在進(jìn)行反射波切除時(shí),選擇的Radon域切除的q值在0值附近,以減少反射波能量的干擾并盡可能完整地保留繞射波信號(hào)。最后采用反Radon變換將分離出的繞射波同相軸恢復(fù)為時(shí)空域數(shù)據(jù),即從全波場(chǎng)中分離出繞射波。

2.1 模擬數(shù)據(jù)測(cè)試

在模擬的傾角域CIGs中,包含兩個(gè)反射波和一個(gè)繞射波響應(yīng)曲線,反射波頂點(diǎn)分別位于1000m、1500m,繞射波則處于地下2000m,如圖3所示。我們對(duì)比分析了最小二乘Radon變換方法、頻率域高分辨率Radon變換方法和本文方法在Radon域中對(duì)繞射波的分離效果。采用最小二乘Radon變換對(duì)傾角域CIGs進(jìn)行繞射波分離,其Radon變換結(jié)果以及繞射波分離后的結(jié)果如圖4所示,基于該方法所得到的Radon域結(jié)果具有明顯的剪刀狀發(fā)散,這些有限孔徑產(chǎn)生的端點(diǎn)假象嚴(yán)重影響了Radon變換的信噪比,進(jìn)而導(dǎo)致在繞射波分離結(jié)果中存在大量殘余的反射能量(圖4b中紅色箭頭所示)。

圖3 模擬數(shù)據(jù)的CIGs記錄

采用頻率域高分辨率Radon變換對(duì)模擬CIGs記錄進(jìn)行繞射波分離,其結(jié)果如圖5所示。從圖5可以看出,頻率域高分辨率Radon變換結(jié)果相較于最小二乘Radon變換結(jié)果,剪刀狀發(fā)散能量減弱,分辨率有了明顯提高,但在繞射波上方的假象仍然存在。我們將內(nèi)循環(huán)迭代次數(shù)分別選取為10次、20次、40次,譜白化因子分別選取為0.01、0.05、0.10,該假象依然存在。該方法得到的繞射波分離后的結(jié)果中同樣存在反射波的殘余。

圖4 最小二乘Radon變換對(duì)繞射波的分離結(jié)果a Radon域結(jié)果; b 繞射波分離后的結(jié)果

圖5 頻率域高分辨率Radon變換對(duì)繞射波的分離結(jié)果a Radon域結(jié)果; b 繞射波分離后的結(jié)果

采用本文方法對(duì)該模擬數(shù)據(jù)進(jìn)行繞射波分離,結(jié)果如圖6所示。從圖6可以看出,基于混合域高分辨率Radon變換的子波幅度無(wú)損失且分辨率很高,繞射波分離后的結(jié)果也優(yōu)于其它兩種方法。

繞射波場(chǎng)分離效果的好壞直接影響繞射波成像效果,對(duì)比3種方法的分離結(jié)果可見(jiàn),基于混合域高分辨率Radon變換的繞射目標(biāo)成像方法相較于其它兩種方法得到的結(jié)果分辨率更高,子波保持得更好,這將有利于后續(xù)小尺度地質(zhì)異常體的準(zhǔn)確識(shí)別和定位。

2.2 實(shí)際資料測(cè)試

選取中國(guó)西部M地區(qū)實(shí)際地震資料進(jìn)行測(cè)試,該地區(qū)溶蝕孔洞和裂縫儲(chǔ)層十分發(fā)育,利用傳統(tǒng)的成像處理流程難以對(duì)儲(chǔ)層內(nèi)的多種繞射目標(biāo)體如速度異常體、河道、斷層、生物礁等進(jìn)行精確成像。首先利用本文方法對(duì)某一CDP點(diǎn)下Kirchhoff疊前深度偏移(PSDM)中產(chǎn)生的傾角域CIGs進(jìn)行反射波和繞射波分離。

圖6 混合域高分辨率Radon變換對(duì)繞射波的分離結(jié)果a Radon域結(jié)果; b 繞射波分離后的結(jié)果

圖7a和圖7b分別為實(shí)際地震數(shù)據(jù)的傾角域CIGs和采用本文方法分離后的繞射波數(shù)據(jù)。從圖7可以看出,原始道集中的強(qiáng)反射能量被明顯壓制,由于大部分反射能量被去除,淹沒(méi)在全波場(chǎng)中的弱繞射能量凸顯出來(lái),位于地下約5.5km處的繞射波同相軸能量明顯加強(qiáng),如圖中紅色箭頭所指部位。

進(jìn)一步將本文方法應(yīng)用于該地區(qū)某條測(cè)線的傾角域CIGs,其波場(chǎng)分離前的PSDM剖面如圖8a所示,由于強(qiáng)反射能量的干擾,繞射點(diǎn)、斷層點(diǎn)的位置難以識(shí)別。采用本文方法分離后的繞射波成像剖面如圖8b所示,巖丘邊界刻畫(huà)更加清晰,并且強(qiáng)反射同相軸被濾除,淹沒(méi)的繞射點(diǎn)、斷層點(diǎn)位置凸顯出來(lái),如圖中紅色箭頭所指之處。圖8a和圖8b中藍(lán)色方框部分的局部放大結(jié)果如圖8c和圖8d所示。從圖8d中能明顯區(qū)分出繞射點(diǎn),如紅色橢圓標(biāo)示,地下小型繞射體也被凸顯出來(lái)。因此,經(jīng)過(guò)混合域高分辨率Radon變換濾波后得到的繞射波成像剖面可以聯(lián)合反射波成像剖面進(jìn)行準(zhǔn)確、高確定性的地震解釋。

圖7 基于混合域高分辨率Radon變換的傾角域CIGs波場(chǎng)分離結(jié)果a 實(shí)際地震數(shù)據(jù)的全波場(chǎng)CIGs; b 波場(chǎng)分離后的繞射波CIGs

圖8 中國(guó)西部M探區(qū)某測(cè)線采用本文方法進(jìn)行波場(chǎng)分離前、后的結(jié)果a 波場(chǎng)分離前的PSDM剖面; b 波場(chǎng)分離后的繞射波成像剖面; c 圖8a的局部放大結(jié)果; d 圖8b的局部放大結(jié)果

3 結(jié)論

本文通過(guò)對(duì)傾角域CIGs中的波場(chǎng)特征進(jìn)行分析得到波場(chǎng)分離的依據(jù),進(jìn)一步利用模擬數(shù)據(jù)對(duì)比分析了常規(guī)最小二乘Radon變換、頻率域高分辨率Radon變換和混合域高分辨率Radon變換方法在波場(chǎng)分離中的效果,并結(jié)合實(shí)際地震資料的應(yīng)用,得到以下結(jié)論。

1) 在傾角域CIGs中,反射波地震響應(yīng)曲線總是呈現(xiàn)出開(kāi)口向上的“笑臉”狀曲線,繞射波同相軸則表現(xiàn)為擬線性。兩者曲線形態(tài)的顯著差異可以實(shí)現(xiàn)波場(chǎng)分離。

2) 與最小二乘Radon變換、頻率域高分辨率Radon變換的傾角域繞射波波場(chǎng)分離方法相比,混合域高分辨率Radon變換方法不僅可以有效保留子波的幅度,而且在提高分辨率方面均優(yōu)于其它兩種方法,在生產(chǎn)中具有更好的應(yīng)用價(jià)值。

3) 本文繞射波分離成像技術(shù)完整保留了地下高分辨率、超高分辨率的繞射波信息,分離的繞射波能夠揭示地下小規(guī)模斷裂構(gòu)造、縫洞儲(chǔ)層,并與反映地下連續(xù)反射層或主要不連續(xù)點(diǎn)(大斷層等)的反射能量聯(lián)合解釋,可有效提高地震解釋的精度。