薛懷艷, 張星宇

(中海石油(中國)有限公司 深圳分公司，廣州 510240)

疊前深度偏移初始速度模型建模方法應用研究
——以JY工區(qū)崎嶇海底為例

薛懷艷, 張星宇

(中海石油(中國)有限公司 深圳分公司，廣州 510240)

Dix公式法和速度約束反演法，是地震資料疊前深度偏移建立初始速度模型的兩種主要方法。首先利用JY工區(qū)深水崎嶇海底的資料對這兩種速度建模方法的成像能力進行了實驗分析，結果表明，這兩種方法都不能得到合理的初始速度模型。此外，通過速度約束反演得到的初始速度模型可以沿層位生成平均層速度模型。再利用垂向上地層埋深與速度關系和橫向上同一套沉積地層與速度關系來共同修正平均層速度模型，這樣可以獲得更符合地質條件的初始速度模型。實驗結果證明，通過引入地質背景的約束有利于建立起更為準確的初始速度模型，為后續(xù)的疊前深度偏移成像提供了有利保障。

崎嶇海底； 速度約束反演； 初始速度模型； 疊前深度偏移

0 引言

珠江口盆地東部，經過多年的勘探開發(fā)，勘探目標已經逐步的由淺水區(qū)向深水區(qū)延伸。由于從潛水到深水的過渡區(qū)域位于陸架破折帶上，海水深度由500 m急劇變?yōu)? 500 m左右，且溝壑縱橫，水道切割關系復雜，造成了海底速度場橫向變化劇烈，利用時間偏移已經不能夠對海底崎嶇區(qū)域的中深層進行準確成像[1]。深度偏移由于考慮了射線偏折和薄透鏡項的影響，依然是現今能夠適應崎嶇海底，對復雜構造進行準確成像最具潛力的技術[2]。對于深度偏移成像方法而言，從Kirchhoff疊前深度偏移、波動方程疊前深度偏移到高斯束疊前深度偏移和逆時偏移，一直吸引著國內、外眾多學者對其進行了深入系統(tǒng)地研究，為復雜構造成像提供了較為完善的技術和方法[3-5]。

關于速度建模地研究也一直是學術界比較活躍的研究領域，由于全波形速度反演還未能實現大規(guī)模工業(yè)化，所以基于層析反演理論的疊前深度偏移速度建模方法，仍將是工業(yè)界速度建模的主流技術[6]。依照目前的常規(guī)技術流程，速度建模主要可以分為2步：①利用疊前時間偏移的層速度模型，提取其低頻分量用來構建深度偏移的初始速度模型；②在初始速度模型的基礎上，利用層析成像與疊前深度偏移進行迭代更新初始速度模型直到速度模型達到較高精度，生成的成像道集拉平為止[7-8]。然而第①步建立的初始速度模型的精度直接決定著后續(xù)層析迭代反演的收斂程度，如果初始速度模型產生了較大的偏差，則通過層析反演修正后的速度模型就很可能會收斂到局部極值，對成像結果造成嚴重的影響[9-10]。

目前，建立深度域初始速度模型的方法主要有3種：①基于Dix公式的方法即直接利用Dix公式，將時間偏移的均方根速度轉換成深度域初始層速度模型[11]；②利用速度約束反演的方法(CVI)，其原理與Dix公式類似只是在反演速度模型的過程中增加了約束條件；③相干反演，該方法的原理是利用射線追蹤，計算每個層位在CMP道集中的時距曲線，時距曲線相關最大值所對應的層速度即為該層的速度，從原理上講相干反演的精度可能較高，但在實際生產中其操作復雜，對資料信噪比要求較高，且在速度橫向變化劇烈的地方多解性強，因此該方法本身實用性不高[12]。目前，在實際生產中用到的主流方法依然是Dix公式和速度約束反演法。這里依據JY地區(qū)實際地震資料，并結合對該地區(qū)已有關于速度研究的認識，總結出了一套更加實用的建立初始速度模型的方法，而且通過與傳統(tǒng)的方法進行分析對比，得出了對實際工作具有參考價值的結論和認識。

1 深度域初始速度模型原理和應用效果分析

1.1 Dix公式計算層速度模型

設有n層水平層狀介質各層層速度為vi，第n層的旅行時為tn，第一層到第層的均方根速度為vRn，則Dix公式為式(1)。

