自动构建大地电磁二维反演的测点中心网格

在大地电磁二维反演中,如果网格不均匀或者测点在网格中的位置不恰当,则将出现测点偏离现象,影响反演的精度.此外,由于正演精度和稳定性的需要,反演中网格间距不能差别太大,需要具有一定的光滑性.本文提出了几种自动构建测点中心网格的方法来消除测点偏离现象,并以正则化约束条件来保证网格的光滑程度.本文还提出了一套网格质量检验指标.算例对比研究表明,对于不同的测点集,最佳网格构建的方法可能不相同,但其中“最平模型约束动态配点法”、“最光模型约束动态配点法”拥有最强的适应能力,它们在各种情况下所构建的网格一般同时具备较光滑的网格间距变化和较小的测点偏离程度.此外,采用高质量（光滑性好、测点偏离程度低）网格在一定程度上不仅可提高反演结果的可靠性,还能提高反演的计算效率.本文的方法为获得高质量的大地电磁二维反演网格提供了简单有效的途径.     

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

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

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

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

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

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

青藏高原东缘川滇构造区深部电性结构特征

本文对位于青藏高原东缘川滇构造区的贡山一绥江大地电磁测深(MT)剖面数据进行反演,获得沿剖面的深部电性结构,为研究喜马拉雅东构造结、川滇菱形地块与华南地块的构造变形特征、壳幔耦合关系、地块间接触关系以及相互作用等问题,提供电性结构的依据.研究发现:(1)电性结构揭示澜沧江断裂带和小金河断裂带为深大断裂带,控制着研究区的深部结构特征和形变机制;(2)澜沧江断裂带和金沙江断裂带之间的高阻体,可能是扬子古地块的残留部分;小金河断裂带和安宁河断裂带之间的高阻体,则是峨眉山大火山省喷发形成的冕宁一越西杂岩带;(3)在滇西地块、川滇地块和大凉山地块均存在低阻层,它们的介质属性有所不同,滇西地块下的低阻层"疑似"高热状态的岩浆囊,主要由缅甸弧向东俯冲运动引起的,中上地壳的高热状态使地块的活动性增强;川滇地块内部的壳内低阻层的成因为:理塘断裂带和小金河断裂带之间的地表低阻层由破碎带充水所致,而金沙江断裂带和理塘断裂带之间的中地壳低阻层可能是由局部熔融物质或含盐流体导致的,其为壳内物质运移的通道.从而在地下物质发生大规模走滑运动的过程中起到引导作用;川滇地块东部和大凉山地块西部的壳内低阻层可能与地慢物质的上涌有关;马边断裂带附近的低阻体可能与破碎带变宽和破碎带内的流体有关.     

误差对大地电磁测深反演的影响

本文讨论了几种不同分布的误差对大地电磁测深反演的影响.首先基于给定的正演模型计算得到相应的理论合成数据,然后在合成数据中分别加入满足不同均值和方差的高斯分布、均匀分布和Gamma分布的误差,最后使用非线性共轭梯度法进行反演.反演结果表明, 较简单的模型在考虑的几种误差分布下可以较好地给出模型的基本信息;但复杂的模型随着误差方差的增大,与已知模型的偏差越明显; 在均值和方差相同时,上述三种不同分布的误差的反演结果相差很小.     

阿尔山火山区地壳上地幔电性结构初探

在阿尔山活火山区一条北北西向测线上进行了7个大地电磁测深点的观测。2-D解释结果表明,研究区内新、老两条火山带可能存在通往深部的岩浆通道。其中,新发现的活火山带地下在10~12km还保持着较高的热状态,很可能富含流体,在30~50km处可能是地幔向上的供热通道;而老火山带30km以上可能存在一个正在冷却的岩浆通道,两条火山条带的深部可能是同源的。     

大别山超高压变质带深部电性结构及其动力学意义初步研究

通过在大别造山带东部横穿超高压变质带的一条NNE向剖面大地电磁测深资料的分析解释,获得了关于沿剖面的地壳上地幔二维电性结构,显示北淮阳与大别地块是电性差异显著的构造单元,它们之间的界面与晓天—磨子潭断裂对应;晓天—磨子潭断裂倾向北,在中上地壳层位出现错动解耦现象;从地表向深处可划分出4个主要电性层：地表风化层、中上地壳高阻层、壳内相对高导层以及上地幔层;大别地块内中、上地壳层位以高阻层为主,与高压-超高压变质岩分布区对应,高阻层最厚处在岳西—英山之间;在大别地块内,推测存在燕山期花岗质岩浆活动的通道,它们造成了超高压变质岩的进一步抬升,同时影响了大别地块内存在的壳内相对高导层的分布,壳内相对高导层在层位上相差较大.     

三维大地电磁自适应L1范数正则化反演

常规三维大地电磁反演的正则项为L₂范数,它以电阻率空间分布函数处处光滑为模型期望,弱化了算法对电性突变界面的分辨能力.本文实现了正则项为L₁范数的三维大地电磁反演算法,让模型空间梯度向量更有机会取得稀疏解,在充分正则的迭代下能够有效突出模型真实电性界面.为避免L₁范数零点不可导带来的求解困难,使用迭代重加权最小二乘法把原问题转换为一系列L₂正则子问题迭代求解.每个子问题的极小方法使用改进型拟牛顿法,其下降方向既能保证正则项海塞矩阵的精确性,又能允许反演过程随迭代灵活更新正则因子.使用比值法或分段衰减法自适应更新正则因子以避免迭代早期陷入奇异解,从而提升反演收敛的稳定性并降低初始模型依赖度.合成的无噪数据反演表明L₁正则算法的模型恢复效果优于L₂正则;不同噪声水平的合成数据反演表明本文的算法具有稳健性;实测数据反演对比表明在合理的正则因子调整策略下,L₁正则反演结果的模型分辨率优于L₂正则.另外,不同初始模型的反演...     

利用大地电磁密集阵列研究蚌埠地震构造环境

本文基于蚌埠及邻区100个高质量测点组成的宽频带大地电磁密集阵列数据,首先利用相位张量和阻抗张量统计成像技术定性分析了研究区构造维性和电性结构特征,然后采用多重网格法结合多次迭代重构的三维反演技术获得了其精细可靠的三维电阻率结构.结果显示蚌埠隆起南、北部电性结构存在显著差异,支持其南部和北部存在不同的构造演化和变形机制的观点.蚌埠隆起内部存在多个切割高阻异常体、与断裂带对应的低阻条带,其中下地壳存在高导异常.综合分析推断,该中下地壳高导异常为流体渗入到先存多条交切断层的破碎带所致,深部物质对蚌埠隆起区构造变形的作用不显著.本文研究认为蚌埠及其邻近地区具备强震孕育的深部结构环境,但该区域构造变形驱动机制有待探讨.