海岸效应对近海地区大地电磁测深数据畸变作用研究

在近海地区采集的大地电磁测深数据通常受到海岸效应的影响,使得大地电磁测深数据发生畸变,因而很难利用大地电磁测深资料较为可靠地获得地下深部的电性结构.本文通过正演模拟方法,分析和总结海水深度变化和海底地形变化对近海地区大地电磁测深数据的畸变影响.当测区与海岸线的距离小于目标频率的大地电磁场趋肤深度时,高导海洋的存在会严重影响测区内电磁场的分布.由于海岸效应的影响,大地电磁测深视电阻率曲线和相位曲线均会发生不同程度的畸变,在低频部分,这种畸变作用尤为明显.大地电磁测深一维Occam反演方法和二维非线性共轭梯度反演方法,对近海地区浅部地层具有较好的反演效果.随着海水深度的增加和海底地形的复杂变化,两种反演方法均会出现不同程度的假异常,为地质解释工作造成了影响.近渤海地区的实测大地电磁测深数据在低频部分可能受到海岸效应的影响而导致视电阻率曲线的严重畸变.     

大地电磁局部畸变异常特征三维数值模拟研究

局部畸变问题曾经困扰大地电磁资料反演解释几十年,大地电磁三维数值模拟技术的发展为剖析局部畸变特点和得到可靠的反演成像结果提供了技术基础。本文采用三维数值模拟成像方法对典型三维局部畸变模型进行模拟分析。三维数值模拟结果显示:电场分量垂直电性分界面的极化模式视电阻率曲线(对应二维情况下TM模式)在穿越低阻异常体界面时,曲线会先上移后下移,而在穿越高阻异常体界面时,曲线会先下移后上移,这与电性分界面处积累面电荷产生的二次电场有关。三维模型中XY模式、YX模式视电阻率和相位在三维异常体附近的水平变化是呈现近似垂向对称的,该现象与电场垂直跨越电性界面时视电阻率的变化规律是吻合的,当测线分别沿X方向和Y方向展布时,三维情况下的XY和YX模式分别对应二维情况下的TM模式。低阻小异常体对区域构造响应的畸变影响比高阻小异常体要严重。低阻小异常体对二维区域响应的两种极化模式视电阻率和相位都有非常明显的畸变影响,相比较而言对TE模式的畸变要大于TM模式,因此我们在做二维反演解释时,可优先考虑拟合TM模式数据。位于小异常体中心上方测点的三维畸变响应虽然与对应真实二维区域响应的差异比较大,但可以等效于某种二维模型响应,这种由局部畸变造成的假二维响应在实际野外数据的解释中是需要注意的。      

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.     

九瑞矿集区三维电性结构研究及找矿意义

九瑞矿集区作为全国重点整装勘查区之一,深部找矿工作一直是研究重点和热点.本文在综合分析研究区区域地质资料、对典型地层和岩石标本进行电阻率参数测试和统计分析的基础上,利用研究区面积性AMT数据进行三维电性结构研究.为获得AMT最佳反演参数,选择一条典型AMT剖面开展二维非线性共轭梯度(NLCG)反演,进而总结出利用TM极化模式数据、拉格朗日乘子等于3的二维反演方案.综合利用地质、钻探和电性特征对全区23条AMT剖面进行二维约束反演,并将所有测线的电阻率反演结果三维网格化,得到了研究区的三维电性结构.最后,我们通过对三维电性结构的综合分析,确定研究区地下地层、岩体展布、断裂构造分布特征以及基底起伏情况.与已知矿床对比,我们发现浅部北西断裂与北东向区域性深部断裂交汇处、基底隆起凹槽边缘以及岩体边缘为成矿有利区.     

松嫩地块东缘和佳木斯地块西缘电性结构

横过松嫩地块东缘和佳木斯地块西缘的大地电磁测深剖面揭示了两块体结合带附近的深部电性结构.本文对剖面测点做了标准化数据处理,并对二维偏离度、构造走向进行了计算和分析,采用非线性共轭梯度(NLCG)二维反演方法对TM模式的数据进行了反演,获得了该剖面的地壳、上地幔电性结构模型,划分出三个典型构造单元:松嫩地块东缘、碰撞拼合带和佳木斯地块西缘.研究结果表明,研究区上地壳基本呈高阻特征,可能为岩浆岩,代表其经历了多期次岩浆作用,而松嫩地块东缘和佳木斯地块西缘的中下地壳的高导体可能与地幔物质的上涌有关;拼合带下方存在西倾的高导体和高阻体,可能是佳木斯地块向西俯冲到松嫩地块下方的构造遗迹;研究区可能发生了拆沉作用,与之伴随的地幔物质上涌可能是后期伸展作用的一个动力.     

