互相关系数自约束的重力三维反演与高效求解

作者：梁生贤

来源：吉林大学学报(地球科学版), 2018, 48(05): 1473-1482.

DOI：10.13278/j.cnki.jjuese.20170167

• 重力三维反演
• 互相关系数
• 自约束反演
• 阻尼LSQR算法
• 云南芦子园铁铅锌铜多金属矿床
• 花岗岩体

本文将拟合残差计算所得的互相关系数作为先验信息,与深度加权函数同时引入到重力正则化反演的模型约束中,以提升反演结果的可靠性。针对三维反演中的大型线性方程组问题,引入阻尼LSQR(最小二乘QR分解)算法,结合等效几何格架技术,将大矩阵按照模型单元划分为若干个子矩阵进行存储与运算。理论模型计算结果表明:同时利用互相关系数和深度加权的模型自约束反演,能较清晰地反映真实异常体;基于分块矩阵的阻尼LSQR算法求解线性方程组较直接法节省了几千甚至上万倍的存储量,且计算速率提高了数倍,可在普通计算机上实现较大规模的反演计算。将其应用于云南芦子园铁铅锌铜多金属矿床隐伏花岗岩体定位,取得了良好的效果。

• 单位
  中国地质调查局成都地质调查中心