Analysis of Long-offset Transient Electromagnetic (LOTEM) Data in Time, Frequency, and Pseudo-seismic Domain

2018 ◽  
Vol 23 (1) ◽  
pp. 15-32
Author(s):  
Muhammad Younis Khan ◽  
Guo-qiang Xue ◽  
Wei-ying Chen ◽  
Hua-sen Zhong

Long-offset transient electromagnetic (LOTEM) has received great attention in mineral, hydrocarbon and hydrogeological investigations for the last several years. Conventionally, TEM soundings have been presented as apparent resistivity curves as function of time. With development of sophisticated computational algorithms, it became possible to extract more realistic geoelectric information by applying inversion programs to 1-D and 3-D problems. Here, we analyze LOTEM data by carrying out analysis in time, frequency and pseudo-seismic domain supported by borehole information. At first, H, K, A & Q type geoelectric models are processed using a proven inversion program (1-D Occam inversion). Second, time-to-frequency transformation is conducted from TEM ρ a ( f) curves for the same models based on all-time apparent resistivity curves. Third, 1-D Bostick's algorithm was applied to the transformed resistivity. Finally, EM diffusion field is transformed into propagating wave field and constructed pseudo-seismic section. The transformed seismic-like wave indicates that some reflection and refraction phenomena appear when the EM wave field interacts with a geoelectric interface at different depth intervals due to contrast in resistivity. In all three cases, synthetic tests showed that conductive anomaly in resistive host environment can be retrieved more clearly than resistive target. A case study illustrates the successful application of proposed approach in recovering a water-filled mined-out area in a coal field located in Ye county, Henan province, China. The results support the introduction of pseudo-seismic imaging technology in long-offset version of TEM which can also be an useful aid if integrated with seismic reflection technique to explore possibilities for high resolution EM imaging in future. [Figure: see text]

Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. B59-B76 ◽  
Author(s):  
Amir Haroon ◽  
Klaus Lippert ◽  
Vladimir Mogilatov ◽  
Bülent Tezkan

The marine differential electric dipole (DED) is applied for the first time to study a subseafloor groundwater body in the coastal region of Bat Yam, Israel. Previous marine long-offset transient electromagnetic applications detected this freshwater body underneath the Mediterranean seafloor. We have applied the novel DED method for the first time in the marine environment to further investigate this natural phenomenon. The main objectives are to locate the freshwater-seawater interface at the western aquifer edge and to identify the mechanism controlling this freshwater occurrence beneath the seafloor. The acquired step-on signals allow one to detect the freshwater body in the vicinity of the Israeli coastline at a depth of approximately 70 m beneath the seafloor. However, aquifer thickness is only poorly determined and may vary between 40 and 100 m. A lateral resistivity contrast is observable between adjacent 1D inversion models and also apparent in data profile curves that constrain the seaward extent of the detected resistive body to a distance of less than 4 km from the coastline. A subsequent 2.5D forward-modeling study aims to find a subseafloor resistivity distribution that adequately explains all measured DED data simultaneously. The results further constrain the lateral extent of the resistive aquifer to approximately 3.6–3.7 km from the Israeli coast. Furthermore, the data indicate that the aquifer system may be susceptible to seawater intrusion, as a superior data fit is achieved if a brackish water zone of approximately [Formula: see text] with a lateral extent of less than 300 m is located at the head of the freshwater body.


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 ◽  
2007 ◽  
Vol 72 (1) ◽  
pp. F25-F34 ◽  
Author(s):  
Benoit Tournerie ◽  
Michel Chouteau ◽  
Denis Marcotte

We present and test a new method to correct for the static shift affecting magnetotelluric (MT) apparent resistivity sounding curves. We use geostatistical analysis of apparent resistivity and phase data for selected periods. For each period, we first estimate and model the experimental variograms and cross variogram between phase and apparent resistivity. We then use the geostatistical model to estimate, by cokriging, the corrected apparent resistivities using the measured phases and apparent resistivities. The static shift factor is obtained as the difference between the logarithm of the corrected and measured apparent resistivities. We retain as final static shift estimates the ones for the period displaying the best correlation with the estimates at all periods. We present a 3D synthetic case study showing that the static shift is retrieved quite precisely when the static shift factors are uniformly distributed around zero. If the static shift distribution has a nonzero mean, we obtained best results when an apparent resistivity data subset can be identified a priori as unaffected by static shift and cokriging is done using only this subset. The method has been successfully tested on the synthetic COPROD-2S2 2D MT data set and on a 3D-survey data set from Las Cañadas Caldera (Tenerife, Canary Islands) severely affected by static shift.


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.


Author(s):  
Jaroslav Smutny ◽  
Viktor Nohal ◽  
Daniela Vukusicova ◽  
Herbert Seelmann

This paper deals with description and application of the Wigner-Ville transformation for vibration analysis. This transformation belongs to the group of non-linear time-frequency processes. Thanks to its properties, it may be successfully used in the area of non-stationary and transitional signals describing various natural processes. The use in the field of the railway constructions testing represents a quite an interesting application area of the transformation. This paper contains mathematical analysis of the transformation, a case study and practical experience obtained and recommendations for its practical use.


Sign in / Sign up

Export Citation Format

Share Document