Controlled Source Radiomagnetotelluric Research Articles

On the Use of Different Sources for Controlled Source Radiomagnetotellurics

On the Use of Different Sources for Controlled Source Radiomagnetotellurics

Study on the detectability of CSRMT in archaeological exploration: taking the ancient city wall remains as an example

Although electromagnetic methods have made significant advances in archaeological exploration, there remains a lack of a method that offers both high operating efficiency and high resolution. This study introduces a novel approach derived from the Controlled Source Radio Magnetotelluric (CSRMT) method. The CSRMT method combines the efficient configuration of the Controlled Source Audio-magnetotelluric (CSAMT) method with enhanced resolution capabilities. To evaluate the efficacy of the new method in archeological exploration, we designed a series of models based on existing references related to ancient city wall remains. Using the 3D contraction integral equation method, we numerically simulated the high-frequency component of the CSRMT method, and then on, the Laterally Constrained Inversion technique is used to invert the data. The effectiveness of the method was assessed qualitatively and quantitatively through analysing the electromagnetic field and apparent resistivity of models with varying widths, heights, top burial depths, and resistivity. Our analysis revealed that models with shallower burial depth and greater height of the underground city wall exhibited higher resolution. The findings demonstrate that the CSRMT method can efficiently and accurately identify buried remains, making it a highly promising tool for archaeological exploration and warranting further development.

Use of Subharmonics of Base Frequencies in the CSRMT Method with Loop Sources

In the controlled source radiomagnetotelluric (CSRMT) sounding method, a horizontal magnetic dipole, HMD (vertical loop) or a horizontal electric dipole, and HED (grounded line) are used as sources. When working with HMD, the source is usually tuned to resonance to increase the current in the loop. However, the disadvantage of this approach is the narrow frequency range realized in the CSRMT method (1–12 kHz) and the short operating distance from the source (600–800 m). The need to tune the source to resonance at each selected frequency reduces the efficiency of the survey. In the case of using HED for sounding, measurements are performed in a wider frequency range of 1 to 1000 kHz, and along with the signal of the base frequency, its subharmonics are measured. In this case, emitted signal measurements are possible at a distance of up to 3–4 km from the source. At the same time, the disadvantage of using HED is that it requires grounding, the arrangement of which requires additional time when working on frozen ground or dry stony soil. We consider the possibilities of generation and registration of signals of subharmonics of base frequencies when applying the CSRMT method with loop sources—HMD and VMD (horizontal loop). A matching unit (MU) based on a step-up transformer was developed, which increases the output voltage of the CSRMT transmitter. In a field test with base frequencies of 20, 40, and 80 kHz, the signal amplitudes increased by a factor of two to four for subharmonics at frequencies of 60–200 kHz and by up to 10–13 times for subharmonics at frequencies of 200–500 kHz due to transformation of signal spectrum provided by MU. The possibility of using odd subharmonics of base frequencies for inversion has been demonstrated in the results of field experiments with different sources (HED, HMD, and VMD). This expands the frequency range of the method when working with loop sources and increases the survey’s effectiveness. The use of loop sources in the CSRMT method is especially advantageous for winter work in Arctic regions.

Tensor CSRMT System with Horizontal Electrical Dipole Sources and Prospects of Its Application in Arctic Permafrost Regions

When studying horizontally-inhomogeneous media, it is necessary to apply tensor modifications of electromagnetic soundings. Use of tensor measurements is of particular relevance in near-surface electrical prospecting because the upper part of the geological section is usually more heterogeneous than the deep strata. In the Enviro-MT system designed for the controlled-source radiomagnetotelluric (CSRMT) sounding method, two mutually perpendicular horizontal magnetic dipoles (two vertical loops) are used for tensor measurements. We propose a variant of the CSRMT method with two horizontal electrical dipole sources (two transmitter lines). The advantage of such sources is an extended frequency range of 1–1000 kHz in comparison with 1–12 kHz of the Enviro-MT system, greater operational distance (up to 3–4 km compared to 600–800 m), and the ability to measure the signal at the fundamental frequency and its subharmonics. To implement tensor measurements with the equipment of the CSRMT method described in the paper, a technique of creating a time-varying polarization of the electromagnetic field (rotating field) has been developed based on the use of two transmitters with slightly different current frequencies and two mutually-perpendicular transmitter lines grounded at the ends. In this way, we made it possible to change the direction of the electrical and magnetic field polarization continuously. This approach allows realization of the technique of tensor measurements using the new modification of the CSRMT method. In permafrost areas, the hydrogenic taliks are widespread. These local objects are important in the context of study of environmental changes in the Arctic and can be successfully explored by the tensor CSRMT method. For the numerical modeling, a 2D model of the talik was used. Results of the interpretation of synthetic data showed the advantage of bimodal inversion using CSRMT curves of both TM and TE modes compared to separate inversion of TM and TE curves. These new data demonstrate the prospects of the tensor CSRMT method in the study of permafrost regions. The problems that can be solved using the CSRMT method in the Arctic permafrost regions are discussed.

