Transient Electromagnetic Research Articles

New fast imaging techniques for electrical source transient electromagnetic data: Approaches and application

The rapid imaging of electrical source transient electromagnetic (TEM) data involves two essential processes: the calculation of apparent resistivity and the conversion of time to depth. Traditionally, the definition of full-time apparent resistivity is defined by considering solely the vertical magnetic field, which is predicated on the monotonic relationship between the resistivity and the electromagnetic field response. Based on the concept of peak time, we have developed distinct methodologies for calculating the apparent resistivity for both the horizontal electric field (ex) and the vertical induced voltage (vz), which demonstrated accuracy across the entire time range examined. We also introduced a formula to address discrepancies in apparent resistivity arising from the non-dipole size effect of the source, thereby ensuring that the algorithm can adapt to any transmitting and receiving configuration. Furthermore, we provided straightforward and precise time-depth conversion equations applicable to both ex and vz, which facilitate the rapid imaging of observational data. Multiple numerical examples were employed to illustrate the effectiveness and robustness of this approach. Finally, we applied this imaging technique to the data processing of actual measured data from a survey area conducted in Ningxia Province, and the imaging results accurately reflected the distribution of the electrical structure of the subsurface strata. The innovative imaging technique presented in this study holds considerable potential for the expedited processing and analysis of ground-based and semi-aerial electrical source transient electromagnetic survey data, which are widely employed in contemporary applications.

Quantitative suppression of unsteady powerline noise in transient electromagnetic surveys: Adjustment of base-frequency and optimal choice of stacking-times.

Powerline noise is a severe interference source in urban or mine transient electromagnetic (TEM) surveys. TEM systems generally adopt the synchronous detection scheme to suppress the powerline noise. However, when the powerline frequency fluctuates differing from its nominal value (50/60Hz), the considerable powerline noise residue still remains even after synchronous detection. To overcome this problem, this paper proposes the quantitative suppression method for powerline noise taking into account the instability of powerline frequency. The method is based on the adjustment of base-frequency and optimal choice of stacking-times. We represent mathematically the sufficient condition for powerline noise suppression by synchronous detection. It consists of two equations, one with respect to base-frequency and other with respect to stacking-times. The base-frequency is adjusted according to powerline frequency estimate. We first derive the mathematical relationship between frequency estimation accuracy, residual noise amplitude, and stacking-times. Based on it, we develop an efficient algorithm to determine the optimal stacking-times. The algorithm takes as an input the powerline noise estimates and the noise tolerance limit. The tolerance limit is set to a certain value desired by the user, such as 10, 4, and 1 µV. By adjusting base-frequency and determining optimal stacking-times, the powerline noise residue after synchronous detection is reduced to below the desired tolerance limit. We verify the effectiveness of the proposed method for both simulated and actual noise. Experimental results show that the method achieves quantitative suppression of unsteady powerline noise without any damage of effective signal and prevents the measurement time loss due to excessive stacking.

Transient Electromagnetic Field Response Analysis Considering Dam Boundary Conditions

Abstract The transient electromagnetic (TEM) method has outstanding advantages in detecting depth and response information and can be used to detect anomalies hidden deep inside dams. However, the current data analysis method based on geological exploration is inapplicable to the interpretation of dam detection data, the main reason being the lack of consideration of the influence of dam boundary conditions on the diffusion of the TEM field. In this work, by establishing a numerical model of a dam for TEM detection, the effects of the dam slope ratio, upstream water body, overburden, and other dam boundary conditions on the TEM field were studied. The results showed that the slope ratio made the propagating TEM field to exhibit a constraining behavior, and a negative deviation could be seen in the attenuation curve, indicating a false high-resistance anomaly. The upstream water body distorted the diffusion of the electromagnetic field, and the low-resistance covering layer caused a low-retardation retention phenomenon in the magnetic field, leading to a positive deviation in the attenuation curve, indicating a false low-resistance anomaly. The influence of the boundary conditions was present throughout the entire observation window in which the detection was performed, and the sensitivity of the magnetic field to these boundary conditions was significantly higher than that of the change rate of the magnetic field. With the increase in the observation time, the main influencing factors were the slope ratio, overburden layer, and upstream water body.

