2020/07/31 by Umair bin Waheed, Ehsan Haghighat, Tariq Alkhalifah +2 · 196 citations
Earth and Planetary Sciences · Mathematics · Physics and Astronomy · #Algorithm #Applied mathematics #Artificial intelligence #Artificial neural network #Computation #Computer science #Eikonal approximation #Eikonal equation #Mathematical analysis #Mathematical optimization #Mathematics #Model Reduction and Neural Networks #Partial differential equation #Physics #Seismic Imaging and Inversion Techniques #Seismic Waves and Analysis #Solver #Weighting #math-ph #math.MP #physics.comp-ph
paper · pdf · doi:10.1016/j.cageo.2021.104833
published in Computers & Geosciences 155, 104833 (Elsevier BV)
arxiv created 2021/05/17 · openalex created_date 2021/05/24 · openalex publication_date 2021/06/04 · arxiv updated 2021/07/07 · openalex updated_date 2026/08/05
The eikonal equation is utilized across a wide spectrum of science and engineering disciplines. In seismology, it regulates seismic wave traveltimes needed for applications like source localization, imaging, and inversion. Several numerical algorithms have been developed over the years to solve the eikonal equation. However, these methods require considerable modifications to incorporate additional physics, such as anisotropy, and may even breakdown for certain complex forms of the eikonal equation, requiring approximation methods. Moreover, they suffer from computational bottleneck when repeated computations are needed for perturbations in the velocity model and/or the source location, particularly in large 3D models. Here, we propose an algorithm to solve the eikonal equation based on the emerging paradigm of physics-informed neural networks (PINNs). By minimizing a loss function formed by imposing the eikonal equation, we train a neural network to output traveltimes that are consistent with the underlying partial differential equation. We observe sufficiently high traveltime accuracy for most applications of interest. We also demonstrate how the proposed algorithm harnesses machine learning techniques like transfer learning and surrogate modeling to speed up traveltime computations for updated velocity models and source locations. Furthermore, we use a locally adaptive activation function and adaptive weighting of the terms in the loss function to improve convergence rate and solution accuracy. We also show the flexibility of the method in incorporating medium anisotropy and free-surface topography compared to conventional methods that require significant algorithmic modifications. These properties of the proposed PINN eikonal solver are highly desirable in obtaining a flexible and efficient forward modeling engine for seismological applications.