Depth of crustal conducting layer and asthenosphere in the Pannonian basin determined by magnetotellurics

On the basis of numerical modelling and statistical investigations of the horizontal electrical inhomogeneities (narrow fractures, change of conductance of the surface sediments) ρ_min magnetotelluric (MT)-sounding curves were found to be less distorted and, therefore, they can be used for the determination of deep conducting formations in the crust and upper mantle.Carefully selected MT curves indicate the crustal conductive layer at depths of 7–13 km and the conductive asthenosphere between 50 and 80 km.The relatively shallow position of the asthenosphere corresponds to the high heat flow (q > 80 mWm^?1) in the Pannonian basin according to the empirical relation derived by Ádám (1978).     

The distribution of deep source rocks in the GS Sag: joint MT–gravity modeling and constrained inversion

The coal-bearing strata of the deep Upper Paleozoic in the GS Sag have high hydrocarbon potential. Because of the absence of seismic data, we use electromagnetic (MT) and gravity data jointly to delineate the distribution of deep targets based on well logging and geological data. First, a preliminary geological model is established by using three-dimensional (3D) MT inversion results. Second, using the formation density and gravity anomalies, the preliminary geological model is modified by interactive inversion of the gravity data. Then, we conduct MT-constrained inversion based on the modified model to obtain an optimal geological model until the deviations at all stations are minimized. Finally, the geological model and a seismic profile in the middle of the sag is analysed. We determine that the deep reflections of the seismic profile correspond to the Upper Paleozoic that reaches thickness up to 800 m. The processing of field data suggests that the joint MT–gravity modeling and constrained inversion can reduce the multiple solutions for single geophysical data and thus improve the recognition of deep formations. The MT-constrained inversion is consistent with the geological features in the seismic section. This suggests that the joint MT and gravity modeling and constrained inversion can be used to delineate deep targets in similar basins.     

Comparative analysis of hypotheses of the geoelectric structure of the Transcaucasian region from magnetotelluric data

Various hypotheses of the geoelectric structure of the Transcaucasian region reflecting a priori geological and geophysical information and data from profile magnetotelluric (MT) soundings are analyzed. These hypotheses are used to construct simplified 3-D models of electrical conductivity differing in the patterns of the shallow and deep structure of the region. Numerical modeling of MT fields is performed. Comparative analysis of its results indicates that the most probable conductivity model consistent with the available data is a model involving a high-conductivity channel connecting the Black and Caspian seas.     

A Coastal Aquifer Study Using Magnetotelluric and Gravity Methods in Abo Zenema, Egypt

Magnetotelluric (MT) soundings and gravity methods were employed to study the deep freshwater aquifer in the area north of Abo Zenema city on the eastern side of the Gulf of Suez, Egypt. Seven MT sites and 48 gravity stations were surveyed along northeast–southwest profiles as close as possible to a line perpendicular to the coast of the Gulf of Suez. The MT survey was conducted using high and low frequencies to investigate shallow and deep areas, respectively. One-dimensional inversion was conducted using a heuristic inversion scheme of the Bostick algorithm. The MT data were also inverted with a 2-D smooth model inversion routine using the nonlinear conjugate gradient method to infer variation in vertical and lateral resistivity inside the Earth. A 100-Ohm-m homogeneous half-space initial model was used to invert the TE mode data only. Then, the inverted model obtained from the TE mode data was used as an initial model for inversion of the TM mode data. The inverted model thus obtained from the TM mode data inversion was used as an initial model for the inversion of the joint TE and TM responses. Two-dimensional (2-D) forward modeling of the gravity data was conducted using the 2-D polygon method of Talwani’s algorithm for an arbitrarily shaped body and was based on the subsurface information from the MT survey and the available information about the geological structure of the study area. This method enabled us to obtain the basement structure of the coastal aquifer in the study area. The results from the analysis and the interpretation of MT and gravity data were used to detect and delineate the groundwater coastal aquifer in the study area.     

