孫 煒,何治亮,李玉鳳,周 雁

(1.中國石油化工股份有限公司石油勘探開發(fā)研究院構造與沉積儲層實驗室,北京100083;2.中國石油大學(北京),北京102249;3.中國石油化工股份有限公司勝利油田物探研究院,山東東營257022)

預測儲層裂縫分布特征的最終目的是得到裂縫的發(fā)育程度和方向,從而指導后續(xù)的勘探開發(fā)工作。目前,基于縱波方位各向異性的裂縫預測方法是能夠解決裂縫預測難題的技術手段之一。國內學者對此開展了較多的應用研究，如：曲壽利等[1]較早地論述了縱波方位各向異性裂縫預測在實際應用中的諸多技術細節(jié);樂紹東[2]利用縱波方位各向異性裂縫檢測技術較為準確地預測了川西坳陷上三疊統(tǒng)須家河組的裂縫分布特征;喻岳鈺等[3]利用瞬時頻域衰減屬性的方位各向異性,在碳酸鹽巖裂縫預測中取得成功;姜傳金等[4]基于縱波阻抗、頻率衰減的方位各向異性信息,準確預測了松遼盆地北部徐家圍子斷陷營城組火山巖的裂縫發(fā)育情況。此外,尹志恒等[5]和劉振峰等[6]等對國內外縱波方位各向異性裂縫預測的科研實踐做了詳細的調研?？傮w來看,縱波方位各向異性裂縫預測技術在不同巖性地層的裂縫預測研究中得到越來越廣泛的應用[7-8]。

前人應用研究成果表明,利用縱波方位各向異性預測裂縫的發(fā)育特征是行之有效的。但是,在選取用于各向異性橢圓擬合的地震屬性時,研究人員往往采用兩種思路:一是經(jīng)驗性的地震屬性優(yōu)選;二是計算出若干種地震屬性數(shù)據(jù)體,分別進行各向異性橢圓擬合并求取出裂縫分布特征,再通過對比各種計算結果與測井裂縫信息的吻合度來確定合適的地震屬性。這兩種方法中,第一種方法帶有人為主觀因素,第二種方法則會極大地增加科研工作量。那么,如何合理、簡便、有效地進行地震屬性優(yōu)選,成為影響方法應用效果的實際問題。

在地震儲層預測研究中,人們經(jīng)常通過正演模擬的方式來確定對儲層最為敏感的地震屬性。相似地,針對裂縫的各向異性正演應當也同樣可以進行裂縫敏感屬性的優(yōu)選。但是,目前將各向異性正演與方位各向異性裂縫預測方法相結合的已發(fā)表文獻尚不多見。

我們將各向異性正演引入方位各向異性裂縫預測中,解決該方法在實際應用中的地震屬性優(yōu)選問題。首先,結合測井、試油信息進行基于HTI介質的正演模擬,得到含有方位各向異性信息的井旁地震道集;然后,對正演地震道集進行多種地震屬性的計算和裂縫敏感性分析;最后,將優(yōu)選出的最佳敏感屬性用于各向異性橢圓擬合和研究區(qū)儲層的裂縫預測。實際數(shù)據(jù)的應用試驗結果表明,基于HTI介質的各向異性正演能夠有效指導用于各向異性橢圓擬合的地震屬性優(yōu)選,優(yōu)化方位各向異性裂縫預測方法的應用過程。

1 原理及思路

1.1 方位各向異性裂縫預測原理

地層重力壓實作用的存在,使高角度裂縫比低角度裂縫更容易保存下來,同時,高角度裂縫對裂縫型油氣藏的貢獻也遠遠大于低角度裂縫。因此,與表征高角度裂縫的HTI介質的相關研究就變得意義重大。Ruger[9-11]在HTI介質Thomsen系數(shù)的基礎上,得到圖1所示模型中縱波速度以及反射系數(shù)隨方位角和入射角變化的公式:

(1)

圖1 Ruger[9-11]的公式(1)對應的地質模型

當入射角較小時,(1)式中的最后一項可以略去。略去最后一項并令

(2)

則(1)式變?yōu)?/p>

(3)

式中:Biso為各向同性項;Bani為各向異性項。

當?shù)卣饠?shù)據(jù)包含3個或3個以上不同入射角、方位角振幅信息時,結合子波可以得到對應的不同入射角、方位角的反射系數(shù)。此時,通過求解超定方程組可以得到Biso和Bani,然后利用Biso和Bani這兩項實現(xiàn)方位各向異性橢圓的擬合。Mallick等[12]認為,由各向同性項和各向異性項擬合出來的橢圓可以用來表征裂縫的發(fā)育特征,裂縫的發(fā)育程度由擬合橢圓的長、短軸之比來量化,裂縫的方向則由橢圓的長軸或短軸予以指示。

