魏欣偉,薛 姣,羅 霞

(1.中國石油化工股份有限公司勝利油田分公司物探研究院,山東東營257000;2.中國地質(zhì)大學(xué)(武漢)地球物理與空間信息學(xué)院,湖北武漢430074;3.地球內(nèi)部多尺度成像湖北省重點(diǎn)實(shí)驗(yàn)室,湖北武漢430074)

裂縫是油氣儲(chǔ)層的儲(chǔ)集空間和滲濾通道,非常規(guī)致密儲(chǔ)層存在孔隙度低、非均質(zhì)性強(qiáng)等特點(diǎn),尋找有效的裂縫發(fā)育帶對致密儲(chǔ)層的勘探和開發(fā)尤為重要。由于地層上覆載荷的壓實(shí)作用,儲(chǔ)層中廣泛發(fā)育中、高角度裂縫,造成了巖石的方位各向異性。通常把垂直定向排列裂縫儲(chǔ)層等效為HTI(具有水平對稱軸的橫向各向同性)介質(zhì),利用地震波在HTI介質(zhì)中傳播的方位各向異性信息可以進(jìn)行裂縫預(yù)測。

地震縱波在HTI介質(zhì)中傳播的方向特性表現(xiàn)為:對于固定偏移距,地震屬性F與測線方位角φ的關(guān)系為F(φ)=A+Bcos[2(φ-φ0)],其中系數(shù)B反映了裂縫發(fā)育強(qiáng)度,角度φ0包含了裂縫走向信息[1]。RüGER等[2]在HTI介質(zhì)線性近似反射系數(shù)研究的基礎(chǔ)上,提出利用AVO(反射振幅隨偏移距變化)梯度、動(dòng)校正(NMO)速度和橫波分裂分析計(jì)算裂縫走向和各向異性參數(shù)。對裂縫方位各向異性比較敏感的地震屬性有NMO速度[2]、反射振幅[3]、縱波旅行時(shí)[4]、頻率和衰減[5]等。國內(nèi)的多個(gè)碳酸鹽巖、火山巖裂縫儲(chǔ)層綜合分析實(shí)例證明了利用P波地震屬性隨方位角變化的橢圓擬合方法進(jìn)行裂縫檢測的有效性[6-8]。宋維琪等[9]利用傅里葉級數(shù)形式的P波方位各向異性反射系數(shù)擬合進(jìn)行裂縫各向異性估計(jì),進(jìn)一步提高了裂縫檢測的準(zhǔn)確性。

隨著裂縫等效介質(zhì)模型各向異性理論研究的發(fā)展,以及寬方位地震采集和處理技術(shù)研究的深入,基于AVOA數(shù)據(jù)反演裂縫參數(shù)成為裂縫儲(chǔ)層預(yù)測研究的重點(diǎn)之一。AVOA裂縫參數(shù)反演中常用的裂縫等效介質(zhì)模型是Schoenberg線性滑動(dòng)模型[10],利用法向和切向裂縫弱度來表征裂縫,具有類似于HTI介質(zhì)的彈性剛度矩陣形式。SHAW等[11]基于線性滑動(dòng)裂縫等效介質(zhì)模型推導(dǎo)了反射系數(shù)與裂縫弱度之間的關(guān)系,提出利用AVOA數(shù)據(jù)反演法向和切向裂縫弱度的方法。基于裂縫介質(zhì)巖石物理模型,利用疊前地震數(shù)據(jù)AVOA進(jìn)行裂縫參數(shù)反演能夠充分利用疊前方位地震信息,獲得裂縫等效介質(zhì)模型參數(shù),實(shí)現(xiàn)裂縫儲(chǔ)層定量描述[12-14]。XUE等[15]提出了一種新流體因子表達(dá)式,以及方位各向異性反射振幅與法向和切向裂縫弱度的矩陣表達(dá)式,利用疊前AVOA技術(shù)反演了裂縫密度和新的流體因子。PAN等[16]和CHEN等[17]利用貝葉斯理論McMc(馬爾科夫鏈蒙特卡洛)算法反演法向和切向裂縫弱度。

疊前AVOA數(shù)據(jù)反演在實(shí)際地震資料中的應(yīng)用得益于寬方位地震數(shù)據(jù)采集和處理的發(fā)展。OVT技術(shù)[18-20]是一種有效的寬方位三維地震數(shù)據(jù)處理技術(shù),OVT域道集具有豐富的炮檢距和方位角信息,基于OVT域地震數(shù)據(jù)可以有效地實(shí)現(xiàn)方位各向異性分析與反演研究,提高裂縫儲(chǔ)層預(yù)測的精度與可靠性[18]。目前,OVT域地震數(shù)據(jù)主要用于分方位地震屬性橢圓擬合[7,19]或Thomsen各向異性參數(shù)反演[20]。

