戴海濤，成劍冰，王紅博，劉懷山，劉志兵

（1.中國(guó)海洋大學(xué)海洋地球科學(xué)學(xué)院，山東青島266100；2.中國(guó)石油東方地球物理公司研究院，河北涿州072750）

1 真地表疊前深度偏移技術(shù)存在的問(wèn)題

1.1 近地表圓滑偏移面分析

中國(guó)西部地區(qū)地表起伏劇烈，低降速帶巨厚，高速地層和低速覆蓋層分布無(wú)規(guī)律，地震波射線的入射角和出射角變化范圍大，靜校正會(huì)造成地震波射線在走時(shí)計(jì)算時(shí)出現(xiàn)誤差，進(jìn)而造成波場(chǎng)失真[1-2]。

常規(guī)時(shí)間域疊前道集的靜校正處理是將靜校正和剩余靜校正的總量在CMP道集域進(jìn)行高頻、低頻分離，分離出的低頻分量是以CMP為單元統(tǒng)計(jì)出的時(shí)間域地形平均面，即CMP參考面。早期的疊前深度偏移處理通常通過(guò)替換速度將時(shí)間域的CMP參考面換算到深度域，作為深度偏移基準(zhǔn)面。復(fù)雜山地地表高程變化大，造成基于真實(shí)地表與CMP參考面計(jì)算得到的深度偏移基準(zhǔn)面差異大，因而射線路徑的走時(shí)誤差大，波場(chǎng)畸變嚴(yán)重。為了解決波場(chǎng)畸變問(wèn)題，目前深度偏移成像時(shí)通常采用對(duì)地表高程進(jìn)行小尺度平滑，選取地表圓滑面作為深度域偏移基準(zhǔn)面的處理方法。地表圓滑面反映地表高程的變化形態(tài)，應(yīng)盡可能減少其與地表的高差。該處理思路是為了消除地表圓滑面與CMP參考面之間的低頻時(shí)差，將激發(fā)點(diǎn)和接收點(diǎn)的真實(shí)位置垂直校正到地表圓滑面上，對(duì)偏移前的道集數(shù)據(jù)應(yīng)用高頻靜校正，以減小真實(shí)激發(fā)點(diǎn)、接收點(diǎn)與地表圓滑面之間的射線旅行時(shí)誤差。

以上方法技術(shù)雖然削弱了道集時(shí)差對(duì)偏移成像的影響，但仍然存在以下問(wèn)題：①雖然地震數(shù)據(jù)中引入了高頻靜校時(shí)差，但道集時(shí)差對(duì)地震波場(chǎng)走時(shí)的影響并未完全消除；②深度偏移基準(zhǔn)面是根據(jù)時(shí)間域基準(zhǔn)面靜校正量求得的，它與后續(xù)建立的深度域近地表速度模型不一定匹配。

1.2 淺層速度建模問(wèn)題分析

疊前深度偏移的核心是獲得準(zhǔn)確的速度模型。中、深層速度模型的準(zhǔn)確程度不僅取決于道集數(shù)據(jù)信噪比，還依賴于淺表層速度模型的精度。淺層速度誤差會(huì)傳遞至下伏地層，進(jìn)而影響深層速度模型的精度，因此建立準(zhǔn)確的淺層速度模型對(duì)深度偏移成像至關(guān)重要。利用反射波進(jìn)行速度層析反演，要求反射波具有高信噪比、大排列長(zhǎng)度、可準(zhǔn)確拾取共反射點(diǎn)（common reflection point，CRP）道 集 的 剩 余 延 遲（residual moveout，RMO）。但因三維地震淺層資料存在覆蓋道數(shù)少、信噪比低的缺點(diǎn)，使得共成像點(diǎn)道集淺層的反射波同相軸很難被識(shí)別，剩余延遲拾取不準(zhǔn)確。因此用于評(píng)價(jià)偏移速度正確與否的判別準(zhǔn)則（同相軸校平）失去其依據(jù)，利用反射波層析成像速度優(yōu)化方法無(wú)法求取準(zhǔn)確的近地表及淺層地層速度。圖1 是我國(guó)西部某復(fù)雜山地的三維地震速度譜及CRP道集，可以看出中深層的反射信號(hào)具有較高信噪比以及較大排列長(zhǎng)度，淺層地震資料信噪比低、覆蓋道數(shù)少，難以識(shí)別有效反射信息，因此利用反射波層析建立準(zhǔn)確的淺層速度模型存在困難。

圖1 西部地區(qū)某復(fù)雜山地三維地震速度譜（a）及CRP道集（b）

