Abstract：Based on the viscoelastic finite element model, combined with the active fault detection in Handan area and historical earthquake data and distribution of the lithosphere medium in the studied area, we construct a three dimensional viscoelastic finite element fault model. We substitute the kinematics characteristic boundary conditions into the model, and simulate the changes of stress and strain field caused by the earthquake on the active faults. Then, we caculate the Coulomb failure stress changes of the rupture surface caused by the 1830 Cixian earthquake and the sliding direction. At last, the risk of strong earthquake in the target area is analyzed.According to the simulation model based on actual geology and current kinematic conditions, the Coulomb stress near the epicenter of the original strong earthquake of the Cixian fault is reduced which indicates that the activity of the Cixian fault reduces seismic risk near the epicenter of the original strong earthquake. But Coulomb stress is strengthened near the junction of the middle and west sections of the Cixian fault and the two northeast trending faults (Zishanxi fault and Taihang mountain piedmont fault), and the seismic risk in the future should be monitored.