Next Article in Journal
Removal Mechanism of Mineral Impurities in Molybdenum Concentrate Treatment Process
Next Article in Special Issue
Three-Dimensional Inversion of Multi-Component Semi-Airborne Electromagnetic Data in an Undulating Terrain for Mineral Exploration
Previous Article in Journal
Organic Geochemistry and Hydrocarbon Generation Characteristics of Shale of the Fourth Member of Yingcheng Formation in the South Shuangcheng Fault Depression, Songliao Basin
Previous Article in Special Issue
Subsurface Insights of the Maricunga Gold Belt through Local Earthquake Tomography
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interpreting Surface Large-Loop Time-Domain Electromagnetic Data for Deep Mineral Exploration Using 3D Forward Modeling and Inversion

1
Department of Earth and Space Sciences, Southern University of Science and Technology, Shenzhen 518055, China
2
Guangdong Provincial Key Laboratory of Geophysical High-Resolution Imaging Technology, Southern University of Science and Technology, Shenzhen 518055, China
3
Guangdong Provincial Geophysical Prospecting Team, Guangzhou 510800, China
*
Author to whom correspondence should be addressed.
Minerals 2023, 13(1), 34; https://doi.org/10.3390/min13010034
Submission received: 9 November 2022 / Revised: 14 December 2022 / Accepted: 23 December 2022 / Published: 26 December 2022

Abstract

:
Surface transient electromagnetic (TEM) data with large transmitter loops for deep mineral exploration are often complicated by the non-trivial coupling between extended sources and arbitrarily oriented geological targets. This case study reports a TEM field data set acquired across terranes with strong lateral inhomogeneity, which is responsible for the high inconsistency in TEM data patterns along the survey line, as well as for the negative TEM transients (sign reversal) at some near-central loop stations. 3D forward modeling and inversion, as maturing tools in recent years, offer unique opportunities to extract as much geological information from such data as possible. 3D forward simulations of representative synthetic models found that the phenomenon of sign reversal at some TEM stations is associated with compact conductors enclosed by the transmitter loop and receivers that are in the loop, but off the conductor–a situation that is common in large-loop TEM and can only be explained by 3D models. However, 3D inversion of the field data with a uniform subspace as the initial and reference model fails to converge, another point of evidence that 3D inversions of large-loop TEM data are more likely to be subject to stability issues. Our solution is to warm-start the inversion with the representative model in the forward simulation experiments as the initial model, so the ill-posed 3D inversion can escape from local minima. Finally, the vertical contact structure in our 3D-inversion model is verified by a resistivity cross section of the CSAMT method. Our case study demonstrates the demand and capability of 3D electromagnetic modeling and inversion for high-resolution deep mineral exploration. It also provides an easy-to-follow template for carrying out 3D interpretation for complex geology in practice.

1. Introduction

Deep mineral exploration relies on cost-effective geophysical methods to obtain the distribution and geometry of ore bodies or ore-bearing structures at the deposit scale [1,2,3,4]. Electromagnetic (EM) methods are tools that have proven to be particularly useful for characterizing the subsurface with electrical conductivity (or its reciprocal, resistivity), as the mineralization is often associated with anomalies in the conductivity property [5,6,7,8]. Time-domain electromagnetic (or transient electromagnetic, TEM) is one of the most commonly used EM methods for minerals due to its high efficiency and high signal–noise ratio compared to other EM methods [9,10,11,12].
A typical TEM survey turns off the transmitter current in a closed wire loop or in a grounded wire line laid on the surface; the abrupt change of current in the wire induces secondary currents in the subsurface, and such secondary fields decay with time; receivers on the surface measure the EM fields as a function of time; qualitative and quantitative analysis of the decay pattern of EM fields provide information about conductivity distribution—a slower decay may indicate the presence of high-conductivity ore bodies or geological objects, and the early and late-time measurements correspond to shallow and deep structures, respectively. In mineral exploration, every conceivable configuration of transmitter and receiver has been tried for TEM to achieve the optimal coupling with the target. There are two configurations that are more frequently used in practice. The first is the central loop configuration, which utilizes a small transmitter loop and a receiver coil in the center (or close proximity) of the loop; both the transmitter and receiver move together and take measurements at a series of locations [13]. Such a configuration can illuminate a geological target from multiple locations, yielding datasets with high spatial resolution. However, a small transmitter loop may not carry a sufficient source excitation moment (the product of loop area and current) for deep penetration, and rough topography or thick vegetation can significantly restrict the deployment and movement of the system. As a result, a large and fixed loop configuration has been developed so there would be only one (or a few) stationary transmitter loop(s), yet multiple receivers could measure in or around the loop(s) [14,15].
The interpretation method of such large and fixed loop TEM data has evolved. Conventionally, data are only acquired in the central region of the loop, where the primary magnetic field at the exploration depth is supposed to be purely vertical; thus, under the assumption of 1D layered earth, the induced current is purely horizontal and the secondary magnetic field is the strongest in the vertical component on the surface. This assumption significantly simplifies the data interpretation, since the concept of apparent conductivity is well-defined and the time channels can be straightforwardly transformed to the depths [16,17]. However, significant lateral variation in geology, as one often encounters in mineral exploration, can violate the assumption and complicate the interpretation. In addition, when the receivers are outside of the source loop, TEM data may have negative responses in the decay curves because of the non-trivial geometric coupling; as a result, the decay rate of the transient signals after source turn-off is no longer indicative of the subsurface conductivity [18]. A common practice in the mining industry is to utilize plate modeling (or other simple geometric objects) to roughly reproduce the TEM responses in 3D [19]. This approach is effective in capturing the position, orientation, and geometry of compact ore bodies that are regularly shaped, but cannot properly handle complex objects or multiple objects simultaneously.
Rigorous 3D modeling and inversion that discretize the subsurface space into many small voxels have been developed to simulate arbitrarily complex geology [20,21,22,23]. The theories for quantitatively interpreting TEM data have been well established, as previous publications in the last decades have proposed a number of numerical approaches to solve both the forward and inverse problems [24,25,26,27]. Despite the advancements in theoretical studies, unfortunately, 3D modeling and inversion are still mostly absent in industrial production. One of the most obvious reasons is that the computational cost, particularly the time required by running the modeling programs, is often prohibitively high [28,29,30]. 3D inversion is also subject to severer non-uniqueness and convergence problems as a result of the drastically increased number of unknowns compared to 1D or 2D [31]. Such problems become much more challenging for the large-loop TEM, in which the subsurface is only illuminated in a fixed way [32]; the negative transient data (sign reversal) because of the out-loop receivers and increased likelihood of complex source–target–receiver coupling relations add the difficulty of fitting zeros-crossing data, which, according to our experiences, requires careful tuning of the initial model and data noise floor to ensure the convergence of inversion. The time cost, stability issue, and few successful know-how cases are all responsible for the slow adoption of 3D-voxel modeling and inversion by the mining industry.
In the present work, we aim to advocate for the application of 3D forward and inverse modeling technology to the practice of mineral exploration by documenting an inspiring case study of large-loop TEM. In particular, we overcome the difficulties in the computing time with a novel parallel solving algorithm called “survey decomposition” [18]. Then, the accelerated 3D modeling facilitates our investigation into the phenomenon of negative transient TEM data, which are deemed to be the cause of instability in the 3D-voxel inversion, and helps us find a more numerically stable initial conductivity model for the inversion. In the end, we demonstrate how a fine-tuned 3D TEM inversion reveals the mineralizing geology characterized by high lateral inhomogeneity that conventional TEM data interpretation methods cannot properly image.

