重力归一化总梯度的稳定解法

重力归一化总梯度法是20世纪80年代兴起的一种求解重力反演的方法,其关键运算是向下延拓及求导。原有的向下延拓方法大多数是在波数域实现的,由于稳定性差,影响了归一化总梯度法的效果,从而难以推广应用。本文以原先的希尔伯特变换法求解重力归一化总梯度值为基础,采用了一种新的、稳定的向下延拓算法——迭代法,替代了原有波数域的向下延拓算子,提高了计算的精度和稳定性。模型实验和实测资料应用显示效果明显改善。

基于向下延拓Milne法的重力归一化总梯度法

重力归一化总梯度法需进行向下延拓和垂向导数的计算,这两个过程都会放大高波数成分,影响计算结果的稳定性。为此,将向下延拓Milne法及积分垂向二阶导数法引入重力归一化总梯度的计算。将该方法应用于估算无限长水平圆柱体模型的中心位置,计算结果准确。然后针对实测重力数据,利用本方法得到的矿洞中心埋深较基于泰勒级数展开的重力归一化总梯度法更准确,因此本文方法具有一定的实用性。