论文部分内容阅读
重力异常界面反演是重力勘探中的一个重要研究课题,它是一个求解第一类Fredholm非线性积分的问题。目前情况下直接求解这类问题比较困难,一般采用拟合法或逐步线性化办法实现反演计算。而在计算过程中,由于反演的不适定性,往往不能很好地获得反演解。
本文先从二度体重力异常正演公式出发,研究了重力异常密度界面反演问题,导出地下界面迭代反演的计算公式;然后利用广义似然函数对模型参数进行优化选择,从而获得模型参数修正量迭代计算方法,实现二度体单-密度界面的稳健迭代(最小二乘加权迭代)反演。模型计算表明方法可以很好反演出地下界面。由于实际资料不能完全分离出正常场,为了消除正常场的残余量对密度界面反演的影响,又用重力异常水平梯度来反演密度界面。先取重力异常的水平梯度,然后利用广义似然函数对模型参数进行优化选择,从而得到地下模型参数修正量的迭代计算公式,将LSQR方法引入到该方法中来,用它求解大型方程组,得到模型参数的修正量,从而实现了二度体单-密度界面的反演,模型计算表明此方法是行之有效的。
最后讨论了正常场对重力异常密度界面反演和重力异常水平梯度密度界面反演的影响。