印興耀,馬正乾,向 偉,宗兆云

(1.中國(guó)石油大學(xué)(華東)地球科學(xué)與技術(shù)學(xué)院,山東青島266580;2.海洋國(guó)家實(shí)驗(yàn)室海洋礦產(chǎn)資源評(píng)價(jià)與探測(cè)技術(shù)功能實(shí)驗(yàn)室,山東青島266071)

地層裂縫參數(shù)研究在油氣田勘探開(kāi)發(fā)、水利工程、土木工程等領(lǐng)域具有重要意義。隨著全球油氣資源勘探程度的不斷提高以及需求量的不斷增大,油氣勘探開(kāi)發(fā)從常規(guī)油氣儲(chǔ)層轉(zhuǎn)向致密砂巖、復(fù)雜碳酸鹽巖、火山巖、頁(yè)巖等復(fù)雜油氣藏,這些油氣藏都涉及到裂縫參數(shù)預(yù)測(cè)。在裂縫型油氣藏中,由于定向構(gòu)造應(yīng)力作用,自然裂縫通常呈現(xiàn)一種定向排列現(xiàn)象,裂縫體的走向、發(fā)育密度、充填物等裂縫參數(shù)與油氣運(yùn)移、聚集、成藏、儲(chǔ)量有著密切的聯(lián)系,另外裂縫參數(shù)也是儲(chǔ)層工程“甜點(diǎn)”的主要評(píng)價(jià)參數(shù)之一。所以儲(chǔ)層裂縫參數(shù)的定量預(yù)測(cè)和評(píng)價(jià)是油氣勘探開(kāi)發(fā)的重要環(huán)節(jié)之一。

目前,裂縫型儲(chǔ)層特征分析和預(yù)測(cè)技術(shù)可以分為三大類(lèi),分別是地質(zhì)巖心綜合分析預(yù)測(cè)法、測(cè)井識(shí)別技術(shù)以及地震預(yù)測(cè)技術(shù)。地質(zhì)巖心綜合分析法主要是通過(guò)地表露頭、巖心直接觀察裂縫的發(fā)育狀態(tài),統(tǒng)計(jì)裂縫參數(shù),進(jìn)而預(yù)測(cè)、標(biāo)定與之相對(duì)應(yīng)地層的裂縫分布[1-5]。這種技術(shù)簡(jiǎn)單、直接、高效,但露頭容易受風(fēng)化等外動(dòng)力地質(zhì)作用的影響而改變其原始地層賦存狀態(tài),且露出地表,原始地層應(yīng)力已經(jīng)釋放,裂縫狀態(tài)也相應(yīng)發(fā)生了變化。巖心裂縫狀態(tài)會(huì)受到取心機(jī)械外力的影響,且隨賦存應(yīng)力變化而發(fā)生變化,不能直接反映地層裂縫分布,而且取心率低,難以覆蓋全井段。測(cè)井技術(shù)識(shí)別裂縫分布,最初是利用常規(guī)測(cè)井曲線在裂縫層的異常響應(yīng)或者基于非線性理論的多曲線聯(lián)合組構(gòu)新特征參數(shù)曲線來(lái)識(shí)別井周裂縫參數(shù)[3,6],后來(lái)隨著成像測(cè)井技術(shù)的發(fā)展,可以獲取能夠直觀、真實(shí)地刻畫(huà)井壁周?chē)貙恿芽p產(chǎn)狀、開(kāi)度、密度等參數(shù)的圖像[7-9]。近些年發(fā)展的偶極聲波遠(yuǎn)探測(cè)技術(shù)在裂縫、斷層的探測(cè)中也有著良好的應(yīng)用效果,使測(cè)井探測(cè)實(shí)現(xiàn)了從近井筒到遠(yuǎn)井筒的突破,探測(cè)直徑提高到近百米[10-13]。測(cè)井技術(shù)具有較高的分辨率,是原始地層狀態(tài)的真實(shí)反映,但測(cè)井?dāng)?shù)據(jù)易受到鉆井液、井壁坍塌等因素的干擾,預(yù)測(cè)精度受到了一定的影響,成像測(cè)井成本高、解釋結(jié)果的主觀性大,且測(cè)井識(shí)別技術(shù)難以準(zhǔn)確地預(yù)測(cè)整個(gè)工區(qū)的裂縫分布。地震預(yù)測(cè)雖然精度不如測(cè)井技術(shù),但在探測(cè)整個(gè)工區(qū)或者深部地層裂縫賦存狀態(tài)方面,有著先天的優(yōu)勢(shì)[14-18]。

