鄂尔多斯盆地西缘构造带北段深部电性结构

在横跨鄂尔多斯盆地西缘构造带北段的查甘池—银川—五湖洞约200 km长的东西向剖面上,进行了67个测点的大地电磁探测.使用“远参考道”和Robust技术处理数据.分析了各测点视电阻率、阻抗相位、二维偏离度、电性主轴方位角、磁实感应矢量等参数,采用NLCG二维反演方法对TE和TM两种模式的数据进行了二维反演.得到的二维电性结构表明,沿剖面查汗断裂带、贺兰山东缘断裂带和黄河断裂带是明显较大型电性边界,为超壳断裂带,而三关口断裂带深部延深不大.沿剖面阿拉善地块、贺兰山褶皱带、银川断陷盆地和鄂尔多斯地块具有明显不同的深部电性结构特征.阿拉善地块内部除浅表电阻率较低外,以下到深度约50 km都表现为高电阻特性.贺兰山褶皱带电性结构复杂,电阻率高低相间.银川盆地具有上宽下窄最深达约8 km低阻层,具有断陷盆地特征.鄂尔多斯地块具有低-高-低的深部电性结构特征,成层性较明显.     

青藏高原东北缘海原构造带马东山阶区深部电性结构特征及其构造意义

海原一六盘山构造带是青藏高原东北缘地区的一条重要边界,在海原断裂带和六盘山断裂带接触区形成了特殊的马东山挤压阶区,本文对跨过该挤压阶区一条密集测点大地电磁剖面数据进行了处理和二维反演,获得的深部电性结构图像揭示在马东山挤压阶区深部电性结构表现为在高阻背景下镶嵌多个向西南倾斜的低阻条带电阻率结构样式,并在深度约25 km汇聚到中下地壳低阻层内,共同组成"正花状"结构;海原一六盘山构造带西南侧到陇中盆地区间呈现高、低阻相互"楔合"的深部结构特征,而其东北侧的鄂尔多斯西缘带自地表到中下地壳为较完整的高阻块体.另外结合跨过海原断裂带中段和西秦岭造山带的大地电磁探测结果,对海原一六盘山构造带分段性及其两侧的陇中盆地和鄂尔多斯地块的接触关系进行了研究分析.大地电磁探测成果佐证了在海原断裂带中段为具有走滑特点的断裂,而其尾端与六盘山断裂带斜交区域的马东山地区发生了强烈的逆冲推覆与褶皱变形;活动构造研究发现沿海原断裂带所产生的左旋走滑位移被其尾端的马东山、六盘山以东西向的地壳缩短调节吸收,GPS观测表明青藏高原东北缘地区现今构造变形分布在海原一六盘山构造带以西上百公里的范围内,陇中盆地一海原一六盘山构造带和鄂尔多斯地块一线的深部电性结构图像也很好地解释了该区变形状态:海原一六盘山构造带带及西南盘的陇中盆地的中下地壳非常破碎,在青藏高原向北东方向的推挤下容易发生变形,而北东盘鄂尔多斯地块地壳结构完整,很难发生构造变形.对海原一六盘山构造带马东山阶区和龙门山构造带的深部电性结构及变形特征等进行了比较分析,发现该区有与2008年汶川地震相似的深部构造背景,应重视该区强震孕育环境的探测研究.     

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

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

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.     

Geophysical Images of the North Anatolian Fault Zone in the Erzincan Basin, Eastern Turkey, and their Tectonic Implications

The collision between the Arabian and Eurasian plates in eastern Turkey causes the Anatolian block to move westward. The North Anatolian Fault (NAF) is a major strike-slip fault that forms the northern boundary of the Anatolian block, and the Erzincan Basin is the largest sedimentary basin on the NAF. In the last century, two large earthquakes have ruptured the NAF within the Erzincan Basin and caused major damage (M _s = 8.0 in 1939 and M _s = 6.8 in 1992). The seismic hazard in Erzincan from future earthquakes on the NAF is significant because the unconsolidated sedimentary basin can amplify the ground motion during an earthquake. The amount of amplification depends on the thickness and geometry of the basin. Geophysical constraints can be used to image basin depth and predict the amount of seismic amplification. In this study, the basin geometry and fault zone structure were investigated using broadband magnetotelluric (MT) data collected on two profiles crossing the Erzincan Basin. A total of 24 broadband MT stations were acquired with 1–2 km spacing in 2005. Inversion of the MT data with 1D, 2D and 3D algorithms showed that the maximum thickness of the unconsolidated sediments is ~3 km in the Erzincan Basin. The MT resistivity models show that the northern flanks of the basin have a steeper dip than the southern flanks, and the basin deepens towards the east where it has a depth of 3.5 km. The MT models also show that the structure of the NAF may vary from east to west along the Erzincan Basin.     

