本文基于三维孔隙弹性理论,建立了紫坪铺水库及周边地区的有限元模型.根据紫坪铺水库开始蓄水到汶川地震发震时刻的水位变化情况,计算了整个区域的孔隙压力和库仑应力.详细讨论了断层及周围地层的弹性模量和扩散系数对计算结果的影响.计算结果表明:从弹性角度看,断层的弹性模量对汶川地震震源处的库仑应力影响很小;震源处的库仑应力随着断层和周围地层的扩散系数增大而增大.当给定弹性模量和扩散系数代表性值的时候,计算结果表明在汶川地震发震时刻,震源处的库仑应力变化量为+1 kPa左右,这表明紫坪铺水库使得汶川地震发震断层更加危险.是否这个量级的库仑应力就能够触发汶川地震还需要进一步探讨.通过分析库区周边小震的分布,发现小震分布区域均是库仑应力增加的地区,因此紫坪铺水库周边的小震应该与紫坪铺水库蓄水有直接关系.

Some large reservoirs are known to induce earthquakes. After occurrence of the Wenchuan earthquake, some studies suggest that it was induced by the Zipingpu reservoir which is not far from the Longmenshan fault. In this paper, based on the poroelastic theory, a 3D finite element model of the Zipingpu area is constructed. According to the water level of the Zipingpu reservoir from the beginning of the impoundment to the occurrence of the 2008 Wenchuan earthquake, we calculate the pore pressure and the change of Coulomb stress with the model. Furthermore, we discuss the effects of elastic modulus and diffusivity of faults and around medium on the results. The results show that the elastic modulus of faults and around medium has little effect on the Coulomb stress in the hypocenter of the Wenchuan earthquake. With the increase of diffusivity of the fault or around medium, the Coulomb stress in the hypocenter increases. If we give the representative values of elastic modulus and diffusivity, the change of Coulomb stress in the hypocenter reaches 1 kPa when the Wenchuan earthquake occurred. The Zipingpu reservoir enhanced the seismic activity of the Longmenshan fault. But it needs further investigation whether increasing 1 kPa of Coulomb stress can trigger the occurrence of the Wenchuan earthquake. The analysis of the distribution of small earthquakes in the Zipingpu area before the Wenchuan earthquake reveals that most of the small earthquakes are located in the increasing region of Coulomb stress. So it indicates that the small earthquakes in the Zipingpu area should have a direct relationship with the impoundment of the Zipingpu reservoir.