Abstract:ABSTRACTComplex topography, the free-surface boundary condition, and anelastic properties of media should be accounted for in the frame of onshore geophysical prospecting imaging, such as full-waveform inversion (FWI). In this context, an accurate and efficient forward-modeling engine is mandatory. We have performed 3D frequency-domain anisotropic elastic wave modeling by using the highly accurate spectral element method and a sparse multifrontal direct solver. An efficient approach similar to computing the matrix-vector product in the time domain is used to build the matrix. We validate the numerical results by comparing with analytical solutions. A parallel direct solver, the sparse direct multifrontal massively parallel solver (MUMPS), is used to solve the linear system. We find that a hybrid implementation of message passing interface and open multiprocessing is more efficient in flops and memory cost. The influence of the deformed mesh, free-surface boundary condition, and heterogeneity of media on MUMPS performance is negligible. Complexity analysis suggests that the memory complexity of MUMPS agrees with the theoretical order [Formula: see text] (or [Formula: see text] with an efficient matrix reordering method) for an [Formula: see text] grid when nontrivial topography is considered. With the available resources, we conduct a moderate scale modeling with a subset of the SEAM Phase II Foothills model, where 60 wavelengths in the [Formula: see text]-axis are propagated. Computing one gradient of FWI based on this model using the frequency-domain modeling is shown to require similar or fewer computational resources than what would be required for a time-domain solver, depending on the number of sources, while larger memory is necessary. An estimation of the increasing trend indicates that approximately 20 Tb of memory would be required for a [Formula: see text] wavelength modeling. The limit of MUMPS scalability hinders the application to larger scale applications.

Three-dimensional preconditioned FM-IBEM solution to broadband-frequency seismic wave scattering in a layered sedimentary basin

A Refinement and High-Degree Stable Marching-on-in-Degree Method for Wide-Band Scattering from Perfect Conductors

Three-dimensional Broadband Simulation of Near-Fault Seismic Ground Motion Considering Complete Source-Path-site Effect: Modeled by Fast Multipole Indirect Boundary Element Method

Application of Multilevel Inverse-Based ILU Preconditioning Approach to Hybrid FE/BI/MLFMA for 3D Scattering

A Three-dimensional Immersed Boundary Method for Accurate Simulation of Acoustic Wavefields with Complex Surface Topography

1-D Finite Element Artificial Boundary Method for Saturated Layered Half-Space Site Response under Obliquely Incident P-SV Waves

Broadband seismic isolation effect of three-dimensional seismic metamaterials in a semi-infinite foundation: modelled by fast multipole indirect boundary element method

3D Frequency-Domain Elastic Wave Modeling with the Spectral Element Method Using a Massively Parallel Direct Solver

A Fast Multipole Boundary Element Method for Three-Dimensional Half-Space Acoustic Wave Problems over an Impedance Plane

A Fast Multipole Boundary Element Method for 3D Multi-Domain Acoustic Scattering Problems Based on the Burton–Miller Formulation

Three-dimensional parallel hybrid PSM/FDM simulation of seismic wave propagation

Fem/wideband fmbem coupling analysis for two dimensional acoustic fluid-structure interaction problems

Seismic wave-field modeling of 3D irregularities in layered half-space under P-, S-, and Rayleigh waves with arbitrary incident azimuth

A Hierarchical Preconditioner For Fe-Bi Matrix Equations Of Three-Dimensional Scattering

An Effective Domain-Decomposition-Based Preconditioner for the FE-BI-MLFMA Method for 3D Scattering Problems

Multilevel Fast Multipole Boundary Element Method for 3D Acoustic Problems and Its Applications

An efficient approach for computing scattering by inhomogeneous objects with thin dielectric structures

3-D Marine CSEM Forward Modeling in Anisotropic Media Using a Space-Wavenumber Domain Integral Equation Method

Preconditioned Multilevel Fast Inhomogeneous Plane Wave Algorithm for Solving Electromagnetic Scattering Problems

Three-dimensional element-by-element parallel spectral-element method for seismic wave modeling

Mixed Spectral-Element Method for Overcoming the Low-Frequency Breakdown Problem in Subsurface EM Exploration