李春鵬,王 哲,李愛山,袁 峰

(1.中海石油 國(guó)際能源服務(wù)(北京)有限公司，北京 100028；2.中國(guó)石油大學(xué)(華東)，青島 266580；3.中海油能源發(fā)展股份有限公司 工程技術(shù)分公司，天津 300452)

0 前言

裂縫作為地殼中一種普遍的構(gòu)造現(xiàn)象，廣泛存在于各類巖石中。到目前為止，已在砂巖、泥頁(yè)巖和碳酸鹽巖，甚至火成巖等各類巖石的裂縫性儲(chǔ)層中獲得了大量的工業(yè)油氣流。實(shí)際地層中裂縫呈一種至多種形態(tài)，包括近水平定向分布裂縫(VTI介質(zhì))、近垂直定向分布裂縫(HTI介質(zhì))、傾斜定向分布裂縫(TTI介質(zhì))和正交裂縫(OA介質(zhì))等。對(duì)于近垂直定向分布裂縫(HTI介質(zhì))，可以利用各向異性梯度進(jìn)行預(yù)測(cè)，各向異性梯度越大，裂縫密度也越大[1-3]。

各向異性梯度可以通過方位地震數(shù)據(jù)反演得到[4]，寬方位比窄方位包含了更多的各向異性信息，實(shí)際生產(chǎn)中一般都利用寬方位地震數(shù)據(jù)反演各向異性參數(shù)[5-6]。不同于三維的疊后數(shù)據(jù)和四維的疊前數(shù)據(jù)，方位地震數(shù)據(jù)是五維數(shù)據(jù)，目前主要依據(jù)炮檢距向量片(OVT)技術(shù)處理寬方位地震數(shù)據(jù)，OVT技術(shù)是一種新穎的疊前數(shù)據(jù)的編排方式，基于OVT數(shù)據(jù)域處理可有效改善寬方位數(shù)據(jù)處理效果，且OVT域偏移結(jié)果含有豐富的方位各向異性信息，能實(shí)現(xiàn)寬方位數(shù)據(jù)的保真處理并提高疊前地震道集的質(zhì)量，是寬方位三維地震數(shù)據(jù)的有效處理技術(shù)[7]。

地震反演的本質(zhì)是求解線性方程組[8-9]，線性方程組能否得到良態(tài)的解，取決于它的系數(shù)矩陣的條件數(shù)。條件數(shù)越小，線性方程組的解越趨于良態(tài)；反之，越趨于病態(tài)。對(duì)于無約束地震反演方程，其系數(shù)矩陣是子波矩陣，然而子波矩陣的條件數(shù)特別大，所得到的反演結(jié)果是病態(tài)的。

為了解決這一問題，一般會(huì)對(duì)反演方程做加法，即加入高斯分布、柯西分布等數(shù)學(xué)約束，同時(shí)也會(huì)加入測(cè)井低頻約束[10-13]。事實(shí)證明，地震反演方程系數(shù)矩陣的條件數(shù)在加入數(shù)學(xué)約束和測(cè)井低頻約束之后顯著降低了，反演結(jié)果由病態(tài)變?yōu)榱紤B(tài)。

實(shí)際工區(qū)一般都有測(cè)井縱、橫波數(shù)據(jù)，然而各向異性數(shù)據(jù)卻不常有。沒有各向異性數(shù)據(jù)也就沒有各向異性低頻約束，這樣就難以對(duì)各向異性反演方程做加法。各向異性梯度反演方程會(huì)因缺少測(cè)井各向異性低頻約束而難以降低反演方程條件數(shù)，從而難以得到準(zhǔn)確的反演結(jié)果[9]。

各向異性梯度反演方程因缺少各向異性測(cè)井信息難以做加法，應(yīng)對(duì)各向異性梯度反演方程進(jìn)行其他改造。筆者基于這種思路，提出了對(duì)各向異性梯度反演方程做減法的方法，即穩(wěn)定的各向異性梯度無井約束反演方法。該方法通過在反演方程中消除子波矩陣和降低反演方程的維數(shù)兩個(gè)減法操作，來降低反演方程系數(shù)矩陣的條件數(shù)，從而得到穩(wěn)定的反演結(jié)果。

1 常規(guī)各向異性梯度反演方法

1.1 各向異性梯度的數(shù)學(xué)表征

高角度定向排列裂縫地層可以等效為HTI介質(zhì)，Ruger[4]推導(dǎo)了上、下HTI介質(zhì)對(duì)稱軸一致的方位縱波反射系數(shù)近似公式，Downton[14]在Ruger基礎(chǔ)上將該近似公式進(jìn)行了線性化，即

