青藏高原东边缘冕宁-宜宾剖面电性结构及其构造意义

在青藏高原东边缘沿冕宁—宜宾进行了大地电磁探测研究,剖面西起康滇地轴,向东穿过大凉山地块,终止于四川盆地.利用带地形的NLCG(非线性共轭梯度)方法对资料进行了反演,得到沿剖面的二维电性结构.康滇地轴和大凉山地块地壳中存在向上拱起的高导层(HCL),顶面埋深为10～15 km,最浅处不足10 km,厚度大约15～25 km,最小电阻率小于10 Ωm.四川盆地中下地壳不存在高导层.和该剖面北侧的石棉—乐山剖面的地壳电性结构对比分析表明,高导层在南北方向上可能连续延伸,长度大于100 km.壳内高导层的高导电性与岩石的部分熔融有关,并可能含有百分之几的含盐流体,易于流动和变形.青藏高原东部地壳内的可流动层在向东或东南方向流动过程中,由于受到四川盆地的阻挡,转向南或南南东方向,大体沿着大凉山地块的走向.在东西方向,壳内高导层自川滇地块向东运动,穿过大凉山地块西边界的安宁河断裂和则木河断裂,在大凉山地块东部,向四川盆地深部倾俯.本文对于壳内可流动层的存在及其与青藏高原东边缘的变形和地震活动性的关系进行了探讨.     

有限元直接迭代算法及其在线源频率域电磁响应计算中的应用

采用有限元直接迭代算法实现了线源频率域测深电磁响应的二维正演计算. 首先给出了线源正演问题的有限元直接迭代格式,然后由迭代法进行求解. 在处理奇异源问题上,采用向内递推的组合网格技巧,在源点附近可进行局部加密,并实现粗细网格的对接,从而较好地解决了奇异源附近的计算问题. 还提出一种迭代求取全区视电阻率的方法,避免了远近区的划分. 通过对均匀半空间、层状介质和二维模型电磁响应的计算,获得了与大地电磁测深相似的视电阻率曲线,验证了算法的正确性;通过对计算结果的分析,在理论上说明了线源频率域近区测深的可行性.     

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

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

黑云斜长片麻岩电导率研究

在压力为1.0GPa、温度为400~1 073K条件下,用交流阻抗谱法研究了采自山东泰山的黑云斜长片麻岩平行及垂直面理方向的电导率。在实验温度范围内,黑云斜长片麻岩实验样品电导率的对数值为-6.0~0.5S/m,满足Arrhenius方程。电导率在平行面理方向比垂直面理方向高出约1个数量级。平行面理方向样品电导率分别在第3轮升、降温过程中的881~1 040K之间出现1个明显的电导率突变过程,这可能与黑云母的脱水有关。在低温段,垂直及平行面理方向样品的活化焓分别为0.43和0.49eV,在高温段则分别为3.40和1.53eV。将该实验电导率结果与华东地区大地电磁结果进行对比发现,在中下地壳范围内,该实验结果部分位于大地电磁结果范围内,说明黑云斜长片麻岩可能是组成这一地区中下地壳的候选岩石之一。     

频率域线源近区(过渡区)测深理论初步研究

对频率域无限长线源近区(过渡区)测深的理论基础进行了探索性研究。通过计算均匀半空间、层状介质等多个模型的线源响应,进而定义并计算全区视电阻率。结果发现,利用水平电场、垂直磁场分别计算的全区视电阻率均能较好地反映地下电性结构的变化。把水平磁场与水平电场联合通过二分搜索算法计算的全区视电阻率也能较好地反映地下电性结构的变化。在近区通过单分量定义的全区视电阻率与一维MT曲线吻合得很好,因此可以用成熟的MT反演技术来进行线源近区电磁资料的反演,从而把有源的问题转化到无源的问题。这初步表明频率域线源近区(过渡区)测深是可行的,但还有一些问题需要解决。     

线源时间域电测深全期全区视电阻率求解

由于线源一维的正演国内外研究尚少,文中主要利用水平电场分量求取全区视电阻率。通过计算发现水平电场与电阻率之间是一一对应的关系,这使得定义全区视电阻率具有惟一性和可靠性。文中以水平电场定义视电阻率,然后利用二分搜索算法进行全区视电阻率求解。模型计算表明,该方法定义的全区视电阻率能较准确地反映地下电阻率信息。进一步研究发现:对时间域线源而言,在近区算不好晚期的数据,因为电磁波的传播、衰减很快,同样在极远区算不好早期的数据,而且探测深度是随收发距的增加而增大。最后,加入5%的随机噪声后与理论结果进行对比,得到了收发距与最大探测深度的关系。     

