物探化探计算技术2017，Vol.39Issue(2)：176-182,7.DOI:10.3969/j.issn.1001-1749.2017.02.04
任意密度分布复杂地质体重力异常快速三维正演方法
Three-dimensional forward modeling method for gravity anomalies of complex geological bodies with arbitrary density distribution
摘要
Abstract
Rapid three-dimensional (3D) forward modeling of gravity anomalies of complex geologic bodies with arbitrary density distribution plays an important role on 3D inversion of gravity and human-computer interaction modeling and interpretation.This paper proposes a wavenumber domain 3D forward modeling method.This new method has three key aspects: (1) prism based subdivision strategy.The research area is divided into many small regular prisms, density for each small prism can be given arbitrarily, in order to simulate complex geological body with arbitrary density distribution and undulating topography;(2) new wavenumber domain formula for computing gravity anomaly of combined prism model;(3) Gauss-FFT method.This paper uses Gauss-FFT method to transform the gravity anomaly from the wave number domain to the spatial domain.The Gauss-FFT method not only has high computational efficiency, but also suppresses the boundary effect and avoids the singular problem at zero wavenumber caused by the traditional FFT method.Numerical tests show that the proposed method is fast and accurate.It takes only a few seconds to calculate gravity anomaly for the combined 3D model with millions of prisms.关键词
重力异常/三维正演/Gauss-FFT/积分方程Key words
gravity anomalies/three-dimensional forward modeling/Gauss-FFT/integral equation分类
天文与地球科学引用本文复制引用
陈轻蕊,陈龙伟,戴世坤,张钱江,欧阳芳..任意密度分布复杂地质体重力异常快速三维正演方法[J].物探化探计算技术,2017,39(2):176-182,7.基金项目
国家自然科学基金项目(41404106) （41404106）