(1)

圖1 地震反演流程圖Fig.1 Work flow of seismic inversion

(2)

通過式(1)可以建立反演方程，據(jù)此可以反演得到C和D，然后根據(jù)式(2)可以計(jì)算各向異性梯度。

1.2 常規(guī)各向異性梯度反演方程

方位地震數(shù)據(jù)的最小二乘反演方程為式(3)。

(WG)T(WG)m=(WG)TS

(3)

式中：S為方位地震數(shù)據(jù);W表為子波矩陣;m=[ABCD]T為待反演參數(shù);G為反射系數(shù)正演算子矩陣;R為反射系數(shù)。

(WG)T(WG)是反演方程的系數(shù)矩陣，該矩陣條件數(shù)極大，反演結(jié)果不穩(wěn)定。一般都是對(duì)該系數(shù)矩陣做加法，通過加入高斯分布、柯西分布等數(shù)學(xué)約束和測(cè)井低頻約束來降低系數(shù)矩陣的條件數(shù)，以得到穩(wěn)定的反演結(jié)果。但是實(shí)際工區(qū)有時(shí)會(huì)缺少各向異性梯度先驗(yàn)分布信息和測(cè)井約束信息，因此難以通過常規(guī)反演方程得到真實(shí)的地下地層各向異性梯度。

2 各向異性梯度無約束反演方法

常規(guī)反演方程做加法的約束項(xiàng)的數(shù)學(xué)作用在于減小反演方程系數(shù)矩陣的條件數(shù)，以得到良態(tài)的反演結(jié)果。其實(shí)除了增加約束項(xiàng)，還有消除子波影響和降低反演方程維度的方法來降低反演方程系數(shù)矩陣的條件數(shù)，該方法分為兩步(圖1)。

1)改造反演方程，方法是子波移位到反演結(jié)果中。假設(shè)方位地震數(shù)據(jù)是方位反射系數(shù)，則方位地震數(shù)據(jù)正演方程如式(5)所示，式(3)和式(5)的最小二乘反演方程如式(6)和式(7)所示。

S=Gm'

(5)

m=[(WG)TWG]-1(WG)TS

(6)

m'=[GTG]-1GTS

(7)

2)疊后去子波，通過地震反褶積將第一步得到的wC和wD數(shù)據(jù)體轉(zhuǎn)換為C和D數(shù)據(jù)體，然后通過式(1)～式(5)可以計(jì)算各向異性梯度大小和各向異性梯度方向。但是常規(guī)反褶積方法是針對(duì)每一個(gè)采樣點(diǎn)計(jì)算的，計(jì)算量較大且不穩(wěn)定。一般采用正則化算法改善反褶積算子，但是這種常規(guī)反褶積方法得到的反射系數(shù)不是稀疏的，不能真正地反映下地層界面的真實(shí)情況。筆者依據(jù)最少反射層地震反演方法[15]，先搜索目的層地下界面位置，只在搜索的界面位置處建立反褶積方程，這樣可以降低各向異性梯度地震記錄反褶積方程維數(shù)，進(jìn)一步也會(huì)降低反褶積方程條件數(shù)，得到穩(wěn)定的反演結(jié)果。

最少反射層稀疏脈沖反演的關(guān)鍵性技術(shù)是:首先利用模擬退火法尋找反射層位置，然后利用線性優(yōu)化方法解目標(biāo)函數(shù)，最終得到準(zhǔn)確的反演結(jié)果[1,6]。

無約束各向異性梯度反演方法的關(guān)鍵是第一步子波移位，該步驟已經(jīng)能得到待反演參數(shù)，即使反演參數(shù)含有子波項(xiàng)，也能體現(xiàn)地下地層的真實(shí)情況。另外該步驟簡(jiǎn)單易操作，只需要把反演方程進(jìn)行相應(yīng)改造，即能得到穩(wěn)定的反演結(jié)果。第二步是對(duì)第一步的補(bǔ)充，其關(guān)鍵是要能找到地下地層界面的位置，它受限于地震分辨率，一般適用于較厚的地層。實(shí)際操作中定義界面數(shù)目是同相軸的兩倍，能得到較好的結(jié)果。

