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.

Hybrid matrix compression for high-frequency problems

Memory-efficient compression of $\mathcal{DH}^2$-matrices for high-frequency problems

A Pre-splitting Green's Function Based Hybrid Fast Algorithm for Multiscale Problems

Uniform $\mathcal{H}$-matrix Compression with Applications to Boundary Integral Equations

A hybrid stochastic interpolation and compression method for kernel matrices

Code Generation and Performance Engineering for Matrix-Free Finite Element Methods on Hybrid Tetrahedral Grids

A matrix-free parallel solution method for the three-dimensional heterogeneous Helmholtz equation

Performance of H-Matrix-Vector Multiplication with Floating Point Compression

Efficient Solution to the Large Scattering and Radiation Problem Using the Improved Finite-Element Fast Multipole Method

A hybrid interpolation ACA accelerated method for parabolic boundary integral operators

Towards High-Performance Linear Potential Flow BEM Solver with Low-Rank Compressions

Matrix-Free Higher-Order Finite Element Methods for Hyperelasticity

Compression of boundary integral operators discretized by anisotropic wavelet bases

Matrix-free GPU-accelerated saddle-point solvers for high-order problems in $H(\mathrm{div})$

Parallel Skeletonization for Integral Equations in Evolving Multiply-Connected Domains

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

Matrix-free multigrid block-preconditioners for higher order Discontinuous Galerkin discretisations

A compression scheme for domain decomposition method in solving electromagnetic problems

Stand-alone Multigrid for Helmholtz Revisited: Towards Convergence Using Standard Components

Exponentially convergent multiscale methods for high frequency heterogeneous Helmholtz equations