常規(guī)AVOA裂縫參數(shù)反演方法先反演法向和切向裂縫弱度,然后根據(jù)線性滑動(dòng)模型與Hudson扁圓幣狀模型[21]的等價(jià)關(guān)系計(jì)算裂縫密度。本文基于含油和含氣條件下的裂縫等效介質(zhì)模型開展AVOA裂縫密度直接反演方法研究,建立基于OVT域疊前地震數(shù)據(jù)的AVOA裂縫密度反演技術(shù)流程,最后基于SK地區(qū)陸上勘探OVT域地震數(shù)據(jù)開展了AVOA裂縫密度反演和裂縫儲(chǔ)層預(yù)測,以驗(yàn)證方法的有效性。

1 基本原理

1.1 裂縫等效介質(zhì)模型

裂縫等效介質(zhì)模型建立了裂縫參數(shù)與彈性參數(shù)之間的關(guān)系。最簡單、最常用的裂縫等效介質(zhì)模型是Schoenberg線性滑動(dòng)模型,忽略裂縫的形狀和結(jié)構(gòu),把裂縫當(dāng)作無限薄的柔軟平面,或者以線性滑動(dòng)邊界條件表示的弱度平面。線性滑動(dòng)模型的基礎(chǔ)是Backus平均[22-23],即裂縫介質(zhì)的等效柔度矩陣等于各向同性背景介質(zhì)的柔度矩陣與裂縫柔度矩陣之和。在線性滑動(dòng)模型中裂縫被當(dāng)作一個(gè)應(yīng)力連續(xù)的位移間斷面,位移間斷和連續(xù)的應(yīng)力之間的關(guān)系是線性的,因此稱為線性滑動(dòng)模型。線性滑動(dòng)模型的剛度矩陣可以用巖石骨架的拉梅常數(shù)λ和μ,法向裂縫弱度ΔN和切向裂縫弱度ΔT來表示[11]:

C=Cb-(λ+2μ)·

(1)

式中:Cb是各向同性背景介質(zhì)彈性剛度矩陣;g是背景介質(zhì)橫波與縱波速度比的平方,g=μ(λ+2μ);r=λ/(λ+2μ)。常規(guī)的HTI介質(zhì)具有3個(gè)各向異性參數(shù)ε,γ和δ,其中ε和γ分別用于表征縱波和橫波的方位各向異性。線性滑動(dòng)模型是一種特殊的HTI介質(zhì)模型,其中,切向裂縫弱度代表橫波各向異性,與裂縫密度呈正比,法向裂縫弱度代表縱波各向異性,與裂縫密度呈正比,同時(shí)也與流體充填情況有關(guān),當(dāng)裂縫中飽和流體且流體無法流動(dòng)時(shí),法向裂縫弱度為零。

為了更加直觀地表征裂縫發(fā)育程度,引入Hudson扁圓幣狀裂縫介質(zhì)模型[21]。Hudson模型假設(shè)各向同性彈性介質(zhì)中嵌入稀疏定向排列的扁圓幣狀(或稱扁橢球狀)裂縫,模擬的是超聲波實(shí)驗(yàn)室條件下飽和巖石的性質(zhì)[24],裂縫密度e可以表示為:

(2)

式中:N/V表示單位體積內(nèi)的裂縫個(gè)數(shù);a是裂縫半徑(即橢球長軸半徑);φ是扁圓幣狀裂縫孔隙度;α是裂縫扁度(即橢球短軸與長軸之間的比值)。Hudson模型的一階等效彈性剛度矩陣為[23]:

(3)

式中:U11和U33是與裂縫性質(zhì)有關(guān)的無量綱常量。線性滑動(dòng)模型與Hudson扁圓幣狀裂縫模型對裂縫的描述有很大差別,但是又存在等價(jià)性。線性滑動(dòng)模型的彈性剛度矩陣與扁圓幣狀裂縫介質(zhì)模型的一階近似具有相同的形式[23]。

(4)

對于流體充填裂縫,流體剪切模量μf=0,切向裂縫弱度簡化為[23]:

(5)

特別地,當(dāng)裂縫中飽含油、水時(shí),流體體積模量Kf≠0,裂縫的扁度α→0,法向裂縫弱度為[23]:

