Abstract:The high-resolution laterolog array (HRLA) tool can simultaneously provide six curves of logging responses with different depths of investigation, which give abundant information for the inversion of geological formation parameters. This paper proposes a trust region algorithm to recover several parameters in the axisymmetric formation from HRLA data, including the mud-invasion radius, mud resistivity and original formation resistivity.Firstly, we propose an accurate and efficient method to solve the forward problem in order to obtain the logging responses of HRLA tool in the axisymmetric formation, where the superposition principle is used to simplify the calculation of resistivity logging responses, and the parallel direct solver PARDISO is employed to solve the linear system with multiple right-hand sides resulting from finite element approximation. The decomposition of the stiffness matrix is needed only once when using PARDISO to solve these linear systems, which is cheap in terms of work and fast in time. Secondly, the problem of identifying formation parameters is transformed into a nonlinear least squares problem that can be solved by the trust region algorithm. In order to speed up the optimization process, we derive an explicit formula for the partial derivative of the objective function with respect to each resistivity parameter. Finally, for a typical formation model we present resistivity four-parameter and traditional three-parameter inversion results by using different initial guesses.In order to test the validity and effectiveness of the proposed method, we design a typical formation model with borehole radius r(h) = 0.1016 m, invasion radius r(i) = 0.4016 m, permeable formation thickness H = 2 m, borehole mud resistivity R-m = 1 Omega m, surrounding rock resistivity = 10 Omega m and original formation resistivity R-t = 500 Omega m. Comparing the forward modeling result with the result in the existing literature, we find that the maximum relative error between the two results is only 0.01%, which verifies the correctness of the forward algorithm program. It takes less than half a second to perform forward modeling with 24881 grid points using PARDISO on a personal laptop. The computational time of the Jacobi conjugate gradient (JCG) method is more than ten times of that of PARDISO, and PARDISO provides a more accurate solution than the JCG. We use the trust region algorithm to invert formation parameters with four different initial values, one of which is far away from the real model. The inversion results show that the method can converge to the true model quickly even with very poor initial guesses. When using noisy data mixed with Gaussian white noise with signal-to-noise ratio of 10dB or even 5dB, the trust inversion algorithm still gives satisfactory results.Trust region algorithm combined with PARDISO is successfully used to multi-parameter inversion of HRLA data. Numerical results show that the trust region inversion algorithm converges fast, and has a global convergence property. Another advantage of the proposed method is the robustness to noise.

Modified Levenberg–Marquardt Inversion for High-Resolution Resistivity Distribution Reconstruction of Multilayered Formation

Parameter Inversion Based on Levenberg-Marquardt Algorithm for Layered Formation Using Electromagnetic Wave Resistivity Tool

Trust Region Inversion Algorithm of High-Resolution Array Lateral Logging in Axisymmetric Formation

Efficient Inversion of 3D IP Data for Continuous Model with Complex Geometry

Real-time forward modeling and inversion of logging-while-drilling electromagnetic measurements in horizontal wells

Inversion of Resistivity and Permittivity Using Phase Induction Logging Data

Joint and laterally constrained inversion of surface MRS and TEM data

Fast hierarchical inversion for borehole resistivity measurements in high-angle and horizontal wells using ADNN-AMLM

Three-Dimensional Inversion of Borehole-Surface Resistivity Method Based on the Unstructured Finite Element

Inversion of Anisotropic Resistivity in Three-dimensional Asymmetric Media and Its Application

Three-Dimensional Anisotropic Inversion for Direct Current Resistivity Method Using Unstructured Tetrahedral Discretization

2.5D Direct‐current Resistivity Forward Modelling and Inversion by Finite‐element–infinite‐element Coupled Method

Self-adaptive inversion method of electromagnetic-propagation resistivity logging while drilling data

An Intelligent Inversion Method for Azimuth Electromagnetic Logging While Drilling Measurements

Numerical modeling of 3D DC resistivity method in the mixed space-wavenumber domain

Multiscale Resistivity Inversion Based on Convolutional Wavelet Transform

3D modeling and inversion of the electrical resistivity tomography using steel cased boreholes as long electrodes

An Inversion Method for Calculating Formation and Borehole Parameters in X-ray Lithology Density Logging

Multiscale Constrained Inversion Method For Direct Current Resistivity Tomography Based On Haar Wavelet Transform

3-D Inversion of High Density Resistivity Method Based on 2-D Exploration Data and Its Application

Inversion of the Wideband and High Frequency Electromagnetic Data