强地震形变过程的有限元数值模拟及其动力学环境的研究

【摘要】： 通过全球定位系统(GPS)和卫星孔径干涉雷达(InSAR)等大地测量观测强地震引起的同震和震后形变,结合地质构造条件分析,利用有限元数值模拟,使得我们可以对地震区域地球动力学环境有更加详细地了解,进而对地震孕育和地壳动力学演化等地质过程有更深入的认识和理解。本文利用固体力学理论,针对不同粘弹介质、孔隙介质以及粘弹—孔隙等介质的本构方程,建立和导出了其适用于FEPG(有限元程序生成系统)的虚功方程。以GPS观测所得地表形变做约束,用FEPG(有限元程序生成系统)实现了强地震引起的同震和震后形变有限元数值模拟,并掘此分析了以介质物性为中心的强地震区域的环境动力学特征。 本文在前人震后形变分析和模拟基础之上,根据粘弹松弛和流体调整所引起的震后形变的物理含义,建立孔隙—粘弹本构方程以及适应于FEPG平台的有限元数值模拟所必需的虚功方程。同时针对震后形变所表现出的粘弹松弛特点,以Burgers体粘弹模型特点以及本构方程出发,建立相应虚功方程,实现了以Burgers体为粘弹介质模型的震后形变有限元数值模拟。以Burgers体为粘弹介质模型探讨了强地震引起的三个阶段的形变效应,同震的瞬时弹性瞬态响应、呈指数衰减的短期形变响应和线性增加的长期稳态形变响应。 在理论模拟和分析的基础上,本文以实际地震为例详细讨论地震引起的同震和震后形变过程。 本文以1999年Mw7.6级集集地震为例,讨论不同的位错模型和介质模型(均匀、层状均匀及横向不均匀等)对同震形变的影响。通过数值模拟分析本文认为,在条件允许的情况,在讨论地震形变时有必要考虑小断层位错模型和更加接近实际的非均匀介质模型,计算和探讨了集集地震引起的库仑应力变化及其与集集地震余震对应关系。 本文以2001年Ms8.1级地震为例,以GPS观测形变为约束,通过有限元数值模拟,解释GPS所观测的震后形变表现出较为奇特的区域特征。用数值模拟结合网格搜索法确定了羌塘块体和柴达木盆地下地壳的粘滞系数分别约为5.0×10~(17)Pa.s,9.0×10~(18)pa.s相差十余倍。正是由于昆仑断裂带南北两侧下地壳粘滞系数的这十余倍的差异引起了断层两侧震后形变格局明显不同,这也可能是羌塘块体和柴达木盆地表现为不同现代地壳运动的直接原因。与此同时,本文还通过数值模拟定性地讨论了断层北侧地表形变在震后短期内转向的现象。我们发现,对于靠近断裂附近的测点上的形变转向也许是粘弹性松弛和孔隙流体调整共同作用的结果,而单独的粘弹性松弛或者孔隙流体调整不能引起这种转向。所以在分析短期震后形变时综合考虑粘弹—松弛介质物性和孔隙流体调整很有必要。 本文探讨了震后不同时间尺度的形变机理问题。研究表明,简单的粘弹介质模型(Maxwell体)不能同时解释震后短期和长期的形变过程,该模型表现出对同一地震同一地区而言,使用震后不同时期内形变资料,将得到不同的介质粘滞性质。本文证明了。以Burgers体为粘弹介质模型可以较好地解决使用Maxwell体等简单的粘弹介质模型所不能解释的震后形变过程随时间变化的问题,即模拟预测表现出来的震后短期和长期形变不协调的现象。通过数值模拟分析了1960年Mw9.5级智利地震震后形变过程,较成功地解释了该强烈地震震后短期和长期形变观测数据。 本文的研究表明,以观测地表形变为约束,结合实际地质条件,选取合理的介质物性模型,构建有限元模型的框架,利用有限元数值模拟的方法,研究强地震形变过程和规律,是分析了以介质物性为中心的强地震区域的环境动力学特征的有效方法。