ΔN=0

(6)

對于含氣裂縫介質(zhì),流體體積模量Kf=0,法向裂縫弱度簡化為[23]:

(7)

垂直定向排列裂縫造成了巖石性質(zhì)的各向異性,利用常規(guī)橢圓擬合方法得到的橢圓扁率是一種定性的裂縫發(fā)育強(qiáng)度表示方法,基于裂縫巖石物理模型進(jìn)行裂縫密度和裂縫弱度反演則是一種定量裂縫參數(shù)反演方法。切向裂縫弱度與裂縫密度呈正相關(guān)關(guān)系,法向裂縫弱度與裂縫密度呈正相關(guān),與充填流體的體積模量呈負(fù)相關(guān)關(guān)系。

1.2 裂縫介質(zhì)反射系數(shù)

裂縫等效介質(zhì)分界面PP波反射系數(shù)可以分為與反射界面兩側(cè)各向同性背景介質(zhì)物性差異有關(guān)的反射系數(shù)各向同性部分Riso(θ)以及與定向排列裂縫有關(guān)的反射系數(shù)各向異性部分Rani(θ,φ)[14]:

RPP(θ,φ)=Riso(θ)+Rani(θ,φ)

(8a)

(8b)

Rani(θ,φ)=-[g(1-2g)cos2(φ-φ0)sin2θ(1+tan2θ)+g2cos4(φ-φ0)sin2θtan2θ]δΔN+gcos2(φ-φ0)sin2θ[1-sin2(φ-φ0)tan2θ]δΔT

(8c)

式中:θ和φ分別表示入射角和測線方位角;φ0表示裂縫對稱軸方向;vP是縱波速度;Z和G分別表示縱波阻抗和橫波模量;符號“-”和δ分別表示反射界面兩側(cè)參數(shù)的平均值和差值,例如δΔN=ΔN2-ΔN1。當(dāng)入射角θ<30°時(shí),忽略高階項(xiàng)sin2θtan2θ,則有:

(9a)

Rani(θ,φ)=-g(1-2g)cos2(φ-φ0)·sin2θδΔN+gcos2(φ-φ0)sin2θδΔT

(9b)

對于含油裂縫介質(zhì),法向裂縫弱度ΔN=0,反射系數(shù)公式簡化為:

(10)

對于含氣裂縫介質(zhì),根據(jù)法向裂縫弱度與裂縫密度之間的關(guān)系(公式(7)),反射系數(shù)公式簡化為:

(11)

1.3 裂縫參數(shù)AVOA反演原理

假設(shè)已知裂縫對稱軸方向φ0,兩個(gè)不同測線方位φ1和φ2的反射系數(shù)之差為:

ΔR(θ)=R(θ,φ2)-R(θ,φ1)=Ce(θ)Δe

(12)

對于含油裂縫,模型參數(shù)正演系數(shù)為:

(13)

對于含氣裂縫,模型參數(shù)正演系數(shù)為:

(14)

對于K+1層裂縫介質(zhì),引入地震子波矩陣,將N個(gè)入射角不同方位的地震道集差值ΔS=S(φ2)-S(φ1)與裂縫參數(shù)之間的關(guān)系表示成矩陣的形式:

(15)

式中:S(θi)=[S(θi,t1)S(θi,t2)S(θi,tK)]T表示第i個(gè)入射角對應(yīng)的地震數(shù)據(jù);W(θi)是子波褶積矩陣;Ce(θi)是由裂縫參數(shù)正演系數(shù)組成的對角矩陣:

Ce(θi)=

(16)

假設(shè)已知裂縫走向φ0,對于含油裂縫介質(zhì),利用公式(13)計(jì)算Ce(θi)矩陣對角線上的正演系數(shù);對于含氣裂縫,利用公式(14)計(jì)算Ce(θi)矩陣對角線上的正演系數(shù)。正演方程(公式(15))是一個(gè)線性方程,可以簡寫為:

d=Gm

(17)

式中:G是由子波矩陣,模型參數(shù)正演系數(shù)矩陣和一階差分算子組成的線性正演算子;m=[e(t0),e(t1),,e(tK)]T表示裂縫密度向量。

假設(shè)地震數(shù)據(jù)噪聲相互獨(dú)立且滿足高斯分布,具有相同的方差,數(shù)據(jù)的似然函數(shù)也滿足高斯分布,將數(shù)據(jù)似然函數(shù)與模型先驗(yàn)分布結(jié)合起來得到后驗(yàn)概率分布,將后驗(yàn)概率最大化,得到目標(biāo)函數(shù):

