管文勝 段文勝 查 明 孫 琦 宋小平 蘇 洲

（①中國(guó)石油大學(xué)（華東）地球科學(xué)與技術(shù)學(xué)院，山東青島266580；②中國(guó)石油塔里木油田公司勘探開發(fā)研究院，新疆庫爾勒841000；③帕拉代姆技術(shù)（北京）有限公司，北京100083）

利用基于模型的層析速度反演進(jìn)行低幅度構(gòu)造成像

管文勝*①②段文勝②查 明①孫 琦②宋小平③蘇 洲②

（①中國(guó)石油大學(xué)（華東）地球科學(xué)與技術(shù)學(xué)院，山東青島266580；②中國(guó)石油塔里木油田公司勘探開發(fā)研究院，新疆庫爾勒841000；③帕拉代姆技術(shù)（北京）有限公司，北京100083）

基于層位的層析反演在低信噪比地區(qū)無法進(jìn)行層位解釋，在劃分速度界面時(shí)也無法保證模型精度，局限性較大；基于網(wǎng)格的層析反演不容易收斂到實(shí)際速度模型，且模型沒有地質(zhì)意義的約束，有時(shí)會(huì)產(chǎn)生不符合地質(zhì)規(guī)律的反演結(jié)果。為此，基于模型的層析反演方法，借鑒了層位約束和全局反演優(yōu)點(diǎn)，以地質(zhì)層位為基礎(chǔ)建立模型，在橫向上沿地質(zhì)層位網(wǎng)格進(jìn)行參數(shù)更新，在縱向上層位間網(wǎng)格更新尺度隨著模型的地質(zhì)層位而變化，既解決了薄層速度建模問題，又能得到精準(zhǔn)的高、低頻速度分量和各向異性參數(shù)。其次，基于模型的層析反演將模型中地質(zhì)層位的深度作為一個(gè)變量納入層析反演數(shù)據(jù)空間，可以同時(shí)反演速度結(jié)構(gòu)與反射點(diǎn)位置（或反射層深度），具有更高的反演精度。將基于模型的層析反演方法應(yīng)用于YM工區(qū)實(shí)際資料，并通過各向異性疊前深度偏移處理，提高了速度模型精度，解決了該區(qū)低幅度構(gòu)造成像問題。

低幅度構(gòu)造 基于模型的層析反演 層析成像矩陣 VTI介質(zhì) 疊前深度偏移

1 引言

基于射線追蹤的層析反演技術(shù)是目前疊前深度偏移最常用的速度模型建立方法，主要利用偏移和層析交替迭代的方法進(jìn)行速度反演，能夠恢復(fù)速度場(chǎng)中的高、低波數(shù)信息，被廣泛用于深度域速度模型的建立［1-5］。在工業(yè)界又細(xì)分為兩種具體方法：①基于層位的層析反演，主要通過在簡(jiǎn)單層狀模型的多個(gè)速度界面求取射線路徑和剩余延遲時(shí)進(jìn)行速度反演［6］。其優(yōu)點(diǎn)是能夠表達(dá)復(fù)雜的地質(zhì)模型，模型受地質(zhì)構(gòu)造的約束，具有地質(zhì)意義；缺點(diǎn)是需要輸入構(gòu)造解釋模型，在低信噪比地區(qū)無法進(jìn)行層位解釋，在劃分速度界面時(shí)也無法保證模型精度，而且這種簡(jiǎn)單模型簡(jiǎn)化了速度模型的細(xì)節(jié)，局限性較大；②基于網(wǎng)格的層析反演，通過在全局求取射線路徑和剩余延遲時(shí)進(jìn)行速度反演，無需輸入解釋模型，克服了層狀模型在層間速度細(xì)節(jié)刻畫上的缺陷，反演精度較高，可以獲得較為準(zhǔn)確的速度場(chǎng)高頻分量；但不容易收斂到實(shí)際速度模型［7-9］，且模型沒有地質(zhì)意義的約束，有時(shí)會(huì)產(chǎn)生不符合地質(zhì)規(guī)律的反演結(jié)果。

