[發明專利]一種基于積分方程的退磁效應下磁測數據的正反演方法

有效

• 申請號: 202110443165.6
• 申請日: 2021-04-23
• 公開（公告）號: CN113139289B
• 公開（公告）日: 2022-06-07
• 發明（設計）人: 左博新
• 申請（專利權）人: 中國地質大學（武漢）
• 主分類號: G06F30/20
• 分類號: G06F30/20;G06F17/11
• 代理公司: 武漢知產時代知識產權代理有限公司 42238
• 代理人: 曹雄

1.一種基于積分方程的退磁效應下磁測數據的正反演方法，其特征在于：

S1、根據地下目標區域，確定網格剖分的空間范圍，使用正交網格對所述空間范圍進行均勻網格剖分，得到剖分后的網格；

S2、設定地磁場的磁傾角、磁偏角和磁感應強度；

S3、在剖分后的網格中設定觀測位置，利用基于積分方程的退磁效應下的磁測數據正演計算方法，獲得在退磁效應下觀測位置的磁測數據；

S4、根據磁測數據，通過基于積分方程的退磁效應下磁測數據的反演方法，結合正則化項進行模型最優化計算，對地下目標進行反演計算，得到反演結果；

步驟S3中，基于積分方程的退磁效應下的磁測數據正演計算方法，其計算公式為：

M＝K(H₀+TM)

AM＝b

A＝I-KT

b＝KH₀

B_ij＝μ₀T_ijM,i∈{x,y,z},j∈{x,y,z}

其中，M表示考慮退磁效應后的磁化強度矩陣；K表示磁化率矩陣，H₀表示磁場強度矢量,B_ij異常場的磁梯度張量數據,μ₀是真空中的磁導率；I為單位陣；p的觀測點位置位于每個正交網格單元的中心；

步驟S4中，基于積分方程的退磁效應下磁測數據的反演方法的目標函數為：

F(κ)＝Gκ

G＝T_ij(B₀+μ₀TM)

其中，F(κ)表示正演過程計算出的磁測數據，d^obs表示實際觀測磁測數據,β代表規整化因子，Δκ表示模型規整化項，G為靈敏度矩陣，κ表示模型磁化率矩陣，B₀為地磁場強度，μ₀表示真空磁導率。