Misidentified subglacial lake beneath the Devon Ice Cap, Canadian Arctic: a new interpretation from seismic and electromagnetic data

Abstract. In 2018 the first subglacial lake in the Canadian Arctic was proposed to exist beneath the Devon Ice Cap, based on the analysis of airborne radar data. Here, we report a new interpretation of the subglacial material beneath the Devon Ice Cap, supported by data acquired from multiple surface-based geophysical methods in 2022. The geophysical data recorded included 9 km of active-source seismic-reflection profiles, seven transient electromagnetic (TEM) soundings, and 17 magnetotellurics (MT) stations. These surface-based geophysical datasets were collected above the inferred locations of the subglacial lakes and show no evidence for the presence of subglacial water. The acoustic impedance of the subglacial material, estimated from the seismic data, is 9.49 ± 1.92 × 106 kg m−2 s−1, comparable to consolidated or frozen sediment. The resistivity models obtained by inversion of both the TEM and MT measurements show the presence of highly resistive rock layers (1000–100 000 Ω m) directly beneath the ice. Re-evaluation of the airborne reflectivity data shows that the radar attenuation rates were likely overestimated, leading to an overestimation of the basal reflectivity in the original radar studies. Here, we derive new radar attenuation rates using the temperature- and chemistry-dependent Arrhenius equation, and when applied to correct the returned bed power, the bed power does not meet the basal reflectivity threshold expected over subglacial water. Thus, the radar interpretation is now consistent with the seismic and electromagnetic observations of dry or frozen, non-conductive basal material.

A rapid imaging approach for transient electromagnetic data based on gradient features

ABSTRACT Data processing techniques in transient electromagnetic (TEM) methods face numerous challenges, such as the problem of 3-D imaging and the demand for new data processing techniques arising from the development of small-loop TEM devices. We have developed an imaging algorithm based on gradient features derived from TEM data to address these challenges. We found that normalising the gradient field of TEM responses under appropriate background conditions enables imaging of the target. Numerical experiments conducted under theoretical models have validated the feasibility of this approach and highlighted the significance of background fields. We further designed more complex models and achieved 3-D imaging of synthetic data within seconds, demonstrating the method's high efficiency and applicability to high-dimensional data. Subsequently, we applied the algorithm to challenging small-loop TEM data and demonstrated its superior performance compared to the inversion method, which presents a novel approach for processing data from small-loop TEM devices. This study demonstrates that the algorithm possesses characteristics of high efficiency and wide applicability, suitable for TEM data from various dimensions and types of devices. Furthermore, it holds promising potential in assisting the enhancement of the efficiency of the three-dimensional inversion of TEM data.

Inversion of TEM Responses in Tunnel with Steel Infrastructure and Its Application

Steel infrastructure in tunnels significantly affects the accuracy of transient electromagnetic (TEM) data inversion. This study proposes and compares two methods to effectively mitigate this interference: the correction factor method and the prior information constraint method. The correction factor method constructs correction factors using interference data obtained in two ways. The first method, numerical simulation, is theoretically precise but relies heavily on accurate model parameters and requires significant computational resources. The second is conducted by measuring the interference data in the working space. While this approach is simple and operationally convenient, its effectiveness in correction is limited by the measured data quality. The a priori information constraint method improves the inversion by introducing a more accurate initial model. The inversion results of synthetic data indicate that the inversion based on the correction factor method and the prior information constraint method can effectively invert the anomaly. However, the delay effect remains a challenge for the correction factor method. Finally, we invert the field data measured in a mine tunnel by combining the interference processing methods.

A small-loop double-coil decoupling TEM device for high sensitivity near surface detection