Deep resistivity structure of the Dikili-Bergama region, west Anatolia, revealed by two- dimensional inversion of vertical electrical sounding data

Within the framework of the National Marine Geological and Geophysical Program, we re‐examined deep vertical electrical sounding (VES) data. The data, measured in 1968 by the General Directorate of Mineral Research and Exploration (MTA) of Turkey with the aim of exploring the deep resistivity structure of the Dikili–Bergama region, focus on the geothermal potential. The geoelectrical resistivity survey was conducted using a Schlumberger array with a maximum electrode half‐spacing of 4.5 km. The two‐dimensional (2D) inversion was utilized to interpret the VES data that were collected along 15‐ to 30‐km profiles. The 2D resistivity–depth cross‐sections obtained show very low resistivity values near the Dikili and Kaynarca hot springs. The 2D inversion results also indicate the presence of fault zones striking nearly N–S and E–W, and fault‐bounded graben‐horst structures that show promising potential for geothermal field resources. The 2D gravity model, which is in good agreement with the density variation of the region, supports the resistivity structure revealed by 2D inversion. The lithology information obtained from the borehole near Kaynarca also confirms the results of the resistivity interpretation and the density model.     

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

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

滇西牟定—香格里拉电性剖面及深部构造

通过对滇西NW—SE向牟定—香格里拉大地电磁测深剖面的反演并结合地质构造、小震精定位资料,综合分析了剖面经过地区的深部构造.在滇西北地区深部发育一规模巨大的近水平产出地壳高导层,此高导层从香格里拉一直延伸到永胜盆地以东,沿剖面水平延伸超过200km.地壳高导层又可分为两段,两段高导层在剖面上形成两个香蕉型相连的复合形态,高导层之上的高阻体则形成两个碗型相连的复合结构,碗型内部有向东倾斜的相对低阻带存在.滇西北地区整体构造格架可以解译为一套两个主推覆面构成的叠瓦式推覆系统,两个主推覆面高低起伏形成总长度超过250km的底部滑脱带.第一个主推覆面的长度超过150km,并在丽江东部上翘到接近地表,小金河断裂是其延伸到地表的主要破裂带.第二个主推覆面长度超过120km,过程海断裂后上翘并在永胜以东地区出露地表.根据反向低阻带影像,可以推断推覆体反冲构造发育,在玉龙山至丽江以东和永胜盆地西缘至平川盆地东缘形成两处冲起构造.推覆面在剖面上呈现出平缓-陡峭-再平缓-再陡峭的断坪和断坡相间的阶梯状结构.     

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.     

南海地区地震背景噪声成像与壳幔深部结构

利用南海地区28个陆地地震台站和2个布设于太平岛和东沙岛的新增海岛地震台站2011—2016年间的连续地震背景噪声波形数据,使用互相关方法计算得到了台站间的互相关函数,并提取出Rayleigh面波群速度和相速度频散曲线.采用快速行进和子空间方法反演获得了南海及周边地区12～40s周期的Rayleigh面波群速度和相速度图像,并联合反演得到了研究区深至60km的三维S波速度结构.考虑到南海数千米厚海水层对于面波频散反演的严重影响,本文在反演模型中加入了水层,显著提高了反演结果的可靠性.成像结果表明:南海及周边地区地壳上地幔顶部S波速度结构存在显著的横向不均匀性,并与这一区域的主要构造单元具有较好的空间对应关系.在5～10km深度,莺歌海—宋红盆地区的低速异常特征可能与盆地较厚的沉积层有关.在5～15km深度,海域高速异常区与海盆空间位置具有高度一致性,推测与海盆区地壳厚度相对陆缘区明显偏薄有关.当深度从20km增加至30km,海盆区的高速特征扩展至了陆缘地区,反映了地壳厚度从海盆至陆缘逐渐增厚的趋势,与OBS(海底地震仪)深地震剖面给出的地壳精细结构结果一致.至35～60km深度,海盆的高速异常特征依然明显,且速度值随深度增加整体呈现上升的趋势,推测南海海盆区的岩石圈厚度应该大于60km.     

Joint inversion of long‐offset and central‐loop transient electromagnetic data: Application to a mud volcano exploration in Perekishkul,Azerbaijan

