Solving the frequency-domain acoustic VTI wave equation using physics-informed neural networks. Issue 2 (11th January 2021)

• Record Type: Journal Article
• Authors: Song, Chao; Alkhalifah, Tariq; Waheed, Umair Bin
• Is Part Of: Geophysical journal international. Volume 225:Issue 2(2021)
• Page Start: 846
• Page End: 859
• Publication Date: 2021-01-11
• Subjects: Neural networks, fuzzy logic -- Numerical modelling -- Seismic anisotropy -- Wave propagation; Geophysics -- Periodicals; 550
• DOI: 10.1093/gji/ggab010
• Languages: English
• ISSNs: 0956-540X

Abstract

SUMMARY: Frequency-domain wavefield solutions corresponding to the anisotropic acoustic wave equation can be used to describe the anisotropic nature of the Earth. To solve a frequency-domain wave equation, we often need to invert the impedance matrix. This results in a dramatic increase in computational cost as the model size increases. It is even a bigger challenge for anisotropic media, where the impedance matrix is far more complex. In addition, the conventional finite-difference method produces numerical dispersion artefacts in solving acoustic wave equations for anisotropic media. To address these issues, we use the emerging paradigm of physics-informed neural networks (PINNs) to obtain wavefield solutions for an acoustic wave equation for transversely isotropic (TI) media with a vertical axis of symmetry (VTI). PINNs utilize the concept of automatic differentiation to calculate their partial derivatives, which are free of numerical dispersion artefacts. Thus, we use the wave equation as a loss function to train a neural network to provide functional solutions to the acoustic VTI form of the wave equation. Instead of predicting the pressure wavefields directly, we solve for the scattered pressure wavefields to avoid dealing with the point-source singularity. We use the spatial coordinates as input data to the network, which outputs the real and imaginary parts of the scattered wavefields and auxiliary function. After training a deep neural network, we can evaluate the wavefield at any point in space almost instantly using this trained neural network without calculating the impedance matrix inverse. We demonstrate these features on a simple 2-D anomaly model and a 2-D layered model. Additional tests on a modified 3-D Overthrust model and a 2-D model with irregular topography further validate the effectiveness of the proposed method.