The long-offset transient electromagnetic method (LOTEM) is widely employed in geophysical exploration, including environmental investigation, mineral exploration, and geothermal resource exploration. However, most interpretations of LOTEM data assume a flat Earth, and the commonly used one-dimensional (1D) interpretation encounters significant challenges in achieving reliable geological interpretations when topography is ignored. To address these challenges, this study presents an effective three-dimensional (3D) LOTEM inversion method. In this study, we discretize the simulation domain using unstructured tetrahedra to accurately simulate complex geological structures. The finite-element time-domain (FETD) method is utilized to calculate the LOTEM forward responses, and the limited-memory BFGS (L-BFGS) optimization method is employed for 3D LOTEM inversion. To avoid explicit calculation of sensitivity, we obtain the product of the transposed sensitivity matrix and the vector through adjoint forward modeling. Several synthetic models are used to verify the developed program, and the influence of topography on LOTEM inversion is examined. The numerical results demonstrate that topography can significantly impact the inversion result, potentially leading to incorrect geological interpretations. Finally, the developed inversion algorithm is applied to a realistic ore model from Voisey’s Bay, Labrador, Canada. The 3D inversion successfully reconstructs the spatial distribution of the ore body, further confirming the effectiveness of the developed algorithm.
When encountering sedimentary rocks with obvious laminations or fracture development zones, the conductivity of the conductive medium in different directions will change significantly, and the subsurface medium will exhibit macroscopic conductivity anisotropy. To analyze the impact of electrical anisotropy on the surface–borehole transient electromagnetic exploration method, we used the finite element method to investigate the electrical anisotropy surface–borehole transient electromagnetic three-dimensional (3D) forward algorithm, in which we used a tetrahedral mesh to spatially discretize the time–domain Maxwell equation. Then, we discretized it using the second-order backward Eulerian difference method, and we obtained the fields through the PARDISO solver. The validity and correctness of the algorithm were verified through comparison of a one-dimensional (1D) anisotropic model, a complex three-dimensional (3D) isotropic model, and a three-dimensional (3D) anisotropic half-space model. A typical anisotropic geological model was constructed to analyze the effects of anisotropic strata and anomalies in the different principal axis directions on the surface–borehole transient electromagnetic response. The results show that the response of the anisotropic medium is related to the direction of the transmitting source, and the response pattern is complex and volatile. The electrical anisotropic anomaly does affect the amplitude, which should be given special attention when performing surface–borehole transient electromagnetic inversion interpretation.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
customersupport@researchsolutions.com
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.