姜曉宇 張 研 甘利燈 宋 濤 杜文輝 周曉越

(中國(guó)石油勘探開發(fā)研究院，北京 100083)

0 花崗巖潛山油氣藏勘探現(xiàn)狀

目前世界上已經(jīng)發(fā)現(xiàn)的基巖油氣藏中，花崗巖油氣藏約占40%，儲(chǔ)量占75%[1]。在花崗巖潛山中發(fā)現(xiàn)油氣藏的有：利比亞錫爾特盆地的Nafoora油田；委內(nèi)瑞拉馬拉開波盆地的La Paz油田；智利西麥哲倫盆地的Lago_Mercedes油田；美國(guó)堪薩斯中部的Orth油田；越南湄公河三角洲的白虎油田；乍得Bongor盆地的Ronier油田等[2-4]。

另外，還有很多基底為花崗巖的含油氣盆地，如：中東波斯灣油氣區(qū)；俄羅斯、哈薩克斯坦的濱里海盆地；蘇丹的穆格萊特盆地；俄羅斯的伏爾加—烏拉爾盆地；澳大利亞的吉普斯蘭盆地；加拿大西部盆地；美國(guó)的洛杉磯盆地、阿巴拉契亞盆地；委內(nèi)瑞拉的東委內(nèi)瑞拉盆地等[5-6]。

中國(guó)花崗巖油氣藏主要分布在渤海灣盆地，如：遼東灣遼西潛山帶的錦州25-1S油氣田；渤海廟西凸起至渤東凹陷東的蓬萊9-1油田；濟(jì)陽(yáng)坳陷陳家莊凸起西南部的王莊油田；渤海沙壘田凸起西北傾沒端的曹妃甸1-6油田等。此外，酒西盆地的鴨兒峽油田和珠江口凹陷的惠州凹陷、松遼盆地[7-9]、柴達(dá)木盆地的局部發(fā)育花崗巖基底[5]。

受地震資料品質(zhì)和地質(zhì)認(rèn)識(shí)的限制，一般情況下，花崗巖潛山不是盆地勘探初期的主要目標(biāo)，潛山油氣藏往往是淺層碎屑巖油氣藏加深兼探而發(fā)現(xiàn)的。以乍得Bongor盆地花崗巖潛山油藏的發(fā)現(xiàn)為例，2007年發(fā)現(xiàn)下白堊統(tǒng)砂巖油氣田，隨著地震資料采集、處理技術(shù)的提高和地質(zhì)認(rèn)識(shí)的深入，系統(tǒng)復(fù)查已鉆遇的測(cè)井和巖心資料以分析基巖潛山類型，結(jié)合鉆井過程中的泥漿漏失情況，確定基巖存在以裂縫為主的儲(chǔ)集空間，可成為油氣成藏的有利場(chǎng)所。2012年針對(duì)潛山進(jìn)行鉆探，獲得了油氣發(fā)現(xiàn)。

隨著花崗巖潛山油氣藏的持續(xù)開發(fā)，潛山儲(chǔ)層出現(xiàn)非均質(zhì)性強(qiáng)、單井產(chǎn)能差異大、部分油井之間存在井間干擾、含水上升快等問題。裂縫識(shí)別與預(yù)測(cè)是解決上述問題的關(guān)鍵，也一直是具有挑戰(zhàn)性的世界性難題[10-17]。目前國(guó)內(nèi)采用的裂縫預(yù)測(cè)方法主要有：基于構(gòu)造應(yīng)力場(chǎng)或應(yīng)變場(chǎng)數(shù)值模擬；疊后地震資料相干和曲率等屬性分析。二者均存在尺度大、精度低等問題，而基于各向異性的疊前地震資料裂縫預(yù)測(cè)技術(shù)尚處于探索階段。本文分析花崗巖潛山裂縫地震預(yù)測(cè)技術(shù)的前提條件、應(yīng)用效果及優(yōu)缺點(diǎn)，以期為花崗巖潛山油氣藏的勘探提供指導(dǎo)作用。

1 花崗巖潛山裂縫地震預(yù)測(cè)技術(shù)

