朱相羽,潘成磊,湯國松,劉立民,熊 燁

(1.中國石油化工股份有限公司江蘇油田分公司勘探管理部,江蘇揚(yáng)州225000;2.中國石油化工股份有限公司江蘇油田分公司物探研究院,江蘇南京210046;3.中國石油化工股份有限公司石油物探技術(shù)研究院,江蘇南京211103)

深度域速度模型的精度決定了深度偏移成像的精度。層析建模技術(shù)一般基于偏移距域共成像點(diǎn)道集,利用旅行時(shí)信息進(jìn)行速度更新。層析建模技術(shù)具有計(jì)算效率高的優(yōu)點(diǎn),目前仍然是工業(yè)界深度域速度建模的主要手段。最早利用DIX公式將時(shí)間域均方根速度轉(zhuǎn)換成深度域?qū)铀俣?但是DIX公式存在水平層狀介質(zhì)的假設(shè)導(dǎo)致結(jié)果往往存在較大誤差。自20世紀(jì)70年代以來,地震層析成像技術(shù)廣泛應(yīng)用于深度域速度建模。ENGDAHL[1]提出采用打靶法進(jìn)行射線追蹤以計(jì)算射線路徑和旅行時(shí)間,但是在復(fù)雜構(gòu)造區(qū)打靶法往往存在陰影區(qū);JULIAN等[2]提出彎曲射線法,該方法在走時(shí)層析應(yīng)用比較廣泛;LANDA等[3]建立了基于偏移前共中心點(diǎn)數(shù)據(jù)的射線追蹤方法,通過偏移前道集的相干性大小來計(jì)算深度域?qū)铀俣?AL-YAHYA[4]基于偏移后共炮檢距地震數(shù)據(jù)校平的原則建立了基于剩余時(shí)差的層析建模方法;STOLK等[5]和STORK等[6]將層析技術(shù)應(yīng)用于疊前深度偏移;隨著寬方位地震受到重視,BUIA等[7]提出了一套基于多方位的層析建模技術(shù)流程;段文勝等[8-9]利用寬方位網(wǎng)格層析建模技術(shù)進(jìn)行西部地震資料串珠成像,取得較好效果,但他們研究的是西部全方位地震資料的層析策略,沒有涉及到東部寬方位地震資料的層析策略。蘇北盆地寬方位地震資料橫縱比一般為0.5～0.7,方位角寬度沒有到達(dá)全方位的程度。地震資料存在主測(cè)線方向覆蓋次數(shù)明顯高于聯(lián)絡(luò)線方向覆蓋次數(shù)的問題,會(huì)導(dǎo)致剩余曲率拾取的密度和質(zhì)量在不同方位存在較大差異,因而影響寬方位層析反演精度。為了解決這個(gè)問題,本文提出一套基于OVT域方位角分組的寬方位層析速度建模技術(shù)流程,目的是解決蘇北盆地地震資料聯(lián)絡(luò)線方向覆蓋次數(shù)不足影響到剩余曲率拾取的質(zhì)量和密度的問題,實(shí)現(xiàn)提高寬方位層析反演精度的目標(biāo)。

1 寬方位層析建模技術(shù)

目前業(yè)內(nèi)深度域速度建模大多采用網(wǎng)格層析建模技術(shù),采用計(jì)算效率較高的克?；舴虔B前深度偏移方法得到共偏移距域共成像點(diǎn)道集(ODCIG),以共成像點(diǎn)道集拉平程度作為判斷疊前深度偏移速度場(chǎng)準(zhǔn)確程度的標(biāo)準(zhǔn)。如果偏移速度偏高,共成像點(diǎn)道集上同相軸下彎,反之同相軸上翹。因此共成像點(diǎn)道集上不同偏移距同相軸深度的差異提供了速度更新的信息。采用網(wǎng)格層析建模技術(shù)實(shí)現(xiàn)速度模型更新[10],即:

Δze(x,y,h,mc)=ze(x,y,h,mc)-
zr(x,y,h=0,mc)

(1)

式中:mc代表當(dāng)前速度場(chǎng);ze代表共成像點(diǎn)道集上偏移距為h的深度剖面;zr代表零偏移距的深度剖面;Δze代表道集的未校平程度。假設(shè)成像點(diǎn)地震道偏移距為0,成像深度接近真實(shí)深度,即:

ze(x,y,h=0,mc)≈zd

(2)

式中:zd代表地層真實(shí)深度。

方程(1)可近似表達(dá)為:

(3)