Mud volcanism is commonly observed in Azerbaijan and the surrounding South Caspian Basin. This natural phenomenon is very similar to magmatic volcanoes but differs in one considerable aspect: Magmatic volcanoes are generally the result of ascending molten rock within the Earth's crust, whereas mud volcanoes are characterised by expelling mixtures of water, mud, and gas. The majority of mud volcanoes have been observed on ocean floors or in deep sedimentary basins, such as those found in Azerbaijan. Furthermore, their occurrences in Azerbaijan are generally closely associated with hydrocarbon reservoirs and are therefore of immense economic and geological interest. The broadside long‐offset transient electromagnetic method and the central‐loop transient electromagnetic method were applied to study the inner structure of such mud volcanoes and to determine the depth of a resistive geological formation that is predicted to contain the majority of the hydrocarbon reservoirs in the survey area. One‐dimensional joint inversion of central‐loop and long‐offset transient electromagnetic data was performed using the inversion schemes of Occam and Marquardt. By using the joint inversion models, a subsurface resistivity structure ranging from the surface to a depth of approximately 7 km was determined. Along a profile running perpendicular to the assumed strike direction, lateral resistivity variations could only be determined in the shallow depth range using the transient electromagnetic data. An attempt to resolve further two‐dimensional/three‐dimensional resistivity structures, representing possible mud migration paths at large depths using the long‐offset transient electromagnetic data, failed. Moreover, the joint inversion models led to ambiguous results regarding the depth and resistivity of the hydrocarbon target formation due to poor resolution at great depths (>5 km). Thus, 1D/2D modelling studies were subsequently performed to investigate the influence of the resistive terminating half‐space on the measured long‐offset transient electromagnetic data. The 1D joint inversion models were utilised as starting models for both the 1D and 2D modelling studies. The results tend to show that a resistive terminating half‐space, implying the presence of the target formation, is the favourable geological setting. Furthermore, the 2D modelling study aimed to fit all measured long‐offset transient electromagnetic Ex transients along the profile simultaneously. Consequently, 3125 2D forward calculations were necessary to determine the best‐fit resistivity model. The results are consistent with the 1D inversion, indicating that the data are best described by a resistive terminating half‐space, although the resistivity and depth cannot be determined clearly.     

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

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

利用密集台阵近震层析成像研究云南宾川上地壳速度结构

滇西北宾川盆地是发育于红河断裂和程海断裂交汇处的晚新生代张扭断陷盆地,该区活动断裂发育且历史地震比较活跃.对宾川盆地及邻区进行高精度浅层和上地壳精细结构研究,有助于深入认识该区主要发震构造的深浅部特征.基于2017年在宾川盆地及其附近开展的为期2个月的密集台阵观测数据,我们对该区96个小震共拾取了117221条初至P波和5475条初至S波震相,并利用simul2000开展了地震重定位和体波层析成像研究.结果表明：（1）小震活动主要集中在宾川盆地东缘断裂的弧形转折部位,并在洱海南侧呈现沿北东向断裂的条带状分布现象,反映了区域上近南北向至北东向断裂是主要控震构造,其次是北西向断裂带.（2）0 km的速度分布与区域地形有很好的对应关系.山地呈现高速异常,宾川盆地呈现低速异常.从3 km至9 km,高低速分界与断层有很好的对应,并且越往深部,近南北向至北北东向的宾川盆地东缘断裂在上地壳构造的控制作用越明显.（3）上地壳层析成像结果同时揭示了宾川盆地东缘断裂的三维形态变化在空间上呈现出南北部倾角大、中部倾角缓的变化特征,可能与区域地块的旋转变形过程有关.（4）综合高精度浅层速度结构和地震重定位结果可知,区域上的近南北向至北北东向断裂正逐步取代北西向构造,成为主要的区域分界断裂和控震构造.新的研究结果为深入理解该区的主要控震构造及其深部结构特征提供了重要依据.     

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

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

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

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

西秦岭与南北地震构造带交汇区深部电性结构特征

西秦岭造山带与南北地震构造带接触区是中国大陆最重要的南北向和东西向构造转化的接合部位之一.本文介绍了分别位于该区106°E东、西两侧的LMS-L3和DBS-L1两条大地电磁剖面的探测结果,两条剖面分别跨过了龙门山构造带东北段的青川段和宁强段.采用大地电磁相位张量分解技术对两条剖面上各测点的电性走向、二维偏离度等进行了计算和分析,采用NLCG二维反演方法对TE+TM模式的视电阻率和阻抗相位数据进行了二维联合反演.反演得到二维电性结构,在经度106°西侧LMS-L3剖面的深部电性结构自北向南揭示出,西秦岭北缘、成县盆地北缘、康县（即勉略构造带）和平武-青川断裂带都表现为明显的电性梯度带,深部延伸可达几十公里;西秦岭造山带、碧口地块与龙门山构造带东北段3个构造单元整体表现为高电阻体、呈现往南叠合且角度逐渐变陡的趋势.在106°E西侧西秦岭造山带区域的深部存在壳内低阻层,而东侧区域表现为高电阻体,深部电性结构在106°E东、西两侧的差异与该区深部速度结构特征一致,东、西两侧深部结构差异可能是该区中强地震分布差异的深层原因.LMS-L3和DBS-L1两条剖面南段的深部电性结构图像揭示出龙门山构造带东北部的青川段和宁强段内的平武-青川断裂带具有明显不同的深部结构特征,平武-青川断裂带在青川段为明显的电性梯度带,在宁强段不再表现为电性梯度带,而是完整的高电阻块体.汶川强余震向东北发展止于青川青木川附近,与平武-青川断裂带延伸深度和向北东方向的延伸长度密切相关,同时高电阻块体的宁强段对汶川强余震东北发展起到了阻挡作用.     