Small loop transient electromagnetic (TEM) devices, several meters in size, are increasingly applied to near-surface geophysical exploration in space-limited environment. But the transmitter (TX) coil significantly affects the early time data acquisition of the receiver (RX) coil during current turn-off, resulting in the loss of shallow signals and creating a significant blind zone problem. We report a double-coil decoupling TEM (DCDTEM) device to solve this problem. The novel DCDTEM device ensures a weak primary magnetic flux region between the two coplanar transmitter (TX) coils with the same current direction and different radii. Simulation results indicate that the transmitting field of the DCDTEM device is more focused compared to the traditional central loop device, ensuring a stronger coupling with the target. 3-D TEM forward modeling results show that the DCDTEM device has the highest detection sensitivity, together with a relatively high lateral resolution compared to other small-loop TEM devices. We further carry out several physical experiments to test the response characteristics of the DCDTEM device with reference to both the offset and the ratio of turns between the TX and RX coils. The experimental results ensure that the position of the zero magnetic flux point can be accurately calibrated and mitigate interference from TX coils. The proposed DCDTEM device is also applied in a field test on water pipe detection, demonstrating its practicality in near-surface detection.

Recovering induced polarization effects from 1-D coupled inversion of transient electromagnetic data

SUMMARY Induced polarization (IP) effects can significantly affect and superimpose the inductive earth response, leading to heavily distorted data and, if overlooked, false geological interpretation. In this paper, we implemented the Levenberg–Marquardt (LM) and very fast simulated annealing (VFSA) algorithms to recover IP effects from central loop transient electromagnetic (TEM) data. To incorporate the IP effect in the TEM response, we used the Cole–Cole parametrization, maximum phase angle (MPA), maximum imaginary conductivity (MIC) and Jeffrey transform of Cole–Cole parameters. The result of 1-D forward calculation and inversion of synthetic TEM data revealed that the Cole–Cole parametrization is more robust and reliable than MPA, MIC and Jeffrey transform, and that the synthetic data were well fitted and IP parameters well recovered using this model. However, the incorporation of the IP effect leads to a highly nonlinear and non-unique inverse problem which requires an accurate starting model, especially for LM inversion. To evaluate the performance of our algorithm using field data, we carried out a 1-D inversion of TEM data acquired along a profile that traverses a waste site located near Cologne, Germany. Furthermore, to obtain a priori information and validate the result of TEM data modelling, we conducted an electrical resistivity tomography (ERT) and time-domain IP (TDIP) survey along the TEM profile. A 2-D inversion was used to retrieve the Cole–Cole parameters as input for TEM interpretation. By including the IP information, the TEM field data can be explained quantitively, and a consistent and improved interpretation of the waste body is achieved.

Surface geometry inversion of transient electromagnetic data

We investigate an emerging method called surface geometry inversion (SGI) for the inversion of transient electromagnetic (TEM) data. Conventional minimum-structure inversion methods parameterize the earth model with many mesh cells within which the physical properties are constant, and construct a physical property model that is usually smoothly varying as well as fits the observations. With these smooth models, it is difficult to extract the interface between different geologic units, and it can be especially difficult to target drillholes for thin, plate-like targets, which are frequently encountered in mineral exploration. Our SGI parameterizes the model in terms of the coordinates of the nodes (vertices) used to connect together the surfaces that define the geologic interfaces. The algorithm then inverts for the locations of these nodes, which directly provides geometric information about the target. This can be more useful than a fuzzy image of conductivity, especially for an exploration project. A genetic algorithm (GA) is used to solve the nonlinear, overdetermined optimization problem. We use a finite-element solver to solve the TEM forward-modeling problem of each candidate model in the GA population. Because forward modeling is independent for each model, we implement a hybrid message passing interface (MPI) + OpenMP parallel method to improve computational efficiency. We investigate a new parameterization method specifically designed for thin, plate-like structures, which is more efficient and can effectively avoid self-intersection. We first illustrate the effectiveness of our SGI algorithm on a synthetic block model before testing the new parameterization method on a synthetic thin plate model. Finally, we apply our SGI to a real data set collected for the exploration of thin graphitic faults in a uranium exploration project in Canada. The constructed model from our SGI corresponds well with the drilling data.