2. Geological Setting

The Shizui polymetallic mining district in Guangning, Guangdong is situated at the southeastern edge of the Guangdong–Guilin Uplift, northwest of the Wuchuan–Sihui Fault Zone, and north of the Guangning Anticlinorium. The area is composed of siliceous rocks, quartz sandstones, and Cambrian black shales, while being intruded by granite at both ends. Two major fault zones have been discovered in the district: the Luoding–Guangning Fault Zone, with strong lithification and silicification, as well as gold and rare earth ore deposits; and the Wuchuan–Sihui Fault Zone, consisting of several north–east-oriented granitic mylonite (Figure 1).
The north–east-oriented Muge Fault, Gushui Fault, and north–south-oriented Chenhuai Fault manifest the tectonic history in the survey area as the primary faults, while the faults in north–east, north–northeast, northwest, and northeast–east directions are the secondary faults of higher maturity. Faults are the main ore-bearing structures in the region, among which the northeast and north–northeast (nearly north–south) oriented faults are the primary hosts of the polymetallic ores. The north–east oriented faults are composed of silicified rocks, silicified faulted rocks, phyllite, and a small number of broken quartz dikes and carbonaceous breccias. The vital mineral alteration process includes limonitization, and secondly, silicification, pyritization, and chalcopyritization, which all contribute to the formation of the gold deposits in the mining area. The northeast–east (nearly east–west) faults developed in the silicified rock belt and deposited gold and tungsten.
Extensively exposed in the area, the Cambrian Gaotan Formation is made of fine-grained quartz sandstone interbedded with siltstone and slate, and a layer of thick laminated grayish-white sandy quartz sandstone at the bottom. At the bottom lamination’s center is a ‘silica cap’ of low conductivity being developed. The Cambrian Niujiaohe Formation is mostly exposed on the northwestern side, composed of siltstone and metamorphic sandstone with a black slate of high conductivity at the bottom. Rock samples were collected from the area for the measurement of electrical properties. The Gaotan Formation generally has high resistivities, ranging from 8000 Ω·m to 20,000 Ω·m, with an average of about 13,000 Ω·m (8 pieces of Ordovician granites). The Niujiaohe Formation’s conductivity ranges from 46 Ω·m to 200 Ω·m, with an average of about 100 Ω·m (16 sandstones and charcoal-bearing slates).

3. Field Data Acquisition

The TEM survey, carried out by the Guangdong Provincial Geophysical Prospecting in 2019, consists of multiple NW–SE parallel lines across the major regional faults and the boundary between the Niujiaohe and Gaotan Formations. The detection area is located in a mountainous area with complex terrain, the available space on both sides of the measurement line is limited, and the field crew implemented a TEM configuration that combines the large fixed loop and moving central loop. The new configuration deploys a rectangular loop (750 m × 150 m) elongated and arranged along the line direction; multiple vertical coil receivers separated by 10 m are placed in the central region of the source loop to measure the time derivative of the vertical magnetic field (dBz/dt) as TEM data; once the measurements are finished at those receiver positions, the rectangular loop moves to the next section until the entire line of 1820 m runs out (Figure 2). The main reason for having so many receiving points within the transmitting coil is to increase the lateral resolution, provide spatial coverage along the survey line, and avoid errors in the acquisition of observations at some measurement points. The measurement points are at the central region of the loop, where the vertical component of the primary field is mostly constant. Thus, the loop is so long that data are not significantly affected by being close to the wire of the loop.
The instrument used in the survey was a GDP-32II system manufactured by Zonge International in Arizona. The source loop carries a 6-A current in a ramp-off waveform with a ramp length of 0.12 ms. The receiver coil has an effective area of 104 m2, and measures dBz/dt data at 28 time channels from 0.05262 to 24.291 ms. In this paper, we focus on Line 406, for which, six large rectangular loops were consecutively deployed to cover the range from Station 340 to 2160 (in meter). The green long rectangles in Figure 2 show the second and the third transmitter loops used in the survey. Some stations are repeatedly measured for different source loops to ensure continuity. The raw data are normalized by the receiver’s effective area and current for plotting and analysis.
The acquired TEM data along Line 406 are characterized by two distinct patterns of decay curve (Figure 3). The TEM data at receivers from the beginning of the line to Station 980 on the side of Niujiaohe Formation behave normally with a smooth and relatively slow decay all above the noise level, indicating a more conductive subsurface (Figure 4a). In contrast, the data from Station 990 to the end of the line on the side of the Gaotan Formation contains a large number of negative transient data and is approximately 7% below the noise level. The negative data generally appear as two types. The first type has a rapid decay to weak signals below the noise level of about 10−10 V/(Am2) at the late time channels (Figure 4c); the second type, as is the most evident near Station 1000, exhibits a transient from positive to negative with reasonably strong signal strengths (Figure 4b). The study area is in a mountainous area with minimum cultural noise and interference from infrastructures. Repeated measurements have confirmed that the negative transient data of the second type are not caused by instrument malfunctioning or environmental noise, and should have geological significance. The variation of decay rate from slow to fast, as observed in the data plot, is consistent with the change of rock type from the more conductive Niujiaohe Formation to the more resistive Gaotan Formation. Next, we quantitatively interpret the data using rigorous 3D forward modeling and inversion, with an emphasis on the transition zone near Station 1000, where the negative transients are the most prominent.

4. 3D Modeling and Inversion Algorithms

4.1. 3D Forward Modeling

