[發明專利]三維重力快速反演優化方法、系統、存儲介質和電子設備

在審

• 申請號: 202210166718.2
• 申請日: 2022-02-23
• 公開（公告）號: CN114611062A
• 公開（公告）日: 2022-06-10
• 發明（設計）人: 荊磊;杜炳銳;徐夢龍;張光之;施蘇利
• 申請（專利權）人: 中國地質科學院地球物理地球化學勘查研究所
• 主分類號: G06F17/12
• 分類號: G06F17/12;G06F17/15;G06F17/14;G06F17/16;G06F30/23
• 代理公司: 北京輕創知識產權代理有限公司 11212
• 代理人: 孟仕杰

1.一種三維重力快速反演優化方法，其特征在于，包括：

計算三維模型中每個深度層的正演核矩陣并將每兩個正演核矩陣進行組合，得到多個組合正演核矩陣

對每個組合正演核矩陣進行快速傅里葉變換，得到所述每個組合正演核矩陣對應的組合正演核矩陣頻譜

基于所述正演核矩陣的大小，對二維矩陣依次進行加零擴展與快速傅里葉變換，得到二維矩陣頻譜

依次將所述每個組合正演核矩陣頻譜與所述二維矩陣頻譜代入預設復合頻譜計算公式，獲取并根據每兩個深度層的和得到G^Tu，其中，G^Tu為三維重力反演過程的中間參數，k為所述三維模型的第k深度層，l為所述三維模型的第l深度層，j²＝-1。

2.根據權利要求1所述的三維重力快速反演優化方法，其特征在于，所述預設復合頻譜公式為：

其中，為所述三維模型中第k深度層的正演核矩陣，為所述三維模型中第l深度層的正演核矩陣，(1≤k,l≤p)，p為所述三維模型中的總深度層數，F表示快速傅里葉變換的正變換，F^-1表示快速傅里葉變換的逆變換，u為所述單層模型，j²＝-1，m為所述三維模型中x(北)方向的網格剖分數量，n為所述三維模型中y(東)方向的網格剖分數量，H為復合頻譜信號。

3.根據權利要求2所述的三維重力快速反演優化方法，其特征在于，所述每兩個深度層的和的計算過程為：

其中，Imag[H]為所述復合頻譜信號的虛部，Real[H]為所述復合頻譜信號的實部。

4.根據權利要求3所述的三維重力快速反演優化方法，其特征在于，所述G^Tu的具體形式為：

5.根據權利要求1-4任一項所述的三維重力快速反演優化方法，其特征在于，在所述對每個組合正演核矩陣進行快速傅里葉變換，得到所述每個組合正演核矩陣對應的組合正演核矩陣頻譜之后，還包括：

對所述每個組合正演核矩陣頻譜進行共軛對稱變換，得到所述每個組合正演核矩陣頻譜對應的變換組合正演核矩陣頻譜

基于所述每個深度層的正演核矩陣的大小，依次對所述每個深度層的密度矩陣進行加零擴展，得到每個密度矩陣對應的擴展密度矩陣，并基于所述每兩個正演核矩陣的組合方式，對每兩個擴展模型矩陣進行組合并經過快速傅里葉變換，得到多個組合密度矩陣頻譜

依次累加每兩個相同深度層組合的變換組合正演核矩陣頻譜與對應的組合密度矩陣頻譜的乘積，得到第一計算結果

對所述第一計算結果進行共軛對稱變換，得到第二計算結果并將所述第一計算結果和所述第二計算結果代入預設正演頻譜公式，得到并對所述三維模型的每個深度層正演頻譜之和進行反傅里葉變換，得到所述三維模型的正演結果。

6.根據權利要求5所述的三維重力快速反演優化方法，其特征在于，所述預設正演頻譜公式為：

其中，為所述三維模型的每個深度層正演頻譜之和，k為所述三維模型的第k深度層，l為所述三維模型的第l深度層，j²＝-1，(1≤k,l≤p)，p為所述三維模型中的總深度層數，F表示快速傅里葉變換的正變換，m為所述三維模型中x(北)方向的網格剖分數量，n為所述三維模型中y(東)方向的網格剖分數量。