胡光輝,杜澤源,何兵紅,孫思宇

(中國石油化工股份有限公司石油物探技術(shù)研究院,江蘇南京211103)

精細油氣勘探開發(fā)對高精度速度建模成像技術(shù)提出了更高的要求。而近地表速度建模精度直接影響地震資料靜校正及成像效果。特別是在我國西部山地、盆嶺及黃土塬等地區(qū),地表條件復(fù)雜,高差變化劇烈,巖性變化大,導(dǎo)致中淺層速度橫向變化大,近地表低降速帶發(fā)育,這些因素對近地表速度建模帶來了嚴峻的挑戰(zhàn)。因此,迫切需要有效進行近地表速度建模方法的研究以精確刻畫近地表地質(zhì)結(jié)構(gòu),進而消除近地表帶來的影響。

經(jīng)過多年研究,形成了很多針對性的方法,一定程度上解決了復(fù)雜近地表建模的很多問題,比如微測井、折射波法、面波法、層析反演方法等[1]。最常用的近地表建模方法主要利用地震波中的直達波、折射波、回轉(zhuǎn)波等最早被檢波器接收到的波,通過基于射線理論或波形理論的層析或波形反演方法進行建模。早期,由于更偏向于大尺度地質(zhì)體的成像以及計算硬件條件的限制,主要使用基于射線理論的旅行時層析方法[2]。隨后專家學者針對走時層析射線追蹤[3-4]、程函方程數(shù)值求解[5]及目標泛函選取[6-7]等方法都提出了許多改進之處,并成功解決了近地表建模的問題[8-9]。但是,如果近地表存在低降速區(qū),將嚴重影響地震記錄,產(chǎn)生波形畸變,使走時層析反演方法得到的結(jié)果不夠理想;如果近地表存在高速異常體,會屏蔽地震波的透射,大部分地震波能量被反射回去,從而在射線追蹤時產(chǎn)生全反射現(xiàn)象,使異常體內(nèi)部缺失射線,嚴重影響此區(qū)域的反演精度。針對復(fù)雜近地表速度異常及復(fù)雜介質(zhì)情況,波形反演方法可以彌補傳統(tǒng)射線理論建模的“盲區(qū)”。

近年來,基于波形類的反演方法發(fā)展迅速。全波形反演(full waveform inversion,FWI)基于全波動方程,綜合考慮運動學和動力學特征,地震波模擬更加符合實際地震波場,而且受低降速帶或高速異常體的影響較小,是目前公認的建模精度最高的方法之一。LAILLY[10]和TARANTOLA[11]提出了基于波場誤差反傳構(gòu)建梯度算子的思想,建立了全波形反演的基礎(chǔ)理論框架。全波形反演理論已經(jīng)出現(xiàn)近40年,但是由于計算條件的限制,直到近年來才成為研究與生產(chǎn)應(yīng)用的熱點[12],并且陸續(xù)出現(xiàn)了實際資料成功應(yīng)用的實例,比如針對海上資料[13-14]、陸上資料[15]及鹽丘模型的反演應(yīng)用[16]。全波形反演在實際應(yīng)用過程中,發(fā)現(xiàn)了許多挑戰(zhàn)性問題,尤其是在陸上三維實際地震資料應(yīng)用中,仍面臨許多實際問題,如全波形反演因正、反演理論及地震資料與模型參數(shù)之間的關(guān)系等因素會產(chǎn)生較強的非線性性;全波形反演在尋優(yōu)過程中,由于全波場數(shù)據(jù)與模型之間的耦合,很難找到全局最優(yōu)解。對于復(fù)雜近地表全波形反演速度建模也存在許多問題。利用全波形反演解決近地表問題時,如果仍然使用全波場信息,很容易使全波形反演陷入局部極值問題,難以得到最優(yōu)的速度模型。因此,必須提取對近地表速度建模最有用的地震數(shù)據(jù)信息進行反演。SHENG等[17]提出了近地表折射波數(shù)據(jù)的早至波層析方法;HANAFY等[18]研究了陸上利用淺層地震數(shù)據(jù)早至波波形反演方法確定淺層地下水情況。胡光輝等[19]研究了采用基于早至波的特征波波形反演方法解決近地表速度建模問題。早至波是指初至以及初至到達之后一段時間內(nèi)被檢波器記錄下來的波的集合[17,19],包括直達波、首波、折射波、回轉(zhuǎn)波等。這些波在近地表傳播,基本不受反射界面影響,攜帶了豐富的近地表速度異常體的信息,因此,我們試圖利用這些早至波進行全波形反演,實現(xiàn)對近地表速度模型的建構(gòu)。ZHOU等[20]研究了基于早至波和反射波波形反演的氣云區(qū)速度建模。早至波全波形反演基于波場誤差反傳構(gòu)建梯度算子的全波形反演理論框架,利用早至波的走時和波形信息,進行數(shù)據(jù)擬合匹配,利用迭代優(yōu)化求解算法,逐步優(yōu)化速度模型。筆者通過模型測試驗證了早至波全波形反演對射線類層析建模“盲區(qū)”的改善,以及早至波全波形反演速度建模的精度[19]。本文分析了全波形反演中非線性來源,研究了早至波提取及早至波正演模擬方法,然后介紹了早至波全波形反演在實際資料應(yīng)用中的關(guān)鍵技術(shù)流程,最后給出了該技術(shù)在實際三維地震資料中的應(yīng)用結(jié)果。

