[發(fā)明專利]一種基于積分方程的退磁效應(yīng)下磁測(cè)數(shù)據(jù)的正反演方法

有效

• 申請(qǐng)?zhí)枺?/td>: 202110443165.6
• 申請(qǐng)日: 2021-04-23
• 公開(kāi)（公告）號(hào): CN113139289B
• 公開(kāi)（公告）日: 2022-06-07
• 發(fā)明（設(shè)計(jì)）人: 左博新
• 申請(qǐng)（專利權(quán)）人: 中國(guó)地質(zhì)大學(xué)（武漢）
• 主分類號(hào): G06F30/20
• 分類號(hào): G06F30/20;G06F17/11
• 代理公司: 武漢知產(chǎn)時(shí)代知識(shí)產(chǎn)權(quán)代理有限公司 42238
• 代理人: 曹雄

1.一種基于積分方程的退磁效應(yīng)下磁測(cè)數(shù)據(jù)的正反演方法，其特征在于：

S1、根據(jù)地下目標(biāo)區(qū)域，確定網(wǎng)格剖分的空間范圍，使用正交網(wǎng)格對(duì)所述空間范圍進(jìn)行均勻網(wǎng)格剖分，得到剖分后的網(wǎng)格；

S2、設(shè)定地磁場(chǎng)的磁傾角、磁偏角和磁感應(yīng)強(qiáng)度；

S3、在剖分后的網(wǎng)格中設(shè)定觀測(cè)位置，利用基于積分方程的退磁效應(yīng)下的磁測(cè)數(shù)據(jù)正演計(jì)算方法，獲得在退磁效應(yīng)下觀測(cè)位置的磁測(cè)數(shù)據(jù)；

S4、根據(jù)磁測(cè)數(shù)據(jù)，通過(guò)基于積分方程的退磁效應(yīng)下磁測(cè)數(shù)據(jù)的反演方法，結(jié)合正則化項(xiàng)進(jìn)行模型最優(yōu)化計(jì)算，對(duì)地下目標(biāo)進(jìn)行反演計(jì)算，得到反演結(jié)果；

步驟S3中，基于積分方程的退磁效應(yīng)下的磁測(cè)數(shù)據(jù)正演計(jì)算方法，其計(jì)算公式為：

M＝K(H₀+TM)

AM＝b

A＝I-KT

b＝KH₀

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

其中，M表示考慮退磁效應(yīng)后的磁化強(qiáng)度矩陣；K表示磁化率矩陣，H₀表示磁場(chǎng)強(qiáng)度矢量,B_ij異常場(chǎng)的磁梯度張量數(shù)據(jù),μ₀是真空中的磁導(dǎo)率；I為單位陣；p的觀測(cè)點(diǎn)位置位于每個(gè)正交網(wǎng)格單元的中心；

步驟S4中，基于積分方程的退磁效應(yīng)下磁測(cè)數(shù)據(jù)的反演方法的目標(biāo)函數(shù)為：

F(κ)＝Gκ

G＝T_ij(B₀+μ₀TM)

其中，F(xiàn)(κ)表示正演過(guò)程計(jì)算出的磁測(cè)數(shù)據(jù)，d^obs表示實(shí)際觀測(cè)磁測(cè)數(shù)據(jù),β代表規(guī)整化因子，Δκ表示模型規(guī)整化項(xiàng)，G為靈敏度矩陣，κ表示模型磁化率矩陣，B₀為地磁場(chǎng)強(qiáng)度，μ₀表示真空磁導(dǎo)率。