DIRECT IMAGING USING PHYSICS INFORMED NEURAL NETWORKS

Imaging is a crucial inversion-based task in fields ranging from medical to structure investigations, and Earth discovery. The exploding reflector assumption provides a direct imaging approach for zero-offset (coincident source-receiver) data, like ground penetrating radar (GPR) data. In imaging, however, we face aliasing problems when the data are coarsely sampled. Formulating the corresponding frequency-domain wavefield as a neural network (NN) function of the lateral and depth coordinates, as well as frequency, allows us to use the physics-informed neural network (PINN) framework to obtain images of the subsurface. In this case, we use a modified Helmholtz equation that incorporates the data on the Earth surface (hard constraint) as the loss function to optimize the NN function. This modified Helmholtz formulation allows us to avoid the inherent weaknesses that PINN has in handling boundary conditions, like the data on the surface as an additional loss term. The frequency dimension allows for image reconstruction by directly summing the wavefield over frequencies (the zero-time imaging condition). This zero-offset implementation serves as a proof of concept for later extensions to prestack data using the double square-root equation.