Exploration of geothermal structure in Puga geothermal field, Ladakh Himalayas, India by magnetotelluric studies

To understand the crustal electric structure of the Puga geothermal field located in the Ladakh Himalayas, wide band (1000 Hz–0.001 Hz) magnetotelluric (MT) study have been carried out in the Puga area. Thirty-five MT sites were occupied with site spacing varying from 0.4 to 1 km. The measurements were carried out along three profiles oriented in east–west direction. After the preliminary analysis, the MT data were subjected to decomposition techniques. The one-dimensional inversion of the effective impedance data and the two-dimensional inversion of the TE (transverse electric) and TM (transverse magnetic) data confirm the presence of low resistive (5–25 Ω m) near surface region of 200–300 m thick in the anomalous geothermal part of the area related to the shallow geothermal reservoir. Additionally, the present study delineated an anomalous conductive zone (resistivity less than 10 Ω m) at a depth of about 2 km which is possibly related to the geothermal source in the area. A highly resistive basement layer separates the surface low resistive region and anomalous conductive part. The estimated minimum temperature at the top of conductive part is about 250 °C. The significance of the deeper conductive zone and its relation to the geothermal anomaly in the area is discussed.     

三河-平谷8级大震区地壳上地幔电性结构特征研究

用电磁阵列剖面法 (EMAP)、大地电磁测深方法 (MT) ,沿三河 -平谷 8级大震震源区 ,作了 31 8km长的EMAP探测和两条总长 150 0 5km共 36个点的MT探测。获得了研究范围内的地壳上地幔电性结构、高导层特征和陡变带、高导异常体、断裂展布、岩石圈结构等结果 ,为搞清地震危险区的深浅构造关系、从电性结构特征推测发震模式和预测未来强震的可能地点提供了介质电性的多种参数     

大兴安岭与两侧盆地结合地带深部电性结构与岩石圈尺度构造关系

横跨大兴安岭与海拉尔盆地和松辽盆地结合地带的大地电磁测深剖面揭示了盆山构造的深部电性结构.剖面西起海拉尔盆地东缘,向东延伸穿过大兴安岭中部,一直到达松辽盆地西缘.本文对剖面测点的二维偏离度、构造走向等进行了计算和分析,采用非线性共轭梯度(NLCG)二维反演方法对TM模式的数据进行了反演,获得了该剖面的地壳、上地幔电性结构模型,划分出三个典型构造单元:海拉尔盆地、大兴安岭和松辽盆地.研究结果表明,海拉尔盆地东缘和松辽盆地西缘浅部都呈低阻特征,但松辽盆地西缘深部电性结构比较复杂,而大兴安岭整体呈高阻特征.海拉尔盆地东缘可能属于兴安块体,松辽盆地西缘与大兴安岭接触关系复杂.海拉尔盆地东缘岩石圈厚度约为110km,大兴安岭岩石圈厚度约为110~150km.大兴安岭上地壳基本呈高阻特征,可能为多次叠置的岩浆岩,代表大兴安岭经历了多期次岩浆作用;中下地壳横向存在较大范围低阻体,可能反映了大兴安岭地壳内部非刚性的特点;残存在岩石圈地幔的高阻异常,说明其下地壳可能发生过拆沉作用.大兴安岭与松辽盆地结合带存在一个岩石圈尺度的西倾低阻带,向下延伸到岩石圈底部,可能是早期松嫩地块向兴安地块俯冲并以软碰撞形式拼合的构造遗迹.     

青藏高原东北缘至鄂尔多斯地块壳幔电性结构及构造变形研究