2 真地表淺層速度建模技術(shù)

2.1 思路及策略

利用初至波層析反演出高精度的淺表層速度模型，并與偏移基準(zhǔn)面配套應(yīng)用，可建立合理的深度偏移淺層速度模型，實(shí)現(xiàn)真地表疊前深度偏移。本文研究并提出了面向復(fù)雜山地低信噪比的真地表淺層速度建模技術(shù)，核心技術(shù)思路包括如下內(nèi)容：

1）在微測(cè)井[3]、鉆測(cè)井等速度信息約束的基礎(chǔ)上，利用回折波層析反演[4-7]方法建立近地表及淺層速度模型，以提高淺層速度模型的精度。

2）靜校正直接面向深度域成像，即將時(shí)間域大炮初至數(shù)據(jù)的回折波反演得到的近地表模型直接作為深度域淺層速度模型，將時(shí)間域靜校正僅作為疊前去噪等處理方法的輔助手段，因此需要在疊前深度偏移前對(duì)應(yīng)用了基準(zhǔn)面靜校正量的地震數(shù)據(jù)進(jìn)行靜校正“反應(yīng)用”，時(shí)間域靜校正引起的時(shí)差問(wèn)題可采用真地表偏移面＋淺層建模＋深度偏移的處理方法解決，根據(jù)深度偏移成像質(zhì)量決定靜校正方法及優(yōu)選參數(shù)。

3）在表層速度準(zhǔn)確的情況下，選取包含高精度旅行時(shí)算法的層析成像和偏移方法，完成整體速度模型的建立以實(shí)現(xiàn)真地表疊前深度偏移。

2.2 微測(cè)井約束回折波分層層析反演

復(fù)雜山地采集的淺層地震資料信噪比低、覆蓋次數(shù)少、缺乏遠(yuǎn)炮檢距波場(chǎng)信息，反射波層析成像無(wú)法有效求取準(zhǔn)確的近地表速度，因此大炮初至數(shù)據(jù)的回折波層析反演技術(shù)[8-11]成為時(shí)間域靜校正處理和深度偏移近地表速度建模的有效技術(shù)手段。

回折波層析反演考慮了地震波在近地表介質(zhì)中傳播的各種現(xiàn)象，對(duì)回折波或連續(xù)折射直達(dá)波通過(guò)高精度的交互反演得到近地表介質(zhì)速度變化情況，從而建立復(fù)雜的近地表模型。但受采集資料的限制，當(dāng)最小炮檢距、低降速層厚度較小時(shí)，利用初至波信息無(wú)法準(zhǔn)確刻畫(huà)表層極低速度的信息，需要聯(lián)合微測(cè)井和近炮檢距初至波信息反演低速層速度，再以反演速度為約束，聯(lián)合微測(cè)井和中、遠(yuǎn)炮檢距初至波信息反演高精度的降速層速度，從而建立精確的近地表速度模型。

大炮初至數(shù)據(jù)的回折波層析反演技術(shù)對(duì)近地表結(jié)構(gòu)進(jìn)行高密度速度單元?jiǎng)澐郑瑢⒒卣鄄ǖ淖邥r(shí)描述為對(duì)介質(zhì)慢度函數(shù)沿射線路徑的線積分，回折波射線追蹤利用初始速度模型和炮檢點(diǎn)位置計(jì)算旅行時(shí)，初至?xí)r間與正演旅行時(shí)之差可被反演成射線經(jīng)過(guò)的每個(gè)速度單元慢度的擾動(dòng)。

為了使反演結(jié)果與實(shí)際模型更吻合并使反演結(jié)果穩(wěn)定合理，本文采用了微測(cè)井約束回折波層析反演方法。該方法的目標(biāo)函數(shù)包含了速度約束項(xiàng)，先將微測(cè)井解釋的近地表模型作為約束條件[11]，再與旅行時(shí)殘差共同構(gòu)成目標(biāo)函數(shù)，最后求解拉格朗日約束條件下的目標(biāo)函數(shù)最小二乘解。在層析反演中，由M條射線和N個(gè)未知數(shù)建立的層析矩陣方程為[12-15]：

式中：A為射線路徑集合的Jacobi矩陣；ΔS為慢度修正向量；ΔT為旅行時(shí)殘差矩陣?？煞謩e表示為：

假設(shè)通過(guò)微測(cè)井獲得了L個(gè)準(zhǔn)確的速度信息，即方程（1）中的N個(gè)速度單元中有L個(gè)已知量，用L維方程可表示為：

式中：C為已知的L個(gè)約束條件的Jacobi矩陣；F為約束旅行時(shí)殘差矩陣。