1.2 各向異性正演原理

裂縫的地震特性受裂縫的傾角、開度、延伸和充填物等眾多因素影響。因此,當巖石含裂縫后其物理特性將變得更為復雜,選擇合適的等效介質模型是進行各向異性正演的基礎。目前,關于裂縫的等效介質模型有Hudson模型[13]、Thomsen模型[14]以及Eshelby-Cheng模型[15]。本文各向異性正演選用的是Hudson模型,通過對含有薄硬幣狀的橢球縫隙或包含物的彈性固體中的平均波長的散射理論分析得到。其模型基于如下假設[16]:①從形狀上,理想化的裂隙形狀為硬幣狀,即高寬比要小,裂隙彼此之間是隔離的,隙間不存在流體流動;②介質包含裂隙半徑及縫隙間距遠小于地震波長的定向疏排列裂隙;③裂縫包體內所含氣體、液體或其它物質的體積模量和剪切模量比圍體的小。

其等效的剛度系數(shù)為

(4)

(5)

式中:V為介質基質的體積;N為體積V內的裂縫總數(shù),即體密度;a為裂縫半徑;φ為裂縫孔隙度;α為裂縫扁率(縱橫比)。

1.3 各向異性正演的思路與過程

由于各向異性正演提供了裂縫在井旁地震道的地震響應,包括疊前各方位角的AVO特征及在裂縫影響下的AVO特征隨方位角的變化規(guī)律,因此,它可以為縱波方位各向異性裂縫分布特征預測提供基礎信息。由正演模擬得到的地震合成記錄獲得振幅、頻率和相位等疊前地震動力學屬性,通過計算這些動力學屬性對于裂縫的敏感性,可以簡便、有效地實現(xiàn)各向異性橢圓擬合的地震屬性優(yōu)選。圖2給出了本文各向異性正演的研究思路。

圖2 各向異性正演研究思路

基于HTI介質Hudson模型各向異性正演過程為:假設在井口布置如圖3所示的觀測系統(tǒng),通過測井資料得到巖石的縱波速度、橫波速度和密度,在給出裂縫所含流體性質及對應彈性模量的情況下,確定出Hudson等效介質模型的剛度系數(shù)及對應的各向異性系數(shù);再根據(jù)Ruger的公式(1)計算出與入射角、方位角有關的反射系數(shù);最后,通過與測井標定子波的褶積,得到不同入射角、方位角的地震波振幅信息。

圖3 各向異性正演模擬觀測系統(tǒng)圖示

2 實際應用試驗

2.1 研究區(qū)地質背景

研究區(qū)位于松遼盆地南部北端某地區(qū),目的層為白堊系營城組火山巖,地層厚度變化較大(鉆井顯示目的層厚度為350～900m)。區(qū)內最早的探井為A井,于2007年開鉆,并于同年在營城組目的層火山巖進行試油,使用50.8mm油嘴獲產(chǎn)天然氣592m3/d,水62m3/d。一年后,在A井附近部署了另外一口探井B,在相近目的層使用12mm油嘴獲天然氣產(chǎn)量20.2464×104m3/d,展現(xiàn)了該地區(qū)良好的勘探潛力。

根據(jù)鉆井資料分析,研究區(qū)火山巖儲層由孔、洞、縫組成,裂縫類型主要包括構造縫和成巖縫,裂縫是保持該地區(qū)儲層連通性和流體運移的主要通道。如何利用地震數(shù)據(jù)準確有效地預測該地區(qū)火山巖裂縫的分布特征是亟待解決的問題。

2.2 正演井段選取

為了檢驗用各向異性正演分析裂縫方位各向異性方法的有效性,需要選取合理的正演參數(shù)來分析不同流體和不同裂縫密度情況下的地層方位各向異性情況,并計算多種地震屬性,分析各種地震屬性的裂縫敏感性,從而指導后續(xù)的疊前方位各向異性裂縫預測。

研究區(qū)內有4口鉆井,分別為A,B,C,D井,其中A井和B井資料較為完全,且有FMI(地層微電阻率成像)測井資料。根據(jù)對FMI裂縫測井資料和試油結果的綜合分析,選取A井的3個井段(自上而下依次命名為A1,A2,A3,其FMI井壁圖像見圖4)進行各向異性正演,具體參數(shù)見表1。由表1 可以看出,所選井段的FMI裂縫密度、試氣結果均有一定差異,可以進行正演模擬結果的對比分析。

