岩石圈介质结构的不均匀性对地震同震、震后位移影响的三维数值模拟.pdf

摘要 摘要 地震同震位移场的模拟受很多参数控制，从震源角度看：诸如断层的几何 结构，倾角的大小，位错的分布等；从介质结构角度看：如半空间均匀模型、 横向分层以及纵向分层介质模型：从弹性参数角度看：如拉梅常数，弹性模量 以及泊松比等。定量的分析这些参数对地震同震位移的影响有助于我们根据地 震破裂过程了解地震造成的地表形变特征或利用观测的形变资料识别断层的活 动特征。在粘弹性松弛的震后形变理论中，震后形变则受控于同震形变和介质 的粘性结构。探讨各种粘性结构介质中震后位移场的特征，有助于我们利用震 后形变资料研究岩石的粘性。 本文首先回顾了同震形变、震后形变理论研究的背景。然后介绍了有限元 方法的基本理论。在此基础之上，定量探讨了介质的不均匀性对地震同震、震 后位移场的影响。最后分别以大陆内部和俯冲区域发生的逆冲大地震：汶川地 震、日本地震为例，构建符合当地发震区域的介质结构，在观测资料的约束下， 定量的模拟介质的不均匀性对地震同震位移以及震后位移的影响。论文主要包 括以下内容： (I)介质的不均匀性和断层倾角对同震位移场的影响 影响同震位移分布的主要因素包括断层的滑动分布、断层的几何参数以及 介质的弹性参数。本项研究主要利用三维有限元数值模拟的方法，定量研究介 质不均匀性和断层倾角变化对同震位移场的影响。而滑动分布则选用简单的有 限矩形走滑、逆冲或正断层滑动。数值模拟显示：(1)在均匀及纵向分层介质 模型中，对于垂直走滑断层地震，水平位移场分量与介质的泊松比呈较弱负相 关，垂直位移场与介质的泊松比呈正相关，水平位移场与剪切模量呈正相关， 垂直位移场与剪切模量呈负相关：(2)模型介质的横向变化对同震位移场具有 较大影响，剪切模量的横向变化对地震的位移场影响最大，断层左侧块体地震 位移场与左侧块体介质的剪切量呈负相关，在垂直走滑断层地震中，断层右侧 块体弹性参数不变、左侧块体剪切模量减半可使左侧块体的最大垂直位移分量 幅度增大55．6％；(3)倾角对同震位移场有重要影响，在断层附近，倾角对位 移场有主要的控制作用，对于高倾角逆冲断层(60--90度)，上盘断层附近区 域的水平位移场出现反向区，反向区域的范围随倾角的增大而增大，当倾角为 90度时，上盘水平位移场全部反向；断层倾角增大时，在断层附近下盘的水平 位移场增大幅度较大。 摘要 (II)介质粘弹性结构的不均匀性对震后位移场的影响 在粘弹性松弛震后形变理论中，介质的粘弹性结构在震后位移场的特征中 起着重要的控制作用，传统的半解析解无法考虑介质粘性的横向差异。我们构 建简单的走滑断层地震，采用有限元的数值方法，定量的分析不同的粘弹性介 质及粘性结构不均匀介质中地震震后形变的特征。数值模拟显示：(1)麦克斯 韦介质在模拟地震震后形变中存在不协调问题，当粘性较低(如1018Pa．s)时， 粘弹性松弛会很快达到平衡，无法模拟时间跨度较长的震后形变，当粘性较高 时，震后位移随时间成线性增大，不符合一般观测的震后位移呈指数变化特点， 采用广义麦克斯韦体介质，能较好的改善这一问题；(2)纵向分层(果冻三明 治模型、奶油焦糖布丁模型)的粘弹性介质模型中，震后位移随时间松弛趋势 与半空间均匀介质中的震后位移松弛趋势一致。相对于上地幔，下地壳的粘性 对震后位移场的影响较大，在相同时间内，较弱的下地壳模型要比较弱的上地 幔模型中震后位移松弛幅度大；(3)在横向不均匀的粘弹性模型中，粘性相对 较强的块体中断层附近震后位移分量呈现与同震位移运动方向相反的特点，在 相同的时间内粘性较弱的块体中震后位移量大于粘性较强块体中的震后位移。 (4)断层一侧块体含有粘性较弱的流动通道时，其震后位移的松弛特点与下地 壳上地幔粘性横向不均匀模型中的震后位移松弛特征相同。 (IⅡ)汶川Ms8．0地震的同震、震后位移场 汶川Ms8．0地震发震区域的介质背景以及断层破裂滑动分布比较复杂，我们 采用有限元算法，分别利用由Jiand and Hayes以及NishimuruYaji给出的滑动分 布模型，参考地震波反演给出的发震区域的介质参数，计算汶川地震的同震位移 场，比较断层滑动分布对地震同震位移场的影响，并将地表同震位移场的计算值 与GPS观测值对比。此外参考CRUST2．0构建汶川地震局部区域的纵向分层介质 模型，比较