3D controlled-source electromagnetic inversion in the radio-frequency band

The classical radio-magnetotelluric (RMT) method is nowadays routinely applied to various environmental, engineering, and exploration problems. The technique uses radio transmitters broadcasting in the frequency range of 10 kHz to 1 MHz, and the measurements are carried out in the far field. The well-known disadvantages of RMT are a lack of robust radio transmitters in remote areas; the absence of transmitters broadcasting below 10 kHz, which limits the penetration depth; and a possible low signal-to-noise ratio. To overcome these difficulties, controlled sources can be used (controlled-source RMT [CSRMT]). We extend the CSRMT method to perform measurements not only in the far field but also in the transition zone. In CSRMT practice, it often is challenging to maintain far-field conditions for logistical reasons. Therefore, part of the measured data contains signatures of the source field, which cannot be interpreted with magnetotelluric software. In addition, the source placed directly in the survey area allows us to increase the signal-to-noise ratio and resolution. Such CSRMT in the transition zone is, in fact, a controlled-source electromagnetic method but with full impedance tensor and tipper vector transfer functions. We develop new procedures for the 3D modeling and inversion of the tensor radio-frequency data measured in the transition zone of two perpendicular horizontal electric dipole sources. In this case, the geometry of the source must be considered in the forward modeling. The developed modeling and inversion software is tested on a synthetic 3D model. The 3D resistivity models derived from the real data confirm the geologic settings and are consistent with the available borehole information. Therefore, we conclude that the CSRMT approach extended to include the source field is feasible and that the developed procedures are reliable.

Application of the Controlled Source Radiomagnetotellurics (CSRMT) in the Study of Rocks Overlying Kimberlite Pipes in Yakutia/Siberia

Тhe task of searching for kimberlite pipes in covered areas of the Yakutia kimberlite province is very difficult due to the significant heterogeneity of the rocks overlying kimberlite pipes. The overlying strata of terrigenous sediments contain rocks of the trap complex (dolerite sills, tuff bodies). We consider the results of the controlled source radiomagnetotelluric (CSRMT) soundings in Yakutia/Siberia. Due to the great thickness of the overlying rocks (near 100 m) and the relatively small horizontal sizes of kimberlite pipes (80–200 m), they cannot confidently be detected directly. An additional difficulty in identifying pipe anomalies is the presence of a layer of conductive carbonaceous siltstones in the overlying strata. Therefore, the main aim of the CSRMT surveys was the study of overlying rocks and the search for indirect indications of the presence of pipes. Possibilities to study the structure of dolerite sills located within overlying sediments and to map the top edge of hosting carbonate rocks are demonstrated using the CSRMT data. The pinching out of dolerite sills above pipes («windows in traps») and the lowering of the top edge of hosting rocks at pipes can be considered as indirect indications of the presence of pipes.

One-dimensional Laterally Constrained Joint Anisotropic Inversion of CSRMT and ERT Data

In this paper, we discuss several approaches for a joint inversion of controlled source radiomagnetotelluric (CSRMT) and electrical resistivity tomography (ERT) data observed over anisotropic media. We compare results of 2D isotropic joint inversion with results of a newly developed joint 1D anisotropic inversion algorithm. The developed algorithm involves the full controlled source high frequency forward and inversion formulations without the plane wave assumption. We demonstrate that for measurements on an anisotropic subsurface the isotropic joint inversion cannot fit both datasets properly due to a high anisotropy of shallow horizons of quaternary sands and loams. The joint anisotropic inversion helps to solve this problem and highlights the advantages of a joint inversion of CSRMT and ERT data. We also demonstrate application of the laterally constrain algorithm for the anisotropic inversion. Results of the joint 1D anisotropic inversion of CSRMT and ERT data were successfully compared with existing borehole data.