(1)

Dix公式由于簡單高效，能夠快速的利用均方根速度得到層速度，所以在工業(yè)生產中經常被使用。但算法本身只適用于水平層狀介質，在速度橫向變化大的區(qū)域利用Dix公式直接去反演層速度很容易造成高頻震蕩的假象，而且均方根速度較小的誤差就會對層速度的計算結果產生很大地影響[13-14]。由圖1(a)可以看出,均方根速度場中深層的速度變化已經受到了崎嶇海底地影響，從圖1(b)中可以看出，在海底溝壑較深的區(qū)域構造產狀普遍下傾，可能是由于海水的速度較低導致的假象，且中深層成像效果欠佳。從圖2可以看出，由Dix公式計算得到的層速度存在嚴重的震蕩，而這種震蕩與實際地層的速度變化趨勢不符。CRP成像道集的剩余曲率比較大，拉伸嚴重，因此可以認為這種簡單的基于垂向單道速度反演的方法，在速度橫向變化劇烈的區(qū)域不能正確地反演出層速度場的變化趨勢。

圖1 均方根速度場及其成像結果Fig.1 RMS velocity model and imaging result(a)均方根速度場；(b)Kirchhoff疊前時間偏移結果

圖2 Dix公式轉換得到的層速度模型及其CRP道集Fig.2 Interval velocity model gained through Dix method and CRP gathers(a)Dix公式轉換得到的層速度模型; (b)Dix公式轉換得到層速度進行Kirchhoff疊前深度偏移生成CRP道集

1.2 速度約束反演(CVI)計算層速度模型

速度約束反演實際上就是在Dix公式的基礎上進行了改進[15]，其核心思想是所反演得到的層速度不僅盡量地滿足Dix公式，而且在層內是連續(xù)變化的，滿足這樣的趨勢函數：

(2)

(3)

圖3是速度約束反演得到的相同測線的層速度模型及其Kirchhoff疊前深度偏移生成的CRP道集。圖2與圖3相比，兩者低頻分量具有一定相似性，但CVI反演速度模型的結果橫向連續(xù)性比Dix公式好，采用這種方法可以有效地消除部分速度異常點，減小由Dix公式反演產生的高頻振蕩對速度模型的影響。從Kirchhoff疊前深度偏移生成的同一成像位置的CRP道集，可以看到道集的成像質量相對于Dix公式在中深層3 000 m～5 000 m范圍內有所改善，但道集整體的能量聚焦性依然欠佳。通過這個實驗可以得到2個基本認識: ①速度約束反演相對于Dix公式具有更好的穩(wěn)定性; ②其反演速度場精度有限依然不能適應速度橫向變化大的情況。

圖3 CVI公式反演出的層速度及其CRP道集Fig.3 Interval velocity model gained through constrained velocity inversion and CRP gathers(a)CVI反演得到的層速度模型; (b)利用CVI公式反演出的層速度模型進行Kirchhoff疊前深度偏移生成CRP道集

2 深度域初始速度模型

上面提到的Dix公式法和速度約束反演法，都不能克服速度橫向變化大的問題。原因主要在于利用時間偏移迭代分析速度場的過程中，并沒有考慮地震射線路徑上速度場的變化。某一個具體的成像點的成像質量，只與該成像位置的均方根速度有關，而與上覆地層和相鄰位置的速度無關，這一點雖然使得疊前時間偏移速度分析和迭代更新變得簡單方便，但同時也嚴重降低了速度場分析的精度，降低了成像質量。甚至在某些速度橫向變化較大的區(qū)域，利用這種速度分析方法會產生較大的誤差，其生成的層速度模型尤其是低頻分量也會出現較為嚴重的偏差。利用這種層速度模型直接作為初始速度模型，有可能導致后續(xù)的疊前深度偏移速度迭代不收斂。因此，在建立初始速度模型的過程中應當引入地質背景的約束，利用埋深與速度在垂向上的變化關系和沉積地層與速度在橫向上的變化關系，合理地約束速度場的宏觀變化趨勢，從而獲得理想的初始速度模型。

圖4 沿層提取的平均層速度模型Fig.4 Average velocity of layer model generated along the horizon