圖4 各向異性正演所選取A1，A2和A3井段的FMI井壁圖像

表1 各向異性正演選取井段參數(shù)及試氣結果

2.3 多屬性敏感性分析優(yōu)選

在儲層預測研究中,疊后地震屬性是一種較常規(guī)的預測儲層及流體分布特征的技術手段。在火山巖地層中,流體及裂縫易引起地震數(shù)據(jù)的頻譜發(fā)生變化,本次研究選取的屬性為與地震頻譜相關的一類屬性,主要包括總能量、最大能量、65%衰減頻率、85%衰減頻率、35Hz能量比、起始衰減頻率和相對波阻抗等7種地震屬性。

研究中,首先利用正演出的道集在給定時窗內分別計算上述7種地震屬性,并將計算出的數(shù)據(jù)導出為ASCII格式;然后以方位角0,入射角0時的屬性值為參考基準值,針對所選取的A1,A2和A3 3個井段,結合矩陣實驗室(Matlab)三維成圖,分析各地震屬性在入射角0≤θ≤30°,方位角0≤φ≤90°對裂縫的敏感性,如圖5至圖8所示。通過對A1井段(無裂縫無流體)、A2井段(裂縫發(fā)育、高產(chǎn)氣層)和A3井段(裂縫發(fā)育、氣水同層)各地震屬性的敏感性分析和結果對比,尋找最適合進行各向異性橢圓擬合的地震屬性。

綜合對比圖5至圖8中7種屬性的裂縫敏感性分析結果,得出以下結論。

圖5 正演模擬道集所計算能量類地震屬性的裂縫敏感性分析a 總能量屬性; b 最大能量屬性

圖6 正演模擬道集所計算頻率類地震屬性的裂縫敏感性分析a 65%衰減對應頻率屬性; b 85%衰減對應頻率屬性; c 起始衰減頻率屬性

圖7 正演模擬道集所計算相對波阻抗屬性的裂縫敏感性分析

圖8 正演模擬道集所計算35Hz低頻能量比屬性的裂縫敏感性分析

1) 能量類地震屬性(圖5)在A1,A2井段表現(xiàn)為明顯的負異常,但在A3井段卻表現(xiàn)為正異常。由于A1井段無裂縫,而A2井段裂縫發(fā)育,那么利用該類地震屬性無法區(qū)分裂縫是否發(fā)育。因此,能量類地震屬性不適用于本區(qū)各向異性橢圓擬合的裂縫預測。

2) 在頻率類屬性(圖6)中,百分比衰減對應頻率屬性(圖6a和圖6b)在A1井段為幅值不大的正異常,在A2井段為幅值不大的負異常,在A3井段為幅值略大的負異常。該類屬性在A1,A2和A3井段的值有一定差異,但幅值變化不夠明顯,在地震數(shù)據(jù)信噪比不高時,無法保證計算結果的準確性。而起始衰減頻率屬性(圖6c)在A1,A2井段表現(xiàn)為正異常,在A3井段為負異常,也不適用于區(qū)分裂縫是否發(fā)育。

3) 相對波阻抗屬性(圖7)在A2井段表現(xiàn)為較高的正異常,在A1井段為中等正異常,在A3井段為較低正異常。該屬性在不同裂縫發(fā)育程度、不同流體性質的地層存在一定區(qū)分性。

4) 低頻能量比屬性(圖8)在裂縫含氣井段A2表現(xiàn)為較高的正異常,在裂縫氣水同層井段A3表現(xiàn)為中等正異常,在無裂縫干層井段A1表現(xiàn)為低的正異常甚至負異常。該屬性對不同井段展現(xiàn)出的這種區(qū)分程度非常有利于進行研究區(qū)各向異性橢圓的擬合。

根據(jù)以上分析結果可知,低頻能量比屬性是研究區(qū)裂縫預測的最敏感屬性。優(yōu)選低頻能量比屬性進行研究區(qū)各向異性橢圓的擬合,預測裂縫的平面分布特征。

2.4 利用優(yōu)選屬性預測裂縫分布特征

利用正演模擬不同方位角道集計算的低頻能量比屬性進行各向異性橢圓的擬合,根據(jù)擬合橢圓的長、短軸之比得到了研究區(qū)營城組裂縫的空間分布預測結果?；诰狭芽p信息,分別從剖面和平面上探討本文方法預測的研究區(qū)裂縫分布特征。

