Update and Test of a Seismicity Simulation Program Based on the Coulomb-Navier Failure Criterion and the Stress Transfer Mechanism

Using a seismicity simulation program based on the Coulomb-Navier failure criterion and the stress transfer mechanism, we added extensions to make the model compatible with different input files and to increase computational efficiency. After testing with a two-vertical-parallel-fault model, we applied this model to the Taiyuan region. The frequency–magnitude data of the simulation results fit well with the observation catalogue, which at the same time tends to follow the characteristic earthquake theory instead of the Gutenberg-Richter law. In addition, a single fault region with a higher slip rate seems to have lower b-values, suggesting a higher risk of strong earthquakes. We have also performed data fitting between the simulation results and the Poisson process of corresponding parameters, along with testing for energy conservation.