基于模型的層析反演方法，借鑒了層位約束和全局反演優(yōu)點(diǎn)，以地質(zhì)層位為基礎(chǔ)建立模型，在橫向上沿地質(zhì)層位網(wǎng)格進(jìn)行參數(shù)更新，在縱向上層位間網(wǎng)格更新尺度隨著模型的地質(zhì)層位而變化，既解決了薄層速度建模問題，又能得到精準(zhǔn)的高、低頻速度分量和各向異性參數(shù)。其次，基于模型的層析反演將模型中地質(zhì)層位的深度作為一個(gè)變量納入層析反演數(shù)據(jù)空間，可以同時(shí)反演速度結(jié)構(gòu)與反射點(diǎn)位置（或反射層深度），具有更高的反演精度。本文將基于模型的層析反演方法應(yīng)用于YM工區(qū)實(shí)際資料，并通過各向異性疊前深度偏移處理，提高了速度模型精度，解決了該區(qū)低幅度構(gòu)造的成像問題。

2 基于模型的層析基本原理

2.1 基于模型的層析速度反演原理

基于模型的層析速度反演包括一個(gè)正演過程（計(jì)算炮點(diǎn)、反射點(diǎn)、檢波點(diǎn)的旅行時(shí)）和一個(gè)反演過程（根據(jù)旅行時(shí)和剩余時(shí)差更新速度模型）。首先建立層析成像矩陣，即通過對(duì)模型中的每一個(gè)反射點(diǎn)做射線追蹤完成正演，一條射線路徑對(duì)應(yīng)一個(gè)等式，其合集即為層析成像矩陣

式中：tm為模型中的地震波旅行時(shí)；Sijk為第ijk網(wǎng)格中的慢度；d ijk為第ijk網(wǎng)格中的射線路徑長(zhǎng)度。

而實(shí)際模型中地震波會(huì)出現(xiàn)旅行時(shí)差，可表示為

式中：Δtm為旅行時(shí)誤差；ΔSijk為第ijk網(wǎng)格中的慢度誤差。

旅行時(shí)差Δtm可以表示為

式中：A v為d ijk組成的系數(shù)矩陣，又稱靈敏因子；Δv為未知變量ΔSijk的向量。

根據(jù)偏移CRP道集中反射波同相軸的上翹或下垂，通過人工或自動(dòng)拾取的方式得到剩余旅行時(shí)Δt。以共軛梯度法求解最小二乘解即可獲得Δv，從而得到更新后的速度體。

在式（2）中只考慮了慢度誤差對(duì)射線旅行時(shí)的影響，而忽略了速度更新對(duì)射線路徑的影響。假設(shè)速度更新量為正值，反射點(diǎn)位置在新速度模型中的深度位置下移，導(dǎo)致射線路徑變長(zhǎng)，求取的速度比真實(shí)值偏高，反之亦然。因此，為了消除速度更新對(duì)射線路徑長(zhǎng)度的影響，將模型的反射點(diǎn)位置的深度變化代入式（1）和式（2）之和，即

式中Δd ijk為模型中射線路徑長(zhǎng)度變化分?jǐn)偟降趇jk網(wǎng)格中的長(zhǎng)度。略去高階項(xiàng)ΔSijkΔd ijk，則式（4）變?yōu)?/p>

式（5）既考慮了速度誤差的影響，也考慮了反射點(diǎn)位置的深度變化，因此更加合理。式（5）又可表示為

在VTI介質(zhì)中，式（6）表示為

式中：ε、δ為VTI介質(zhì)各向異性參數(shù)；Aδ、Aε分別為δ、ε的系數(shù)矩陣；Δδ、Δε分別為δ、ε的誤差向量；A z為由Sijk組成的系數(shù)矩陣；Δz為Δd ijk的向量，即模型中反射點(diǎn)深度變化向量，通過求解該變量，即可反演模型中地質(zhì)層位的深度（圖1）。

2.2 各向異性速度分析

