模拟地震波传播的三维逐元并行谱元法

高效地震波场正演模拟对于复杂模型中地震波传播与成像研究至关重要.本文在谱元法原理框架内,对已有逐元谱元法改进,提出一种新的逐元并行谱元法求解三维地震波运动方程,并得到地震波场.逐元并行谱元法的核心思想在于在单元上进行质量矩阵与解向量的乘积运算,并将此运算平均分配至每一个CPU计算核心,此处理有利提升谱元法的并行计算效率.同时,根据Gauss-Lobatto-Legendre (GLL)数值积分点与插值点重合的特点,将稠密单元刚度矩阵的存储转化成单元雅克比矩阵行列式的值及其逆的存储,大幅减少谱元法计算内存开销.此外,在模型边界上利用逐元并行谱元法求解二阶位移形式完美匹配层(PML)吸收边界条件,消除边界截断而引入的虚假反射.通过逐元并行谱元法得到的数值解与解析解对比,以及实际地震波场模拟,数值结果证实了逐元并行谱元法用于地震波场模拟的高效性.

     

基于迭代反演的层间多次波压制方法

层间多次波与有效波在走时、频率和叠加速度上差异较小,因此层间多次波压制常难以获得理想效果.本文提出一种基于迭代反演的层间多次波压制方法（MSI,Multiple Suppression Inversion）,该方法以共聚焦点（CFP,Common Focus Point）层间多次波压制理论为基础,通过构建卷积因子,将层间多次波压制转变为迭代反演的问题,直接利用观测的地震数据进行迭代反演计算,进而完成多次波压制.MSI方法避免了共聚焦点方法中构建CFP道集的聚焦运算,大幅降低层间多次波预测的计算成本;同时该方法为全数据驱动方法,无需地下介质任何先验信息,算法容易实现.模型数据测试表明,本文提出的方法可有效压制层间多次波,而且对有效波也具有很好的保幅性.

     

面波频散谱多模式高分辨率成像的多道信号比较法

高分辨率的面波频散谱成像是浅层地震勘探领域基于频散性质反演横波速度结构中的一个关键步骤.在天然地震探测领域,仅利用两个台站记录的线性信号比较法(LSC),被广泛用来计算面波的频散谱,并用于大尺度的面波层析成像.然而互相关的成像方式会造成频散谱在低频端较低的分辨率.非线性信号比较法(NLSC)利用指数函数克服了这个问题,...     

利用接收函数散射核方法探测中国东北地区地幔转换带界面三维形态

地幔转换带是上下地幔物质运移和能量交换的必经通道,其速度结构和上下界面的起伏能够为认识地幔内部的温度和物质变化、地幔对流模式等地球演化相关的科学问题提供关键约束.本文利用布设在中国东北地区的高密度固定台网和流动台阵所记录的远震体波接收函数,采用Ps散射核叠前深度偏移成像方法,获得了台站下方地幔转换带界面及其内部速度间断面的三维精细图像.研究结果表明：西北太平洋俯冲板片上下界面在地幔转换带内清晰可见,在高纬度（44°N）区域存在约30°的倾角;660-km间断面的深度起伏具有明显的分区性,在与俯冲板片相交处以西200~300 km,界面出现约20~40 km的下沉,而长白山和龙岗火山的西北区域存在约5~15 km的抬升,分别与板片滞留引起的低温异常和局部热物质上涌相对应;410-km间断面的起伏形态较复杂,在大部分区域观测到大于10 km的下沉,且表现出明显的区域性横向变化,与深俯冲动力学背景下冷的温度异常造成的影响不一致.我们认为板片俯冲、停滞和海沟后撤过程中引起的地幔转换带物质异常、含水状态及分布的变化是显著改变410-km间断面形态的主要原因.本文获得的高精度地幔转换带界面三维形态为更好地认识东北亚地区俯冲板片在地幔物质分布和能量交换中的作用提供了重要参考.

     

基于耐震时程法的进水塔结构抗震性能分析

进水塔结构作为水电站的进水建筑物,其抗震性能对水电站的安全运行有重要意义.以沙牌水电站的高耸独立进水塔结构为研究对象,应用3条耐震时程曲线对其进行抗震性能分析与评估,并将计算结果与天然地震动下的增量动力分析结果进行比较,对比分析进水塔塔顶位移、基底剪力及地震损伤历程.实例研究表明:耐震时程法可以较好地分析结构的抗震响应...     