为了获取青藏高原东北缘至鄂尔多斯地块的壳幔电性结构,研究祁连造山带、鄂尔多斯地块及六盘山构造带的构造变形,布设一条甘肃陇西至陕西黄陵的近东西向大地电磁测深剖面,获取了91个大地电磁测深点的响应.经过对全剖面观测资料的数据处理、分析及二维反演,获得了剖面壳幔电性结构模型.研究结果表明:剖面横向可划分为三个区块,分别对应祁连造山带、六盘山构造带与鄂尔多斯地块;祁连造山带东段可能残存沟弧盆体系的构造格架,青藏高原北东向生长可能是在这一先存格架上的叠加与改造;六盘山构造带壳幔结构复杂,以中地壳拆离断层为界,上地壳发育拆离断层系统而下地壳挤压缩短增厚;鄂尔多斯地块成层性较好,地块总体较为稳定,但局部经历了与地幔上涌相关的物质与结构再造.     

大地电磁资料精细处理和二维反演解释技术研究(七)——云南盈江——龙陵地震区深部电性结构及孕震环境

本文对一条布设在滇西盈江—龙陵地区的大地电磁剖面(苏典—中山剖面)数据进行了精细处理和二维反演解释,得到了测区较高置信度的二维电性结构.该电性模型纵向上表现为高阻-低阻-高阻的"三明治"式岩石圈电性结构,上地壳为平均厚度约为10km的高阻地层,在约6～16km地壳深度范围发育有电阻率为几欧姆米的显著高导层,下地壳底部和上地幔顶部表现为电性较为均匀的相对高阻层.横向上自西向东划分出以大盈江断裂带、龙陵—瑞丽断裂带为限的3个主要构造区域.壳内分布的高导层沿剖面表现出一定的横向不均匀性,其在龙陵—瑞丽断裂带下方消失,在该处形成了腾冲地块和保山地块的电性构造边界.电性结构表明,大盈江断裂附近高导层顶界面浅,两侧高阻体厚度小,因此难以形成较大规模的相互作用,致其附近浅震源、小震级的地震活跃;龙陵—瑞丽断裂两侧的高阻体较厚,易积累较大的应力,具有大震的深部孕震环境,故其附近发生过多次7级以上强震.     

河北石家庄地区深部结构大地电磁探测

石家庄位于太行山隆起带和华北平原冀中坳陷盆地接触区,西邻太行山山前断裂带,1966年在其东南曾发生邢台7.2级强震.为研究该区的深部构造背景,并为分析地震活动性趋势提供基础资料,2010年10月采用宽频带大地电磁法对该区的深部结构进行探测研究.大地电磁剖面穿过石家庄南侧区域,长约167km,获得了64个测点数据.采用远...     

Structural controls on the 1998 volcanic unrest at Iwate volcano: Relationship between a shallow, electrically resistive body and the possible ascent route of magmatic fluid

Magnetotelluric (MT) measurements were conducted at Iwate volcano, across the entirety of the mountain, in 1997, 1999, 2003, 2006, and 2007. The survey line was 18 km in length and oriented E–W, comprising 38 measurements sites. Following 2D inversion, we obtained the resistivity structure to a depth of 4 km. The surface resistive layer (~ several hundreds of meters thick) is underlain by extensive highly conductive zones. Based on drilling data, the bottom of the highly conductive zone is interpreted to represent the 200 °C isotherm, below which (i.e., at higher temperatures) conductive clay minerals (smectite) are rare. The high conductivity is therefore mainly attributed to the presence of hydrothermally altered clay. The focus of this study is a resistive body beneath the Onigajo (West-Iwate) caldera at depths of 0.5–3 km. This body appears to have impeded magmatic fluid ascent during the 1998 volcanic unrest, as inferred from geodetic data. Both tectonic and low-frequency earthquakes are sparsely distributed throughout this resistive body. We interpret this resistive body as a zone of old, solidified intrusive magma with temperatures in excess of 200 °C. Given that a similar relationship between a resistive body and subsurface volcanic activity has been suggested for Asama volcano, structural controls on subsurface magmatic fluid movement may be a common phenomenon at shallow levels beneath volcanoes.     

Papua New Guinea MT: looking where seismic is blind1

