基于三维层状介质模型地震层析成像正演网格研究

地震层析成像技术是研究地层速度结构的一个重要手段,该技术是利用震源到接收点的地震波走时或波形的观测数据,结合建立的数学物理模型来重建地层速度结构。而层析成像正演模拟的精度将直接影响反演出的速度与地层真实速度的拟合度,因此找出一种精度较高的层析成像正演模拟算法是非常重要的。本文针对三维层状介质,通过网格化模型,采用最小走时射线追踪方法展开层析成像正演研究,分别采用三角形网格、矩形网格和六边形网格进行模型参数化,并对比分析各自正演结果的精度,总结出针对三维层状介质模型的最佳网格划分方式。      

一种锥束CT动态Bowtie的蒙特卡罗模拟研究

在CT中,前置于患者的衰减器（"Bowtie滤波器"或"Bowtie"）被用来减少患者所受的辐射剂量和帮助满足探测器的动态范围,同时Bowtie也能提高CT图像质量。传统的静态滤波器对调制X射线能力有限,我们之前提出了一种能够有效调制X射线的动态Bowtie滤波器。本文使用蒙特卡罗模拟方法设计和研究提议动态Bowtie滤波器。首先,建立了一个仿真框架,包括一个目标试样和一个由反射靶生成的多能谱射线源。接着,一个设计试样被用来辅助设计针对于目标试样的动态Bowtie滤波器。最后,我们分别从没有Bowtie滤波器和带有动态Bowtie滤波器两种模式中,模拟获得两个试样的投影数据并进行分析研究。结果显示,在设计的动态Bowtie作用下,到达探测器的X射线动态范围明显减小、X射线光子通量一致性很好,并且动态Bowtie可以减小辐射剂量和提高CT图像质量。      

西太平洋俯冲板块对中国东北构造演化的影响及其动力学意义

中国东北地区在古生代期间以众多微陆块的拼合以及古亚洲洋的闭合为特征,其后又经历了中-新生代太平洋构造域及中生代蒙古—鄂霍茨克构造域的叠加与改造,以致东北地区的构造行迹显得极为复杂,而大兴安岭重力梯级带及其西部地区构造演化是否与西太平洋俯冲有关仍然存在争议.本研究利用分布于中国东北、华北地区以及韩国、日本等部分台网所接收的近震与远震走时数据获得了中国东北地区壳幔精细的三维P波速度结构.成像结果显示,太平洋板块持续西向俯冲,俯冲板片的前缘停滞在大兴安岭—太行山重力梯度带以东区域的地幔转换带之中;长白山火山区上地幔存在着显著的低速异常体,推测西太平洋板块的深俯冲脱水导致了上地幔底部岩石的熔点降低,从而形成了大范围的部分熔融物质上涌.通过分析上地幔的速度结构,我们认为由于太平洋板块的大规模西向深俯冲,在大地幔楔中发生板片脱水、低速热物质上涌等复杂的地球动力学过程;俯冲板片前缘带动上地幔中不均匀分布的地幔流强烈作用于上部的岩石圈,这对东北地区深部壳幔结构乃至大兴安岭重力梯级带的形成、演化有着重要的影响.     

2017年米林6.9级地震震源区速度结构与余震重定位

利用西藏自治区林芝地区的固定地震台站与南迦巴瓦流动测震台站在2017年11月18日至2017年11月24日记录到的430个余震的直达波走时数据反演得到了震源区的三维P波速度、S波速度结构,并利用三维速度结构对余震进行了重定位.成像结果显示,米林地震震源区在0~5km深度内存在低地震波速度异常;在5~15km深度内,存在高地震波速度异常,该高速异常致使震源区西南侧的地震波速度高于东北侧.重定位结果中,余震呈条带状以NW-SE走向展布,震源深度具有西南方向深、东北方向浅的特征.主震位于11km深度处、高地震波速异常体顶部,余震主要分布在高地震波速度与低地震波速度过渡的区域.对成像结果的分析表明,震源区浅部的低速异常具有低泊松比的特性,与富石英的沉积变质杂岩体-东久杂岩单元的岩性特征有关;深部的速度结构特征则可能反映了发震断层上盘地震波速度高,下盘地震波速度低的介质特性.余震重定位结果与成像结果联合表明:此次地震发震断层从11km深度处,东久杂岩体下方的高地震波速度异常顶部开始破裂,继而在5~15km深度内发生后续破裂,后续破裂的发生区域正处于喜马拉雅构造单元与冈底斯构造单元接触的形变区内.此外,根据地震波速度计算的泊松比反映了震源区持续的低泊松比特征,暗示此次地震与流体活动并无直接关系.     

龙门山断裂带精细速度结构的双差层析成像研究