花崗巖潛山裂縫地震預(yù)測(cè)技術(shù)大體上可分為疊前和疊后兩類，目前常用的裂縫地震預(yù)測(cè)技術(shù)的優(yōu)缺點(diǎn)如表1所示，不同技術(shù)適用性及應(yīng)用效果不同。

由于潛山內(nèi)幕地震資料成像差，信噪比低，一般可在裂縫預(yù)測(cè)前采用構(gòu)造導(dǎo)向?yàn)V波處理技術(shù)提高地震資料品質(zhì)。

構(gòu)造導(dǎo)向?yàn)V波技術(shù)是利用地層傾角和方位角沿地層定向?yàn)V波，再利用曲率和相干屬性描述地層不連續(xù)性，并對(duì)無意義的不連續(xù)性做平滑處理，從而達(dá)到邊緣保護(hù)的目的，使地震數(shù)據(jù)同相軸的連續(xù)性和錯(cuò)斷特征更明顯。地震資料經(jīng)構(gòu)造導(dǎo)向?yàn)V波技術(shù)處理后，同相軸能量增強(qiáng)，斷點(diǎn)清晰，潛山內(nèi)幕信息更加明確(圖1)。

表1 花崗巖潛山裂縫地震預(yù)測(cè)技術(shù)對(duì)比

圖1 構(gòu)造導(dǎo)向?yàn)V波處理前(a)、后(b)效果對(duì)比[18]箭頭所指為花崗巖潛山頂面

1.1 裂縫疊后地震預(yù)測(cè)技術(shù)

1.1.1 相干體屬性技術(shù)

相干體屬性技術(shù)是基于地震反射同相軸的不連續(xù)性預(yù)測(cè)斷裂分布。相干類屬性有以下優(yōu)點(diǎn)：

(1)相干體屬性技術(shù)不受任何解釋誤差的影響，可直接從三維地震數(shù)據(jù)體中獲取斷層和地層信息，用于斷層、巖性體邊界的不連續(xù)性檢測(cè)，極大地提高了解釋精度；

(2)對(duì)橫向一致的地層、構(gòu)造進(jìn)行了壓縮，提取的水平時(shí)間切片可以顯示任意方向的斷層，解決了平行同相軸的斷層難以解釋的問題；

(3)由于裂縫性儲(chǔ)層的非均質(zhì)性會(huì)引起地震反射特征的不規(guī)則變化，在相干體切片上表現(xiàn)為低相干的特征，因此相干體屬性技術(shù)也可應(yīng)用于裂縫性儲(chǔ)層發(fā)育區(qū)的預(yù)測(cè)[19]。

相干體屬性技術(shù)大致經(jīng)歷了三代發(fā)展。第一代是基于歸一化互相關(guān)算法的相干體屬性技術(shù)(Correlation)，計(jì)算速度快，分辨率高，但只適用于信噪比高的地震數(shù)據(jù)。第二代是基于相似性算法的相干體屬性技術(shù)(Semblance)，可以對(duì)任意多道地震數(shù)據(jù)計(jì)算相似性，抗噪性強(qiáng)，但分辨率低。第三代是基于本征值算法的相干體屬性技術(shù)(Eigenstructure)，將多道地震數(shù)據(jù)組成協(xié)方差矩陣，并計(jì)算特征值，提高了抗噪能力和分辨率，但計(jì)算量大，而且計(jì)算中采用了以計(jì)算點(diǎn)為中心的固定時(shí)窗長(zhǎng)度，無法兼顧具有寬頻帶的淺層數(shù)據(jù)和窄頻帶的深層數(shù)據(jù)，導(dǎo)致計(jì)算的相干體屬性不能真實(shí)反映波組特征[20-21]。

