CSEM Optimization Using the Correspondence Principle

Traditionally, 3D modeling of marine controlled-source electromagnetic (CSEM) data (in the frequency domain) involves high-memory demand, requiring solving a large linear system for each frequency. To address this problem, we propose to solve Maxwell's equations in a fictitious dielectric medium with time-domain finite-difference methods, with the support of the correspondence principle. As an advantage of this approach, we highlight the possibility of its implementation for execution with GPU accelerators, in addition to multi-frequency data modeling with a single simulation. Furthermore, we explore using the correspondence principle to the inversion of CSEM data by calculating the gradient of the least-squares objective function employing the adjoint-state method to establish the relationship between adjoint fields in a conductive medium and their counterparts in the fictitious dielectric medium, similar to the approach used in forward modeling. We validate this method through 2D inversions of three synthetic CSEM datasets, computed for a simple model consisting of two resistors in a conductive medium, a model adapted from a CSEM modeling and inversion package, and the last one based on a reference model of turbidite reservoirs on the Brazilian continental margin. We also evaluate the differences between the results of inversions using the steepest descent method and our proposed momentum method, comparing them with the limited-memory BFGS (Broyden-Fletcher-Goldfarb-Shanno) algorithm (L-BFGS-B). In all experiments, we use smoothing by model reparameterization as a strategy for regularizing and stabilizing the iterations throughout the inversions. The results indicate that, although it requires more iterations, our modified momentum method produces the best models, which are consistent with results from the L-BFGS-B algorithm and require less storage per iteration.