圖9給出了低頻能量比屬性方位各向異性預測的過A井和B井目的層段裂縫密度剖面(圖中色標為各向異性橢圓的長、短軸之比,即各向異性強度值;無量綱),可以看出,在A井試氣段和B井產(chǎn)氣段預測的裂縫密度均為高值。可見利用低頻能量比屬性各向異性預測的裂縫發(fā)育情況不僅與FMI測井資料吻合較好,也符合井上已知的地質認識。

研究區(qū)內發(fā)育若干條南北走向的斷層。從預測的裂縫平面分布特征(圖10)來看,研究區(qū)營二段、營三段和營四段(T4a至T4c)裂縫方向以近南北向為主,北東—南西向為輔(圖10a),營一段(T4c至T41)裂縫方向以南北向為主,東西向為輔(圖10b),裂縫方向與斷層走向較為相似(玫瑰圖中色標代表某一方向的裂縫占該位置目的層段所有裂縫方向的百分比，紅色指示預測出的裂縫優(yōu)勢方向)。預測的裂縫發(fā)育區(qū)主要分布在斷層附近,同時在斷層周邊存在若干成巖縫發(fā)育帶。對于有FMI資料的A井和B井,在T4a-T4c層段,A井和B井裂縫均發(fā)育;在T4c-T41層段,A井裂縫發(fā)育而B井裂縫不發(fā)育。從裂縫分布特征(圖10)來看,本文方法預測結果與單井的裂縫發(fā)育情況較為一致。

圖9 低頻能量比屬性方位各向異性預測的過A井和B井目的層段裂縫密度剖面

圖10 低頻能量比屬性方位各向異性預測的不同層段裂縫平面分布a T4a-T4c層段; b T4c-T41層段

將預測結果與FMI解釋結果進行對比,進一步驗證本文方法的有效性。圖11給出了FMI裂縫方向(圖11a，圖11b，圖11c)與低頻能量比屬性方位各向異性預測的A井裂縫方向(圖11d)的對比結果(白底玫瑰圖為FMI裂縫方向,黑底玫瑰圖為預測方向)。從圖11中可以看出,預測的A井裂縫方向幾乎為FMI測井資料顯示的多種裂縫類型裂縫方向的綜合顯示(B井的對比結果也相似),說明本文方法預測的裂縫方向較為可信。

圖11 低頻能量比屬性方位各向異性預測的A井裂縫方向與FMI對比a 高導縫的走向為北東—南西向; b 高阻縫的走向為南北向; c 微裂縫的走向為北東—南西向; d 本文方法預測的A井裂縫方向

3 結論與認識

各向異性橢圓擬合是縱波方位各向異性裂縫預測方法中重要的一步,其關鍵在于如何合理有效地進行地震屬性的優(yōu)選。我們提出了利用各向異性正演進行方位各向異性橢圓擬合地震屬性優(yōu)選的新思路。通過基于HTI介質Hudson模型的正演模擬得到含有各向異性信息的井旁地震道集,并對地震道集進行多種地震屬性的計算,對比分析各地震屬性的裂縫敏感性,從而優(yōu)選出最適合進行方位各向異性橢圓擬合的地震屬性。松遼盆地南部某地區(qū)白堊系營城組火山巖儲層的實際應用試驗取得了較好的裂縫預測效果,并得到以下認識。

1) 能量類地震屬性、頻率類屬性和相對波阻抗屬性不適用于研究區(qū)各向異性橢圓擬合的裂縫預測,低頻能量比是進行該地區(qū)各向異性橢圓擬合的優(yōu)選敏感屬性。但這一結果并不具有普遍性,在不同地區(qū)對不同儲層進行各向異性正演時,地震屬性的優(yōu)選結果也會不同。

2) 研究區(qū)裂縫分布特征與斷層展布較為一致,營二段、營三段和營四段裂縫方向為以近南北向為主,北東—南西向為輔;營一段裂縫方向為南北向為主,東西向為輔。此外,預測的裂縫分布特征在單井FMI資料上得到了較好的驗證。

3) 多屬性裂縫敏感性分析是基于Hudson模型的高角度裂縫,但以低角度裂縫、水平裂縫為主的儲層是客觀存在的,在面對這一類儲層的裂縫預測問題時,文中涉及的多屬性裂縫敏感性分析技術不適用。

致謝:在本文編寫過程中得到中國地質大學(北京)王彥春教授、劉學清博士和劉志偉博士的指導和啟發(fā),在此表示衷心感謝！

參 考 文 獻

[1] 曲壽利,季玉新,王鑫,等.全方位P波屬性裂縫檢測方法研究[J].石油地球物理勘探,2001,36(4):390-397