Use of remote sensing, spatial and geophysical modeling, and real recharging capabilities to identify suitable areas for groundwater exploitation in dry coastal areas

Accurate identification of groundwater potential areas in arid regions is an important task for groundwater management and sustainability. As a result, this study used the innovative integration of remote sensing (RS), geographic information system (GIS), watershed modeling system (WMS), geophysical survey, and water mass balance equation to identify potential groundwater areas in the W. Dara, Eastern Desert, Egypt. A weighted spatial probability model (WSPM) of groundwater potential based on eight regulatory factors was implemented within ArcGIS software. Drainage density (DD), precipitation (P), net groundwater recharge (NGR), terrain slope (TS), lineament density (LD), lithologic group (LG), water quality (TDS), and depth to groundwater level (DGW) are the aspects considered. The Analytical hierarchy process (AHP) method was used to assign weights to these parameters, and their accuracy was estimated using the consistency ratio (CR). The resulting groundwater potential map classified W. Dara study area into five categories, ranging from very low to very high potential. A geophysical survey, in the form of Vertical Electrical Sounding (VES) and Transient Electromagnetic (TEM), was conducted along W. Dara to validate the results of the WSPM, which identified areas of high groundwater potential. The 1D inversion of VES/TEM shows that the central and western parts of W. Dara are considered the most promising areas for groundwater occurrence, and are located in areas of high and very high potential classes derived from WSPM. Moreover, the results of VES and TEM surveys showed that the proposed aquifers (Nubian Sandstone, Miocene, and Quaternary) in the study area are horizontally and vertically connected through a set of normal faults traversing NW-SE. Ten sites have been proposed for drilling additional exploitative wells in W. Dara area based on the WSPM and geophysical survey with the aim of sustainable development. Thus, the integrated techniques applied in this study proved effective in accurately determining the development strategy for arid and semi-arid coastal areas, especially those that suffer from scarcity of rainfall and increased agricultural reclamation requirements in remote areas.

Evaluating the applicability of magnetic resonance sounding for groundwater exploration at selected locations in Aswan, Egypt

ABSTRACT Magnetic resonance sounding (MRS) was performed to evaluate its applicability for characterising groundwater at selected locations in Aswan Governorate, Egypt. The MRS mainly detects a weak but noticeable alternating magnetic field as MRS signals when the proton spin in groundwater produces non-magnetic moments and rotates around the geomagnetic field, indicating the detection of groundwater. The MRS system was first correlated at a piezometer with a known water table depth (~80 m). Subsequently, three MRS sounding stations with unknown water tables were tested at Qism Aswan at ~ 18 km from the piezometer using NUMIS Auto magnetic resonance sounding system following the MRS site common practices. The MRS data were subjected to adaptive notch filtering to remove spiky noise and potential power line harmonics, and finally enhance the MRS signals at the selected stations. The water depth (~75–95 m) estimated from the MRS data around the piezometer largely coincides with the actual depth (~80 m) at the piezometer itself. A transient electromagnetic (TEM) profile performed near the measured MRS stations indicates no water content. Unlike the MRS data for NNL-1, the data for NNL-2 and NNL-3 match well with TEM data at shallow depths. The evaluation generally indicates that the MRS method is applicable for groundwater exploration at the selected piezometer and the other three stations. However, the method necessarily requires additional denoising efforts to comprehensively characterise groundwater occurrence and aquifer parameters.

Effects of base frequency, duty cycle, and waveform repetition on transient electromagnetic responses: Insights from models of a deep-buried conductor