Study of a Permafrost Area in the Northern Part of Siberia Using Controlled Source Radiomagnetotellurics

Controlled source radiomagnetotelluric (CSRMT) measurements with a horizontal electrical dipole were carried out for the planned construction of a railway bridge in the Northern part of Siberia (the Yamalo-Nenets region). The aim of this survey was the delineation of horizons of loams, which are mostly thixotropic, and the estimation of the depth of solid sands serving as a basement for bridge piles. The remote area of the investigation is characterized by a few number of broadcasting radio transmitters. Therefore, the possibilities of the conventional radiomagnetotelluric (RMT) method to study the geological structure are quite limited. Because of swamps, the construction area is only accessible in winter or early spring periods. Our CSRMT survey was carried out by using ungrounded (capacitive) electric lines. The measurements were conducted in the broadside configuration using a 580 m long grounded electrical line at a distance of about 640 m from the transmitter. At the first stage of the CSRMT data inversion, we used the 1D controlled source inversion code considering the geometry of the source and the effect of the transition zone. In addition, we selected frequencies associated with the far-field response of the electromagnetic field. To demonstrate the advantages of the controlled source, we compared the responses and inversion results obtained from both CSRMT and RMT data. The geoelectrical section obtained from CSRMT data is more detailed in the shallow part and allows to resolve the deep structures because of an extended frequency range and an increased depth of investigation. The conductivity model derived from the 2D inversion of the CSRMT data showed a good correlation with the borehole data. The CSRMT results were very useful for the construction of the railway bridge because of the reliable delineation of loams and sands horizons in the subsurface.

On the study of industrial waste sites on the Karelian Isthmus/Russia using the RMT and CSRMT methods

Radiomagnetotelluric (RMT) and Controlled Source Radiomagnetotelluric (CSRMT) measurements were carried out on landfills of industrial wastes located on the Karelian Isthmus/Russia. The RMT and CSRMT data were interpreted using the 2D inversion techniques, and the 2D conductivity structures of the investigated industrial waste sites were derived. The RMT method is based on measurements of electromagnetic fields of remote radio transmitters in a frequency range of 10–1000 kHz or in a frequency range of 1–1000 kHz by using a controlled source (CSRMT method). A landfill in the settlement Morozova near St. Petersburg was formed as a result of a long-term warehousing of toxic wastes by the sulfuric acid manufacture. Soils and ground waters were polluted with salts of heavy metals. In the adjacent to the landfill territory, investigations were realized by the RMT method and a polluted zone was determined. This zone is traced up to the coast of the Ladoga Lake and the Neva River. It constitutes a potential danger because the Neva River is used for the water supply of St. Petersburg. Another landfill of debris in a former sandy career is located in the Vyborg area. The high conductivity and a relative large thickness of the wastes did not allow us to map the bottom of the landfill by the RMT method using the frequency range of 10–1000 kHz. Therefore, CSRMT measurements with an extended frequency range up to 1 kHz were carried out on this landfill, and the 2D inversion of the CSRMT data clearly showed the lateral and vertical extension of it.

Mapping of buried faults using the 2D modelling of far-field controlled source radiomagnetotelluric data

Controlled source radiomagnetotellurics (CSRMT) is a relatively new geophysical method for near-surface applications. A rectangular signal with base frequencies between 0.1 and 150 kHz is injected through a grounded electric dipole which is used as a transmitter. Electric and magnetic field components are observed at these frequencies and at their subharmonics, usually in the far-field zone so that apparent resistivities and impedance phases can be obtained in a broad frequency range between 1 and 1000 kHz. Inline or broadside configuration can be used for measurements. Similar to the controlled source audiomagnetotelluric method, tensor measurements are also possible when locating two transmitters perpendicular to each other. A scalar CSRMT survey was carried out on the buried faults in the Vuoksa region, 110 km north of St. Petersburg to test the applicability of this method to the mapping of near-surface faults. A 700 m electric dipole with base frequencies of 0.5, 11.3, 30 and 105 kHz was used as a transmitter. Smooth apparent resistivity and phase values as a function of frequency from 1 kHz to 1 MHz were observed in the far-field zone for the inline configuration at 57 stations using a station distance of 20 m. Electric fields observed in the direction of the transmitter were perpendicular to the assumed strike direction of the buried faults so that they could be associated with the TM mode. The observed apparent resistivity and phase TM mode data were interpreted using the 2D inversion algorithm, and a good data fitting could be obtained. The resistivity structure beneath the survey area (down to a depth of 80 m) could be derived and the buried faults could be mapped successfully. In addition to the CSRMT observations, a conventional radiomagnetotelluric (RMT) survey was also carried out on the same profile. An excellent correlation of the observed RMT and CSRMT transfer functions and 2D conductivity models was achieved.

