一种重力梯度的自适应剖分反演方法与流程

技术特征：
1.一种重力梯度的自适应剖分反演方法，其特征在于，包括以下步骤：步骤1：获得重力梯度异常数据或正演出模型的重力梯度异常数据；步骤2：采用倾斜角边界识别法，在结果大于等于0的范围内进行四面体剖分；步骤3：构建非全局核函数矩阵，并实现非全局四面体剖分反演计算，得到空间密度分布图。2.根据权利要求1所述的重力梯度的自适应剖分反演方法，其特征在于，所述倾斜角边界识别法，公式如下：式中：其中g
z
为重力梯度异常。3.根据权利要求1所述的重力梯度的自适应剖分反演方法，其特征在于，所述构建非全局核函数矩阵的过程具体为：将任意多面体重力异常表达式应用于所述四面体，任意观测点(x
n
,y
n
,z
n
)处的重力异常计算公式如下：重力梯度为重力异常的一阶垂直导：其中：g为核函数矩阵；γ为万有引力常数6.672
×
10
‑
11
m3/(kg
·
s2)；ρ为四面体密度；φ
i
为四面体第i个三角形与z轴的夹角；j
j,i
由如下公式得到：其中：i为三角形序号，j为第i个三角形中的边序号；其中：i为三角形序号，j为第i个三角形中的边序号；其中：i为三角形序号，j为第i个三角形中的边序号；θ
i
为四面体第i个面外法向在水平方向的投影与x方向的夹角，式中0≤θ
i
≤2π，0≤φ
i
≤π。
4.根据权利要求1所述的重力梯度的自适应剖分反演方法，其特征在于，在实现非全局四面体剖分反演计算时，所求解的方程解为欠定方程的解，在求解欠定方程时，构造下式目标函数：m是密度参数的向量；m
ref
是参考模型，一般取值为0；λ为正则化参数，起到平衡数据拟合函数与权函数的作用；w
d
是数据加权矩阵；为深度权，用于约束反演结果使结果更加符合实际克服趋肤效应，且我们在四面体剖分反演中加上与体积有关的权：v代表四面体的体积，z0是由深度权函数与观测点下的核函数进行匹配得到，α是一个经验值，通常取0.5，β是一个经验值。

技术总结
本发明公开了一种重力梯度的自适应剖分反演方法，包括如下过程：获得重力梯度异常数据或正演出模型的重力梯度异常数据；采用倾斜角边界识别法，在结果大于等于0的范围内进行四面体剖分；构建非全局核函数矩阵，并实现非全局四面体剖分反演计算，得到空间密度分布图。本发明通过边界识别圈定出目标矿体的范围，在该范围内进行网格剖分；并采用非全局四面体剖分来拟合起伏地形和不规则地质体，更加贴合地形形状；相比传统的物性反演方法，本发明能够实现在同样的区域内通过缩小剖分范围来减少剖分网格，以达到加快反演速度，提高反演精度的目的，同时降低对计算机性能的要求。同时降低对计算机性能的要求。同时降低对计算机性能的要求。


技术研发人员：牛润馨 马国庆 李丽丽 高桐
受保护的技术使用者：吉林大学
技术研发日：2021.07.12
技术公布日：2021/10/19