Hydrocarbon exploration in the Papuan fold belt is made extremely difficult by mountainous terrain, equatorial jungle and thick karstified Miocene limestones at the surface. The high-velocity karstified limestones at or near the surface often render the seismic technique useless for imaging the subsurface. In such areas magnetotellurics (MT) provides a valuable capability for mapping subsurface structure. The main structural interface which can be mapped with MT, due to the large electrical contrast, is the contact between the resistive Darai limestone and the underlying conductive sediments of the Ieru Formation. In some areas the base of the Darai can be mapped with reasonable accuracy by fitting 1D models to the observed MT data. However, in many cases where 2D and 3D effects are severe, 1D interpretations can yield dramatically incorrect results. Numerical and field data examples are presented which demonstrate the severity of the 1D errors and the improvements in accuracy which can be achieved using a 2D inverse solution. Two MT lines over adjacent anticlines, both with well control and seismic data, are used to demonstrate the application of 1D and 2D inversions for structural models. In both cases the seismic data provide no aid in the interpretations. The example over the Hides anticline illustrates a situation where 1D inversion of either TE or TM mode provides essentially the same depth to base of Darai as 2D inversion of both TE and TM. Both models provide base Darai depth estimates which are within 10% of that measured in the well. The example over the Angore anticline illustrates the inadequacy of 1D inversion in structurally complex geology complicated by electrical statics. The TE mode fits a 1D Darai thickness of 800 metres while the TM mode fits a 1D Darai thickness of 3500 metres, bracketing the thickness of 2450 metres observed in the well. The final 2D inversion model provides a depth estimate of 2250 metres. Four MT lines along the Angore anticline have been interpreted using 2D inversion. A high degree of correlation exists between lineaments observed on an airborne radar image and zones of low resistivity within the high-resistivity material interpreted as Darai limestone. These low-resistivity zones are interpreted as fault zones. Three-dimensional modelling has been used to simulate 3D statics in an otherwise 2D earth. These data were used to test the Groom-Bailey (GB) decomposition for possible benefits in reducing static effects and estimating geoelectric strike in the Papua New Guinea (PNG) field data. It has been found that the GB decomposition can provide improved regional 2D strike estimates in 3D contaminated data. However, in situations such as PNG, where the regional 2D strike is well established and hence can be fixed, the GB decomposition provides apparent resistivities identical to those simply rotated to strike.     

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.     

STUDY ON THE ELECTRICAL STRUCTURE OF THE ANQIU AND JUXIAN ELECTROMAGNETIC STATIONS IN THE TANLU FAULT ZONE

The Yishu fault zone is one of the branch faults of the Tanlu fault zone in its central part. Moderate and strong earthquakes occurred in the Yishu fault zone repeatedly. Due to its complex structure, the Yishu fault zone attracts much attention from earthquake researches. The Anqiu and Juxian electromagnetic stations in Shandong Province locate near the Anqiu-Juxian Fault and Changyi-Dadian Fault, which are branches of the Yishu fault zone, respectively. Geoelectric field and geomagnetic field observation were carried out in these two stations. The Wudi electromagnetic station is in the west of Tanlu fault zone in the Jidong-Bohai block and 230km from Anqiu electromagnetic station. This paper firstly describes the crustal structure near the electromagnetic stations by using magnetotelluric(MT)method. By processing the data carefully, we obtain the MT data in good quality near the stations. The MT data of each electromagnetic station and its nearby area suggests that the electrical structure and geological structure of the station are comparable. This paper applied 1-D and 2-D inversion for MT data and obtained the crustal electrical structure model beneath the Anqiu and Juxian seismic station. The shallow electrical structure from the MT method was compared with the results of symmetrical quadrupole electrical sounding. The model suggests that the electrical structure beneath the Anqiu and Juxian electromagnetic stations is complex and shows the feature of block boundary. The Wudi electromagnetic station is located inside a basin, the crustal structure shows layered feature typical for the stable blocks. Beneath the Anqiu electromagnetic station, there is a 1km-thick relative low resistivity layer in the shallow crust and a high resistivity body beneath it with a depth of 13km. There is a high resistivity structure in the crust beneath the Juxian electromagnetic station. The crustal structures are divided into two different parts by Anqiu-Juxian Fault and Changyi-Dadian Fault, respectively. More conductive layers appear to the west of the two faults. Plenty of fluid possibly exists within the conductive body to the west of Changyi-Dadian Fault, which plays important role in the earthquake generation. There is a relative low resistivity layer in the crust within 1~2km beneath the Wudi electromagnetic station. Beneath the relatively low resistivity layer, a relatively high resistivity layer extends to a depth of around 15km, and the resistivity value decreases with the increase of depth. The electrical resistivity model suggests the seismic activity of the Yishu fault zone around the Anqiu and Juxian electromagnetic stations should be taken into account seriously, and monitoring and research on it need to be strengthened. The results of this paper provide a certain reference value for the crustal structure research to similar stations.     

青藏高原东缘及四川盆地的壳幔导电性结构研究

自从2008年M_S8.0级汶川大地震发生以来,青藏高原东缘便成为地质与地球物理研究的热点区域.该区域的龙门山断裂带标志着青藏高原东缘与四川盆地的边界.汶川地震即发生于龙门山断裂带内的映秀-北川断裂上.该地区现有的研究工作多集中于青藏高原东缘及四川盆地的西部,对四川盆地东部构造情况的研究目前较少.在SinoProbe项目的资助下,完成了一条跨越青藏高原东缘及整个四川盆地的大地电磁测深剖面.该剖面自西北始于青藏高原内部的松潘-甘孜地块,向东南延伸穿过龙门山断裂带、四川盆地内部及四川盆地东部的华蓥山断裂,最终止于重庆东南的川东滑脱褶皱带附近.维性分析表明剖面数据整体二维性较好,通过二维反演得到了最终的电性结构模型.该模型表明,从电性结构上看,沿剖面可分为三个主要的电性结构单元,分别为:浅部高阻、中下地壳低阻的松潘-甘孜地块,浅部低阻、中下地壳相对高阻的四川盆地,以及华蓥山以东整体为高阻特征的扬子克拉通地块.龙门山断裂带在电性结构上表现为倾角较缓、北西倾向的逆冲低阻体,反映了青藏高原东缘相对四川盆地的推覆作用.其在地下向青藏高原内部延伸,深度约为20 km左右.在标志逆冲推覆滑脱面的低阻层下存在一电性梯度带,表征着低阻的青藏高原中下地壳与高阻的扬子地壳之间的电性转换.位于四川盆地东边界的华蓥山断裂在电性结构上表现为一倾向为南东向的低阻体插入高阻的扬子克拉通结晶基底,切割深度约为30 km左右.这一结构反映出华蓥山向西的推覆作用.在电性结构模型的基础上,进一步讨论了青藏高原东缘的壳内物质流、青藏块体与扬子块体的深部关系以及青藏高原东部的隆升机制等构造问题.     

A preparation zone for volcanic explosions beneath Naka-dake crater,Aso volcano,as inferred from magnetotelluric surveys

The 1st crater of Naka-dake, Aso volcano, is one of the most active craters in Japan, and known to have a characteristic cycle of activity that consists of the formation of a crater lake, drying-up of the lake water, and finally a Strombolian-type eruption. Recent observations indicate an increase in eruptive activity including a decrease in the level of the lake water, mud eruptions, and red hot glows on the crater wall. Temporal variations in the geomagnetic field observed around the craters of Naka-dake also indicate that thermal demagnetization of the subsurface rocks has been occurring in shallow subsurface areas around the 1st crater. Volcanic explosions act to release the energy transferred from magma or volcanic fluids. Measurement of the subsurface electrical resistivity is a promising method in investigating the shallow structure of the volcanic edifices, where energy from various sources accumulates, and in investigating the behaviors of magma and volcanic fluids. We carried out audio-frequency magnetotelluric surveys around the craters of Naka-dake in 2004 and 2005 to determine the detailed electrical structure down to a depth of around 1 km. The main objective of this study is to identify the specific subsurface structure that acts to store energy as a preparation zone for volcanic eruption. Two-dimensional inversions were applied to four profiles across the craters, revealing a strongly conductive zone at several hundred meters depth beneath the 1st crater and surrounding area. In contrast, we found no such remarkable conductor at shallow depths beneath the 4th crater, which has been inactive for 70 years, finding instead a relatively resistive body. The distribution of the rotational invariant of the magnetotelluric impedance tensor is consistent with the inversion results. This unusual shallow structure probably reflects the existence of a supply path of high-temperature volcanic gases to the crater bottom. We propose that the upper part of the conductor identified beneath the 1st crater is mainly composed of hydrothermally altered zone that acts both as a cap to upwelling fluids supplied from deep-level magma and as a floor to infiltrating fluid from the crater lake. The relatively resistive body found beneath the 4th crater represents consolidated magma. These results suggest that the shallow conductor beneath the active crater is closely related to a component of the mechanism that controls volcanic activity within Naka-dake.     

大地电磁测深二维反演方法求解复杂电性结构问题的适应性研究