近年來，基于幾何結(jié)構(gòu)張量的相干體屬性技術(shù)、基于高階統(tǒng)計(jì)量的相干體屬性技術(shù)、基于小波變換的多尺度相干體屬性計(jì)算技術(shù)等對(duì)第三代相干體屬性技術(shù)進(jìn)行了改進(jìn)，也有人稱為新一代(或第四代)相干體屬性技術(shù)。但是，相干體屬性技術(shù)的應(yīng)用仍有一定的局限性：首先，要求地震數(shù)據(jù)的信噪比高；其次，相干體屬性技術(shù)對(duì)地層的細(xì)微彎曲不敏感，相干體屬性無法識(shí)別反射波同相軸不斷開的斷裂(圖2a)。對(duì)于斷裂構(gòu)造復(fù)雜的花崗巖潛山裂縫型儲(chǔ)層，廣泛發(fā)育著不同級(jí)次的斷層和裂縫，地層產(chǎn)狀(傾角和方位角)變化大，尤其是潛山內(nèi)幕地震資料存在照明不均、成像差、信噪比低等問題，即使采用最先進(jìn)的相干體算法，結(jié)果也不盡理想。

1.1.2 曲率屬性技術(shù)

曲率屬性屬于地震幾何屬性的一種，與地震反射體的彎曲程度相對(duì)應(yīng)，用于識(shí)別裂縫、斷層和褶皺等地質(zhì)構(gòu)造。一般曲率越大，裂縫就越發(fā)育。

曲率屬性技術(shù)是基于垂直作用下的彎曲薄板理論，應(yīng)用時(shí)須滿足以下前提條件: ①巖石是脆性的，裂縫形成主要是由于巖石破裂作用，不考慮巖層的塑性變形作用; ②只能用于預(yù)測(cè)彎曲巖層面上抗張應(yīng)力形成的張性縫，即斷裂兩側(cè)若均不存在地層的褶曲現(xiàn)象，則曲率屬性無法識(shí)別(圖2b)； ③巖層受力變形而彎曲[22-23]。

對(duì)于相干體屬性無法預(yù)測(cè)的反射波同相軸不斷開現(xiàn)象(圖2a)，由于地層發(fā)生了彎曲，曲率屬性可以識(shí)別。地層既發(fā)生錯(cuò)斷又產(chǎn)生彎曲現(xiàn)象，相干和曲率屬性都可以識(shí)別 (圖2c)。乍得Bongor盆地某區(qū)花崗巖潛山的曲率屬性切片(圖3)和基于曲率體計(jì)算的裂縫玫瑰花圖(圖4)均可以清楚地展現(xiàn)潛山的裂縫發(fā)育方向與程度。

圖2 曲率屬性、相干屬性識(shí)別斷裂比較(a)相干屬性不能識(shí)別，曲率屬性能識(shí)別； (b)曲率屬性不能識(shí)別，相干屬性可識(shí)別； (c)曲率、相干屬性都可識(shí)別。黑線表示斷層，其他色線表示層位

圖3 某區(qū)曲率屬性平面分布圖深色表示裂縫發(fā)育，藍(lán)色虛線為北東—南西走向斷裂[24]

圖4 某區(qū)裂縫走向玫瑰花屬性水平切片[24]黑線表示裂縫方向

1.1.3 譜分解技術(shù)

譜分解技術(shù)是通過提高地震資料分辨率增強(qiáng)儲(chǔ)層的橫向不連續(xù)性，從而描述微斷裂信息[25]。Stockwell等[26]提出的S變換結(jié)合了小波變換和短時(shí)傅里葉變換的優(yōu)點(diǎn)，在信號(hào)分析中得到了廣泛的應(yīng)用。自該方法提出以來，為了突破S變換中固定的基本小波的限制，許多學(xué)者對(duì)S變換算法提出了改進(jìn)方案[27]，并通過改變S變換中窗函數(shù)隨頻率變化的趨勢(shì)，提出了廣義S變換的概念。應(yīng)用廣義S變換對(duì)地震資料進(jìn)行頻譜分解可實(shí)現(xiàn)分頻顯示，并可獲取豐富的地質(zhì)信息，對(duì)于識(shí)別花崗巖潛山發(fā)育的小斷層具有良好的效果[28-34]。

1.1.4 螞蟻體追蹤技術(shù)

