弹性波的三维有限元模拟

本文讨论了弹性波瞬态传播问题的三维有限元计算方法及当前存在的实际困难.针对要求计算机内存大和计算时间长的问题,采取了改进措施.由于采用了集中质量矩阵和修正的中心差分时间积分显格式相结合的方法,可以使计算机内存和计算时间大为减少;由于采用结点定位法,最适合用于目前发展的并行计算机系统,可使计算速度大大增快;还采用了有效激发震源法,有效激发区是随时间步进的增加而逐步增大,这不仅能节省计算时间,而且使波场的传播过程一目了然,本文计算了由两种介质组成的三维楔形问题,得到若干典型剖面的瞬时波场图及随时间变化的合成地震图.

近场地震快慢横波到时差测量离散辨析和改正

快慢横波到时差是地震物理预报的重要线索之一.长期以来由于快慢横波到时差的较大离散,横波观测到时差变化与地震发生过程中应力变化之间无法建立能被绝大多数观测结果证实的物理联系.快慢横波测量离散既可能来自观测台站下方地层结构和介质的各向异性差异,也可能来自数据处理过程.本文利用新西兰南岛北部KHZ台站从2013年1月到2017年1月的单台近震数据,在观测台站底部介质物性和结构基本稳定的前提下,基于横波偏振矢量分析和可视化方法对横波数据进行处理,对横波处理过程中引起快慢横波倒时差测量离散的频率滤波、距离归一化因素进行了逐一的对比,并对台站下方介质各向异性分层可能引起的多次分裂所表现出来的"假离散"现象进行了分析,由此给出了改正横波测量离散的方法.用"有效震源距"作为到时差归一化路径长度得到的横波归一化到时差随震源距变化的关系图清晰地展示出7.8级地震前后应力变化的空间范围和影响程度,是较大地震前后观测地区存在应力积累的明确标识.给出了确定横波分裂界面深度的方法,得到了新西兰南岛北部KHZ台站下3个横波分裂界面深度.