亚洲狠狠干,亚洲国产福利精品一区二区,国产八区,激情文学亚洲色图

一種基于hti介質(zhì)的各向異性單參數(shù)反演方法

文檔序號:9563616閱讀:807來源:國知局
一種基于hti介質(zhì)的各向異性單參數(shù)反演方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明屬于地球物理反演領(lǐng)域,具體涉及一種基于HTI介質(zhì)的各向異性單參數(shù)反 演方法,適用于HTI (橫向各向異性介質(zhì))。
【背景技術(shù)】
[0002] 隨著地質(zhì)勘探的不斷發(fā)展,對裂縫性油氣藏的研究也越來越廣泛,在裂縫預(yù)測中 對裂縫參數(shù)的定量預(yù)測主要通過各向異性參數(shù)得到。通過疊前方位道集獲取各向異性參數(shù) 目前主要分為三種方法:
[0003] 1.通過分析NMO速度的方位變化計算裂縫參數(shù),當(dāng)速度分析精度比較高時,分析 NMO速度的方位變化就能夠做方位各向異性分析,但在實際問題中,很難確定速度的方位各 向異性;
[0004] 2.通過疊前方位道集中提取剩余方位時差計算,疊前道集中剩余時差比較容易拾 取,該方法易于實現(xiàn),但是由于拾取剩余時差時存在較多誤差,無法區(qū)分產(chǎn)生剩余方位時差 的原因是否是各向異性引起的,導(dǎo)致該方法計算結(jié)果準(zhǔn)確性較低;
[0005] 3.利用縱波反射振幅隨偏移距(入射角)和方位變化的信息,以Ruger (1998)方 程為基礎(chǔ),通過推導(dǎo)近似等轉(zhuǎn)換來化簡Ruger方程,計算出與各向異性參數(shù)相關(guān)的綜合參 數(shù),這些方法大多相對復(fù)雜,假設(shè)條件多并且計算效率不高。

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

[0006] 本發(fā)明的目的在于解決上述現(xiàn)有技術(shù)中存在的難題,提供一種基于HTI介質(zhì)的各 向異性單參數(shù)反演方法,基于簡化的Ruger方程,通過疊后反演的方法解決疊前各向異性 參數(shù)反演問題。
[0007] 本發(fā)明是通過以下技術(shù)方案實現(xiàn)的:
[0008] -種基于HTI介質(zhì)的各向異性單參數(shù)反演方法,包括:
[0009] (1)選取任意三個不同方位角4?,,中3對應(yīng)的疊后地震數(shù)據(jù)體,通過稀疏脈 沖的反演方法計算得到對應(yīng)該三個不同方位角的反射系數(shù)
[0010] (2)分析疊前方位道集獲得入射角信息,所述入射角信息包括最大入射角Θ 2、最 小入射角Θ i和間隔,然后通過公式
>計算出B的值;
[0011] (3)利用方位角信息%?物、反射系數(shù)〗
和參 數(shù)Β,計算得到各向異性參數(shù)△ ε (ν),完成反演。
[0012] 所述方法進一步包括:
[0013] (4)將步驟(3)反演得到的Λ ε (ν)與井旁道或模型數(shù)據(jù)的Λ ε (ν)值進行對比,分 析反演得到的Λ ε (ν)值與井旁道或模型數(shù)據(jù)的Λ ε (ν)值之間的誤差。
[0014] 所述步驟(2)中,0彡Θ i < Θ 2彡30。
[0015] 所述步驟(3)是利用下面的公式計算得到各向異性參數(shù)Λ ε (v):
[0017] 與現(xiàn)有技術(shù)相比,本發(fā)明的有益效果是:
[0018] (1)目前疊前反演受疊前道集的影響,穩(wěn)定性不高,本方法將疊前道集在入射角域 進行疊加,將疊加后的具有方位信息的疊后方位體作為反演的輸入,使反演方法具有疊后 地震反演穩(wěn)定性高、信噪比高的優(yōu)點;
[0019] (2)目前疊前各向異性參數(shù)反演的理論要求輸入數(shù)據(jù)體的個數(shù)盡量多,反演方法 則更加穩(wěn)定,而實際應(yīng)用過程中由于計算效率及數(shù)據(jù)品質(zhì)的限制難以得到多個輸入數(shù)據(jù) 體,本反演方法從理論上只需要3個疊前方位體數(shù)據(jù),在實際應(yīng)用過程中易于實現(xiàn);
[0020] (3)目前由于觀測系統(tǒng)的影響,絕大多數(shù)疊前角度域道集存在,近、遠(yuǎn)道能量弱,中 入射角道集能量強的缺陷,本反演方法從理論上實現(xiàn)了疊前角度域道集的數(shù)據(jù)篩選功能, 可直接選取地震數(shù)據(jù)品質(zhì)高的部分進行疊加作為反演輸入數(shù)據(jù);
[0021] (4)反演參數(shù)少,計算方法簡潔。
【附圖說明】
[0022] 圖1本方法的步驟框圖
[0023] 圖2理論模型圖
[0024] 圖3正演得到的0°方位的疊后數(shù)據(jù)體
[0025] 圖4正演得到的30°方位的疊后數(shù)據(jù)體
[0026] 圖5正演得到的60°方位的疊后數(shù)據(jù)體
[0027] 圖6反演得到的0°方位疊后數(shù)據(jù)體對應(yīng)的反射系數(shù)
[0028] 圖7反演得到的30°方位疊后數(shù)據(jù)體對應(yīng)的反射系數(shù)
[0029] 圖8反演得到的60°方位疊后數(shù)據(jù)體對應(yīng)的反射系數(shù)
[0030] 圖9計算得到的各向異性參數(shù)Λ ε (v)值
[0031] 圖10模型中各向異性參數(shù)Λ ε (v)值。
【具體實施方式】
[0032] 下面結(jié)合附圖對本發(fā)明作進一步詳細(xì)描述:
[0033] 直接進行疊前各向異性參數(shù)反演時,Ruger方程內(nèi)同時含有縱波速度、橫波速度、 密度和Λ S (ν),Λ ε (v)和Λ Y六個未知變量,這就使方程求解時引入多種近似關(guān)系且解法 復(fù)雜,穩(wěn)定性較低。
[0034] 本發(fā)明通過計算各個方位疊后地震數(shù)據(jù)的反射系數(shù),以化簡后的Ruger方程為基 礎(chǔ)計算出各向異性參數(shù)△ ε (ν),將疊后反演的穩(wěn)定性優(yōu)勢引入到解決疊前反演問題中,該 反演算法得到各向異性參數(shù)可以用于裂縫型儲層的精細(xì)預(yù)測。
[0035] Ruger(1998)借助弱各向異性的概念,推導(dǎo)出HTI介質(zhì)的縱波近似反射系數(shù)關(guān)系 式:
[0036]
[0037] 式中:S, P和p分別為HTI介質(zhì)上下兩層的縱波、橫波速度平均值以及密度平均 值;Λ δ (ν),Λ ε (v)和Λ Y分別為上下兩層各向異性參數(shù)差值,與裂縫密度和縫隙充填流 體有關(guān);Θ和平分別代表入射角和方位角,與裂縫的發(fā)育方向有關(guān)。
[0038] 根據(jù)Thomsen弱各向異性近似將均勻各向同性介質(zhì)的反射系數(shù)描述為HTI介質(zhì)的 各向同性背景將各向異性參數(shù)S (v),^"和Y看成各向異性擾動。因此各向異性介質(zhì)的 反射系數(shù)為均勻各向同性背景反射系數(shù)R PP_IS_}和各向異性中)反射系數(shù)之和。
[0043] 由于實際地震資料中疊前角度域道集的入射角范圍大多數(shù)在0-30°之間,因此令
各向異性項_Rp:P_.ANI(?> Φ)化簡為
[0045]為了提高反演的穩(wěn)定性,將各向異性反射系數(shù)按照入射角方向疊加, 方程化簡為
[0048]由于各向同性反射系數(shù)Rppii^e)不隨方位角變化,令
[0050] 其中 0 彡 Θ i < Θ 2 彡 30,
[0051] 則將公式6, 7中帶入公式2轉(zhuǎn)換為,
[0057] 解由方程(9) (10) (11)組成的方程組,得
[0059] 其中
分別代表
時的反射系數(shù),上述為HTI介質(zhì)中近似各向異性參數(shù)計算公式,在此基礎(chǔ) 上進行各向異性參數(shù)Λ ε (v)的反演研究。
[0060] 反射系數(shù)
通過稀疏脈沖反演計算,常規(guī)的稀疏脈沖 流程可分為三步:
[0061] (1)反射系數(shù)反演
[0062] 采用最大似然反褶積進行反射系數(shù)的反演,最大似然反褶積對地層的假設(shè)認(rèn)為: 地層的反射系數(shù)是由較大的反射界面的反射和具有高斯背景的小反射疊
當(dāng)前第1頁1 2 
網(wǎng)友詢問留言 已有0條留言
  • 還沒有人留言評論。精彩留言會獲得點贊!
1