The Development and Applications of the Helicopter-borne Transient Electromagnetic System CAS-HTEM

2019 ◽  
Vol 24 (4) ◽  
pp. 653-663 ◽  
Author(s):  
Xin Wu ◽  
Guangyou Fang ◽  
Guoqiang Xue ◽  
Lihua Liu ◽  
Leisong Liu ◽  
...  

Over the past decade, helicopter-borne transient electromagnetic (HTEM) systems have been rapidly developed. A new HTEM prototype (referred to as a CAS-HTEM) has been developed by the Chinese Academy of Sciences. In terms of hardware, the CAS-HTEM system uses an inflatable structure to carry the transmitting loop, which significantly reduces the weight of the system and makes it easier to transport. A dual gain receiver was innovated to extend the dynamic range of the system. In addition, an observation circuit for transmitting voltage waveform is introduced, so that the derivative waveform of transmitting current with higher SNR could be calculated. In terms of data processing, more reliable early time data could be obtained by band-limited effect removal; a weighted stacking algorithm is introduced to reduce the narrow band noise more effectively and increase the sensitivity of data to the anomaly location; a method based on τ-domain transform is used for late time signal processing. The results of the field test which was carried out in Inner Mongolia were found to be consistent with the drill data, which effectively verified the performance of this HTEM prototype.

Geophysics ◽  
1992 ◽  
Vol 57 (9) ◽  
pp. 1116-1126 ◽  
Author(s):  
James E. Hanneson

An algorithm for computing the transient electromagnetic (TEM) response of a dipping plate in a conductive half‐space has been developed. For a stationary [Formula: see text] current loop source, calculated profiles simulate the response of the University of Toronto electromagnetic system (UTEM) over a plate in a 1000 Ω ⋅ m half‐space. The objective is to add to knowledge of the galvanic process (causing poloidal plate currents) and the local induction process (causing toroidal currents) by studying host and plate currents with respect to surface profiles. Both processes can occur during TEM surveys. Plates are all [Formula: see text] thick with various depths, dips, and conductances. Calculated host and plate currents provide quantitative examples of several effects. For sufficiently conductive plates, the late time currents are toroidal as for a free‐space host. At earlier times, or at all times for poorly conducting plates, the plate currents are poloidal, and the transitions to toroidal currents, if they occur, are gradual. At very late times, poloidal currents again dominate any toroidal currents but this effect is rarely observed. Stripped, point‐normalized profiles, which reflect secondary fields caused by the anomalous plate currents, illustrate effects such as early time blanking (caused by noninstantaneous diffusion of fields into the target), mid‐time anomaly enhancement (caused by galvanic currents), and late time plate‐in‐free‐space asymptotic behavior.


Geophysics ◽  
2001 ◽  
Vol 66 (4) ◽  
pp. 1090-1096 ◽  
Author(s):  
Z. Zhang ◽  
J. Xiao

Large‐loop electromagnetic (EM) systems that measure transient EM (TEM) data on the surface or in boreholes have shown increased application in exploration geophysics. Accurate interpretation of borehole TEM data is necessary to discover deep hidden targets that cannot be detected with surface systems. However, the inversion of borehole TEM data has not been fully addressed. In this paper, we study the propagation of the TEM field from a large‐loop EM borehole system inside a layered earth and develop a new inversion algorithm to reconstruct layered conductivity structures from large‐loop TEM data measured with both surface and borehole configurations. The magnetic field and sensitivities are first computed in the frequency domain and are then transformed into the time domain where the inversion is performed. The surface data have a higher S/N ratio at early time channels, while the borehole data have a higher S/N ratio at late time channels. Consequently, the surface data can be inverted to better resolve shallow structures, and the borehole data can be used to better detect deep structures. The merits of joint inversions of borehole and surface data are explored. We test our inversion algorithm using numeric examples.


Geophysics ◽  
1986 ◽  
Vol 51 (7) ◽  
pp. 1462-1471 ◽  
Author(s):  
Brian R. Spies ◽  
Dwight E. Eggers

