Three‐dimensional modelling of magnetotelluric data to image Sehqanat hydrocarbon reservoir in southwestern Iran

A detailed magnetotelluric survey was conducted in 2013 in the Sehqanat oil field, southwestern Iran to map the geoelectrical structures of the sedimentary Zagros zone, particularly the boundary between the Gachsaran Formation acting as cap rock and the Asmari Formation as the reservoir. According to the electrical well logs, a large resistivity contrast exists between the two formations. The Gachsaran Formation is formed by tens to hundreds of metres of evaporites and it is highly conductive (ca. 1 Ωm–10 Ωm), and the Asmari Formation consists of dense carbonates, which are considerably more resistive (more than 100 Ωm). Broadband magnetotelluric data were collected along five southwest–northeast directed parallel lines with more than 600 stations crossing the main geological trend. Although dimensionality and strike analysis of the magnetotelluric transfer functions showed that overall they satisfied local 2D conditions, there were also strong 3D conditions found in some of the sites. Therefore, in order to obtain a more reliable image of the resistivity distribution in the Sehqanat oil field, in addition to standard 2D inversion, we investigated to what extent 3D inversion of the data was feasible and what improvements in the resistivity image could be obtained. The 2D inversion models using the determinant average of the impedance tensor depict the main resistivity structures well, whereas the estimated 3D model shows significantly more details although problems were encountered in fitting the data with the latter. Both approaches resolved the Gachsaran–Asmari transition from high conductivity to moderate conductivity. The well‐known Sehqanat anticline could also be delineated throughout the 2D and 3D resistivity models as a resistive dome‐shaped body in the middle parts of the magnetotelluric profiles.     

冻土区天然气水合物各向异性储层的AMT响应特征

陆域天然气水合物通常发育于冻土层下方破碎带和岩层裂隙处,其储层会表现明显电性各向异性特征.音频大地电磁法（AMT）能有效探测陆域水合物储藏范围,且不易受高阻冻土层的压制和干扰,可用于陆域天然气水合物探测研究.本文采用AMT对角各向异性二维正演方法,对多种陆域天然气水合物各向异性储层模型及相关参数进行了模拟试算,分析其视电阻率和相位响应特征.结果表明AMT能清晰显示各向异性储层空间位置和分布情况,并对储层的水合物饱和度变化以及储层数量等特征都有所反映.天然气水合物电性各向异性模拟研究为这一新型潜力能源的勘查工作提供了新的思路.     

基于三维正演的音频大地电磁阻抗相位不变量校正技术

复杂地形、地质条件的大地电磁数据解释容易出现假象,采用三维正演技术模拟地形和地表不均匀体的背景响应,对实测数据阻抗相位不变量进行校正,实现更准确的定性分析;对三维异常体模型的合成数据进行一维、二维多参数反演试算,以确定地形剧变区选择反演技术的最佳方案.合成数据的试反演结果显示一维反演水平切片假异常较多,二维反演能压制测向假异常,但不能压制走向的假异常,水平切片多出现测向条带.使用本文提出的阻抗相位不变量校正法扣除地形、地表背景响应,结合一维、二维反演,能使实际资料解释成果更加可靠.     

The potential of audiomagnetotellurics in the study of geothermal fields: A case study from the northern segment of the La Candelaria Range,northwestern Argentina

Despite its reduced penetration depth, audiomagnetotelluric (AMT) studies can be used to determine a broad range of features related to little studied geothermal fields. This technique requires a stepwise interpretation of results taking into consideration diverse information (e.g. topographic, hydrological, geological and/or structural data) to constrain the characteristics of the study area. In this work, an AMT study was performed at the hot springs in the northern segment of the La Candelaria Range in order to characterize the area at depth. Geometric aspects of the shallow subsurface were determined based on the dimensional and distortion analysis of the impedance tensors. Also, the correlation between structural features and regional strikes allowed us to define two geoelectric domains, useful to determine the controls on fluid circulation. The subsurface resistivity distribution was determined through 1D and 2D models. The patterns of the 1D models were compared with the morpho-structure of the range. Shallow and deep conductive zones were defined and a possible shallow geothermal system scheme proposed. A strong correlation was found between the AMT results and the geological framework of the region, showing the relevance of using AMT in geothermal areas during the early stages of subsurface prospecting.     