綜上所述，無約束各向異性梯度反演方法和常規(guī)反演方法,對(duì)反演方程的改造存在較大的不同。無約束反演方法對(duì)反演方程做“減法”，減去子波項(xiàng)和非界面采樣點(diǎn)數(shù)，以降低反演方程條件數(shù)。而常規(guī)反演方法是對(duì)反演方程做“加法”，加入統(tǒng)計(jì)約束和測(cè)井低頻約束。雖然二者手段不同，但是目的都是增強(qiáng)反演方程的穩(wěn)定性，以獲得準(zhǔn)確的反演結(jié)果。但是這并不代表無約束反演方法可以代替常規(guī)反演方法，無約束反演得到的是相對(duì)結(jié)果，常規(guī)反演得到的是絕對(duì)結(jié)果且有更多的約束更符合地質(zhì)規(guī)律。無約束反演更適用于缺少測(cè)井資料而難以提供低頻約束的情況，是地震反演的權(quán)宜之計(jì)，若資料允許，還是建議優(yōu)先使用常規(guī)反演方法。

3 模型試算

為了驗(yàn)證以上HTI介質(zhì)各向異性梯度反演方法的有效性，引用四層三維模型[1-2,9](圖2(a))，其中第三層介質(zhì)是近垂直定向分布的裂縫介質(zhì)，第一層、第二層和第四層是各向同性介質(zhì)。第三層裂縫介質(zhì)的背景彈性參數(shù)和物性參數(shù)取自于A井的測(cè)井信息(圖3)，據(jù)此可以得到裂縫密度平面分布、裂縫走向平面分布和裂縫走向玫瑰圖(圖4)。根據(jù)以上模型數(shù)據(jù)可以建立HTI介質(zhì)彈性矩陣，再根據(jù)HTI介質(zhì)彈性波反射透射方程[16]計(jì)算方位反射系數(shù)[16]，模型正演采用相移30°后的50 Hz雷克子波，入射角分別是15°、25°和35°，方位角分別是0°、30°、60°、90°、120°、150°和180°。圖2(b)顯示了inline75，xline50處，入射角是35°，信噪比是10的方位縱波地震記錄，圖2(b)中最上面的同相軸表示模型第一層和第二層的各向同性界面的方位地震記錄，該地震記錄不隨方位角變化而變化；圖2(b)中最下面的兩個(gè)同相軸分別表示第三層高角度裂縫介質(zhì)的頂和底的方位地震記錄，該地震記錄隨方位角變化而變化。

各向異性梯度“減法”反演可以分成兩步。

1)假設(shè)方位地震數(shù)據(jù)是方位反射系數(shù)，使得反演方程隱去子波項(xiàng)，則式(7)反演矩陣GTG條件數(shù)約為172。但若采用式(6)的反演方程，則反演矩陣(WG)T(WG)的條件數(shù)約為7.3×1017，遠(yuǎn)大于式(7)計(jì)算的反演矩陣條件數(shù)，可見式(7)的反演方程穩(wěn)定性更高。單道反演結(jié)果見圖5(a)，圖5中反演的自激自收反射地震記錄wA和各向異性梯度地震記錄wB、wC、wD與真實(shí)值吻合的較好，說明反演結(jié)果比較可信。另外，wC和wD平方根可以表示各向異性梯度地震記錄wBani大小，不難發(fā)現(xiàn)wC和wD平方根在各向同性界面處幾乎為零，在裂縫介質(zhì)頂、底界面處有值，說明wBani能夠指示地層裂縫發(fā)育情況。

2)通過地震反褶積將第一步得到的wC和wD數(shù)據(jù)體轉(zhuǎn)換為C和D數(shù)據(jù)體，但是以上模型的反褶積算子wTw的條件數(shù)約為1.6×1018，反褶積方程極不穩(wěn)定。因此筆者采用最少反射層假設(shè)反演C和D數(shù)據(jù)體，單道反演結(jié)果見圖5(b)，其中C和D的平方根可以表示各向異性梯度Bani大小，不難發(fā)現(xiàn)C和D的平方根在裂縫地層頂、底界面處有值，其余地層分界面處幾乎為零，說明最少反射層假設(shè)反演方法能準(zhǔn)確體現(xiàn)裂縫地層的發(fā)育情況。

將穩(wěn)定各向異性梯度反演方法推廣到整個(gè)四層介質(zhì)模型中，比較圖6(a)和圖4(a)可以發(fā)現(xiàn)，反演的各向異性梯度大小和裂縫密度比較吻合，說明反演的各向異性梯度基本可以反映裂縫地層裂縫密度的平面展布；比較圖6(b)和圖4(b)可以發(fā)現(xiàn),反演的各向異性梯度方向與裂縫地層裂縫走向基本一致或垂直，另外圖6(c)中各向異性梯度方向玫瑰圖中有兩個(gè)互相垂直的主方向，其中18°左右的主方向與圖4(c)中的裂縫走向主方向一致，因此各向異性梯度方向可以指示裂縫走向，但是存在90°不確定性，實(shí)際裂縫地層預(yù)測(cè)應(yīng)用中,需要加入裂縫走向先驗(yàn)信息以輔助識(shí)別地層真實(shí)裂縫走向。