J(m)=(d-Gm)T(d-Gm)+μR(m)

(18)

式中:R(m)是由模型參數(shù)先驗(yàn)分布得到的正則化項(xiàng),通過改變系數(shù)μ改變模型約束條件在反演求解過程中占的比重。目標(biāo)函數(shù)對模型參數(shù)m的導(dǎo)數(shù)為0即可求得目標(biāo)函數(shù)的最小值,得到反演問題的解。本文采用Tikhonov正則化,并加入非負(fù)約束進(jìn)行迭代反演。

1.4 基于OVT域數(shù)據(jù)的AVOA反演流程

常規(guī)AVAZ(振幅隨方位角變化)分析或裂縫介質(zhì)各向異性橢圓擬合方法中,通常通過定義疊加量版,對OVT道集數(shù)據(jù)部分偏移距、部分方位角進(jìn)行疊加,從而進(jìn)一步提高數(shù)據(jù)信噪比。而在AVOA反演中,需要利用偏移速度從OVT道集中抽取分方位-入射角道集,即從偏移距轉(zhuǎn)換到(入射角)角度域,加強(qiáng)振幅能量和穩(wěn)定性。

AVOA反演利用的是相同反射點(diǎn)的振幅隨方位角和入射角的變化,而疊前道集中可能存在方位各向異性時(shí)差,以及遠(yuǎn)偏移距數(shù)據(jù)仍然存在動(dòng)校正未校平的殘余時(shí)差。方位各向異性和動(dòng)校正誤差引起的方位時(shí)差會(huì)影響反演結(jié)果的準(zhǔn)確性,因此,需要進(jìn)一步對AVOA道集做時(shí)差校正。利用滑動(dòng)時(shí)窗波形匹配剩余時(shí)差校正技術(shù),將疊后數(shù)據(jù)作為參考道,選定時(shí)窗內(nèi)某一道集的所有道與參考道進(jìn)行互相關(guān),求取每道最大相關(guān)系數(shù)對應(yīng)的時(shí)移量用于時(shí)差校正。

AVOA反演中需要已知裂縫走向或裂縫對稱軸方向,首先利用井信息或者地震數(shù)據(jù)的橢圓擬合或AVAZ反演獲取裂縫走向,繼而進(jìn)行AVOA裂縫參數(shù)反演[15-16,25-26]。

基于OVT道集數(shù)據(jù)的AVOA裂縫參數(shù)反演的流程如圖1所示。

圖1 基于OVT數(shù)據(jù)的AVOA反演流程

2 模型試算

建立裂縫等效介質(zhì)模型,其中各向同性背景介質(zhì)縱橫波速度、密度和裂縫密度參數(shù)分別為R1,R2和R3,如圖2所示。該模型中存在3段裂縫儲(chǔ)層,裂縫走向?yàn)楸睎|向45°,裂縫密度均為0.1。假設(shè)裂縫扁度α=10-6,法向裂縫弱度隨充填流體體積模量的增大而減小(圖3),當(dāng)裂縫中充填油或水時(shí)可忽略法向弱度,切向裂縫弱度不隨充填流體變化。

圖2 裂縫介質(zhì)模型縱橫波速度(a)、密度(b)和裂縫密度(c)

圖3 裂縫弱度隨充填物體積模量的變化(紅圈表示含氣、含油和含水裂縫介質(zhì)的法向裂縫弱度)

利用本文提出的AVOA裂縫密度反演算法對裂縫儲(chǔ)層模型的合成地震記錄進(jìn)行試算。假設(shè)裂縫儲(chǔ)層中飽和充填油,利用Zoeppritz方程計(jì)算反射系數(shù),并將反射系數(shù)與30Hz雷克子波進(jìn)行褶積得到合成疊前地震記錄,然后分別添加20%,50%和100%的隨機(jī)噪聲,得到的合成地震記錄信噪比分別為RS/N=5,RS/N=2和RS/N=1。首先,對合成疊前地震數(shù)據(jù)進(jìn)行AVAZ裂縫走向反演,圖4顯示,當(dāng)信噪比大于等于2時(shí),儲(chǔ)層段裂縫走向反演結(jié)果與裂縫模型走向相吻合,非裂縫儲(chǔ)層段(裂縫密度為0的位置)裂縫走向反演受噪聲影響較大。利用與裂縫對稱軸夾角分別為0和90°的疊前地震數(shù)據(jù)(圖5和圖6)進(jìn)行AVOA裂縫密度反演,其結(jié)果如圖7所示。由于裂縫密度不可能是負(fù)值,因此在反演過程中加入非負(fù)約束條件。圖7顯示,當(dāng)信噪比大于等于2時(shí)滿足反演要求,當(dāng)信噪比為1時(shí)反演結(jié)果與裂縫模型真實(shí)值有較大差異。