新疆喀什地区1996年11月震前超低频电磁辐射异常与震源方向的关系

研究了1996年11月19-27日在新疆喀什地区连续发生6次4.5级以上地震前记录到的超低频电磁辐射异常现象的偏振特性,发现震前记录到的该异常用0.0-0.16Hz的带通滤波器滤波后,水平磁场分量的偏振近似为线偏振,且震前1天的偏振主轴方向日均值与喀什台到强震震源的方向之间有近乎垂直的关系.     

标定大陆科学钻探孔区地震反射体

在中国大陆科学钻探孔区(江苏省东海县南部)进行了系统的地球物理调查,包括二维地震测网和专门的地震剖面,大地电磁法和位场方法等.地震调查表明,在超高压变质岩出露区上地壳充满了反射体,包括倾斜反射体与上拱的反射弧. 本文介绍大陆科学钻探先行研究中地震调查的成果.根据大陆科学钻探预先导孔的岩芯和测井资料、井旁VSP和数值模拟结果证实,高波速的榴辉岩体、破碎断裂带和大型韧性剪切带都可引起倾斜的地震反射,而上拱的弧形反射体则是由近似直立的榴辉岩体和其中的破碎带的综合反映.由于地壳深部广泛分布着经受变质的岩石,上述研究结果对标定地壳中的反射地震信号具有一定意义.     

双相各向异性介质中弹性波方程的有限元解法及波场模拟

基于双相各向异性介质模型,首先推导了双相各向异性介质中弹性波传播的动力学方程及其Galerkin变分方程和有限元运动方程,然后给出了孔隙弹性波方程的有限元数值解法以及二维双相PTL介质中波场模拟的人为吸收边界条件. 最后,利用本文给出的有限元方法对双相PTL介质和双相各向同性介质中的弹性波传播进行了数值模拟. 结果表明:有限元方法和吸收边界条件有效、可行,在理想相界条件下,不论是从固体位移,还是从流体位移的波场快照都能看到明显的慢速拟P波;在黏滞相界情况下,能否观察到慢速拟P波,与含流体地层介质的耗散性质有关.对实际含流体介质,从流体位移分量的波场快照比从固体位移波场快照更容易观察到慢速拟P波.     

青藏高原北部异常SKS分裂成因的初步探讨--被熔体强化的岩石圈各向异性

当人们试图解释青藏高原异常的剪切波分裂成因时,以下的问题让人们感到困惑:(1)为什么异常大的SKS分裂延时(1.91-2.4s)出现在青藏高原北部Sn波缺失区;(2)为什么分裂延时突变(1.47s和1.09s)出现在Sn波缺失区的边缘;(3)为什么快波极化方向(FPD)与地表大规模的构造走向之间存在约20°-30°的偏差. 本文在综合分析流变学实验和岩石物理学实验研究成果、青藏高原地质和地球物理资料的基础上,提出青藏高原北部地震波各向异性受岩石圈地幔主要矿物的晶格优选方位(LPO)和熔体的定向分布(MPO)的双重控制,并模拟计算了MPO对青藏高原北部岩石圈地幔各向异性强度的贡献. 研究结果表明,由MPO强化的青藏高原北部岩石圈地幔各向异性强度可达10髎,相应的各向异性层厚度平均为94km. 该结果为研究区SKS分裂的成因解释以及造山带深部地质过程的研究提供了新的约束条件.     

重庆云阳新津乡中侏罗统新田沟组恐龙的发现及沉积环境

在重庆市云阳县新津乡中侏罗统新田沟组中发现了大量恐龙化石,且恐龙化石产出层位呈现出与双壳类、鱼类、龟鳖类和蛇颈龙类化石整体混杂埋藏的特征,这是新田沟组中首次发现。通过野外实测剖面、岩石特征、沉积构造、古生物特征、粒度特征,初步判定研究区新田沟组为滨-浅湖亚相沉积,局部为洪泛亚相。重庆地区新田沟组古生物化石组合的首次发现为四川盆地中生代古生物化石及其古环境的研究提供了更多依据。