[發明專利]基于AVO特征的剩余振幅補償方法

有效

• 申請號: 201410194584.0
• 申請日: 2014-05-09
• 公開（公告）號: CN103954998B
• 公開（公告）日: 2017-01-18
• 發明（設計）人: 劉力輝
• 申請（專利權）人: 成都晶石石油科技有限公司
• 主分類號: G01V1/36
• 分類號: G01V1/36
• 代理公司: 成都行之專利代理事務所(普通合伙)51220
• 代理人: 王記明

1.基于AVO特征的剩余振幅補償方法，其特征在于包括以下步驟：

（a）疊前CRP道集，計算參數；

（b）在滑動時窗計算每個時窗內的振幅值；

（c）計算CRP道集內（2j+1）道信號的時變振幅衰減曲線；

（d）給一定道間隔增量判斷是否在偏移距范圍內，如果是，則重復進行步驟（c），如果不是則進行下一步驟（e）；

（e）計算不同偏移距時空變振幅衰減曲線，得到處理后的道集振幅；

（f）將步驟（e）得到的道集振幅與補償目標道的振幅進行計算，利用均方根速度公式，求取偏移距對應入射角及相應的振幅值；

（g）將步驟（f）得到的結果進行計算并比較結果：如果剔除道的百分比小于設定值，則保留誤差小的振幅數據，并利用最小二乘擬合求系數R、W、V；如果剔除道的百分比不小于設定值，則利用最小二乘擬合求系數R、W、V；

（h）根據步驟（g）得到的擬合求系數R、W、V，重新確定偏移距與振幅的關系，補償修正目標道。

2.根據權利要求1所述的基于AVO特征的剩余振幅補償方法，其特征在于，所述步驟（b）在滑動時窗計算每個時窗內的振幅值是指：假設道集內有N道地震信號，計算以第i道為中心，左右相鄰j道共計（2j+1）道，在時窗tu-td內的振幅值；具體過程如下：

（b1）統計時窗內（2j+1）道信號的波峰和波谷總數和對應的振幅極值，即波谷取絕對值；

（b2）對振幅極值進行排序；

（b3）利用線性擬合法或中值計算法計算代表該時窗的振幅值。

3.根據權利要求1所述的基于AVO特征的剩余振幅補償方法，其特征在于，所述步驟（c）計算CRP道集內（2j+1）道信號的時變振幅衰減曲線包括以下步驟：

（c1）給定計算時變振幅衰減曲線的起始時間、結束時間及時窗長度，采用滑動時窗法計算每個時窗的極值振幅值或平均振幅值；

（c2）根據每個時窗的極值振幅值或平均振幅值進行插值、平滑處理，計算出每個時間采樣點的極值振幅值或平均振幅值，從而得到時變振幅衰減曲線。

4.根據權利要求1所述的基于AVO特征的剩余振幅補償方法，其特征在于，（e）計算不同偏移距時空變振幅衰減曲線，得到處理后的道集振幅步驟在于：在偏移距方向，每隔一定道增量計算一條時變振幅衰減曲線，相鄰曲線間通過插值就可以得到每道的時變振幅衰減曲線。

5.根據權利要求1所述的基于AVO特征的剩余振幅補償方法，其特征在于，所述步驟（g）中最小二乘擬合求系數的步驟如下：

（g1）根據以下公式描述的振幅隨偏移距變化滿足拋物線的規律???????????????????????????????????????????????、已知的反射系數項y和x，用最小二乘逐點進行AVO拋物線擬合，

（g2）加入剔除擬合方法，求得拋物線的系數R、W、V值，

（g3）根據求得的拋物線系數以及已知的x值再代步驟（g1）中重新求得新的地震信號y值，并根據用戶需求將近道、遠端能量弱的道進行補償或者進行近道多次波的剔除。

6.根據權利要求5所述的基于AVO特征的剩余振幅補償方法，其特征在于，所述步驟（g2）中的剔除擬合方法包括以下步驟：

（g21）優先剔除小炮檢距處的大誤差點子，若炮檢距x處的誤差ex，則加權誤差為：E_x?=?W_xex；

權系數W_x為W_x=1+Aexp[-B(x/x_max)]

其中：A及B為兩個正的常數；

（g22）利用中值或均值計算各道誤差：先求取所有道的均值，然后計算各道的振幅與均值之差e_x，并把差值當作誤差e_x，重新帶入步驟（g21）計算，得到加權后的誤差E_x。

7.根據權利要求1至6中任意一項所述的基于AVO特征的剩余振幅補償方法，其特征在于：所述步驟（g）中，剔除道的百分比設定值為20%。