Problems and misunderstandings arise with the concept of apparent resistivity when the analogy between an apparent resistivity computed from geophysical observations and the true resistivity structure of the subsurface is drawn too tightly. Several definitions of apparent resistivity are available for use in electromagnetic methods; however, those most commonly used do not always exhibit the best behavior. Many of the features of the apparent resistivity curve which have been interpreted as physically significant with one definition disappear when alternative definitions are used. It is misleading to compare the detection or resolution capabilities of different field systems or configurations solely on the basis of the apparent resistivity curve. For the in‐loop transient electromagnetic (TEM) method, apparent resistivity computed from the magnetic field response displays much better behavior than that computed from the induced voltage response. A comparison of “exact” and “asymptotic” formulas for the TEM method reveals that automated schemes for distinguishing early‐time and late‐time branches are at best tenuous, and those schemes are doomed to failure for a certain class of resistivity structures (e.g., the loop size is large compared to the layer thickness). For the magnetotelluric (MT) method, apparent resistivity curves defined from the real part of the impedance exhibit much better behavior than curves based on the conventional definition that uses the magnitude of the impedance. Results of using this new definition have characteristics similar to apparent resistivity obtained from time‐domain processing.


1999 ◽  
Vol 2 (03) ◽  
pp. 271-280 ◽  
Author(s):  
Ekrem Kasap ◽  
Kun Huang ◽  
Than Shwe ◽  
Dan Georgi

Summary The formation-rate-analysis (FRASM) technique is introduced. The technique is based on the calculated formation rate by correcting the piston rate with fluid compressibility. A geometric factor is used to account for irregular flow geometry caused by probe drawdown. The technique focuses on the flow from formation, is applicable to both drawdown and buildup data simultaneously, does not require long buildup periods, and can be implemented with a multilinear regression, from which near-wellbore permeability, p * and formation fluid compressibility are readily determined. The field data applications indicate that FRA is much less amenable to data quality because it utilizes the entire data set. Introduction A wireline formation test (WFT) is initiated when a probe from the tool is set against the formation. A measured volume of fluid is then withdrawn from the formation through the probe. The test continues with a buildup period until pressure in the tool reaches formation pressure. WFTs provide formation fluid samples and produce high-precision vertical pressure profiles, which, in turn, can be used to identify formation fluid types and locate fluid contacts. Wireline formation testing is much faster compared with the regular pressure transient testing. Total drawdown time for a formation test is just a few seconds and buildup times vary from less than a second (for permeability of hundreds of millidarcy) to half a minute (for permeability of less than 0.1 md), depending on system volume, drawdown rate, and formation permeability. Because WFT tested volume can be small (a few cubic centimeters), the details of reservoir heterogeneity on a fine scale are given with better spatial resolution than is possible with conventional pressure transient tests. Furthermore, WFTs may be preferable to laboratory core permeability measurements since WFTs are conducted at in-situ reservoir stress and temperature. Various conventional analysis techniques are used in the industry. Spherical-flow analysis utilizes early-time buildup data and usually gives permeability that is within an order of magnitude of the true permeability. For p* determination, cylindrical-flow analysis is preferred because it focuses on late-time buildup data. However, both the cylindrical- and spherical-flow analyses have their drawbacks. Early-time data in spherical-flow analysis results in erroneous p* estimation. Late-time data are obtained after long testing times, especially in low-permeability formations; however, long testing periods are not desirable because of potential tool "sticking" problems. Even after extended testing times, the cylindrical-flow period may not occur or may not be detectable on WFTs. When it does occur, permeability estimates derived from the cylindrical-flow period may be incorrect and their validity is difficult to judge. New concepts and analysis techniques, combined with 3-D numerical studies, have recently been reported in the literature.1–7 Three-dimensional numerical simulation studies1–6 have contributed to the diagnosis of WFT-related problems and the improved analysis of WFT data. The experimental studies7 showed that the geometric factor concept is valid for unsteady state probe pressure tests. This study presents the FRA technique8 that can be applied to the entire WFT where a plot for both drawdown and buildup periods renders straight lines with identical slopes. Numerical simulation studies were used to generate data to test both the conventional and the FRA techniques. The numerical simulation data are ideally suited for such studies because the correct answer is known (e.g., the input data). The new technique and the conventional analysis techniques are also applied to the field data and the results are compared. We first review the theory of conventional analysis techniques, then present the FRA technique for combined drawdown and buildup data. A discussion of the numerical results and the field data applications are followed by the conclusions. Analysis Techniques It has been industry practice to use three conventional techniques, i.e., pseudo-steady-state drawdown (PSSDD), spherical and cylindrical-flow analyses, to calculate permeability and p* Conventional Techniques Pseudo-Steady-State Drawdown (PSSDD). When drawdown data are analyzed, it is assumed that late in the drawdown period the pressure drop stabilizes and the system approaches to a pseudo-steady state when the formation flow rate is equal to the drawdown rate. PSSDD permeability is calculated from Darcy's equation with the stabilized (maximum) pressure drop and the flowrate resulting from the piston withdrawal:9–11 $$k {d}=1754.5\left({q\mu \over r {i}\Delta p {{\rm max}}}\right),\eqno ({\rm 1})$$where kd=PSSDD permeability, md. The other parameters are given in Nomenclature.