The field data at Shizui were acquired along lines, but this does not invalidate the 3D modeling, since the effect of source geometry and topography can be fully accounted for in 3D. The numerical approaches for the forward problem mostly follow the computational framework for the time-domain EM in Oldenburg et al. (2013) [27], Yang et al. (2014) [33] and Yang and Oldenburg (2016) [18], with a coding implementation in C++. In the quasi-static regime neglecting the displacement current, the time-domain Faraday’s law, and Ampere’s law are written as
× E = B t
× μ 1 B = σ E + J s
where B is the magnetic flux intensity in Tesla, E is the electric field in V/m, Js is the current density in A/m2 imposed by external sources, μ is the magnetic permeability in H/m, and σ is the electric conductivity in S/m. Approximating the time derivative with finite differences using the backward Euler method, we obtain a second-order system for the magnetic flux
× σ 1 × μ 1 B i + 1 + B i + 1 δ t = × σ 1 J s i + 1 + B i δ t
where the superscript i denotes the magnetic flux at the time i.
The spatial domain of modeling is discretized by a 3D rectilinear mesh for its ease of use and numerical stability in practice. The differential operators in Equation (3) are approximated by the finite volume method (FVM) on a staggered grid, in which the electric field and the magnetic flux are designated to the edges and faces, respectively. Our algorithm discretizes arbitrary line sources specified by a series of waypoints to many small electric dipoles, then redistributes the dipole moments to the surrounding mesh edges using trilinear interpolation projection. The similar projection is also used to calculate the magnetic flux at arbitrary positions in the mesh. A natural boundary condition is applied to all six outer boundaries.
Discretizing Equation (3) yields a linear system of equations at every time step,
A i + 1 ( σ , δ t i + 1 ) B i + 1 = u ( J s i + 1 , B i , σ , δ t i + 1 )
The large computational costs of 3D modeling come from two sources. First, the coefficient matrix A in Equation (4) is large, sparse, and full of null spaces, making the solution of Equation (4) numerically cumbersome. Second, Equation (4) must be solved at every time step throughout the entire transient process. As a result, a single 3D TEM forward modeling can take tens of minutes or hours, slowing down the turnaround in data acquisition, quality control, and analysis. To promote the practical application of 3D exploration technology, we adopt the survey decomposition (SD) method, a framework to massively parallelize and speed up the 3D EM modeling using the concepts of “scale matching” and “local discretization”. Yang and Oldenburg (2016) [18] presented the basic theory of survey decomposition (SD) along with examples that validate the correctness of the framework. In time, SD abandons the sequential time stepping from early to late time channels; instead, it finds a characteristic time step size scale-adapted for each time channel, then steps from the beginning to the desired time channel using its characteristic step size (local discretization in time), so multiple time channels at the same receiver can be computed independently in parallel. In space, extensive sources like the large loop require a large volume of modeling, resulting in a large number of cells. Thus, SD treats a large loop as a linear combination of many magnetic dipoles (Figure 5a); the TEM responses at both in-loop (Figure 2) and out-of-loop (Figure 5a) receivers due to a single magnetic dipole source can be more efficiently computed using a local mesh that is oriented and optimized only for that particular source-receiver pair at that time channel (Figure 5b), also known as a S-R-T (Source–Receiver–Time) subproblem [18]. The final modeling results are trivial ensemble of many independent and easy-to-solve S-R-T subproblems (e.g., 16 subproblems in Figure 5a). With those measures in SD and modern parallel direct solvers, for example, PARDISO, as we use here, a typical 3D modeling can be finished in high precision in seconds. Our forward code, along with the SD approach, has been compared with the analytic solution of the uniform half-space and the semi-analytic solution of the layer model, and it can meet the computational accuracy requirements. We also carried out accuracy tests against some other 3D solutions, including FDTD [20], BEDS-FDTD [23], and shift-and-invert Krylov subspace (SAI) [34] on 3D models with a conductivity contrast up to 1:10,000 in Árnason (1999) [35], and our results were in high agreement with others.

4.2. 3D Inversion Algorithms

Our inversion is built on the Tikhonov regularization technique. The objective function to be minimized consists of two components: the data misfit term and the model norm term,
ϕ = 1 2 W d [ F ( m ) d ] 2 2 + β 1 2 W m ( m m 0 ) 2 2
where Wd is the data weighting matrix specifying the relative importance of each datum, Wm is the model weighting matrix penalizing the model deviation from a reference model m0 and the model variations in x-, y-, and z-direction, F(m) is the simulated response of the current model m, and β is the regularization parameter balancing the relative importance of data misfit and model norm. We note that in this study, the use of L-2 based measures in Equation (5) implies the noise in the data are in a normal distribution, and the recovered model is expected to be as smooth as possible while honoring the data.
We calculate the gradient of Equation (5) and set it to be zero,
g ( m ) = J ( m ) T W d T W d [ F ( m ) d ] + β W m T W m ( m m 0 ) = 0
where J(m) is the Jacobian or the sensitivity matrix dependent on the model m in nonlinear problems. Because of the nonlinearity, the solution of Equation (6) must be obtained iteratively. Expressing the current model as mk and using the Gauss–Newton method, we can compute the model update δmk+1 at the k + 1-th iteration by solving
[ J ( m k ) T W d T W d J ( m k ) + β W m T W m ] δ m k + 1 = g ( m k )
Equation (7) is commonly solved iteratively by conjugate gradient (CG) solvers that only require J and JT times a vector. The implicit method has been preferred because the operation of Jv and JTv is readily available after the forward solutions are obtained and one does not need to work with dense matrix J. In practice, the inversion often starts with a large β, so the model is simple and smooth and the data are poorly fitted at early iterations; at later iterations, β is reduced to allow more structure to be constructed for better data fit until the misfit is below a prescribed target level.

5. 3D Forward Modeling Analysis

Many studies have shown that central-loop TEM data over a 1D layered conductivity model must not show negative transients or sign reversals, unless a strong induced polarization (IP) effect is present [36]. In the field data presented in (Marchant et al., 2014) [37], the central-loop TEM data with the IP effect have two features. First, the data along the survey line may show a symmetrical distribution when crossing the chargeable object. Second, the TEM data above the chargeable object have higher values at early times and lower (even negative) values at late times, sometimes referred to as “shape reversal”. The field data in this paper does not look similar to any of the known IP features. While admitting that mineralization can often be associated with the IP effect, we consider the non-1D effect to be mostly responsible for the negative transients, because of the strong lateral inhomogeneity evidenced by the geological mapping, rock sample resistivity data and the abrupt change in TEM decay patterns. Fast 3D TEM modeling using survey decomposition effectively reduces the computing time to seconds or minutes, so exploration geophysicists can carry out quick modeling tests on such hypophysis, even during the survey operation if connected to the computational resources. In the following, we demonstrate how we advance our understanding about the unusual decays in the TEM data of Shizui with synthetic conceptual models and 3D forward modeling exercises, and how forward analysis can be critical to successful inverse modeling.

5.1. The Infinite Vertical Contact Model