裂縫參數(shù)地震預(yù)測(cè)技術(shù)主要包括三維地震屬性預(yù)測(cè)技術(shù)[19-25]、多波多分量預(yù)測(cè)技術(shù)[26-27]、縱波方位各向異性預(yù)測(cè)技術(shù)[28-35]。利用三維地震屬性預(yù)測(cè)地層裂縫參數(shù),即利用地震數(shù)據(jù)體的一階導(dǎo)數(shù)——傾角[19]、二階導(dǎo)數(shù)——曲率[21,25,36-37]以及地震波相鄰道相似性——相干體[38-40]等幾何屬性實(shí)現(xiàn)地層裂縫刻畫(huà)。反射面曲率代表對(duì)應(yīng)地層的變形程度,某一層位曲率值越大,代表該層位褶皺程度越大,通常裂縫越發(fā)育,即曲率預(yù)測(cè)裂縫的基本假設(shè)是認(rèn)為地層構(gòu)造程度與裂縫發(fā)育程度之間存在單調(diào)映射關(guān)系,因此對(duì)于裂縫發(fā)育的構(gòu)造空白區(qū),曲率屬性則無(wú)法應(yīng)用。相干屬性發(fā)展主要經(jīng)歷了三代,基于三道互相關(guān)的C1算法[38],基于多道相似性的C2算法[38],基于本征結(jié)構(gòu)的相似算法C3[40]。C3算法較前兩種算法無(wú)論在分辨率上,還是抗噪性上都有著明顯的優(yōu)勢(shì),但其計(jì)算速度較慢。后來(lái)許多學(xué)者提出了一系列的改進(jìn)算法,如基于梯度結(jié)構(gòu)張量?jī)A角估計(jì)算法[41]、基于方位疊前地震數(shù)據(jù)的相干算法[42]等。多波多分量預(yù)測(cè)技術(shù),主要利用裂縫儲(chǔ)層中的橫波分裂現(xiàn)象[27]實(shí)現(xiàn)裂縫參數(shù)的預(yù)測(cè),但高昂的采集處理成本,純橫波勘探激發(fā)難等問(wèn)題極大地限制了其推廣應(yīng)用。轉(zhuǎn)換PS波(縱-橫轉(zhuǎn)換波)勘探使得多分量勘探有了進(jìn)一步的發(fā)展,但相較于縱波勘探,依舊存在高成本、對(duì)地震數(shù)據(jù)質(zhì)量要求高、對(duì)裂縫充填物不敏感等問(wèn)題。P波(縱波)方位各向異性裂縫預(yù)測(cè)技術(shù)主要利用裂縫誘導(dǎo)的縱波速度、振幅、頻率等屬性在觀測(cè)方位角域的變化規(guī)律實(shí)現(xiàn)裂縫密度、方位、填充物等參數(shù)預(yù)測(cè),包括基于縱波運(yùn)動(dòng)學(xué)特性的預(yù)測(cè)技術(shù)和基于縱波動(dòng)力學(xué)特性的預(yù)測(cè)技術(shù)。前者有NMO(normal moveout)速度方位變化裂縫預(yù)測(cè)技術(shù)[43-46]、正交地震測(cè)線縱波時(shí)差預(yù)測(cè)技術(shù)[47-49]和基于VSP數(shù)據(jù)的各向異性參數(shù)預(yù)測(cè)[50]等,這3項(xiàng)技術(shù)均是基于速度在方位角域和偏移距(入射角)域的規(guī)律變化實(shí)現(xiàn)裂縫參數(shù)預(yù)測(cè),但分辨率不足是其最難解決的問(wèn)題。后者主要利用地震P波振幅[51]、阻抗[52-54]、衰減[55-56]等屬性隨方位角的變化實(shí)現(xiàn)儲(chǔ)層預(yù)測(cè),具有更高的分辨率。整體來(lái)講,地震預(yù)測(cè)技術(shù)的探測(cè)深度大、范圍廣等特征是地質(zhì)巖心技術(shù)和測(cè)井技術(shù)所不具備的,而且基于P波屬性隨入射角和方位角變化(五維地震數(shù)據(jù))的技術(shù)較其它裂縫地震預(yù)測(cè)技術(shù)具有更高的分辨率和準(zhǔn)確度,可獲得更多的儲(chǔ)層信息且縱波數(shù)據(jù)更易獲取。與此同時(shí),裂縫儲(chǔ)層的等效地震巖石物理理論是連接儲(chǔ)層微觀參數(shù)與宏觀地震響應(yīng)之間的橋梁,也是反射系數(shù)參數(shù)化和裂縫參數(shù)穩(wěn)定預(yù)測(cè)的理論基礎(chǔ)。所以我們將從裂縫型儲(chǔ)層等效巖石物理理論和五維地震數(shù)據(jù)裂縫預(yù)測(cè)技術(shù)兩個(gè)方面介紹巖石物理驅(qū)動(dòng)的裂縫預(yù)測(cè)技術(shù)研究現(xiàn)狀與進(jìn)展。本文主要綜述裂縫型儲(chǔ)層等效各向異性巖石物理理論的研究現(xiàn)狀與進(jìn)展。五維地震數(shù)據(jù)裂縫預(yù)測(cè)技術(shù)研究現(xiàn)狀與進(jìn)展將在第二部分(Ⅱ)中詳細(xì)介紹。

1 裂縫儲(chǔ)層等效各向異性巖石物理理論

地震巖石物理模型建立了地層巖石物性參數(shù)和彈性參數(shù)之間的關(guān)系,是地震反演和裂縫參數(shù)預(yù)測(cè)的橋梁,所以了解巖石物理裂縫模型,厘清各模型的理論前提和物理機(jī)理是利用地震資料預(yù)測(cè)裂縫參數(shù)的基礎(chǔ)[57]。所以,本節(jié)將從3個(gè)方面綜述裂縫型儲(chǔ)層地震巖石物理理論:①考慮孤立孔縫的地震巖石物理理論,包括Hudson理論[58-62]、Schoenberg線性滑動(dòng)理論[63-68]、Eshelby-Cheng模型[69]、Aniso-DEM模型(anisotropic differential effective model)[70]、Aniso-SCA模型(anisotropic self-consistent approximation model)[70]等;②衰減各向異性地震巖石物理理論,包括修正Shoenberg線性滑動(dòng)模型[71-75]、Chapman多尺度裂隙模型[76-77]等;③極限頻率下考慮孔縫連通性的地震巖石物理理論,包括Thomsen含等徑孔隙裂隙模型[78]、基于廣義Gassmann方程的飽和裂縫巖石模型[79-80]等。

1.1 考慮孤立孔縫的地震巖石物理理論

以靜力學(xué)為基礎(chǔ)導(dǎo)出的裂縫型儲(chǔ)層地震巖石物理理論,主要研究不考慮孔隙-裂隙或裂隙-裂隙之間流體流動(dòng)的含孤立孔縫巖石。該理論將裂縫作為巖石的包含物,采用等效平均思想,推導(dǎo)參考元的等效巖石模量。下面重點(diǎn)介紹巖石彈性模量解析表達(dá)的Hudson理論和Schoenberg線性滑動(dòng)理論。

1.1.1 Hudson理論

如圖1所示,HUDSON[58-59]基于平均波場(chǎng)散射理論研究了各向同性背景中包含一組水平硬幣型裂縫介質(zhì)的等效彈性剛度矩陣C,提出了Hudson理論(本文坐標(biāo)系統(tǒng)均采用Cartesian坐標(biāo)系,并規(guī)定z或x3軸正方向垂直向下,n為裂隙的法向單位向量):

圖1 Hudson模型裂縫分布[81]

C=C0+C1+C2+o(e2)

(1)

式中:C為裂縫介質(zhì)等效剛度張量;C0為各向同性背景剛度張量;C1,C2分別是由裂縫引起的剛度張量的一階校正項(xiàng)和二階校正項(xiàng);o(e2)為裂隙密度e2的高階無(wú)窮小量。剛度矩陣C的非零元素分別為:

(2a)

(2b)

(2c)

(2d)

(2e)

其中,Cij代表等效介質(zhì)剛度系數(shù);λb,μb為各向同性背景拉梅參數(shù);e=υca3為裂縫密度且υc為單位體積巖石中裂縫的條數(shù),a為硬幣型裂縫的半長(zhǎng)軸長(zhǎng)度。

Uii的物理意義是硬幣型水平裂縫表面受到i方向的單位應(yīng)力時(shí),在i方向產(chǎn)生的位移不連續(xù)性,其值由裂縫自身的參數(shù)以及內(nèi)部充填物的性質(zhì)決定。由公式(2)可知,Uii真正獨(dú)立的分量為U11和U22,可表達(dá)為[59]U11=(16/3)[(λb+2μb)/(3λb+4μb)](1+M)-1和U33=(4/3)[(λb+2μb)/(λb+μb)](1+K)-1。當(dāng)裂縫中充填空氣時(shí),M=0,K=0;當(dāng)裂縫中充填無(wú)粘性液體時(shí),M=0,K→∝?U33=0;當(dāng)裂縫中充填弱介質(zhì)或粘性流體時(shí),M=[4μ′(λb+2μb)]/[παratioμb(λb+μb)],K={[k′+(4/3)μ′)](λb+2μb)}/[παratioμb(λb+μb)],其中αratio代表裂隙橫縱比,k′和μ′分別為充填物的體積模量和剪切模量。

