地球动力学若干问题的计算仿真研究

论文摘要

运用三维有限元动力学分析方法,结合地球物理观测数据,研究了地壳介质中应力的传播特性,得出了大陆地壳对海洋潮汐的响应、仿真地震图、全球模型、断裂面对地震波的影响、动应力对断层错动的影响等研究结果。基于后两项内容,探讨了地震触发的动力学机理。在此基础上,进一步研究了由此形成的地震波与断层活动之间的正反馈关系,解释了大地震发生有一个长时间酝酿过程这一现象。通过研究,对地壳中的力学现象从应力的传播和影响角度作出物理解释。同时探讨了计算力学与地球物理观测相结合的可行性。本文的主要内容包括:（1）将三维动态有限元方法应用于地球动力学分析,研究了大尺度的地球物理对象,本文致力于深入地进行三维动力学分析与固体地球的动力学过程的联合研究,在方法上,着眼于地球物理的动力学现象,对地球物理现象作动力学机理的探讨。地球物理的研究对象具有大尺度特征,通常的研究,往往着眼于缩小的模型,再外推到大尺度的实际情况,外推必须非常谨慎。研究小尺度模型而外推到大尺度,涉及到相似性分析,其正确性还需要论证,比较繁琐。本文利用了动态有限元建模,直接研究大尺度的对象。在大陆地壳对海洋潮汐的响应、仿真地震图、全球模型、断裂面对地震波的影响、动应力对断层错动的影响等问题的研究中,均采用了真实尺寸地壳建模,进行仿真计算。得到了直观的结果。（2）利用反演思想实现了仿真地震图研究,基于传统的理论地震图的研究方法,本文提出仿真地震图方法。针对地壳的有限元模型,做了实际的仿真地震图分析,与地震观察结果进行对比研究。理论地震图是一种解析方法,而仿真地震图基于数值方法。前者采用一种解析的格林函数,代表波在地壳中的传播特性,而后者对地壳介质及其结构具有明确的定义,对波的传播过程有具体的描述。仿真地震图属于反演方法,通过假设介质和震源参数,得到理论上的仿真地震图,将实际观察到的地震图与之比较,用试错法调整介质和震源参数,最终使仿真地震图与观察到的地震图符合较好,这时的介质和震源参数就是所求出的未知量。（3）研究了地震触发的动力学机理:地震机理有很多模式,本文提出了一种可能的模式。这个模式指出,库伦应力的变化可能触发地震。动态应力的传播即地震波,波在自由面上反射会发生半波损失,也就是说,压缩波会变成拉伸波,因此,无论是压缩波还是拉伸波,都有可能激发应力降。而应力下降将会减小（内）摩擦,导致库伦-莫尔断裂准则中的库伦等效应力增大,从而触发岩石破裂现象。岩石的破裂会导致地震波的辐射,因此,应力的下降和岩石的破裂会相互激发,形成一个正反馈环路,最终触发大规模的岩石破裂事件,这就是大地震。本项研究从冲击载荷下的断裂现象入手,论述拉伸波的作用,考虑到围压的影响,指出拉伸波虽然可以在爆破中直接引起岩石的断裂,但是在地层深处,不能直接引起断裂事件。因此引出了内摩擦和库伦应力的影响机理。在此基础上,理解到大地震的触发,为一个不断通过内摩擦改变形成的正反馈过程,解释了大地震发生有一个漫长的地应力下降过程和局部变形过程这一实际现象,从而对已有的地震的触发机理给出了合理的补充。

论文目录

• 摘要
• Abstract
• 1 引言
• 1.1 国内外研究现状
• 1.2 地震波研究
• 1.3 研究内容和主要思路
• 2 非线性动态有限元原理
• 2.1 非线性动态有限元的特点
• 2.2 动态有限元原理的数学表述
• 2.3 总体控制方程
• 2.4 计算方案
• 2.5 软件功能介绍
• 2.6 小结
• 3 中国中北部地壳对潮汐的响应
• 3.1 研究内容和方案
• 3.2 建立模型
• 3.3 载荷条件
• 3.4 计算结果及分析
• 3.5 海潮应力与地震的关系
• 3.6 模拟结果与孕震体应力环境的关系
• 3.7 海洋潮汐载荷与地震孕育过程的关系
• 3.8 与孕震介质的关系
• 3.9 小结
• 4 仿真地震图
• 4.1 研究背景和研究方案
• 4.2 有限元仿真研究
• 4.3 仿真地震图实施过程
• 4.4 小结
• 5 断裂面对动态应力传播的影响
• 5.1 概述
• 5.2 断层面反射模型
• 5.3 模拟结果
• 5.4 一个形状不规则的模型
• 5.5 结合观测的分析
• 5.6 小结
• 6 动应力对断层摩擦力的影响
• 6.1 研究内容和技术路线
• 6.2 建模和计算仿真
• 6.3 对模拟结果的物理解释
• 6.4 小结
• 7 展望
• 7.1 环球剪切带
• 7.2 动量矩迁移
• 7.3 全球模型
• 7.4 小结
• 8 全文总结
• 8.1 完成的研究工作
• 8.2 创新之处
• 8.3 展望
• 致谢
• 参考文献
• 附录 攻读博士学位期间已发表的论文