隨著寬方位地震數(shù)據(jù)越來越受到重視,1988年,VERMEER[11]提出了炮檢距向量片(OVT)的概念。OVT的優(yōu)勢(shì)在于保留數(shù)據(jù)的方位角信息,有利于利用方位各向異性信息進(jìn)行裂縫預(yù)測(cè),OVT域技術(shù)在裂縫預(yù)測(cè)等方面取得較大的進(jìn)展[12-16]。隨著OVT域研究的深入,提出了針對(duì)寬(全)方位地震數(shù)據(jù)的速度建模技術(shù)[17]。寬方位層析建模過程主要包括以下3個(gè)方面。①將地震數(shù)據(jù)由共中心點(diǎn)域分選到OVT域,在一個(gè)十字排列內(nèi)按照檢波線距和炮線距將地震數(shù)據(jù)劃分成多個(gè)小矩形,每個(gè)小矩形代表一個(gè)OVT,對(duì)每一個(gè)OVT進(jìn)行編號(hào)。將整個(gè)三維工區(qū)所有的十字排列具有相同OVT編號(hào)的地震數(shù)據(jù)集合在一起構(gòu)成一個(gè)OVT片,理論上每一個(gè)OVT片是覆蓋整個(gè)工區(qū)的單次覆蓋數(shù)據(jù),OVT片個(gè)數(shù)等于工區(qū)的覆蓋次數(shù)。同一個(gè)OVT片的地震數(shù)據(jù)具有相近的偏移距和方位角。②使用同一個(gè)偏移速度對(duì)每一個(gè)OVT片分別進(jìn)行疊前偏移得到OVT域共成像點(diǎn)道集(OVG),OVG攜帶共成像點(diǎn)的三維坐標(biāo)、偏移距和方位角信息。③剩余曲率拾取和寬方位層析反演,一般在不同方位地震道分布比較均勻的全方位地震資料可以直接在OVG上進(jìn)行剩余曲率(深度時(shí)差)拾取和寬方位層析反演。但是,蘇北盆地寬方位地震資料橫縱比一般為0.5～0.7,主測(cè)線方向覆蓋次數(shù)是聯(lián)絡(luò)線方向覆蓋次數(shù)的2倍左右,存在主測(cè)線方向射線密度明顯高于聯(lián)絡(luò)線方向射線密度的問題,這會(huì)導(dǎo)致剩余曲率拾取的密度和質(zhì)量在不同方位存在較大差異,影響寬方位層析結(jié)果。為了解決這個(gè)問題,本文根據(jù)蘇北盆地寬方位地震資料方位角分布特點(diǎn)對(duì)寬方位層析流程進(jìn)行優(yōu)化,對(duì)OVG進(jìn)行方位角分組,目的是保證每個(gè)方位角組剩余曲率拾取的密度和質(zhì)量基本相當(dāng)。具體步驟是將OVG按照方位角劃分成幾組。假設(shè)每組數(shù)據(jù)對(duì)應(yīng)相同的方位角,將每組的中心方位角作為該組數(shù)據(jù)共同的方位角。對(duì)每組數(shù)據(jù)分別進(jìn)行剩余曲率拾取,盡可能保證每組數(shù)據(jù)剩余曲率拾取的準(zhǔn)確性,而不需要考慮不同方位角組之間剩余曲率拾取結(jié)果是否相同。分別對(duì)每個(gè)方位角組進(jìn)行射線追蹤,然后建立線性方程組。將線性方程組及其對(duì)應(yīng)的方位角作為輸入,進(jìn)行寬方位網(wǎng)格層析反演,得到一個(gè)速度更新量,將此速度更新量加到原速度模型得到更新后的速度模型。圖1展示了常規(guī)層析建模技術(shù)(圖1a)與寬方位層析建模技術(shù)(圖1b)思路上的差別。

圖1 常規(guī)層析建模技術(shù)(a)與寬方位層析建模技術(shù)(b)思路

1) 用于偏移的數(shù)據(jù)域不同。常規(guī)層析建模技術(shù)多基于共偏移距成像道集,進(jìn)行共偏移距域克希霍夫深度偏移,產(chǎn)生的共成像點(diǎn)道集按照偏移距由小到大進(jìn)行排序。共成像點(diǎn)道集保留偏移距信息,但不保留方位角信息。寬方位層析建模技術(shù)基于OVG,OVG保留偏移距和方位角信息。

