[發(fā)明專利]弱化正則化參數(shù)的三維重力反演方法、系統(tǒng)和存儲(chǔ)介質(zhì)

在審

• 申請(qǐng)?zhí)枺?/td>: 202210188458.9
• 申請(qǐng)日: 2022-02-28
• 公開(kāi)（公告）號(hào): CN114662052A
• 公開(kāi)（公告）日: 2022-06-24
• 發(fā)明（設(shè)計(jì)）人: 荊磊;徐夢(mèng)龍;施蘇利;孫誠(chéng)業(yè);邱隆君
• 申請(qǐng)（專利權(quán)）人: 中國(guó)地質(zhì)科學(xué)院地球物理地球化學(xué)勘查研究所
• 主分類號(hào): G06F17/16
• 分類號(hào): G06F17/16;G06F17/18;G06F30/20;G06F111/04
• 代理公司: 北京輕創(chuàng)知識(shí)產(chǎn)權(quán)代理有限公司 11212
• 代理人: 趙秀斌

1.一種弱化正則化參數(shù)的三維重力反演方法，其特征在于，包括：

根據(jù)三維模型中的模型加權(quán)矩陣和所述模型加權(quán)矩陣的逆矩陣，生成用于弱化正則化參數(shù)的反演替換公式；

根據(jù)所述反演替換公式和所述三維模型的原始反演公式中，得到所述三維模型的第一反演公式；

根據(jù)所述第一反演公式對(duì)所述三維模型進(jìn)行三維重力反演。

2.根據(jù)權(quán)利要求1所述的弱化正則化參數(shù)的三維重力反演方法，其特征在于，所述反演替換公式為：W^-1W＝I，其中，W為所述模型加權(quán)矩陣，W^-1為所述模型加權(quán)矩陣的逆矩陣，I為單位矩陣。

3.根據(jù)權(quán)利要求2所述的弱化正則化參數(shù)的三維重力反演方法，其特征在于，所述原始反演公式為：ρ＝ρ₀+(G^TD^-1G+αW^-1)^-1G^TD^-1(d-Gρ₀)，所述第一反演公式為：ρ＝ρ₀+(WG^TD^-1G+αI)^-1WG^TD^-1(d-Gρ₀)，其中，ρ為所述三維模型的模型空間，ρ₀為參考模型，D^-1為對(duì)角數(shù)據(jù)協(xié)方差矩陣，α為所述正則化參數(shù)。

4.根據(jù)權(quán)利要求3所述的弱化正則化參數(shù)的三維重力反演方法，其特征在于，所述根據(jù)所述第一反演公式對(duì)所述三維模型進(jìn)行三維重力反演，具體包括：

當(dāng)所述第一反演公式中的正則化參數(shù)取值為零時(shí)，則所述第一反演公式演變?yōu)榈诙囱莨剑龅诙囱莨綖棣眩溅?Sub>0+(WG^TD^-1G)^-1WG^TD^-1(d-Gρ₀)；

根據(jù)所述第二反演公式對(duì)所述三維模型進(jìn)行三維重力反演。

5.根據(jù)權(quán)利要求1-4任一項(xiàng)所述的弱化正則化參數(shù)的三維重力反演方法，其特征在于，所述根據(jù)三維模型中的模型加權(quán)矩陣和所述模型加權(quán)矩陣的逆矩陣，生成用于弱化正則化參數(shù)的反演替換公式，具體包括：

當(dāng)所述模型加權(quán)矩陣的逆矩陣的對(duì)角元素不為零時(shí)，基于所述模型加權(quán)矩陣和所述模型加權(quán)矩陣的逆矩陣，生成用于弱化正則化參數(shù)的所述反演替換公式。

6.一種弱化正則化參數(shù)的三維重力反演系統(tǒng)，其特征在于，包括：生成模塊、處理模塊和運(yùn)行模塊；

所述生成模塊用于：根據(jù)三維模型中的模型加權(quán)矩陣和所述模型加權(quán)矩陣的逆矩陣，生成用于弱化正則化參數(shù)的反演替換公式；

所述處理模塊用于：根據(jù)所述反演替換公式和所述三維模型的原始反演公式中，得到所述三維模型的第一反演公式；

所述運(yùn)行模塊用于：根據(jù)所述第一反演公式對(duì)所述三維模型進(jìn)行三維重力反演。

7.根據(jù)權(quán)利要求6所述的弱化正則化參數(shù)的三維重力反演系統(tǒng)，其特征在于，所述反演替換公式為：W^-1W＝I，其中，W為所述模型加權(quán)矩陣，W^-1為所述模型加權(quán)矩陣的逆矩陣，I為單位矩陣。