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.

Robust Inversion Scheme for Logging Responses Interpretation of Micro-Cylindrically Focused Logging.

A Robust Inversion of Induction Logging Responses in Anisotropic Formation Based on Supervised Descent Method

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

Pixel-based Inversion of Induction Logging Data in Anisotropic Media with Supervised Descent Method

Geological structure-guided hybrid MCMC and Bayesian linearized inversion methodology

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

Inversion of Resistivity and Permittivity Using Phase Induction Logging Data

Modeling and Inversions of Acoustic Reflection Logging Imaging Using the Combined Monopole–dipole Measurement Mode

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

A Supervised Descent Learning Technique for Solving Directional Electromagnetic Logging-While-Drilling Inverse Problems

Theory and Application of Geostatistical Inversion: A Facies-Constrained MCMC Algorithm

The Discrete Newton Method of the Three Parameters Inversion in Electrical Logging

FMG_INV, a Fast Multi-Gaussian Inversion Method Integrating Well-Log and Seismic Data

Layered multi-exponential inversion method on surface magnetic resonance sounding dataset

Three-dimensional Regularized Inversion of Magnetotelluric Data with a Minimum Support Gradient Constraint

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

Research on Quantitative Inversion Characterization of High-Definition Electrical Imaging Logging in Oil-Based Mud Based on BPNN and MPGA-LM Algorithm

An Intelligent Inversion Method for Azimuth Electromagnetic Logging While Drilling Measurements

Near-Borehole Formation Acoustic Logging Imaging: A Full Waveform Inversion Algorithm in Cylindrical Coordinates

Joint Data and Physics Model Driven Full-Waveform Inversion Using CMP Gathers and Well-Logging Data

A Semianalytic Solution to the Response of A Triaxial Induction Logging Tool in A Layered Biaxial Anisotropic Formation