張世鑫,杜向東,韓文明,孫林潔

(中海油研究總院,北京100027)

?

適于Russell流體因子提取的孔隙彈性介質(zhì)反射系數(shù)近似方程

張世鑫,杜向東,韓文明,孫林潔

(中海油研究總院,北京100027)

摘要:Russell流體因子是表征儲(chǔ)層孔隙流體特征的重要敏感流體因子,其提取結(jié)果的可靠性是決定流體判識(shí)質(zhì)量的核心問題之一。疊前地震反演是提取流體因子的主要途徑,而反射系數(shù)近似方程是疊前反演的基礎(chǔ)。結(jié)合Biot-Gassmann多孔彈性介質(zhì)理論對(duì)反射系數(shù)近似方程進(jìn)行了重新推導(dǎo),得到了包含Russell流體因子的孔隙彈性介質(zhì)反射系數(shù)線性近似方程;通過分析方程的精度和反演敏感特征,驗(yàn)證了新反射系數(shù)近似方程的可靠性和Russell流體因子項(xiàng)的高反演敏感度;最終基于貝葉斯反演框架的AVO反演實(shí)現(xiàn)了Russell流體因子的定量提取。模型試算和實(shí)際應(yīng)用結(jié)果表明,基于新反射系數(shù)近似方程的疊前反演方法能夠穩(wěn)定可靠地從地震資料中提取Russell流體因子,為儲(chǔ)層流體識(shí)別提供了一種新的流體因子計(jì)算方法。

關(guān)鍵詞:Russell流體因子;反射系數(shù)近似方程;疊前地震反演;流體識(shí)別

地震波所蘊(yùn)含的豐富的動(dòng)力學(xué)與運(yùn)動(dòng)學(xué)信息有助于揭示地下巖層及其孔隙流體特征。隨著物探技術(shù)的進(jìn)步,以地震資料為主體的地下巖層孔隙流體識(shí)別技術(shù)已成為現(xiàn)階段儲(chǔ)層精細(xì)描述的關(guān)鍵技術(shù)之一。在保證地震資料品質(zhì)的前提下,儲(chǔ)層流體識(shí)別精度主要取決于所選取的流體因子是否敏感及其計(jì)算方式是否可靠[1-5]。

現(xiàn)今,單純依靠振幅異常的“亮點(diǎn)”等定性檢測(cè)技術(shù)已無法滿足復(fù)雜巖性油氣藏的流體識(shí)別需求,SMITH等[6]率先提出了利用縱、橫波速度相對(duì)變化率的加權(quán)疊加檢測(cè)烴類異常的方法,并首次提出了流體因子的概念;隨后,GIDLOW等[7]和FATTI等[8]利用縱、橫波阻抗相對(duì)變化率構(gòu)建了新流體因子;GRAY等[9]則使用拉梅參數(shù)的相對(duì)變化率為流體因子進(jìn)行流體判識(shí)。流體識(shí)別的研究重點(diǎn)也逐漸從定性檢測(cè)過渡到利用流體因子實(shí)現(xiàn)儲(chǔ)層含流體異常的定量判識(shí),并且隨著地震反演技術(shù)的發(fā)展,流體因子類型也逐步從反射系數(shù)域發(fā)展到阻抗域。GOODWAY等[10]提出了拉梅常數(shù)-密度(Lambda-Mu-Rho,簡(jiǎn)稱LMR)技術(shù),利用儲(chǔ)層的拉伸剪切特性區(qū)分儲(chǔ)層含流體類型;QUAKENBUSH等[11]提出了泊松阻抗的概念;RUSSELL等[12-13]基于多孔彈性介質(zhì)巖石物理理論,提出了可以反映孔隙流體彈性效應(yīng)的Russell流體因子和Gassmann流體/孔隙項(xiàng),通過研究干巖石縱橫波速度比的取值,認(rèn)為LMR技術(shù)只是Russell流體因子的特例,并且與其它利用權(quán)差運(yùn)算得到的流體因子相比,Russell流體因子和Gassmann流體/孔隙項(xiàng)具有更明確的物理意義和更敏感的流體判識(shí)能力。目前,二者在碎屑巖儲(chǔ)層流體判識(shí)中占據(jù)重要位置,且研究發(fā)現(xiàn),雖然兩者對(duì)碎屑巖儲(chǔ)層流體判識(shí)的敏感度基本一致,但是Russell流體因子綜合了Gassmann流體/孔隙項(xiàng)和密度項(xiàng)參數(shù),對(duì)一些固結(jié)成熟的碎屑巖儲(chǔ)層,表現(xiàn)為更高的流體判識(shí)敏感性[1,14]。