Imaging a 3D geological structure from HEM, airborne magnetic and ground ERT data in Kalat-e-Reshm area, Iran

A set of geophysical data collected in an area in Iran are analyzed to check the validity of a geological map that was prepared in connection to a mineral prospecting project and also to image the spatial electrical resistivity distribution. The data set includes helicopter electromagnetic (HEM), airborne magnetic and ground electrical resistivity measurement. Occam approach was used to invert the HEM data to model the resistivity using a layered earth model with fixed thicknesses. The algorithm is based on a nonlinear inverse problem in a least-squares sense.The algorithm was tested on a part of an HEM dataset acquired with a DIGHEM helicopter EM system at Kalat-e-Reshm, Semnan in Iran. The area contains a resistive porphyry andesite that is covered by Eocene sedimentary units. The results are shown as resistivity sections and maps confirming the existence of an arc like resistive structure in the survey area. The resistive andesite seems to be thicker than it is indicated in the geological maps. The results are compared with the reduced to the pole (RTP) airborne magnetic anomaly field data as well as with two ground resistivity profiles. We found reasonable correlations between the HEM 1D resistivity models and 2D models from electrical resistivity tomography (ERT) inversions. A 3D visualization of the 1D models along all flight lines provided a useful tool for the study of spatial variations of the resistivity structure in the investigation area.     

The main structure and tectonic features of the Chernorud-Mukhor site on the western shore of Lake Baikal from TEM and SP measurements

Transient electromagnetic (TEM), self-potential (SP) and geoelectrical mapping measurements were carried out at the Chernorud-Mukhor site in the Priolkhonje area on the western shore of the Lake Baikal. All measurements were made along several profiles across the main strike of the regional Primorsky fault. TEM measurements were carried out in a time range from a few tens of microseconds to several tens of milliseconds. The most important result of the 1D modelling of TEM soundings is the discovery of nearly horizontal boundaries that divide high resistive overlying and well conducting underlying rocks. The resistivity of the former is in the range from 100 Ωm to 1000 Ωm, while the resistivity of the latter varies from less than 1 Ωm to several tens of Ωm. This good conductive zone could also be verified by geoelectrical mapping using Schlumberger array (AB/2=100 m). Due to high conductivity of the underlying rocks only the upper boundary of the conductive layer could be determined by TEM soundings. A regional SP anomaly with amplitude of about −450 mV has also been observed above the low resistivity zone indicating the electron nature of its conductance. Geologically, the conductive zone is represented by a graphite-bearing layer within the region of archean rocks. Since that layer extends over a large area, it may be used as a key in studying structures and tectonics of the Priolkhonje area. A 1D TEM geoelectric section shows a wide, gently sloping syncline as a probable base structure of the Chernorud-Mukhor site. Neotectonic faults divide the syncline into vertically displaced blocks that form a wide complicated graben with a total amplitude of about 250 m.     

Geoelectrical structure and hydrogeological investigations of the southern Rif Cordillera (Morocco)