2) 剩余曲率的拾取模式不同。網(wǎng)格層析建模技術(shù)根據(jù)同相軸的校平程度計(jì)算速度更新量,因此道集上剩余曲率拾取模式影響層析反演結(jié)果。常規(guī)層析建模技術(shù)認(rèn)為共成像點(diǎn)道集上同相軸在偏移距方向上表現(xiàn)為下彎或者上翹反映了速度模型偏高或者偏低。剩余曲率曲線形式上與偏移距有關(guān),與方位角無關(guān)(圖2a)。而寬方位層析建模技術(shù)認(rèn)為速度場(chǎng)不準(zhǔn)確時(shí),OVG同相軸剩余曲率曲線形式上不僅同偏移距有關(guān),同時(shí)還與方位角有關(guān)(圖2b)。

圖2 常規(guī)層析剩余曲率拾取(a)和寬方位層析剩余曲率拾取(b)示意

2 應(yīng)用實(shí)例分析

2019年,江蘇油田實(shí)施了YA高密度三維地震數(shù)據(jù)采集。YA高密度三維橫縱比為0.72。從方位角寬度來講YA高密度三維屬于寬方位三維。

在對(duì)YA高密度三維以復(fù)雜小斷塊成像為目標(biāo)進(jìn)行處理的過程中,通過井震對(duì)比確定Thomsen各向異性參數(shù)場(chǎng),在此基礎(chǔ)上按照常規(guī)層析建模思路進(jìn)行了多輪速度迭代(圖1a),應(yīng)用此速度模型進(jìn)行OVT域疊前深度域偏移得到OVG,部分OVG上存在一定的同相軸隨方位周期性抖動(dòng)的現(xiàn)象。以往認(rèn)為這種現(xiàn)象可能是地下裂縫引起的速度隨方位的變化。但是從YA地區(qū)地質(zhì)情況看,YA地區(qū)目的層三垛組到戴南組是砂泥巖互層沉積,砂巖百分含量為30%～70%,本區(qū)泥巖裂縫普遍不發(fā)育,我們認(rèn)為這種現(xiàn)象不太可能是由地下裂縫引起的。YA地區(qū)速度場(chǎng)較復(fù)雜,地震剖面上顯示中部是HL斷裂帶,阜寧組—戴南組有高角度火成巖侵入,地層傾角變化大,中淺層小斷塊發(fā)育。因此認(rèn)為OVG上同相軸不能較好校平的主要原因是速度精度不夠。常規(guī)層析建模技術(shù)假設(shè)所有的地震道方位都是主方位,反演過程中未考慮射線方位信息,導(dǎo)致所建立的速度模型不能滿足所有方位地震數(shù)據(jù)的要求,僅能使某些方位的道校平,另外一些方位的道上翹或者下拉,從而影響斷層成像精度。為了解決上述問題,充分利用YA高密度三維數(shù)據(jù)方位寬的優(yōu)勢(shì),提高復(fù)雜小斷層成像精度,本文采用寬方位層析建模技術(shù)(圖1b)進(jìn)行速度建模。

制定寬方位層析建模策略要考慮地震資料方位角特點(diǎn)。YA高密度三維縱向覆蓋次數(shù)為20次,橫向覆蓋次數(shù)為9次,主測(cè)線方向射線密度明顯高于聯(lián)絡(luò)線方向射線密度,導(dǎo)致主測(cè)線方向和聯(lián)絡(luò)線方向的剩余曲率拾取密度和質(zhì)量存在較大的差異。為了解決這個(gè)問題,將OVG進(jìn)行方位角分組來達(dá)到每個(gè)方位角組拾取權(quán)重基本相同的目的。方位角分組既要保證每一個(gè)方位角組內(nèi)有足夠的覆蓋次數(shù)以確保剩余曲率拾取可靠,又要確保不同方位角組之間覆蓋次數(shù)基本相當(dāng),使每個(gè)方位角組數(shù)據(jù)在反演過程中共同發(fā)揮作用。如果將地震數(shù)據(jù)方位角按照每組30°等分為6組(圖3),會(huì)造成聯(lián)絡(luò)線方向的覆蓋次數(shù)(圖3c 和圖3d)明顯低于主測(cè)線方向的覆蓋次數(shù)(圖3a 和圖3f)。本文根據(jù)YA高密度三維地震資料方位角特點(diǎn)對(duì)方位角分組參數(shù)進(jìn)行適當(dāng)調(diào)整,一方面減少主測(cè)線方向的方位角寬度,增加聯(lián)絡(luò)線方向的方位角寬度,另一方面將方位角組的個(gè)數(shù)由6個(gè)調(diào)整成5個(gè)(圖4)。調(diào)整后,實(shí)現(xiàn)了同組內(nèi)覆蓋次數(shù)較高并且不同組之間覆蓋次數(shù)相當(dāng)?shù)哪康?。將同一個(gè)OVG上剩余曲率拾取結(jié)果按照方位角組和偏移距兩個(gè)關(guān)鍵字進(jìn)行排序顯示(圖5),同一個(gè)反射層位不同方位角組拾取的剩余曲率有一定的差異,這反映了同一個(gè)成像點(diǎn)上同一地層來自不同方位的射線旅行時(shí)的差異。