In transient electromagnetic (TEM) methods, the recorded data for the magnetic field and its time derivative are influenced by the waveforms of the transmitting current. These waveforms are characterized by parameters such as the durations of the turn-on, steady, and turn-off stages as well as the base frequency, duty cycle, and waveform repetition. The full-waveform effects encompass all counteracting effects observed in TEM responses that arise from using a waveform other than a step-off or step-on waveform. To address these complexities, a 3D forward-modeling solver is developed for TEM methods, capable of calculating TEM responses while considering realistic waveforms with considerations for base frequency, duty cycle, and waveform repetition. The effects of these parameters on the magnetic field and its time derivative are investigated through numerical experiments using models involving deep-buried horizontal and vertical conductors. The results indicate that as the waveform period increases, there is a significant improvement in the detection and discrimination capability of the magnetic field for perfect conductors with a conductivity of 1000 S/m or higher in the models presented in this study. However, the improvement in the time derivative of the magnetic field is minimal. Interestingly, the improvement in the magnetic field and its time derivative does not indicate a consistent trend as the conductor conductivity increases. In addition, a novel equivalence phenomenon is uncovered in the magnetic field and its time derivative data. The results also suggest that as the number of waveform repetitions increases or the duty cycle decreases, the discrepancies between the recorded TEM responses at the two off-time stages within a single period gradually diminish. However, the influence of the full waveform on the TEM responses becomes more pronounced with an increasing number of waveform repetitions or a decrease in the duty cycle, especially in terms of the magnetic field.

Exploring the challenges of interpreting near-surface towed transient electromagnetic data on saline permafrost

We analyze a continuous near-surface transient electromagnetic (TEM) survey on permafrost around Ilulissat in western Greenland, an area characterized by continuous saline permafrost. The TEM data are severely affected by induced polarization (IP) effects, causing a large range of decay shapes. We identify seven unique decay shapes: oversteepened, sign change, all negative, double sign change, no apparent IP, flat spot, and positive nonmonotonic, the last two of which have not previously been identified in the scientific literature. A clear spatial dependency of the decay shapes is observed. Inversion of the data is carried out using a Cole-Cole model that poses a highly nonunique inversion problem with an extreme starting-model dependency. A series of inversion and forward-modeling experiments demonstrate these challenges and indicate that a low-resistivity and highly chargeable layer with time constant ([Formula: see text]) values between [Formula: see text] and [Formula: see text] and frequency exponent ( C) values above 0.74 is needed to fit the data used in the inversion experiment. Forward modeling further indicates that low [Formula: see text] and high C values are needed to reproduce the observed flat-spot and positive-nonmonotonic decay shapes. Based on these observations, we attribute the IP effects to the orientational polarization of ice in the soil column. This mechanism allows for low-resistivity, high-chargeability layers due to partially frozen saline sediments, a combination that is difficult to explain by the IP mechanisms traditionally considered. Based on the forward modeling of a realistic structural model of a layered sediment pack over dipping bedrock, we are able to reproduce all the observed decay shapes. This model provides a consistent framework for qualitative interpretation of the entire data set and evidence for the presence of saline deposits in the central parts of the sedimentary basins.

Bayesian inversion and uncertainty analysis of transient electromagnetic data

AbstractQuantification of non‐uniqueness and uncertainty is important for transient electromagnetism (TEM). To address this issue, we develop a trans‐dimensional Bayesian inversion schema for TEM data interpretation. The trans‐dimensional posterior probability density (PPD) offers a solution to model selection and quantifies parameter uncertainty resulting from the model selection from all possible models rather than determining a single model. We use the reversible‐jump Markov chain Monte Carlo sampler to draw ensembles of models to approximate PPD. In addition to providing reasonable model selection, we address the reliability of the inversion results for uncertainty analysis. This strategy offers reasonable guidance when interpreting the inversion results. We make the following improvements in this paper. First, in terms of algorithmic acceleration, we use the nonlinear optimization inversion results as the initial model and implement the multi‐chain parallel method. Second, we develop double factors to control the sampling step size of the proposed distribution, so that the sampling models cover the high‐probability region of the parameter space as much as possible. Finally, we provide the potential scale reduction factor‐η convergence criteria to assess the convergence of the samples and ensure the rationality of the output models. The proposed methodology is first tested on synthetic data and subsequently applied to a field dataset. The TEM inversion results show that probability inversion can provide reliable references for data interpretation through uncertainty analysis.

Near-Surface Geophysical Characterization of a Marble Deposit to Promote a Sustainable Small-Scale Mining

