1976年唐山大震发生对华北地区各地块运动与变形影响的数值模拟研究.pdfVIP

( D 辑) 第33 卷 增 刊 SCIENCE IN CHINA ( Series D ) 2003 年4 月 1976 年唐山大震发生对华北地区各地块运动 与变形影响的数值模拟研究* 白武明 ** 林邦慧 陈祖安 ( 中国科学院地质与地球物理研究所, 北京 100101; 中国地震局地球物理研究所, 北京 100081) 摘要 用非连续变形分析方法(DDA+FEM) 数值模拟在华北地区各地块相互制约的 构造环境中, 发生 1976 年唐山地震的动力学过程; 研究大震引起的华北地区各地块特 别是鄂尔多斯地块运动变形及各地块边界断层上应力状态变化的特征. 采用的 1976 年 唐山地震复杂震源模型由 3 个子事件组成: 北北东走向的主破裂, 与主破裂北东端部 相连接的北东走向子事件以及位于主破裂南东侧的北西走向的子事件. 同前人的结果 比较, 这个模型与该震早期余震分布图像的位置更为相符. 计算得到了地震释放的主 应力场, 最大剪应力变化等值线图, 大震邻近地区地块滑移变形图, 发震断层滑移距离 随时间的变化以及最大走滑错距和应力降等震源参数. 这些结果与该震的震源机制 解, 震源参数, 宏观等震线, 地表观测的水平位移矢量图基本一致. 唐山地震的发生使 得与其邻近的地块及边界断层受到不同程度的影响, 发生了不同规模的移动和变 形. 鄂尔多斯地块似乎发生了整体性小规模的逆时针转动及变形, 其中鄂尔多斯地块 东北部与太行山, 燕山, 阴山地块交接带的相对运动变形较大. 唐山地震的发生也引起 了华北地区各地块边界断层上的应力状态发生了一定程度的变化. 其中通过唐山震区 的张家口-蓬莱北西走向断裂带以及鄂尔多斯地块北部边界断层的主要地带上剪应力增 加而法向应力减小, 使这些断裂带的地震危险度增加. 这与唐山地震后鄂尔多斯地块 北部边界断层上发生一系列6 级强震的事实基本相符. 关键词 非连续变形分析方法 数值模拟 唐山地震 震源过程 地块移动与变形 危险度 回顾20 世纪60~70 年代在华北地区地块系统的某些边界断层上发生一系列6, 7 级强震, 包括 1966 年邢台地震 1969 年渤海地震 1975 年海城地震及 1976 年和林格尔地震 唐山 地震等的过程, 可以发现, 地块边界断层的失稳似乎与地块系统中各地块的运动与变形的相 互影响有一定关系. 某个地震的发生, 不仅仅是震源自身的发震断层的应力积累和释放, 而且 [1~7] 与整个地块系统中地块运动与变形的相互影响以及引起的应力状态的改变有关 . 某个地块 边界断层上失稳发生地震, 不仅仅在该发震断层上产生应力降, 释放应变能, 同时由于各地块 的运动与变形的相互影响, 有可能向其周围的地块传递应力, 调整和改变其他地块边界断层 上的应力状态, 加速或延缓其他地块边界断层的失稳. 1976 年唐山地震发生后, 华北地区至今 再也没有发生 7 级以上大震. 但在地质上比较稳定的鄂尔多斯地块的北部边界断层上却发生 了一系列的6 级强震, 同时北京西北3 个省(河北, 山西, 内蒙古) 与鄂尔多斯地块交界地区震 情不断, 似乎显示了唐山大震对鄂尔多斯地块运动, 变形的影响. 然而这只是