大地电磁测控深二维地形影响及其校正方法研究

采用有限元法实现ＭＴ二维模型的正演模拟,将地形条件下的二维地电模型的ＭＴ响应分解为地形引起的畸变场和与测点高程位置有关的稳定场两部进行分析。用有限元法模拟二维地形条件下地电结构的ＭＴ响应和地形的ＭＴ响应,采用幽会法消除地形引起的干扰场,然后,将消除地形影响后的视电阻率值通过延拓方法获水平地形条件下二维地电结构的视电阻率分布,完善在地电磁测深视电阻率的地形校正。     

青藏高原东北缘祁连山造山带至阿拉善地块壳幔电性结构研究

在青藏高原东北缘祁连山造山带至阿拉善地块之间完成了一条372km的大地电磁剖面,通过二维反演计算,获得了沿剖面180km深的壳幔电性结构模型,结合研究区地质和地球物理资料开展综合分析,研究结果表明:(1)剖面自南向北所经过的祁连山造山带、走廊过渡带和阿拉善地块对应3种壳幔电性结构模型:东祁连壳幔高-低-高阻似层状电性结构、河西走廊壳幔低阻带状电性结构和阿拉善南缘壳幔高-低-高阻层状电性结构.(2)剖面所经过的主要断裂带在电性结构上表现为低阻异常带或电性梯度带,并且止于中上地壳或消失于下地壳低阻层中.除这些分布于中上地壳的断裂系统以外,在下地壳至上地幔顶部还存在两条切割莫霍面的壳幔韧性剪切带:西华山北缘壳幔韧性剪切带和阿拉善南缘壳幔韧性剪切带.其中,西华山北缘壳幔韧性剪切带可能是1920年海原8.6级地震发生的深部背景之一;而阿拉善南缘壳幔剪切带可能是卫宁北山燕山晚期和喜山期幔源岩浆上升到地壳浅部或喷出到地表的通道,为在该区域寻找晚中生代至新生代含矿隐伏岩体提供了深部电性结构依据.(3)由若干形状不规则、彼此不相连的"碎块状"极高阻块体组成的中上地壳与"似层状"的中下地壳低阻层共同构成的地壳电性结构,是引起青藏高原东北缘强烈破坏性地震最佳的地壳电性结构组合之一.印度板块向欧亚板块俯冲碰撞楔入引起青藏高原块体向北东方向运移与阿拉善地块向南的俯冲碰撞楔入,是青藏高原东北缘强震活动带产生的动力学背景.     

地震放大效应与地下地质结构——龙门山山前玉皇观区域观测数据分析

通过布置于龙门山断裂带中段、龙门山山前玉皇观区域的地震观测台站阵列接收地震数据,研究该区域的地震动放大效应和地下地质结构.观测阵列共10台宽频带地震仪,分布在玉皇观河口冲积扇区域.分别采用参考场址谱比法(RSSR)和HV谱比法(HVSR)计算64个高信噪比近震数据的振幅谱比函数,结果显示在玉皇观区域具有较明显的地震动放大效应,并且局部场址效应显著.以S06场址为例,建立近地表地震地质模型,通过SH波放大效应正演模拟研究该场址的地震动放大模式.RSSR与HVSR的结果表明,两者所计算的场址放大效应主频一致,但是HVSR的放大峰值却比RSSR的放大峰值大一倍左右,表明HVSR的结果可能包含了波场在近地表低速层之下传播路径的改造作用.另外,采用27个远震P波的接收函数计算了该区域地壳上地幔S波速度结构.接收函数研究结果显示玉皇观地区的莫霍面深度为44 km,沉积盖层、结晶地壳和上地幔的S波速度分别为2.5 km·s~(-1)、3.5 km·s~(-1)和4.5 km·s~(-1).观测阵列台站之间的接收函数反演结果一致性较好,说明本研究区域范围内地形地貌等近地表结构因素的相对变化对接收函数的影响不大.     

龙门山构造带深部动力学过程与地表地质过程的耦合关系