Hudson理論有著較為苛刻的假設(shè)條件:①裂縫為硬幣型,縱橫比較小;②裂縫單組、定向排列;③裂縫稀疏分布于各向同性背景中且裂縫之間孤立,即忽略裂縫之間的彈性互動(dòng)及流體流動(dòng);④裂縫尺度遠(yuǎn)小于地震波長(zhǎng)。在此基礎(chǔ)上,HUDSON等[60-61,82]繼續(xù)發(fā)展該模型,不斷突破假設(shè)條件的限制:針對(duì)假設(shè)②,提出了各向同性背景中發(fā)育多組定向裂縫系統(tǒng)的等效模型[60];針對(duì)假設(shè)③,推導(dǎo)了飽和流體的連通裂隙-裂隙和孔隙-裂隙介質(zhì)模型[82]和各向異性背景中包含定向裂縫模型[61],后者利用各向同性介質(zhì)包含虛擬定向裂縫體來(lái)模擬各向異性背景,這是基于弱各向異性背景假設(shè),本質(zhì)上并沒(méi)有考慮真實(shí)裂縫系統(tǒng)與各向異性背景之間的作用,所以對(duì)于裂縫型頁(yè)巖,由于強(qiáng)各向異性頁(yè)巖背景,該模型很難適用。另外,HUDSON[83]還研究了不飽和流體填充時(shí),孤立裂縫介質(zhì)的彈性性質(zhì);LIU等[84]進(jìn)一步考慮裂縫表面形態(tài),提出了孤立滑移面裂隙模型、非連續(xù)接觸裂隙模型和液體薄層裂隙模型;TOD等[85]針對(duì)裂隙面與裂縫面不共面但裂隙中心位于裂縫表面的情況,提出了裂縫法向與裂隙法向不平行的裂縫模型。

1.1.2 Schoenberg線性滑動(dòng)理論

SCHOENBERG忽略了裂縫的形狀和微觀結(jié)構(gòu),將其看作一個(gè)非常薄的滑動(dòng)面且表面位移不連續(xù),這種位移不連續(xù)性與穿過(guò)裂縫的應(yīng)力呈線性關(guān)系?；瑒?dòng)面與背景介質(zhì)按照體積比例構(gòu)成層狀介質(zhì),并假設(shè)滑動(dòng)面的剛度模量與其厚度為等價(jià)無(wú)窮小,然后利用平滑方法實(shí)現(xiàn)裂縫介質(zhì)等效模量的求取[63,65-66,86-87],形成Schoenberg線性滑動(dòng)理論。由該理論發(fā)展的模型是一種適用于具有線性連續(xù)邊界、充滿弱強(qiáng)度(充填物模量小)介質(zhì)的平行層模型。按照Schoenberg線性滑動(dòng)理論,含單組平行定向裂縫巖石(如圖2a所示)的有效柔度張量可寫(xiě)為[67]:

(3)

式中:下標(biāo)b代表背景介質(zhì);Sijkl表示巖石等效柔度張量的第ijkl個(gè)元素;Zik是裂縫柔度矩陣的第ik個(gè)元素;ni是裂縫法向矢量的第i個(gè)分量。介質(zhì)的柔度張量即為其剛度張量的逆S=C-1。對(duì)于包含多組定向裂縫系統(tǒng)的巖石,如圖2b所示,其等效柔度張量可進(jìn)一步表示為[67]:

圖2 Schoenberg線性滑動(dòng)模型及其推廣

(4)

假設(shè)一組稀疏、旋轉(zhuǎn)不變的垂直裂縫鑲嵌在各向同性背景巖石中,可用法向、切向柔度ZN,ZT對(duì)裂縫進(jìn)行完全表征,當(dāng)裂縫法向與x軸一致時(shí),裂縫巖石等效剛度矩陣CHTI可用Voigt縮寫(xiě)形式表達(dá)為[67]:

(5)

式中:Mb,λb和μb分別為各向同性背景巖石的縱波模量、第一拉梅參數(shù)和剪切模量;χb=1-2γb,γb=μb/Mb代表背景介質(zhì)橫縱波速度比的平方;δN=ZNMb/(1+ZNMb)和δT=ZTμb/(1+ZTμb)分別為裂縫法向和切向弱度參數(shù)。

綜上,Schoenberg線性滑動(dòng)理論的基本假設(shè)可以總結(jié)為:①長(zhǎng)波長(zhǎng)限制——長(zhǎng)波長(zhǎng)限制是指裂縫的大小遠(yuǎn)小于地震波的波長(zhǎng),相比裂縫間距,裂縫分布細(xì)節(jié)及張開(kāi)度可以忽略,也就是說(shuō),Schoenberg模型將裂縫看成滿足線性光滑邊界條件的無(wú)限薄且非常松軟的地層或平面,忽略裂縫的形狀和微結(jié)構(gòu);②對(duì)于有多套裂縫存在的情況,忽略裂縫間的相互作用;③Schoenberg模型中假設(shè)裂縫面為位移間斷,但是界面兩側(cè)應(yīng)力保持連續(xù)。

實(shí)驗(yàn)證明Schoenberg理論中的裂縫和Hudson理論中的裂隙在地震響應(yīng)上是分不開(kāi)的[59],(1)式和(5)式是等效的,可得[65]:

(6a)

(6b)

式中:αratio=c/a表示裂縫橫縱比。設(shè)置μ′=0,可研究飽和孤立裂縫巖石弱度參數(shù)特征;設(shè)置μ′=κ′=0,可研究干裂縫巖石弱度參數(shù)特征。

1.2 衰減各向異性地震巖石物理理論

對(duì)于包含飽和、相互連通的定向裂縫及背景孔隙的地下巖石,當(dāng)?shù)卣鸩ǖ竭_(dá)時(shí),會(huì)導(dǎo)致孔縫中的流體產(chǎn)生振蕩擴(kuò)散[88],這種波致流將誘導(dǎo)地震波衰減各向異性,發(fā)生頻散現(xiàn)象,表現(xiàn)出動(dòng)力學(xué)特征。所以考慮這種波致流衰減的巖石物理模型在本文中稱為衰減各向異性地震巖石物理理論,該類(lèi)理論已經(jīng)被多位學(xué)者研究,并給出了各種情況下的巖石彈性模量復(fù)值表達(dá)式。本節(jié)將主要綜述修正的Shoenberg線性滑動(dòng)模型及相關(guān)理論和Chapman多尺度裂隙模型。

1.2.1 修正的Schoenberg線性滑動(dòng)模型

(7a)

(7b)

(7c)

(8a)

(8b)

① 對(duì)于機(jī)制1,如圖3a所示:

(9a)

(9b)

式中:ηf和κf分別表示裂縫流體的粘滯系數(shù)和體積模量;Kr代表巖石滲透率;ω和k表示角頻率和波數(shù)。

② 對(duì)于機(jī)制2,如圖3b所示:

(10a)

(10b)

式中:J2=φmκfKm/2ωηf(J>0);φm和Km分別為無(wú)裂縫背景基質(zhì)的孔隙度和滲透率。

③ 對(duì)于機(jī)制3,如圖3c所示:

圖3 裂縫型巖石3種不同流體流動(dòng)模型

(11a)

式中:qm和ηm分別代表部分飽和裂縫m相流體的飽和度和粘滯系數(shù),m=l或g,表示液態(tài)相或氣態(tài)相,ql+qg=1。此時(shí),公式(8b)可以改寫(xiě)為:

(11b)

其中,

(11c)

(11d)

