詹自敏 陈连旺 陆远忠

(中国地震局地壳应力研究所，北京 100085)

大华北地壳动力学参数及三维有限元模型的建立*

詹自敏 陈连旺 陆远忠

(中国地震局地壳应力研究所，北京 100085)

利用华北地区的断层数据、高程数据、波速结构数据、弹性常数以及多种流变机制参数，建立大华北地区三维有限元数值模型，并通过给每个单元赋属性的方法，实现物性参数的连续变化。以模型边界附近实测GPS数据插值作为边界位移条件，进行数值模拟。模拟结果表明，大华北区域位移场与实测位移场基本一致。

华北;有限元模型;数值模拟;GPS;地壳动力学

1 引言

大华北地区多次发生7级以上强震，给人民的生命财产安全造成了极大的威胁。

有限元的方法是在有限差分的基础上发展起来的数值模拟方法，运用有限元方法建立的地块实体模型，可以模拟地块深部应力场的变化和强震的孕育的过程。建立大华北地区三位实体模型，需要比较详尽的参数，而华北地区大量的监测数据，对于建立华北地区有限元模型提供了良好的条件。

早在1980年王仁等［1］便对华北地区进行了二维有限元数值模拟，采用弹性材料，并用初始应力场和应力释放的方法模拟了邢台、河间、渤海、海城、唐山5次地震的发生;2003年，文献［2，3］用非连续变形和有限元相结合的方法对华北地区进行了二维数值模拟，同样采用了弹性的本构关系;1998年，刘昌铨等［4］建立了华北的三维有限元弹性模型;1999年及其后，陈连……