Thomsen［10］于1986年給出了表征VTI介質(zhì)彈性性質(zhì)的5個(gè)參數(shù)，并對(duì)其進(jìn)行定義和描述，其中描述P波各向異性的參數(shù)為VP、ε和δ，其中VP表示P波的垂直方向傳播速度，ε為P波各向異性的度量，δ值決定了P波在垂直方向附近相速度的各向異性程度。隨后，人們對(duì)VTI介質(zhì)的各向異性速度進(jìn)行了研究［11，12］。

對(duì)某一套地層，通過合成記錄標(biāo)定可以得到頂、底反射時(shí)差，用單程時(shí)差與已知的地震層速度相乘得到層厚度ΔZI（未考慮各向異性），而實(shí)際的測(cè)井記錄則準(zhǔn)確記錄了該套地層的頂、底深度，考慮了各向異性影響的真實(shí)層厚度表示為ΔZA，則

在分析地震資料的各向異性時(shí)，一般基于弱各向異性假設(shè)，即縱波速度取決于參數(shù)δ和ε，當(dāng)ε=δ時(shí)，簡(jiǎn)化為橢圓各向異性。因此在獲取了各向異性參數(shù)δ后，先假設(shè)ε=δ，進(jìn)行基于VTI介質(zhì)的初始疊前深度偏移，然后進(jìn)行ε迭代，以求取準(zhǔn)確的ε參數(shù)。

3 基于模型的層析成像速度模型建立方法

3.1 初始層速度模型的建立

3.1.1 淺層速度建模

由于覆蓋次數(shù)低、炮檢距較小的原因，在YM工區(qū)近地表難以獲得準(zhǔn)確、可靠的剩余延遲，也就難以求取層速度。淺層速度模型的誤差會(huì)累積影響到深層目的層的成像。為此，采用微測(cè)井信息約束大炮初至的層析反演，彌補(bǔ)了應(yīng)用反射波方法求取近地表速度的缺陷，得到的淺層速度場(chǎng)更加可靠。

圖1 網(wǎng)格層析反演（a）與模型層析反演（b）網(wǎng)格尺度對(duì)比

3.1.2 中深層速度建模

YM工區(qū)的中深層構(gòu)造相對(duì)簡(jiǎn)單，可采用模型約束反演方法得到初始層速度。通過對(duì)疊前時(shí)間偏移處理得到的較準(zhǔn)確的均方根速度場(chǎng)進(jìn)行約束，建立一個(gè)模型約束的瞬時(shí)速度場(chǎng)。具體步驟為：首先建立初始的低頻趨勢(shì)模型速度場(chǎng)，然后對(duì)每一個(gè)反演的垂向函數(shù)，假設(shè)局部變化是一維模型，再使用最小二乘基本原理求解反問題。

3.2 基于模型的層析速度更新與優(yōu)化

結(jié)合前面介紹的基本原理，基于模型的層析速度建模步驟如下。

（1）以近地表速度高速頂為分界面，將由大炮初至層析速度模型得到的淺層速度和由模型約束反演得到的中深層速度進(jìn)行拼接，并作為初始速度進(jìn)行速度迭代。

（2）計(jì)算旅行時(shí)并利用Kirchhoff積分法進(jìn)行疊前深度偏移，通過大量CMP道集振幅樣點(diǎn)的不斷歸位，使真正反射點(diǎn)位置的振幅不斷疊加，從而提高信噪比。

（3）根據(jù)區(qū)域地質(zhì)情況，在深度域地震數(shù)據(jù)體或時(shí)間域偏移數(shù)據(jù)體上選取測(cè)井聲波速度變化的界面進(jìn)行地震解釋層位拾取，這樣的層位通常對(duì)應(yīng)地震強(qiáng)反射同相軸，信噪比較高，獲得的剩余延遲可靠。在模型的地質(zhì)層位上按固定間隔選取反射點(diǎn)（圖1藍(lán)色圓點(diǎn)）進(jìn)行射線追蹤，可以大量減少反射點(diǎn)數(shù)量，節(jié)約計(jì)算量。

（4）通過交互方式和試驗(yàn)確定合適的參數(shù)，如深度誤差值（深度偏移資料深度與實(shí)鉆井深度差）等，在研究區(qū)CRP偏移道集上進(jìn)行批量應(yīng)用，以獲得剩余延遲。

