基于入射角的avo近似公式進行屬性提取的方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明設(shè)及地震資料處理領(lǐng)域，是一種利用地震數(shù)據(jù)進行儲層識別的方法。
【背景技術(shù)】
[0002] AV0技術(shù)是根據(jù)振幅隨炮檢距的變化規(guī)律所反映出的地下巖性及其孔隙流體的性 質(zhì)來直接預(yù)測油氣和估計地殼巖性參數(shù)的一項技術(shù)。其理論基礎(chǔ)是描述平面波在水平分界 面上反射和折射的Zo巧pritz方程。完整的Zo巧pritz方程形式復(fù)雜，物理意義不明確， 無法直接應(yīng)用于地震數(shù)據(jù)的處理和參數(shù)反演。盡管人們早就認(rèn)識到了該一點，但在早期的 研究中，人們僅從數(shù)值計算的角度對Zoeppritz描述的反射系數(shù)進行了各種討論和近似。 Muskat系統(tǒng)地計算了彈性介質(zhì)中平面波的反射和透射系數(shù)；Koe化ed利用17組不同彈性參 數(shù)的模型，詳細地研究了泊松比對兩個各向同性介質(zhì)之間的反射界面所產(chǎn)生的產(chǎn)生的反射 系數(shù)的影響，并第一次給出了將泊松比與反射系數(shù)直接聯(lián)系起來的Zoeppritz近似方程。 Bcxrtfeld第一個給出了物理意義比較明確的反射系數(shù)近似計算公式，使得基于Zo巧pritz 方程的反射和透射問題真正應(yīng)用到了勘探地球物理領(lǐng)域。在W上工作的基礎(chǔ)上，不少人做 了大量的工作，分別從不同的方面對Zoeppritz方程進行了詳細的討論和簡化，提出了不 同假設(shè)條件下的縱波反射振幅的近似表達式。但是所有的近似公式中角度都是應(yīng)用的上層 入射角和下層透射角的平均，給進一步的參數(shù)估算帶來一定的誤差。

【發(fā)明內(nèi)容】