螞蟻體追蹤技術(shù)基于蟻群算法實(shí)現(xiàn)對(duì)斷裂的追蹤和識(shí)別。該算法原理為模擬螞蟻在食物與巢穴之間根據(jù)可吸引螞蟻的信息素濃度尋求最短路徑。在地震數(shù)據(jù)中，“螞蟻”根據(jù)振幅及相位之間的差異，沿著可能的斷層和裂縫移動(dòng)完成對(duì)二者的刻畫。

與相干屬性相比(圖5)，螞蟻體屬性凸顯了斷裂線狀構(gòu)造特征，去除了與斷裂無關(guān)的信息，提高了斷裂解釋精度。缺點(diǎn)是平面預(yù)測(cè)結(jié)果往往過于雜亂，無規(guī)律。原因之一是控制螞蟻?zhàn)粉櫧Y(jié)果的參數(shù)太多，調(diào)節(jié)困難。

圖5 相干屬性(a)與螞蟻體屬性(b)平面切片對(duì)比深色表示斷裂或裂縫位置

螞蟻體追蹤技術(shù)是基于疊后地震數(shù)據(jù)運(yùn)算的，雖然其精度比相干等屬性高，但也只適用于對(duì)小斷層和大尺度裂縫的預(yù)測(cè)?？深A(yù)測(cè)裂縫發(fā)育的方向，但難以定量化表征裂縫發(fā)育密度。

王軍等[35]、馬曉宇等[36]將螞蟻體追蹤技術(shù)與疊前各向異性技術(shù)相結(jié)合，將疊前分方位各向異性信息融入到螞蟻體追蹤技術(shù)中，得到裂縫的方向與密度，使裂縫預(yù)測(cè)由定性轉(zhuǎn)向定量化，可預(yù)測(cè)大尺度及中等尺度裂縫發(fā)育帶。圖6為與基于螞蟻體各向異性與基于振幅各向異性預(yù)測(cè)的裂縫密度切片對(duì)比，圖中顏色的變化表示裂縫密度的變化，紅色區(qū)域?yàn)榱芽p集中發(fā)育區(qū)域?？梢钥闯?，基于螞蟻體各向異性方法不僅可以清晰地預(yù)測(cè)裂縫發(fā)育帶展布特征，還可以定量計(jì)算裂縫密度，對(duì)于交叉裂縫也能有效識(shí)別，比基于振幅各向異性刻畫的裂縫更為精細(xì)和全面。

1.1.5 最大似然屬性技術(shù)

Hale[37]在研究斷面提取和斷距估算時(shí)提出了最大似然屬性。最大似然屬性是通過對(duì)整個(gè)地震數(shù)據(jù)體掃描、計(jì)算數(shù)據(jù)樣點(diǎn)之間的相似性，獲得斷裂發(fā)育的最可能位置及概率。結(jié)合傾角掃描和構(gòu)造導(dǎo)向?yàn)V波等技術(shù)，可以增強(qiáng)斷裂的地震成像效果，提升斷裂刻畫精度。

最大似然屬性主要包括Likelihood和Thinlikelihood屬性計(jì)算兩個(gè)方面。Likelihood屬性的計(jì)算基于地震相似性，是為了壓制噪聲從而突出斷裂的成像。與相干體技術(shù)中的Semblance屬性不同，Likelihood屬性是Semblance屬性8次冪與1的差值，因此放大了相鄰樣點(diǎn)間相似性的對(duì)比關(guān)系，進(jìn)一步凸顯斷裂的成像。Thinlikelihood屬性是在Likelihood屬性體計(jì)算的基礎(chǔ)上進(jìn)行樣點(diǎn)掃描，對(duì)于每個(gè)樣點(diǎn)附近一定步長(zhǎng)范圍內(nèi)不同走向、不同傾角處Likelihood屬性值只保留最大值及其對(duì)應(yīng)的傾角、方位角信息，其余值設(shè)置為0，以此記錄斷裂最有可能發(fā)育的位置及該位置斷裂發(fā)育的概率。Thinlikelihood屬性比Likelihood對(duì)斷裂的刻畫更加準(zhǔn)確。

