摘要
实用化的重力数据三维反演过程中需要大量的计算时间和存储空间,针对此问题本文提出一种组合优化算法:(1)利用空间域的精确性和波数域的快速性,提出一种高精度快速三维正演算法;(2)利用反演矩阵中的对称性,将重力共轭梯度反演(CG)中的主要计算分解成两次正演计算,从而优化了三维重力反演的计算效率。通过不同网格剖分数量模型的数值实验,验证了本文优化算法的计算精度和计算效率。
The practical application of 3D inversion of gravity data requires a lot of computation time and storage space.To solve this problem,we present an integrated optimization algorithm with the following components:(1)targeting high accuracy in the space domain and fast computation in the wavenumber domain,we design a fast 3D forward algorithm with high precision;and(2)taking advantage of the symmetry of the inversion matrix,the main calculation in gravity conjugate gradient inversion is decomposed into two forward calculations,thus optimizing the computational efficiency of 3D gravity inversion.We verify the calculation accuracy and efficiency of the optimization algorithm by testing various grid-number models through numerical simulation experiments.
作者
荆磊
姚长利
杨亚斌
徐梦龙
张光之
纪若野
Jing Lei;Yao Chang-Li;Yang Ya-Bin;Xu Meng-Long;Zhang Guang-Zhi;Ji Ruo-Ye(School of geophysics and information technology,China University of Geosciences(Beijing),Beijing 100083,China;The National Center for Geological Exploration Technology,Langfang 065000,China;Institute of Geophysical and Geochemical Exploration,Chinese Academy of Geological Sciences,Langfang 065000,China.)
基金
Financial support by the China Geological Survey Project(Nos.DD20190030,DD20190032)
关键词
重力
三维反演
优化算法
高效率
海量数据
Gravity
3D inversion
optimization algorithm
high efficiency
massive data