疊前地震反演作為提取儲(chǔ)層彈性參數(shù)的主要途徑,其反演結(jié)果的可靠性直接影響流體因子的可靠性。在疊前反演中,反射系數(shù)近似方程是反演框架構(gòu)建的基礎(chǔ),如包含縱、橫波速度和密度的Aki-Richards近似等[15]。由于反射系數(shù)近似方程中的彈性參數(shù)信息權(quán)重存在較大差異,即反射系數(shù)近似方程中的彈性參數(shù)具有不同的反演敏感性,且受地震反演本身的“病態(tài)解”影響,疊前反演得到的彈性參數(shù)不可避免地存在多解性[16],而通過間接代數(shù)運(yùn)算得到的流體因子則會(huì)進(jìn)一步放大累計(jì)誤差,影響流體識(shí)別的可靠性[17]。為了提高流體因子的估算精度,學(xué)者們進(jìn)行了一系列研究,RUSSELL等[13]基于多孔彈性介質(zhì)巖石物理理論首次建立了包含Gassmann流體/孔隙項(xiàng)、剪切模量和密度項(xiàng)的反射系數(shù)近似方程,提出了從疊前道集資料中直接估算Gassmann流體/孔隙項(xiàng),并以此作為流體因子直接進(jìn)行儲(chǔ)層流體識(shí)別的研究思路。隨后,敏感流體因子的直接反演逐步成為學(xué)者們研究的熱點(diǎn),印興耀等[18]利用彈性阻抗反演實(shí)現(xiàn)了Gassmann流體/孔隙項(xiàng)的估算;楊培杰等[19]基于Russell反射系數(shù)近似公式提出了可變點(diǎn)約束AVO疊前反演算法,提高了Gassmann流體/孔隙項(xiàng)的估算精度。考慮到流體因子的敏感性具有區(qū)域適用性,圍繞Russell流體因子的直接反演,學(xué)者們也進(jìn)行了一些研究,印興耀等[1]基于包含Russell流體因子的兩項(xiàng)彈性阻抗方程實(shí)現(xiàn)了Russell流體因子的直接提取,取得了較好的深層流體識(shí)別應(yīng)用效果,同時(shí)也指出在大角度疊前地震資料可用的情況下,要盡可能借助大角度地震資料來進(jìn)一步提高Russell流體因子計(jì)算的可靠性。本文為充分利用疊前地震資料包含的有效信息,在多孔彈性介質(zhì)巖石物理理論指導(dǎo)下,推導(dǎo)了包含Russell流體因子的3項(xiàng)反射系數(shù)近似方程,在貝葉斯反演理論框架下利用疊前AVO反演實(shí)現(xiàn)了Russell流體因子的直接反演,為流體因子的可靠計(jì)算提供了一種新的方法。

1新反射系數(shù)近似方程

在地震波引起孔隙流體壓力變化符合Biot-Gassmann多孔彈性介質(zhì)理論應(yīng)用條件的前提下[20-21],RUSSELL等[12]利用孔隙彈性介質(zhì)的力學(xué)特征對(duì)含流體多孔介質(zhì)的縱波速度與橫波速度進(jìn)行了改寫,并進(jìn)一步推導(dǎo)得到了Russell流體因子F,即:

(1)

式中:f為Gassmann流體/孔隙項(xiàng);ρ為多孔流體飽和彈性介質(zhì)的密度;vP為縱波速度;vS為橫波速度;γdry為干巖石縱橫波速度比。

為得到包含Russell流體因子的反射系數(shù)近似方程,對(duì)Aki-Richards反射系數(shù)近似方程進(jìn)行重新改寫。Aki-Richards反射系數(shù)近似方程為:

(2)

對(duì)(2)式兩邊同時(shí)乘以(vPρ)2,得到:

(3)

對(duì)(1)式兩邊微分,則有:

(4)

整理后可以得到:

(5)

將(5)式代入(3)式得:

(6)

又因?yàn)榱黧w因子可以表示為:

(7)

兩邊微分并變形可以得到:

(8)

