Seismic wave propagation and inversion with Neural Operators

Abstract

Seismic wave propagation forms the basis for most aspects of seismological research, yet solving the wave equation is a major computational burden that inhibits the progress of research. This is exacerbated by the fact that new simulations must be performed when the velocity structure or source location is perturbed. Here, we explore a prototype framework for learning general solutions using a recently developed machine learning paradigm called Neural Operator. A trained Neural Operator can compute a solution in negligible time for any velocity structure or source location. We develop a scheme to train Neural Operators on an ensemble of simulations performed with random velocity models and source locations. As Neural Operators are grid-free, it is possible to evaluate solutions on higher resolution velocity models than trained on, providing additional computational efficiency. We illustrate the method with the 2D acoustic wave equation and demonstrate the method's applicability to seismic tomography, using reverse mode automatic differentiation to compute gradients of the wavefield with respect to the velocity structure. The developed procedure is nearly an order of magnitude faster than using conventional numerical methods for full waveform inversion.

References in corpus (9)

• Adam: A Method for Stochastic Optimization
• Fourier Neural Operator for Parametric Partial Differential Equations
• Finite Basis Physics-Informed Neural Networks (FBPINNs): a scalable domain decomposition approach for solving differential equations
• Multipole Graph Neural Operator for Parametric Partial Differential Equations
• Neural Operator: Graph Kernel Network for Partial Differential Equations
• Fast approximate simulation of seismic waves with deep learning
• Solving the wave equation with physics-informed deep learning
• A General Approach to Seismic Inversion with Automatic Differentiation
• Integrating Deep Neural Networks with Full-waveform Inversion: Reparametrization, Regularization, and Uncertainty Quantification