圖6 基于振幅各向異性(a)與基于螞蟻體各向異性(b)預(yù)測(cè)的裂縫密度切片對(duì)比[36]

最大似然屬性技術(shù)關(guān)鍵步驟包括： ①分析斷裂的地震反射特征，便于最大似然屬性計(jì)算參數(shù)的提取以及切片的地質(zhì)解釋； ②傾角控制下斷裂成像加強(qiáng)，目的是去除隨機(jī)噪聲、提高信噪比、突出斷裂位置同相軸錯(cuò)斷的成像效果； ③提取最大似然屬性;④解釋屬性切片，在斷裂地質(zhì)模式的指導(dǎo)下進(jìn)行斷裂刻畫[38]。

圖7是花崗巖潛山最大似然體(Thinlikelihood屬性)與相干體屬性平面對(duì)比圖。由圖可見： 相干體屬性比較清楚地刻畫了主干斷裂，但是無法刻畫斷裂帶內(nèi)部結(jié)構(gòu)以及裂縫發(fā)育區(qū)； 最大似然體屬性比相干體精度高，比較清楚地刻畫了斷裂帶內(nèi)部結(jié)構(gòu)，并且對(duì)于裂縫密集發(fā)育區(qū)的預(yù)測(cè)效果較好，但對(duì)主干斷裂整體展布的刻畫效果不如相干屬性。

1.2 裂縫疊前地震預(yù)測(cè)技術(shù)

疊后方法預(yù)測(cè)裂縫速度快，但精度低，難以實(shí)現(xiàn)裂縫密度與方位信息的定量預(yù)測(cè)。近年來，關(guān)于裂縫疊前地震預(yù)測(cè)技術(shù)的研究越來越多。地震波在裂縫型介質(zhì)中傳播時(shí)表現(xiàn)為旅行時(shí)、速度、振幅、衰減、頻率、相位等方位各向異性。雖然橫波比縱波的方位各向異性更強(qiáng)，但由于橫波的勘探成本高且信噪比低，限制了橫波裂縫預(yù)測(cè)技術(shù)的發(fā)展。多分量轉(zhuǎn)換橫波裂縫檢測(cè)技術(shù)尚不完善。VSP地震探測(cè)范圍有限而無法對(duì)裂縫進(jìn)行區(qū)域檢測(cè)。目前，常用的裂縫疊前地震預(yù)測(cè)技術(shù)還是基于縱波振幅方位各向異性[39-47]。

基于縱波振幅方位各向異性裂縫預(yù)測(cè)技術(shù)是利用縱波振幅隨方位角變化規(guī)律預(yù)測(cè)裂縫走向和密度。郝守玲等[48]對(duì)裂縫介質(zhì)的縱波方位各向異性特征進(jìn)行物理模型試驗(yàn)時(shí)發(fā)現(xiàn)，當(dāng)裂縫走向與觀測(cè)方位平行時(shí)(夾角為0°)，反射波振幅和速度最大；隨著裂縫走向與測(cè)線方位之間夾角的增大，反射波的振幅和速度逐漸減小；當(dāng)夾角為90°時(shí)，反射波振幅和速度最小。Vladimir等[49]將振幅或旅行時(shí)等隨著裂縫走向與測(cè)線方位之間夾角的變化規(guī)律總結(jié)為橢圓方程

R(φ)=A+Bsin2φ

(1)

圖7 過花崗巖潛山屬性平面圖(a)相干體； (b)最大似然體。深色表示斷裂或裂縫發(fā)育位置

式中:R為振幅、旅行時(shí)或速度的幅值；φ為觀測(cè)方位與裂縫的夾角，當(dāng)φ=0°時(shí)，表示裂縫走向，當(dāng)φ=90°時(shí)，表示垂直裂縫走向；A表示各向同性振幅、旅行時(shí)或速度；B表示振幅、旅行時(shí)或速度隨反射角的變化量。

