Efficient 3D large-scale forward modeling of gravity anomaly in space-wavenumber mixing domain

• DOI: 10.1190/gem2019-018.1
• Publication Date: 2019-09-11T15:25:26Z
• Authors: Dongdong Zhao; Xiaoying Hu; Shikun Dai*; Deqiang Tao; Yinming Zhou; Kun Li; Xulong Wang
• Keywords: Wavenumber; Anomaly (physics)

A novel 3D numerical modeling method with topography for gravity anomaly in space-wavenumber mixing domain is proposed to reduce the computation time and memory requirements. By performing a two-dimensional Fourier transform along horizontal directions, three-dimensional partial differential equations spatial are transformed into group of independent one-dimensional engaged different wave numbers. method, highly parallel among The preserves vertical component space domain, so mesh can be fine at shallow depths coarse deep depth. finite element used solve numbers, efficiency solving fixed-bandwidth linear further improved by chasing method. In synthetic test, using 1 processor, more than 40,401 observation points 4,040,100 source elements computed 18.9 s. traditional Gauss-FFT requires 81.4 Hence, our new algorithm substantially accelerate solution inversion problems.