Deep structure of the northeastern margin of the Parnaiba Basin, Brazil, from magnetotelluric imaging

The magnetotelluric (MT) method has been applied to the determination of the deep resistivity structure of the northeastern margin of the Parnaiba Basin. Transient electromagnetic (TEM) and MT data were collected in early 1999 along a 95 km long N–S line, extending from the coast across the projected subcrop position of a discontinuous fault found to the west of the study area that is believed to be a possible basin‐bounding fault. The MT data were processed to yield the TE‐ and TM‐mode responses and then corrected for static shift using central‐loop and single‐loop TEM data, respectively. Regularized 2D MT inversion was subsequently undertaken using a structured initial model with the near‐surface constrained by TEM inversion results. As a consistency check, we performed another set of 2D inversions using different smooth initial models. The various optimal 2D inversion models show clearly the presence of a major basement trough, over 2 km deep, located about 70 km from the coast. We interpret it as possibly marking the main basin margin and suggest that it may have implications for groundwater resource development in the area.     

大地电磁资料精细处理和二维反演解释技术研究(三)——构建二维反演初始模型的印模法

在主流的线性最优化大地电磁二维反演中,如何合理构建初始模型是一个亟待解决的问题.常用的是采用均匀半空间或一维反演结果构建初始模型,不易获得稳定可靠的反演效果.实践表明,尽管基于不同初始模型的大地电磁二维反演结果差别较大,但均较初始模型更为接近真实模型.基于这样一种认识,经过反复的理论和实践探索,我们提出构建大地电磁二维反演初始模型的印模法.印模法的基本思想是依据已有反演结果和均匀半空间模型之间的加权来确定下一步二维反演的初始模型,它一方面保留了已有反演结果中关于真实模型的宏观轮廓信息,另一方面,保证了深部电性结构的均匀性,从而满足大地电磁二维正演所要求的底边界条件.基于印模法,本文进一步提出了迭代重构的反演思想.通过多个理论模型和实测数据的反演计算,验证了上述方法可在很大程度上压制初始模型对反演结果的影响.     

Mapping deep targets based on integrated 3D MT-gravity interpretation: a case study

Mapping deep geological hydrocarbon targets is of significant importance in basin exploration. In areas lacking reliable seismic data, magnetotelluric (MT) and gravity explorations are helpful to delineate the distribution of potential deep geological hydrocarbon targets. Here we investigate the effectiveness of the integrated 3D MT and gravity explorations for mapping the potential deep hydrocarbon source rocks. The result based on the data from the W Basin (part of the Ordes Basin) of China demonstrates that the method is efficient and economical for basin exploration. The method is particularly useful in target areas which are of great interest for oil and gas exploration but lack high quality seismic data. In our method, we first use the high-precision 3D small-bin MT data acquisition to improve the data accuracy. Then we perform datum static correction method and apply 3D inversion to obtain the3D resistivity distribution. We also develop a layered resistivity model based on resistivity logging to assist the interpretation of the inverted 3D resistivity data so as to derive an initial 3D geological model. Starting from the initial model, we use 2D gravity data to update the model via 2D inversion line by line, and then pass the updated model for the next round of the 3D MT inversion. The integrated inversion is implemented iteratively so the model converges to satisfy the need of final geological analysis. The application to the W Basin shows that we could successfully delineate the geological distribution of the potential deep hydrocarbon source rocks within the basin and map the thickness of the upper Paleozoic.     

APPLYING 3D INVERSION OF SINGLE-PROFILE MAGNETOTELLURIC DATA TO IDENTIFY THE SHADE AND YUNONGXI FAULTS