圖3 偏移后OVT道集按方位角等分6組時(shí)覆蓋次數(shù)屬性a 方位角15°覆蓋次數(shù)30次; b 方位角45°覆蓋次數(shù)30次; c 方位角75°覆蓋次數(shù)15次; d 方位角105°覆蓋次數(shù)15次; e方位角135°覆蓋次數(shù)30次; f 方位角165°覆蓋次數(shù)30次

圖4 偏移后OVT道集按方位角分5組時(shí)覆蓋次數(shù)屬性a 方位角18°覆蓋次數(shù)40次; b 方位角54°覆蓋次數(shù)40次; c 方位角90°覆蓋次數(shù)40次; d 方位角126°覆蓋次數(shù)40次; e 方位角162°覆蓋次數(shù)40次

圖5 OVG上分5個(gè)方位角組分別進(jìn)行剩余曲率拾取a 中心方位角18°; b 中心方位角54°; c 中心方位角90°; d 中心方位角126°; e 中心方位角162°

采用相同的Thomsen各向異性參數(shù)場(chǎng)、地質(zhì)信息約束條件和反演參數(shù),分別應(yīng)用常規(guī)層析建模技術(shù)(圖6a)和寬方位層析建模技術(shù)(圖6b)反演得到各自的速度模型,將各自的速度模型和地震數(shù)據(jù)疊合顯示,寬方位層析建模技術(shù)反演速度細(xì)節(jié)刻畫更清楚,速度變化趨勢(shì)和地震剖面吻合程度更高。寬方位層析建模技術(shù)一定程度上消除了OVG上旅行時(shí)隨方位周期性抖動(dòng)的現(xiàn)象(圖7b),提高了同相軸校平程度。說明提高速度模型精度可以較好解決YA高密度三維OVG上出現(xiàn)的旅行時(shí)周期性抖動(dòng)的問題。

圖6 常規(guī)層析建模技術(shù)(a)和寬方位層析建模技術(shù)(b)反演得到的速度模型對(duì)比

圖7 應(yīng)用常規(guī)層析建模技術(shù)(a)和寬方位層析建模技術(shù)(b)得到的OVG道集對(duì)比(粉紅線代表方位角)

寬方位層析技術(shù)反演的速度精度高,提高了同相軸校平質(zhì)量和小斷層成像精度。偏移成果剖面上紅圈內(nèi)HL斷裂帶附近的小斷塊細(xì)節(jié)刻畫更清晰(圖8),小斷層的識(shí)別精度提高。應(yīng)用YA高密度三維處理成果資料落實(shí)紅圈內(nèi)多個(gè)復(fù)雜小斷塊(圖9),部署的2口探井獲得平均15t/d以上的高產(chǎn)油流,先后成功實(shí)施多口高效開發(fā)井。依托永安高密度三維處理成果資料新增商業(yè)可采儲(chǔ)量占2020年江蘇油田新增商業(yè)可采儲(chǔ)量的16.3%。

圖8 應(yīng)用常規(guī)層析建模技術(shù)(a)和寬方位層析建模技術(shù)(b)偏移成像效果對(duì)比

圖9 老解釋成果(a)和新解釋成果(b)對(duì)比

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

寬方位層析建模技術(shù)是針對(duì)寬(全)方位地震數(shù)據(jù)的建模技術(shù),與常規(guī)層析建模技術(shù)相比,寬方位層析建模技術(shù)在反演過程中考慮了方位角信息,所建的速度模型精度更高,可以在一定程度上解決OVG道集上出現(xiàn)的旅行時(shí)周期性抖動(dòng)現(xiàn)象,提高復(fù)雜小斷塊成像精度。寬方位層析建模技術(shù)在YA高密度三維工區(qū)的應(yīng)用驗(yàn)證了該技術(shù)在提高寬方位地震資料復(fù)雜小斷塊成像精度方面的有效性,可在類似地區(qū)推廣應(yīng)用。同時(shí)需要注意,應(yīng)用該技術(shù)之前需要考慮工區(qū)地質(zhì)特點(diǎn),如果OVT域道集上的同相軸變化反映的是由地下裂縫存在引起的速度隨方位角的變化,采用本文方法是不合理的。