位场向下延拓系数矩阵性质及Barzilai-Borwein向下延拓法  

Coefficient Matrix Properties of Downward Continuation for Potential Fields and Barzilai-Borwein Downward Continuation Method

作　　者：张志厚 廖晓龙 姚禹 范祥泰 路润琪 ZHANG Zhihou;LIAO Xiaolong;YAO Yu;FAN Xiangtai;LU Runqi(Faculty of Geosciences and Environmental Engineering,Southwest Jiaotong University,Chengdu 611756,China)

机构地区：[1]西南交通大学地球科学与环境工程学院,四川成都611756

出　　处：《西南交通大学学报》2021年第2期323-330,362,共9页Journal of Southwest Jiaotong University

基　　金：国家自然科学基金(41672295);四川省科技计划(2019YFG0460,2019YFG0001,2019YFG0047,21YYJC3115);中国中铁股份有限公司科技研究开发计划项目(CZ01-重点-05)。

摘　　要：位场的向下延拓不仅仅能够提高地球物理数据解释的可靠性,在导航方面也有着重要的作用.为了进一步提高计算精度和速度,提出了位场向下延拓的Barzilai-Borwein(BB)法.首先证明了位场向下延拓的系数矩阵为对称的双重Toeplitz系统矩阵(block-Toeplitz-Toeplitz-block,BTTB);其次,假定该系数矩阵为正定的条件下,采用BB法迭代求解下延方程组,并约束其迭代步长确保算法收敛;最后,分别通过理论模型无噪声数据和实际资料对BB法进行检验,并与积分迭代法进行对比.结果表明:理论模型验证时,同一收敛精度条件下,BB法的计算速度是积分迭代法的2倍以上;实际资料检验时,在相同计算次数下,BB法与积分迭代法的平均相对误差分别为6.1%与7.7%.Downward continuation of potential fields can not only improve the reliability of geophysical data interpretation,but also plays an important role in navigation.In order to further enhance the calculation accuracy and speed of downward continuation,the Barzilai-Borwein(BB)downward continuation method is proposed.First,the coefficient matrix of downward continuation for potential fields is proved to be a symmetric block-Toeplitz-Toeplitz-block matrix(BTTB).Next,assuming that the coefficient matrix is positive definite,the BB method is used to solve the equations of downward continuation,and the iterative step is limited to ensure the convergence of the algorithm.Finally,the BB method is validated by noise-free data of the theoretical model and practical data,and compared with the integral iterative method.The results show that under the same convergence accuracy,in the case of the theoretical model,the calculation speed of BB method is more than 2 times that of integral iteration method.In the case of practical data,the average relative errors of BB method and integral iteration method are 6.1%and 7.7%,respectively.

关 键 词：向下延拓 系数矩阵 Barzilai-Borwein法 收敛性 

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