兰州地区深部电性结构的初步研究

利用兰州地区现有的大地电磁资料，对该区的深部电性结构进行了初步研究。结果表明，兰州地区各地质构造单元的电性差异较明显，与区域构造有较好的对应性。     

南北地震带南段川滇黔接壤区电性结构特征和孕震环境

本文利用大地电磁测深数据,对穿过兰坪—思茅地块和川滇菱形地块以及进入扬子地块的云南兰坪—贵州贵阳大地电磁测深剖面展开了深部电性结构研究.采用大地电磁数据处理分析以及反演技术,对观测资料进行了由定性到定量全面地分析,通过二维非线性共轭梯度反演得到了沿剖面的较为详细的地壳上地幔电性结构,结合其他地质和地球物理资料的分析,对该剖面的二维电性结构进行解释,确定了主要断裂带和边界带的位置以及深部延伸情况,同时确定了壳内低阻层的分布位置,最后进行了区域动力学和孕震构造环境的探讨.研究表明:剖面壳幔电性结构分块性特征与区域地质构造分布特征基本一致,不同地块的电性结构存在显著差异,其中川滇菱形地块的结构相对复杂,上地壳的电性结构为高低阻相间分布特征,电阻率的突变带与地表断裂具有很好的对应关系;兰坪—思茅地块存在中上地壳低阻层,川滇菱形地块中西部存在下地壳低阻层,川滇菱形地块东部和华南地块西部存在中上地壳的低阻层;川滇菱形地块中部攀枝花附近的低阻层埋深最深,而华南地块西部会泽附近的低阻层埋深则最浅;兰坪—思茅地块和川滇菱形地块的中下地壳的低阻层可能与青藏高原物质的东南逃逸有关;华南块体的宣威以东的下地壳不存在低阻层,华南块体下地壳和上地幔的电阻率较高;攀枝花附近的高阻体可能是峨眉山玄武岩喷发导致底侵作用及幔源物质上侵的结果.     

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

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

基于电性结构模型的青藏高原东缘上地幔热结构研究

研究青藏高原东缘地区的深部物质结构对于理解青藏高原的隆升及扩张机制具有重要的科学意义.本文将青藏高原东缘实测大地电磁测深剖面反演所得的岩石圈电性结构模型与高温高压岩石物理实验测得的上地幔矿物和熔融体导电性定量关系相结合,通过Hashin-Shtrikman(HS)边界条件建立上地幔电导率与温度、熔融百分比等参数的定量关系,在此基础上计算得到了青藏高原东缘上地幔热结构及熔融百分比分布模型.研究结果表明在青藏高原东缘地区通过大地电磁测深方法所探测到的上地幔低阻体可以解释为由高温作用所产生的局部熔融区域.其中,松潘—甘孜地块上地幔高导体对应的温度介于1300～1500℃之间,熔融百分比可高达10%,支持前人将松潘—甘孜地块内部的低阻体解释为局部熔融的观点.龙门山断裂带以东、四川盆地西缘的上地幔高导体温度介于1200～1400℃之间,熔融百分比介于1%～5%左右,表明扬子克拉通的西缘可能正在经历一定程度的活化作用.龙门山断裂带下方的上地幔高阻体温度介于1100℃附近,基本没有发生局部熔融,具有较冷的刚性块体特征,与该区域频发的地震活动相吻合.四川盆地东部的扬子上地幔温度介于800～900℃之间,没有发生局部熔融,符合古老稳定的克拉通块体的基本特征.     

Numerical simulation and inversion of MT fields in the 3D electric conductivity model of the Vesuvius volcano

Possibilities for three-dimensional (3D) magnetotelluric (MT) sounding of local objects contained in the Earth’s crust are estimated in a case study of the magma chamber of the Vesuvius volcano. Stochastic inversion of the model MT data by the Markov Chain Monte Carlo (MCMC) method has shown that the most efficient approach is not simultaneous but successive estimation of the geometry and the depth of the anomaly and the assessment of the conductivity distribution within the anomalous region. A zone of equivalence is revealed between the a priori estimate of the depth of the anomalous zone and the a posteriori distribution of electric conductivity within it. Based on the present estimation and previous results, an algorithm for determination of the parameters of local crustal anomaly is proposed.     

On the Coupling of Geodynamic and Resistivity Models: A Progress Report and the Way Forward

Magnetotelluric (MT) studies represent the structure of crust and mantle in terms of conductivity anomalies, while geodynamic modelling predicts the deformation and evolution of crust and mantle subject to plate tectonic processes. Here, we review the first attempts to link MT models with geodynamic models. An integration of MT with geodynamic modelling requires the use of relationships between conductivity and rheological parameters such as viscosity and melt fraction, which are provided by laboratory measurements of rock properties. Owing to present limitations in our understanding of these relationships, and in interpreting the trade-off between scale and magnitude of conductivity anomalies from MT inversions, most studies linking MT and geodynamic models are qualitative rather than providing hard constraints. Some recent examples attempt a more quantitative comparison, such as a study from the Himalayan continental collision zone, where rheological parameters have been calculated from a resistivity model and compared to predictions from geodynamic modelling. We conclude by demonstrating the potential in combining MT results and geodynamic modelling with examples that directly use MT results as constraints within geodynamic models of ore bodies and studies of an active volcano-tectonic rift.     

2D modelling of resistivity and magnetotelluric data from the Belvedere Spinello salt mine,Italy1

The Belvedere Spinello salt mine is located in the Catanzaro Province of Calabria in Southern Italy. An extensive mining program has caused the development of Underground cavities filled with brine and the migration of this brine has been of great environmental concern to the mine owners. This paper presents the results of a multidimensional interpretation of a two-phase resistivity and magnetotelluric (MT) survey that was performed in an attempt to determine the complex conductivity structure of the mine area and to gain information on brine development and migration pathways. Key resistivity soundings were interpreted using a 2.5D algorithm based on the Polozhii decomposition method. The MT data were interpreted using a 2D finite-element code. A conductivity model was developed, integrating available geological and drill-hole information. The interpretation of the MT data, collected five years after the acquisition of the resistivity data, shows a conductive feature of depth that is not resolved in the resistivity interpretation. This feature has been interpreted as a thick brine zone that has developed as a result of mining during the interval between the resistivity and the MT measurements.     

郯庐断裂带中段(沂沭断裂带)电性结构研究与孕震环境

穿过郯庐断裂带中段（沂沭断裂带,36°N）所做的大地电磁测深（MT）剖面长约150km.使用Robust技术和远参考道大地电磁方法处理观测数据.通过分析视电阻率、阻抗相位、Swift二维偏离度和区域走向,定性确定测区的电性结构.二维反演解释中选择非线性共轭梯度（NLCG）方法,使用TE、TM两种模式资料联合反演,沿剖面的二维电性结构显示：自西向东,鲁西隆起、郯庐断裂带、胶莱坳陷及鲁东隆起4个电性区块分别对应,鲁东和鲁西隆起区为高阻,郯庐断裂带电性结构复杂,高、低阻相间,胶莱坳陷为低阻（高导）区.沿MT剖面附近曾发生3个地震,其震源区处在电性变化剧烈部位,并在震源区附近存在高导体.     

Large-scale three-dimensional inversion of EarthScope MT data using the integral equation method

In this paper we apply 3D inversion to MT data collected in the Northwestern United States as a part of the EarthScope project. By the end of 2009 MT data had been collected from 262 stations located throughout Oregon, Washington, Idaho, and most of Montana and Wyoming. We used data from 139 MT stations in this analysis. We developed fully parallelized rigorous 3D MT inversion software based on the integral equation method with variable background conductivity. We also implemented a receiver footprint approach which considerably reduced the computational resources needed to invert the large volumes of data covering vast areas. The data set used in the inversion was obtained through the Incorporated Research Institutions for Seismology (IRIS). The inversion domain was divided into 2.7 M cells. The inverted electrical conductivity distribution agrees reasonably well with geological features of the region.     

Deep structure of Western Yakutia from magnetotelluric soundings

Analysis of magnetotelluric (MT) sounding studies in western Yakutia showed that amplitude curves are strongly affected by the regional distribution of electrical conductivity heterogeneities of the sedimentary cover, while phase curves are distorted to a lesser extent. The distribution of the longitudinal conductance in the sedimentary cover is constructed, and MT fields and resistivity and phase curves are modeled. Based on combined analysis of experimental and model data and results of inversion of corrected curves, an interpretation of main anomalies of MT fields is given and a regional crustal anomaly of the electrical conductivity is discovered.     

带地形的大地电磁各向异性二维模拟及实例对比分析

带地形的大地电磁二维正演数值模拟多数基于电性各向同性理论,由于地球内部电性各向异性现象的普遍存在,基于电性各向异性理论研究地形起伏情况下大地电磁二维正演数值模拟就显得非常迫切.本文首先由麦克斯韦方程出发,引入张量电导率,求得一组关于平行走向的电场分量Ex和磁场分量Hx的二阶偏微分方程,使用有限差分法求解出Ex和Hx的近似解,并以此求得其他场分量;其次,引入地形因素,改变变量在网格节点中的排列方式,选择交错排列方式从而给有限差分系数矩阵的最大带宽分配合理的存储空间;最后,使用Weaver的方法解决TM模式下,在地-空分界面垂直于构造走向的一些区域存在不同电导率的问题.通过对带地形的二维电性各向异性结构做正演模拟,研究地形因素对大地电磁响应的影响;以电性各向异性理论为基础,将地形因素引入对实测大地电磁资料的处理中,通过做二维正演拟合和未引入地形因素的结果做对比,说明电性各向异性现象的普遍存在,认识地形因素对观测大地电磁场的影响,为今后分析解释实测大地电磁资料包含地形因素和电性各向异性情况提供理论基础和技术指导.     

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.     

3-D inversion of MT data from the Sabalan geothermal field,Ardabil, Iran

Since the true Earth is 3-D in nature, a three-dimensional (3-D) inversion has clear advantages over lower dimensional inversions. We utilized a 3-D magnetotelluric (MT) inversion code, WSINV3DMT, to obtain a realistic resistivity model using a long period MT data set collected in the Northwest Sabalan geothermal field in Ardabil, Iran. The apparent resistivity and phase curves, the magnetic induction vectors, the impedance polar diagrams and the rotational invariant of impedance tensor, indicate a complex 3-D conductivity structure. After setting up the model parameters and designing the appropriate block discretization, we performed the 3-D inversions for two sets of observed data; one set includes the full MT impedance tensor and another set contains only off-diagonal elements of the MT impedance. The final model was selected according to the relative magnitude of the data misfit and the model norm with respect to various Lagrangian multipliers. The results of this study illustrate the 3-D inversion of the off-diagonal elements of MT impedance tensor is precisely enough to explain the structures related to the geothermal source. The obtained results were compared with the results of available 2-D models and they are then interpreted using all of the geological and drilling data of the area. The main outcome of this study is the precise delineation of the geometry of geothermal source that is located at the center of the study area with a surface coverage of about 7 km².     

On the Causes of Electrical Conductivity Anomalies in Tectonically Stable Lithosphere

Magnetotelluric (MT) data can image the electrical resistivity of the entire lithospheric column and are therefore one of the most important data sources for understanding the structure, composition and evolution of the lithosphere. However, interpretations of MT data from stable lithosphere are often ambiguous. Recent results from mineral physics studies show that, from the mid-crust to the base of the lithosphere, temperature and the hydrogen content of nominally anhydrous minerals are the two most important controls on electrical conductivity. Graphite films on mineral grain boundaries also enhance conductivity but are stable only to the uppermost mantle. The thermal profile of most stable lithosphere can be well constrained, so the two important unknowns that can affect the conductivity of a lithospheric section are hydrogen content and graphite films. The presence of both of these factors is controlled by the geological history of the lithosphere. Hydrogen in nominally anhydrous minerals behaves as an incompatible element and is preferentially removed during melting or high-temperature tectonothermal events. Grain-boundary graphite films are only stable to ~900 °C so they are also destroyed by high-temperature events. Conversely, tectonic events that enrich the lithosphere in incompatible elements, such as interaction with fluids from a subducting slab or a plume, can introduce both hydrogen and carbon into the lithosphere and therefore increase its electrical conductivity. Case studies of MT results from central Australia and the Slave Craton in Canada suggest that electrical conductivity can act as a proxy for the level of enrichment in incompatible elements of the lithosphere.     

喜马拉雅西北部逆冲带的地壳电性结构

印度板块北部地形起伏较大的喜马拉雅山地区由几个构造互异的地质单元组成,依地形高、低把喜马拉雅碰撞带分成低喜马拉雅和高喜马拉雅.为了研究与主要逆冲带（含主缝合带MCT和主边界带MBT）有关的地壳电性结构,沿Rohtangpass (海拔4000 m) 到Mandi (海拔400 m)剖面进行了MT探测.通过对16个测点观测资料的分析和考虑地形的二维反演,获得了沿剖面的二维电性结构.电性结构显示,在Chail和主逆冲边界带下方,东西走向的缝合带突然转向北.在下喜马拉雅的Rampur 区段的元古代基底为范围较大的高阻体,而浅部地壳被逆冲带分成向北倾的电导性块体和电阻性块体.Chail 逆冲带东侧低喜马拉雅Rampur 区段的推挤和它西侧的基底脊柱体导致主边界带及相关的逆冲带（Kangra 拐角）向北转弯,Kangra拐角处的应力可能是由于西侧基底脊柱体进入到Kangra 区引起的.     

起伏地形下大地电磁L-BFGS三维反演方法

为推进大地电磁三维反演的实用化,本文实现了基于L-BFGS算法的带地形大地电磁三维反演.首先推导了大地电磁法三维反演的Tikhonov正则化目标函数以及Hessian矩阵逆矩阵近似表达式和计算方法,然后设计了一种既能保证空气电阻率固定不变又能保证模型平滑约束的协方差矩阵统一表达式,解决带地形反演问题.在反演算法中采用正则化因子冷却法以及基于Wolf条件的步长搜索策略,提升了反演的稳定性.利用开发的算法对多个带地形地电模型（山峰地形下的单个异常模型、峰-谷地形下的棋盘模型）的合成数据进行了三维反演,并与已有大地电磁三维反演程序（ModEM）进行对比,验证了本文开发的三维反演算法的正确性和可靠性.最后,利用该算法反演了华南某山区大地电磁实测数据,得到该区三维电性结构,揭示了研究区以高阻介质为基底,中间以低阻不整合面和相对低阻介质连续分布,浅部覆盖高阻介质的电性结构特征,进一步验证了本文算法的实用性.     

ELECTRICAL STRUCTURE OF THE 2017 MS7.0 JIUZHAIGOU EARTHQUAKE REGION AND THE EASTERN TERMINUS OF THE EAST KUNLUN FAULT

The East Kunlun Fault is a giant fault in northern Tibetan, extending eastward and a boundary between the Songpan-Ganzi block and the West Qinling orogenic zone. The East Kunlun Fault branches out into a horsetail structure which is formed by several branch faults. The 2017 Jiuzhaigou M_S7.0 earthquake occurred in the horsetail structure of the East Kunlun Fault and caused huge casualties. As one of several major faults that regulate the expansion of the Tibetan plateau, the complexity of the deep extension geometry of the East Kunlun Fault has also attracted a large number of geophysical exploration studies in this area, but only a few are across the Jiuzhaigou earthquake region. Changes in pressure or slip caused by the fluid can cause changes in fault activity. The presence of fluid can cause the conductivity of the rock mass inside the fault zone to increase significantly. MT method is the most sensitive geophysical method to reflect the conductivity of the rock mass. Thus MT is often used to study the segmented structure of active fault zones. In recent years MT exploration has been carried out in several earthquake regions and the results suggest that the location of main shock and aftershocks are controlled by the resistivity structure. In order to study the deep extension characteristics of the East Kunlun Fault and the distribution of the medium properties within the fault zone, we carried out a MT exploration study across the Tazang section of the East Kunlun Fault in 2016. The profile in this study crosses the Jiuzhaigou earthquake region. Other two MT profiles that cross the Maqu section of East Kunlun Fault performed by previous researches are also collected. Phase tensor decomposition is used in this paper to analyze the dimensionality and the change in resistivity with depth. The structure of Songpan-Ganzi block is simple from deep to shallow. The structure of West Qinlin orogenic zone is complex in the east and simple in the west. The structure near the East Kunlun Fault is complex. We use 3D inversion to image the three MT profiles and obtained 3D electrical structure along three profiles. The root-mean-square misfit of inversions is 2.60 and 2.70. Our results reveal that in the tightened northwest part of the horsetail structure, the East Kunlun Fault, the Bailongjiang Fault, and the Guanggaishan-Dieshan Fault are electrical boundaries that dip to the southwest. The three faults combine in the mid-lower crust to form a "flower structure" that expands from south to north. In the southeastward spreading part of the horsetail structure, the north section of the Huya Fault is an electrical boundary that extends deep. The Tazang Fault has obvious smaller scale than the Huya Fault. The Minjiang Fault is an electrical boundary in the upper crust. The Huya Fault and the Tazang Fault form a one-side flower structure. The Bailongjiang and the Guanggaishan-Dieshan Fault form a "flower structure" that expands from south to north too. The two "flower structures" combine in the high conductivity layer of mid-lower crust. In Songpan-Ganzi block, there is a three-layer structure where the second layer is a high conductivity layer. In the West Qinling orogenic zone, there is a similar structure with the Songpan-Ganzi block, but the high conductivity layer in the West Qinling orogenic zone is shallower than the high conductivity layer in the Songpan-Ganzi block. The hypocenter of 2017 M_S7.0 Jiuzhaigou earthquake is between the high and low resistivity bodies at the shallow northeastern boundary of the high conductivity layer. The low resistivity body is prone to move and deform. The high resistivity body blocked the movement of low resistivity body. Such a structure and the movement mode cause the uplift near the East Kunlun Fault. The electrical structure and rheological structure of Jiuzhaigou earthquake region suggest that the focal depth of the earthquake is less than 11km. The Huya Fault extends deeper than the Tazang Fault. The seismogenic fault of the 2017 Jiuzhaigou earthquake is the Huya Fault. The high conductivity layer is deep in the southwest and shallow in the northeast, which indicates that the northeast movement of Tibetan plateau is the cause of the 2017 Jiuzhaigou earthquake.     

用大地电磁法研究构造走向及维性特征

介绍了大地电磁GB张量分解法及其对它的改进法 ,可确定出更可靠更真实的区域构造走向 .将分解结果结合传统的座标旋转法所确定的视电阻率、相位、走向、偏离度等响应函数及维权参数进行分析 ,可得到更详细的电性结构维性质信息 .对兰州地区的实测资料研究表明 ,区域电性结构主体呈 2 D结构 ,走向方向大致为南北或东西向     

高温高压下麻粒岩电导率研究

麻粒岩是下地壳的重要组成物质,模拟其在下地壳温压条件下的电导率对于认识下地壳的电导率分布具有十分重要的意义.本文用交流法在1 GPa下,373~1002 K范围内研究了麻粒岩样品的复阻抗,并且对测量结果进行了阻抗谱分析.研究结果表明,复阻抗对频率具有依赖性,且随温度的升高,复阻抗的实部、虚部均变小.在实验给出的温度范围内,电导率结果符合Arrhenius关系式.当温度在373~663 K范围内时,实验所获得的激化焓为0.31 eV,表明样品的电导率由低能带杂质离子所控制;当温度在673~1002 K范围内时,激化焓为0.67 eV,此时可能为小极化子导电.将所得电导率结果与西南峨边-马边地区以及华北应县―商河地区的大地电磁结果进行了对比,发现在所模拟的下地壳温压范围内,实验室测得的电导率位于野外MT数据范围内.     

The North American Central Plains conductivity anomaly and its correlation with gravity, magnetic, seismic, and heat flow data in Saskatchewan, Canada

The North American Central Plains conductivity anomaly (NACP) lies, virtually in its entirety, within the Trans-Hudson Orogen. Accordingly, should these two features prove to be contemporaneous, then the geometrical relationship between these two is of foremost importance to any evolutionary tectonic model proposed to explain the collision of the Superior and Churchill Provinces in the Hudsonian, and a model that does not include a mechanism for the generation of this anomaly is obviously untenable.

In order to map better the trend of the NACP in the Province of Saskatchewan, Canada, two magnetotelluric (MT) profiles were conducted over the NACP previously defined by magnetometer arrays and profiles. Data from these two profiles, along with an earlier MT profile just north of the U.S./Canadian border, suggest that the NACP is not a continuous feature, but rather that it exhibits a definite break at latitude 51° N. Other geophysical evidence examined herein is concordant with this characteristic. If a second MT anomaly mapped to the northwest of this break is, in fact, a manifestation of the same geological structure, then one possible interpretation is of a major NW-SE trending sinistral fault in the deep crust, previously undetected, with a movement of some 100–150 km along strike.

The MT data from the southernmost profile, after correction for static shift, are modelled in a 2D manner, and it is shown that the NACP is consistent with an arcuate structure in vertical section of high conductivity (> 2 S m⁻¹) beginning at a depth of 10 km centred on 103° W dipping down to the west reaching possibly the base of the crust. It is also shown that such an arcuate shape in section can explain an observed gravity high of 40 mgal. Such high conductivities cannot be explained in terms of connected fluids, as it would require implausibly high porosities (12–20%). A comprehensive interpretation of all the geophysical data would have to account for the observed characteristics of the anomalous region; which are, high electrical conductivity, positive density contrast, no magnetization, and high heat flow in the basement, in terms of a single causative body. One possible explanation is presented in terms of a zone of lithospheric weakness along which was emplaced low-density differentia from a mantle-derived body. An alternative explanation could be the phase transformation of obducted material into eclogite or serpentinite. However, difficulties exist with both of these explanations and further data are required.