Many synthetic model studies suggested that the best way to obtain good 3D interpretation results is to distribute the MT sites at a 2D grid array with regular site spacing over the target area. However, MT 3D inversion was very difficult about 10 years ago. A lot of MT data were collected along one profile and then interpreted with 2D inversion. How to apply the state-of-the-art 3D inversion technique to interpret the accumulated mass MT profiles data is an important topic. Some studies on 3D inversion of measured MT profile data suggested that 2D inversions usually had higher resolution for the subsurface than 3D inversions. Meanwhile, they often made their interpretation based on 2D inversion results, and 3D inversion results were only used to evaluate whether the overall resistivity structures were correct. Some researchers thought that 3D inversions could not resolute the local structure well, while 2D inversion results could agree with the surface geologic features much well and interpret the geologic structures easily. But in the present paper, we find that the result of 3D inversion is better than that of 2D inversion in identifying the location of the two local faults, the Shade Fault(SDF)and the Yunongxi Fault(YNXF), and the deep structures. In this paper, we first studied the electrical structure of SDF and YNXF based on a measured magnetotelluric(MT) profile data. Besides, from the point of identifying active faults, we compared the capacity of identifying deep existing faults between 2D inversion models and 3D models with different inversion parameters. The results show that both 2D and 3D inversion of the single-profile data could obtain reasonable and reliable electrical structures on a regional scale. Combining 2D and 3D models, and according to our present data, we find that both SDF and YNXF probably have cut completely the high resistivity layer in the upper crust and extended to the high conductivity layer in the middle crust. In terms of the deep geometry of the faults, at the profile's location, the SDF dips nearly vertically or dips southeast with high dip angle, and the YNXF dips southeast at depth. In addition, according to the results from our measured MT profile, we find that the 3D inversion of single-profile MT data has the capacity of identifying the location and deep geometry of local faults under present computing ability. Finally, this research suggests that appropriate cell size and reasonable smoothing parameters are important factors for the 3D inversion of single-profile MT data, more specifically, too coarse meshes or too large smoothing parameters on horizontal direction of 3D inversion may result in low resolution of 3D inversions that cannot identify the structure of faults. While, for vertical mesh size and data error thresholds, they have limited effect on identifying shallow tectonics as long as their changes are within a reasonable range. 3D inversion results also indicate that, to some extent, adding tippers to the 3D inversion of a MT profile can improve the model's constraint on the deep geometry of the outcropped faults.     

标定大陆科学钻探孔区地震反射体

在中国大陆科学钻探孔区 (江苏省东海县南部 )进行了系统的地球物理调查 ,包括二维地震测网和专门的地震剖面 ,大地电磁法和位场方法等 .地震调查表明 ,在超高压变质岩出露区上地壳充满了反射体 ,包括倾斜反射体与上拱的反射弧 .本文介绍大陆科学钻探先行研究中地震调查的成果 .根据大陆科学钻探预先导孔的岩芯和测井资料、井旁VSP和数值模拟结果证实 ,高波速的榴辉岩体、破碎断裂带和大型韧性剪切带都可引起倾斜的地震反射 ,而上拱的弧形反射体则是由近似直立的榴辉岩体和其中的破碎带的综合反映 .由于地壳深部广泛分布着经受变质的岩石 ,上述研究结果对标定地壳中的反射地震信号具有一定意义     

A comparative study of results obtained in magnetotelluric deep soundings in Villarrica active volcano zone (Chile) with gravity investigations, distribution of earthquake foci, heat flow empirical relationships, isotopic geochemistry Sr/Sr and SB systematics

The first magnetotelluric deep soundings in Chile were carried out during 1986 in the Villarrica active volcano zone (39°25′S, 71°57′W). In the TM mode of polarization, the curves show a distorted segment with dispersion. A static distortion at long periods is observed in curves in the TE mode of polarization; the segment was shifted vertically to fit the geomagnetic global model values at daily periods. This modified curve was used for 1D modelling to determine the electrical structure in the study area. The upper level of the intermediate conducting layer of resistivity 20–60 ω m is found to be at 35–50 km depth. A higher resistivity layer (600 ω m), starting at 100 km depth, may be resolved in the intermediate conducting layer. A sharp decrease in the resistivity is shown by the model at 500 km.

Large heterogeneities at the level of the conducting layer encountered in the 1D modelling, and increased resistivity of the ultimate layer, may account for distortion observed at long periods. Two-dimensional test models show that the conducting layer in the area of Villarrica volcano may be an anomalous heated layer surrounded by rocks of higher resistivity of about 2 × 10³ ω m. These features correspond to the interaction with a subsiding oceanic lithosphere resulting in a complex thermal structure and perturbed resistivity distribution in transition zones of the Pacific type such as Chile, and to the existence of a megafault and a system of fractures in the sounding area. These facts make it difficult to determine the conductance of the electrical asthenosphere.

The parameters of the model structure correlate well with geophysical and geochemical results obtained in the area by other workers. Gravity studies indicate a maximum crustal thickness of about 37 km, which implies a non-full compensation according to the Airy hypothesis. The morphology of the Wadati-Benioff zone clearly shows a sharp decrease of earthquake foci at 50 km depth, and a reinitiation of seismic activity from a depth of 100 km down to 160 km. Kaufman and Keller and Levi and Lysak empirical relationships between heat flow and thickness of the lithosphere are somewhat consistent when the parameters of the model structure are assumed. The strontium isotope ratios are indicative of minimal crustal contamination of mantle-derived magmas, thus allowing a relatively thin crust under the area. Furthermore, the SB index of partial melting of mantle peridotite may suggest the contribution of uprising material from the deep asthenosphere through fractures extending to depth.     

Three-dimensional interpretation considering the static and the sea-effects of magnetotelluric data obtained in Jeju,Korea

Three-dimensional (3D) magnetotelluric (MT) surveys have been performed in Jeju, the largest volcanic island in Korea to figure out any possible structures or potential anomaly for remnant deep geothermal resources. Various approaches have been applied to interpret MT data observed in Jeju. MT dataset shows generally simple stratigraphy of four layers, though contains the severe static and the sea-effects. In our previous works, the induction vectors and 3D inversion results have commonly indicated the existence of a conductive anomaly in central parts of the island, beneath Mt. Halla. The 3D inversion dealt the static shifts as inversion parameters. The Jeju MT dataset, however, still contains the effect of conductive sea water surrounding the island.The sea-effect on MT impedance can be represented as a distortion tensor and excluded from the Jeju MT dataset by an iterative sea-effect correction. In this study, 3D inversion incorporating static shift parameterization was conducted using MT dataset corrected using 1D resistivity model obtained from the iterative scheme. Reasonably reconstructed images are obtained through the 3D inversion and using the MT dataset with sea-effect correction. The inversion result still shows the conductive anomaly in a similar depth. RMS misfits converged to a lower value than that of inversion using MT data before the sea-effect correction. From the fact, it is highly possible that the conductive anomaly is not an artifact but a real underground structure. Further investigation about the anomaly including exploration drilling is needed to see if it is from a fracture containing conductive sea water or related to the old volcanic activities.     

The changes in the electric conductivity of the lithosphere in the source region of the strongest Olyutora earthquake in the Koryak highlands

The results of the magnetotelluric (MT) soundings before and after an earthquake are analyzed. The interpretation is based on the longitudinal and transverse MT curves (along and across the strike of the main tectonic elements, respectively). The MT curves are distorted by the ρ- and coast effect. The distortions due to the coast effect are estimated by the testing three-dimensional (3D) model. It is established that the coast effect distortion at the periods up to 1000 s is small and can be disregarded. The divergence of the longitudinal and transverse MT curves, which points to the presence of the deep faults, is thoroughly studied. The inversion of the MTsounding curves is carried out by the REBOCC program of the numerical two-dimensional modeling. This program implements the procedures of elimination of the ρ-effect and the joint inversion of the longitudinal and transverse MT sounding curves. The obtained geoelectrical cross sections provide an insight into the structure of electrical conductivity of the lithosphere before and after the earthquake. The more intense variations in the electric conductivity are observed in the zone of the deep faults. These variations are related to the changes in the porosity and saturation of the rocks by the highly mineralized fluids.     

宁波深厚软土动力特性研究及其地震响应分析

利用Davidenkov模型对宁波深厚软土地区土动力试验数据进行拟合分析,确定相关拟合参数,结果表明该模型能较好地描述软土的动剪切模量和动阻尼比随动剪应变的变化规律。然后,基于等价线性化方法以及ANSYS软件的APDL语言进行简单的二次开发,编写Davidenkov模型计算程序,并利用经典的SHAKE软件进行验证,结果表明该程序计算结果很好,能解决目前商业软件中缺乏合适的土体非线性动力本构模型的不足。最后,利用开发的模型计算分析宁波深厚软土的地震响应规律,并与规范推荐的简化计算公式进行比较分析,结果表明规范推荐的公式得到的土体响应偏大,对于地下结构的抗震设计偏于保守。研究成果对确定深厚软土地区土体动力特性及其地震响应研究提供了合理的计算方法。     

长周期大地电磁测深研究——青藏高原东部LMT响应函数及应用

长周期大地电磁(LMT)是基于常规MT理论发展起来的电磁测深技术,青藏高原东部岩石圈较厚、视电阻率较低,应用LMT方法能够弥补常规MT仪器对低频信号响应的不足,获得有用的深部信息。文中概述了LMT的数据采集、处理及反演解释,并介绍了将LMT应用于东喜马拉雅构造结及其周围地区完成的长周期观测剖面——下察隅—玉树LMT剖面的岩石圈结构探测实例。实际观测表明,在重叠频段内LMT结果与常规MT具有良好的一致性;LMT数据处理实现的阻抗和倾子估计,在深部信息上具有更高的分辨率。LMT剖面结果揭示了青藏高原东部普遍存在地下低阻体,为局部地区部分熔融和地下流体存在的可能性提供了依据。LMT测深较好地弥补了常规MT方法在巨厚岩石圈和遇到低阻体时探测深度的局限性。     

面向目标自适应有限元法的带地形三维大地电磁各向异性正演模拟

传统三维大地电磁各向异性模拟均是基于规则六面体网格,计算精度有限且较难拟合复杂地质条件.本文采用面向目标自适应非结构矢量有限元法,对三维大地电磁各向异性介质进行模拟.首先从电场双旋度方程出发,利用伽辽金方法建立变分方程;然后利用电流密度连续性条件构建适合大地电磁各向异性问题的加权后验误差估计方法,实现面向目标的网格自适应正演;最后通过典型算例分析各向异性对网格自适应和大地电磁响应的影响特征以及各向异性的识别方法.本文算法能够高精度地拟合起伏地表和任意各向异性介质,适用于分析复杂地电条件大地电磁响应特征,为提高大地电磁资料解释水平提供了理论基础.     

The MT inversion for conductivity anisotropy and EDA precursor, stress field and deformation band in the Earth’s deep crust

The conductivity anisotropy behaviour is described for certain environment in the Earth’s crust and the MT inversion method for a layered symmetrically anisotropic model is presented. The inversion interpretations of the anisotropic model from the observational data are helpful to identify the earthquake precusors as indicated by the deep conductivity anisotropic variations, and also provide some useful information to investigate the stress states and deformation bands in the deep crust of the Earth. This project is sponsored by the Chinese Joint Seismological Science Foundation.     

广西龙滩库区深部孕震结构大地电磁探测研究

在广西龙滩库区布置2条北北西-南东东向大地电磁探测剖面进行了32个测点的观测,观测频段为320 Hz~1000 s.对2条剖面上测量和计算得到的各种电磁参数、深部电性结构图像以及库区地质构造、岩性和水域分布等进行分析研究.结果显示库区发育的6条断裂在电性上表现为不同程度的电性差异边界,其中拉浪-达良断裂和拉色-辉马断裂为库区主要断裂带.库区深部电性结构总体表现为纵向上自地表到20 km深度呈现高-低-高的波浪起伏的三层结构,横向上呈明显的高低电阻块状组合样式.龙滩库区4个地震丛集区的震源都位于高、低电阻接触带附近,其中3个地震丛集区的震源位于低阻特性的石炭系地层的下部,为岩溶水体诱发地震;发生最大震级的地震丛集区的震源位于上宽下窄似"铆钉状"的高阻体下部,推测是因水库蓄水后水体压力增大和库水渗透作用下在聚集高变形能的脆性高电阻体内部发生的地震.     