圖4是在CVI反演獲得的層速度模型基礎上沿層提取的同一測線的平均層速度模型。根據對該區(qū)域以往地認識，該區(qū)域地層速度整體垂向上主要受埋深壓實控制，速度在橫向上主要受沉積地層的影響，速度沿同一套沉積地層變化相對穩(wěn)定。因此可以首先利用已有的地質解釋層位作為框架，提取出每一層的平均速度作為初始背景速度；然后統(tǒng)計不同層位的速度控制點與埋深的梯度關系，剔除異常值，在選取速度控制點的過程中，要以平面上沿層平均速度值作為參考值對所選的種子點進行質量控制；最后利用統(tǒng)計出的速度控制點與埋深的變化關系分別沿層位垂向和橫向進行插值，就重新構建起了新的初始層速度模型。

圖5是通過利用埋深和沉積地層關系共同約束建立起的初始速度模型及其Kirchhoff疊前深度偏移生成的CRP道集。從圖5(a)可以看出，新生成的速度模型不但保留了原始平均層速度的部分低頻分量，而且在垂向上層間層厚度內速度趨勢基本上符合埋深的梯度變化，在橫向上速度沿層位的變化也比較平緩，符合地質沉積規(guī)律。從圖5(b)可以看出，無論是成像質量還是能量聚焦性都明顯優(yōu)于前兩種方法，證明采用的這種在地學信息約束下建立起的初始速度模型是可靠的。

圖5 引入地學信息約束建立的初始層速度模型及其CRP道集Fig.5 Initial velocity model gained through introducing the constraints of geological information and CRP gathers(a)初始速度模型;(b)Kirchhoff疊前深度偏移生成CRP道集

3 疊前深度偏移成像效果

圖6(a)是利用高精度網格層析對圖5(a)地質信息約束下建立的初始速度模型迭代更新后的結果，以及在同一測線位置的最終疊前深度偏移結果和CRP成像道集。通過對比更新前后的速度模型可以發(fā)現，網格層析更新后的速度模型與初始速度模型依然具有較高的相似性，表明初始速度模型的低頻分量是比較穩(wěn)定的，同時也可以看出，在低頻分量穩(wěn)定性較好的前提下網格層析，可以進一步反演出速度模型的高頻分量，豐富速度模型細節(jié)信息。這一點不但在速度模型上有體現，而且還可以從CRP道集上得到映證,圖6(c)與圖5(b)相比，在深度3 800 m的位置CRP道集成像質量明顯改善，道集能量聚焦性增強，出現了較強的反射界面。從成像剖面上來看相比于疊前時間偏移，在海底溝壑深的位置CDP在1 400～1 500附近地層構造受海水的影響被消除，構造形態(tài)不在下傾并且真實形態(tài)得到恢復。此外，對于中深層的成像質量明顯優(yōu)于疊前時間偏移。

4 結論和認識

1)對JY區(qū)域而言，由于崎嶇海底導致海底速度場橫向變化較大，疊前時間偏移已經不能夠滿足對構造形態(tài)進行準確成像的要求，而疊前深度偏移可以取得比較理想的成像結果。

2) Dix公式和速度約束反演兩種方法在速度橫向變化劇烈的情況下,都不能夠反演出較為理想的初始速度模型，其Kirchhoff疊前深度偏移生成的CRP道集質量欠佳。

圖6 網格層析反演得到的速度模型及其Kirchhoff疊前深度偏移結果和CRP道集Fig.6 Interval velocity model gained through grid-based tomography inversion， Kirchhoff prestack depth migration result， CRP gathers(a)層速度模型;(b)Kirchhoff疊前深度偏移結果;(c)CRP道集

3)以速度約束反演獲得的平均層速度模型作為參考，再結合埋深與速度在垂向上的關系和同一套沉積地層與速度在橫向上的關系共同約束建立起的初始速度模型，能夠較準確地反映出速度模型的低頻趨勢，其疊前深度偏移生成的CRP道集質量較好。

4)相對于速度模型的高頻分量，速度模型的低頻分量對成像質量具有更大地影響。

5)在初始速度模型的低頻分量比較可靠的前提下，網格層析可以反演出速度模型的高頻分量，進一步提高成像質量。

[1] ZIQIANG XU,ZHONGYU FANG,ZHIWEI DAN,et al.The anisotropic prestack depth migration in deepwater rough bottom area[C].Beijing 2014 International Geophysical Conference & Exposition,2014:419-422.