Qu S L,Ji Y X,Wang X,et al.Seismic method fro using full-azimuth P-wave attributes to detect fracture [J].Oil Geophysical Prospecting,2001,36(4):390-397

[2] 樂紹東.AVA裂縫檢測技術在川西JM構造的應用[J].天然氣工業(yè),2004,24(4):22-24

Yue S D.Application of AVA fracture detection technique in JM structure in west Sichuan [J].Natural gas industry,2004,24(4):22-24

[3] 喻岳鈺,楊長春,王彥飛,等.瞬時頻域衰減屬性及其在碳酸鹽巖裂縫檢測中的應用[J].地球物理學進展,2009,24(5):1717-1722

Yu Y Y,Yang C C,Wang Y F,et al.P-wave azimuthal attenuation attributes in wavelet-scale domain and its application to fracture detection in carbonate [J].Progress in Geophysics,2009,24(5):1717-1722

[4] 姜傳金,鞠林波,張廣穎,等.利用地震疊前數(shù)據(jù)預測火山巖裂縫的方法和效果分析——以松遼盆地北部徐家圍子斷陷營城組火山巖為例[J].地球物理學報,2011,54(2):515-523

Jiang C J,Ju L B,Zhang G Y,et al.The method and effect analysis of volcanic fracture prediction with prestack seismic data——an example from the volcanic rocks of Yingcheng formation in Xujiaweizi fault depression,north of Songliao Basin [J].Chinese Journal of Geophysics,2011,54(2):515-523

[5] 尹志恒,狄?guī)妥?李向陽,等.國外應用縱波各向異性技術檢測裂縫的研究進展[J].科技導報,2011,29(30):73-78

Yin Z H,Di B R,Li X Y,et al.Progress in P-wave anisotropy technology for fracture detection [J].Science and Technology Review,2011,29(30):73-78

[6] 劉振峰,曲壽利,孫建國.地震裂縫預測技術研究進展[J].石油物探,2012,51(2):191-198

Liu Z F,Qu S L,Sun J G.Progress of seismic fracture characterization technology[J].Geophysical Prospecting for Petroleum,2012,51(2):191-198

[7] 孫煒,王彥春,李梅,等.利用疊前地震數(shù)據(jù)預測火山巖儲層裂縫[J].物探與化探,2010,34(2):229-232

Sun W,Wang Y C,Li M,et al.The detection of fractures in volcanic reservoir with pre-stack seismic data [J].Geophysical and Geochemical Exploration,2010,34(2):229-232

[8] 孫煒,李玉鳳,何巍巍,等.P波各向異性裂縫預測技術在ZY區(qū)碳酸鹽儲層中的應用[J].石油與天然氣地質,2013,34(2):137-144

Sun W,Li Y F,He W W,et al.Sun Using P-wave azimutal anisotropy to predict fractures in carbonate reservoirs of the ZY block[J].Oil & Gas Geology,2013,34(2):137-144

[9] Ruger A.Reflection coefficients and azimuthal AVO analysis in anisotropic media[D].Colorado:Colorado School of Mines,1996

[10] Ruger A.P-wave reflection coefficients for transversely isotropic models with vertical and horizontal axis of symmetry [J].Geophysics,1997,62(3):713-722

[11] Ruger A.Using AVO for fracture detection:analytic basis and practical solutions [J].The Leading Edge,1997,16(10):1429-1430,1432-1434

[12] Mallick S,Craft K L,Meister L,et al.Determination of the principle direction of azimuthal anisotropy from P-wave seismic data [J].Geophysics,1998,63(2):692-706

[13] Hudson J A.Wave speeds and attenuation of elastic waves in material containing cracks[J].Geophysical Journal of the Royal Astronomical Society,1981,64:133-150

[14] Thomsen L.Elastic anisotropy due to aligned cracks in porous rock[J].Geophysical Prospecting,1995,43(6):805-829

[15] Cheng C H.Crack models for a transversely anisotropic medium[J].Journal of Geophysical Research,1993,98:675-684

[16] 葛瑞·馬沃可,塔潘·木克基,杰克·徳沃金.巖石物理手冊:孔隙介質中地震分析工具[M].徐海濱,戴建春，譯.合肥:中國科學技術大學出版社,2008:113-118

Mavko G,Mukerji T,Dvorikin J.The rock physics handbook:tools for seismic analysis in porous media[M].Xu H B, Dai J C,translators.Hefei:University of Science and Technology of China Press,2008:113-118

[17] Crampin S.Effective anisotropic elastic constants for wave propagation through cracked solids[J].Geophysical Journal of the Royal Astronomical Society,1984,76(1):135-145