Small-scale mining (SSM) is responsible for almost all the production of non-metallic minerals in the world and represents around 80% of the mining in Brazil. The lack of direct geological information increases the level of uncertainty associated with the exploratory process, compromises mine planning, limits mineral extraction, and contributes to maximizing environmental issues. In this research, near-surface geophysical methods, including Electrical Resistivity, Capacitive Resistivity, Ground Penetrating Radar (GPR), and Transient Electromagnetic (TEM), were applied to characterize a marble deposit in an SSM located in the Campos do Jordão region, São Paulo state, southeast Brazil. The geophysical methods used provide indirect information about the subsurface geology based on the contrast in electrical and electromagnetic properties. Resistivity results show the efficiency of locating marble deposits, as well as fracture zones. GPR profiles allowed for the investigation of the structural heterogeneities in the subsurface. Geophysical data and lithological information from drill holes were integrated into Micromine software and guided the development of a geological model and a conceptual pit model. The information inferred from the pit modeling allowed us to analyze the potential of the deposit and should be used to assist in developing sustainable mining planning. The results of this work demonstrate that the investment in geophysical research can support the modernization of an SSM and contribute to more sustainable and productive mining.

A conditioned Latin hypercube sampling design methodology for ground-truthing transient EM resistivity models

Translating 3D resistivity models of the subsurface to relevant hydrogeological parameters such as sediment texture requires applying a rock physics transform that is ideally calibrated using co-located 1D sediment texture logs. It is difficult to select optimal logging sites in the model that span the range of observed resistivity values using expert opinion alone. Instead, we demonstrate a methodology for selecting logging sites using conditioned Latin hypercube sampling (cLHS), a maximally stratified sampling design methodology. cLHS employs simulated annealing to select sampling sites that minimize an objective function by reproducing the statistics of the resistivity model. We implement numerous preprocessing steps to isolate desirable regions within the model for logging, and we use principal component analysis to reduce the dimensionality of the model to aid the convergence of the simulated annealing procedure while preserving the majority of the model variance. Additionally, we impose an additional objective function condition to penalize the selection of sites in the resistivity model that exhibit a high degree of lateral variance. This condition is used to increase the likelihood that the sediment texture logs acquired at the selected sites correlate well with the resistivity model. We tested our approach on a resistivity model obtained using a towed transient EM system in an almond orchard in California’s Central Valley, comparing our selected sets of logging sites against the five cone penetrometer testing logs originally acquired at the site, as well as randomly selected logging sites. The sample sets selected with cLHS were consistently able to reproduce the statistics of the resistivity model and did so more effectively than the true sampling locations or the random sampling locations.

Sensitivity analysis of inverted model parameters from transient electromagnetic measurements affected by induced polarization effects

We investigate the application of the distance-based global sensitivity analysis (DGSA) to evaluate the sensitivity of electrical model parameters obtained from transient electromagnetic (TEM) data including induced polarization (IP) effects. We propose novel open-source forward modeling and inversion routines for single-loop TEM data including IP effects with the maximum phase angle model to model the frequency dependence of the complex resistivity. In a first step, we evaluate the accuracy of our forward modeling and inversion routines using numerical studies, where the actual variations of layer thicknesses and resistivities, as well as the frequency dependence of the complex resistivity is known. In a second step, we extend our investigation to field data and apply our approach to three distinct case studies in layered media: 1) a confined aquifer corresponding to conductive non-polarizable media, 2) a graphite deposit corresponding to highly conductive and polarizable anomalies in a resistive host rock and 3) an ice glacier corresponding to highly resistive polarizable media. Our DGSA results reveal that standard depth of investigation (DOI) approaches may overestimate the true sensitivity of the model obtained from the inversion. TEM data collected in conductive layered media without IP effects show reduced sensitivity above the predicted DOI. The case studies in polarizable media demonstrate that the maximum phase angle is more influential on the TEM model response than the relaxation time and dispersion coefficient. Our DGSA results for polarizable media reveal that TEM field data collected at the graphite deposit and at the ice glacier are sensitive to the geometry of the polarizable layer.