发生在地球浅层的2008年汶川地震驱动了龙门山及前陆地区的地表同震垂向位移.根据冲断带-前陆盆地弹性挠曲模型理论,在进行弹性挠曲模拟反演的基础上,结合对深部地球物理特征(泊松比、电性结构)的分析,发现龙门山前陆盆地现今岩石圈有效弹性厚度(T_e)具有自东向西逐渐减薄的趋势,自川中地区的30~40 km减至龙门山地区的10~20 km.在对晚三叠世以来前陆盆地各阶段盆地结构进行刻画的基础上,进行弹性挠曲模拟反演,推断龙门山前陆盆地的前渊地区(四川盆地西部)岩石圈的T_e值自晚三叠世以来具有逐渐减薄的趋势.这可能与松潘一甘孜地块下方广泛存在的软流圈热物质对四川盆地西部岩石圈下部的长期加热而导致的熔融有关,反映了地球深部动力学过程与地球表层盆地演化之间的耦合关系.     

基于VTI各向异性介质的频率域海洋可控源电磁三维约束反演

近年来,海洋可控源电磁法(MCSEM)被引入油气勘探领域以降低勘探风险.在海洋环境中,受沉积因素所造成的电阻率各向异性的影响,地电模型往往会非常复杂.为更好地反映地下电性结构,本文实现了基于VTI各向异性介质的频率域海洋可控源电磁三维反演.其中,正演采用基于Yee氏交错网格的三维有限差分算法,所形成的离散线性系统通过大规模并行矩阵直接求解器(MUMPS)进行求解.反演采用基于不等式约束的有限内存BFGS(L-BFGS)算法.最后,利用VTI各向异性介质合成数据,分别进行了电阻率各向异性覆盖层和电阻率各向异性高阻层的三维反演,结果表明:(1)基于并行直接法的MCSEM非常适用于海洋电磁所特有的多场源问题;(2)针对各向异性覆盖层模型进行三维各向异性约束反演,提高了解的可靠性;(3)针对电阻率各向异性高阻层,Inline和broadside数据覆盖的反演结果对异常体位置有很好的反映.     

龙门山构造带壳幔电性结构特征及其与汶川、芦山强震关系

青藏高原东缘龙门山构造带是研究青藏高原地壳物质向东侧向挤出的焦点地区.为探索龙门山构造带活动构造特征及其与发震构造的关系,本文通过布置垂直龙门山构造带南段芦山地震震源区的大地电磁测深剖面,运用多种数据处理手段,得到研究区可靠的电性结构,并通过与已有龙门山中段和北段剖面进行对比分析.研究表明：（1）青藏高原东缘岩石圈存在明显的低阻异常带--松潘岩石圈低阻带,该低阻异常带沿龙日坝断裂-岷山断裂-龙门山后山断裂分布,形成松潘-甘孜地块向扬子地块俯冲的深部动力学模式,通过统计研究区的历史强震,发现震源主要沿低阻异常带东侧分布,同时,低阻异常带也是低速度、低密度异常带,松潘岩石圈低阻带可能是扬子地块的西缘边界;（2）青藏高原物质东移过程中,受到克拉通型四川盆地的强烈阻挡,龙门山构造带表层岩块和物质发生仰冲推覆,表现为逆冲推覆特征的薄皮构造,中下地壳和上地幔顶部物质向龙门山构造带岩石圈深部俯冲,印支运动晚期,扬子古板块持续向华北板块俯冲,在上述构造运动作用下,呈现出刚性的上扬子地块西缘高阻楔形体向西插入柔性青藏块体的楔状构造;（3）根据电性结构推断,芦山地震受到深部上里隐伏壳幔韧性剪切带向上扩展的影响,构成芦山地震的深部主要动力来源;汶川地震的发生,在龙门山南段形成应力加载区,是触发或加快芦山地震孕育发生的另一个动力来源.     

特殊区域旋翼无人机航磁测量研究

为解决地面人员难以到达的特殊区域大比例尺磁测数据获取问题,本文提出一种基于三轴磁通门传感器,采用旋翼无人机进行航空磁测的方法.项目组研制适合旋翼无人机搭载的航磁数据记录仪,并完成旋翼无人机、三轴磁通门传感器、航磁数据记录仪的集成工作,开展了三轴磁通门传感器安装点位噪声水平试验、地面空中对比试验、磁性体识别试验、动态飞行一致性试验、区域飞行试验.旋翼无人机航磁系统具有操作简单、设备轻便、效率高的特点,能在地面人员难以达到的区域开展低空航磁测量.中国四川省黑竹沟地区是地磁异常区域,项目组在该区域内的沟壑、峡谷、湖泊等实地开展了野外应用测试,完成了这些特殊区域的航磁调查工作.