共轭梯度法向残差法的重力异常向下延拓研究  

Downward continuation of gravity anomaly by conjugate gradient normal residual method

作　　者：赵予菲 柯宝贵[2] ZHAOYufei;KE Baogui(Faculty of Geomatics,Lanzhou Jiaotong University,Lanzhou 730070,China;Chinese Academy of Surveying and Mapping,Beijing 100036,China;National-Local Joint Engineering Research Center of Technologies and Applications for National Geographic State Monitoring,Lanzhou 730070,China;Gansu Provincial Engineering Laboratory for National Geographic State Monitoring,Lanzhou 730070,China)

机构地区：[1]兰州交通大学,测绘与地理信息学院,兰州730070 [2]中国测绘科学研究院,北京100036 [3]地理国情监测技术应用国家地方联合工程研究中心,兰州730070 [4]甘肃省地理国情监测工程实验室,兰州730070

出　　处：《测绘科学》2022年第10期29-38,共10页Science of Surveying and Mapping

基　　金：国家自然科学基金项目(41874022,41930535);青岛海洋科学与技术试点国家实验室“问海计划”项目(2021WHZZB1005);兰州交通大学优秀平台支持项目(201806)

摘　　要：针对重力异常向下延拓的不适定性导致求解难度大的问题,该文采用了共轭梯度法向残差(CGNR)法与频率域积分迭代法,使用模拟重力异常数据进行向下延拓,并对两种方法的延拓效果进行对比。实验结果表明:二者均需对变量的数据进行扩边,以更好地消除边缘效应;CGNR法结合快速傅里叶变换(FFT)可以提高计算效率,有利于抑制噪声的影响,算法稳健性强且相比频率域积分迭代法具有较高的延拓精度,向下延拓深度可达20倍点距,且保持良好的延拓精度。In view of the difficulty of solving the problem caused by the ill-posed of downward continuation of gravity anomaly,conjugate gradient normal residual(CGNR)method and frequency domain integral iterative method were adopted in this paper to downward continuation of simulated gravity anomaly data,and the continuation effects of the two methods were compared.The experimental results showed that both of them needed to expand the edge of variable data to better eliminate the edge effect.The CGNR method combined with fast Fourier transform(FFT)could improve the calculation efficiency and suppressed the influence of noise.Compared with frequency domain integral iterative method,the CGNR algorithm had strong robustness and higher continuation accuracy,and the downward continuation depth could reach 20 times the point distance,and maintained good continuation accuracy.

关 键 词：向下延拓 频率域 积分迭代法 共轭梯度法向残差法 

分 类 号：P223[天文地球—大地测量学与测量工程]