[发明专利]一种重力多界面反演方法

摘要

本发明公开了一种重力多界面反演方法，该方法包括以下步骤：获取观测点的信息、重力观测值、初始模型范围、地下介质信息、待反演层位信息、及模型泛函各部分的正则化权重；对地下介质进行网格剖分；计算初始模型网格剖分后的各柱状体对观测点的重力效应并进行累加，得到重力正演值；根据重力观测值与重力正演值，计算失配泛函F₁(x)；计算全变差函数F₂(x)；计算模型界面反演偏差函数F₃(x)；根据计算结果、模型泛函各部分的正则化权重λ₂～λ₃与失配泛函F₁(x)构成目标函数：F(x)＝F₁(x)+λ₂F₂(x)+λ₃F₃(x)；采用迭代方法使目标函数最小。通过上述方式，本发明能够具有较好的反演精度和效率。

主权项

一种重力多界面反演方法，其特征在于，该方法包括以下步骤：(1)获取观测点的信息、重力观测值、初始模型范围、地下介质信息、待反演层位信息、及模型泛函各部分的正则化权重；(2)对地下介质进行网格剖分；(3)计算初始模型网格剖分后的各柱状体对观测点的重力效应并进行累加，得到重力正演值；(4)根据所述重力观测值与重力正演值，计算失配泛函F₁(x)，其中，x为待求的各柱状体底界面深度；(5)计算全变差函数F₂(x)；(6)计算模型界面反演偏差函数F₃(x)；(7)根据步骤(5)～(6)的计算结果构成模型泛函，并根据所述模型泛函各部分的正则化权重λ₂～λ₃与所述失配泛函F₁(x)一起构成目标函数：F(x)＝F₁(x)+λ₂F₂(x)+λ₃F₃(x)；(8)采用迭代方法使所述目标函数最小，期间根据迭代的结果与给定的不等约束，不断修改模型以使反演结果趋向真实值。