利用川西流动地震台阵、汶川地震震后应急台网记录到的P波到时资料,对2008年5月至2008年10月期间发生的汶川地震余震序列应用双差层析成像方法进行了地震震源和三维P波速度结构的联合反演.结果显示,联合反演获得的地震重定位结果与基于一维地壳参考模型的双差定位方法结果相近;研究区15 km以上速度结构与地表断裂分布密切相关,20 km以下深度呈现北东向和北西向交错结构.汶川地震破裂带南段龙门山断裂带之间上地壳呈现高速异常,速度结构的非均匀变化是控制余震分布和主震破裂传播的主要因素;联合反演结果给出了小鱼洞-理县方向存在隐伏断裂的速度结构证据,同时发现,破裂带北东段可能沿新发断裂扩展;结果确认了汶川地震起始段的高角度逆冲断裂特征,也确认了前山断裂和中央断裂在约20 km深度合并到脆韧转换带的特征.     

苏鲁造山带地壳上地幔结构层析成像研究

扬子与华北板块在三叠纪的俯冲碰撞形成了著名的苏鲁超高压造山带,其板块碰撞接触关系一直是热点问题.利用国家台网中心64个省台记录的1 079个近震事件的10 922个P波到时和251个远震事件的11 931个P波到时数据,采用远近震联合反演的层析成像方法对苏鲁地区进行了地壳上地幔速度结构反演.结果显示,研究区内两个低速异常区分别对应山东半岛西部的华北板块地幔上隆区和壳幔相互作用强烈的长江中下游成矿带地区.在地幔300 km深度之下出现的高速异常体可能代表了早中生代扬子与华北板块碰撞之前俯冲拆沉的古特提斯洋板块.传统观点的扬子板块岩石圈向北俯冲不明显,华北板块表现为向东南俯冲的高速特征.华北板块俯冲以苏鲁造山带中部的北纬35°为界,分为南北两种俯冲样式.北部俯冲不明显,华北板块停滞在郯庐断裂带以西;南部则表现华北板块向东南陡倾俯冲到苏鲁造山带之下.      

鄂尔多斯地块北部及邻区Pn波速度结构与各向异性

利用鄂尔多斯地块北部及其邻区2008—2018年期间固定台网的地震波形记录,手动拾取出高质量的Pn波到时资料,反演获得了研究区上地幔顶部的Pn波速度结构及各向异性。结果表明:鄂尔多斯地块北部及其邻区上地幔顶部的Pn波速度存在明显的横向不均匀性,与区域地质构造和地震活动相关;研究区内平均Pn波速为8.18 km/s,鄂尔多斯地块内部表现出大范围的高速异常,阿拉善地块的高速异常体中存在低速异常现象,河套断陷带、阴山—燕山造山带、银川—吉兰泰断陷带和海原—六盘山弧形断裂带区域均表现为显著的低速异常,河套断陷带下方存在向鄂尔多斯地块内部延伸的明显低速异常条带,大同火山群下方存在强低速异常;多数历史强震均发生在低速异常区或高低速异常过渡带上;鄂尔多斯地块内部Pn波各向异性快波方向西部为近NE?SW向,而东部为近NW?SE向,河套断陷带和鄂尔多斯地块西缘、青藏高原东北缘与阿拉善地块的交界带以及阴山—燕山造山带的各向异性快波方向总体均呈现为NW?SE向,而阴山—燕山造山带东部则呈NE?SW向。      

爆炸当量估算的初步研究

爆炸当量是评估爆炸能量的重要参数,快速估算爆炸当量对于有效开展应急救援工作至关重要。收集攀枝花铁矿爆破、宁夏煤矿爆破、河南汝阳钼矿爆破等矿山爆破资料,得到爆破的地方性震级M_L与爆炸当量Y之间的经验关系,并据此估算2019年7月19日河南省三门峡市气化厂爆炸事故的爆炸当量为16 t左右。     

弹性波逆散射微扰论

本文将普遍声逆散射微扰论应用于弹性波层析成像问题,在Born变换下推出了以旋转角为补偿参数的各阶微扰重建公式,实现了对非均匀各向同性散射体内3个参数（质量密度ρ和两个Lamé系数λ,μ）的同时重建. 对于层析成像问题,在弹性波的传播过程中P波与SV波有耦合,但它们不会和SH波发生耦合,于是可以得到3个形式相对简单的标量方程. 在Born变换下,在散射波中引入微扰参数,将散射体的3个参数分别按该微扰参数展开,然后利用二维自由空间的Green函数分别得到散射的P波、SV波和SH波的积分表示. 最后,经一维傅氏变换后,得到Born变换下散射体3个参数的各阶微扰重建公式.