Robust data driven discovery of a seismic wave equation

TL;DR

This paper presents D-WE, a machine learning framework that automatically discovers the wave equation from observed wavefield data, even with noise and sparse measurements, by combining neural networks, genetic algorithms, and physics-informed evaluation.

Contribution

The paper introduces a novel data-driven method, D-WE, for discovering physical wave equations directly from observed data using neural networks and genetic algorithms.

Findings

Successfully discovers 2D acoustic wave equation

Demonstrates robustness to noisy data

Validates effectiveness with sparse measurements

Abstract

Despite the fact that our physical observations can often be described by derived physical laws, such as the wave equation, in many cases, we observe data that do not match the laws or have not been described physically yet. Therefore recently, a branch of machine learning has been devoted to the discovery of physical laws from data. We test such discovery algorithms, with our own flavor of implementation D-WE, in discovering the wave equation from the observed spatial-temporal wavefields. D-WE first pretrains a neural network (NN) in a supervised fashion to establish the mapping between the spatial-temporal locations (x,y,z,t) and the observation displacement wavefield function u(x,y,z,t). The trained NN serves to generate meta-data and provide the time and spatial derivatives of the wavefield (e.g., u_tt and u_xx) by automatic differentiation. Then, a preliminary library of potential terms for the wave equation is optimized from an overcomplete library by using a genetic algorithm. We, then, use a physics-informed information criterion to evaluate the precision and parsimony of potential equations in the preliminary library and determine the best structure of the wave equation. Finally, we train the "physics-informed" neural network to identify the corresponding coefficients of each functional term. Examples in discovering the 2D acoustic wave equation validate the feasibility and effectiveness of D-WE. We also verify the robustness of this method by testing it on noisy and sparsely acquired wavefield data.