其中，

在初至波層析中，假如[～S1，～S2，…，～SL]T為采用微測(cè)井或其它方法獲得的部分已知近似速度，那么（4）式C和F的矩陣元素分別為：

1）利用微測(cè)井約束回折波分層層析反演求取高精度近地表速度模型。

2）中深層的速度模型的建立：①利用時(shí)間域均方根速度場(chǎng)，通過(guò)約束速度反演得到時(shí)間域?qū)铀俣龋虎谠跇?gòu)造模型約束條件下，通過(guò)時(shí)深轉(zhuǎn)換獲取深度域?qū)铀俣?；③利用VSP等鉆、測(cè)井資料對(duì)速度場(chǎng)進(jìn)行修正得到中、深層速度模型。

3）在近地表速度模型中選取一個(gè)平滑的高速層頂界面作為速度融合界面，在設(shè)定的時(shí)窗內(nèi)將近地表速度鑲嵌到中、深層速度模型中，從而獲得用于深度偏移的初始速度模型，在后續(xù)速度模型迭代修正過(guò)程中，保持近地表速度模型不變，僅對(duì)高速層頂界面以下的地層進(jìn)行速度迭代更新。

解決起伏地表的波場(chǎng)畸變問(wèn)題，不僅需要高精度的近地表速度模型，還必須將偏移基準(zhǔn)面選為真實(shí)地表面，從而實(shí)現(xiàn)真地表疊前深度偏移。為此需要正確處理時(shí)間域靜校正與偏移基準(zhǔn)面的關(guān)系。既然深度域淺層速度模型與時(shí)間域靜校正是一致的，即已經(jīng)在深度域利用近地表速度解決了時(shí)間域靜校正問(wèn)題，那么應(yīng)用了靜校正的時(shí)間域偏移前道集數(shù)據(jù)應(yīng)該在深度域偏移前對(duì)應(yīng)用的靜校正展開(kāi)“反應(yīng)用”，使深度偏移前的數(shù)據(jù)無(wú)需進(jìn)行基準(zhǔn)面靜校正。

圖3 采用不同迭代次數(shù)得到的淺層速度模型及對(duì)應(yīng)的疊前深度偏移剖面

對(duì)柴達(dá)木盆地西部某工區(qū)回折波進(jìn)行層析反演迭代，得到的淺層速度模型及疊前深度偏移剖面如圖3所示，本工區(qū)共進(jìn)行了4輪迭代，可以看出隨著迭代次數(shù)的增加，速度場(chǎng)細(xì)節(jié)逐漸清晰，最終得到了精確的淺層速度模型，疊前深度偏移剖面成像質(zhì)量也得到了大幅改善。對(duì)該工區(qū)應(yīng)用微測(cè)井等資料約束前、后的回折波層析反演結(jié)果如圖4 所示，圖4c和圖4d中藍(lán)線代表VSP測(cè)井速度，紅線代表層析成像反演的層速度?？梢钥闯?，經(jīng)微測(cè)井資料約束后，回折波層析反演的表層速度橫向刻畫(huà)更精細(xì)，整體速度與VSP測(cè)井速度在淺、中層一致性更好，淺、中層速度很好地反映出逆沖斷層速度上沖的趨勢(shì)，符合目前的地質(zhì)認(rèn)識(shí)，疊前深度偏移剖面上反射波同相軸連續(xù)性明顯增強(qiáng)，斷面波特征更清楚，斷層下部成像質(zhì)量也得到明顯改善。

2.4 真地表疊前深度偏移

目前制約復(fù)雜山地真地表疊前深度偏移實(shí)際應(yīng)用的因素主要包括高精度淺表層速度建模和高精度偏移算法。提高旅行時(shí)計(jì)算精度，使旅行時(shí)計(jì)算結(jié)果能夠描述真實(shí)近地表地層速度橫向變化特點(diǎn)及規(guī)律，即讓深度偏移解決時(shí)間域的靜校正問(wèn)題，有利于深度域高精度淺層速度模型建立。在建立淺層速度模型的基礎(chǔ)上，實(shí)現(xiàn)真地表疊前深度偏移還需要選取合適的層析成像和偏移方法，但無(wú)論如何，旅行時(shí)計(jì)算都是重要步驟之一，在條件允許的情況下，須盡可能采用波前構(gòu)建法和頻率域波動(dòng)方程數(shù)值解法進(jìn)行旅行時(shí)計(jì)算。建立高精度的整體速度模型后，可采用波場(chǎng)外推的偏移方法實(shí)現(xiàn)真地表疊前深度偏移。