An electrical prospecting survey is conducted in the Rharb basin, a semi‐arid region in the southern part of the Rifean Cordillera (Morocco) to delineate characteristics of the aquifer and the groundwater affected by the marine intrusion related to Atlantic Ocean. Analysis and interpretations of electrical soundings, bi‐logarithmic diagrams and the geoelectrical sections highlight a monolayer aquifer in the southern part, a multilayer system in the northern part of the Rharb basin and lenticular semi‐permeable formations. Several electrical layers have been deduced from the analysis of bi‐logarithmic diagrams: resistant superficial level (R₀), conducting superficial level (C₀), resistant level (R), intermediary resistant level (R′), conducting level (Cp) and intermediary layer of resistivity (AT). Spatial distribution of the resistivity deduced from the interpretation of apparent resistivity maps (AB = 400 and 1000 m) and the decreasing of resistivity values (35–10 Ωm), in particular in the coastal zone show that this heterogeneity is related to several anomalies identified in the coastal area, which result from hydraulic and geological processes: (i) heterogeneous hydraulic conductivity in particular in the southern part of the Rharb; (ii) lateral facies and synsedimentary faulting and (iii) the relationship between the electrical conductivity and chloride concentration of groundwater shows that salinity is the most important factor controlling resistivity. The distribution of fresh/salt‐water zones and their variations in space along geoelectrical sections are established through converting subsurface depth‐resistivity models. Copyright © 2010 John Wiley & Sons, Ltd.     

New insights into the hydrogeology of a basaltic shield volcano from a comparison between self-potential and electromagnetic data: Piton de la Fournaise, Indian ocean

In order to investigate aquifers, several geophysical surveys have been carried out in the Baril area of the southern flank of Piton de la Fournaise volcano on Reunion in the Indian Ocean using audiomagnetotelluric (AMT), very-low-frequency (VLF) and self-potential (SP) methods. We present the results with emphasis on a comparison between SP data and the findings of geoelectric surveys. AMT soundings have indicated, from the surface downward, three layers: (i) resistive volcanic rocks, (ii) an intermediate resistivity layer, and (iii) a conductive basement attributed to a seawater-bearing aquifer. VLF measurements allow the mapping of the first layer apparent resistivity, and therefore its bottom, when the true resistivity is supposed to be isotropic and homogenous. When this assumption does not hold, only the SP method permits the mapping of this bottom. Because of the good agreement between the SP and electromagnetic results, we propose the SP method as the first tool that should be used in studying shallow hydrogeological structures in volcanic areas.     

Integrated geophysical and chemical study of saline water intrusion

Surface geophysical surveys provide an effective way to image the subsurface and the ground water zone without a large number of observation wells. DC resistivity sounding generally identifies the subsurface formations-the aquifer zone as well as the formations saturated with saline/brackish water. However, the method has serious ambiguities in distinguishing the geological formations of similar resistivities such as saline sand and saline clay, or water quality such as fresh or saline, in a low resistivity formation. In order to minimize the ambiguity and ascertain the efficacy of data integration techniques in ground water and saline contamination studies, a combined geophysical survey and periodic chemical analysis of ground water were carried out employing DC resistivity profiling, resistivity sounding, and shallow seismic refraction methods. By constraining resistivity interpretation with inputs from seismic refraction and chemical analysis, the data integration study proved to be a powerful method for identification of the subsurface formations, ground water zones, the subsurface saline/brackish water zones, and the probable mode and cause of saline water intrusion in an inland aquifer. A case study presented here illustrates these principles. Resistivity sounding alone had earlier failed to identify the different formations in the saline environment. Data integration and resistivity interpretation constrained by water quality analysis led to a new concept of minimum resistivity for ground water-bearing zones, which is the optimum value of resistivity of a subsurface formation in an area below which ground water contained in it is saline/brackish and unsuitable for drinking.     

金属护栏对洛阳地震台地电阻率干扰定量分析

自2018年12月起,因洛阳地震台测区内安装金属护栏,该台地电阻率出现异常变化。为正确评估此次影响,结合洛阳地震台电测深曲线、岩性资料建立了三维水平层状结构模型,利用有限元法计算了测区内金属护栏对洛阳地震台地电阻率观测的影响。同时,建立典型“H”型电性结构模型,分析了不同属性的金属导线对地电阻率观测的影响。结果显示：①金属护栏是洛阳地震台此次地电阻率年变异常的主要影响因素。当位于测区内影响系数为负的区域时,金属护栏引起地电阻率的趋势性上升;当位于影响系数为正的区域时,金属护栏引起地电阻率的趋势性下降。②不同属性的金属导线对地电阻率观测的影响不同。导线与测道夹角减小、导线长度增加、导线电阻率减小以及导线横截面积增加等均会导致干扰源对地电阻率年变影响作用增强。     