Controlled source radiomagnetotellurics: A tool for near surface investigations in remote regions

In the course of near surface measurements, we investigated the application of the controlled source radiomagnetotelluric (CSRMT) sounding method. The used equipment RMT-C includes a five-channel recorder (frequency range 1–1000kHz), two electric and three magnetic sensors, and a transmitter (1kW) with a horizontal electric dipole as a source. The electric field measurements are carried out with ungrounded lines thus enabling surveys in the wintertime on snow and ice, and in the summertime in conditions adverse to the groundings. The use of electric dipole as a source ensures operating in a wide frequency range at a significant distance (up to 3–4km) from the source. Measurements of basic signals and their odd sub harmonics raise the efficiency of the field surveys. For the covering of the full frequency range 1–1000kHz, three basic frequencies are usually transmitted; each one is accompanied by 8–12 sub harmonics. This provides a high rate of measurements, i.e. about 70 sounding stations per day - 10 times faster than by the vertical electric sounding (VES) method. The CSRMT method offers the possibility to work in the far-field zone of the controlled source where the plane wave approximation of the primary field can be used. For the analysis of the high-frequency field of the electric dipole, an approach is used considering quasi-stationary and wave zones of the source. Features of the electromagnetic field are analyzed by taking into account displacement currents in the air. In the wave zone, at a certain distance from the source (tens to hundreds of meters), field components are characterized by a number of differences compared to a quasi-stationary zone (slower drop of the field's amplitudes, change of the directional diagram of the source, appearance of the ellipticity of polarization of electric and magnetic fields). The CSRMT method is intended for the application in remote regions where surveys by the standard RMT method could only use Very Low Frequency (VLF) radio transmitter signals. Examples for the application of the CSRMT method are presented from the Arkhangelsk region (geophysical support of an onshore gas pipeline construction), and the Chukotka Peninsula (survey at a site of the ore-dressing factory construction). The CSRMT results were confirmed by the subsequent drilling.

Research on technology of Controlled Source Radio MagnetoTelluric system’s transmitter

The basic theory and principle of the artificial Controlled Source Radio MagnetoTelluric (CSRMT) method is studied and a novel CSRMT transmitter in kHz frequency band is designed. The specific circuit and measured results of some key modules in transmitter are presented, and some outdoor experimental tests have been carried out, which shows that the completed prototype of transmitter can generate a continuous sine current with frequency up to 35.33 kHz, peak-to-peak amplitude up to 40 A, and realize a reliable transmitting mode with multi-frequency and high-current. The transmitter has a wide operating band, large magnetic moment and high waveform fidelity, and can meet the requirements of shallow geological exploration with high-resolution.

Integration of controlled-source and radio magnetotellurics, electric resistivity tomography, and reflection seismics to delineate 3D structures of a quick-clay landslide site in southwest of Sweden

Quick clay, which is the main cause of landslides that occur in the northern countries, liquefies easily, and its presence implies an increased risk of landslide. Geophysical methods have been increasingly used in landslide investigations. Three-dimensional electric resistivity tomography, radio magnetotelluric (RMT), controlled-source RMT (CSRMT), and high-resolution reflection seismic data were acquired at a quick-clay landslide site in the southwest of Sweden. The main objectives were to evaluate the capability of each method in delineating different subsurface geologic structures that controlled a peculiar and hazardous retrogressive-type landslide in the study area. A 3D resistivity model from the inversion of CSRMT data showed the best correlation with the reflection seismic data and borehole information, thanks to the broad frequency range of the data set. It better imaged the resistive crystalline bedrock underlying the marine conductive clays and showed considerable correlations with the 3D reflection seismic data in resolving a coarse-grained layer that was interpreted to act as a conduit directing freshwater into the clays under a confined pressure, leaching their salt and forming quick clays. The 3D CSRMT resistivity model and 3D reflection seismic data showed that the coarse-grained layer has a varying thickness. At some locations, it was too thin to be resolved by the methods used here. Combination of the CSRMT model, reflection seismic data, and the borehole data suggested that a layer with thickness of approximately 5 m and resistivity between [Formula: see text] was potentially quick clay, which probably extended laterally in the entire study area. These observations suggested that future developments should focus on joint inversion of such 3D data sets incorporating sharp boundaries as constraints in the inversion and particularly when quick clays were studied.

CSRMT measurements in the frequency range of 1–250kHz to map a normal fault in the Volvi basin, Greece

Abstract In order to gain a better understanding of the geometry of surface faults, five Controlled Source/Radio Magnetotelluric (CSRMT) profiles were measured across the Volvi basin, 45 km northeast of the city of Thessaloniki in Greece. The data were collected in two frequency ranges: a) 1–12.5 kHz using a remotely controlled double horizontal magnetic dipole transmitter (CSAMT measurements), and b) 15–250 kHz using the signal from distant radio transmitters (RMT measurements). The transition from the RMT band to the CSAMT band was smooth and continuous allowing us to combine both datasets for plane-wave modeling. The surface geology shows a predominantly 2D structure, and therefore we planned the survey into profiles perpendicular to the geological strike. We have used a 2D interpretation tool to model the data in TE, TM, TE + TM and determinant modes. Using a 4% error floor on the impedance, 2D resistivity models from inversion of the determinant data provide lower RMS data fits (4.2 and 1.2 for resistivity and phase, respectively) compared to the combined TE + TM data (4.4, 2.8, overall resistivity and phase, respectively). 2D inversion of the measured tensor data shows a sharp change in the depth to the top of resistive gneiss–schist basement that is overlain by a less resistive overburden at southern basin flanks. The change in depth to the bedrock is clearly seen in all 2D models along the measured profiles suggesting the existence of normal faults with strike directions of NE–SW to E–W. The 2D electrical resistivity models suggest that the bedrock deepens towards south-west. The resistivity models are also compared with the existing borehole information in the area and show a reasonable correlation. For example the sharp change of depth to the bedrock towards the center of the basin as seen in the resistivity models are also confirmed by the borehole data.

