重力梯度全张量数据三维共轭梯度聚焦反演  

The 3D focusing inversion of full tensor gravity gradient data based on conjugate gradient

作　　者：陈闫[1] 李桐林[1] 范翠松[2] 许博文 孙业发 

机构地区：[1]吉林大学地球探测科学与技术学院,长春130026 [2]中国地质调查局天津地质调查中心,天津300170 [3]山东省第一地质矿产勘察院,济南250014

出　　处：《地球物理学进展》2014年第3期1133-1142,共10页Progress in Geophysics

基　　金：国家重大科学仪器设备开发专项"大深度三维电磁探测技术工程化开发"(2011YQ05006009);吉黑东部综合找矿方法研究(201301-060-004)联合资助

摘　　要：随着全张量重力梯度测量技术的日趋成熟和应用领域的不断扩大,重力梯度全张量数据的三维反演技术越来越受重视.本文利用剖分单元之间几何架构等效性,实现了重力梯度全张量场三维正演快速计算和导数矩阵优化存储.并将积分灵敏度、粗糙度和最小支撑泛函约束以及参考模型和模型参数界限约束引入到目标函数中,采用共轭梯度法进行反演迭代,实现了重力梯度全张量数据三维快速正反演计算.多种模型的反演试算表明:本文提出的反演算法的可靠性和稳定性较好,并且算法速度快、占用内存低且易于并行化.As the full tensor gravity gradient measurement technology matures and its applications continues to expand,the 3Dinversion of full tensor gravity gradient data is paid more and more attention.Using the split unit between geometric structures equivalence, this paper realizes the threedimensional fast forward calculation of gravity gradient full tensor field and its optimization of derivative matrix storage.Taking integral sensitivity,roughness, minimum support functional constraints and reference model,model parameter range constraints into the objective function,and using the conjugate gradient method for inversion iteration,realizing the full tensor gravity gradient data 3Drapid inversion calculation.A variety of the model inversion trial indicates that the proposed has good reliability and stability of the inversion algorithm,and has high speed,low memory and easy to parallelize.

关 键 词：重力梯度全张量 三维反演 共轭梯度 聚焦反演 多重模型约束 

分 类 号：P631[天文地球—地质矿产勘探]