3D electrical resistivity inversion using prior spatial shape constraints

To minimize the number of solutions in 3D resistivity inversion, an inherent problem in inversion, the amount of data considered have to be large and prior constraints need to be applied. Geological and geophysical data regarding the extent of a geological anomaly are important prior information. We propose the use of shape constraints in 3D electrical resistivity inversion, Three weighted orthogonal vectors （a normal and two tangent vectors） were used to control the resistivity differences at the boundaries of the anomaly. The spatial shape of the anomaly and the constraints on the boundaries of the anomaly are thus established. We incorporated the spatial shape constraints in the objective function of the 3D resistivity inversion and constructed the 3D resistivity inversion equation with spatial shape constraints. Subsequently, we used numerical modeling based on prior spatial shape data to constrain the direction vectors and weights of the 3D resistivity inversion. We established a reasonable range between the direction vectors and weights, and verified the feasibility and effectiveness of using spatial shape prior constraints in reducing excessive structures and the number of solutions. We applied the prior spatially shape-constrained inversion method to locate the aquifer at the Guangzhou subway. The spatial shape constraints were taken from ground penetrating radar data. The inversion results for the location and shape of the aquifer agree well with drilling data, and the number of inversion solutions is significantly reduced.     

Mapping Coastal Aquifers by Joint Inversion of DC and TEM Soundings-Three Case Histories

Electrical and electromagnetic methods are well suited for coastal aquifer studies because of the large contrast in resistivity between fresh water-bearing and salt water-bearing formations. Interpretation models for these aquifers typically contain four layers: a highly resistive unsaturated zone; a surficial fresh water aquifer of intermediate resistivity; an underlying conductive, salt water saturated aquifer; and resistive substratum. Additional layers may be added to allow for variations in lithology within the fresh water and salt water layers. Two methods are evaluated: direct current resistivity and time domain electromagnetic soundings. Use of each method alone produces nonunique solutions for resistivities and/or thicknesses of the different layers. We show that joint inversion of vertical electric and time domain electromagnetic soundings produces a more tightly constrained interpretation model at three test sites than is produced by inversion methods applied to each data set independently.     

3D electrical structure of porphyry copper deposit: A case study of Shaxi copper deposit

Located in Lu-Zong ore concentration area, middle-lower Yangtze metallogenic belt, ShaXi porphyry copper deposit is a typical hydrothermal deposit. To investigate the distribution of deep ore bodies and spatial characteristics of host structures, an AMT survey was conducted in mining area. Eighteen pseudo-2D resistivity sections were constructed through careful processing and inversion. These sections clearly show resistivity difference between the Silurian sandstones formation and quartz diorite porphyry and this porphyry copper formation was controlled by the highly resistive anticlines. Using 3D block Kriging interpolation method and 3D visualization techniques, we constructed a detailed 3D resistivity model of quartz diorite porphyry which shows the shape and spatial distribution of deep ore bodies. This case study can serve as a good example for future ore prospecting in and around this mining area.     

Geoelectrical and geological structure of the crust in Western Slovakia

Electrical resistivity of the Earth’s crust is sensitive to a wide range of petrological and physical parameters, and it particularly clearly indicates crustal zones that have been tectonically or thermodynamically disturbed. A complex geological structure of the Alpine nappe system, remnants of older Hercynian units and Neogene block tectonics in Western Slovakia has been a target of recent magnetotelluric investigations which made a new and more precise identification of the crustal structural elements of the Western Carpathians possible. A NW-SE magnetotelluric profile, 150 km long, with 30 broad-band and 3 long-period magnetotelluric sites, was deployed, crossing the major regional tectonic elements listed from the north: Brunia (as a part of the European platform), Outer Carpathian Flysch, Klippen Belt, blocks of Penninic or Oravicum crust, Tatricum and Veporicum. Magnetotelluric models were combined with previous seismic and gravimetric results and jointly interpreted in the final integrated geological model. The magnetotelluric models of geoelectrical structures exhibit strong correlation with the geological structures of the crust in this part of the Western Carpathians. The significant resemblance in geoelectrical and crustal geological structures are highlighted in shallow resistive structures of the covering formations represented by mainly Tertiary sediments and volcanics. Also in the deeper parts of the crust highly resistive and conductive structures are shown, which reflect the original building Hercynian crust, with superposition of granitoids or granitised complexes and lower metamorphosed complexes. Another important typical feature in the construction of the Western Carpathians is the existence of young Neogene steep fault zones exhibited by conductive zones within the whole crust. The most significant fault zones separate individual blocks of the Western Carpathians and the Western Carpathians itself from the European Platform.     