The decay pattern of TEM data changes abruptly at about Station 1000. Thus, as our best initial guess, the first group of synthetic models consists of two quarter-spaces, one conductive and the other resistive, that together make a half-space representing the subsurface, and the geological contact of the quarter-spaces is located at X = 1000 m (Figure 6). We fix the conductivity of the resistive terrane (Gaotan Formation) to be 10−4 S/m and alter the conductive terrane (Niujiaohe Formation) to be 10−2, 10−1, 100, and 101 S/m as Models 1, 2, 3, and 4, respectively.
A rectangular source loop of 750 m × 150 m is laid across the geological contact (Figure 6). Five receivers, distributed symmetrically about the contact, are placed along the central line of the loop at X = 900, 950, 1000, 1050, and 1100 m; the receiver in the middle rides on the contact boundary at X = 1000 m. We adopt the same survey parameters in the field data for our synthetic studies. The 3D global mesh used to contain the model has a refined grid around the source and receiver locations with a smallest cell size of 50 m × 50 m × 50 m; the cell size geometrically expands from the refined grid to about 20,000 m away; a flat surface with no topography is assumed, but the air is modelled by some layers of cells that are assigned a conductivity of 10−8 S/m.
For the 3D forward modeling, we adopt the framework of survey decomposition (SD) to speed up the computing and return the results in seconds or minutes. For large-loop TEM, the decomposition takes place for both time and space. For example, at the receiver at X = 1000 m, 25 time channels are modelled independently and in parallel with the characteristic time step (Δt) 7.21 × 10−6 and 6.0728 × 10−4 s for the earliest and latest time channels, respectively. For the modeling of each time channel, the rectangular source loop is further split to smaller vertical magnetic dipoles in order to minimize the local mesh. The decomposition of the loop is also time-dependent; in this example, the earliest and latest time channels request the 750 m × 150 m loop to be considered as the superposition of 45 and 2 vertical magnetic dipoles respectively. Each magnetic dipole source is then bundled with the receiver and solved for that particular time channel on a local mesh, whose cell and domain size are also optimized for that particular scale of EM induction. Table 1 summarizes the decomposition of survey data associated with the receiver at X = 1000 m. The entire modeling results in 1755 independent subproblems, each of which requires an average computing time of about 10 s on a 40-core server node. When solved in parallel on 50 nodes of the “Taiyi” HPC platform at SUSTech, the total time is about 3.82 min, fast enough for us to analyze the model and its responses in a trial-and-error manner. Employing more than a thousand nodes, if available, can accelerate the computing of 3D models to near real-time.
Simulations of the first group of models confirm that the large lateral contrast in conductivity can lead to sign reversals in the observed TEM data when the source loop crosses the contact interface. When the conductive terrane is 10−2 S/m in Model 1, which has the smallest conductivity contrast, all five receivers observe positive decays (Figure 7a). However, if the conductivity continues to increase, for example, in Models 2, 3, and 4, we start to see the sign reversals at the receivers not directly over the conductive terrane (Figure 7b–d). The receivers above the conductive terrane always have positive decays and the decay rates reflect the change of conductivity, as expected.
Further inspection of the decay patterns for different conductivities reveals complex coupling and inductive interaction between the instrument and geological features, which give rise to diverse sign-reversed transients, including negative–positive (Model 2), positive–negative–positive (Model 3) and positive–negative (Model 4), depending on the specific geo-electric conditions (Figure 7). These numerical simulation experiments powered by fast 3D modeling profoundly imply that the data acquired at the stations off the good conductor can also bear critical information about its conductivity, sometimes not as decay rate, but in the form of sign reversal. Such sensitivity is desirable from the perspective of pursuing high detectability; data change drastically if the model changes, but it brings serious stability issues for recoverability; the inversion may be trapped in local minimum if initiated from a background model not sufficiently close to the true answer, as shown later in this paper.

5.2. The Vertical Contact Model with a Block

Resistivity measurements of the rock samples suggest the relatively conductive Niujiaohe Formation has an average conductivity of about 0.01 S/m. However, our numerical simulations in the previous subsection require the conductive quarter-space to be about 10 S/m to have a positive–negative transient. Such high conductivity is very rare in natural rocks and in large volume. Therefore, we reasonably infer that Niujiaohe Formation may contain isolated spots of high conductivity (10 S/m), while the majority of that formation is still intermediately conductive (0.01 S/m). In order to verify the plausibility of this conceptual model, the second group of synthetic models is set up by adding a very conductive block to the quarter-space of 0.01 S/m representing the conductive terrane. The block, directly touching the resistive terrane and outcropping on the surface, measures 500, 400, and 400 m in the x-, y-, and z-directions, respectively (Figure 8). The survey configurations and parameters remain the same as in the previous subsection.
To reproduce the positive–negative transient pattern observed in the field data with a reasonable conductivity model, for the second group of synthetic models, the block’s conductivity is assigned to be 10 and 100 S/m as Models 5 and 6, respectively. The numerical setup and computing times required by these models are similar to those of Models 1 to 4.
The 3D simulation results with a block added to the contact model show positive–negative transients at the off-conductor receiver locations, whereas the data directly above the conductive block have all-positive decay, consistent with the patterns in the field data (Figure 9). The timing of sign reversal is found to be related to the receiver position and block’s conductivity; the moment of positive–negative transition comes earlier if the receiver is closer to the conductive terrane or if the block is more conductive. In the field data, the sign reversal is usually before 10−3 s, so Model 6 with a block’s conductivity of 100 S/m can be considered as a more representative model of the reality (Figure 9b), although the exact boundaries of the rock units must be fine-tuned using inversions.
We acknowledge that the conductivity and the size of the block may be unrealistic. The high-conductivity block is required to qualitatively fit the field data, whereas in reality, the ground truth model can be in moderate conductivities with a complex 3D geometry. Additional simulations with buried, deeper, and smaller conductive blocks all generate positive-only transients at all of the receiver locations, which are inconsistent with the patterns in the field data. For these reasons, we consider Model 6 as the optimal candidate for the purpose of interpreting the sign reversed TEM data.

6. 3D Inversion of Field Data

