Magnetotellurics (MT) is an electromagneticgeophysical method for inferring the earth's subsurface electrical conductivity from measurements of natural geomagnetic and geoelectric field variation at the Earth's surface.
Investigation depth ranges from 100 m below ground by recording higher frequencies down to 200 km or deeper with long-period soundings. Proposed in Japan in the 1940s, and France and the USSR during the early 1950s, MT is now an international academic discipline and is used in exploration surveys around the world.
Commercial uses include hydrocarbon (oil and gas) exploration, geothermal exploration, carbon sequestration, mining exploration, as well as hydrocarbon and groundwater monitoring. Research applications include experimentation to further develop the MT technique, long-period deep crustal exploration, deep mantle probing, sub-glacial water flow mapping, and earthquake precursor research.
History
The magnetotelluric technique was introduced independently by Japanese scientists in 1948[1] (Hirayama, Rikitake), Soviet geophysicist Andrey Nikolayevich Tikhonov in 1950[2] and the French geophysicist Louis Cagniard in 1953.[3] With advances in instrumentation, processing and modelling, magnetotellurics has become one of the most important tools in deep Earth research.
Since first being created in the 1950s, magnetotelluric sensors, receivers and data processing techniques have followed the general trends in electronics, becoming less expensive and more capable with each generation. Major advances in MT instrumentation and technique include the shift from analog to digital hardware, the advent of remote referencing, GPS time-based synchronization, and 3D data acquisition and processing.
Commercial applications
Hydrocarbon exploration
For hydrocarbon exploration, MT is mainly used as a complement to the primary technique of reflection seismology exploration.[4][5][6][7] While seismic imaging is able to image subsurface structure, it cannot detect the changes in resistivity associated with hydrocarbons and hydrocarbon-bearing formations. MT does detect resistivity variations in subsurface structures, which can differentiate between structures bearing hydrocarbons and those that do not.[8]
At a basic level of interpretation, resistivity is correlated with different rock types. High-velocity layers are typically highly resistive, whereas sediments – porous and permeable – are typically much less resistive. While high-velocity layers are an acoustic barrier and make seismic ineffective, their electrical resistivity means the magnetic signal passes through almost unimpeded. This allows MT to see deep beneath these acoustic barrier layers, complementing the seismic data and assisting interpretation.[9] 3-D MT survey results in Uzbekistan (32 x 32 grid of soundings) have guided further seismic mapping of a large known gas-bearing formation with complex subsurface geology.[10][11]
INCO's 1991 proof-of-concept study in Sudbury, Ontario, Canada sensed a 1750-meter-deep nickel deposit. Falconbridge followed with a feasibility study in 1996 that accurately located two Ni-Cu mineralized zones at about 800 m and 1350 m depth. Since then, both major and junior mining companies are increasingly using MT and audio-magnetotellurics (AMT) for both brownfields (near known deposits) and greenfields (uncharted land) exploration. Significant MT mapping work has been done on areas of the Canadian Shield.[13]
Diamond exploration, by detecting kimberlites, is also a proven application.[14]
Geothermal exploration
MT geothermal exploration measurements allow detection of resistivity anomalies associated with productive geothermal structures, including faults and the presence of a cap rock, and allow for estimation of geothermal reservoir temperatures at various depths.[15][16][17] Dozens of MT geothermal exploration surveys have been completed in Japan and the Philippines since the early 1980s, helping to identify several hundred megawatts of renewable power at places such as the Hatchobaru plant on Kyushu[18][19] and the Togonang plant on Leyte.[20][21][22] Geothermal exploration with MT has also been done extensively in the United States, Iceland,[23] New Zealand, Hungary,[16] China,[24]Ethiopia, Indonesia, Peru,[25] Australia, and India.[26]
Other
MT is also used for groundwater exploration and mapping research, hydrocarbon reservoir monitoring, deep investigation (100 km) of the electrical properties of the bedrock for high-voltage direct current (HVDC) transmission systems,[27]carbon dioxide sequestration,[28][29] and other environmental engineering applications (e.g. nuclear blast site monitoring[30] and nuclear waste disposal site monitoring).
Research applications
Crust and mantle
Since the MT is highly sensitive to the composition and temperature of the Earth, it has been widely used to understand numerous geological phenomena in the Earth's mantle and crust. These include investigating the composition and distribution of melts,[31][32] understanding fault mechanics and earthquake generation,[33] imaging deep lithospheric architecture and composition, which can be tied to many geodynamic processes.[34][35] Large investigations have focused on the conterminous US (e.g. the National Science Foundation EarthScope MT Program and its successor NASA and USGS MTArray[36]), the East Pacific Rise, Australia (AusLAMP MT Program[37]), Southern Africa (SAMTEX MT Project[38]), China (Part of the Sinoprobe project[39]) and the Tibetan Plateau.
Earthquake precursor prediction
Fluctuations in the MT signal may be able to predict the onset of seismic events.[40][41][42] Stationary MT monitoring systems have been installed in Japan since April 1996, providing a continuous recording of MT signals at the Wakuya Station (previously at the Mizusawa Geodetic Observatory) and the Esashi Station of the Geographical Survey Institute of Japan (GSIJ). These stations measure fluctuations in the Earth's electromagnetic field that correspond with seismic activity.[43] The raw geophysical time-series data from these monitoring stations is freely available to the scientific community, enabling further study of the interaction between electromagnetic events and earthquake activity.[44]
Additional MT earthquake precursor monitoring stations in Japan are located in Kagoshima, in Sawauchi, and on Shikoku. Similar stations are also deployed in Taiwan on Penghu Island, as well as in the Fushan Reserve on the island of Taiwan proper.[45]
POLARIS is a Canadian research program investigating the structure and dynamics of the Earth's lithosphere and the prediction of earthquake ground motion.[46]
Different rocks, sediments and geological structures have a wide range of different electrical conductivities. Measuring electrical resistivity allows different materials and structures to be distinguished from one another and can improve knowledge of tectonic processes and geologic structures.
The Earth's naturally varying electric and magnetic fields are measured over a wide range of magnetotelluric frequencies from 10,000 Hz to 0.0001 Hz (10,000 s). These fields are due to electric currents flowing in the Earth and the magnetic fields that induce these currents. The magnetic fields are produced mainly by the interaction between the solar wind and the magnetosphere. In addition, worldwide thunderstorm activity causes magnetic fields at frequencies above 1 Hz. Combined, these natural phenomena create strong MT source signals over the entire frequency spectrum.
The ratio of the electric field to magnetic field provides simple information about subsurface conductivity. Because the skin effect phenomenon affects the electromagnetic fields, the ratio at higher frequency ranges gives information on the shallow Earth, whereas deeper information is provided by the low-frequency range. The ratio is usually represented as both apparent resistivity as a function of frequency and phase as a function of frequency.
MT measurements can investigate depths from about 300 m down to hundreds of kilometers, though investigations in the range of 500 m to 10,000 m are typical. Greater depth requires measuring lower frequencies, which in turn requires longer recording times. Very deep, very long-period measurements (mid-crust through upper mantle depths), may require recordings of several days to weeks or more to obtain satisfactory data quality.
Horizontal resolution of MT mainly depends on the distance between sounding locations- closer sounding locations increase the horizontal resolution. Continuous profiling (known as Emap) has been used, with only meters between the edges of each telluric dipole.
Vertical resolution of MT mainly depends on the frequency being measured, as lower frequencies have greater depths of penetration. Accordingly, vertical resolution decreases as depth of investigation increases.
Signal strength and recording times
Magnetic fields in the frequency range of 1 Hz to approximately 20 kHz are part of the audio-magnetotelluric (AMT) range. These are parallel to the Earth surface and move towards the Earth's centre. This large frequency band allows for a range of depth penetration from several metres to several kilometres below the Earth's surface. Due to the nature of magnetotelluric source, the waves generally fluctuate in amplitude height. Long recording times are needed to ascertain usable reading due to the fluctuations and the low signal strength. Generally, the signal is weak between 1 and 5 kHz, which is a crucial range in detecting the top 100 m of geology. The magnetotelluric method is also used in marine environments for hydrocarbon exploration and lithospheric studies.[49] Due to the screening effect of the electrically conductive sea water, a usable upper limit of the spectrum is around 1 Hz.
2D and 3D magnetotellurics
Two-dimensional surveys consist of a longitudinal profile of MT soundings over the area of interest, providing two-dimensional "slices" of subsurface resistivity.
Three-dimensional surveys consist of a loose grid pattern of MT soundings over the area of interest, providing a more sophisticated three-dimensional model of subsurface resistivity.
Variants
Audio-magnetotellurics
Audio-magnetotellurics (AMT) is a higher-frequency magnetotelluric technique for shallower investigations. While AMT has less depth penetration than MT, AMT measurements often take only about one hour to perform (but deep AMT measurements during low-signal strength periods may take up to 24 hours) and use smaller and lighter magnetic sensors. Transient AMT is an AMT variant that records only temporarily during periods of more intense natural signal (transient impulses), improving signal-to-noise-ratio at the expense of strong linear polarization.[50]
Controlled source electromagnetics
CSEM controlled source electro-magnetic is a deep-water offshore variant of controlled source audio magnetotellurics; CSEM is the name used in the offshore oil and gas industry.[51] and for onshore exploration mostly Lotem is used in Russia, China the USA and Europe[52][53]
Onshore CSEM / CSAMT may be effective where electromagnetic cultural noise (e.g. power lines, electric fences) present interference problems for natural-source geophysical methods. An extensive grounded wire (2 km or more) has currents at a range of frequencies (0.1 Hz to 100 kHz) passed through it. The electric field parallel to the source and the magnetic field which is at right angles are measured. The resistivity is then calculated, and the lower the resistivity, the more likely there is a conductive target (graphite, nickel ore or iron ore). CSAMT is also known in the oil and gas industry as onshore controlled source electromagnetics (Onshore CSEM).
An offshore variant of MT, the marine magnetotelluric (MMT) method,[54][page needed] uses instruments and sensors in pressure housings deployed by ship into shallow coastal areas where water is less than 300 m deep.[6][55][56][57][58] A derivative of MMT is offshore single-channel measurement of the vertical magnetic field only (the Hz, or "tipper"), which eliminates the need for telluric measurements and horizontal magnetic measurements.[59]
Exploration surveys
MT exploration surveys are done to acquire resistivity data which can be interpreted to create a model of the subsurface. Data is acquired at each sounding location for a period of time (overnight soundings are common), with physical spacing between soundings dependent on the target size and geometry, local terrain constraints and financial cost. Reconnaissance surveys can have spacings of several kilometres, while more detailed work can have 200 m spacings, or even adjacent soundings (dipole-to-dipole).
The HSE impact of MT exploration is relatively low because of light-weight equipment, natural signal sources, and reduced hazards compared to other types of exploration (e.g. no drills, no explosives, and no high currents).
Remote reference soundings
Remote Reference is an MT technique used to account for cultural electrical noise by acquiring simultaneous data at more than one MT station. This greatly improves data quality, and may allow acquisition in areas where the natural MT signal is difficult to detect because of man-made EM interference.
Equipment
A typical full suite of MT equipment (for a "five component" sounding) consists of a receiver instrument with five sensors: three magnetic sensors (typically induction coil sensors), and two telluric (electric) sensors. For exclusively long-period MT (frequencies below approximately 0.1 Hz), the three discrete broadband magnetic field sensors may be substituted by a single compact triaxial fluxgate magnetometer. In many situations, only the telluric sensors will be used, and magnetic data borrowed from other nearby soundings to reduce acquisition costs.
A complete five-component set of MT equipment can be backpack-carried by a small field team (2 to 4 persons) or carried by a light helicopter, allowing deployment in remote and rugged areas. Most MT equipment is capable of reliable operation over a wide range of environmental conditions, with ratings of typically −25 °C to +55 °C, from dry desert to high-humidity (condensing) and temporary full immersion.
Data processing and interpretation
Post-acquisition processing is required to transform raw time-series data into frequency-based inversions. The resulting output of the processing program is used as the input for subsequent interpretation. Processing may include the use of remote reference data or local data only.
Processed MT data is modelled using various techniques to create a subsurface resistivity map, with lower frequencies generally corresponding to greater depth below ground. Anomalies such as faults, hydrocarbons, and conductive mineralization appear as areas of higher or lower resistivity from surrounding structures. Various software packages are used for interpretation (inversion) of magnetotelluric data, where apparent resistivity is used to create a model of the subsurface.
Instrument and sensor manufacturers
Four companies supply most of the commercial-use world market: one in the United States (Zonge International, Inc.[60]), one in Canada; (Phoenix Geophysics, Ltd.[61]); one in Germany (Metronix Messgeraete und Elektronik GmbH).[62]) and One in Russia (Vega Geophysics, LLC).[63]
^Rikitake, T. (1948). "Notes on electromagnetic induction within the Earth". Bull. Earthq. Res. Inst. 24 (1): 4.
^Tikhonov, A.N. (1950). "On determining electrical characteristics of the deep layers of the Earth's crust". Doklady. 73 (2): 295–297. CiteSeerX10.1.1.462.5940. NAID10004593077.
^Sinharay, Rajib K; Bhattacharya, Bimalendu B. (2001). "An analysis of magnetotelluric (MT) data over geothermal region of Bakreshwar, West Bengal". Journal of Geophysics. 22 (1). Hyderabad: 31–39. INIST1145977.
^Thunehed, H.; Astrom, U.; Westman, B. (2007). "Geophysical and geological pre-investigations for HVDC-electrodes". 2007 IEEE Power Engineering Society Conference and Exposition in Africa - Power Africa. pp. 1–3. doi:10.1109/PESAFR.2007.4498123. ISBN978-1-4244-1477-2. S2CID7541303.
^Unsworth, Martyn; Soyer, Wolfgang; Tuncer, Volkan; Wagner, Anna; Barnes, David (1 May 2007). "Hydrogeologic assessment of the Amchitka Island nuclear test site (Alaska) with magnetotellurics". Geophysics. 72 (3): B47–B57. Bibcode:2007Geop...72B..47U. doi:10.1190/1.2539353. S2CID2288777. INIST18787533.
^Uyeda, Seiya; Nagao, Toshiyasu; Kamogawa, Masashi (2009). "Short-term earthquake prediction: Current status of seismo-electromagnetics". Tectonophysics. 470 (3–4): 205–213. Bibcode:2009Tectp.470..205U. doi:10.1016/j.tecto.2008.07.019.
^Tzanis, A.; Vallianatos, F. (April 2001). "A critical review of Electric Earthquake Precursors". Annali di Geofisica. 44 (2): 429–460. hdl:2122/1213. OSTI20222530.
^Cantwell, T. (1960) Detection and Analysis of Low-Frequency Magnetotelluric Signals, PhD Thesis, Massachusetts Institute of Technology, Cambridge, Massachusetts
^Strack, Kurt (1992). Exploration with Deep Transient Elctromagnetics. Elsevier. ISBN0444895418.
^Strack, K. -M.; Pandey, P. B. (2007). "Exploration with controlled-source electromagnetics under basalt cover in India". The Leading Edge. 26 (3): 360. Bibcode:2007LeaEd..26..360S. doi:10.1190/1.2715056.
^Stéphane, Sainson (2017). Electromagnetic seabed logging : a new tool for geoscientists. Springer. ISBN978-3-319-45355-2.