基于縱波屬性方位各向異性的裂縫疊前地震預(yù)測(cè)技術(shù)對(duì)地震資料品質(zhì)要求較高。傳統(tǒng)窄方位采集的地震資料已不能滿足花崗巖潛山裂縫預(yù)測(cè)精度要求。隨著資料采集、處理技術(shù)的發(fā)展，“兩寬一高”(寬頻、寬方位、高密度)地震勘探技術(shù)逐漸興起。與常規(guī)地震資料相比，高密度、寬頻、寬方位地震資料能顯著改善潛山的成像質(zhì)量、拓寬頻帶和提高信噪比，可提高潛山頂面及內(nèi)幕成像精度(圖8)。

利用“兩寬一高”地震資料提供的OVT域五維道集資料能提高各向異性分析的精度[50]。夏亞良等[51]以中非花崗巖潛山裂縫儲(chǔ)層為例，以花崗巖裂縫地質(zhì)特征為指導(dǎo)，以式(1)為理論基礎(chǔ)，利用橢圓擬合法預(yù)測(cè)裂縫，不僅能預(yù)測(cè)裂縫發(fā)育的走向(圖9a黑色線)，還能定量預(yù)測(cè)裂縫發(fā)育的強(qiáng)度(圖9a彩色部分)，預(yù)測(cè)結(jié)果與測(cè)井成像結(jié)果一致。但是，用該方法預(yù)測(cè)裂縫時(shí)，需要進(jìn)行局部炮檢距和方位角疊加，可能會(huì)影響預(yù)測(cè)結(jié)果的準(zhǔn)確性。另外，OVT處理后的道集資料數(shù)據(jù)量大，運(yùn)算時(shí)間長(zhǎng)，如果擬合點(diǎn)過多，運(yùn)算時(shí)間將會(huì)更長(zhǎng)。

圖8 常規(guī)地震資料(a)與“兩寬一高”地震資料(b)[50]對(duì)比

圖9 花崗巖基底裂縫強(qiáng)度及走向平面預(yù)測(cè)結(jié)果(a)裂縫預(yù)測(cè)圖黑色桿狀表示走向，彩色表示裂縫強(qiáng)度； (b)測(cè)井裂縫走向玫瑰圖。從左至右分別為W1、W2、W3井[51]

花崗巖潛山可能發(fā)育單組或多組高角度縫和網(wǎng)狀縫?；跈E圓擬合法的疊前縱波方位各向異性裂縫預(yù)測(cè)技術(shù)更適用于預(yù)測(cè)單組高角度裂縫發(fā)育帶，不太適用于多組裂縫發(fā)育帶(多個(gè)方位的裂縫)。多個(gè)方位裂縫發(fā)育時(shí)，反射系數(shù)隨方位變化特征更為復(fù)雜，導(dǎo)致預(yù)測(cè)結(jié)果不太準(zhǔn)確。如圖10所示，在單組裂縫發(fā)育情況下，反射系數(shù)隨方位呈近似橢圓變化(圖10a)； 在正交裂縫情況下，反射系數(shù)隨方位變化規(guī)律復(fù)雜(圖10b)。利用不同炮檢距的振幅進(jìn)行橢圓擬合法預(yù)測(cè)裂縫存在多解性。以圖11為例，當(dāng)發(fā)育有兩組走向不同的裂縫時(shí)，小炮檢、中炮檢距、大炮檢距擬合的裂縫走向均不同。同時(shí)，擬合橢圓的離心率會(huì)隨著裂縫組數(shù)的增多而變小，各向異性強(qiáng)度也隨之減弱，導(dǎo)致部分裂縫無法識(shí)別。風(fēng)化殼和基底下方可能發(fā)育的網(wǎng)狀縫，由于不存在各向異性響應(yīng)，疊前裂縫預(yù)測(cè)技術(shù)不適用。

針對(duì)上述問題，王霞等[52]、陳志剛等[53]提出了基于統(tǒng)計(jì)法的各向異性強(qiáng)度預(yù)測(cè)技術(shù)。與橢圓擬合法不同的是，統(tǒng)計(jì)法各向異性強(qiáng)度預(yù)測(cè)雖然也是基于縱波方位各向異性原理，但沒有進(jìn)行局部炮檢距與方位角疊加，而是首先通過道集規(guī)則化，將不規(guī)則的道集處理成螺旋道集后，直接在OVT道集中統(tǒng)計(jì)振幅或雙程走時(shí)的方差以指示各向異性的強(qiáng)度。計(jì)算方差的離散公式