[0003] 本發(fā)明的目的在于避免現(xiàn)有技術(shù)的不足而提供基于入射角的AV0近似公式進行 屬性提取的方法。
[0004] 本發(fā)明的技術(shù)方案是：基于入射角的AV0近似公式進行屬性提取的 方法，實現(xiàn)過程包括；1)基于入射角的AV0近似公式的提出；2)基于入射角 的AV0近似公式的理論推導(dǎo)；3)應(yīng)用常規(guī)方法提取P、G，Rs的屬性；4)將屬性歸 一化；5)由井?dāng)?shù)據(jù)或者研究區(qū)域的經(jīng)驗公式確定系數(shù)b值，然后應(yīng)用公式（14)
對常規(guī)方法提取的屬性進行修正；6)將 修正結(jié)果與井?dāng)?shù)據(jù)進行對比分析，根據(jù)分析結(jié)果對修正結(jié)果進行校正，得到最終的屬性；其 中，P是由零炮檢距構(gòu)成的地震道，即縱波疊加道，代表對反射界面兩側(cè)的波阻抗變化的響 應(yīng)；斜率G為梯度疊加道，代表對橫波速度、縱波速度和體密度變化的響應(yīng)，是由入射角近 似公式計算得到的，GO代表對橫波速度、縱波速度和體密度變化的響應(yīng)，由常規(guī)近似公式計 算的得到；Rs表示由入射角近似公式計算的橫波阻抗反射系數(shù)，馬》為由常規(guī)近似公式計算 的橫波阻抗反射系數(shù)； 具體的求取方法為： 1)由Aki和化U巧近似可得；
R( 0 )表示縱波反射系數(shù)，妍和巧分別為反射界面上、下地層介質(zhì)的密度，和分 別是反射界面上下上、下地層的縱波速度，分別是上、下地層的橫波速度，巧和巧分 別是地層的入射角和透射角 2)根據(jù)Snell定律有；
對公式（1)左邊項進行泰勒展開，及
將公式（3)代入公式（2)可得：
令
4)Garner等人于1974提出的密度與縱波速度之間的經(jīng)驗公式的形式為[6];
其中，a和b為系數(shù)（一般取a=0. 31，b=0. 25) 由上式可得
即：
5) -般由化uey近似計算截距和梯度屬性P和G(進一步可計算橫波阻抗反射系數(shù)Rs 等屬性）的公式為：
其中，巧是根據(jù)平均角計算出的梯度 6) 在實際應(yīng)用公式（11)計算屬性時，由于求取角度A不方便，常常用入射角S代替6進 行估算，即應(yīng)用下式計算P和Rs:
比較式（11)和式（13)可得：
所W，計算橫波阻抗反射系數(shù)Rs的公式為：
該里的P和Rs。由公式（13 )計算得到。
[0005] 與現(xiàn)有技術(shù)相比，本發(fā)明的有益效果是：計算方便、精確，與實際的參數(shù)誤差小。
【附圖說明】
[0006] 圖 1 為Ostrander模型。
[0007] 圖 2 為Goodway模型。
[0008] 圖3為MarmousiII模型中的氣層。
[0009] 圖4為修改前后誤差對比分析。
[0010] 圖5為從mar2模型中抽取的一道數(shù)據(jù)。
[0011] 圖6為MarmousiII模型的屬性對比。
[0012] 圖7MarmousiII模型橫波合成地震記錄的比較。其中（a)為垂直入射計算的 Rs與子波權(quán)積得到的合成地震記錄?；槌Ｒ?guī)方法計算的Rs與子波權(quán)積得到的合成地 震記錄。（C)為修正之后的Rs與子波權(quán)積得到的合成地震記錄。
[0013] 圖8為實施例的疊前角度道集。其中（a)為某地區(qū)8°角度疊加道集。（b)為某 地區(qū)16°角度疊加道集。（C)為某地區(qū)24°角度疊加道集。
[0014] 圖9為實施例基于修正公式計算的Rs屬性。
【具體實施方式】
[0015] 基于入射角的AV0近似公式進行屬性提取的方法，實現(xiàn)過程包括；1)基于入射角 的AV0近似公式的提出；2)基于入射角的AV0近似公式的理論推導(dǎo)；3)應(yīng)用常規(guī)方法提取 P、G，Rs的屬性；4)將屬性歸一化；5)由井?dāng)?shù)據(jù)或者研究區(qū)域的經(jīng)驗公式確定系數(shù)b值，然 后應(yīng)用公式（14)
對常規(guī)方法提取的屬性 進行修正；6)將修正結(jié)果與井?dāng)?shù)據(jù)進行對比分析，根據(jù)分析結(jié)果對修正結(jié)果進行校正，得 到最終的屬性；其中，P是由零炮檢距構(gòu)成的地震道，即縱波疊加道，代表對反射界面兩側(cè) 的波阻抗變化的響應(yīng)；斜率G為梯度疊加道，代表對橫波速度、縱波速度和體密度變化的響 應(yīng)，是由入射角近似公式計算得到的，GO代表對橫波速度、縱波速度和體密度變化的響應(yīng)， 由常規(guī)近似公式計算的得到；Rs表示由入射角近似公式計算的橫波阻抗反射系數(shù)，馬C為由 常規(guī)近似公式計算的橫波阻抗反射系數(shù)； 具體的求取方法為： 1) 由Aki和化U巧近似可得；
R( 0 )表示縱波反射系數(shù)，A和A分別為反射界面上、下地層介質(zhì)的密度，?和分 別是反射界面上下上、下地層的縱波速度Ki和吃分別是上、下地層的橫波速度，巧和巧分 別是地層的入射角和透射角 2) 根據(jù)Snell定律有；
對公式（1)左邊項進行泰勒展開，及
3)該樣可W得到基于入射角的AVO近似公式：
4) Ga