1 早至波全波形反演

早至波與初至波一樣,并不是特指某一類型的波,而是包含了很多波的信息。近偏移距早至波主要包括淺層的直達波、潛波及小折射波等。這些波沒有經(jīng)過反射界面改造,攜帶了豐富的近地表速度異常體的信息,因此,可以利用這些早至波,通過全波形反演方法,構(gòu)建含速度異常體的近地表模型[17]。大偏移距早至波主要包括中深層的折射波、反射波以及潛波等,這些波的信息對中深層速度模型的構(gòu)建有著重要意義。因此,對于近地表速度建模來說,可以僅利用小偏移距的早至波信息來實現(xiàn)中淺層的速度建模,避免使用大偏移距數(shù)據(jù)可以減少累積誤差,降低問題的非線性性。為了獲得理想的近地表速度模型,提高對近地表速度異常體的認識,將早至波引入到基于最小二乘理念的全波形反演框架中,利用早至波信息,完成基于早至波的全波形反演,實現(xiàn)近地表復(fù)雜區(qū)速度建模。

1.1 敏感核函數(shù)分析

在全波形反演的目標泛函里,反射波是主要的非線性來源。在沒有準確初始速度模型的情況下,無法將反射點準確歸位,因此在進行早至波全波形反演時,將反射信息切除,利用早至波信息是降低反演非線性性的有效手段,尤其是在近地表建模方面。圖1所示的敏感核函數(shù)充分說明了這一點。

從圖1中的敏感核函數(shù)可以看出,其有效區(qū)域為炮點與檢波點連線的菲涅爾帶。反射響應(yīng)與其射線路徑并不相似,如果背景速度不準確,則更新過程中不能校正速度模型,這也是反射波強非線性的根本所在。全波形反演中,由于采用地面觀測方式,模型的淺層分辨率要比深層更高。因此,早至波部分往往具有更高的分辨率,擬合具有較高分辨率的模型淺層部分數(shù)據(jù)要容易得多。利用早至波進行近地表建模,甚至中深層建模(大偏移距早至波信息)是行之有效的全波形反演手段。

1.2 早至波正演模擬方法

常規(guī)地震數(shù)據(jù)處理中,早至波往往被當作噪聲予以切除。而我們在需要利用早至波信息時,就要首先對地震數(shù)據(jù)進行預(yù)處理,然后提取保幅、保真的早至波數(shù)據(jù)。近地表速度建模的早至波全波形反演主要利用近偏移距的早至波信息,因此,在數(shù)據(jù)準備階段,需要對早至波進行準確識別與提取。本文采用切除手段,根據(jù)地質(zhì)背景優(yōu)化選取切除函數(shù),移除觀測數(shù)據(jù)中反射波同相軸,保留數(shù)據(jù)中的直達波、回折波和折射波信息。圖2給出了某三維工區(qū)用于近地表建模的地震資料中早至波識別與提取實例。

圖1 敏感核函數(shù)

圖2 某三維工區(qū)早至波識別與提取實例

準確模擬出與實際資料相吻合的早至波是早至波全波形反演的基礎(chǔ)。反演過程中,采用基于阻尼的波動方程正演方法實現(xiàn)早至波模擬。根據(jù)觀測記錄的切除函數(shù),設(shè)定阻尼區(qū)域,模擬對應(yīng)的早至波信息。設(shè)置阻尼因子可以直接獲得早至波,不用模擬出全波場數(shù)據(jù)后再進行切除,從而避免了冗余過程及誤差累積,同時保證在反演中只利用早至波信息,消除反射波信息,降低反演非線性程度。

(1)

式中:v為模型速度;ρ為模型密度;x,y和z分別為空間坐標方向。在正演模擬時對阻尼因子加一時移來衰減早至波之后波形,以保證早至波信息的完整。

對圖2中所示三維工區(qū)資料進行早至波正演模擬,得到的模擬地震記錄如圖3所示?？梢钥闯?正演模擬得到的早至波與實際資料切除得到的早至波基本吻合。這也說明了本文基于阻尼的早至波波動方程正演模擬方法的有效性。

圖3 正演模擬的早至波場

1.3 早至波全波形反演方法