(2)

圖10 反射系數(shù)隨裂縫方位變化特征[36](a)單組裂縫; (b)正交裂縫

圖11 橢圓擬合法的多解性分析[52]a)振幅切片，黃綠色為強(qiáng)振幅; (b)小炮檢距擬合;(c)中炮檢距擬合; (d)大炮檢距擬合

式中：β為方差；N為離散樣點(diǎn)個(gè)數(shù)；Ri為道集中某個(gè)點(diǎn)的振幅或雙程旅行時(shí)；R′為振幅或雙程旅行時(shí)的平均值。

與橢圓擬合法相比，方位統(tǒng)計(jì)法對(duì)方位角進(jìn)行離散，考慮了每組裂縫導(dǎo)致的方位屬性值的偏離程度。當(dāng)發(fā)育多組裂縫時(shí)，各向異性強(qiáng)度值增大，減少了多組裂縫預(yù)測(cè)的多解性。

圖12所示的橢圓擬合法預(yù)測(cè)結(jié)果符合率雖達(dá)到80%，但從平面圖上看，裂縫在全區(qū)普遍發(fā)育，規(guī)律性不明顯。統(tǒng)計(jì)法各向異性強(qiáng)度預(yù)測(cè)結(jié)果與鉆井符合率高且地質(zhì)規(guī)律明顯(圖13)。

圖12 橢圓擬合法各向異性強(qiáng)度預(yù)測(cè)結(jié)果(a)各向異性強(qiáng)度平面圖； (b)PS-3井(上)、RS-10井(中)和ME-2井(下)成像測(cè)井玫瑰圖(左)與預(yù)測(cè)玫瑰圖(右)對(duì)比[53]

圖13 統(tǒng)計(jì)法各向異性強(qiáng)度預(yù)測(cè)結(jié)果(a)各向異性強(qiáng)度平面圖； (b)PS-3井(上)、RS-10井(中)和ME-2井(下)成像測(cè)井玫瑰圖(左)與預(yù)測(cè)玫瑰圖(右)對(duì)比[53]

3 認(rèn)識(shí)與展望

(1)構(gòu)造導(dǎo)向?yàn)V波處理技術(shù)可改善潛山內(nèi)幕地震資料成像品質(zhì)。

(2)花崗巖裂縫疊后地震預(yù)測(cè)技術(shù)主要有相干、曲率、螞蟻?zhàn)粉?、最大似然屬性等。每種技術(shù)都有其優(yōu)缺點(diǎn)及適用條件。其中螞蟻體追蹤與最大似然屬性比相干等屬性精度高，但也只能實(shí)現(xiàn)對(duì)大尺度及中等尺度裂縫發(fā)育帶的定性預(yù)測(cè)。

(3)裂縫疊前地震預(yù)測(cè)技術(shù)能定量預(yù)測(cè)裂縫發(fā)育的密度和方位，預(yù)測(cè)效果較疊后方法精度更高；橢圓擬合法適用于單組裂縫預(yù)測(cè)，統(tǒng)計(jì)法適用于多組裂縫預(yù)測(cè)；疊后預(yù)測(cè)方法更適用于風(fēng)化殼儲(chǔ)層和基底下方發(fā)育的網(wǎng)狀縫。

(4)裂縫疊前地震預(yù)測(cè)技術(shù)對(duì)資料品質(zhì)要求較高?；凇皟蓪捯桓摺辈杉卣鹳Y料可提高預(yù)測(cè)精度，但面對(duì)其產(chǎn)生的海量地震數(shù)據(jù)還需發(fā)展新的解釋技術(shù)和數(shù)據(jù)壓縮方法。目前大數(shù)據(jù)時(shí)代已經(jīng)來臨，期待更多的人工智能方法能與地震處理、解釋方法相結(jié)合，解決五維地震數(shù)據(jù)量大的問題。