利用大地电磁技术揭示2016年1月21日青海门源MS6.4地震隐伏地震构造和孕震环境

2016年1月21日01时13分在青海省海北州门源县发生了M_S6.4地震,震中位置位于青藏高原东北缘地区祁连造山带内的祁连—海原断裂带冷龙岭断裂部分附近,震源深度约11.4 km,震源机制解显示该次地震为一次纯逆冲型地震.我们于2015年7—8月期间完成了跨过祁连造山带紧邻穿过2016年1月21日青海门源M_S6.4地震震中区的大地电磁探测剖面（DKLB-M）和古浪地震大地电磁加密测量剖面（HYFP）.本文对所采集到的数据进行了先进的数据处理和反演工作,获得了二维电性结构图.结合青藏高原东北缘地区最新获得的相对于欧亚板块2009—2015年GPS速度场分布特征,1月21日门源M_S6.4地震主震与余震分布特征以及其他地质与地球物理资料等,探讨了门源M_S6.4地震的发震断裂,断裂带空间展布、延伸位置,分析了门源M_S6.4地震孕震环境与地震动力学背景等以及祁连山地区深部构造特征等相关问题.所获结论如下：2016年门源M_S6.4地震震源区下存在较宽的SW向低阻体,推测冷龙岭断裂下方可能形成了明显的力学强度软弱区,这种力学强度软弱区的存在反映了介质的力学性质并促进了地震蠕动、滑移和发生;冷龙岭北侧断裂可能对门源M_S6.4地震主震和余震的发生起控制作用,而该断裂为冷龙岭断裂在青藏高原北东向拓展过程中产生的伴生断裂,表现出逆冲特征;现今水准场、重力场、GPS速度场分布特征以及大地电磁探测结果均表明祁连—海原断裂带冷龙岭断裂部分为青藏高原东北缘地区最为明显的一条边界断裂,受控于青藏高原北东向拓展和阿拉善地块的阻挡作用,冷龙岭断裂附近目前正处于青藏高原北东向拓展作用最强烈、构造转化最剧烈的地区,这种动力学环境可能是门源M_S6.4地震发生的最主要原因,与1927年古浪M_S8.0地震和1954年民勤M_S7.0地震相似,2016年门源M_S6.4地震的发生同样是青藏高原北东向拓展过程中的一次地震事件.     

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.     

华北地区上部地壳结构的三维有限差分层析成像

采用正方形网格对华北地区速度模型进行参数化、网格维数为649×90×16,网格间距为1km. 通过对华北地区30条深地震测深剖面的Pg波三维有限差分速度层析成像反演,获得了该地区上部地壳的三维速度结构. 反演结果显示,低速的"棱形"华北盆地被高速的隆起所包围;而华北盆地又由两个大的断陷带中间夹一个隆起所组成.两个大的断陷带又被分割为多个小的断陷盆地.研究区内的断层多为铲式正断层,结晶基底深度变化较大. 所有这些表明华北地区具有复杂的地质构造和发生强震的深部背景.     

利用远震接收函数探测四川盆地及周边地区的地壳结构

青藏高原东南缘的龙门山断裂两侧具有陡峭的地形特征,在约50~100 km的水平距离内,地形高程从2000 m增加到4000 m,该区强烈的壳幔变形特征及地球动力学模式一直是研究的热点问题.本文从四川地区49个固定台站记录的远震资料提取了P波接收函数,获得了四川盆地及周边的地壳厚度和泊松比,并以此构建反演的初始模型.在线性反演的基础上,引入了分别拟合低频和高频接收函数的两步反演技术,用以反演台站下方的地壳S波速度结构.数字试验表明,该方法可以有效抑制接收函数反演的不唯一性,为了得到最优解,最后用Bootstrap重采样技术估计解的不确定性.结果表明,四川盆地的地壳厚度在40~46 km,松潘-甘孜块体北部的地壳厚度为46~52 km,而南部增厚到50~60 km.从四川盆地向西跨过龙门山断裂,地壳厚度增加了10~15 km.在四川盆地及周边地区,地壳泊松比在0.26~0.32之间,呈块体分布特征,高泊松比（0.28~0.32）主要沿龙门山断裂以及安宁河-小江断裂分布.地壳S波速度结构表明,来自青藏高原中部的中下地壳低速层可能受到了坚硬的四川盆地阻挡,改变原来的运动方向并沿龙门山断裂展布,由于低速层的囤积导致该区地形陡峭和下地壳增厚.