圖4 應(yīng)用微測(cè)井等資料約束前、后的回折波層析反演結(jié)果

采用多方位各向異性層析成像或者更為先進(jìn)的繞射層析成像方法可提高速度建模精度。通常采用由淺及深、逐步遞進(jìn)的建模方法，在獲得高精度的淺層速度模型后，采用基于井控地質(zhì)導(dǎo)向約束的TTI網(wǎng)格層析成像迭代技術(shù)完成整體速度建模。該技術(shù)是在數(shù)字化露頭、地層傾角測(cè)井、三維地震構(gòu)造解釋模型等多信息的約束下，通過(guò)聯(lián)合反演減少層析解的不確定性，并確保速度模型在符合地質(zhì)及構(gòu)造認(rèn)識(shí)規(guī)律的條件下迭代更新[18-19]。實(shí)際資料處理結(jié)果表明，網(wǎng)格尺度影響走時(shí)計(jì)算精度，進(jìn)而對(duì)速度模型的細(xì)節(jié)刻畫(huà)造成影響，因此在層析反演過(guò)程中，可通過(guò)逐步細(xì)分層析反演的網(wǎng)格，提高反演結(jié)果的分辨率。

3 應(yīng)用效果分析

圖5和圖6是西部某工區(qū)近地表圓滑面疊前深度偏移剖面和真地表疊前深度偏移剖面，該區(qū)開(kāi)展過(guò)多輪次的處理攻關(guān)，以往速度建模都是在近地表圓滑面上開(kāi)展的，存在斷層刻畫(huà)不清楚，構(gòu)造細(xì)節(jié)未落實(shí)，部分地區(qū)井震深度誤差大的問(wèn)題。此次采用了真地表疊前深度偏移技術(shù)，首先通過(guò)微測(cè)井約束初至回轉(zhuǎn)波分層層析反演建立深度0～1500 m 的速度模型，然后采用基于井控地質(zhì)導(dǎo)向約束的多方位網(wǎng)格層析[20]成像及TTI網(wǎng)格層析成像建立并優(yōu)化深度大于1500 m 的速度模型，最后利用TTI真地表逆時(shí)疊前深度偏移實(shí)現(xiàn)高精度偏移成像。由圖5可見(jiàn)，真地表疊前深度偏移剖面上逆沖斷層下盤(pán)鹽內(nèi)斷背斜構(gòu)造形態(tài)清晰可靠，成像效果較近地表圓滑面的疊前深度偏移成果明顯改善，主要表現(xiàn)在以下方面：

1）真地表疊前深度偏移剖面，反射波特征突出，波組特征清晰，偏移歸位準(zhǔn)確，斷點(diǎn)、斷面位置更清楚，有利于精細(xì)構(gòu)造解釋。

2）真地表疊前深度偏移較好地恢復(fù)了地下界面的真實(shí)形態(tài)。

如圖6所示，與以往疊前深度偏移的結(jié)果相比，真地表疊前深度偏移得到的剖面構(gòu)造形態(tài)變化明顯，利用本工區(qū)已經(jīng)完鉆的兩口井資料驗(yàn)證可知，真地表疊前深度偏移的結(jié)果與井吻合程度高，證明真地表疊前深度偏移的結(jié)果真實(shí)可靠。

圖5 西部某工區(qū)近地表圓滑面疊前深度偏移剖面（a）及真地表疊前深度偏移剖面（b）

圖6 西部某工區(qū)以往疊前深度偏移剖面（a）及真地表疊前深度偏移剖面（b）

4 結(jié)論與認(rèn)識(shí)

總結(jié)真地表淺層速度建模的研究及其實(shí)際應(yīng)用成果得出以下結(jié)論：

1）真地表疊前深度偏移淺層建模技術(shù)在很大程度上解決了近地表圓滑偏移面與時(shí)間域靜校正引起的波場(chǎng)畸變問(wèn)題。

2）微測(cè)井約束回折波分層層析反演是一種可靠地建立高精度淺層速度模型的方法，在多種信息的約束下，只有在獲得高精度淺層速度模型的基礎(chǔ)上，才能夠建立準(zhǔn)確的中、深層模型。

3）復(fù)雜山地的真地表疊前深度偏移成像精度受多種因素的影響，只有從旅行時(shí)算法、層析成像、偏移算法等多方面深入開(kāi)展理論研究與實(shí)際資料試驗(yàn)，提高其中每一個(gè)環(huán)節(jié)的計(jì)算精度，才能夠解決山前帶的復(fù)雜構(gòu)造成像問(wèn)題。