Improved adaptive regularization for simulated annealing inversion of transient electromagnetic

Geophysical inversion usually involves ill-posed problem. Regularization is the most commonly used method to mitigate this problem. There are many regularization parameter selection methods, among which the adaptive regularization method can automatically update parameters during iteration, reducing the difficulty of parameter selection. Therefore, it is widely used in linear inversion. However, there are very few studies on the use of adaptive regularization methods in stochastic optimization algorithms. The biggest difficulty is that in stochastic optimization algorithms, the search direction of any iteration is completely random. Data fitting term and stabilizing term vary in a wide range, making it difficult for traditional methods to work. In this paper, we consider the contributions of the data fitting term and the stabilizing term in the objective function and give an improved adaptive regularization method for very fast simulated annealing (VFSA) inversion for transient electromagnetic (TEM) data. The optimized method adjusts the two terms dynamically to make them in balance. We have designed several numerical experiments, and the experimental results demonstrate that the method in this paper not only accelerates the convergence, but also the inversion results are very little affected by the initial regularization parameter. Finally, we apply this method to field data, and the inversion results show very good agreements with nearby borehole data.

Quasi-2D inversion of surface large fixed-loop transient electromagnetic sounding data

Abstract In many cases, 1D inversion is still an important step in transient electromagnetic data processing. Potential issues may arise in the calculation of apparent resistivity using induced electromotive force (EMF) due to overshoot and the presence of multi-valued functions. Obtaining reliable and consistent inversion results using a uniform half-space as the initial model is challenging, especially when aiming for efficient inversion. Focusing on these problems, we use the land-based transient electromagnetic (TEM) sounding data, which was acquired by using a large fixed-loop transmitter, and adopt a quasi-2D inversion scheme to generate improved images of the subsurface resistivity structure. First, we have considered directly using magnetic field data or converting induced EMF into magnetic field, and then calculating the apparent resistivity over the whole zone. Next, a resistivity profile that varies with depth is obtained through fast smoke ring imaging. This profile serves as the initial model for the subsequent optimal inversion. The inversion scheme uses a nonlinear least-squares method, incorporating lateral and vertical constraints, to produce a quasi-2D subsurface image. The potentiality of the proposed methodology has been exemplified through the interpretation of synthetic data derived from a 3D intricate resistivity model, as well as field data obtained from a TEM survey conducted in a coalmine field. In both cases, the inversion process yields quasi-2D subsurface images that exhibit a reasonable level of accuracy. These images appear to be less moulded by 3D effects and demonstrate a satisfactory level of agreement with the known target area.

An optimized and hybrid gating scheme for the suppression of very low-frequency radios in transient electromagnetic systems

Abstract. One of the most widely used approaches for measuring the earth's subsurface resistivity is the transient electromagnetic (TEM) method. However, noise and interference from different sources, e.g., radio communication, the instrument, the atmosphere, and power lines, severely taint these types of signals. In particular, radio transmission in the very low-frequency (VLF) range between 3 and 30 kHz is one of the most prominent sources of noise. Transient electromagnetic signals are normally gated to increase the signal-to-noise ratio. A precise selection of gate shapes is required to suppress undesired noise while allowing the TEM signal to pass unaltered. We employ the multi-objective particle swarm optimization technique to choose optimal gate shapes and placements by minimizing an objective function composed of standard error bars, the covariance between gates, and the distortion of the gated signal. The proposed method is applied to both fully sampled synthetic TEM data and to boxcar-gated field data. The best output from the search space of gate shapes was found to be a hybrid combination of boxcar and Hamming gates. The effectiveness of hybrid gating over traditional boxcar and semi-tapered gating is confirmed by an analysis of covariance matrices and error bars. The results show that the developed method effectively suppresses VLF noise in the middle gates, which are gates with center times spanning 30 to 200 µs , and in the late gates, which are gates with center times spanning 200 to 1130 µs. The analysis shows that the average improvement in standard errors obtained for the hybrid gating scheme over boxcar gating is 1.719 and 1.717 for synthetic and field data, respectively.