（5）在層位模型、偏移數(shù)據(jù)體、剩余延遲體、層速度體、δ體、ε體以及地震反射層屬性數(shù)據(jù)體（傾角體、方位角體、連續(xù)性體等）的基礎(chǔ)上建立數(shù)據(jù)庫。

（6）對(duì)每一個(gè)反射點(diǎn)進(jìn)行射線追蹤完成正演，其中每一條射線路徑對(duì)應(yīng)一個(gè)等式，從而建立層析反演的層析成像矩陣（式（5））。

（7）求解層析成像矩陣，更新每個(gè)網(wǎng)格的速度（式（7）、式（8）），得到初步迭代的各向異性速度體、δ體和ε體。

（8）經(jīng)過層析速度反演的多次迭代，完成層析正演和反演過程，直至近道集基本拉平、剩余延遲為零，并輸出最終的各向異性速度體、δ體、ε體以及各向異性疊前深度偏移數(shù)據(jù)體。

結(jié)合YM工區(qū)實(shí)際情況，確定了基于模型層析速度反演的VTI各向異性疊前深度偏移流程（圖2）。

圖2 基于模型層析速度反演的VTI各向異性疊前深度偏移流程

4 應(yīng)用實(shí)例及效果分析

TLM盆地YM工區(qū)三維地震資料于2006年采集，覆蓋次數(shù)為48，目的層為白堊系巴西改組。疊前時(shí)間偏移資料表明，該區(qū)存在一個(gè)北東向的低速帶，發(fā)育一個(gè)大型低幅度斷背斜構(gòu)造（圖3a），在該構(gòu)造有利部位部署的第一口井（1號(hào)井）即獲得突破。為加快該構(gòu)造的評(píng)價(jià)，在該構(gòu)造的3個(gè)高部位同時(shí)部署了3口評(píng)價(jià)井（2、3、4號(hào)井），評(píng)價(jià)井的設(shè)計(jì)高點(diǎn)海拔比1號(hào)井高20～30m，完鉆后3口井的實(shí)際海拔比1號(hào)井低20～30m，因而相繼失利（圖3b）。經(jīng)綜合研究認(rèn)為，由于工區(qū)速度復(fù)雜，沒有考慮各向異性，因而造成假構(gòu)造的低速帶，導(dǎo)致鉆探失利。因此，針對(duì)速度變化大、存在低速陷阱等各向異性突出的地區(qū)，需采用基于模型層析速度反演的VTI介質(zhì)疊前深度偏移技術(shù)，以提高低幅度構(gòu)造成像質(zhì)量和精度，并發(fā)現(xiàn)和落實(shí)與斷層相關(guān)的低幅度構(gòu)造。

圖3 YM工區(qū)各向同性速度平面圖（a）和構(gòu)造圖（b）圖a中紅色區(qū)域?yàn)榈退賲^(qū)，圖b中紅色區(qū)域?yàn)楸承睒?gòu)造區(qū)

從CRP道集上看，通過模型層析反演反復(fù)迭代后，CRP遠(yuǎn)道集被明顯拉平，目的層段能量聚焦，信噪比得到提高（圖4b）。

基于模型的層析反演在求取速度模型時(shí)，考慮到每個(gè)深度成像點(diǎn)，同時(shí)速度的更新以地質(zhì)層位為對(duì)象，在不同層位依據(jù)剩余延遲情況獲得不同的速度變化值，速度精度得以提高，地質(zhì)界面速度發(fā)生了明顯變化（圖5）。

式（8）表明：當(dāng)各向同性速度大于各向異性速度時(shí)，δ＞0；當(dāng)各向同性速度小于各向異性速度時(shí)，δ＜0。分析YM工區(qū)更新迭代后的δ體及ε體數(shù)據(jù)表明（圖6），目的層段的δ、ε值小于零（黃色至紅色范圍），表明各向同性速度小于各向異性速度，可能會(huì)造成異常低速區(qū)，使深度域構(gòu)造圖的深度比實(shí)際深度小，從而出現(xiàn)假構(gòu)造（圖3）。