Geophysics ◽  
1986 ◽  
Vol 51 (6) ◽  
pp. 1291-1297 ◽  
Author(s):  
Yang Sheng

Early‐time and late‐time apparent resistivity approximations have been widely used for interpretation of long‐offset transient electromagnetic (LOTEM) measurements because it is difficult to find a single apparent resistivity over the whole time range. From a physical point of view, Dr. C. H. Stoyer defined an apparent resistivity for the whole time range. However, there are two problems which hinder its use: one is that there is no explicit formula to calculate the apparent resistivity, and the other is that the apparent resistivity has no single solution. A careful study of the two problems shows that a numerical method can be used to calculate a single apparent resistivity. A formula for the maximum receiver voltage over a uniform earth, when compared with the receiver voltage for a layered earth, leads to the conclusion that, in some cases, a layered earth can produce a larger voltage than any uniform earth can produce. Therefore, our apparent resistivity definition cannot be applied to those cases. In some other cases, the two possible solutions from our definition do not merge, so that neither of them is meaningful for the whole time range.


Geophysics ◽  
2012 ◽  
Vol 77 (3) ◽  
pp. E203-E212 ◽  
Author(s):  
Nicklas Skovgaard Nyboe ◽  
Kurt Sørensen

The transient electromagnetic method (TEM) is a recognized tool for determining the subsurface resistivity structure over a wide depth interval. A key requirement for an accurate characterization of the shallow part of this interval is a sufficiently wideband receiver system. The maximum depth of investigation, on the other hand, is determined by the late-time signal-to-noise ratio (S/N). It has been demonstrated that the use of compact wideband receiver coils tends to deteriorate the late-time S/N due to their inherent low sensitivity to time-varying magnetic fields. To overcome this problem, we used a bandwidth- and sensitivity-optimized parallel recording setup using two separated receiver coils optimized for measuring the early- and late-time part of the transient signal, respectively. Using this setup, we experienced substantial improvements in the achievable depth of investigation while retaining shallow resolution. Further noise reduction may be obtained by adapting the applied synchronous detection measuring scheme to the local noise conditions. This is especially important when TEM are carried out in culturally developed areas, which is often the case for environmental investigations. Such soundings frequently suffer from inferior signal quality due to EM noise from local cultural sources. Comparing the frequency-domain filtering effect of the applied synchronous detection measuring scheme with the actual noise frequency content allowed us to evaluate the relative importance of the many sources of noise encountered in the field. Further improvements could be obtained by tailoring the applied synchronous detection to better reject specific noise frequencies through modification of the properties of time-gate integration and gate stacking.


Geophysics ◽  
2020 ◽  
Vol 86 (1) ◽  
pp. E1-E11
Author(s):  
Peter K. Fullagar ◽  
Ralf Schaa

In the resistive limit, discrete conductors give rise to magnetic dipole fields. Magnetostatic modeling of time integrals, or moments, of transient electromagnetic (TEM) data therefore offers a means for fast approximate 3D modeling and inversion of TEM data sets. In our approximate inversion scheme, the net TEM moment response is represented as the combination of discrete conductor and uniform host responses. The inversion algorithm first estimates a homogeneous host conductivity, and then it subtracts the host response and fits the residual moment data by adjusting the conductivities of cells comprising a 3D rectangular mesh. To expedite calculation of the host response, we have derived analytic formulas for first-order TEM moments produced on and under the ground by a horizontal electric dipole on the surface of a homogeneous conductive half-space. We present analytic expressions for idealized all-time “complete” moments, or resistive limits, as well as for realizable finite-time “incomplete” moments. The moments produced by an arbitrary horizontal polygonal loop are determined by combining contributions from appropriately oriented electric dipoles. Downhole TEM moments computed with the new expressions reveal substantial differences between incomplete and complete moments when early time data are excluded and between step and impulse response incomplete moments. The role of the formulas in the first stage of our moment-based 3D inversion scheme is illustrated via analysis of downhole TEM data recorded at Santander, Peru. The host conductivity of best fit for early time B-field moments is 2.40 mS/m, consistent with apparent conductivities derived from ground TEM data recorded in the same area.