將(8)式代入(6)式,并進(jìn)行變換可以得到:

(9)

又因?yàn)?

(10)

(11)

(9)式最終化為:

(12)

(12)式為包含Russell流體因子F,剪切模量μ和密度ρ的孔隙彈性介質(zhì)反射系數(shù)近似方程。與常規(guī)反射系數(shù)近似方程不同的是該方程考慮了孔隙流體的彈性效應(yīng),不僅更直觀地表征了地震波中蘊(yùn)含的流體信息,而且有利于流體因子的定量反演。

2近似精度與反演敏感性分析

2.1近似精度分析

利用RUSSELL等[13]提出的砂巖模型來驗(yàn)證公式(12)的精度。砂巖模型為兩種雙層砂巖模型,模型參數(shù)如表1所示,模型一為上層砂巖含氣(含氣飽和度為90%,含水飽和度為10%),下層砂巖含水,兩層孔隙度均為25%;模型二的兩層砂巖均含水,但孔隙度存在差異,上層孔隙度為25%,下層孔隙度為15%。分別采用精確Zoeppritz方程、Aki-Richard近似方程和本文新推導(dǎo)的反射系數(shù)近似方程計(jì)算這兩種模型中砂巖界面的反射系數(shù),并比較了兩種近似方程與精確Zoeppritz方程之間的相對(duì)誤差,其結(jié)果如圖1和圖2所示。

表1　兩種砂巖模型參數(shù)[13]

圖1　模型一的反射系數(shù)計(jì)算結(jié)果a 反射系數(shù); b 反射系數(shù)相對(duì)誤差

從圖1a可以看到,當(dāng)界面兩側(cè)介質(zhì)孔隙度沒有變化時(shí),(12)式與Zoeppritz精確方程的結(jié)果匹配較好;從圖1b可以看到,當(dāng)入射角為0時(shí),公式(12)與Zoeppritz方程的相對(duì)誤差為-2.5%,且隨著入射角的增加,兩者之間的相對(duì)誤差逐漸減小。從圖2a則可以看到,當(dāng)界面兩側(cè)砂巖孔隙度不同時(shí),公式(12)在小角度入射時(shí),與精確解相差不大,隨著入射角的增大,公式(12)的計(jì)算數(shù)值逐漸偏離精確Zoeppritz結(jié)果;由圖2b可知,公式(12)與Zoeppritz精確方程的最大誤差在入射角為40°時(shí)出現(xiàn),相對(duì)誤差為-8%,但是公式(12)的計(jì)算值與Aki-Richards近似的結(jié)果基本一致,這是因?yàn)樾路瓷湎禂?shù)近似方程本身是從Aki-Richards近似推導(dǎo)出來的,兩者精度基本一致,誤差僅在于干巖樣縱橫波速度比值的選取;考慮到反射系數(shù)近似是建立在界面兩側(cè)彈性差異較小的前提下,而孔隙度的變化對(duì)巖石模量的影響要大于孔隙流體變化造成的影響,因此界面兩側(cè)孔隙度的變化會(huì)產(chǎn)生更大的大角度近似誤差,但從總體上說,公式(12)在小角度入射情況下仍然可以滿足反射系數(shù)的近似精度?？紤]到實(shí)際應(yīng)用中的角度道集一般不超過35°,即在誤差允許的范圍內(nèi)使用新的反射系數(shù)近似方程進(jìn)行疊前反演是可行的。

圖2　模型二的反射系數(shù)計(jì)算結(jié)果a 反射系數(shù); b 反射系數(shù)相對(duì)誤差

2.2反演敏感性分析

由URSIN等[22],NICOLAO等[23]的研究得出,不同反射系數(shù)線性近似方程所表達(dá)的彈性參數(shù)信息是不同的,進(jìn)而影響了相應(yīng)彈性參數(shù)反演的可靠性。為檢驗(yàn)基于公式(12)的Russell流體因子的反演可行性,在此采用NICOLAO等[23]提出的信息敏感分析方法對(duì)公式(12)包含的彈性參數(shù)信息進(jìn)行定量分析。

將公式(12)以矩陣形式表示為:

(13)

NICOLAO等[23]指出,通過分析模型空間特征值及特征向量,可以定量分析反射系數(shù)近似方程中蘊(yùn)含的彈性參數(shù)的信息敏感性。其中,正演映射矩陣的大特征值及對(duì)應(yīng)的特征向量表示信息敏感的彈性參數(shù),而小特征值則表示信息不敏感、在反演中不易得到的彈性參數(shù)。