Audio frequency magneto-telluric survey of Norikura Volcano in central Japan

Norikura Volcano has not been active during the last 10,000 years in spite of the activity of the surrounding volcanic mountains. To study past volcanic activities, geological studies were carried out extensively. However, quite a few geophysical investigations were conducted to contribute to volcanology. Our objective is to detect the present subsurface structure of Norikura Volcano and to define volcanic stratifications. In the vicinity of Norikura Volcano, geothermal fields are still active. Subsurface volcanic rocks in this area have been exposed to geothermal activity and altered. To comprehend volcanic stratifications of Norikura and geothermal activity, we conducted audio frequency magneto-telluric (AMT) surveys around Norikura Volcano. AMT survey is useful in clearly defining the resistivity structure related to volcanic regions. The AMT data were acquired over a frequency range 10 Hz–10 kHz. Decomposition analysis was applied to the tensor impedance data. Subsequently, apparent resistivity and phase data were inverted using a two-dimensional magneto-telluric (MT) inversion and a model of Norikura was derived. The final model manifests that the surface resistors are in agreement with andesite lava or dacite lava. As for the deeper structure, a horizontal conductor is situated above resistive basements. The alteration of the conductor was weak, while basement rocks were strongly altered and/or heated through the thermal activity. The existence of these layers seems to indicate the degree of thermal activity of Norikura Volcano.     

Study of the Santa Catarina aquifer system (Mexico Basin) using magnetotelluric soundings

A tensor magnetotelluric test survey was carried out in the region of Santa Catarina, located in the Chalco sub-basin of the Mexico Basin. The objective was to define the stratification at depth with an emphasis on the geometry of the main aquifer of that region which is partially known from DC resistivity soundings and drilling. High-quality magnetotelluric soundings could be recorded in the immediate vicinity of large urban zones because the sub-surface is very conductive. Interpretation shows that the solid bedrock is located at a depth of at least 800 m to the south and 1300 m to the north; it could, however, be much deeper. Using complementary DC resistivity sounding and well-logging data, three main layers have been defined overlying the bedrock. These layers are, from surface to bottom, an unsaturated zone of sand, volcanic ash and clay about 10 m thick, followed by a very conductive (1.5 ohm·m) 200 m thick layer of sand and ash with intercalated clay, saturated with highly mineralized water, and finally a zone with resistivity increasing gradually to 60 ohm·m. The investigated deep aquifer constitutes most of this third layer. It consists of a sequence of sand, gravel, pyroclastites and mainly fractured basalts. MT resistivity soundings and magnetic transfer functions also indicate that a shallow resistive structure is dipping, from the northwest, into the lacustrine deposits of the basin. This geologic feature is likely to be highly permeable fractured basaltic flows, which provide a channel by which water contaminated by the Santa Catarina landfill may leak into the basin.     

Multiple caldera collapses inferred from the shallow electrical resistivity signature of the Las Cañadas caldera, Tenerife, Canary Islands

