[發明專利]一種基于HTI介質的各向異性單參數反演方法

有效

• 申請號: 201410369859.X
• 申請日: 2014-07-30
• 公開（公告）號: CN105319590B
• 公開（公告）日: 2018-03-06
• 發明（設計）人: 張克非;肖鵬飛;李呈呈;胡華鋒;林正良;朱童
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油物探技術研究院
• 主分類號: G01V1/30
• 分類號: G01V1/30;G06F19/00
• 代理公司: 北京思創畢升專利事務所11218
• 代理人: 郭韞

技術領域

本發明屬于地球物理反演領域，具體涉及一種基于HTI介質的各向異性單參數反演方法，適用于HTI(橫向各向異性介質)。

背景技術

隨著地質勘探的不斷發展，對裂縫性油氣藏的研究也越來越廣泛，在裂縫預測中對裂縫參數的定量預測主要通過各向異性參數得到。通過疊前方位道集獲取各向異性參數目前主要分為三種方法：

1.通過分析NMO速度的方位變化計算裂縫參數，當速度分析精度比較高時，分析NMO速度的方位變化就能夠做方位各向異性分析，但在實際問題中，很難確定速度的方位各向異性；

2.通過疊前方位道集中提取剩余方位時差計算，疊前道集中剩余時差比較容易拾取，該方法易于實現，但是由于拾取剩余時差時存在較多誤差，無法區分產生剩余方位時差的原因是否是各向異性引起的，導致該方法計算結果準確性較低；

3.利用縱波反射振幅隨偏移距(入射角)和方位變化的信息，以Ruger(1998)方程為基礎，通過推導近似等轉換來化簡Ruger方程，計算出與各向異性參數相關的綜合參數，這些方法大多相對復雜，假設條件多并且計算效率不高。

發明內容

本發明的目的在于解決上述現有技術中存在的難題，提供一種基于HTI介質的各向異性單參數反演方法，基于簡化的Ruger方程，通過疊后反演的方法解決疊前各向異性參數反演問題。

本發明是通過以下技術方案實現的：

一種基于HTI介質的各向異性單參數反演方法，包括：

(1)選取任意三個不同方位角對應的疊后地震數據體，通過稀疏脈沖的反演方法計算得到對應該三個不同方位角的反射系數

(2)分析疊前方位道集獲得入射角信息，所述入射角信息包括最大入射角θ₂、最小入射角θ₁和間隔，然后通過公式計算出B的值；

(3)利用方位角信息反射系數和參數B，計算得到各向異性參數Δε^(V)，完成反演。

所述方法進一步包括：

(4)將步驟(3)反演得到的Δε^(V)與井旁道或模型數據的Δε^(V)值進行對比，分析反演得到的Δε^(V)值與井旁道或模型數據的Δε^(V)值之間的誤差。

所述步驟(2)中，0≤θ₁＜θ₂≤30。

所述步驟(3)是利用下面的公式計算得到各向異性參數Δε^(V)：

與現有技術相比，本發明的有益效果是：

(1)目前疊前反演受疊前道集的影響，穩定性不高，本方法將疊前道集在入射角域進行疊加，將疊加后的具有方位信息的疊后方位體作為反演的輸入，使反演方法具有疊后地震反演穩定性高、信噪比高的優點；

(2)目前疊前各向異性參數反演的理論要求輸入數據體的個數盡量多，反演方法則更加穩定，而實際應用過程中由于計算效率及數據品質的限制難以得到多個輸入數據體，本反演方法從理論上只需要3個疊前方位體數據，在實際應用過程中易于實現；

(3)目前由于觀測系統的影響，絕大多數疊前角度域道集存在，近、遠道能量弱，中入射角道集能量強的缺陷，本反演方法從理論上實現了疊前角度域道集的數據篩選功能，可直接選取地震數據品質高的部分進行疊加作為反演輸入數據；

(4)反演參數少，計算方法簡潔。

附圖說明

圖1本方法的步驟框圖

圖2理論模型圖

圖3正演得到的0°方位的疊后數據體

圖4正演得到的30°方位的疊后數據體

圖5正演得到的60°方位的疊后數據體

圖6反演得到的0°方位疊后數據體對應的反射系數

圖7反演得到的30°方位疊后數據體對應的反射系數

圖8反演得到的60°方位疊后數據體對應的反射系數

圖9計算得到的各向異性參數Δε^(V)值

圖10模型中各向異性參數Δε^(V)值。

具體實施方式

下面結合附圖對本發明作進一步詳細描述：

直接進行疊前各向異性參數反演時，Ruger方程內同時含有縱波速度、橫波速度、密度和Δδ^(V)，Δε^(V)和Δγ六個未知變量，這就使方程求解時引入多種近似關系且解法復雜，穩定性較低。