Abstract:Among offshore hydrocarbon exploration technologies, the controlled source electromagnetic (CSEM) method have gained a lot of interest in both academia and industry because of its ability to detect hydrocarbon reservoirs [1]. In order to maximally extract information from the data, a full nonlinear inversion approach is employed [2]. In such an approach, the investigation domain is subdivided into pixels and by using an optimization process, the conductivity distribution in the domain can be generated and hence the location, the shape and the conductivity of the reservoir can be inferred. Most of these nonlinear methods use iterative schemes where the conductivity distribution is updated in each iteration based on a search direction computed from a gradient of the cost function. Therefore, in these derivative-based approaches the Jacobian matrix plays a key role. The elements of this matrix are the derivative of the simulated data with respect to the pixel conductivities. Its size is equal to the number of measurement data times the number of unknown pixels. In the CSEM data inversion, the size of the data set and the inversion region can be very large. Hence, the storage of the Jacobian matrix requires a huge amount of memory. This is one of the bottlenecks of using gradient-type nonlinear inversion approaches. Moreover, because the Jacobian matrix is a dense matrix, the arithmetic operation of a matrix-vector multiplication can be very expensive as the size of the Jacobian matrix increases. To reduce the size of the Jacobian matrix, we usually invert a subset of the data at the risk of missing important data points. An alternative way is to compress the Jacobian matrix based on the fact that the electromagnetic field has a limited spatial bandwidth. We can either use physics-based techniques that rely on the field kernel or use pure numerical methods that tend to be very expensive to compute. In this work, we use the adaptive cross approximation (ACA) technique introduced by Bebendorf [3] to compress the Jacobian matrix. The ACA technique converts the Jacobian matrix into two smaller rectangular matrices. This approach reduces both memory usage and CPU time of the Gauss-Newton inversion approach as well as stabilizes the inversion process. To demonstrate the Gauss-Newton inversion using the compressed Jacobian matrix, we employ inversion examples in two-and-half dimensional (2.5D) geometries.

Computing sparse derivatives and consecutive zeros problem

Computing a sparse Jacobian matrix by rows and columns

Applying Analytical Derivative and Sparse Matrix Techniques to Large-Scale Process Optimization Problems

A new approximate descent derivative-free algorithm for large-scale nonlinear symmetric equations

Computing a Sparse Approximate Inverse on Quantum Annealing Machines

Keeping it sparse: Computing Persistent Homology revisited

A structured diagonal Hessian approximation method with evaluation complexity analysis for nonlinear least squares

A Matrix-Free Newton Method

Dictionary Learning for Sparse Coding: Algorithms and Convergence Analysis

Sparse Tensors and Subdivision Methods for Finding the Zero Set of Polynomial Equations

An Alternating Manifold Proximal Gradient Method for Sparse Principal Component Analysis and Sparse Canonical Correlation Analysis

SPARSE AND LOW-RANK MATRIX DECOMPOSITION VIA ALTERNATING DIRECTION METHOD

A Proximal-Gradient Homotopy Method for the Sparse Least-Squares Problem

A scalar Hessian estimation with a sparse nonmonotone line search technique for the sparse recovery problem

An Efficient Gauss--Newton Algorithm for Symmetric Low-Rank Product Matrix Approximations

What is the gradient of a scalar function defined on a subspace of square matrices ?

Computing Symmetric Positive Definite Solutions of Three Types of Nonlinear Matrix Equations

Scalable Subspace Methods for Derivative-Free Nonlinear Least-Squares Optimization

Low-complexity subspace-descent over symmetric positive definite manifold

An Efficient Modified Residual-Based Algorithm for Large Scale Symmetric Nonlinear Equations by Approximating Successive Iterated Gradients

Adaptive cross approximation for compressing the jacobian matrix in the gauss-newton inversion