The numerical simulation analyses using fast 3D forward modeling have revealed the complexity of large-loop TEM data in the presence of strong lateral variation in conductivity. Conventional interpretation techniques using apparent conductivity, decay constant or 1D layered earth inversion will certainly fail in such circumstance. Plate or block modeling can capture the first-order 3D geometric coupling between the instruments and the geological targets, sometimes correctly reproducing the sign reversal, but it lacks the flexibility of adding fine-scale and irregular structure to the model for an acceptable data fit. In this study, we demonstrate how 3D voxel inversions can play a key role in quantitatively interpreting the complex TEM field data, and how the inversion result offers critical information of geological significance.
Here, we concentrate on the data from the second and third rectangular source loops, whose data cover a range from Station 730 to 1300 at a 10-m spacing (Figure 2). The entire dataset is assumed a 5% relative error due to random noise; the absolute noise floor is 10−10 V/(Am2). The TEM data over the resistive terrane decays rapidly to below the noise floor, so some late channel data, even with an accidentally high magnitude, are deemed unreliable and excluded from the inversion. Realistic digital elevation model of the survey area is obtained to build the topography in the 3D model. In general, topographical factors have a relatively significant effect on the electric field data but not so much on the dB/dt data, as the magnetic field is determined by the electric current that is continuous across lateral conductivity interfaces. As such, the topography is approximated by staircasing mesh cells; the cells in the air are assigned a constant 10−8 S/m and fixed in the inversion.
Our 3D inversion uses a regularization of smoothness constraint based on Gauss–Newton’s method. Following the common procedure and assuming no a prior knowledge about the subsurface, we first carry out the inversion using a resistive homogeneous medium of 1.25 × 10−4 S/m as the initial and reference model. The general expectation on regularized inversions solved iteratively is that the data misfit gradually reduces towards the convergence target with an increase of model complexity (model norm), reflecting a continuous transition from "poor data fit and simple model" to "good data fit and complex model". Ideally, the desired final data misfit should be approximately unity, after normalized by the total number of data. Unfortunately, in our first inversion, the convergence of data misfit stalls after 16 iterations, as the total reduction of data misfit is only 7%, although the model norm has increased from 0 to beyond 106 (Figure 10a). The resultant conductivity model, with an unsatisfactory data fit, depicts a thick resistive formation overlain some conductive bulges at depth, a structure bearing no geological meaning and no indication of consistency with the expected contact between the Niujiaohe and Gaotan Formations (Figure 10b). At such a high data misfit, the inversion completely fails to fit the decay patterns in the TEM field data. We have also investigated different conductivities (e.g. 0.01 S/m) for the initial model and found that the inversion always tends to fall into local minima.
Initiating an inversion using a homogeneous model is a reasonable and convenient choice in practice. However, we have shown that the combination of large loop and strong lateral contrast significantly complicate the data pattern, as well as the linearity and the numerical stability of the inverse problem. For example, in our first group of synthetic modeling, the time channel of 3 × 10−3 s at Station 1100 can go from positive to negative, then to positive as the conductivity increases (Figure 7). A gradient-based inversion does not know that keeping increasing the conductivity can lead to an even better data fit outside of the local minimum, although that means a temporarily worse fit.
Our analysis and previous experiences suggest that 3D EM inversions often need to be warm-started; starting the inversion with some pre-existing structure in the background model so the model search starts somewhere in the correct convergence domain. Such a strategy has been proven necessary and effective in 3D airborne and semi-airborne TEM inversion [31,32]. Our experience also shows that the warm-start model does not need to be too sophisticated, and even “about-right” models as simple as a few blocks can work well.
For the TEM field data at Shizui, to start the 3D inversion in the correct convergence domain, Model 6 in Figure 8 is used as the initial and reference model, because it has the positive–negative transient pattern that is most consistent with that of the field data. The warm-start model consists of a resistive terrane of 10−4 S/m and a conductive terrane made of a highly conductive block of 100 S/m embedded in an intermediately conductive background of 10−2 S/m. The second inversion using the warm-start strategy converges after 18 iterations with a reduction of normalized data misfit from 78.56 to 1.06 and an increase of model norm by no more than two orders of magnitude (Figure 11a). The data misfit has achieved an acceptable level, as most of the decay patterns from both the conductive and resistive terranes have been well-fitted (Figure 12). In particular, we note that the positive–negative decay at Station 1020 is precisely reproduced with the correct time of sign reversal. The possibility of inverting the field data with both conductivity and chargeability models will be explored in future.
The recovered 3D model, as shown on the cross section below the survey line, clearly distinguishes two stratigraphic units of significantly different conductivities and their contact interfaces (Figure 11b). The first-order structure of the inversion model is consistent with the warm-starting model, but with local adjustment and smoothing required by fitting the field data and by the L2-norm in the regularization term. The final model, compared with the starting blocky model conceptually constructed, is much more natural and realistic. The geological contact characterized by two conductivities is interpreted as the boundary between the Niujiaohe and Gaotan Formations, and is believed to be associated with the mineralization processes of hydrothermal alteration. The conductivities in our model are sometimes unrealistically high, and do not match the values of rock samples. There can be many possible reasons. For example, the available data are only from a single survey line, so the recovered model can be an equivalent representation of the true model with much more complicated 3D features. Additional survey data and geological constraints can help refine the model in future.
Our delineation of the 3D structure of the contact a few hundred meters below the surface provides critical morphometric information about the alteration for understanding the mineralization and its economic potentials (Figure 13). As an independent verification, we obtain another geo-electric model reconstructed using CSAMT (controlled source audio magnetotelluric) data (Figure 14). The CSAMT cross section suggests a similar conductivity structure, since the vertical contact near Station 1000 is clearly shown from near surface to depth, and the vertical and horizontal ranges of the high conductor are close to those in the TEM model. The CSAMT model is relatively rich in small-scale variations in conductivity, partly because CSAMT data involve measuring the electric field data, which is known to be very sensitive to the local near-surface inhomogeneity. But the rectangular-loop TEM generally has higher efficiency in field work. The recovered conductivity contrast is located close to the F1 fault on the geological map (Figure 2). As far as we know, the F1 fault is mainly inferred from the surface outcropping, and there is no accurate geological information about its geometry and extension. Detailed study on the F1 fault would require additional geophysical and geological information in the future.

7. Conclusions

3D forward modeling and inversion are maturing and ready for use in the practice of mineral exploration. However, the computational costs and the ill-posedness of the inverse problem still make the rigorous 3D interpretation less accessible for practitioners. This paper contributes to the literature of 3D inversion application by presenting an inspiring case study on a large-loop TEM data from Shizui polymetallic mining area in Guangdong Province, China. The TEM data, acquired using large fixed rectangular loops and near-central-loop configuration, have complex and unusual transient decay patterns, including positive–negative sign reversals near the known boundary of two geological formations. Such a decay pattern cannot be fully explained by conventional interpretation methods, so we demonstrate how 3D TEM forward and inversion techniques can be applied here.
In particular, we develop a numerical framework, called survey decomposition, to carry out fast and parallel 3D simulations of arbitrary models and configurations in seconds or minutes. The fast 3D modeling algorithm facilitates our investigation of TEM sign reversal by quickly simulating and analyzing a number of representative synthetic models made based on prior geological and geophysical knowledges. Our 3D modeling finds that when the source loop crosses a contact of drastically different conductivities or encircles a good conductor, the receivers off the conductive terrane can observe the phenomenon of negative transients or sign reversals. In addition, the pattern of sign reversals can be very complex and nonlinear, depending mostly on how conductive the conductive terrane is. Such sensitivity of the off-target data, however, brings significant difficulties in convergence for 3D inversions, as we have shown that starting with a uniform medium subsurface can lead the inversion to a stuck in local minima and erroneous models. Because the gradient-based 3D inversion needs strong a priori information to initiate the inversion in the right convergence domain, we adopt one of the conceptual blocky models built during the forward modeling exercise as the initial and reference model to “warm-start” the inversion. The warm-start model, despite its simplicity, successfully captures the critical structure of the vertical contact between two formations and leads to an inversion with a good performance of convergence and a 3D conductivity model consistent with most known geological and geophysical information, including another conductivity model independently obtained from CSAMT data. The vertical contact structure delineated by our 3D inversion, especially at large depth, provides mineral exploration with critical information about the process and extent of hydrothermal alteration during mineralization. This application shows that 3D forward and inverse modeling, being the state-of-arts, can nowadays effectively and efficiently solve exploration problems in geologically complex areas and offer accurate quantitative evidences for the search of minerals. The workflow described in this paper can also serve as a template to further promote the application of 3D exploration techniques in the mining industry.

