## 1 Introduction
[17,18] used FLAC3D software(2018) to calculate the accumulated energy of the disturbed surrounding rock and used seismology theory to predict the degree of rockburst risk at each mining level. [PERSON] et al. [19,20] constructed a real three-dimensional geo-mechanical model of adiamond mine in Canada, and they used the tangential stress criterion and the Burst Potential Index (BPI) based on the energy principle to evaluate the degree of rockburst risk in the stope and the roadway and identify potential rockburst risk areas. Based on theoretical analysis and numerical simulation, [PERSON] et al. [21] studied the distribution characteristics of high-stress-concentration zones and elastic strain energy of surrounding rock around sharply inclined thin vein mining, and they analyzed the rockburst risks of stope with different mining sequences. [PERSON] et al. [22] predicted the rockburst of the surrounding rock in the roadway of the Linglong Gold Mine based on strain energy analysis using three-dimensional numerical simulation. [PERSON] et al. [23] introduced the strain energy density index (SED) in the numerical simulation to research the characteristics of energy accumulation and dissipation in the process of dynamic instability, and the simulation results were in better agreement with the field observation. [PERSON] et al. [24] introduced the rock ultimate energy storage rate in conventional triaxial loading and unloading tests and proposed a new Rockburst Energy Release Rate index (RBERR) for assessing rockburst risk. Adopting a three-dimensional discrete element method and an independently developed multi-factor rockburst criterion, [PERSON] et al. [25] researched and analyzed the distribution of stress field, energy characteristics, damage characteristics, and the degree of rockburst risk of the surrounding rock of the tunnel under the conditions of in situ stress at different buried depths. [PERSON] et al. [26] considered that the rockburst risk is directly proportional to the maximum energy density and inversely proportional to the distance between the location of the maximum energy density and the free face, proposing a new energy density criterion, and combined with the numerical simulation method and the independently proposed energy density criterion, the effect of the tunnel mining sequence on the rockburst risk of the pillar was examined. As mentioned above, most research has focused on the rockburst risk assessment of a single roadway or stope with consistent lithology. However, few efforts have been made to evaluate the rockburst risk of deep-buried roadways under multi-stratigraphic tectonics. There are differences in the rockburst characteristics that occur when constructing in hard, harder, and soft rock, and the rockburst risk varies when excavating from soft rock towards hard rock or in the opposite direction. [PERSON] et al. [3] described the situation of the rockburst that occurred in interbedded siltstone and sandstone formations in the N-J hydroelectric power station in Pakistan, and the rockburst risk from sandstone to siltstone is higher than during excavation from siltstone to sandstone. Therefore, rockburst risk assessment for deep-buried roadways traversing the vicinity of the stratigraphic boundary under multi-stratigraphic tectonics is of great engineering significance. This study is based on a developing roadway near the 14# exploration line in the middle section of 1104 of the Huize lead-zinc mine, which crosses different formations as the engineering background and establishes a three-dimensional geo-mechanical model through the actual engineering geological investigation on site.