Welcome to EARTHQUAKE,

EARTHQUAKE ›› 2013, Vol. 33 ›› Issue (4): 162-171.

Previous Articles     Next Articles

Comparison of 2D and 3D Finite Element Simulations in the Research of Reservoir Induced Earthquakes

ZHENG Liang1,2, ZHANG Huai1,2, SUN Yu-jun1,2, CHENG Hui-hong1,2, ZHANG Bei1,2, SHI Yao-lin1,2   

  1. 1. Key Laboratory of Computational Geodynamics, CAS, Beijing 100049, China;
    2. College of Earth Science, University of CAS, Beijing 100049, China
  • Received:2013-01-11 Revised:2013-03-18 Published:2020-09-27

Abstract: It has been widely discussed in geoscience community whether Coulomb stress change caused by Zipingpu reservoir loading activities has triggered the 2008 Wenchuan MS8.0 earthquake. Many conclusions have come from 2D finite element simulations, which are quite different with the 3D results. In this paper, we quantitatively discussed the difference of the 2D and the 3D results. First, a comparison of the 3D finite element solutions and the analytic solutions was given to ensure the enough precision of the 3D finite element simulation with 1 million grid points. Then, the 2D and the 3D numerical solutions of the elastic effect and the seepage effect caused by Zipingpu reservoir were compared. The simulation results show that elastic load caused the thrust fault to be more secure, while the water pressure caused it to be more dangerous. The superimposed effect, however, didn't result an increase in the risk of Wenchuan earthquake. The comparison of the 2D and the 3D results show the 2D result exaggerated the truth. At the hypocenter of the Wenchuan earthquake, the computed stress of the 3D result is only 1/3 to 1/4 of the 2D result, thus in a serious study, 3D modeling is necessary.

Key words: Finite element numerical simulation, Coulomb stress, Wenchuan Earthquake, Zipinpu reservoir earthquake

CLC Number: