[發明專利]Radon域能量掃描疊加的微地震事件定位方法

有效

• 申請號: 201811292105.3
• 申請日: 2018-10-31
• 公開（公告）號: CN109188516B
• 公開（公告）日: 2021-07-20
• 發明（設計）人: 王東凱;王延光;苗永康;尚新民;李繼光;劉鵬翔;徐強;隋旭強
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司勝利油田分公司物探研究院
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 濟南日新專利代理事務所(普通合伙) 37224
• 代理人: 崔曉艷

1. Radon域能量掃描疊加的微地震事件定位方法，其特征在于，該Radon域能量掃描疊加的微地震事件定位方法包括：

步驟1，輸入微地震數據及信息，并進行空間網格化；

步驟2，計算單個網格點Radon域數據投影；

步驟3，搜索Radon域能量疊加極大值E_i；

步驟4，重復步驟2-3,計算完成所有網格對應的Radon域能量疊加極大值；

步驟5，確定疊加能量最大值，并進行最值點三維空間切片；

步驟6，輸出震源空間位置，完成微地震事件定位。

2.根據權利要求1所述的Radon域能量掃描疊加的微地震事件定位方法，其特征在于，在步驟1中，輸入微地震監測觀測系統采集到的微地震數據。

3.根據權利要求1所述的Radon域能量掃描疊加的微地震事件定位方法，其特征在于，在步驟1中，綜合考慮精度需求與硬件條件進行地下空間網格化。

4.根據權利要求1所述的Radon域能量掃描疊加的微地震事件定位方法，其特征在于，在步驟2中，選定單個網格點i作為假設震源位置，計算該震源與檢波器相對空間位置關系下的微地震記錄Radon域數據投影；根據微地震記錄同向軸的不同空間形態，選取線性和非線性不同的積分路徑進行Radon變換，實現Radon域能量聚焦。

5.根據權利要求4所述的Radon域能量掃描疊加的微地震事件定位方法，其特征在于，在步驟2中，采用的線性Radon變換過程為：

（1）

其中S為投影變換前的時間空間域微地震數據，為與網格點對應的變換后的Radon域數據，M為需要遍歷的震源網格總數，為截距時間，p為射線參數，x為假設震源與檢波器間的位置關系；

采用的非線性Radon變換過程為：

（2）

其中S為投影變換前的時間空間域微地震數據，為與網格點對應的變換后的Radon域數據，M為需要遍歷的震源總數，為截距時間，q為曲率參數，x為假設震源與檢波器間的位置關系。

6.根據權利要求5所述的Radon域能量掃描疊加的微地震事件定位方法，其特征在于，在步驟3中，對應單一網格點i，掃描并記錄已得到的Radon域數據的能量疊加極大值E_i，表示為：

（3）。

7.根據權利要求1所述的Radon域能量掃描疊加的微地震事件定位方法，其特征在于，在步驟4中，重復步驟2-3，計算對應下一個網格點的Radon域數據并記錄能量疊加極大值，直至計算完成所有假設的震源點對應的Radon域數據。

8.根據權利要求1所述的Radon域能量掃描疊加的微地震事件定位方法，其特征在于，在步驟5中，掃描對應所有M個震源網格點的Radon域數據能量疊加極大值并確定最大值E_max，表示為：

（4）。

9.根據權利要求8所述的Radon域能量掃描疊加的微地震事件定位方法，其特征在于，在步驟5中，提取對應E_max的三維空間能量切片，并通過三維切片的疊加能量最大值確定對應的空間位置。

10.根據權利要求1所述的Radon域能量掃描疊加的微地震事件定位方法，其特征在于，在步驟6中，當需進一步提高精度時，在該位置處繼續剖分加密網格并重復步驟2-6。