Journal article
Fast 3D gravity and magnetic modelling using midpoint quadrature and 2D FFT
Scientific reports, Vol.13(1), pp.9304-9
06/08/2023
DOI: 10.1038/s41598-023-36525-2
PMCID: PMC10250347
PMID: 37291183
Abstract
To avoid the problem of the traditional methods consuming large computational resources to calculate the kernel matrix and 2D discrete convolution, we present a novel approach for 3D gravity and magnetic modelling. This method combines the midpoint quadrature method with a 2D fast Fourier transform (FFT) to calculate the gravity and magnetic anomalies with arbitrary density or magnetic susceptibility distribution. In this scheme, we apply the midpoint quadrature method to calculate the volume element of the integral. Then, the convolution of the weight coefficient matrix with density or magnetization is efficiently computed via the 2D FFT. Finally, the accuracy and efficiency of the proposed algorithm are validated by using an artificial model and a real topography model. The numerical results demonstrate that the proposed algorithm's computation time and the memory requirement are decreased by approximately two orders of magnitude compared with the space-wavenumber domain method.
Details
- Title: Subtitle
- Fast 3D gravity and magnetic modelling using midpoint quadrature and 2D FFT
- Creators
- Xulong Wang - Central South UniversityJianxin Liu - Central South UniversityJian Li - Central South UniversityHang Chen - Boise State University
- Resource Type
- Journal article
- Publication Details
- Scientific reports, Vol.13(1), pp.9304-9
- DOI
- 10.1038/s41598-023-36525-2
- PMID
- 37291183
- PMCID
- PMC10250347
- NLM abbreviation
- Sci Rep
- ISSN
- 2045-2322
- eISSN
- 2045-2322
- Publisher
- NATURE PORTFOLIO
- Number of pages
- 9
- Language
- English
- Date published
- 06/08/2023
- Academic Unit
- Earth and Environmental Sciences
- Record Identifier
- 9984962532102771
Metrics
4 Record Views