电阻率横向剖面法在倾斜断层精确探测中的应用研究

在横向剖面法研究的基础上,选取30度和60度倾斜良导体模型,利用有限元三维数值模拟方法,对比研究横向剖面法和常规剖面法的异常特征.研究结果表明：两方法对断层的倾向都有较好反映,对倾角难以判定.为研究倾角定量化判定,本文提出了横向剖面电阻率成像方法与常规的电阻率层析成像方法对比研究,结果表明：两者均较准确的反映出断层的倾角变化,前者效果略优于后者.     

华北地区地壳P波和S波速度结构的双差层析成像

基于华北地区(37°N—42°N,113.5°E—118.5°E)133个固定地震台站收集到的P波和S波震相数据,利用双差层析成像法反演了该地区地壳三维速度结构并对所用地震进行了重定位.结果显示:地震走时残差均方根的平均值由重定位前的0.265 s下降至0.008 s;重定位后的震源主要分布于6—16 km深度范围内;...     

阿尔金断裂带地壳和上地幔结构的P波层析成像

阿尔金断裂是青藏高原西北边缘最主要的断裂．天然地震P波层析成像结果揭示了阿尔金断裂（严格地说这里指阿尔金断裂中部,下同）为一条宽约40km左右的低速带,并以比较直立的产状向下延伸至150km左右深度．结果同时显示塔里木岩石圈曾经挤入到柴达木盆地下面,并受到后来发育起来的阿尔金断裂的切割．本次研究结果支持青藏高原中北部上地幔热物质参与了支撑高原高海拔地形的重力均衡作用的假设。     

中国台湾地区-菲律宾群岛深部纵横波结构成像及其构造意义

使用地震层析成像方法反演了中国台湾地区-菲律宾群岛的1197126条P波震相和1217821条S波震相,首次同时得到了该地区从地表到100km的纵横波速度结构。成像结果揭示了沿着马尼拉海沟向东俯冲的欧亚板块从中国台湾西南部到吕宋岛南端构造形态上的变化;中国台湾中央山脉由于受到造山运动的影响,地壳厚度可达55~60km,而山脉两侧的地壳厚度多为20~35km;此外,成像结果还提供了菲律宾海板块在中国台湾东北部及琉球海沟下方北向俯冲的地球物理学证据;研究区域南部的菲律宾群岛由于同时受到两侧欧亚板块和菲律宾海板块的双向俯冲,使得岛内的岩石圈变形严重,岛弧岩浆和地震活动十分发育;群岛东侧的菲律宾海板块在东吕宋海槽的活动性较弱,很可能是本哈姆海台的碰撞和俯冲所造成的,而菲律宾海板块在菲律宾海沟的俯冲活动则十分强烈,随着深度的增加逐渐向南发展。研究表明,板块俯冲造成了中国台湾地区-菲律宾群岛的地壳及上地幔具有较强的不均一性,这不仅孕育了大量的地震和火山活动,同时也对该地区的地质构造产生了深远的影响。     

大别造山带地壳结构反演及其动力学意义

大别造山带是全球最大的碰撞造山带之一,三叠纪时期,扬子板块深俯冲至地幔的200km处,经历了超高压变质作用。白垩纪早期,该造山带发生了强烈的伸展和垮塌,以及大规模的后造山地幔源岩浆侵入和火山活动。本研究收集了大别造山带及其邻区(29°~34°N、114°~119°E)的震相资料,采用双差层析成像技术,对大别造山带地壳结构进行反演,研究地壳结构与后造山地幔源岩浆侵入和火山活动之间的关系。结果显示,大别造山带中上地壳存在低速结构,该低速结构可能是熔融的幔源侵入物质,由于俯冲板片断裂,或下地壳/岩石圈发生拆沉,导致软流圈物质上涌至地壳底部、侵入地壳中,形成大别造山带地壳中的低速结构;同时,合肥盆地显示为低速区,可能是受浅部沉积层影响。研究中横切大别山的4条剖面显示,该地区下方存在北向倾斜高速结构,该高速结构可能是襄樊-广济断层,或者是扬子板块向华北板块下方俯冲的遗迹。     

青藏高原东北缘上地幔多尺度层析成像