并且:

Fl(q)=F1(q)≈A1(1-q)[1+cos(π-πq)]

(11e)

Fg(q)=F2(q)≈A2(1-q)[1+cos(π-πq)]

(11f)

當(dāng)液態(tài)相占據(jù)裂縫中心位置時(shí),A1=0.053,A2=0.058;當(dāng)液態(tài)相分布在裂縫邊緣時(shí),F1(q)=F2(q),Fg(q)=Fl(q)。

1.2.2Chapman多尺度裂縫模型

不同尺度裂縫的定向排列均能誘導(dǎo)地層呈現(xiàn)各向異性特征,但只有足夠大尺度的裂縫才對(duì)油氣的運(yùn)移、成藏及生產(chǎn)起主要作用[89-90]。常規(guī)等效介質(zhì)理論采用裂縫密度、縱橫比等參數(shù)來(lái)考慮裂縫對(duì)等效彈性模量的影響,并未研究裂縫尺度與彈性模量的關(guān)系。多尺度非均勻性的引入誘導(dǎo)地震反射/透射的頻散現(xiàn)象[91-92],所以多尺度定向裂縫的引入也會(huì)產(chǎn)生各向異性的頻變特征,LIU等[93]和MAULTZSCH等[94]在實(shí)際數(shù)據(jù)處理中已經(jīng)發(fā)現(xiàn)了各向異性與頻率相關(guān)的現(xiàn)象。為了分析多尺度裂縫的地震各向異性響應(yīng)特征,實(shí)現(xiàn)多尺度裂縫的識(shí)別,CHAPMAN[77]提出頻變多尺度裂縫巖石物理模型,接下來(lái)對(duì)該模型進(jìn)行簡(jiǎn)要闡述。

如圖4所示,Chapman多尺度裂縫模型[76-77]假設(shè)巖石中含球狀等徑孔隙、隨機(jī)排列的硬幣型微裂隙以及定向排列的硬幣型裂縫,其中基質(zhì)等徑孔隙和微裂隙與巖石顆粒尺度相當(dāng),而定向排列的裂縫尺度遠(yuǎn)大于顆粒尺度,但小于地震波長(zhǎng)尺度。同時(shí),假設(shè)微裂隙與微裂隙間、背景等徑孔隙與微裂隙間相互連通構(gòu)成單元c1,定向排列的裂縫與多個(gè)等徑孔隙或微裂隙連通構(gòu)成單元c2,但每個(gè)等徑孔隙或微裂隙最多與一條定向排列的裂縫連通,且定向裂縫之間不連通。當(dāng)?shù)卣鸩ǖ竭_(dá)含孔隙裂隙巖石時(shí),認(rèn)為壓力梯度將誘導(dǎo)兩種不同尺度流體擴(kuò)散效應(yīng),即定向排列的裂縫與背景等徑孔隙、微裂隙之間的達(dá)西流效應(yīng),及微裂隙之間或與背景等徑孔隙之間的噴射流效應(yīng)。因此Chapman模型的等效剛度參數(shù)Cchapman可表示為:

圖4 Chapman裂隙模型[95]

(12)

中觀尺度定向裂縫和微觀尺度孔裂隙的存在引入兩個(gè)尺度的弛豫時(shí)間參數(shù)τc和τm,且兩者存在近似函數(shù)關(guān)系:

(13)

式中:ac和ag分別是定向排列裂縫與顆粒尺寸。公式(13)物理含義為:隨著裂縫半徑增大,裂縫表面積與體積的比值變小,達(dá)到流體壓力平衡時(shí)其單位表面積將有更多的流體通過(guò),所需時(shí)間也更長(zhǎng);而與定向裂縫相關(guān)的弛豫時(shí)間τc將導(dǎo)致地震頻率范圍內(nèi)出現(xiàn)速度的頻散與衰減,也即定向裂縫誘導(dǎo)的各向異性具有頻率相關(guān)性。由于實(shí)際應(yīng)用中,無(wú)縫多孔背景的模量較易獲得,為此,CHAPMAN等[96]將公式(12)進(jìn)行擴(kuò)展,得:

(14)

式中:f是頻率;固體顆粒拉梅常數(shù)分別為λg=Γ-Φmc,p(λg,μg,fω)與μg=Z-Φmc,p(λg,μg,fω),其中,Γ和Z是頻率fω下經(jīng)背景等徑孔隙與微裂隙校正后的拉梅常數(shù),Φmc,p是校正項(xiàng)。對(duì)于大多數(shù)的裂縫型儲(chǔ)層,微裂隙孔隙度要遠(yuǎn)遠(yuǎn)小于等徑孔隙度。當(dāng)勘探頻率低于微觀噴射流頻率時(shí),可忽略公式(14)中的微裂隙項(xiàng)。

CHAPMAN[97]在單組定向多尺度裂縫頻變彈性特征的研究基礎(chǔ)上,提出了多組定向裂縫頻率依賴的巖石物理模型。主要解決思路為:首先求出由于等徑孔隙存在而產(chǎn)生的彈性參數(shù)修正量,然后求出本征坐標(biāo)系下,每組裂縫單獨(dú)存在時(shí)產(chǎn)生的彈性修正量,接著通過(guò)坐標(biāo)旋轉(zhuǎn)將裂縫引起的修正量變換至觀測(cè)坐標(biāo)系下,最后將背景基質(zhì)張量、等徑孔隙修正量、多組裂縫修正量求和,便可得到最終的彈性張量,可表達(dá)為:

(15)

1.3 極限頻率下考慮孔縫連通性的地震巖石物理理論

上節(jié)表明:地球巖石受到外界波場(chǎng)擾動(dòng)時(shí),巖石內(nèi)流體在壓力梯度的作用下會(huì)發(fā)生振蕩擴(kuò)散,且這種擴(kuò)散特征與波場(chǎng)頻率有關(guān)。地震頻率相對(duì)較低,一般認(rèn)為在地震波擾動(dòng)下,巖石內(nèi)流體在連通孔縫間有充足時(shí)間流動(dòng)擴(kuò)散以達(dá)到壓力均衡,表現(xiàn)為極限頻率下的動(dòng)力學(xué)特征。所以本節(jié)將以極限頻率下考慮孔縫連通性的地震巖石物理理論為主題,重點(diǎn)介紹Thomsen含等徑孔隙裂隙模型和基于廣義Gassmann方程的飽和裂縫巖石模型。

1.3.1Thomsen含等徑孔隙裂隙模型

Thomsen將裂隙看作一個(gè)橫縱比較小的橢球體,然后假設(shè)一組水平定向裂隙稀疏地分布在各向同性多孔巖石中,裂隙與背景等徑孔隙之間是連通的,在壓力梯度作用下,流體發(fā)生振蕩擴(kuò)散,如圖5所示。在先前弱各向異性參數(shù)的研究基礎(chǔ)上[98],THOMSEN[78]提出了考慮流體壓力均衡的孔隙裂隙介質(zhì)模型:

(16a)

(16b)

(16c)

(16d)

此時(shí)裂縫與背景孔隙之間有充足的時(shí)間通過(guò)流體流動(dòng)達(dá)到局部壓力均衡;在中高頻條件下,背景孔隙與裂縫之間沒(méi)有足夠的時(shí)間達(dá)到流體壓力均衡,則有:

(16e)

式中:flow和fm-h分別代表低頻和中高頻條件;φc為裂隙孔隙度;φp為背景等徑孔隙度。

1.3.2 基于廣義Gassmann方程的飽和裂縫巖石模型

GASSMANN[99]提出了等效的各向異性巖石流體替換方程[100]:

(17a)

i,j=1,2,…,6

(17b)

(17c)

對(duì)于垂直定向裂縫誘導(dǎo)的水平橫向各向同性介質(zhì)(horizontally transverse isotropy,HTI),聯(lián)合Schoenberg線性滑動(dòng)模型和廣義Gassmann方程,公式(17b)和公式(17c)可表達(dá)為[79,102]:

(17d)

(17e)

(17f)

(18a)

(18b)

(18c)

(18d)

(18e)

(18f)

(18g)

(18h)

(18i)

(18j)

其它剛度系數(shù)均為零。由于調(diào)用了廣義Gassmann方程,因此公式(18)是低頻假設(shè)下,含連通孔隙裂縫巖石的等效剛度系數(shù),即要求在地震波的擾動(dòng)下,背景等徑孔隙與定向裂縫中的流體有充足的時(shí)間達(dá)到壓力均衡。GUREVICH[79]研究表明,由于該模型考慮了背景孔隙和裂縫之間的流體流動(dòng),所以不再滿足HTI型巖石總?cè)岫仁歉飨蛲员尘叭岫扰c裂縫柔度的疊加。CIZ等[103]通過(guò)假設(shè)在加載圍壓和孔隙空間應(yīng)力時(shí)巖石骨架和孔隙充填材料變形較小,將各向異性流體替換方程推廣為各向異性固體替換方程,以實(shí)現(xiàn)孔隙空間充填固體或固-液混合物的各向異性儲(chǔ)層剛度參數(shù)的估算,可用柔度張量Sijkl表達(dá)為:

(19)

式中:上標(biāo)sat,dry和0分別代表與飽和巖石、無(wú)充填多孔骨架和背景基質(zhì)相關(guān)的量。上標(biāo)if和φ分別代表與孔隙充填物體平均應(yīng)變有關(guān)的量和與骨架孔隙空間有關(guān)的量,且CIZ等[103]也給出了Sif和Sφ的估計(jì)方式。

巖石物理模型是利用五維地震數(shù)據(jù)預(yù)測(cè)裂縫參數(shù)的理論基礎(chǔ)。本文根據(jù)儲(chǔ)層流體在地震波擾動(dòng)下的響應(yīng)特征,從3個(gè)角度分別介紹了當(dāng)前國(guó)內(nèi)外常用的裂縫型儲(chǔ)層地震巖石物理理論。每種模型都有其適用的地質(zhì)環(huán)境:①對(duì)于喉道不發(fā)育的儲(chǔ)層,當(dāng)微裂隙/裂縫飽和均勻流體時(shí),可以采用孤立孔縫Hudson理論[58-59]、線性滑動(dòng)模型[65,67-68]或高頻Thomsen含等徑孔隙裂隙模型[78]來(lái)模擬裂縫、流體等對(duì)巖石屬性的影響;當(dāng)裂隙/裂縫充填多相流體時(shí),可以采用部分飽和Hudson模型[83]等;②對(duì)于喉道充分發(fā)育的裂縫型地層,在地震數(shù)據(jù)超低頻假設(shè)下,可以采用基于廣義Gassmann方程的飽和裂縫巖石模型[79,102]或低頻Thomsen含等徑孔隙裂隙模型[78]計(jì)算地層彈性參數(shù),但研究巖石的等效彈性模量在全頻帶的特征時(shí),應(yīng)該采用連通孔-裂隙的Hudson模型[82,85]等。此時(shí),我們并沒(méi)有強(qiáng)調(diào)基質(zhì)孔隙(或硬孔隙)的發(fā)育情況,但可以通過(guò)調(diào)節(jié)每種模型的相關(guān)參數(shù)來(lái)改變基質(zhì)孔隙的含量以適應(yīng)更廣泛的地質(zhì)情況。另外,當(dāng)上述兩類(lèi)儲(chǔ)層發(fā)育多組相同尺度裂縫體時(shí),多組裂隙Hudson模型[60]、線性滑動(dòng)模型[104]等將是較好的選擇。但對(duì)于基質(zhì)孔隙-隨機(jī)裂隙-定向裂縫發(fā)育的儲(chǔ)層,地震波經(jīng)過(guò)時(shí)誘導(dǎo)多尺度流體流動(dòng),此時(shí)可采用Chapman波致流衰減理論等多尺度裂縫模型[77,105]來(lái)研究?jī)?chǔ)層宏觀彈性特征?？傊?地下介質(zhì)非均質(zhì)性強(qiáng),孔隙、裂縫結(jié)構(gòu)多樣,油氣水關(guān)系復(fù)雜,一種巖石物理理論難以全面刻畫(huà)地下巖石的彈性特征,應(yīng)該針對(duì)不同的地質(zhì)環(huán)境,研究和構(gòu)建適用的地震巖石物理模型。

2 裂縫型儲(chǔ)層地震響應(yīng)特征分析

研究不同巖石物理模型之間的聯(lián)系和差異,分析裂縫型儲(chǔ)層微觀因素對(duì)宏觀地震響應(yīng)(反射系數(shù))的影響,將為五維地震裂縫參數(shù)預(yù)測(cè)方法及穩(wěn)定性研究提供理論基礎(chǔ)。所以本節(jié)在前文介紹的巖石物理理論基礎(chǔ)上,主要討論Hudson模型、Thomsen含等徑孔隙裂隙模型和基于廣義Gassmann方程的飽和裂縫巖石模型(后面將其簡(jiǎn)稱為Gurevich模型)間的聯(lián)系和差異,并分析裂縫、流體及背景孔隙度等參數(shù)對(duì)地震響應(yīng)特征(反射系數(shù))的影響,從而為裂縫和流體參數(shù)反演策略的制定提供思路。

為了便于分析裂縫誘導(dǎo)各向異性特征,我們?cè)O(shè)計(jì)一個(gè)兩層模型,如圖6所示,如無(wú)特殊說(shuō)明,上層介質(zhì)始終為各向同性致密介質(zhì),下層為發(fā)育定向垂直裂縫的儲(chǔ)層,裂縫法向與x軸相同,上下層介質(zhì)的礦物組成相同且分布均勻。孔-裂隙空間中飽含均勻流體(水的體積模量2.2GPa,油的體積模量1.4GPa,氣的體積模量0.133GPa)。表1展示了兩類(lèi)不同礦物的速度和密度參數(shù),這兩類(lèi)礦物的區(qū)別是橫縱波速度比的平方不一樣,礦物A為0.475,礦物B為0.350。采用RüGER[106]給出的反射系數(shù)方程研究裂縫密度、充填物和孔隙度等參數(shù)對(duì)地震響應(yīng)的影響,方程表達(dá)如下:

表1 兩類(lèi)礦物的速度和密度參數(shù)