为探讨二维反演方法在三维电性结构中的适应性问题,本研究中设计了一系列的二维/三维正演模型进行实验计算,分析了三维高阻/低阻异常体对模型响应的畸变作用,并从反演模式选择和数据旋转方向两个方面进行模型二维反演的对比分析,与三维反演的结果进行了比较,最后采用了实测数据进一步进行了二维和三维反演的比较实验.实验和研究结果表明,在剖面选择方面,在剖面方向与垂直主构造方向相差不大的情况下,截取剖面方向,将电性主轴旋转到垂直剖面方向的二维反演结果与垂直主构造方向的反演结果都可以较好地还原正演模型,在大的构造的反映上并无太大差异.在地下为二维或近三维条件时,正演模型的主要结构都可以较好地被二维和三维反演解析出来.二维的反演结果可能甚至会比三维的反演结果的边界更清晰,更精确.然而,对于具有较强的三维结构的模型而言,其二维反演结果与原始模型可能仍然存在较大差异,其中TM+TP或TM模式的二维反演结果相对更接近原始模型,而TE模式的结果往往会有较大误差,需要在解释时特别注意以免得出错误结论.     

Rapid evaluation of multifrequency EMI data to characterize buried structures at a historical Jesuit Mission in Argentina

We present results from a geophysical prospection at a historical Jesuit Mission, located in San Ignacio (Misiones-Argentina), declared World Monument some years ago. We studied different sectors looking for buried structures; the total area under study covered 36 ha. In this work we will show the results obtained in a sector that at the time of the prospection was at imminent risk of being damaged by a wrong management of the historical resource. To optimize the data acquisition and the preliminary in-situ analysis of the results, we performed an electromagnetic survey using a multifrequency electromagnetic induction device (GEM-2) to have a first insight of the near surface electrical distribution. From 2D and 3D visualization of data, different targets were identified as possible historical structures. Around these anomalous zones, we performed different dipole–dipole profiles, forming high resolution grids for later 3D inversions. Further inversions of the electric and electromagnetic data completed the characterization of the anomalies. One of the main results of this work is that the 3D electric image obtained from 1D inversion of electromagnetic data coincides with the 3D images obtained from resistivity inversion, but at a much less time consuming cost.     

A regional strategy for geothermal exploration with emphasis on gravity and magnetotellurics

As part of the resource evaluation and exploration program conducted by Los Alamos Scientific Laboratory for the national Hot Dry Rock (HDR) Geothermal Program, a regional magnetotelluric (MT) survey of New Mexico and Arizona is being performed. The MT lines are being located in areas where the results of analysis of residual gravity anomaly maps of Arizona and New Mexico, integrated with other geologic and geophysical studies indicate the greatest potential for HDR resources.The residual gravity anomalies are derived by applying the concept of predicting gravity anomalies from topography. This can be accomplished by employing reductions similar to those used in some isostatic investigations, in which a regional topographic surface is used as the Bouguer reduction datum. The datum is derived by comparison of various harmonics of Bouguer anomalies and elevations of stations. Topography can be used to predict Bouguer anomalies because of isostatic compensation; the resultant anomalies can be considered high frequency residual anomalies or isostatic anomalies corrected for regional compensation. Such maps have been produced for Arizona, New Mexico, west Texas, and Chihuahua, Mexico.The main objective of the MT project is to produce a regional geoelectric contour map of the pervasive deep electrical conductor within the crust and/or upper mantle beneath the Colorado Plateau and the adjacent Basin and Range Province and Rio Grande Rift. The MT survey consists of 200 sites along several long profiles with site spacing of 15–20 km. Pre-existing available MT data are being integrated with the new data. After the data are processed, a one-dimensional inversion is applied to the sounding curve and used as a starting point for 2-D modeling. Such a project and ultimate map will be of major value in studying the regional geophysics and tectonics of the southwest United States as they now apply to HDR resources in particular and geothermal resources in general.Electrical conductivity anomalies of large areal extent are of particular interest in geothermal exploration. Correlation analysis of large conductive anomalies with other geophysical, geological, and geotectonic data is being performed. Preliminary analysis of the data has suggested several major regions of anomalously shallow high electrical conductivity. Among these is the Aquarius area of northwest Arizona which is the site of a longwavelength residual anomaly low, which when modeled and correlated with other geophysical data can be shown to be possibly related to low density and high temperature in the crust at depths of 20 km or less. Preliminary analysis of MT data indicates the possible existence of a mid-crustal high electrical conductivity anomaly in this same region.