圖4 不同信噪比疊前地震數(shù)據(jù)裂縫走向AVAZ反演a 不含噪聲; b 含20%噪聲; c 含50%噪聲; d 含100%噪聲

圖5 不同信噪比疊前地震數(shù)據(jù)(Ⅰ)a 不含噪聲; b 含20%噪聲

圖6 不同信噪比的疊前地震數(shù)據(jù)(Ⅱ)a 含50%噪聲; b 含100%噪聲

圖7 不同信噪比疊前地震數(shù)據(jù)裂縫密度反演結(jié)果(紅線)與模型真實(shí)值(藍(lán)線)對比a 不含噪聲; b 含20%噪聲; c 含50%噪聲; d 含100%噪聲

假設(shè)裂縫模型為干裂縫或者含氣裂縫,依然用含油情況下的AVOA裂縫密度反演方法進(jìn)行反演,則其裂縫密度反演結(jié)果小于真實(shí)值(圖8)。

圖8 干裂縫介質(zhì)模型裂縫密度反演結(jié)果(紅線)與模型真實(shí)值(藍(lán)線)對比a 不含噪聲; b 含20%噪聲; c 含50%噪聲; d 含100%噪聲

3 實(shí)際應(yīng)用

中國中東部某探區(qū)致密砂巖油氣資源潛力大,巖性致密、基質(zhì)孔隙結(jié)構(gòu)復(fù)雜,具有較強(qiáng)的非均質(zhì)性,裂縫的發(fā)育程度對致密砂巖儲(chǔ)層的控制作用非常明顯,該探區(qū)中SK工區(qū)多層系含油,裂縫較為發(fā)育,某些井試油日產(chǎn)達(dá)百噸以上。實(shí)際地震數(shù)據(jù)為SK工區(qū)陸上勘探得到的OVT道集。

3.1 OVT道集分析與預(yù)處理

SK工區(qū)OVT道集覆蓋次數(shù)為93～312次,不同位置覆蓋次數(shù)分布不均勻。方位角-偏移距變化關(guān)系分析圖(圖9)顯示OVT道集偏移距范圍為300～6000m,南北方向炮檢距分布均勻,覆蓋次數(shù)較高,然而橫縱比較小,方位分布不均勻;南北方向最大偏移距達(dá)6000m,而東西方向最大偏移距不足2500m,東西方向遠(yuǎn)偏移距信息缺失。偏移距大于2500m的部分,方位角主要分布在0～60°和120°～180°之間,如果使用偏移距大于2500m的疊前地震數(shù)據(jù)進(jìn)行AVOA反演則會(huì)產(chǎn)生較大誤差。經(jīng)過分析,有效的偏移距范圍是300～2500m,其覆蓋次數(shù)和方位角分布較為均衡。

圖9 方位角-偏移距分析

對OVT數(shù)據(jù)進(jìn)行分方位-入射角道集抽取時(shí),需要考慮各面元覆蓋次數(shù)以及方位角和入射角個(gè)數(shù)之間的平衡。理論上講,方位角和入射角個(gè)數(shù)越多,反演的精度越高,但隨著方位角和入射角的增多,覆蓋次數(shù)變少,AVOA道集數(shù)據(jù)的信噪比逐漸降低。

雖然目的層T6和T7界面最大入射角約為50°,但由于不同方位炮檢距分布不均,在抽取分方位-入射角道集時(shí)應(yīng)充分考慮最大有效炮檢距(2500m)對應(yīng)的入射角作為最大入射角(約30°)。經(jīng)過試驗(yàn)最終確定劃分為4個(gè)方位和3個(gè)入射角,生成中心方位角分別為22.5°,67.5°,112.5°和157.5°,中心入射角分別為10°,17°和24°的分方位-入射角道集(圖10a)。方位角劃分范圍分別為0～45°,45°～90°,90°～135°和135°～180°,入射角劃分范圍分別為6°～13°,13°～20°,和20°～27°。選用合適的滑動(dòng)時(shí)窗參數(shù)進(jìn)行剩余時(shí)差校正后(圖10b),方位各向異性和動(dòng)校正誤差引起的時(shí)差被校正,同相軸被拉平。

