Study on tectonic fractures based on the inversion of tectonic stress fields is an effective method. In this study, a geological model was set up based on geological data from the Hudi Coal Mine, Qinshui Basin, a mechanical model was established under the condition of rock mechanics and geostress, and the finite element method was used to simulate the paleotectonic stress field. Based on the Griffith and Mohr‐Coulomb criterion, the distribution of tectonic fractures in the Shanxi Formation during the Indosinian, Yanshanian, and Himalayan period can be predicted with the index of comprehensive rupture rate. The results show that the acting force of the Pacific Plate and the India Plate to the North China Plate formed the direction of principal stress is N‐S, NW ‐ SE, and NE ‐ SW, respectively, in different periods in the study area. Changes in the direction and strength of the acting force led to the regional gradients of tectonic stress magnitude, which resulted in an asymmetrical distribution state of the stress conditions in different periods. It is suggested that the low‐stress areas are mainly located in the fault zones and extend along the direction of the fault zones. Furthermore, the high‐stress areas are located in the junction of fold belts and the binding site of multiple folds. The development of tectonic fractures was affected by the distribution of stress intensity and the tectonic position of folds and faults, which resulted in some developed areas with level I and II. There are obvious differences in the development of tectonic fractures in the fold and fault zones and the anticline and syncline structure at the same fold zones. The tectonic fractures of the Shanxi Formation during the Himalayan period are more developed than those during the Indosinian and Yanshanian period due to the superposition of the late tectonic movement to the early tectonic movement and the differences in the magnitude and direction of stress intensity.