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.

Super Large-Scale Magnetic Data Inversion

Advanced Computational-Time Reduction Technology of Nonlinear Fourier-Based and Magnetic Circuit Hybrid Model

Improved <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M1"><mml:mrow><mml:mi>k</mml:mi></mml:mrow></mml:math>-<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M2"><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:math> PCA Algorithm Using Artificial Sparsity in Dynamic MRI

The Modified Solutoin for Large Sparse Symmetric Linear Systems in Electromagnetic Field Analysis

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

Efficient Solution Scheme for Large-Scale Anisotropic Forward Modeling of 3-D Magnetotelluric Data

Hyperspectral Unmixing Based on Incremental Kernel Nonnegative Matrix Factorization.

Three-Dimensional Magnetotelluic Parallel Inversion Algorithm Using Data Space Method

Reducing computational costs in large scale 3D EIT by using a sparse Jacobian matrix with block-wise CGLS reconstruction

Parallel Inversion of 3D Airborne Transient Electromagnetic Data Using an Approximate Jacobi Matrix

An Efficient Multigrid Solver Based on a Four-Color Cell-Block Gauss-Seidel Smoother for 3D Magnetotelluric Forward Modeling

A compressed implicit Jacobian scheme for 3D electromagnetic data inversion

An efficient three-dimensional (3D) inversion algorithm for recovering magnetic targets when remanence exists

Three-dimensional Magnetotelluric Inversion Based on Multigrid Algorithm

Robust 3D Joint Inversion of Gravity and Magnetic Data: A High-Performance Computing Approach

Acceleration of 2-D Multiplicative Regularized Contrast Source Inversion Algorithm Using Paralleled Computing Architecture

Improving the Use of the Randomized Singular Value Decomposition for the Inversion of Gravity and Magnetic Data

Parallelized 3-D Csem Inversion With Secondary Field Formulation And Hexahedral Mesh

Hierarchical Matrices Method and Its Application in Electromagnetic Integral Equations

A Stochastic Gradient Descent Approach with Partitioned-Truncated Singular Value Decomposition for Large-Scale Inverse Problems of Magnetic Modulus Data

Accelerated MRI Reconstruction With Separable and Enhanced Low-Rank Hankel Regularization