圖10 時(shí)差校正前(a)、后(b)分方位-入射角道集

3.2 AVOA裂縫參數(shù)反演結(jié)果

在疊前OVT道集資料分析的基礎(chǔ)上,抽取分方

位-入射角道集,利用滑動(dòng)時(shí)窗波形匹配時(shí)差校正技術(shù)消除方位各向異性和動(dòng)校正剩余時(shí)差,然后求取相互垂直的角道集之差:

(19)

式中:S(θi,φj)表示第i個(gè)入射角、第j個(gè)方位角對應(yīng)的地震道;對于SK實(shí)際地震數(shù)據(jù),S(θi,φ3)-S(θi,φ1)表示入射角為θi時(shí)第3個(gè)方位(φ3=112.5°)與第1個(gè)方位(φ1=22.5°)地震道之間的差異。首先,利用AVAZ反演估計(jì)裂縫走向,圖11顯示研究區(qū)目的層(T6上10ms到T6下10ms)裂縫走向主要為近東西向,裂縫走向與斷裂帶的延展方向一致。然后,應(yīng)用本文提出的AVOA裂縫參數(shù)反演方法估計(jì)裂縫密度,達(dá)到裂縫儲(chǔ)層預(yù)測的目的。測井、鉆井、巖心和試油資料表明SK工區(qū)多層系含油,為裂縫型含油致密砂巖儲(chǔ)層,因此將含油裂縫的模型參數(shù)正演系數(shù)(公式(13))代入模型正演方程(公式(15))中,然后利用目標(biāo)函數(shù)(公式(18))最小化進(jìn)行AVOA反演計(jì)算。裂縫密度沿T6層切片(圖12)顯示在研究區(qū)西北部yi177井到y(tǒng)i178井之間裂縫密度較大,為裂縫發(fā)育帶。裂縫的分布主要受斷裂帶控制,裂縫的分布與斷裂系統(tǒng)的分布具有一定相關(guān)性,裂縫密度較大的位置位于斷裂帶周圍和斷裂帶交會(huì)處。

圖11 沿T6層裂縫走向切片(藍(lán)線走向表示裂縫走向,藍(lán)線長度表示各向異性程度)

圖12 沿T6層裂縫密度切片

將測井資料解釋成果與AVOA裂縫密度反演結(jié)果進(jìn)行對比分析,驗(yàn)證基于AVOA裂縫密度反演的裂縫儲(chǔ)層預(yù)測效果。yi176井T6～T7含油層段為粉細(xì)砂巖,綜合解釋為油層,常規(guī)測井曲線上具有“高聲波時(shí)差、高中子、低電阻、低密度和自然電位負(fù)異常”的特征。圖13是過yi176井(線號702)裂縫密度反演剖面,其中,粉紫色曲線顯示了正交多極子聲波陣列測井得到的各向異性參數(shù),黑色曲線是井旁地震道的裂縫密度反演結(jié)果。井旁地震道的裂縫密度反演結(jié)果(黑色曲線)與聲波陣列各向異性(粉紫色曲線)異常相吻合,裂縫密度反演結(jié)果顯示在T6到T7層中間及T7層上部裂縫較為發(fā)育,測井解釋為裂縫儲(chǔ)層。裂縫密度反演結(jié)果顯示裂縫儲(chǔ)層在橫向上呈似層狀和團(tuán)塊狀分布。

圖13 過yi176井裂縫密度剖面

4 結(jié)論

本文基于裂縫等效介質(zhì)模型,分別推導(dǎo)了含油和含氣情況下反射系數(shù)與裂縫密度之間的關(guān)系,在此基礎(chǔ)上提出了基于疊前AVOA反演的裂縫密度估計(jì)方法和基于OVT域疊前數(shù)據(jù)的裂縫密度反演技術(shù)流程。SK工區(qū)裂縫參數(shù)反演結(jié)果與測井資料解釋成果相吻合,顯示了OVT域AVOA裂縫參數(shù)反演的有效性及其在SK工區(qū)裂縫儲(chǔ)層預(yù)測中的適用性。基于OVT域地震數(shù)據(jù)的AVOA裂縫參數(shù)反演結(jié)果能夠定量刻畫裂縫分布特征,為SK工區(qū)致密砂巖含油裂縫儲(chǔ)層預(yù)測和開發(fā)提供了重要依據(jù)。