利用布设在青藏高原东北缘地区的甘肃宽频带地震台阵记录到的远震P波走时数据,采用小波域参数化和基于L1范数的稀疏约束反演算法的多尺度层析成像方法,得到了该地区400km深度范围内上地幔的P波速度结构.本文采用的多尺度层析成像方法可以自适应数据非均匀采样的情况,有效降低谱泄漏效应和反问题的多解性,明显提高解的分辨率和可靠性.层析成像结果表明青藏高原东北缘上地幔整体上显示为低速特征,扬子地块上地幔则显示为高速特征,两者之间上地幔存在清晰的块体边界带,该边界带位于东经104°—105°之间并且随深度的增加逐渐东移.该特征暗示了青藏高原上地幔物质向东扩张的机制,但在西秦岭上地幔顶部不存在物质运移的通道.青藏高原东北缘内部也具有明显的分区特征,松潘—甘孜地块上地幔P波速度整体呈低速特征,而柴达木地块的上地幔顶部具有相对高速特征,而在上地幔200km以下这两个地块间的差别逐渐减小.1654年天水地震和1879年武都地震都发生在扬子地块与青藏高原的碰撞交汇区,其震中下方上地幔显示为高低速转换结构.     

一阶Rytov近似有限频走时层析

传统的波动方程走时核函数(或走时Fréchet导数)多基于互相关时差测量方式及地震波场的一阶Born近似导出,其成立条件非常苛刻.然而,地震波走时与大尺度的速度结构具有良好的线性关系,对于小角度的前向散射波场,Rytov近似优于Born近似.因此,本文基于Rytov近似和互相关时差测量方式,导出了基于Rytov近似的有限频走时敏感度核函数的两种等价形式:频率积分和时间积分表达式.在此基础之上,本文提出了一种隐式矩阵向量乘方法,可以直接计算Hessian矩阵或者核函数与向量的乘积,而无需显式计算和存储核函数及Hessian矩阵.基于隐式矩阵向量乘方法,本文利用共轭梯度法求解法方程实现了一种高效的Gauss-Newton反演算法求解走时层析反问题.与传统的敏感度核函数反演方法相比,本文方法在每次迭代过程中,无需显式计算和存储核函数,极大降低了存储需求.与基于Born近似的伴随状态方法走时层析相比,本文方法具有准二阶的收敛速度,且适用范围更广.数值试验证明了本文方法的有效性.     

大地电磁、重力、磁法和地震初至波走时的交叉梯度二维联合反演研究

为了准确的探测和描绘地下复杂的地质结构,同时克服地球物理单一方法反演的多解性和单一参数反演模型的不一致性等问题,近年来基于交叉梯度联合反演的综合地球物理解释已经得到了广泛的关注和应用.本文首先研究了两种地球物理方法的交叉梯度联合反演算法,在此基础上,推导并实现了多种地球物理方法(大地电磁,重力,磁法,地震初至波走时)的多交叉梯度约束的二维联合反演算法;其次,我们设计了结构不一致模型和复杂模型,针对多物性联合反演算法的准确性和有效性进行了模拟试算,并对复杂模型的单独反演结果和联合反演结果进行了交叉梯度值和物性交会图的对比;最后,本文将成熟的卫星资料多光谱综合分析技术应用到联合反演中,将多物性参数反演模型结果图通过RGB(红-绿-蓝)模式进行合成,得到融合的RGB合成图.结果表明:通过对结构不一致模型和复杂模型的联合反演结果和单独反演结果的对比分析,可以得出联合反演得到的结果更接近真实模型,并从得到的交叉梯度值进一步证明了联合反演模型相似度高,也从物性交会图中得到联合反演的物性相关性更好的结论,反向证明了算法的正确性.最终从得到的RGB合成图像,我们可以更直观的分析反演结果,更有利于准确划分地下模型结构.     

一种新的预条件伴随状态法初至波走时层析

伴随状态法初至波走时层析是基于最优化理论的一种层析成像方法,该方法不必进行射线追踪,用两次正演的计算量便可以获得梯度,具有计算效率高、内存占用小等优点,但是其一阶方向在初始模型或观测孔径不理想的情况下往往无法获得正确的反演结果,而二阶方向的实现又比较困难且费时.在伴随状态法的基础上,将走时差替换为定值,再次进行反演,便可以得到类似于射线密度的矩阵,用该矩阵的逆可以方便地进行预条件.基于该方法,本文提出了一种简单易行的预条件伴随状态法初至波走时层析的实现方法.理论模型和实际资料处理结果都表明,该方法既保留了伴随状态法初至波走时层析的优点,又可以克服一阶方向的局限,获得良好的反演效果.     

基于Abel变换的图像重建自适应方法

本文论述了利用轴对称物体的单幅投影信息进行密度重建的一种自适应正则化模型。所提模型基于全变分正则项与高阶全变分正则项的联合使用,主要的优点是在保持清晰的界面及恢复平稳变化区域的同时减弱了阶梯效应。并且使用自适应方法,提高了效果的同时简化了所使用的参数。对于其中涉及的最优化问题,我们采用增广拉格朗日方法来解。数值结果表明,这一模型提高了关于密度界面位置及密度值的准确度,具有较好的抗噪性。