圖6 兩層裂縫模型示意(礦物組成相同)

(20a)

(20b)

(20c)

公式(20)中,G=ρβ2和Z=ρα。R代表反射系數(shù),θ為入射角,φ為方位角。α和β分別為垂直方向的縱、橫波速度,ρ為密度,Δδ(V),Δε(V)和Δγ為上下兩層介質(zhì)Thomsen各向異性參數(shù)差值。接下來(lái),我們將從背景無(wú)等徑孔隙和背景包含等徑孔隙兩種情況展開(kāi)討論。

2.1 背景無(wú)等徑孔隙的巖石(背景孔隙度假設(shè)為0)

假設(shè)背景介質(zhì)僅由均勻的礦物A組成,則上下層介質(zhì)的差異僅由垂直裂縫和充填的流體而產(chǎn)生。采用Hudson孤立裂隙模型、Thomsen含等徑孔隙裂隙模型和Gurevich模型分別計(jì)算干巖石和飽含水巖石的剛度模量,并代入公式(20)可計(jì)算反射系數(shù)。圖7展示了由Gurevich模型和公式(20)計(jì)算的無(wú)背景孔隙的兩層模型反射系數(shù)隨入射角和方位角變化的三維圖。圖7a和圖7b中裂縫分別無(wú)充填(干裂縫)和飽含水,沿紅色箭頭方向,下層裂縫密度依次為0.001,0.010,0.025和0.050。從圖7中可以發(fā)現(xiàn),裂縫引起的方位振幅變化較弱,且隨著充填物模量的增大,方位振幅逐漸降低。圖8和圖9展示了無(wú)背景孔隙的兩層模型反射系數(shù)隨方位角和入射角變化的二維曲線。深粉紅色、青藍(lán)色、綠色和深藍(lán)色曲線分別代表利用低頻Thomsen理論、中高頻Thomsen理論、Gurevich模型和Hudson模型計(jì)算得到的反射系數(shù),圖8和圖9對(duì)應(yīng)的下層裂縫密度分別為0.001和0.050,裂縫縱橫比均為0.005。

圖7 無(wú)背景孔隙的兩層模型反射系數(shù)隨入射角θ和方位角φ變化(背景礦物為礦物A,沿紅色箭頭方向,裂縫密度依次為0.001,0.010,0.025和0.050)

對(duì)于干裂縫情況,從圖8a、圖8b和圖9a、圖9b中我們可以發(fā)現(xiàn),Gurevich模型與Hudson模型的結(jié)果一致,但與高、低頻Thomsen模型的結(jié)果存在差異,同時(shí),高頻與低頻Thomsen模型的結(jié)果相同。這是由于Gurevich模型是在Hudson干裂縫模型基礎(chǔ)上采用各向異性Gassmann流體替換而推導(dǎo)的,所以對(duì)于干裂縫,二者是一致的。但此時(shí),Thomsen理論給出的各向異性參數(shù)表達(dá)是Hudson理論給出的結(jié)果關(guān)于裂縫密度的線性近似[107]。因此小裂縫密度情況下,Hudson理論、Gurevich模型和Thomsen理論比較接近,正如圖8a和圖8b所示;大裂縫密度情況下,Hudson理論、Gurevich模型和Thomsen理論的差異較大,這也被圖9a和圖9b證實(shí)。另外,由于干燥巖石不存在流體流動(dòng)現(xiàn)象,所以高頻與低頻Thomsen模型的結(jié)果相同。

圖8 無(wú)背景孔隙的兩層模型反射系數(shù)對(duì)比(下層介質(zhì)裂縫密度為0.001)

圖9 無(wú)背景孔隙的兩層模型反射系數(shù)對(duì)比(下層介質(zhì)裂縫密度為0.050)

對(duì)于飽含水的巖石,分析圖8d和圖9d可以發(fā)現(xiàn),當(dāng)小裂縫密度時(shí),基于Hudson理論、Gurevich模型和Thomsen理論的反射系數(shù)差異較小;對(duì)于大裂縫密度,Hudson理論和Gurevich模型所得結(jié)果的差異較小,但與Thomsen理論給出的反射系數(shù)有較大差異。這是由于在小裂縫密度情況下,流體的作用大于裂縫的影響,盡管Gurevich模型與低頻Thomsen理論假設(shè)激勵(lì)頻率足夠低,但定向裂隙的形狀相同,所以在地震波長(zhǎng)范圍內(nèi)孔隙壓力均衡,不會(huì)發(fā)生裂隙之間的流體流動(dòng),其表現(xiàn)與孤立裂縫情況相似[71]。因此Gurevich模型、Hudson模型和Thomsen理論在小裂縫密度情況下表現(xiàn)相近。在大裂縫密度時(shí),裂縫影響大于充填水的作用,所以Thomsen理論的線性近似導(dǎo)致其與Gurevich模型和Hudson模型間的差異變大。對(duì)比圖8c和圖9c可以發(fā)現(xiàn),無(wú)論裂縫密度大小,在裂縫走向上(觀測(cè)方位90°),Hudosn模型與Gurevich模型存在差異,且Gurevich模型與Thomsen理論更接近。這是由于裂縫走向的反射系數(shù)由上下層垂直速度和密度的差異決定,與各向異性參數(shù)無(wú)關(guān)(根據(jù)公式(20)可知),流體充填弱化裂縫的作用,Gurevich模型與Thomsen理論考慮流體作用的方式相同,所以導(dǎo)致圖8c和圖9c所示的現(xiàn)象。這種現(xiàn)象是理論模型差異而導(dǎo)致,并非地層屬性所造成。

綜上,對(duì)于背景孔隙度極其小的地層(例如頁(yè)巖地層),地震波的傳播不再引起孔-裂隙流體的流動(dòng),和高頻極限情況一致,所以無(wú)論喉道結(jié)構(gòu)發(fā)育情況,我們均可采用孤立裂隙模型(Hudson模型和高頻Thomsen理論)研究地層的孔彈性特征。且當(dāng)裂縫密度較小時(shí),利用Hudson理論、Gurevich模型和Thomsen理論均可獲得較好的理論模擬,但當(dāng)裂縫密度較大時(shí),Thomsen理論的模擬結(jié)果精度較Hudson理論和Gurevich模型的精度低。

