Abstract:The electromagnetic wave signal from the electromagnetic field source generates induction signals after reaching the target geological body through the underground medium. The time and spatial distribution rules of the artificial or the natural electromagnetic fields are obtained for the exploration of mineral resources of the subsurface and determining the geological structure of the subsurface to solve the geological problems. The goal of electromagnetic data processing is to suppress the noise and improve the signal-to-noise ratio and the inversion of resistivity data. Inversion has always been the focus of research in the field of electromagnetic methods. In this paper, the three-dimensional borehole-surface resistivity method is explored based on the principle of geometric sounding, and the three-dimensional inversion algorithm of the borehole-surface resistivity method in arbitrary surface topography is proposed. The forward simulation and calculation start from the partial differential equation and the boundary conditions of the total potential of the three-dimensional point current source field are satisfied. Then the unstructured tetrahedral grids are used to discretely subdivide the calculation area that can well fit the complex structure of subsurface and undulating surface topography. The accuracy of the numerical solution is low due to the rapid attenuation of the electric field at the point current source and the nearby positions and sharply varying potential gradients. Therefore, the mesh density is defined at the local area, that is, the vicinity of the source electrode and the measuring electrode. The mesh refinement can effectively reduce the influence of the source point and its vicinity and improve the accuracy of the numerical solution. The stiffness matrix is stored with Compressed Row Storage (CSR) format, and the final large linear equations are solved using the Super Symmetric Over Relaxation Preconditioned Conjugate Gradient (SSOR-PCG) method. The quasi-Newton method with limited memory (L_BFGS) is used to optimize the objective function in the inversion calculation, and a double-loop recursive method is used to solve the normal equation obtained at each iteration in order to avoid computing and storing the sensitivity matrix explicitly and reduce the amount of calculation. The comprehensive application of the above methods makes the 3D inversion algorithm efficient, accurate, and stable. The three-dimensional inversion test is performed on the synthetic data of multiple theoretical geoelectric models with topography (a single anomaly model under valley and a single anomaly model under mountain) to verify the effectiveness of the proposed algorithm.

Extraction of High Chargeability Distribution using Joint Inversion of Direct Current Resistivity and Time Domain Induced Polarization Data: A Development and Test Results on Uranium Bearing Target Rock

Reservoir evaluation method for complex resistivity using the borehole–surface electromagnetic method: A case study of an igneous reservoir in the K exploration area, China

Electrode array and data density effects in 3D induced polarization tomography and applications for mineral exploration

Inversion for 3D Conductivity and Chargeability Models Using EM Data Acquired by the New Airborne TargetEM System in Ontario, Canada

Two-dimensional electrical resistivity tomography (ERT) and time-domain-induced polarization (TDIP) study in hard rock for groundwater investigation: a case study at Choutuppal Telangana, India

2-D transdimensional joint inversion of radio magnetotelluric and electrical resistivity tomography data

Inversion Algorithms for Large-Scale Geophysical Electromagnetic Measurements

Joint Inversion of Potential Field Data with Adaptive Unstructured Tetrahedral Mesh

Improved methodology for deep aquifer characterization using hydrogeological, self-potential, and magnetotellurics data

Joint Inversion of DC Resistivity and MT Data using Multi-Objective Grey Wolf Optimization

Pseudo-3D Constrained Inversion of Transient Electromagnetic Data for a Polarizable SMS Hydrothermal System in the Deep Sea

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

The Appication Of Induced Polarization Method In The Uranium Deposit Of The Meng Gongjie

A new automatic geo-electric self-potential imaging technique for diverse sustainable development scenarios

Joint Inversion of DC Resistivity and Gravity Data with Undulating Terrain Based on Deformed Hexahedral Mesh

Stability and uncertainty assessment of geoelectrical resistivity model parameters: a new hybrid metaheuristic algorithm and posterior probability density function approach

Application of Induced Polarization Method in Mineral Resource Exploration

Anisotropic induced polarization modeling with neural networks and effective medium theory

Coupled Hydrogeophysical Inversion to Identify Non-Gaussian Hydraulic Conductivity Field by Jointly Assimilating Geochemical and Time-Lapse Geophysical Data

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

Inversion of downhole resistivity properties through infrared spectroscopy and whole‐rock geochemistry using machine‐learning