Investigation of groundwater resources using controlled-source radio magnetotellurics (CSRMT) in glacial deposits in Heby, Sweden

We have combined tensor radio magnetotelluric- (RMT, 15–250 kHz) and controlled source tensor magnetotelluric (CSTMT, 1–12 kHz) data for the mapping of aquifers in gravel formations lying in between crystalline bedrock and clay rich sediments in the Heby area some 40 km west of Uppsala in Sweden. The estimated transfer functions, the impedance tensor and the tipper vector generally satisfy 1D or 2D necessary conditions except for the lowest CSTMT frequencies where near field effects become more dominant. The data measured from 8 profiles were inverted with the Rebocc code of Siripunvaraporn and Egbert (2000) assuming plane wave conditions. This meant that only 12 frequencies in the range of 4–180 kHz could be used. The four lowest frequencies of CSTMT in the range of 1–2.8 kHz were excluded because of source effects. Data from all profiles were inverted with a starting model of 100 Ω-m and a relative error floor of 0.02 on apparent resistivity, corresponding to less than 1° on phase. Tipper vectors are generally small except when source effects become dominant in the lowest frequencies of CSTMT and were therefore not used for inversion. Comparing with models derived from vertical electrical soundings, refraction and reflection seismic data as well as ground truth from exploration wells assessed the reliability of the deep part of the models. Furthermore we carried out a non-linear resolution analysis to better quantify the depth extent of the aquifers. The inverted models from the Heby area show well the thickness variations of glacial deposits overlying crystalline bedrock. Generally, the upper 20 m of the models have resistivities below 40 Ω-m, taken to represent clay rich formations. Below the clay layer resistivities increase to about 40–400 Ω-m, interpreted to represent sand/gravel formations with a maximum thickness of about 40 m and a width of several hundred metres. This is a potential aquifer that extends in approximately N–S direction for some kilometres.