Geophysics ◽  
1985 ◽  
Vol 50 (6) ◽  
pp. 1002-1019 ◽  
Author(s):  
Peggie R. Gallagher ◽  
Stanley H. Ward ◽  
G. W. Hohmann

The computer program PLATE, developed at the University of Toronto, models the electromagnetic (EM) response of an inductively thin plate in free space. We used PLATE to compute two components of the time derivative of the magnetic field for a range of models for the EM37 fixed‐source transient system ([Formula: see text] loop). Analysis of the response curves produced methods of interpretation for obtaining plate geometry and conductance. The overall width of an anomaly, the distance between peaks and the width of the updip lobes, can provide an estimate of depth. Dip has the dominant effect on the ratio of the peak amplitudes. A rough estimate of plate size and the position in time (early or late) of the currents is essential before proceeding with interpretation. Strike length is not obviously reflected in the shape of the curves, but depth extent is indicated by the rate at which the downdip tail returns to the baseline, except for vertical plates. For vertical plates, curve matching may be the only method of obtaining an estimate of depth extent. Varying conductance for a particular model in free space affects whether a channel represents an early, intermediate, or late time response. The shape of a profile varies with the time of measurement. The estimated time constant can be used to calculate the conductance, provided an estimate of the shortest dimension of the plate is available. Extinction angles appear frequently for plates of small depth extent but do not occur for plates which are of infinite strike and depth extent with respect to the size of the transmitting loop.


2019 ◽  
Vol 24 (3) ◽  
pp. 351-359
Author(s):  
Junjie Xue ◽  
Jiulong Cheng ◽  
Guoqiang Xue ◽  
Hai Li ◽  
Dongyang Hou ◽  
...  

The diffusive electromagnetic field can be transformed into the wave domain by means of mathematical conversion. The transformed field can then be interpreted with the tools in seismic data processing so that the identification to the underground targets can be effectively improved. However, the conversion is typically an ill-posed problem that needs to be solved using regularization tools. Based on the conventional regularization with smooth constraints in the L2 norm, the inversion result is of low resolution, while that obtained using truncated singular value decomposition (TSVD) methods is typically accurate, but has poor stability. To obtain a stable and accurate transformed electromagnetic field value, this study proposed to combine conventional regularization tools and singular value decomposition algorithms by incorporating a set of weighting coefficients. The proposed method is validated on both synthetic and observed data. The results from the proposed method are more accurate at the early time, and at the late time are more stable compared with the other methods. Furthermore, the example of field data shows that the proposed method could potentially further improve the interpretation accuracy of future mining explorations.


SPE Journal ◽  
2015 ◽  
Vol 20 (04) ◽  
pp. 717-728 ◽  
Author(s):  
Francisco J. Pacheco-Roman ◽  
S. Hossein Hejazi

Summary Solubility and diffusivity of gases in heavy oils, quantified by Henry's constant (Hij) and the diffusion coefficient (D), respectively, are essential properties for the design of recovery processes that require the injection of gas or vapor solvents into the reservoir. Data, obtained from various experimental procedures such as the pressure-decay technique (PDT), are used to estimate these two parameters. The PDT uses a pressure/volume/temperature (PVT) cell where the gas-phase pressure declines as gas diffuses into heavy oil following an early- and a late-time regime. Current approaches to analyze data from the conventional PDT are either graphical techniques based on early-time data or full numerical simulation. Early-time data, when the diffusing gas has not reached the bottom of the PVT cell, do not provide enough information to simultaneously estimate both the diffusion coefficient and Henry's constant. Hence, existing graphical procedures are limited to diffusion-coefficient estimation. In this paper, we propose a novel graphical technique to estimate the diffusion coefficient and Henry's constant by use of the late-time data from pressure-decay experiments. Our method is derived from the modeling of gas-phase pressure decay by use of Fick's second law and gas-phase mass-balance equations. We use the integral method to provide an approximate analytical solution to the set of equations. In addition, by use of the resultant solution, we develop a simple graphical method to directly estimate both the diffusion coefficient and Henry's constant. The estimated parameters through the proposed technique for methane/bitumen and carbon dioxide/bitumen experiments are in close agreement with those reported in the literature.


Sign in / Sign up

Export Citation Format

Share Document