安徽贵池矿集区深部精细结构——来自综合地球物理探测结果的认识

安徽贵池矿集区是长江中下游成矿带中重要的铜多金属矿集区,随着浅表找矿难度的加大,急需了解矿集区深部结构为深部找矿提供深部找矿提供指示.为揭示贵池矿集区深部地壳结构格架、控岩、控矿构造和岩浆系统的空间结构,深化对成矿作用的认识,在该区实施了两条长度各约60km垂直于构造走向的反射地震和MT探测剖面,通过对MT数据有针对性的去噪和反演,获得了研究区10km以浅的地电结构;通过对反射数据的层析反演静校正、叠前保幅多域去噪、精细速度分析和叠前时间偏移等处理,获取了该地区地壳结构的地震成像剖面.结合区域地质构造背景和针对性处理的区域重磁数据,对两条综合探测剖面的反演处理结果进行了分析与解释,研究发现了本区地层的构造反射特征、地壳结构特征和Moho特征,厘定了盖层与基底之间的区域滑脱面、中下地壳滑脱面和Moho面的深度和形态.发现了地壳上冲和俯冲的反射特征,可能是华夏块体与扬子块体晋宁期发生板内碰撞的痕迹.高坦断裂是一个深达地壳级的断裂,其南北两侧具有不同的重磁特征、电性结构特征和地震反射特征,该断裂可能是深部流体上涌的通道,岩浆和地幔热液沿断裂上涌并沿地层和裂隙侵入中上地壳,形成岩基或岩体,或与围岩发生强烈矿化作用形成贵池矿集区多金属矿床.     

部分区域约束下的交叉梯度多重地球物理数据联合反演

交叉梯度联合反演方法通过对多种地球物理模型实现结构耦合,在岩石物性关系不确定的情况下,既能提高反演结果的可靠性,又能减少反演的多解性,还能减少不同方法解释结果之间的矛盾.当不同的模型观测数据覆盖范围不一致时,交叉梯度联合反演通常需要取出重叠区域数据进行联合反演,并且建模时还要扩展一些模型范围.本文首先提出并实现了部分区域约束下的交叉梯度多重地球物理数据联合反演算法;接着进行了算法的模型试算;最后,我们将该反演算法用于本溪—集安深部地质调查重磁电综合地质地球物理解释中.结果表明:该算法不但能在重叠区域内很好地恢复结构相似的模型,而且在非重叠区域与重叠区域的边界处仍然可以得到平滑变化的模型;在本溪—集安10号剖面所获得的结构上相似的电阻率、密度及磁化率模型较好地反映了该区的深部地质结构,对于确定深部地质体的性质提供了有力的证据.     

佳木斯地块及东缘岩石圈电性结构特征

佳木斯地块及东缘是中国东北地区的重要地质构造单元. 为探测该地区地壳深部结构与构造关系,沿桦南—饶河实施了240 km的大地电磁测深的探测研究. 采用光滑模型二维反演方法对桦南—饶河大地电磁剖面的探测数据进行了二维反演和综合地球物理解释. 研究结果揭示了研究区断面域的电性结构：(1)剖面西段具有稳定的高阻特征且具有稳定的岩石圈厚度（约90km）,在十几公里深度范围存在壳内高导层,该区段对应佳木斯地块. (2)剖面中部具有明显的电性梯度带,该梯度带反映了佳木斯地块的东界位置及其深部的结构形态. (3)剖面东段电性特征揭示了佳木斯地块以东区段浅部为逆冲推覆体,深部为多个高阻块体和低阻条带相间的电性结构,这些高阻块体可能为早期俯冲的岩石圈残片.