陈 婷，章 阳，王晓山，孙丽娜

（河北省地震局，石家庄 050021）

0 引言

静态库仑应力变化最早应用在分析大地震对余震分布的影响中，研究认为当静态库仑应力改变量超过0.01 MPa（0.1 bar）就可以影响余震的发生位置[1]。强震或大地震引起的库仑应力变化不仅对余震的发生有影响，更重要的是对周边断层上的强震同样能产生触发作用，使处于库仑应力增加区域的地震提前发生[1-2]。近年来，国内众多学者综合分析断裂带上库仑应力动态演化与强震活动的关系，进而分析区域地震危险性[3-7]。

1966年3月邢台震群共发生5次6.0级以上地震，前人对邢台地震的库仑应力变化展开大量研究。万 永革等[8]研 究3月8日MS6.8地 震 对3月22日MS7.2地震的触发作用；陈连旺等[9]研究邢台地震长时间尺度的库仑应力变化；李守勇等[10]分析邢台地震序列2个地震间的相互触发关系，由于使用弹性模型计算产生较大误差，7.2级地震以后的触发关系不好；万永革等[11]研究认为粘弹性模型计算效果优于完全弹性；Wang R等[12]编写了粘弹性应力应变程序PSGRN/PSCMP程序，该程序基于Burgers粘弹性模型[13]，考虑地球重力对形变场的影响，可计算同震和震后长时间尺度的应力变化。

前人对邢台地震库仑应力的计算基于弹性模型假设和粘弹性有限元模型，研究内容也仅限于2个地震间的触发分析。为了更加接近岩石圈的真实地壳结构，本文基于Burgers粘弹性模型，综合考虑同震位错和震后粘滞松弛效应，采用PSGRN/PSCMP程序计算地震破裂产生的库仑应力变化，分析邢台地震序列中6级以上地震受该地震之前所有6级以上地震的综合库仑应力的叠加作用。……