在上述分析基礎(chǔ)上,我們采用Hudson模型研究的裂縫密度、充填物類(lèi)型對(duì)地震反射特征的影響。圖10 和圖11展示了礦物A和礦物B分別作為基質(zhì)礦物時(shí)的反射系數(shù)R(θ,φ)隨方位角φ和入射角θ的變化,其中每種線型代表不同的裂縫密度,依次為0.01,0.02,0.03,0.05,裂縫橫縱比為0.005,圖10a、圖10c和圖10e以及圖11a、圖11c和圖11e代表裂縫走向反射系數(shù)相對(duì)垂直入射反射系數(shù)的改變量(即R(θ,90°)-R(0,90°))隨入射角變化,圖10b、圖10d和圖10f以及圖11b、圖11d和圖11f代表入射角為30°反射系數(shù)相對(duì)裂縫法向反射系數(shù)的改變量(即R(30°,φ)-R(30°,0))隨方位角變化。圖10a、圖10b、圖11a和圖11b為氣充填裂縫,圖10c、圖10d、圖11c和圖11d為油充填,圖10e、圖10f、圖11e 和圖11f為水充填。圖10和圖11 顯示地震反射振幅隨觀測(cè)方位角呈現(xiàn)余弦變化規(guī)律,且隨著裂縫密度的增加,由定向裂縫引起的AVO梯度及振幅方位變化幅度均會(huì)增加。這是由于裂縫的存在會(huì)降低巖石中波的傳播速度,進(jìn)而增大界面上下巖石的彈性差異,從公式(20b)可得,AVO梯度將會(huì)增加;公式(20c)表明反射振幅方位變化幅度是各向異性參數(shù)的線性函數(shù),各向異性參數(shù)是裂縫密度的線性函數(shù)[108],則反射振幅方位變化幅度是裂縫密度的近似線性函數(shù),所以裂縫密度的增加會(huì)提高振幅方位變化幅度。很多油氣勘探工作者已經(jīng)利用這種特性開(kāi)展裂縫密度預(yù)測(cè)方法研究。

圖10 反射系數(shù)改變量隨入射角θ和方位角φ的變化(礦物A組成背景基質(zhì),裂縫密度由小到大依次為0.01,0.02,0.03,0.05,裂縫橫縱比為0.005)

圖11 反射系數(shù)改變量隨入射角θ和方位角φ的變化(礦物B組成背景基質(zhì),各個(gè)圖件的描述與圖10相同)

綜上,對(duì)于背景孔隙度非常小或孔喉結(jié)構(gòu)發(fā)育較差的地層,盡管裂縫密度的增大將會(huì)導(dǎo)致地震反射振幅方位變化幅度的增大,但裂縫充填物類(lèi)型也會(huì)影響反射振幅方位變化幅度,所以直接利用反射振幅方位變化幅度難以準(zhǔn)確預(yù)測(cè)裂縫密度與流體類(lèi)型。為了準(zhǔn)確預(yù)測(cè)裂縫參數(shù)、實(shí)現(xiàn)油氣識(shí)別,需要研究如何從反射振幅方位變化信息中解耦背景基質(zhì)、裂縫及油氣等信息,以提高裂縫型儲(chǔ)層預(yù)測(cè)精度。除此,雖然反射振幅方位變化幅度與充填物類(lèi)型之間的關(guān)系依賴橫縱波速度比等參數(shù),但從圖10和圖11中可以發(fā)現(xiàn)充填物的模量越小——從水充填到氣充填,裂縫走向(90°觀測(cè)方位)反射振幅與裂縫傾向(0觀測(cè)方位)反射振幅的差(即ΔR)將會(huì)越大。這或許為流體識(shí)別指明一條好的方向——在獲得裂縫方位的前提下,可利用五維地震數(shù)據(jù)直接檢測(cè)裂縫流體類(lèi)型。

2.2 背景基質(zhì)發(fā)育等徑孔隙的巖石

假設(shè)上下層基質(zhì)均由礦物A組成,下層背景介質(zhì)包含10%的等徑孔隙,此時(shí)如果孔喉發(fā)育良好,在低頻激勵(lì)下,下層介質(zhì)流體將在定向裂縫和背景孔隙之間流動(dòng)以平衡局部壓力。我們依舊采用Hudson孤立裂隙模型、Thomsen含等徑孔隙裂隙模型和Gurevich模型分別計(jì)算飽含氣和飽含水巖石的剛度模量,并代入公式(20)計(jì)算反射系數(shù)。圖12和圖13展示了反射系數(shù)隨方位角和入射角變化的曲線,圖中不同線型與圖8中線型代表的意義一致,且圖12和圖13對(duì)應(yīng)的下層裂縫密度分別為0.001和0.050,裂縫縱橫比均為0.005。

圖12 含背景等徑孔隙的兩層模型反射系數(shù)對(duì)比(下層介質(zhì)裂縫密度為0.001,背景孔隙度為10%)

Hudson孤立裂隙模型和高頻Thomsen模型均假設(shè)地層孔裂隙之間流體無(wú)法流動(dòng),低頻Thomsen模型和Gurevich模型均考慮了因局部壓力差而導(dǎo)致的孔-裂隙間流體流動(dòng)。比較圖9a和圖13a,可以發(fā)現(xiàn)孤立孔裂隙和連通孔裂隙介質(zhì)的地震響應(yīng)存在差異,即流體流動(dòng)對(duì)地層地震響應(yīng)的影響不可忽視,這也證明了地震方法預(yù)測(cè)儲(chǔ)層流體流動(dòng)性是可能的。圖13a中Gurevich模型和低頻Thomsen模型之間存在差異的原因之一是Thomsen模型對(duì)裂縫密度的線性近似,另一部分原因是Thomsen假設(shè)裂縫與流體引起的柔度擾動(dòng)矩陣為對(duì)角矩陣,而Gurevich模型中裂縫與流體引起的柔度擾動(dòng)矩陣不為對(duì)角矩陣。所以,從圖13a到圖13c,水充填孔裂隙替換氣充填孔裂隙,減弱了裂縫的影響,則圖13c中Gurevich模型和低頻Thomsen模型之間的差異較圖13a更小。圖12對(duì)應(yīng)的裂縫密度非常小,流體流動(dòng)效應(yīng)也非常微弱,所以圖12a和圖12c中,基于這4種模型的地震響應(yīng)差異非常小。

為了分析流體流動(dòng)等因素對(duì)地震方位響應(yīng)的影響,我們繪制地震振幅方位差異(任意觀測(cè)方位地震振幅與0觀測(cè)方位的地震振幅的差異)隨方位角變化的曲線,如圖12b,圖12d和圖13b,圖13d所示。我們可以發(fā)現(xiàn)孤立孔-裂隙介質(zhì)的方位振幅變化規(guī)律(Hudson模型和高頻Thomsen理論)與背景無(wú)等徑孔隙情況(圖10和圖11)相似,這是由于此時(shí)含等徑孔隙的背景可作為一個(gè)整體,作用等效于無(wú)等徑孔隙的背景。除此,流體流動(dòng)效應(yīng)對(duì)地震振幅方位響應(yīng)有著較大的影響。不同于孤立裂縫的情況,飽含流體連通裂縫情況的裂縫走向(90°觀測(cè)方位)反射振幅一直大于裂縫傾向(0觀測(cè)方位)反射振幅。這是因?yàn)榭紤]孔-縫流體流動(dòng)時(shí),γ不被影響,但ε(V)模值增大且達(dá)到主導(dǎo)地位,所以ΔR<0,即R(θ,0)

圖13 含背景等徑孔隙的兩層模型反射系數(shù)對(duì)比(下層介質(zhì)裂縫密度為0.050,其它描述與圖12相同)

仍然假設(shè)上下層基質(zhì)均由礦物A組成,下層介質(zhì)裂縫密度為0.050。首先假設(shè)裂縫橫縱比為0.005,采用Gurevich模型和公式(20)繪制不同等徑孔隙度下反射振幅方位差異隨觀測(cè)方位變化的曲線,如圖14 所示,其中圖14a代表孔裂隙充填氣,圖14b代表孔裂隙充填水。然后假設(shè)背景孔隙度固定,采用Gurevich模型和公式(20)繪制不同裂縫縱橫比下反射振幅方位差異隨觀測(cè)方位變化的曲線,如圖15所示,其中圖15a和圖15b代表背景孔隙度為0.005,圖15c 和圖15d代表背景孔隙度為0.050,圖15a和圖15c 代表孔裂隙充填氣,圖15b和圖15d代表充填水。根據(jù)圖14,我們可以發(fā)現(xiàn),對(duì)于含氣儲(chǔ)層(圖14a),當(dāng)背景孔隙度處于1%～10%時(shí),反射振幅方位變化幾乎獨(dú)立于背景孔隙度,對(duì)于含水儲(chǔ)層(圖14b),這種獨(dú)立性成立的孔隙度范圍大約變?yōu)?%～15%。這指導(dǎo)我們?cè)诎l(fā)育垂直裂縫體的含氣頁(yè)巖或含油水砂巖中,可利用測(cè)井等探測(cè)方法估計(jì)背景孔隙度,然后將估計(jì)值作為先驗(yàn)常數(shù)直接代入反演流程中來(lái)預(yù)測(cè)裂縫密度等參數(shù),以避免背景孔隙度地震難預(yù)測(cè)[110]的問(wèn)題。分析圖15,可以總結(jié)出含氣裂縫儲(chǔ)層反射振幅方位變化對(duì)裂縫橫縱比的依賴度比含水裂縫儲(chǔ)層低(即流體模量的降低將減弱反射振幅方位變化對(duì)裂縫縱橫比的依賴度);背景孔隙度較大的裂縫儲(chǔ)層反射振幅方位變化對(duì)裂縫橫縱比的依賴度比背景孔隙度較小的裂縫儲(chǔ)層低(即背景孔隙度的增大將減弱反射振幅方位變化對(duì)裂縫縱橫比的依賴度)。所以對(duì)于孔隙、喉道發(fā)育良好的含氣儲(chǔ)層,我們可以將裂縫橫縱比設(shè)置為常數(shù),直接利用反射振幅方位變化預(yù)測(cè)裂縫密度等參數(shù)。

圖14 反射振幅方位差異隨觀測(cè)方位φ和孔隙度變化曲線

圖15 反射振幅方位差異隨觀測(cè)方位φ和裂隙縱橫比變化曲線

綜合上述分析,我們可以總結(jié):考慮孔-裂隙間的流體流動(dòng)時(shí),裂縫走向始終對(duì)應(yīng)方位反射振幅的極大值點(diǎn),這為裂縫方位預(yù)測(cè)提供了理論基礎(chǔ);但孔、縫孤立時(shí),裂縫走向與方位反射振幅極值點(diǎn)(極大值點(diǎn)或極小值點(diǎn),這兩個(gè)點(diǎn)的觀測(cè)方位相差90°)的對(duì)應(yīng)關(guān)系受制于充填流體類(lèi)型等因素,直接利用反射振幅方位變化規(guī)律預(yù)測(cè)裂縫方位存在不確定性,此時(shí)需要其它預(yù)測(cè)技術(shù)為地震預(yù)測(cè)裂縫方位提供先驗(yàn)信息。另外,裂縫密度越大,地震振幅方位變化越劇烈;飽含的流體模量越小,裂縫走向反射振幅與裂縫法向反射振幅的差值越大。這些規(guī)律是地球物理學(xué)家利用地震方法預(yù)測(cè)裂縫密度和流體類(lèi)型的理論基礎(chǔ)之一。但地震反射振幅的方位變化規(guī)律是裂縫密度、孔裂隙中充填的流體類(lèi)型、喉道發(fā)育情況,甚至巖石的礦物成分等因素的綜合響應(yīng),難以直接利用反射振幅方位變化實(shí)現(xiàn)裂縫和流體的精準(zhǔn)識(shí)別,所以在分析清楚喉道和等徑孔隙發(fā)育情況的前提下,需要研究如何從反射振幅方位變化信息中解耦基質(zhì)礦物、裂縫及油氣等信息,提出一套“固-液-裂縫”解耦的裂縫儲(chǔ)層地震描述及油氣檢測(cè)的方法。除此,流體流動(dòng)對(duì)地震響應(yīng)影響較大,所以我們應(yīng)該在分析儲(chǔ)層孔喉結(jié)構(gòu)和滲透率的基礎(chǔ)上,準(zhǔn)確選擇巖石物理模型,防止巖石物理模型的誤用帶來(lái)解釋誤差。最后,給出在地震勘探實(shí)踐中或許有用的兩條結(jié)論:對(duì)于垂直裂縫發(fā)育的含氣頁(yè)巖或含油水砂巖地層,可將非地震方法的背景孔隙度估計(jì)值(甚至常數(shù)值)直接代入反演流程中來(lái)預(yù)測(cè)裂縫密度等參數(shù);對(duì)于孔喉發(fā)育良好的含氣儲(chǔ)層,可將裂縫縱橫比設(shè)置為常數(shù),直接利用反射振幅方位變化預(yù)測(cè)儲(chǔ)層參數(shù)。這樣避免了地震預(yù)測(cè)背景孔隙度和裂縫縱橫比難度大的問(wèn)題！

3 挑戰(zhàn)和機(jī)遇

等效各向異性巖石物理理論的快速發(fā)展有效推動(dòng)了裂縫型油氣儲(chǔ)層的勘探發(fā)現(xiàn),但隨著勘探目標(biāo)復(fù)雜程度的不斷提高,以及油氣儲(chǔ)層的開(kāi)發(fā)生產(chǎn)活動(dòng)對(duì)地震預(yù)測(cè)精度的要求不斷提高,裂縫儲(chǔ)層等效各向異性巖石物理理論的發(fā)展也迎來(lái)了新的挑戰(zhàn)。結(jié)合上文的分析和討論,在未來(lái)的研究中,我們應(yīng)聚焦于以下幾方面:

1) 裂縫儲(chǔ)層等效介質(zhì)巖石物理理論目前主要考慮了裂縫夾雜物和各向同性背景之間的彈性互動(dòng)作用且研究成果得到了應(yīng)用。但對(duì)于頁(yè)巖等裂縫型儲(chǔ)層,強(qiáng)各向異性背景與裂縫之間的彈性互動(dòng)鮮有研究;

2) 實(shí)際地層中的裂縫形狀、大小、表面粗糙度等特征具有較高的復(fù)雜性,應(yīng)該探索如何克服當(dāng)前裂縫儲(chǔ)層等效介質(zhì)巖石物理理論的理想假設(shè)帶來(lái)的誤差,建立一套能夠更精確地描述實(shí)際裂縫特征的模型;

3) 多尺度或跨頻帶的裂縫介質(zhì)地震巖石物理理論逐步發(fā)展起來(lái),例如Chapman巖石物理理論,已經(jīng)引起地球物理工作者的研究興趣,但地震波衰減特征或與地震波頻率有關(guān)的地震屬性和復(fù)雜儲(chǔ)層裂縫參數(shù)的量化關(guān)系有待深入研究和應(yīng)用。