對(duì)正演映射矩陣G(θ)進(jìn)行奇異值分解:

(14)

式中:U表示數(shù)據(jù)空間的特征向量;Λ表示奇異值對(duì)角矩陣;V表示模型空間的特征向量。設(shè)γsat=2.000,γdry=1.732,最大入射角范圍為0～50°。圖3 給出了特征值隨最大入射角的變化曲線。其中,特征值的大小用分貝(20lgλ)表示。由圖3可知,所有特征值都隨最大入射角的增加而增大,這說明模型空間包含的彈性參數(shù)信息量隨著入射角的增加而增大,即入射角越大,反演的彈性參數(shù)越可靠。在最大入射角為30°時(shí),第一特征值幾乎是第二特征值所包含能量的100倍,是第三特征值所包含能量的10000倍?？紤]到第三特征值包含能量很少,且隨入射角的增加也很難確定待反演的參數(shù)組合,因此只分析第一、第二特征值對(duì)應(yīng)的特征向量即可確定易于反演的彈性參數(shù)。

圖3　模型空間中的3個(gè)特征值隨最大入射角的變化曲線

圖4a和圖4b分別給出了模型空間第一特征值和第二特征值對(duì)應(yīng)的特征向量的方向余弦隨最大入射角的變化曲線。當(dāng)方向余弦為0時(shí),無法反演該彈性參數(shù),且數(shù)值偏離0的程度越大,反演該彈性參數(shù)的難度越小。從圖4a可知,當(dāng)最大入射角較小時(shí),第一特征值對(duì)應(yīng)的Δμ/μ特征向量方向余弦偏離0的幅度最大,其次是ΔF/F,Δρ/ρ偏離0的幅度最小;隨著最大入射角的增加,ΔF/F和Δρ/ρ的特征向量方向余弦偏離0的幅度逐漸增大,而Δμ/μ則逐漸減小;當(dāng)最大入射角大于27°時(shí),ΔF/F的偏離幅度要高于Δμ/μ,即Russell流體因子信息隨角度的增加而增大,而Δρ/ρ的幅度則始終很小。圖4b中3個(gè)彈性參數(shù)的特征向量方向余弦偏離0的幅度相差不大,相比較而言,Δρ/ρ的偏離幅度最大,其次是ΔF/F和Δμ/μ。

由于第一特征值對(duì)應(yīng)的信息量最大,第二特征值對(duì)應(yīng)的信息只有在入射角很大的情況下才能生效。由圖4可知,基于公式(12)可較好地確定ΔF/F和Δμ/μ,而確定Δρ/ρ的難度相對(duì)較大,且為了保證ΔF/F的準(zhǔn)確性,需選用入射角大于30°的地震道集。因此,當(dāng)最大入射角大于30°時(shí),ΔF/F具有最高的反演敏感性,這與印興耀等[1]的認(rèn)識(shí)也是一致的。

圖4　模型空間中特征向量的方向余弦隨最大入射角的變化曲線a 第一特征值對(duì)應(yīng)的方向余弦; b 第二特征值對(duì)應(yīng)的方向余弦

3Russell流體因子反演

3.1反演方程

貝葉斯反演理論能夠利用必要的先驗(yàn)約束縮小解空間,且保證反演結(jié)果的穩(wěn)定性[24-26]。在貝葉斯反演框架下建立了基于新反射系數(shù)近似方程的反演方程,通過先驗(yàn)地質(zhì)約束最大似然的反演過程,實(shí)現(xiàn)了Russell流體因子的直接反演。

一般情況下,地震記錄中的噪聲信息服從正態(tài)分布,即似然函數(shù)是高斯分布。為了保護(hù)地下的弱反演信息,在此假設(shè)反映地下彈性參數(shù)相對(duì)變化量特征的先驗(yàn)信息服從柯西分布,利用貝葉斯公式將似然函數(shù)與先驗(yàn)函數(shù)結(jié)合起來,得到包含反演參數(shù)m的后驗(yàn)概率函數(shù)P(m|d):

(15)

通過求解后驗(yàn)概率最大化可以得到反演目標(biāo)函數(shù):

(16)

(17)