圖2 四層介質(zhì)模型Fig.2 Four layers model (a)示意圖；(b)信噪比10的方位地震記錄

圖3 L工區(qū)A井的測(cè)井曲線Fig.3 Well A logs of block L

圖4 第三層裂縫介質(zhì)參數(shù)Fig.4 Fracture parameter of the third layer(a)裂縫密度；(b)裂縫走向;(c)走向玫瑰圖

圖5 反演結(jié)果(黑)與真實(shí)值(紅)對(duì)比Fig.5 Comparison between inversion result (red) and real result (black)(a)第一步；(b)第二步

圖6 反演結(jié)果Fig.6 Inversion result(a)各向異性梯度；(b)各向異性梯度方向；(c)方向玫瑰圖

圖7 不同方位疊后地震道集Fig.7 Post stack seismic of different azimuth(a)0°～30°;(b)30°～60°;(c)60°～90°;(d)90°～120°;(e)120°～150°;(f)150°～180°

圖8 裂縫預(yù)測(cè)結(jié)果Fig.8 Fracture prediction results(a)各向異性梯度；(b)各向異性梯度方向玫瑰圖；(c)測(cè)井結(jié)果

4 實(shí)際裂縫地層預(yù)測(cè)應(yīng)用

利用某工區(qū)過A井的實(shí)際地震數(shù)據(jù)參與反演，該工區(qū)發(fā)育大段泥頁(yè)巖裂縫地層，工區(qū)內(nèi)地震、測(cè)井資料比較完備，適合于利用方位道集進(jìn)行裂縫型儲(chǔ)層預(yù)測(cè)研究。該方位地震數(shù)據(jù)的炮檢距和方位角分布比較均勻，屬于寬方位角地震數(shù)據(jù)，因此適合于利用方位道集進(jìn)行裂縫型儲(chǔ)層預(yù)測(cè)研究。通過地震資料處理，得到了6個(gè)部分方位疊加道集，分別是0°～30°、30°～60°、60°～90°、90°～120°、120°～150°、150°～180°(圖7)，圖7中紅框內(nèi)的裂縫發(fā)育區(qū)地震振幅隨著方位角變化而變化，說明方位地震對(duì)裂縫有較好的響應(yīng)。

圖8是該工區(qū)S地層裂縫預(yù)測(cè)結(jié)果，已知A井在S地層附近存在裂縫，圖8中反演的S地層各向異性梯度在A井附近存在異常，這與已知認(rèn)識(shí)吻合。另外，反演的各向異性梯度方向與測(cè)井的裂縫走向結(jié)果一致。以上表明，這里研究的各向異性梯度反演方法能得到準(zhǔn)確的地層各向異性梯度。

筆者發(fā)展的各向異性梯度無約束反演方法，在第一步子波移位中就能得到較好的反演結(jié)果，該步簡(jiǎn)單易操作且對(duì)地震資料分辨率要求不高。第二步去子波，利用模擬退火法搜索反射界面位置，對(duì)于調(diào)諧地層效果較差，因此適用于較厚的地層，實(shí)際應(yīng)用需根據(jù)地震資料分辨率決定反演步驟。

5 結(jié)論

各向異性梯度對(duì)預(yù)測(cè)裂縫型儲(chǔ)層具有十分重要的意義，利用寬方位地震數(shù)據(jù)可以反演各向異性梯度。常規(guī)的地震反演是對(duì)反演方程做加法，但是對(duì)于沒有各向異性測(cè)井?dāng)?shù)據(jù)的工區(qū)，是難以增加測(cè)井低頻約束的。針對(duì)這種情況，可以利用筆者發(fā)展的消除子波和降低反演方程維度對(duì)反演方程做減法的方法，得到穩(wěn)定的反演結(jié)果。該方法分為兩步：①子波移位簡(jiǎn)單易操作能得到相對(duì)裂縫發(fā)育區(qū);②去子波能得到絕對(duì)裂縫發(fā)育區(qū)，但操作較為復(fù)雜且對(duì)地震資料分辨率要求較高。模型試算和實(shí)際工區(qū)應(yīng)用表明，該方法能夠得到準(zhǔn)確的各向異性梯度，據(jù)此可以較為準(zhǔn)確地預(yù)測(cè)高角度裂縫發(fā)育帶，進(jìn)一步可以為鉆探裂縫型地層提供借鑒。