Author Contributions

Conceptualization, D.Y. and Q.L.; methodology, D.Y.; software, M.C.; validation, M.C.; formal analysis, M.C.; investigation, M.C.; resources, D.Y. and Q.L.; data curation, Q.L.; writing—original draft preparation, M.C.; writing—review and editing, D.Y.; visualization, M.C.; supervision, D.Y.; project administration, D.Y.; funding acquisition, D.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the National Natural Science Foundation of China (no. 42274102) and the Guangdong Provincial Key Laboratory of Geophysical High-resolution Imaging Technology (2022B1212010002).

Acknowledgments

Some high-performance computing tasks were supported by the Center for Computational Science and Engineering of Southern University of Science and Technology. Bo Ouyang helped with data archiving, plotting and preliminary 1D inversion. Xuping Feng collected the geological information.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Smith, R.J. Geophysics in Australian mineral exploration. Geophysics 1985, 50, 2637–2665. [Google Scholar] [CrossRef]
  2. Sabins, F.F. Remote sensing for mineral exploration. Ore Geol. Rev. 1999, 14, 157–183. [Google Scholar] [CrossRef]
  3. Hu, X.Y.; Peng, R.H.; Wu, G.J.; Wang, W.P.; Huo, G.P.; Han, B. Mineral exploration using CSAMT data: Application to Longmen region metallogenic belt, Guangdong Province, China. Geophysics 2013, 78, B111–B119. [Google Scholar] [CrossRef]
  4. Dentith, M.; Yuan, H.Y.; Johnson, S.; Murdie, R.; Piña-Varas, P. Application of deep-penetrating geophysical methods to mineral exploration: Examples from Western Australia. Geophysics 2018, 83, WC29–WC41. [Google Scholar] [CrossRef]
  5. Sheard, S.N.; Ritchie, T.J.; Christopherson, K.R.; Brand, E. Mining, environmental, petroleum, and engineering industry applications of electromagnetic techniques in geophysics. Surv. Geophys. 2005, 26, 653–669. [Google Scholar] [CrossRef]
  6. Vallee, M.A.; Smith, R.S.; Keating, P. Metalliferous mining geophysics—State of the art after a decade in the new millennium. Geophysics 2011, 76, W31–W50. [Google Scholar] [CrossRef] [Green Version]
  7. Smith, R. Electromagnetic induction methods in mining geophysics from 2008 to 2012. Surv. Geophys. 2014, 35, 123–156. [Google Scholar] [CrossRef]
  8. Guo, Z.W.; Xue, G.Q.; Liu, J.X.; Wu, X. Electromagnetic methods for mineral exploration in China: A review. Ore Geol. Rev. 2020, 118, 103357. [Google Scholar] [CrossRef]
  9. Rozenberg, G.; Tykajlo, R.; Pesowski, M.; Eaton, D. Combined mode transient electromagnetic survey for mineral exploration: A case history. In SEG Technical Program Expanded Abstracts 1985; Society of Exploration Geophysicists: Houston, TX, USA, 1985; pp. 257–258. [Google Scholar]
  10. Sasaki, Y.; Cho, S.J. 3D finite-difference modeling of time-domain electromagnetic data for mineral exploration. In Proceedings of the 10th SEGJ International Symposium, Kyoto, Japan, 20–22 November 2011; Society of Exploration Geophysicists of Japan: Tokyo, Japan, 2011; pp. 1–4. [Google Scholar]
  11. Asten, M.W.; Duncan, A.C. The quantitative advantages of using B-field sensors in time-domain EM measurement for mineral exploration and unexploded ordnance search. Geophysics 2012, 77, WB137–WB148. [Google Scholar] [CrossRef]
  12. Zeng, S.H.; Hu, X.Y.; Li, J.H.; Farquharson, C.G.; Wood, P.C.; Lu, X.S.; Peng, R.H. Effects of full transmitting-current waveforms on transient electromagnetics: Insights from modeling the Albany graphite deposit. Geophysics 2019, 84, E255–E268. [Google Scholar] [CrossRef]
  13. Swidinsky, A.; Hölz, S.; Jegen, M. On mapping seafloor mineral deposits with central loop transient electromagnetics marine central loop TEM. Geophysics 2012, 77, E171–E184. [Google Scholar] [CrossRef] [Green Version]
  14. Flores, C.; Peralta-Ortega, S.A. Induced polarization with in-loop transient electromagnetic soundings: A case study of mineral discrimination at El Arco porphyry copper, Mexico. J. Appl. Geophys. 2009, 68, 423–436. [Google Scholar] [CrossRef]
  15. Li, D.; Shan, X.; Weng, A. Applying three-dimensional inversion to the frequency-domain response converted from transient electromagnetic data for a rectangular fixed loop. J. Appl. Geophys. 2022, 196, 104489. [Google Scholar] [CrossRef]
  16. Fraser, D.C. Resistivity mapping with an airborne multi-coil electromagnetic system. Geophysics 1978, 43, 144–172. [Google Scholar] [CrossRef]
  17. Eaton, P.A. Application of an improved technique for interpreting transient electromagnetic data. Explor. Geophys. 1998, 29, 175–183. [Google Scholar] [CrossRef]
  18. Yang, D.; Oldenburg, D.W. Survey decomposition: A scalable framework for 3D controlled-source electromagnetic inversion. Geophysics 2016, 81, E69–E87. [Google Scholar] [CrossRef]
  19. Gallagher, P.R.; Ward, S.H.; Hohmann, G.W. A model study of a thin plate in free space for the EM37 transient electromagnetic system. Geophysics 1985, 50, 1002–1019. [Google Scholar] [CrossRef] [Green Version]
  20. Wang, T.; Hohmann, G.W. A finite-difference, time-domain solution for 3-dimensional electromagnetic modeling. Geophysics 1993, 58, 797–809. [Google Scholar] [CrossRef] [Green Version]
  21. Commer, M.; Newman, G. A parallel finite-difference approach for 3D transient electromagnetic modeling with galvanic sources. Geophysics 2004, 69, 1192–1202. [Google Scholar] [CrossRef] [Green Version]
  22. Commer, M.; Newman, G. An accelerated time domain finite difference simulation scheme for three-dimensional transient electromagnetic modeling using geometric multigrid concepts. Radio Sci. 2006, 41, 1–15. [Google Scholar] [CrossRef]
  23. Liu, S.B.; Chen, C.C.; Sun, H.F. Fast 3D transient electromagnetic forward modeling using BEDS-FDTD algorithm and GPU parallelization. Geophysics 2022, 87, E359–E375. [Google Scholar] [CrossRef]
  24. Newman, G.; Commer, M. New advances in three dimensional transient electromagnetic inversion. Geophys. J. Int. 2005, 160, 5–32. [Google Scholar] [CrossRef] [Green Version]
  25. Haber, E.; Oldenburg, D.W.; Shekhtman, R. Inversion of time domain three-dimensional electromagnetic data. Geophys. J. Int. 2007, 171, 550–564. [Google Scholar] [CrossRef] [Green Version]
  26. Cox, L.H.; Wilson, G.A.; Zhdanov, M.S. 3D inversion of airborne electromagnetic data using a moving footprint. Explor. Geophys. 2010, 41, 250–259. [Google Scholar] [CrossRef] [Green Version]
  27. Oldenburg, D.W.; Haber, E.; Shekhtman, R. Three dimensional inversion of multisource time domain electromagnetic data. Geophysics 2013, 78, E47–E57. [Google Scholar] [CrossRef] [Green Version]
  28. Wilson, G.A.; Cox, L.H.; Zhdanov, M.S. Practical 3D inversion of entire airborne electromagnetic surveys. Preview 2010, 146, 29–33. [Google Scholar]
  29. Zaslavsky, M.; Druskin, V.; Knizhnerman, L. Solution of 3D time-domain electromagnetic problems using optimal subspace projection. Geophysics 2011, 76, F339–F351. [Google Scholar] [CrossRef]
  30. Liu, Y.H.; Yin, C.C.; Qiu, C.K.; Hui, Z.J.; Zhang, B.; Ren, X.Y.; Weng, A.H. 3-D inversion of transient EM data with topography using unstructured tetrahedral grids. Geophys. J. Int. 2019, 217, 301–318. [Google Scholar] [CrossRef]
  31. Yang, D.K.; Fournier, D.; Kang, S.G.; Oldenburg, D.W. Recovery of compact conductors in 3D voxel inversion of time-domain EM data. In SEG Technical Program Expanded Abstracts; Society of Exploration Geophysicists: Tulsa, OK, USA, 2018; pp. 1873–1876. [Google Scholar]
  32. Yang, D.; Oldenburg, D.W. 3D inversion of total magnetic intensity data for time-domain EM at the Lalor massive sulphide deposit. Explor. Geophys. 2017, 48, 110–123. [Google Scholar] [CrossRef] [Green Version]
  33. Yang, D.; Oldenburg, D.W.; Haber, E. 3-D inversion of airborne electromagnetic data parallelized and accelerated by local mesh and adaptive soundings. Geophys. J. Int. 2014, 196, 1492–1507. [Google Scholar] [CrossRef] [Green Version]
  34. Zhou, J.M.; Liu, W.T.; Li, X.; Qi, Z.P. 3D transient electromagnetic modeling using a shift-and-invert Krylov subspace method. J. Geophys. Eng. 2018, 15, 1341–1349. [Google Scholar] [CrossRef]
  35. Árnason, K. Consistent Discretization of Electromagnetic Fields and Transient Modeling. Geophys. Dev. Ser. 1999, 7, 103–118. [Google Scholar]
  36. Flis, M.F.; Newman, G.A.; Hohmann, G.W. Induced-polarization effects in time-domain electromagnetic measurements. Geophysics 1989, 54, 514–523. [Google Scholar] [CrossRef] [Green Version]
  37. Marchant, D.; Haber, E.; Oldenburg, D.W. Three-dimensional modeling of IP effects in time-domain electromagnetic data. Geophysics 2014, 79, E303–E314. [Google Scholar] [CrossRef]