在給定初值的前提下,利用共軛梯度方法循環(huán)迭代求解公式(16)即可得到反問題的解。得到的結(jié)果是3個(gè)參數(shù)的相對(duì)變化率mF,mμ和mρ,再通過求解方程(18)可以最終得到Russell流體因子

F,剪切模量μ和密度ρ的數(shù)據(jù)體。

(18)

式中:t是時(shí)間采樣點(diǎn);t0是起始時(shí)間。

3.2模型測(cè)試及分析

為驗(yàn)證反演方法的可靠性,以實(shí)際測(cè)井資料為模型數(shù)據(jù)進(jìn)行測(cè)試。模型如圖5所示,用精確Zoeppritz方程得到不同采樣點(diǎn)的P波反射系數(shù),將其與40Hz零相位Ricker子波褶積生成疊前角度道集(角度范圍為0～45°,采樣間隔2ms),以此為基礎(chǔ)對(duì)不同信噪比條件下的反演結(jié)果進(jìn)行分析。圖6a和圖6b分別為未含噪聲和信噪比為2∶1的角度道集數(shù)據(jù)。針對(duì)這兩套數(shù)據(jù),分別采用直接反演方法和常規(guī)方法(先反演得到縱、橫波速度和密度,再通過間接代數(shù)運(yùn)算得到流體因子等參數(shù))反演得到Russell流體因子、剪切模量和密度,與模型數(shù)據(jù)進(jìn)行對(duì)比,結(jié)果如圖7和圖8所示。圖7和圖8中黑色實(shí)線代表實(shí)際模型的反演結(jié)果,紅色虛線代表利用本文方法直接反演的結(jié)果,藍(lán)色虛線代表利用常規(guī)方法間接計(jì)算的結(jié)果。通過比較發(fā)現(xiàn),無噪聲時(shí),兩種方法得到的Russell流體因子及剪切模量均與模型數(shù)據(jù)吻合較好,只有密度項(xiàng)的反演結(jié)果誤差相對(duì)較大,這與2.2節(jié)中的反演敏感性分析結(jié)論一致;當(dāng)信噪比為2∶1時(shí),兩種方法的估算精度均受到一定的影響,但與常規(guī)方法相比,本文方法直接反演的Russell流體因子以及剪切模量的總體趨勢(shì)更接近實(shí)際模型,但是兩種方法得到的密度項(xiàng)則均存在較大的差異,這是因?yàn)榉瓷湎禂?shù)近似方程本身缺少密度信息,反演得到的密度參數(shù)穩(wěn)定性不高,且易受噪聲影響產(chǎn)生較大的反演誤差。

圖5　地層模型曲線

圖6　模型角度道集a 無噪聲角度道集; b 信噪比為2∶1的角度道集

圖7　無噪聲反演結(jié)果

圖8　信噪比為2∶1情況下的反演結(jié)果

考慮到Russell流體因子是參與儲(chǔ)層流體識(shí)別的主體數(shù)據(jù),為進(jìn)一步比較本文方法與常規(guī)方法提取Russell流體因子的精度,對(duì)比了兩種方法估算的Russell流體因子與模型數(shù)據(jù)的相對(duì)誤差,結(jié)果如圖9所示。圖9中紅色虛線代表直接反演結(jié)果的相對(duì)誤差;藍(lán)色虛線代表常規(guī)方法計(jì)算結(jié)果的相對(duì)誤差。從圖9可見,在無噪聲情況下,兩種方法估算的Russell流體因子的相對(duì)誤差基本一致,但是,當(dāng)有噪聲時(shí),常規(guī)間接計(jì)算方法的相對(duì)誤差更大,這是因?yàn)橛性肼暤那闆r下,密度項(xiàng)的反演結(jié)果會(huì)存在較大的不確定性(圖8),利用這種密度項(xiàng)進(jìn)行Russell流體因子的間接計(jì)算會(huì)進(jìn)一步造成誤差累計(jì),使得間接計(jì)算的Russell流體因子存在較大誤差,而直接反演方法則有效規(guī)避了這種現(xiàn)象。考慮到實(shí)際地震數(shù)據(jù)均會(huì)存在一定的噪聲,所以本文方法直接反演的Russell流體因子具有更高的可靠性。

圖9　無噪聲(a)和信噪比為2∶1(b)情況下反演的Russell流體因子與實(shí)際模型數(shù)據(jù)的相對(duì)誤差

3.3實(shí)際數(shù)據(jù)應(yīng)用

實(shí)際資料來自加拿大Alberta地區(qū)的疊前道集數(shù)據(jù),選取了100個(gè)CMP道集,每個(gè)道集含有8道,偏移距為70～560m。圖10為選用的道集數(shù)據(jù)剖面,在剖面0.64～0.66s用黑色橢圓標(biāo)出了含氣儲(chǔ)層位置。假設(shè)該數(shù)據(jù)不含噪聲,加入白噪聲使得數(shù)據(jù)信噪比為5∶1,利用本文方法直接反演得到Russell流體因子數(shù)據(jù),結(jié)果如圖11a所示;利用常規(guī)AVO三參數(shù)反演方法反演得到縱、橫波速度與密度參數(shù),再根據(jù)公式(1)計(jì)算得到的Russell流體因子如圖11b所示。從圖11中可以看到,雖然兩者在含氣目的儲(chǔ)層位置均表現(xiàn)為較為明顯的低值異常,但是相比較直接反演結(jié)果,間接計(jì)算結(jié)果的信噪比明顯偏低,且產(chǎn)生了較大累計(jì)計(jì)算誤差。

圖10　測(cè)試道集數(shù)據(jù)(黑色橢圓表示含氣儲(chǔ)層位置)

圖11　測(cè)試數(shù)據(jù)的Russell流體因子計(jì)算結(jié)果a 直接定量反演結(jié)果; b 常規(guī)方法間接計(jì)算結(jié)果

為了更加清晰地對(duì)比反演結(jié)果的可靠性,圖12 給出了疊前反演得到的Russell流體因子曲線和實(shí)際測(cè)井曲線(為減小尺度差異,此處將測(cè)井?dāng)?shù)據(jù)進(jìn)行了Backus平均處理)。由圖12可見,直接反演的Russell流體因子與實(shí)際測(cè)井曲線吻合度更高,具有更高的抗噪性,而常規(guī)方法得到的Russell流體因子則存在較大誤差,這與模型測(cè)試的結(jié)果也是吻合的,即不可靠的密度項(xiàng)參數(shù)造成了Russell流體因子的累計(jì)誤差。通過實(shí)際數(shù)據(jù)的應(yīng)用也進(jìn)一步驗(yàn)證了基于新反射系數(shù)近似方程直接反演得到的Russell流體因子更為精確。

圖12　疊前反演得到的Russell流體因子與實(shí)際測(cè)井曲線的對(duì)比a 基于新反射系數(shù)近似方程的反演結(jié)果與測(cè)井曲線的對(duì)比; b 常規(guī)疊前反演結(jié)果與測(cè)井曲線的對(duì)比

4結(jié)束語

包含Russell流體因子的反射系數(shù)近似方程更直觀地表征了孔隙介質(zhì)中的流體彈性效應(yīng)。通過疊前地震反演直接提取Russell流體因子,不僅實(shí)現(xiàn)了疊前反演與儲(chǔ)層流體識(shí)別的有效結(jié)合,而且為流體檢測(cè)提供了更加可靠的數(shù)據(jù)基礎(chǔ)。主要結(jié)論如下:

1) 結(jié)合Biot-Gassmann孔隙彈性介質(zhì)理論推導(dǎo)的包含Russell流體因子的反射系數(shù)近似方程不僅可以滿足精度要求,而且Russell流體因子在公式中的信息量可以保證參數(shù)反演的穩(wěn)定性;

2) 基于貝葉斯理論框架的反演方法實(shí)現(xiàn)了Russell流體因子參數(shù)的可靠反演,有效避免了常規(guī)方法的累計(jì)誤差,提高了Russell流體因子的估算精度,有助于儲(chǔ)層流體類型的可靠判識(shí);

3) 在實(shí)際資料的應(yīng)用中,需要借助各向異性動(dòng)校正等技術(shù)對(duì)疊前道集進(jìn)行處理,在保證遠(yuǎn)道大偏移距數(shù)據(jù)可用的前提下,盡可能利用大角度疊前數(shù)據(jù)提高Russell流體因子反演的可靠性。

參考文獻(xiàn)

[1]印興耀,張世鑫,張峰.針對(duì)深層流體識(shí)別的兩項(xiàng)彈性阻抗反演與Russell流體因子直接估算方法研究[J].地球物理學(xué)報(bào),2013,56(7):2378-2390

YIN X Y,ZHANG S X,ZHANG F.Two term elastic impedance inversion and Russell fluid factor direct estimation method for deep reservoir fluid identification[J].Chinese Journal of Geophysics,2013,56(7):2378-2390

[2]鄭靜靜,印興耀,張廣智.流體因子關(guān)系分析以及新流體因子的構(gòu)建[J].地球物理學(xué)進(jìn)展,2011,26(2):579-587

ZHENG J J,YIN X Y,ZHANG G Z.Fluid factor analysis and the construction of the new fluid factor[J].Progress in Geophysics,2011,26(2):579-587

[3]管路平.地震疊前反演與直接烴類指示的探討[J].石油物探,2008,47(3):228-234

GUAN L P.Pre-stack inversion and direct hydrocarbon indicator[J].Geophysical Prospecting for Petroleum,2008,47(3):228-234

[4]郝前勇,印興耀,王玉梅,等.彈性模量流體因子在永新工區(qū)儲(chǔ)層預(yù)測(cè)中的應(yīng)用[J].石油物探,2012,51(5):502-507

HAO Q Y,YIN X Y,WANG Y M,et al.Application of elastic modulus fluid factor on reservoir prediction at Yongxin Block[J].Geophysical Prospecting for Petroleum,2012,51(5):502-507

[5]李超,印興耀,張廣智,等.基于貝葉斯理論的孔隙流體模量疊前AVA反演[J].石油物探,2015,54(4):467-476

LI C,YIN X Y,ZHANG G Z,et al.Prestack AVA inversion for pore fluid modulus based on the Bayesian theory[J].Geophysical Prospecting for Petroleum,2015,54(4):467-476

[6]SMITH G C,GIDLOW P M.Weighted stacking for rock property estimation and detection of gas[J].Geophysical Prospecting,1987,35(8):993-1014

[7]GIDLOW P M,SMITH G C,VAIL P J.Hydrocarbon detection using fluid factor traces,a case study:how useful is AVO analysis?[J].Technical Program and Abstracts of Joint SEG/EAEG Summer Research Workshop,1992:78-89

[8]FATTI J L,VAIL P,SMITH G C,et al.Detection of gas in sandstone reservoirs using AVO analysis:a 3-D seismic case history using the geostack technique[J].Geophysics,1994,59(9):1362-1376

[9]GRAY D,GOODWAY B,CHEN T.Bridging the gap:using AVO to detect changes in fundamental elastic constants[J].Expanded Abstracts of 69thAnnual Internat SEG Mtg,1999:852-855

[10]GOODWAY B,CHEN T,DOWNTON J.Improved AVO fluid detection and lithology discrimination using Lamé petrophysical parameters;“λρ”,“μρ”,&“λ/μfluid stack”,from P and S inversions[J].Expanded Abstracts of 67thAnnual Internat SEG Mtg,1997:183-186

[11]QUAKENBUSH M,SHANG B,TUTTLE C.Poisson impedance[J].The Leading Edge,2006,25(2):128-138

[12]RUSSELL B H,HEDLIN K,HILTERMAN F J.Fluid-property discrimination with AVO:a Biot-Gassmann perspective[J].Geophysics,2003,68(1):29-39

[13]RUSSELL B H,GRAY D,HAMPSON D P.Linearized AVO and poroelasticity[J].Geophysics,2011,76(3):C19-C29

[14]張世鑫.基于地震信息的流體識(shí)別方法研究與應(yīng)用[D].青島:中國石油大學(xué)(華東),2012

ZHANG S X.Methodology and application of fluid identification with seismic information[D].Qingdao:China University of Petroleum,2012

[15]AKI K,RICHARDS P G.Quantitative seismology[M].2nd ed.San Francisco:W H Freeman and Co,2002:153-154

[16]TARANTOLA A.Inverse problem theory and methods for model parameter estimation[M].Philadelphia:SIAM,2005:5-8

[17]王保麗,印興耀,張繁昌.基于Gray近似的彈性波阻抗方程及反演[J].石油地球物理勘探,2007,42(4):435-439

WANG B L,YIN X Y,ZHANG F C.Gray approximation-based elastic wave impedance equation and inversion[J].Oil Geophysical Prospecting,2007,42(4):435-439