圖7為各向同性、基于模型速度反演的各向異性深度域成像剖面。由圖可見：各向同性深度域剖面（圖7a）與基于模型速度反演的各向異性深度域剖面（圖7b）的整體構(gòu)造形態(tài)基本一致，但后者（圖7b）的信噪比更高，波組特征清楚，斷面成像清晰可靠，目的層段的儲(chǔ)層內(nèi)部反射特征明顯，低幅度構(gòu)造成像質(zhì)量好于前者（圖7a），有利于識(shí)別巖性圈閉；從地震地質(zhì)標(biāo)定來看，各向異性深度域剖面（圖7b）的古近系底（Bot_E）和白堊系巴西改組頂（Bot_K1bs）、底（Bot_K1bx）的層位深度比各向同性深度域剖面（圖7a）的大。

圖4 各向同性（a）和各向異性（b）道集

圖5 初始速度（a）、模型層析迭代更新速度（b）

圖6 各向異性參數(shù)δ體（a）和ε體（b）

圖7 各向同性（a）、基于模型速度反演的各向異性（b）深度域成像剖面

圖8 各向同性、各向異性疊前深度偏移與實(shí)鉆的深度數(shù)據(jù)

圖8為各向同性、各向異性疊前深度偏移與實(shí)鉆的深度數(shù)據(jù)。由圖可見，基于模型層析速度反演的VTI各向異性疊前深度偏移的深度數(shù)據(jù)與鉆井深度數(shù)據(jù)基本一致，而各向同性疊前深度偏移的深度比鉆井深度小。因此，基于模型層析的各向異性偏移可以滿足YM工區(qū)的低幅度構(gòu)造的地震成像精度需求。圖9為YM工區(qū)疊前深度偏移各向異性速度場(chǎng)和構(gòu)造圖。由圖可見，1號(hào)井處于構(gòu)造高點(diǎn)上，2、3、4號(hào)井處于構(gòu)造斜坡部位，證實(shí)北部的低速帶不存在。

5 結(jié)束語

基于模型的層析反演方法具有以下優(yōu)點(diǎn)：

（1）模型受地質(zhì)層位約束，即在橫向上是沿地質(zhì)層位的網(wǎng)格進(jìn)行速度和各向異性參數(shù)更新，尤其在縱向上層間網(wǎng)格尺度隨著地質(zhì)層位而變化，解決了薄層速度建模問題，避免了單一網(wǎng)格尺度層析成像可能產(chǎn)生的射線覆蓋程度過低的問題。同時(shí)這種多尺度層析方法具有粗網(wǎng)格和細(xì)網(wǎng)格的剖分，既兼顧層析成像的分辨率和精度，又能得到精準(zhǔn)的高、低頻速度分量和各向異性參數(shù)；

（2）基于模型的層析反演將模型的地質(zhì)層位深度值作為一個(gè)變量納入層析反演，即對(duì)速度結(jié)構(gòu)與反射點(diǎn)位置（或反射層深度）同時(shí)反演，具有更高的反演精度，接近地下實(shí)際速度模型，與鉆井深度更加吻合，使偏移成像質(zhì)量更高；

（3）在迭代更新過程中，模型考慮了地層中的各向異性特點(diǎn)，如δ、ε等反映速度各向異性的參數(shù)，符合地下介質(zhì)的實(shí)際情況。

總之，基于模型的層析反演方法解決了低速陷阱、速度橫向變化大的各向異性問題，得到的各向異性疊前深度偏移體的地質(zhì)現(xiàn)象清晰、可靠，斷層歸位準(zhǔn)確，其深度與鉆井深度差別小，落實(shí)了與斷層相關(guān)的低幅度構(gòu)造。該方法有具有使用方便、計(jì)算效率高的特點(diǎn)，適合工業(yè)化推廣應(yīng)用。

圖9 YM工區(qū)疊前深度偏移各向異性速度場(chǎng)（a）和構(gòu)造圖（b）