The Las Cañadas caldera of Tenerife (LCC) is a well exposed caldera depression filled with pyroclastic deposits and lava flows from the active Teide–Pico Viejo complex (TPVC). The caldera's origin is controversial as both the formation by huge lateral flank collapse(s) and multiple vertical collapses have been proposed. Although vertical collapses may have facilitated lateral slope failures and thus jointly contribute to the exposed morphology, their joint contribution has not been clearly demonstrated. Using results from 185 audiomagnetotelluric (AMT) soundings carried out between 2004 and 2006 inside the LCC, our study provides consistent geophysical constraints in favour of multiple vertical caldera collapse. One-dimensional modelling reveals a conductive layer at shallow depth (30–1000 m), presumably resulting from hydrothermal alteration and weathering, underlying the infilling resistive top layer. We present the resistivity distribution of both layers (resistivity images), the topography of the conductive layer across the LCC, as well as a cross-section in order to highlight the caldera's evolution, including the distribution of earlier volcanic edifices. The AMT phase anisotropy reveals the structural and radial characteristics of the LCC.     

On deep electric conductivity of the lithosphere in the Tobol-Ishim interfluve (West Siberia)

On the basis of exploration magnetotelluric soundings (MTS) carried out by the Tyumen Geological Survey in 1980–1981 within the Tobol-Ishim interfluve, the top of the conductive horizon with the electric resistivity of dozens of Ω m is identified in the geoelectrical section at a depth of 75–80 km. Two segments of the locally elevated top of this horizon are revealed. One segment, where the depth of the top of the conductor is 60–65 km, corresponds to the region of the Kiselevskii fault-the main fault in this territory; another segment is associated with the edge zone of the Ishim branch of the Triassic rift system within the West Siberian Plate. Within this segment, the top of the conductive layer rises to a depth of 55-60 km. It is assumed that the Ishim geoelectrical anomaly, which projects into the geothermal anomaly and is coherent with the features of the deep geological structure of its display region, is conditioned by the element of the mantle-crustal magmatogene fluid paleosystem.     

Integrated interpretation of the magnetotelluric and magnetic data from Mahallat geothermal field,Iran

Magnetotelluric (MT) and ground magnetic surveys were conducted on the Mahallat geothermal field situated in Markazi province, central Iran, as a primary part of the explorations and developments of a geothermal energy investigation program in the region. Mahallat region has the greatest geothermal fields in Iran. MT survey was performed in November 2011 on an 8 km profile crossing the hot springs with a total of 17 stations. The 2D inversion of the determinant MT data was performed using a 2D inversion routine based on the Occam approach. The 2D resistivity model obtained from the determinant data shows a low resistivity zone at 800-2000 m depth and a higher resistivity zone above the low resistivity zone, interpreted as geothermal reservoir and cap rock, respectively. It also revealed two major concealed faults which are acting as preferential paths for the circulation of hydrothermal fluids. To obtain more geophysical evidence, a ground magnetic survey with 5000 stations was also performed over an area of 200 km² around the MT profile. Magnetic measurements show a main positive anomaly of about +1000 nT over the study area, which could be interpreted as an intrusive body with the high magnetic susceptibility (i.e. mafic and ultramafic rocks) into the sedimentary host rocks. We interpret the body as the heat source of the geothermal system. Structural index and depth estimation of the anomaly indicate that the intrusive body is similar to a cylinder extending from about one kilometer depth down to greater depths. The results of MT and magnetic investigations indicate a geothermal reservoir which proves the preliminary geological observations to a great extent.     

Sensitivity of transversal apparent resistivity to the changes in electrical resistivity of the elements of a 2D geoelectrical structure

The sensitivity of transversal apparent resistivity to the changes in the electrical resistivity of elements of a two-dimensional (2D) geoelectrical structure is studied by numerical modeling. This sensitivity is found to have a series of specific features, due to which the monitoring of transversal apparent resistivity can serve as a helpful addition to the monitoring of longitudinal apparent resistivity to trace the dynamics of subsurface and shallow crustal elements of the medium. It is shown that the method previously suggested for the inversion of relative changes in apparent resistivity into relative changes in electrical resistivities of the elements of the geoelectrical structure is applicable to the transversal electrical resistivity.