[2] XIUMEI CHEN,ZAITIAN MA,HUAZHONG WANG,et al.Curved-wave prestack depth migration[J].GEOPHYSICS,2011,76(3)：115-129.

[3] CHAOSHUN HU, MRINAL K. SEN, PAUL STOFFA, et al. Plane-wave Gaussian-beam prestack depth migration[C]. SEG Technical Program Expanded Abstracts, 2007:2200-2204.

[4] YONGHE SUN, FUHAO QIN, STEVE CHECKLES, et al. Leveille. 3-D prestack Kirchhoff beam migration for depth imaging[J]. GEOPHYSICS, 2000,65(5)：1592-1603.

[5] YAN HUANG, BING BAI, HAIYONG QUAN, et al. Application of RTM 3D angle gathers to wide-azimuth data subsalt imaging[J]. GEOPHYSICS, 2011,76(5)：175-182.

[6] YI TAO, MRINAL K. SEN. Frequency-domain full waveform inversion with plane-wave data[J]. Geophysics,2013,78(1)：13-23.

[7] SCHNEIDER. J. F. Residual moveout and residual radii of curvature from approximate prestack depth migrations[C] . SEG Technical Program Expanded Abstracts, 2009:3605-3609.

[8] RICKETT, J., P. SAVA. Offset and angle-domain common image point gathers for shot-profile migration[J]. Geophysics, 2002,67(3): 883-889.

[9] HUA‐WEI ZHOU. Multiscale traveltime tomography [J]. Geophysics,2003,68(5):1639-1649.

[10]馬彥彥,李國發(fā),張星宇,等.疊前深度偏移速度建模方法分析[J]. 石油地球物理勘探,2014,49(4):687-693. MA Y Y, LI G F, ZHANG X Y, et al. Analysis on the strategy of building velocity model in prestack depth migration[J]. Oil Geophysical Prospecting,2014,49(4):687-693.(In Chinese)

[11]HEWITT DIX. C. Seismic velocities from surface measurements[J]. Geophysics, 1955, 20(1)：68-86.

[12]LANDA,E. THORE,P.SORIN,V.et al.Interpretation of velocity estimates from coherency inversion[J].Geophysics,1991,56(9)：1377-1383.

[13]HSING-ZEN LEE, KARL SCHLEICHER, HANS MEINARDUS. Refined dix approach for seismic velocity[C]. SEG Technical Program Expanded Abstracts, 1992:1206-1209.

[14]ZVI KOREN, IGOR RAVVE.Constrained Dix inversion[J]. Geophysics, 2006,71(6): 113-130.

[15]ZVI KOREN,IGOR RAVVE.Constrained velocity inversion[J].SEG Technical Program Expanded Abstracts,2005:2289-2292.

Application research on the strategy of building initial velocity model in prestack depth migration:A case from rough deep sea bottom in JY area

XUE Huaiyan, ZHANG Xingyu

(China National Offshore Oil Corporation(Shenzhen), Guangzhou 510240，China)

Dix formula method and constrained velocity inversion are two main initial model-building methods of prestack depth migration in seismic data processing. Using seismic data of rough deep sea bottom from JY area, the imaging abilities of the two methods were experimentally analyzed, and thus concluded that both methods could not obtain the reasonable initial velocity model. Furthermore, the average velocity of layer model can be generated along the horizon through velocity model obtained by using constrained velocity inversion. Then by using the relationship between vertical velocity function with depth, as well as the relationship between the same set of sedimentary strata with velocity on horizontal, which could common correct the average velocity of layer model. So we can get more in line with the initial velocity model of geological conditions. Experimental results show it could be helpful to establish a more accurate initial velocity model through introducing the constraints of geological background, which would provide the advantageous safeguard for subsequent imaging quality of prestack depth migration.

rough deep sea bottom; constrained velocity inversion; initial velocity model; prestack depth migration

2015-11-04 改回日期：2016-10-19

薛懷艷(1965-)，女，高級工程師，研究方向為油氣勘探開發(fā)高分辨率地震資料處理與解釋，E-mail：xuehy@cnooc.com.cn。

1001-1749(2017)01-0103-06

P 631.4

A

10.3969/j.issn.1001-1749.2017.01.15