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.

Detection of Boulder Using GEM-2 EM System Based on Three-Dimensional Inversion Algorithm

Detection method for boulders in subway shield zones based on data fusion multi-resistivity three-dimensional tomography

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

Two-dimensional Inversion of CSAMT Data and Three-Dimensional Geological Mapping for Groundwater Exploration in Tongkeng Area, Hunan Province, China

3D Inversion of Multifrequency Controlled-Source Electromagnetic Data and Its Application to Geothermal Exploration in the Tianzhen Region of the Northern Datong Basin, China

Three-Dimensional ERT Advanced Detection Method with Source-Position Electrode Excitation for Tunnel-Boring Machines

Multifrequency 3-D Inversion of GREATEM Data by BCGS-FFT-BIM

2.5D Inversion of Borehole Transient Electromagnetic Method with Scanning Detection Based on RCQPSO-LMO Combined Algorithm

Boulder Detection in Metro Shield Zones Based on 3D Seismic Method, and Its Application

COMPLEX GEOLOGICAL STRUCTURE INDENTIFICATION AND PARA-METER ESTIMATION METHOD BASED ON 3D RESISTIVITY MEASURE-MENT

Borehole Transient Electromagnetic Response Calculation and Experimental Study in Coal Mine Tunnels

3-D imaging of whole-space environments with electromagnetic induction sensors based on linear approximation

A New Method for Abnormal Underground Rocks Identification Using Ground Penetrating Radar

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

3-D Joint Inversion of MT and CSEM Data for Imaging a High-Temperature Geothermal System in Yanggao Region, Shanxi Province, China

Three-dimensional inversion analysis of transient electromagnetic response signals of water-bearing abnormal bodies in tunnels based on numerical characteristic parameters

Identification of Unstable Subsurface Rock Structure Using Ground Penetrating Radar: an EEMD-Based Processing Method

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

3D inversion of audio-magnetotelluric data for mineral exploration: A case study of Layikeleke buried porphyry copper deposit, Xinjiang, China

Geological detection of hard rocks by GPR and signal time-frequency characteristics analysis in urban underground trenchless construction

Application of a New Geophone and Geometry in Tunnel Seismic Detection