电磁卫星和地震预测

地震前几周、几天、几小时甚至几分钟,电磁场的强度、相位或者谱密度等将发生异常变化,电离层的电子、离子浓度和温度等会发生异常扰动,这些现象的出现向人们发出地震即将来临的信号.这些现象不仅在地面可观测到,在卫星上也能观测到,它们在确定异常现象发生的地点和时间等方面具有各自的优势,把两种观测相结合,对异常信号进行立体监测,是实现地震预测预报的发展趋势.至今,对地面观测已经进行了大量的研究和实践,利用卫星进行观测近几年刚刚开始,但已经显示了它监测地震异常的独特能力,引起各方面的高度重视,成为各国建立立体预警系统不可缺少的组成部分.本文重点介绍利用卫星等观测空间电磁异常现象的研究结果.     

东昆仑断裂带东端和2017年九寨沟7.0级地震区深部电性结构探测

东昆仑断裂带是青藏高原北部一条近EW走向的巨型断裂,其东南尾端发生分叉形成了复杂的马尾状断裂系统,2017年在该区域发生了九寨沟M_S7. 0地震。文中对跨过东昆仑断裂带东端和九寨沟地震区的3条剖面上的大地电磁探测数据进行处理分析,采用三维电磁成像反演技术获取了三维深部电性结构图像。所得结果表明,东昆仑断裂带东端及周边区域内的东昆仑断裂、白龙江断裂和光盖山-迭山断裂表现出向SW倾斜的电性差异带,这些断裂向下延伸并合并于中下地壳的低阻层中,共同组成了由南向北扩展的花状构造。在马尾状构造中,隐伏的虎牙断裂带(北段)在深部表现为明显的低阻边界带;塔藏断裂的规模明显小于虎牙断裂(北段),并与虎牙断裂(北段)组成单侧花状结构;白龙江断裂和光盖山-迭山断裂依然表现为由南向北扩展的花状构造,2组花状结构在深部衔接并统一归并于壳内的低阻层中。2017年九寨沟7. 0级地震的震源区位于高、低阻交界区域,处于松潘-甘孜地块壳内低阻层向NE涌动的端点附近,根据震源区的电性结构和流变结构推测震源深度≤11km。虎牙断裂(北段)的延伸深度和规模大于东侧的塔藏断裂,是2017年九寨沟地震的发震...     

大地电磁资料精细处理和二维反演解释技术研究（四）——阻抗张量分解的多测点-多频点统计成像分析

基于大地电磁阻抗张量分解技术,本文提出了两种电性主轴方位的统计描述图像：随频率变化的统计分布成像（频率分布云图）和随测点序列变化的统计分布成像（测点分布云图）.这两种图像与传统的统计玫瑰图一起,较全面地描述了最佳主轴的分布特征.在进行构造维性分析过程中,通过定义二维有效因子e2d,来压制一维结构和三维结构、突出纯二维结构的影响.e2d被用于电性主轴的统计加权,有效地起到了滤波的作用;同时,统计成像中还考虑了数据质量的影响.为了得到稳定、高质量的区域阻抗张量数据,提出并实现了共主轴的多测点-多频点阻抗张量分解新算法.最终,完成了以上各项处理手段的可视化实现.本文通过两个理论模型和一个实测算例,以共轭阻抗法（CCZ法）为基础,展示了这一新技术的有效性.     

汶川MS8.0级地震成因三“层次”分析——基于深部电性结构

2008年5月12日汶川M_S8级地震的发生不是局部地区孤立的构造事件,研究汶川地震的孕震机制,应该把局部分析和区域分析相结合,关注地壳上地幔直至地幔过渡带的深部结构.基于近年来在东北、华北和汶川地震附近地区进行的深部结构电磁探测结果,结合地震学等其他资料,从太平洋板块的俯冲、印度板块的碰撞和松潘甘孜地块的推挤三个“层次”探讨分析汶川特大地震的成因.太平洋板块向亚洲大陆的俯冲作用,导致中国大陆东部地幔过渡带深度较普遍地存在着停滞的板片,它对汶川地震的影响不可忽视.印度板块与青藏高原的碰撞,使组成高原的各地块发生向北和向东的运动,各地块向东的运动作用于南北地震带中南段,影响到该区域的地震活动.松潘甘孜地块向四川地块的推挤,使松潘甘孜地块运动方向和龙门山断裂带形成“丁”字形结构,龙门山断裂带显示为较陡直的电性边界,加剧了汶川地震前的应力积累,可能是汶川地震发生的最直接的诱因.