[發明專利]微震震源定位方法、裝置、電子設備及介質

有效

• 申請號: 202210418079.4
• 申請日: 2022-04-21
• 公開（公告）號: CN114518595B
• 公開（公告）日: 2022-08-02
• 發明（設計）人: 常瑩;張達;楊小聰;冀虎;戴銳;石雅倩;吳彥博
• 申請（專利權）人: 礦冶科技集團有限公司
• 主分類號: G01V1/28
• 分類號: G01V1/28;G01V1/30;G06F17/16
• 代理公司: 北京超凡宏宇專利代理事務所(特殊普通合伙) 11463
• 代理人: 張萌

1.一種微震震源定位方法，其特征在于，包括：

獲取接收到微震事件的傳感器的位置坐標和所述傳感器接收到微震事件的微震到達時間；

基于第一個接收到微震事件的傳感器的位置坐標和所述微震到達時間確定初始震源位置和初始震源發生時間；

基于所述初始震源位置和所述初始震源發生時間對目標函數進行傳感器權重系數配置、系數矩陣變換求解以及迭代，直至滿足預設條件；其中，所述目標函數為地震波走時殘差的二范數；所述傳感器權重系數是基于地震波信噪比、地震波震源距和所述地震波走時殘差確定的，所述系數矩陣變換包括將地震波走時進行參數分離和歸一化處理；

將滿足所述預設條件時所述目標函數的解確定為震源位置和震源發生時間；

基于所述初始震源位置和所述初始震源發生時間對目標函數進行傳感器權重系數配置、系數矩陣變換求解以及迭代，直至滿足預設條件，包括：基于所述初始震源位置和所述傳感器的位置坐標，計算每個所述傳感器的地震波走時，并基于所述初始震源發生時間、微震到達時間與所述傳感器的地震波走時確定每個所述傳感器的地震波走時殘差；基于所述地震波信噪比、所述地震波震源距和所述地震波走時殘差計算所述目標函數的傳感器權重系數；基于所述地震波走時計算所述目標函數的系數矩陣，并對所述系數矩陣進行參數分離和歸一化處理，得到目標系數矩陣；基于所述地震波走時殘差、所述傳感器權重系數和所述目標系數矩陣確定所述目標函數，并基于所述目標函數計算得到修正的震源位置和震源發生時間；將修正后的所述震源位置和所述震源發生時間作為震源位置和震源發生時間，重復對所述目標函數進行迭代，直至滿足預設條件；

基于所述地震波信噪比、所述地震波震源距和所述地震波走時殘差計算所述目標函數的傳感器權重系數，包括：計算所述地震波信噪比，并基于所述地震波信噪比計算第一傳感器權重系數；計算所述地震波震源距，并基于所述地震波震源距計算第二傳感器權重系數；計算所述地震波走時殘差的平均值、中位數和標準差，并基于所述地震波走時殘差的平均值、中位數和標準差確定損失函數，以及基于所述損失函數計算第三傳感器權重系數；將所述第一傳感器權重系數、所述第二傳感器權重系數和所述第三傳感器權重系數的組合，確定為所述目標函數的傳感器權重系數。

2.根據權利要求1所述的方法，其特征在于，對所述系數矩陣進行參數分離和歸一化處理，得到目標系數矩陣，包括：

照以下公式對所述系數矩陣進行參數分離和歸一化處理：

其中，A為系數矩陣，，T_i為第i個傳感器的地震波走時，是傳感器i的走時對未知數的偏導，=(x,y,z)；B為進行參數分離后的系數矩陣，C為歸一化處理后的目標系數矩陣。

3.根據權利要求1所述的方法，其特征在于，將滿足所述預設條件時所述目標函數的解確定為震源位置和震源發生時間，包括：

若所述地震波走時殘差小于第一預設值，則將所述目標函數為最小值時的解確定為震源位置和震源發生時間。

4.根據權利要求1所述的方法，其特征在于，將滿足所述預設條件時所述目標函數的解確定為震源位置和震源發生時間，包括：

若迭代次數等于預設次數，則將歷次迭代結果中使得所述目標函數為最小值時的解確定為震源位置和震源發生時間。