Figure 1. Regional geological map. Legend: 1. Cambrian system; 2. Sinian system; 3. Magmatite; 4. Luoding Formation; 5. Monzonite granite; 6. Binhuai Fault; 7. Muge Fault; 8. Gushui Fault; 9. Bilong Fault; 10. Wuchuan-Sihui Fault; 11. our study area.
Figure 1. Regional geological map. Legend: 1. Cambrian system; 2. Sinian system; 3. Magmatite; 4. Luoding Formation; 5. Monzonite granite; 6. Binhuai Fault; 7. Muge Fault; 8. Gushui Fault; 9. Bilong Fault; 10. Wuchuan-Sihui Fault; 11. our study area.
Minerals 13 00034 g001
Figure 2. Rectangular-loop TEM survey layout. The green lines indicate the two of the six source loops used in our inversion. є1n: Cambrian Niujiaohe Formation; єg: Cambrian Gaotan Formation; O 3 1 γ: magmatite; F1: Binhuai Fault; F2: Muge Fault.
Figure 2. Rectangular-loop TEM survey layout. The green lines indicate the two of the six source loops used in our inversion. є1n: Cambrian Niujiaohe Formation; єg: Cambrian Gaotan Formation; O 3 1 γ: magmatite; F1: Binhuai Fault; F2: Muge Fault.
Minerals 13 00034 g002
Figure 3. Measured TEM data along Line 406. The red and blue dots represent the positive and negative data, respectively.
Figure 3. Measured TEM data along Line 406. The red and blue dots represent the positive and negative data, respectively.
Minerals 13 00034 g003
Figure 4. Decay of TEM data at select receiver stations. (a) Station 960; (b) Station1040; (c) Station 1300.
Figure 4. Decay of TEM data at select receiver stations. (a) Station 960; (b) Station1040; (c) Station 1300.
Minerals 13 00034 g004
Figure 5. Diagram of survey decomposition. (a) Decomposition of a large square loop and pairing of magnetic dipole sources (Tx) and receiver (Rx); (b) diagram of an efficient local mesh customized for one of the source–receiver pair. The grids in gray and green are the global and local mesh, respectively.
Figure 5. Diagram of survey decomposition. (a) Decomposition of a large square loop and pairing of magnetic dipole sources (Tx) and receiver (Rx); (b) diagram of an efficient local mesh customized for one of the source–receiver pair. The grids in gray and green are the global and local mesh, respectively.
Minerals 13 00034 g005
Figure 6. The synthetic model of vertical contact. The quarter-space in dark blue is the resistive terrane of 10−4 S/m, and that in light blue is the conductive terrane, whose conductivity is 0.01, 0.1, 1, and 10 S/m for Model 1, 2, 3, and 4, respectively. The red line outlines the source loop; the red points mark the receiver locations.
Figure 6. The synthetic model of vertical contact. The quarter-space in dark blue is the resistive terrane of 10−4 S/m, and that in light blue is the conductive terrane, whose conductivity is 0.01, 0.1, 1, and 10 S/m for Model 1, 2, 3, and 4, respectively. The red line outlines the source loop; the red points mark the receiver locations.
Minerals 13 00034 g006
Figure 7. 3D simulated data for the first group of synthetic models. (a) Model 1; (b) Model 2; (c) Model 3; (d) Model 4. The solid and dashed lines indicate positive and negative data, respectively.
Figure 7. 3D simulated data for the first group of synthetic models. (a) Model 1; (b) Model 2; (c) Model 3; (d) Model 4. The solid and dashed lines indicate positive and negative data, respectively.
Minerals 13 00034 g007aMinerals 13 00034 g007b
Figure 8. The synthetic model of vertical contact with a block. The quarter-space in dark blue is the resistive terrane of 10−4 S/m, the light blue region is the conductive terrane of 10−2 S/m, and the yellow represents the high-conductivity block, whose conductivity is 10 and 100 S/m for Model 5 and 6, respectively. The red line outlines the source loop; the red points mark the receiver locations.
Figure 8. The synthetic model of vertical contact with a block. The quarter-space in dark blue is the resistive terrane of 10−4 S/m, the light blue region is the conductive terrane of 10−2 S/m, and the yellow represents the high-conductivity block, whose conductivity is 10 and 100 S/m for Model 5 and 6, respectively. The red line outlines the source loop; the red points mark the receiver locations.
Minerals 13 00034 g008
Figure 9. 3D simulated data for the second group of synthetic models. (a) Model 5; (b) Model 6. The solid and dashed lines indicate positive and negative data, respectively.
Figure 9. 3D simulated data for the second group of synthetic models. (a) Model 5; (b) Model 6. The solid and dashed lines indicate positive and negative data, respectively.
Minerals 13 00034 g009
Figure 10. The first inversion using a uniform subsurface as the initial and reference models. (a) Tikhonov curve: model norm versus normalized data misfit using dots to indicate the iterations; (b) recovered conductivity model sliced on the cross section Y = 0 m. The receiver locations are marked in red.
Figure 10. The first inversion using a uniform subsurface as the initial and reference models. (a) Tikhonov curve: model norm versus normalized data misfit using dots to indicate the iterations; (b) recovered conductivity model sliced on the cross section Y = 0 m. The receiver locations are marked in red.
Minerals 13 00034 g010
Figure 11. The second inversion using the warm-start strategy. (a) Tikhonov curve: model norm versus normalized data misfit using dots to indicate the iterations; (b) recovered conductivity model sliced on the cross section at Y = 0 m.
Figure 11. The second inversion using the warm-start strategy. (a) Tikhonov curve: model norm versus normalized data misfit using dots to indicate the iterations; (b) recovered conductivity model sliced on the cross section at Y = 0 m.
Minerals 13 00034 g011
Figure 12. TEM field data used in 3D inversion and the data fitting after inversion. Green lines indicate the estimated noise level.
Figure 12. TEM field data used in 3D inversion and the data fitting after inversion. Green lines indicate the estimated noise level.
Minerals 13 00034 g012
Figure 13. The final 3D conductivity model with realistic topography viewed in 3D. The red dots indicate the receiver locations involved in the inversion.
Figure 13. The final 3D conductivity model with realistic topography viewed in 3D. The red dots indicate the receiver locations involved in the inversion.
Minerals 13 00034 g013
Figure 14. Cross section of CSAMT conductivity model.
Figure 14. Cross section of CSAMT conductivity model.
Minerals 13 00034 g014
Table 1. Survey decomposition of the synthetic forward modeling at one receiver location.
Table 1. Survey decomposition of the synthetic forward modeling at one receiver location.
Time Channel (s)Δt (s)Number of Magnetic Dipole SourcesCells in Local MeshCPU Time of One Subproblem (s)
1.442 × 10−47.21 × 10−64528,800–29,7607
1.747 × 10−48.735 × 10−64228,800–30,7207
2.052 × 10−41.026 × 10−53929,760–31,6208
2.504 × 10−41.252 × 10−53629,760–31,6208
3.116 × 10−41.558 × 10−53330,720–32,6408
3.727 × 10−41.8635 × 10−53030,720–32,6748
4.482 × 10−42.241 × 10−51831,680–33,7289
5.399 × 10−42.6995 × 10−51632,640–34,7829
6.599 × 10−43.2995 × 10−51433,660–35,8369
8.271 × 10−44.1355 × 10−51434,680–35,8369
1.0249 × 10−33.416 × 10−51235,700–36,8909
1.2672 × 10−34.224 × 10−51235,700–36,8909
1.5709 × 10−35.236 × 10−5536,720–37,80010
1.9633 × 10−36.544 × 10−5536,750–37,80010
2.4785 × 10−38.262 × 10−5437,80010
3.1166 × 10−31.039 × 10−4437,80010
3.9047 × 10−31.302 × 10−4337,94410
4.8908 × 10−31.6303 × 10−4339,06011
6.1464 × 10−32.0488 × 10−4340,32011
7.7398 × 10−32.5799 × 10−4341,61611
9.7119 × 10−33.2373 × 10−4242,76812
1.2196 × 10−23.049 × 10−4244,06412
1.5348 × 10−23.8370 × 10−4245,92812
1.9306 × 10−24.8265 × 10−4246,65613
2.4291 × 10−26.0728 × 10−4249,24814
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cheng, M.; Yang, D.; Luo, Q. Interpreting Surface Large-Loop Time-Domain Electromagnetic Data for Deep Mineral Exploration Using 3D Forward Modeling and Inversion. Minerals 2023, 13, 34. https://doi.org/10.3390/min13010034

AMA Style

Cheng M, Yang D, Luo Q. Interpreting Surface Large-Loop Time-Domain Electromagnetic Data for Deep Mineral Exploration Using 3D Forward Modeling and Inversion. Minerals. 2023; 13(1):34. https://doi.org/10.3390/min13010034

Chicago/Turabian Style

Cheng, Ming, Dikun Yang, and Qiang Luo. 2023. "Interpreting Surface Large-Loop Time-Domain Electromagnetic Data for Deep Mineral Exploration Using 3D Forward Modeling and Inversion" Minerals 13, no. 1: 34. https://doi.org/10.3390/min13010034

APA Style

Cheng, M., Yang, D., & Luo, Q. (2023). Interpreting Surface Large-Loop Time-Domain Electromagnetic Data for Deep Mineral Exploration Using 3D Forward Modeling and Inversion. Minerals, 13(1), 34. https://doi.org/10.3390/min13010034

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop