張繁昌,李傳輝,吳國忱,王保麗
(中國石油大學(xué) 地球資源與信息學(xué)院,山東東營 257061)
隨著油氣勘探重點(diǎn)由構(gòu)造型油氣藏向巖性油氣藏、隱蔽油氣藏的轉(zhuǎn)移,地層對地震波的吸收衰減越來越引起地球物理工作者的重視。地震波在地下介質(zhì)中傳播時(shí)的能量衰減分為二類:①與地震波傳播特性有關(guān)的衰減,包括球面擴(kuò)散,散射及透射損失等,這部份能量在資料處理時(shí)已對其進(jìn)行補(bǔ)償;②地層本身的本征衰減,即地層的吸收,這種介質(zhì)本身所固有的吸收特性通常用固有品質(zhì)因子Q來描述[1]。Q值與介質(zhì)的結(jié)構(gòu)、孔隙度,以及孔隙流體的性質(zhì)密切相關(guān),被認(rèn)為是隱蔽油氣藏預(yù)測最為敏感、可靠的依據(jù)之一[2]。
目前,地震資料的解釋已進(jìn)入以巖性解釋為主的階段,而巖石的彈性參數(shù)在巖性解釋中,扮演著非常重要的角色。當(dāng)?shù)貙游諊?yán)重時(shí),就需要考慮吸收對巖石彈性參數(shù)帶來的影響[3]。例如對地層波阻抗、彈性阻抗的影響,以及由于吸收引起的AVO類型的改變等,這些影響將會(huì)直接關(guān)系到巖性解釋的準(zhǔn)確性,因而研究地層吸收對巖石彈性參數(shù)的影響具有十分重要意義。
自Futter man把吸收衰減作為地層的基本特性以來,許多學(xué)者在地層吸收方面作了大量的研究:Tonn[4]曾對十多種地層吸收和品質(zhì)因子計(jì)算方法作了比較;尹陳等人[5]利用二維粘彈性波動(dòng)方程研究了地震波的吸收衰減特性;Koesoemadinata等[6]利用實(shí)驗(yàn)室測量數(shù)據(jù),建立了砂巖儲(chǔ)層物性參數(shù)(孔隙度、含水飽和度、泥質(zhì)含量、有效地層壓力等)與地球物理參數(shù)(縱波速度Vp、橫波速度Vs、密度ρ、品質(zhì)因子Q等)的經(jīng)驗(yàn)關(guān)系,并利用這些經(jīng)驗(yàn)關(guān)系由儲(chǔ)層物性參數(shù)來估算Q值。因此,綜合利用地層的品質(zhì)因子和縱橫波信息進(jìn)行儲(chǔ)層預(yù)測,可以進(jìn)一步降低勘探風(fēng)險(xiǎn)。
作者在本文中從粘彈性介質(zhì)的本構(gòu)方程出發(fā),推導(dǎo)了地層吸收對波阻抗值的影響,進(jìn)而分析并討論了地層吸收對地震數(shù)據(jù)的AVO特征,以及彈性阻抗的影響,最后將地層品質(zhì)因子Q融入到疊前彈性阻抗反演過程中,實(shí)現(xiàn)了地層品質(zhì)因子與縱橫波阻抗參數(shù)的同時(shí)反演。
地下介質(zhì)對地震波的吸收通常用品質(zhì)因子Q來衡量,品質(zhì)因子越小,地層對地震波的吸收越嚴(yán)重。為了量化地層吸收對巖石彈性參數(shù)的影響,下面從粘彈性介質(zhì)本構(gòu)方程出發(fā),來推導(dǎo)品質(zhì)因子與波阻抗之間的關(guān)系。
在線性粘彈性介質(zhì)的基本模型中,Kelvin-Voigt模型(見圖1(a))未考慮應(yīng)力作用下應(yīng)變的突然變化,也不能表示應(yīng)力消失后的剩余應(yīng)變;而Maxwell模型(見圖1(b))不具備蠕變特征;標(biāo)準(zhǔn)線性粘彈性模型(見圖1(c))彌補(bǔ)了前二者的不足[7],作者采用該粘彈性模型來討論吸收特性。
圖1 線性粘彈性介質(zhì)基本模型Fig.1 Basic linear viscoelastic model
標(biāo)準(zhǔn)線性粘彈性模型的本構(gòu)方程為[8]:
其中 E1和E2是彈性模量;η為粘性系數(shù);σ為應(yīng)力;ε為應(yīng)變。
應(yīng)力與應(yīng)變的關(guān)系為σ=M(ω)ε,其中關(guān)于頻率ω的復(fù)模量M(ω)為:
根據(jù)復(fù)模量,地層品質(zhì)因子Q即可表示為:
由式(3)得到最大衰減為:其中 M0、M∞分別是高頻、低頻極限下的模量。
為了使品質(zhì)因子Q與模量之間的關(guān)系更明顯,對式(4)進(jìn)行變換得:
考慮到Q的取值范圍,忽略Qmax-1的平方項(xiàng),則有:
從式(6)看出,地層吸收作用的確對巖石彈性參數(shù)產(chǎn)生了影響,這種影響可能在大多數(shù)情況下是很小的,但是如果介質(zhì)含烴,使介質(zhì)吸收嚴(yán)重,那這種影響將比較大。
當(dāng)不考慮地下介質(zhì)吸收的情況下,隨入射角θ變化的P波反射系數(shù)RPP,可以用Aki-Richard方程計(jì)算[9]:
其中 R(ˉθ)表示隨角度變化的PP波反射系數(shù);ˉα、ˉβ、ˉρ、ˉγ和ˉθ分別表示平均P波速度、平均S波速度、平均密度、ˉβ/ˉα比值,以及分界面的入射角和透射角的平均角度。類似的,Δ α、Δ β、Δ ρ是界面二側(cè)P波速度、S波速度,以及密度的變化量。
由于在小角度情況下,sin2θ≈tan2θ,并取S波和P波的背景速度比值=0.5,則由方程(7)就可以得到式(8)表示的W iggins方程[10]:
其中 RP和RS分別表示法向入射時(shí)縱波和橫波的反射系數(shù)。
利用式(6)對式(8)進(jìn)行改寫得到:
把式(8)、式(9)二式相減,得:
由于 RP=(IP2-IP1)/(IP2+IP1);
于是有:
由式(11)可知,在碎屑巖油氣藏中,如果砂巖中含有了油氣,那么其Q值比較小(一般小于50),而泥巖的Q值比較大(一般大于100)。所以式(11)中/2的值一般為正,也即AVO的截距和梯度由于吸收的影響,都以相同的衰減量減小。
下面通過四類典型的砂泥巖AVO模型來進(jìn)行分析,模型數(shù)據(jù)見表1。
如圖2所示,圖2中深色點(diǎn)線表示不考慮吸收時(shí)的AVO曲線。從圖2(a)、圖2(c)看出,雖然泥巖的Q值從100變化到150,但是它們對應(yīng)的AVO曲線幾乎重合;而從圖2(b)、圖2(d)卻明顯看出,下層砂巖的Q值從50變化到5,AVO曲線的形狀發(fā)生了顯著變化,截距與梯度都隨著Q值的減小而減小,當(dāng)Q值取到5時(shí),第一類AVO甚至轉(zhuǎn)化成了第二類AVO,而第二類AVO轉(zhuǎn)化成了第三類AVO,這說明了砂巖的吸收對AVO的影響比較明顯。
從圖3(見下頁)可以看出,地層吸收對第三類和第四類AVO的影響,與前二類基本相同。曲線的截距和梯度都隨著Q值的減小而減小,且下層砂巖的吸收,對AVO影響比上層泥巖要大。當(dāng)砂巖吸收嚴(yán)重時(shí),AVO的類型甚至發(fā)生了變化。
通過上面對四類AVO的分析可以得出,上層泥巖的吸收對AVO的影響非常小,而下層砂巖的吸收對其影響比較大。如果砂巖含氣,Q值很低時(shí),AVO的類型甚至?xí)l(fā)生變化:第一類AVO會(huì)轉(zhuǎn)化為第二類,第二類和第四類AVO會(huì)轉(zhuǎn)化為第三類,而第三類AVO保持不變。因此,當(dāng)?shù)卣鸩ㄔ诘叵聜鞑ミ^程中遇到儲(chǔ)層時(shí),地層的吸收使地震波能量衰減,使其AVO響應(yīng)的截距和梯度都減小。如果忽略這種影響,有可能使我們在利用AVO識(shí)別儲(chǔ)層時(shí),誤把油層認(rèn)為是氣層,把水層誤認(rèn)為油層,從而造成錯(cuò)誤的判斷。
表1 四類AVO典型界面模型Tab.1 Stratigraphic models for four classes ofAVO
圖2 地層吸收對第一類AVO和第二類AVO的影響Fig.2 The effect of absorption on class 1 and class 2 ofAVO
圖3 地層吸收對第三類AVO和第四類AVO的影響Fig.3 The effect of absorption on class 3 and class 4 of AVO
彈性阻抗是聲阻抗的推廣,它隨入射角的變化而變化,聲阻抗是入射角為零時(shí)彈性阻抗的一個(gè)特例[11]。常規(guī)彈性阻抗方程是基于Aki-Richard方程推導(dǎo)出來的,作者在本文將基于地層吸收情況下的W iggins方程來推導(dǎo)彈性阻抗方程。由式(9)得:
對其積分,最后可以得到彈性阻抗方程式(12)。
表2 用于驗(yàn)證彈性阻抗方程的模型數(shù)據(jù)Tab.2 Model data to test the EI equation
與常規(guī)的彈性阻抗公式類似,EIQ表達(dá)式也存在著彈性阻抗值隨角度變化在量綱上有很大變化的問題,這不利于不同角度彈性阻抗值的對比。為了消除入射角變化對量綱的影響,可引入常量IP0和IS0對式(12)進(jìn)行規(guī)范化,得到式(13):
從式(13)容易看出,地層吸收對彈性阻抗也有影響,且這種影響隨著角度的不同而不同。下面,以Marmousi2模型為例來討論地層吸收對彈性阻抗的影響。在圖5所示的Mar mousi2模型中,第2 200道經(jīng)過二套儲(chǔ)層(圖5中橢圓標(biāo)注區(qū))。從圖6中看出,不同角度的彈性阻抗值,都因地層吸收而減小,尤其是在儲(chǔ)層處更明顯,且隨著角度的增加,彈性阻抗值的衰減量也增加。充分利用這些性質(zhì),將有利于彈性阻抗在識(shí)別儲(chǔ)層方面的應(yīng)用。
圖4 反射系數(shù)對比圖Fig.4 Comparison of reflectivities
圖5 Marmousi2模型Fig.5 Marmousi2 model
彈性阻抗反演技術(shù)充分利用了疊前數(shù)據(jù)豐富的振幅和旅行時(shí)信息,利用了不同入射角的部份疊加數(shù)據(jù)體得到彈性阻抗數(shù)據(jù)體。利用彈性阻抗參數(shù),可以分析地層在近、中、遠(yuǎn)偏移距上的變化規(guī)律,也可以從其中提取出地層的縱橫波速度、密度、泊松比等巖性參數(shù)。用式(13)進(jìn)行彈性阻抗反演與常規(guī)彈性阻抗反演類似,需要地震資料處理,測井資料標(biāo)定,角度子波提取和彈性阻抗反演這幾個(gè)流程[12]。作者在本文中利用不同角度的疊前彈性阻抗結(jié)果,從EIQ數(shù)據(jù)體中同時(shí)提取出地層品質(zhì)因子Q和縱橫波阻抗,來更好地指導(dǎo)儲(chǔ)層的預(yù)測。
對式(13)二邊取對數(shù),可得:
圖6 有吸收和無吸收情況下彈性阻抗曲線的對比Fig.6 Comparison of EI logs under the situation of absorption and non-absorption
為了得到縱橫波阻抗和品質(zhì)因子,需要三個(gè)不同角度的EIQ數(shù)據(jù)體。將三個(gè)角度分別代入上式,并用矩陣形式表示為式(14)。
方程(14)的求解形如對AX=B的求解,由于角度已知,只需知道三個(gè)相互獨(dú)立的EIQ數(shù)據(jù)體,用最小平方法便可以求得縱橫波阻抗和品質(zhì)因子。
但在實(shí)際應(yīng)用中,直接利用方程(14)求解所得到的ln(IP/IP0)、ln(IS/IS0)和ln(1-Q-1P)值極不穩(wěn)定,并且有些值不符合實(shí)際的物理和地質(zhì)意義。因此,在實(shí)際應(yīng)用時(shí),是利用井資料先擬合方程中的系數(shù)。
可以證明,方程(14)可以寫成下面的形式:
其中 a′1、b′1、c′1;a′2、b′2、c′2;a′3、b′3、c′3為待定系數(shù)。
以ln(IP/IP0)為例,對于同一道的不同采樣點(diǎn)(即不同時(shí)間t)有:
由于在計(jì)算系數(shù)a′1、b′1、c′1時(shí),所用的EIQ數(shù)據(jù)體都來自于反演所得結(jié)果,因此,這些系數(shù)可選取反演得到的井旁道EIQ數(shù)據(jù)和測井曲線,通過最小二乘擬合得到。這種巖性參數(shù)與EIQ之間的關(guān)系最密切,得到的a′1、b′1、c′1也最具有代表性。
同理,ln(IS/IS0)和ln(1-Q-1P)對應(yīng)的系數(shù),也可以由該方法獲得,最后得到九個(gè)常系數(shù)a′1、b′1、c′1;a′2、b′2、c′2;a′3、b′3、c′3之后,將它們代入式(15),就可得到方程(20):
將反演所得的各角度EIQ數(shù)據(jù)體帶入方程(20),即可獲得各道任意采樣點(diǎn)處的縱橫波阻抗和地層品質(zhì)因子。
以Marmousi2模型為例,對上述方法進(jìn)行測試。首先利用圖7(見下頁)所示的Marmousi2模型參數(shù),由式(13)求得三個(gè)角度的EIQ數(shù)據(jù)(角度值分別為25°、15°和5°),然后利用方程(20)從中反演出縱橫波阻抗和品質(zhì)因子;圖8(見下頁)為EIQ數(shù)據(jù)中沒有噪音的反演結(jié)果,可以看出,在無噪音時(shí),該方法可以完全恢復(fù)出Mar mousi2模型的縱橫波阻抗和品質(zhì)因子值;圖9(見下頁)為加入5%噪音時(shí)反演的結(jié)果,反演結(jié)果與真實(shí)模型數(shù)據(jù)吻合很好;圖10(見下頁)為加入10%噪音時(shí)的反演結(jié)果,反演出的縱橫波阻抗和真實(shí)模型吻合很好,品質(zhì)因子稍有抖動(dòng),但也有很好的相似性;圖11(見下頁)為從模型中抽出一道的顯示結(jié)果。經(jīng)過測試表明,該方法是一種有效的求取縱橫波阻抗和品質(zhì)因子的方法。
圖7 Mar mousi2模型參數(shù)Fig.7 Marmousi2 model parameters
圖8 無噪時(shí)的反演剖面Fig.8 Inversion section with noise free seismology
圖9 5%噪音時(shí)的反演剖面Fig.9 Inversion section with seismology containing 5%noise
圖10 10%噪音時(shí)的反演剖面Fig.10 Inversion section with seismology containing 10%noise
圖11 反演結(jié)果曲線對比Fig.11 Comparison of inverted logs and model data
在研究區(qū)的中層氣藏具備良好的巖石物理?xiàng)l件,測井資料顯示,氣層縱波速度較低,地震反射特征表現(xiàn)為強(qiáng)振幅、低頻率,具有亮點(diǎn)型反射特征。目的層處的砂巖儲(chǔ)層含氣后縱波阻抗明顯降低,由于孔隙流體類型對橫波的影響很小,所以目的層段的橫波阻抗與圍巖相當(dāng)。
利用QAVO方法,從本區(qū)的疊前部份角度疊加地震數(shù)據(jù)反演得到相應(yīng)的數(shù)據(jù)體,并從中提取出縱橫波阻抗和品質(zhì)因子。圖12(見下頁)是在1.35 s處,穿過含氣儲(chǔ)層位置的縱橫波阻抗和品質(zhì)因子時(shí)間切片,可以看出,縱波阻抗在含氣目標(biāo)區(qū)表現(xiàn)為低值,切片上可明顯看出含氣砂體的大致范圍,而橫波阻抗切片上則看不出儲(chǔ)層范圍。品質(zhì)因子在含氣目標(biāo)區(qū)也表現(xiàn)為明顯低值,從而可以比較準(zhǔn)確地推斷含氣砂體的空間展布。該結(jié)果與驗(yàn)證井Y74的結(jié)論一致,表明了含氣儲(chǔ)層預(yù)測結(jié)果與實(shí)際情況吻合。由此可見,利用基于粘彈性理論的地層品質(zhì)因子,與縱橫波阻抗進(jìn)行QAVO同步反演是可行的。綜合利用品質(zhì)因子和縱橫波信息,可以有效識(shí)別儲(chǔ)層,符合儲(chǔ)層含氣后縱波阻抗降低,橫波阻抗基本不變且吸收嚴(yán)重的規(guī)律。
(1)從粘彈性介質(zhì)的本構(gòu)方程出發(fā),作者在文中討論了地層吸收對波阻抗等巖石彈性參數(shù)的影響,給出了吸收前、后波阻抗之間的定量關(guān)系,證明了地層的吸收作用使波阻抗減小,特別是在含氣儲(chǔ)層段。
(2)作者在文中討論了地層吸收對地震AVO響應(yīng)的影響,得出由于吸收使得P波反射系數(shù)的截距和梯度都發(fā)生減小的結(jié)論。并且通過四類AVO模型驗(yàn)證發(fā)現(xiàn),當(dāng)吸收嚴(yán)重時(shí),AVO類型甚至發(fā)生變化,會(huì)造成AVO假象。并進(jìn)一步從考慮吸收情況下的W iggins方程出發(fā),推導(dǎo)出了帶品質(zhì)因子的彈性阻抗方程。
(3)利用疊前彈性阻抗反演框架,實(shí)現(xiàn)了地層品質(zhì)因子與縱橫波阻抗的同步QAVO反演,并將地層的吸收特征與縱橫波阻抗等彈性參數(shù)相結(jié)合,進(jìn)行綜合儲(chǔ)層表征,提高了儲(chǔ)層預(yù)測和烴類檢測的可靠性。
圖12 1.35 s處的各參數(shù)切片F(xiàn)ig.12 Inverted parameter slices at 1.35 s
[1] FUTTERMAN W I.Dispersive body waves[J].Journal of Geophysical Research,1962,67(12):5279.
[2] 張璐,白曉寅.單程波法地震波衰減特征數(shù)值模擬分析[J].物探化探計(jì)算技術(shù),2009,31(6):541.
[3] UPENDRA K T,MCMECHAN G A.Effects of incomplete parameterization on full-wavefield viscoelastic seismic data for petrophysical reservoir properties[J].Geophysics,2007,72(3):O9.
[4] TONN R.The determination of the seismic quality factor Q from VSP data:A comparison of different computationalmethods[J].Geophysical Prospecting,1991,39(1):1.
[5] 尹陳,賀振華,黃德濟(jì).基于二維粘滯性波動(dòng)方程的地震波衰減理論研究[J].物探化探計(jì)算技術(shù),2008,30(5):353.
[6] KOESOEMAD INATA A P,MCMECHAN GA.Petroseismic inversion for sandstone properties[J].Geophysics,2003,68(5):1611.
[7] JOSEM C.Wave fields in real media:Wave Propagation In Anisotropic,Anelastic And Porous Media[M].Seismic Exploration,2001.
[8] MAVKO G,MUKERJI T,DVORKIN J.Rock physics handbook[M].Cambridge University Press,1998.
[9] AKI K,R ICHARDS P G.Quantitative Seismology:Theory and Methods[J].W.H.Freeman and Co,1980:123.
[10]WIGG INS R,KENNY G S,MCCLURE C.A method for determining and displaying the shear-velocity reflectivities of a geologic formation[M].European patent Application,1983.
[11]CONNOLLY P.Elastic impedance[J].The Leading Edge,1999,18(4):438.
[12]王保麗,印興耀,張繁昌.彈性阻抗反演及應(yīng)用研究[J].地球物理學(xué)進(jìn)展,2005,20(1):89.