[18]印興耀,張世鑫,張繁昌,等.利用基于Russell近似的彈性波阻抗反演進(jìn)行儲(chǔ)層描述和流體識(shí)別[J].石油地球物理勘探,2010,45(3):373-380

YIN X Y,ZHANG S X,ZHANG F C,et al.Utilizing Russell approximation based elastic wave impedance inversion to conduct reservoir description and fluid identification[J].Oil Geophysical Prospecting,2010,45(3):373-380

[19]楊培杰,王長江,畢俊鳳,等.可變點(diǎn)約束疊前流體因子直接提取方法[J].地球物理學(xué)報(bào),2015,58(6):2188-2200

YANG P J,WANG C J,BI J F,et al.Direct extraction of the fluid factor based on variable point-constraint[J].Chinese Journal of Geophysics,2015,58(6):2188-2200

[20]BIOT M A.Mechanics of deformation and acoustic propagation in porous median,low frequency range[J].Journal of the Acoustical Society of America,1956,28(2):168-178

[21]GASSMANN F.über die elastizit?t por?ser median.Vier.Der Natur[J].Gesellschaft in Zürich,1951,96:1-23

[22]URSIN B,TJALAND E.Information content of the elastic reflection matrix[J].Expanded Abstracts of 62ndAnnual Internat SEG Mtg,1992:796-799

[23]NICOLAO A D,DRUFUCA G,ROCCA F.Eigenvalues and eigenvectors of linearized elastic inversion[J].Geophysics,1993,58(5):670-679

[24]WUBSHET M A.Regularization of the AVO inversion problem by means of a multivariate Cauchy probability distribution[D].Alberta:University of Alberta,2010

[25]張世鑫,印興耀,張繁昌.基于三變量柯西分布先驗(yàn)約束的疊前三參數(shù)反演方法[J].石油地球物理勘探,2011,46(5):737-743

ZHANG S X,YIN X Y,ZHANG F C.Prestack three term inversion method based on trivariate Cauchy distribution prior constraint[J].Oil Geophysical Prospecting,2011,46(5):737-743

[26]印興耀,周琪超,宗兆云,等.基于t分布為先驗(yàn)約束的疊前AVO反演[J].石油物探,2014,53(1):84-92

YIN X Y,ZHOU Q C,ZONG Z Y,et al.AVO inversion witht-distribution as priori constraint[J].Geophysical Prospecting for Petroleum,2014,53(1):84-92

(編輯:陳杰)

A reflection coefficient approximation equation of poroelastic media for Russell fluid factor estimation

ZHANG Shixin,DU Xiangdong,HAN Wenming,SUN Linjie

(CNOOCResearchInstitute,Beijing100027,China)

Abstract:The Russell fluid factor is an important sensitive parameter for reservoir pore-fluid identification and the estimation techniques for Russell fluid factor plays a key role to the quality of reservoir fluid identification.The prestack seismic inversion based on the reflection coefficient approximation equation is the main approach to estimate fluid factors.With the help of Biot-Gassmann poroelastic theory,the Aki-Richards reflection coefficient approximation is rearranged to get a new linear approximation equation containing the Russell fluid factor parameter.Analytical analysis is carried out to check the accuracy and inversion sensitivity of the new approximation equation,which shows the reliability of the new reflection coefficient approximation and high inversion sensitivity of Russell fluid factor term.Finally,the directly quantitative inversion of pore-fluid impedance is realized based on the Bayesian AVO inversion.The model test and real application denote the feasibility and stability of this Russell fluid factor estimation technique.This is proposed as a new fluid estimation approach for reservoir fluid identification.

Keywords:Russell fluid factor,reflection coefficient approximation equation,prestack seismic inversion,fluid identification

文章編號(hào):1000-1441(2016)02-0178-10

DOI:10.3969/j.issn.1000-1441.2016.02.003

中圖分類號(hào):P631

文獻(xiàn)標(biāo)識(shí)碼:A

基金項(xiàng)目:國家科技重大專項(xiàng)(2011ZX05030-004)資助。

作者簡(jiǎn)介:張世鑫(1985—),男,博士,主要從事地震資料解釋、儲(chǔ)層油氣預(yù)測(cè)等方面的研究工作。

收稿日期:2015-03-24;改回日期:2015-12-25。

This research is financially supported by the National Science and Technology Major Project of China (Grant No.2011ZX05030-004).