全波形反演是一個非線性問題,其目標泛函一般定義為正演波場與觀測波場誤差的L2范數(shù)。本文早至波全波形反演基于傳統(tǒng)的L2范數(shù)全波形反演框架,定義目標泛函為:

(2)

式中:‖·‖2表示L2范數(shù);Δd為觀測波場與模擬波場之間的殘差。在提取觀測數(shù)據(jù)早至波時,需要根據(jù)不同的地質(zhì)任務(wù)和目的層,只保留適當?shù)男∑凭嘣缰敛ㄐ畔ⅰ?/p>

公式(2)在形式上與常規(guī)全波形反演類似,不同的是,公式中誤差僅為觀測早至波與模擬早至波之間的誤差。目標泛函的求解采用梯度類算法,通過梯度求解得到迭代更新模型,為方便理解,寫成頻率域的形式:

(3)

早至波全波形反演同樣是一個非線性問題,需要多次迭代直到誤差收斂。本文采用L-BFGS算法求取迭代步長,可在時間域引入海森矩陣進行照明補償,提高反演的穩(wěn)定性。

本文早至波全波形反演的技術(shù)流程如圖4所示。主要過程是:首先利用初至走時層析方法恢復(fù)近地表模型低波數(shù)信息,并將其作為早至波全波形反演的初始模型,然后利用早至波信息進行反演,恢復(fù)模型的高波數(shù)成分,從而建立高精度近地表速度模型。相對于常規(guī)全波形反演,早至波全波形反演選取近偏移距淺層的地震波信息,降低了反演的非線性性,反演過程更加穩(wěn)定,收斂效果更好[19]。本文初至波走時層析與早至波全波形反演逐級速度建模流程引入了多尺度反演策略,不依賴于初始模型的精度,可大幅提高近地表速度建模的精度和效率。

圖4 早至波全波形反演流程

2 實際資料應(yīng)用

2.1 實際數(shù)據(jù)1

將本文方法應(yīng)用于中國東部某陸上工區(qū)三維實際地震資料,工區(qū)近地表部分存在明顯的低降速帶,近地表速度建模不準對后續(xù)資料處理影響較大。工區(qū)地形較平緩,高差相對較小,受地震地質(zhì)條件的影響,資料信噪比不同,干擾波的類型及特點也不同,面波、折射波、直達波干擾較強,高頻干擾普遍存在于整個工區(qū)。早至波全波形反演對資料的品質(zhì)十分敏感,因此在進行早至波全波形反演之前需要對原始地震資料進行預(yù)處理,主要包括:①去面波、去壞道、去壞炮、去異常振幅及線性噪聲等,保證在進行預(yù)處理時既能最大限度壓制隨機干擾,又能不損壞有效信號;②球面擴散補償、地表一致性振幅補償,削弱大地的衰減效應(yīng)和球面擴散效應(yīng)以及由于不同的激發(fā)條件和地面檢波器因地表條件引起的不耦合因素造成的振幅不一致影響。本測試工區(qū)共選取445炮地震數(shù)據(jù),震源距為200m,檢波點距為50m,原始最大偏移距為5500m,經(jīng)處理后的偏移距范圍是500～3000m。工區(qū)大小為Lx=8375m,Ly=4775m,早至波全波形反演模型深度Lz=1250m。

首先提取觀測地震數(shù)據(jù)早至波信息,圖5為其中一炮數(shù)據(jù)提取的早至波,保留了比較完整的早至波信息。試驗中所用地震子波為通過原始地震數(shù)據(jù)估計的子波,圖6展示了估計的地震子波,子波的能量主要集中在4～15Hz。

圖5 提取的早至波

圖6 原始地震記錄估計的地震子波

然后采用商業(yè)軟件構(gòu)建初始層析速度模型,如圖7 所示。再以此速度模型作為初始模型進行早至波全波形反演,結(jié)果如圖8所示。由圖8可以看出,相較于商業(yè)軟件層析速度模型,早至波全波形反演結(jié)果細節(jié)更加豐富,成層性更好,速度異常體空間展布更加清晰。

全波形反演的質(zhì)量監(jiān)控可以通過對比正演模擬記錄和觀測記錄來實現(xiàn)。為驗證本文早至波全波形反演效果,從數(shù)據(jù)域的角度進行互相關(guān)結(jié)果對比。圖9 為第40炮觀測地震記錄的歸一化非零延遲自相關(guān)結(jié)果。從圖9可以看出,當炮記錄延遲為0時相關(guān)的能量達到最大。圖10為第40炮觀測地震記錄和初始模型正演模擬地震記錄的互相關(guān)結(jié)果,顯示模擬結(jié)果和觀測記錄在零延遲處能量沒有聚焦,說明初始模型和真實模型差距較大。圖11為第40炮觀測地震記錄和反演結(jié)果正演模擬的地震記錄的互相關(guān)結(jié)果,在零延遲處,互相關(guān)值達到最大,說明在一定程度上模擬記錄的相位已校正到正確位置。