［1］ 李振春.地震偏移成像技術(shù)研究現(xiàn)狀與發(fā)展趨勢(shì).石油地球物理勘探，2014，49（1）：1-21.Li Zhenchun.Research status and development trends for seismic migration technology.OGP，2014，49（1）：1-21.

［2］ Stork C.Reflection tomography in the postmigrated domain.Geophysics，1992，57（5）：680-692.

［3］ Wang B，Pann K，Meek R A.Macro velocity model estimation through model based globally-optimized residual-curvature analysis.SEG Technical Program Expanded Abstracts，1995，14：1084-1087.

［4］ Bloor R.Building velocity models using tomography.SEG Technical Program Expanded Abstracts，1998，17：1195-1198.

［5］ Woodward M，F(xiàn)armer P，Nichols D et al.Automated 3D tomographic velocity analysis of residual moveout in prestack depth migrated common image point gathers.SEG Technical Program Expanded Abstracts，1998，17：1218-1221.

［6］ Kosloff D，Sherwood J，Koren Z et al.Velocity and interface depth determination by tomography of depth migrated gathers.Geophysics，1996，61（5）：1511-1523.

［7］ 馬彥彥，李國(guó)發(fā)，張星宇等.疊前深度偏移速度建模方法分析.石油地球物理勘探，2014，49（4）：687-693.Ma Yanyan，Li Guofa，Zhang Xingyu et al.Strategy of velocity model building in prestack depth migration.OGP，2014，49（4）：687-693.

［8］ 李慧，成德安，金婧.網(wǎng)格層析成像速度建模方法與應(yīng)用.石油地球物理勘探，2013，48（增刊1）：12-16.Li Hui，Cheng Dean，Jin Jing.Velocity model building based on grid tomography.OGP，2013，48（S1）：12-16.

［9］ Robein E.Seismic Imaging：A Review of the Techniques，Their Principles，Merits and Limitations.Hutten，Netherlands，2010.

［10］ Thomesen L.Weak elastic anisotropy.Geophysics，1986，51（10）：1954-1966.

［11］ Tsvankin I D.橫向各向同性介質(zhì)速度分析和時(shí)深轉(zhuǎn)換（張青偉譯）.國(guó)外油氣勘探，1997，9（4）：504-513，517.

［12］ 李源，劉偉，劉微等.各向異性全速度建模技術(shù)在山地地震成像中的應(yīng)用.石油物探，2015，54（2）：157-164.Li Yuan，Liu Wei，Liu Wei et al.Application of anisotropic full velocity modeling in the mountainous seismic imaging.GPP，2015，54（2）：157-164.

P631

A

10.13810／j.cnki.issn.1000-7210.2017.01.013

管文勝，段文勝，查明，孫琦，宋小平，蘇洲.利用基于模型的層析速度反演進(jìn)行低幅度構(gòu)造成像.石油地球物理勘探，2017，52（1）：87-93.

1000-7210（2017）01-0087-07

*山東省青島市黃島區(qū)長(zhǎng)江西路66號(hào)中國(guó)石油大學(xué)（華東）地球科學(xué)與技術(shù)學(xué)院，266580。Email：guanws-tlm＠petrochina.com.cn

本文于2015年12月30日收到，最終修改稿于2016年11月21日收到。

本項(xiàng)研究受國(guó)家重大科技專項(xiàng)課題（2011ZX05001-002-003、2011ZX05001-001-006）聯(lián)合資助。

（本文編輯：劉勇）

管文勝 高級(jí)工程師，1968年生；1991年獲江漢石油學(xué)院（現(xiàn)長(zhǎng)江大學(xué)）勘查地球物理專業(yè)學(xué)士學(xué)位；2009年獲中國(guó)地質(zhì)大學(xué)（武漢）構(gòu)造地質(zhì)學(xué)專業(yè)碩士學(xué)位；現(xiàn)為中國(guó)石油大學(xué)（華東）地球科學(xué)與技術(shù)學(xué)院博士研究生。目前在塔里木油田分公司勘探開發(fā)研究院從事地震處理、解釋一體化及構(gòu)造、圈閉評(píng)價(jià)等方面的研究。