圖7 商業(yè)軟件構(gòu)建的初始層析速度模型

圖8 早至波全波形反演得到的速度模型

2.2 實際數(shù)據(jù)2

將本文方法應(yīng)用于中國西部某工區(qū)三維實際地震資料。該工區(qū)淺層存在較明顯的速度異常體。在進行早至波全波形反演之前同樣對地震資料進行相應(yīng)的預(yù)處理。本測試工區(qū)共選取635炮地震數(shù)據(jù),震源距為200m,檢波點距為50m,原始炮集最大偏移距為6000m,經(jīng)處理后的偏移距范圍是500～3000m。工區(qū)大小為Lx=10250m,Ly=6250m,早至波全波形反演模型深度Lz=1750m。試驗中所用地震子波為通過原始地震數(shù)據(jù)估計的子波,子波的能量主要集中在4～15Hz。首先使用商業(yè)軟件走時層析方法對工區(qū)進行反演,建立大尺度背景速度場,然后使用本文所述早至波全波形反演方法建立高精度速度模型,速度建模結(jié)果如圖12所示。

圖9 第40炮觀測地震記錄的歸一化非零延遲自相關(guān)結(jié)果

圖10 第40炮觀測地震記錄和初始模型模擬地震記錄歸一化非零延遲互相關(guān)結(jié)果

圖11 第40炮觀測地震記錄和反演結(jié)果模擬的地震記錄歸一化非零延遲互相關(guān)結(jié)果

圖12a為商業(yè)軟件走時層析結(jié)果;圖12b為早至波全波形反演結(jié)果。由圖12可見,使用早至波全波形反演方法得到的速度場在走時層析結(jié)果的基礎(chǔ)上更加精細,可以將速度異常體清晰地描繪出來。

分別使用圖12中兩個速度模型進行疊后深度偏移,結(jié)果如圖13所示。其中,圖13a為使用商業(yè)軟件層析結(jié)果進行偏移得到的剖面;圖13b為使用早至波全波形反演結(jié)果進行偏移得到的剖面。對比圖13a 和圖13b可以看出,在橢圓標識處,使用早至波全波形反演建模結(jié)果得到的偏移剖面斷面更加清晰,同相軸橫向連續(xù)性更好,分辨率更高,而且中深層成像質(zhì)量有較大改善,更符合地質(zhì)認識成果。

圖12 不同方法反演的速度模型a 商業(yè)軟件走時層析; b 早至波全波形反演

圖13 疊后深度偏移結(jié)果a 使用商業(yè)軟件走時層析模型; b 使用早至波全波形反演模型

3 結(jié)論與認識

本文將早至波全波形反演方法應(yīng)用于實際工區(qū)三維地震資料近地表速度建模,建立了走時層析與早至波全波形反演相結(jié)合的多尺度速度建模技術(shù)流程。分析表明,反射波是全波形反演強非線性的主要來源,選取合適的近偏移距早至波資料進行反演,可有效降低反演的非線性,實現(xiàn)高精度近地表速度建模。早至波全波形反演方法的關(guān)鍵是地震資料早至波識別與選取、早至波正演模擬、早至波全波形反演算法。早至波近地表速度建模并不能完全使用所有的早至波信息,正確選擇早至波以及高精度早至波正演模擬方法是早至波全波形反演近地表速度建模的基礎(chǔ)。早至波全波形反演方法與常規(guī)速度建模方法相比,更加適應(yīng)近地表低速異常,縱、橫向速度變化以及低降速帶發(fā)育等情況。本文方法在三維地震資料中的實際應(yīng)用情況表明,早至波近地表速度建模技術(shù)得到的建模結(jié)果優(yōu)于常規(guī)射線層析類建模方法,建模精度明顯提高。相對于商業(yè)軟件走時層析結(jié)果,早至波全波形反演得到的速度細節(jié)更加豐富,速度異常體刻畫更加清晰。從疊后深度偏移剖面可以看出,使用早至波全波形反演速度的成像結(jié)果斷面更加清楚,深層成像質(zhì)量明顯改善。

對于早至波全波形反演,地震子波的提取,起伏地表情況下自由界面的準確加載、轉(zhuǎn)換波的準確識別及模擬等都是亟待解決的問題,影響早至波全波形反演的建模精度,且往往地表起伏較大等近地表復(fù)雜地區(qū)的實際地震資料信噪比較低,這些都是進一步研究的重點?？偟膩碚f,早至波是地震炮集里信噪比最高、最容易識別的有效信號。早至波全波形反演無需人工拾取,自動化程度較高,建模精度優(yōu)于傳統(tǒng)射線類方法,具有廣闊的應(yīng)用前景。