Mechanical Response and Stability Optimization of Shallow-Buried Tunnel Excavation Method Conversion Process Based on Numerical Investigation

: Different excavation methods lead to substantial changes in the relaxation zone of the surrounding rock. The stress characteristics of the support structure become more complex during the process of excavation method conversion. It is essential to design a well-considered construction plan to minimize the disruptions caused by excavation method conversion. This work takes a tunnel in Jiangsu Province, China, as the engineering background and establishes a numerical model for the conversion from the double-side-drift method to the three-bench method. Based on a numerical investigation, this study analyzed the deformation of surrounding rock and the stress state of the support structure during the excavation method conversion. The results showed that excavation method conversion can accelerate the deformation rate of the surrounding rock, causing the support structure to bear greater pressure. This leads to a significant increase in the circumferential stress at the junction of two excavation methods. Setting up a reinforced area during the process of excavation method conversion can improve the stress state of the support structure. Finally, an analysis of the parameters of the reinforced area was conducted, and a reasonable construction plan was proposed. This study can provide guidance for subsequent construction projects.


Introduction
In tunnel engineering, geological conditions such as soil properties have a significant impact on the mechanical response of support structures [1,2].Therefore, it is necessary to choose appropriate excavation methods, according to different tunnel sizes and geological conditions [3].When the surrounding rock conditions are good, the bench method and full-section method can be employed for rapid construction [4,5].If the conditions are poor or the tunnel span is large, the subsection excavation methods can be used to divide the tunnel into several parts, which can enhance the stability of the tunnel [6].The subsection excavation methods have been widely applied in practical engineering projects.These methods include the double-side-drift method (DSD), the center diaphragm method (CD) and the cross diaphragm method (CRD) [7,8].
Many scholars have conducted research on the selection of excavation methods.Some studies focused on specific engineering aspects and compared tunnel safety when employing different excavation methods [9][10][11].A small number of studies have summarized extensive engineering experience and proposed a series of criteria for selecting excavation methods [12][13][14].During the construction of long tunnels, various excavation methods are often combined due to factors such as design changes, tight schedules, and unexpected changes in surrounding rock conditions [15].Different excavation methods will lead to substantial changes in the relaxation zone of the surrounding rock.The stress characteristics of the support structure become more complex during the process of excavation method conversion.Therefore, improper conversion schemes can easily lead to tunnel instability, causing unfavorable impacts on the construction [16].It is essential to design a well-considered construction plan to effectively minimize the disruptions caused by the excavation method conversion.
However, the current research on the excavation method conversion is limited.Zhang et al. [17] used the Jiaoweiqin Tunnel as the research subject and provided a clear exposition of the excavation method transition plan.The analysis of tunnel deformation confirmed the safety of the method transition.This study primarily focused on whether the excavation method conversion is feasible under current conditions and did not investigate the mechanical response of support structures.In order to effectively prevent tunnel instability during the process of excavation method conversion, it is necessary to conduct a systematic analysis of the support stress characteristics, with the proposal of reasonable reinforcement measures.Advanced support has been often employed in tunnel engineering to address issues such as weak surrounding rock, support deformation, and face instability caused by tunnel construction [18][19][20].However, there is still a lack of comprehensive research on advanced support schemes for the excavation method conversion process.
To address the deficiencies of existing research, this study takes a tunnel project within Jiangsu Province as a case study and investigates the deformation and stress characteristics of the support structure during the conversion from the DSD method to the three-bench method.Based on numerical results and field monitoring, the response mechanisms of the tunnel were summarized, and reasonable reinforcement measures were proposed.The research findings can provide guidance for future projects of a similar nature.

Project Overview and Model Establishment 2.1. Tunnel Geological Conditions
The tunnel is situated at Jiangsu Province, China, with a total length of 660 m.The maximum burial depth of the tunnel is 75 m.According to the tunnel design, the excavation span is 17.39 m, and the excavation area is 163.34 m 2 .As displayed in Figure 1, the tunnel exit section (from ZK10+660 to ZK10+780) has a relatively shallow burial depth.The maximum burial depth is 17.8 m, the minimum burial depth is 4.6 m.The main components of the surrounding rock are shale and limestone.According to the BQ rock mass classification system (China) [21], it is classified as Class V.As displayed in Figure 2, the middle section and the portion above the tunnel crown are composed of highly weathered mud shale, which exhibits overall fragmentation.The lower section consists of moderately weathered mud shale with thick bedding planes.will lead to substantial changes in the relaxation zone of the surrounding rock.The stress characteristics of the support structure become more complex during the process of excavation method conversion.Therefore, improper conversion schemes can easily lead to tunnel instability, causing unfavorable impacts on the construction [16].It is essential to design a well-considered construction plan to effectively minimize the disruptions caused by the excavation method conversion.However, the current research on the excavation method conversion is limited.Zhang et al. [17] used the Jiaoweiqin Tunnel as the research subject and provided a clear exposition of the excavation method transition plan.The analysis of tunnel deformation confirmed the safety of the method transition.This study primarily focused on whether the excavation method conversion is feasible under current conditions and did not investigate the mechanical response of support structures.In order to effectively prevent tunnel instability during the process of excavation method conversion, it is necessary to conduct a systematic analysis of the support stress characteristics, with the proposal of reasonable reinforcement measures.Advanced support has been often employed in tunnel engineering to address issues such as weak surrounding rock, support deformation, and face instability caused by tunnel construction [18][19][20].However, there is still a lack of comprehensive research on advanced support schemes for the excavation method conversion process.
To address the deficiencies of existing research, this study takes a tunnel project within Jiangsu Province as a case study and investigates the deformation and stress characteristics of the support structure during the conversion from the DSD method to the three-bench method.Based on numerical results and field monitoring, the response mechanisms of the tunnel were summarized, and reasonable reinforcement measures were proposed.The research findings can provide guidance for future projects of a similar nature.

Tunnel Geological Conditions
The tunnel is situated at Jiangsu Province, China, with a total length of 660 m.The maximum burial depth of the tunnel is 75 m.According to the tunnel design, the excavation span is 17.39 m, and the excavation area is 163.34 m 2 .As displayed in Figure 1, the tunnel exit section (from ZK10+660 to ZK10+780) has a relatively shallow burial depth.The maximum burial depth is 17.8 m, the minimum burial depth is 4.6 m.The main components of the surrounding rock are shale and limestone.According to the BQ rock mass classification system (China) [21], it is classified as Class Ⅴ.As displayed in Figure 2, the middle section and the portion above the tunnel crown are composed of highly weathered mud shale, which exhibits overall fragmentation.The lower section consists of moderately weathered mud shale with thick bedding planes.

Excavation Methods and Support Structure
The original construction plan for the tunnel was to use the DSD method or the CRD method for the Class V rock mass, and the three-bench method for the Class IV rock mass.The tunnel exit section was originally constructed using the DSD method.After excavation to ZK10+727, it was converted to the three-bench method to accelerate the excavation speed.As shown in Figure 3, the tunnel face is segmented into seven parts using the DSD method.The excavation sequence is ➀→➁→➂→➃→➄→➅; the tunnel face is segmented into three parts using the three-bench method.The excavation sequence is ➀→➁→➂.The excavation footage is 2 m.When using the DSD method, the upper steps of the left and right sides lead the lower steps by 4 m, the upper steps in the middle lead the lower steps by 2 m, the left side leads the right side by 6 m, and the right side leads the middle step by 6 m.When using the three-bench method, the upper step leads the middle step by 4 m, and the middle step leads the lower step by 4 m.
The parameters of the support structure are as follows: C25 shotcrete with a thickness of 28 cm; φ8 steel mesh with a spacing of 20 cm × 20 cm; I22a steel arch with a longitudinal spacing is 70 cm; D25 hollow grouting bolt with a length of 4.0 m and a spacing of 70 cm × 100 cm.

Model and Parameters
This paper adopts the birth-death element within the finite element method to simulate the processes of excavation of rock mass and the installation and removal of support structures.Select ZK10+745-ZK10+705 of the tunnel for numerical simulation.According

Excavation Methods and Support Structure
The original construction plan for the tunnel was to use the DSD method or the CRD method for the Class V rock mass, and the three-bench method for the Class IV rock mass.The tunnel exit section was originally constructed using the DSD method.After excavation to ZK10+727, it was converted to the three-bench method to accelerate the excavation speed.As shown in Figure 3, the tunnel face is segmented into seven parts using the DSD method.The excavation sequence is 1  ⃝→ 2 ⃝→ 3 ⃝→ 4 ⃝→ 5 ⃝→ 6 ⃝; the tunnel face is segmented into three parts using the three-bench method.The excavation sequence is 1  ⃝→ 2 ⃝→ 3 ⃝.

Excavation Methods and Support Structure
The original construction plan for the tunnel was to use the DSD method or the CRD method for the Class V rock mass, and the three-bench method for the Class IV rock mass.The tunnel exit section was originally constructed using the DSD method.After excavation to ZK10+727, it was converted to the three-bench method to accelerate the excavation speed.As shown in Figure 3, the tunnel face is segmented into seven parts using the DSD method.The excavation sequence is ➀→➁→➂→➃→➄→➅; the tunnel face is segmented into three parts using the three-bench method.The excavation sequence is ➀→➁→➂.The excavation footage is 2 m.When using the DSD method, the upper steps of the left and right sides lead the lower steps by 4 m, the upper steps in the middle lead the lower steps by 2 m, the left side leads the right side by 6 m, and the right side leads the middle step by 6 m.When using the three-bench method, the upper step leads the middle step by 4 m, and the middle step leads the lower step by 4 m.
The parameters of the support structure are as follows: C25 shotcrete with a thickness of 28 cm; φ8 steel mesh with a spacing of 20 cm × 20 cm; I22a steel arch with a longitudinal spacing is 70 cm; D25 hollow grouting bolt with a length of 4.0 m and a spacing of 70 cm × 100 cm.

Model and Parameters
This paper adopts the birth-death element within the finite element method to simulate the processes of excavation of rock mass and the installation and removal of support structures.Select ZK10+745-ZK10+705 of the tunnel for numerical simulation.According The excavation footage is 2 m.When using the DSD method, the upper steps of the left and right sides lead the lower steps by 4 m, the upper steps in the middle lead the lower steps by 2 m, the left side leads the right side by 6 m, and the right side leads the middle step by 6 m.When using the three-bench method, the upper step leads the middle step by 4 m, and the middle step leads the lower step by 4 m.
The parameters of the support structure are as follows: C25 shotcrete with a thickness of 28 cm; φ8 steel mesh with a spacing of 20 cm × 20 cm; I22a steel arch with a longitudinal spacing is 70 cm; D25 hollow grouting bolt with a length of 4.0 m and a spacing of 70 cm × 100 cm.

Model and Parameters
This paper adopts the birth-death element within the finite element method to simulate the processes of excavation of rock mass and the installation and removal of support structures.Select ZK10+745-ZK10+705 of the tunnel for numerical simulation.According to Saint-Venant's principle, the dimensions of the numerical model should be at least three times the diameter of the excavated tunnel in order to eliminate boundary effects [22].
Consequently, the width, length and height of the numerical model were 120, 40 and 80 m.The thickness of the surrounding rock above the tunnel is 15.9 m, based on the average burial depth.Considering the actual constraint effect of surrounding rock, we apply normal constraints to the four sides of the model, fix constraints at the bottom and leave the top surface unconstrained [23].This model assumes homogeneity for both the rock mass and the support structure.The rock mass is simulated by 8-node hexahedron elements and follows the Mohr Coulomb criterion.Meanwhile, the support structure adopts a linear elastic model and is simulated by 4-node quadrilateral elements.The finite element model is illustrated in Figure 4.The focus of this study is on the variation patterns of deformation and stress characteristics of supporting structures.Therefore, the geological conditions in the model have been simplified, disregarding the effects of tectonic stress and groundwater.
Buildings 2024, 14, x FOR PEER REVIEW 4 of 16 to Saint-Venant's principle, the dimensions of the numerical model should be at least three times the diameter of the excavated tunnel in order to eliminate boundary effects [22].Consequently, the width, length and height of the numerical model were 120, 40 and 80 m.The thickness of the surrounding rock above the tunnel is 15.9 m, based on the average burial depth.Considering the actual constraint effect of surrounding rock, we apply normal constraints to the four sides of the model, fix constraints at the bo om and leave the top surface unconstrained [23].This model assumes homogeneity for both the rock mass and the support structure.The rock mass is simulated by 8-node hexahedron elements and follows the Mohr Coulomb criterion.Meanwhile, the support structure adopts a linear elastic model and is simulated by 4-node quadrilateral elements.The finite element model is illustrated in Figure 4.The focus of this study is on the variation pa erns of deformation and stress characteristics of supporting structures.Therefore, the geological conditions in the model have been simplified, disregarding the effects of tectonic stress and groundwater.The mechanical parameters of the rock mass were obtained from the 'Geological Engineering Handbook' [24].The mechanical parameters of the support structure were obtained from the 'Specifications for Design of Highway Tunnels' (JTG 3370.1-2018)[25].To simplify the numerical model, an equivalent substitution is made by converting the elastic modulus of the steel arch and steel mesh into shotcrete.The equivalent formula is given in Equation ( 1) [26].The mechanical parameters of the rock mass and the support structure are listed in Table 1.
where Ec and E represent the initial and modified elastic modulus of shotcrete; Ea and Em are the elastic modulus of the steel arch and steel mesh; Sc, Sa, Sm are the cross-sectional areas of the shotcrete, steel arch and steel mesh, respectively.The mechanical parameters of the rock mass were obtained from the 'Geological Engineering Handbook' [24].The mechanical parameters of the support structure were obtained from the 'Specifications for Design of Highway Tunnels' (JTG 3370.1-2018)[25].
To simplify the numerical model, an equivalent substitution is made by converting the elastic modulus of the steel arch and steel mesh into shotcrete.The equivalent formula is given in Equation ( 1) [26].The mechanical parameters of the rock mass and the support structure are listed in Table 1.
where E c and E represent the initial and modified elastic modulus of shotcrete; E a and E m are the elastic modulus of the steel arch and steel mesh; S c , S a , S m are the cross-sectional areas of the shotcrete, steel arch and steel mesh, respectively.

The Impact of Excavation Method Conversion on Tunnel Response
To investigate the impact of the excavation method conversion on tunnel safety, numerical simulations were conducted for the two construction conditions: (1) The DSD method is employed for excavation within the range of ZK10+745-ZK10+705 (no conversion).( 2) The DSD method is employed for excavation within the range of ZK10+745-ZK10+727; the three-bench method is employed for excavation within the range of ZK10+727-ZK10+705 (conversion).
ZK10+727 is the junction point between the DCD method and the three-bench method.Therefore, the numerical results at ZK10+727 are selected for analysis, including tunnel deformation and the stress state of the support structure.

Deformation Analysis of the Surrounding Rock
During the excavation process of a shallow-buried tunnel, the pressure arch cannot naturally form in the rock mass above the tunnel.The self-bearing capacity of the surrounding rock is not fully utilized.The surrounding rock pressure is approximately equal to the weight of the rock mass above the tunnel.Therefore, vertical displacement is one of the crucial indicators reflecting the stability of the tunnel.The isoline map of vertical displacement at the completion of the construction simulation is illustrated in Figure 5. Vertical displacement significantly increases under the influence of the excavation method conversion.The vault settlement increases from 7.64 mm to 8.56 mm, with a difference of 0.91 mm; the invert uplift increases from 6.65 mm to 7.53 mm, with a difference of 0.88 mm.The evolution process of vault settlement and invert uplift at ZK10+727 is shown in Figure 6.Since there are significant differences in the construction procedures between the DSD method and the three-bench method, it is difficult to compare each excavation step with one another.Therefore, the distance between the ZK10+727 section and the excavation face is used as the abscissa to uniformly describe the two different excavation processes.It can be observed that both the vault settlement and invert uplift undergo three stages: rapid deformation, continuous deformation and slow deformation.At the end of the continuous deformation stage, the difference in vault settlement between the two construction conditions is 0.67 mm, accounting for 73.6% of the final difference; the difference in invert uplift is 0.87 mm, representing 98.9% of the final difference.This indicates that the impact of excavation method conversion on tunnel response is mainly evident in the early stages of construction and gradually increases as construction progresses.At this time, the initial support structure of the tunnel has just been constructed, and excessive displacement development will cause the initial support to bear greater pressure, thus causing hidden dangers to the safety of tunnel construction.

Stress Analysis of the Support Structure
The circumferential stress nephogram of shotcrete is shown in Figure 7.It can be observed that there is an obvious stress increase at ZK10+727.The distribution of the circumferential stress in shotcrete at ZK10+727 is illustrated in Figure 8. Stress is defined as positive in compression and negative in tension.When the DCD method is employed for the This indicates that the impact of excavation method conversion on tunnel response is mainly evident in the early stages of construction and gradually increases as construction progresses.At this time, the initial support structure of the tunnel has just been constructed, and excessive displacement development will cause the initial support to bear greater pressure, thus causing hidden dangers to the safety of tunnel construction.

Stress Analysis of the Support Structure
The circumferential stress nephogram of shotcrete is shown in Figure 7.It can be observed that there is an obvious stress increase at ZK10+727.The distribution of the circumferential stress in shotcrete at ZK10+727 is illustrated in Figure 8. Stress is defined as positive in compression and negative in tension.When the DCD method is employed for the entire excavation process, the shotcrete is under compression, and the circumferential stress is uniform and symmetrically distributed.After the excavation method conversion, the circumferential stress of the shotcrete at the arch and sidewalls tends to increase, with the maximum increase of 1.90 MPa occurring at the left wall waist (E).The maximum stress appears at the left arch foot (D), which is 1.66 MPa higher than that without conversion.It is noteworthy that tensile stress occurs at the middle invert (H), with a value of 0.22 MPa.
During the construction process, the stress state of the support structure continuously changes.In order to propose effective reinforcement measures to enhance the stability of the tunnel, it is essential to analyze the evolution of circumferential stress in shotcrete.Due to the symmetrical distribution of circumferential stress, an analysis at positions A-H is shown in Figure 9.
In these two conditions, the trend of circumferential stress development in shotcrete is significantly different.When excavation is carried out using the DCD method, the evolution process of stress can be divided into two stages.The first stage involves excavation of the left and right sides.As shown in Figure 9a, during this stage, the excavation volume of the rock for each excavation step is small, resulting in a slow increase in stress.The second stage involves the excavation of the surrounding rock in the middle, during which the stress at each point increases rapidly.Point C is located at the arch shoulder, where it bears the pressure from the surrounding rock.Therefore, the stress increase is most significant at this juncture.Before the start of the second stage, the support on both sides has already completed closure, allowing it to maintain stability well in the face of stress increase.Overall, stress development is gradual.
entire excavation process, the shotcrete is under compression, and the circumferential stress is uniform and symmetrically distributed.After the excavation method conversion, the circumferential stress of the shotcrete at the arch and sidewalls tends to increase, with the maximum increase of 1.90 MPa occurring at the left wall waist (E).The maximum stress appears at the left arch foot (D), which is 1.66 MPa higher than that without conversion.It is noteworthy that tensile stress occurs at the middle invert (H), with a value of 0.22 MPa.During the construction process, the stress state of the support structure continuously changes.In order to propose effective reinforcement measures to enhance the stability of the tunnel, it is essential to analyze the evolution of circumferential stress in shotcrete.Due to the symmetrical distribution of circumferential stress, an analysis at positions A-H is shown in Figure 9. entire excavation process, the shotcrete is under compression, and the circumferential stress is uniform and symmetrically distributed.After the excavation method conversion, the circumferential stress of the shotcrete at the arch and sidewalls tends to increase, with the maximum increase of 1.90 MPa occurring at the left wall waist (E).The maximum stress appears at the left arch foot (D), which is 1.66 MPa higher than that without conversion.It is noteworthy that tensile stress occurs at the middle invert (H), with a value of 0.22 MPa.
(a) (b)  During the construction process, the stress state of the support structure continuously changes.In order to propose effective reinforcement measures to enhance the stability of the tunnel, it is essential to analyze the evolution of circumferential stress in shotcrete.Due to the symmetrical distribution of circumferential stress, an analysis at positions A-H is shown in Figure 9. Due to the accelerated deformation rate of the surrounding rock, the support structure bears greater vertical pressure.As shown in Figure 9b, the development of stress is highly uneven after the conversion.The circumferential stress at each point is already significant from the beginning and increases rapidly.With just one excavation step, the circumferential stress at point D (left arch foot) has already increased to 7.00 MPa, accounting for 77.95% of the final value.After the support is closed, the increase in stress slows down significantly.This is unfavorable for construction stability.
In addition, tensile stress occurs at point H (invert) and continues to decrease within the two excavation steps after the support closure.In shallow-buried tunnels, the lateral pressure from the surrounding rock is relatively low, meaning the horizontal force is insufficient to counteract the vertical squeezing, which causes tension at the invert.However, the tensile stress is minor and well below the tensile strength of the shotcrete.Safety can be ensured by enhancing monitoring and timely implementing secondary support for the invert.
In summary, the conversion of excavation methods can adversely affect the stress state of the support structure.The main risk is the excessive early stress at point D. It is necessary to take measures to enhance construction stability.In these two conditions, the trend of circumferential stress development in shotcrete is significantly different.When excavation is carried out using the DCD method, the evolution process of stress can be divided into two stages.The first stage involves excavation of the left and right sides.As shown in Figure 9a, during this stage, the excavation volume of the rock for each excavation step is small, resulting in a slow increase in stress.The second stage involves the excavation of the surrounding rock in the middle, during which the stress at each point increases rapidly.Point C is located at the arch shoulder, where it bears the pressure from the surrounding rock.Therefore, the stress increase is most significant at this juncture.Before the start of the second stage, the support on both sides has

Leading Conduit Grouting Reinforcement
Leading conduit grouting is a commonly used reinforcement measure in tunnel excavation.As shown in Figure 10, the cement grout, surrounding rock, and small conduits are closely integrated, forming a reinforced area.The mechanical properties of the surrounding rock within the reinforced area are enhanced.

Leading Conduit Grouting Reinforcement
Leading conduit grouting is a commonly used reinforcement measure in tunnel excavation.As shown in Figure 10, the cement grout, surrounding rock, and small conduits are closely integrated, forming a reinforced area.The mechanical properties of the surrounding rock within the reinforced area are enhanced.In order to optimize the stress state of the support structure of after ZK10+727, leading conduits are set within the range of ZK10+729-ZK10+725.According to the 'Guidelines for Design of Highway Tunnel' (JTG/T D70-2010) [27], φ42 × 4 mm leading conduits are selected with a length of 3.5 m; the spacing is set at 40 cm × 240 cm, with a layout scope of 120°.The thickness of the reinforced area is 1 m, and the se ing angle of the small conduit is 9.2°.The finite element model of the reinforcement area is shown in Figure 11.Furthermore, the mechanical parameters of the reinforced area are adjusted to match those of the rock mass to simulate the grouting of the leading conduit [28].The mechanical parameters of the leading conduit and the reinforced area are listed in Table 2.In order to optimize the stress state of the support structure of after ZK10+727, leading conduits are set within the range of ZK10+729-ZK10+725.According to the 'Guidelines for Design of Highway Tunnel' (JTG/T D70-2010) [27], φ42 × 4 mm leading conduits are selected with a length of 3.5 m; the spacing is set at 40 cm × 240 cm, with a layout scope of 120 • .The thickness of the reinforced area is 1 m, and the setting angle of the small conduit is 9.2 • .The finite element model of the reinforcement area is shown in Figure 11.Furthermore, the mechanical parameters of the reinforced area are adjusted to match those of the rock mass to simulate the grouting of the leading conduit [28].The mechanical parameters of the leading conduit and the reinforced area are listed in Table 2.The distribution of circumferential stress in shotcrete at ZK10+727 is shown in Figure 12.It can be observed that the circumferential stress of the arch and sidewall decreases.This indicates that the reinforced area can bear a portion of the surrounding rock weight, reducing the rock pressure acting on the support structure.The distribution of circumferential stress in shotcrete at ZK10+727 is shown in Figure 12.It can be observed that the circumferential stress of the arch and sidewall decreases.This indicates that the reinforced area can bear a portion of the surrounding rock weight, reducing the rock pressure acting on the support structure.
Figure 13 illustrates the process of circumferential stress development at point.The final stress decreases from 8.98 MPa to 8.57 MPa, representing a reduction of 4.57%; early stress decreases from 7.00 MPa to 6.61 MPa, representing a reduction of 5.57%.This indicates that the optimization effect is limited under the current design, and it is necessary to further study the parameters of the reinforced area.The distribution of circumferential stress in shotcrete at ZK10+727 is shown in Figure 12.It can be observed that the circumferential stress of the arch and sidewall decreases.This indicates that the reinforced area can bear a portion of the surrounding rock weight, reducing the rock pressure acting on the support structure.

Parameter Studies of the Reinforced Area
The grouting volume, grouting radius, range and spacing of the leading conduit, among other parameters, will influence the reinforcement effect [29].The above parameters are simplified to match the length, thickness, and strength of the reinforced area.The length of the reinforced area is as shown in Table 3.The thickness of the reinforced area is correlated with the se ing angle of the small conduit, as shown in Table 4.According to the research in reference [30], the elastic modulus has the greatest influence on the strength of the reinforced area.To simplify the calculations, the Poisson's ratio, cohesion and internal friction angle remain constant.The values of elastic modulus and gravity density are given in Table 5, taking the early stress and final stress at point D as the criteria for evaluation.

Parameter Studies of the Reinforced Area
The grouting volume, grouting radius, range and spacing of the leading conduit, among other parameters, will influence the reinforcement effect [29].The above parameters are simplified to match the length, thickness, and strength of the reinforced area.The length of the reinforced area is as shown in Table 3.The thickness of the reinforced area is correlated with the setting angle of the small conduit, as shown in Table 4.According to the research in reference [30], the elastic modulus has the greatest influence on the strength of the reinforced area.To simplify the calculations, the Poisson's ratio, cohesion and internal friction angle remain constant.The values of elastic modulus and gravity density are given in Table 5, taking the early stress and final stress at point D as the criteria for evaluation.
As shown in Figure 14, with the increase in the length of the reinforced area, both the final stress and the early stress show a trend of first increasing and then stabilizing.When the length increases from 4 m to 6 m and 8 m, the final stress decrease by 0.07 MPa and 0.031 MPa, respectively, while the early stress decrease by 0.118 MPa and 0.039 MPa, respectively.This indicates that increasing the length of the reinforced area has a limited effect on optimizing the stress state of the support structure.When the length reaches 8 m, further increasing the length does not further reduce the circumferential stress.This suggests that the influence range of the reinforced area is limited.In subsequent calculations, the length of the reinforced area is set to 8 m.As shown in Figure 14, with the increase in the length of the reinforced area, both the final stress and the early stress show a trend of first increasing and then stabilizing.When the length increases from 4 m to 6 m and 8 m, the final stress decrease by 0.07 MPa and 0.031 MPa, respectively, while the early stress decrease by 0.118 MPa and 0.039 MPa, respectively.This indicates that increasing the length of the reinforced area has a limited effect on optimizing the stress state of the support structure.When the length reaches 8 m, further increasing the length does not further reduce the circumferential stress.This suggests that the influence range of the reinforced area is limited.In subsequent calculations, the length of the reinforced area is set to 8 m.From Figure 15a, it is observed that when the thickness of the reinforced area increases from 1 m to 1.5 m, the final stress decreases from 8.47 MPa to 8.33 MPa; however, when the thickness is further increased, the circumferential stress changes li le, with only a 0.007 MPa reduction when the thickness increases from 1.5 m to 2 m.From Figure 15b, it can be seen that the trend of change in the early stress is different from that of the final stress.When the thickness increases from 1 m to 1.5 m, the early stress decreases by 0.249 MPa; as the thickness continues to increase, the trend of reduction in early stress significantly slows down, with only a 0.095 MPa reduction when the thickness increases from 1.5 m to 2 m.This suggests that the circumferential stress will not indefinitely decrease with the increase in the thickness of the reinforced area.From Figure 15a, it is observed that when the thickness of the reinforced area increases from 1 m to 1.5 m, the final stress decreases from 8.47 MPa to 8.33 MPa; however, when the thickness is further increased, the circumferential stress changes little, with only a 0.007 MPa reduction when the thickness increases from 1.5 m to 2 m.From Figure 15b, it can be seen that the trend of change in the early stress is different from that of the final stress.When the thickness increases from 1 m to 1.5 m, the early stress decreases by 0.249 MPa; as the thickness continues to increase, the trend of reduction in early stress significantly slows down, with only a 0.095 MPa reduction when the thickness increases from 1.5 m to 2 m.This suggests that the circumferential stress will not indefinitely decrease with the increase in the thickness of the reinforced area.Based on the mechanism of the leading conduit, increasing the thickness of the reinforced area enhances the bearing capacity of the surrounding rock over a larger zone.At the same time, the se ing angle of the small conduit increases with the thickness.When the se ing angle is less than 15°, the beam-arching effect of the small conduit can be effectively utilized; when the angle exceeds 15°, the small conduit primarily functions as anchors.Additionally, considering the site construction, an excessive angle of the small conduit is not conducive to grouting control.Therefore, the thickness of the reinforced area should not exceed 1.5 m.
As shown in Figure 16a, when the elastic modulus of the reinforced area increases from 1 GPa to 1.5 GPa, the final stress decreases by 0.159 MPa; when it increases from 1.5 GPa to 2 GPa, the final stress only decreases by 0.013 MPa.This is because, as the strength of the reinforced area increases, its bearing capacity correspondingly improves; at the same time, the weight of the rock mass above is also increasing, ultimately acting on the support structure.As can be seen from Figure 16b, the early stress uniformly decreases with the increase in the strength of the reinforced area.When the elastic modulus increases from 1 GPa to 2 GPa, the early stress decreases by a total of 0.337 MPa.Based on the mechanism of the leading conduit, increasing the thickness of the reinforced area enhances the bearing capacity of the surrounding rock over a larger zone.At the same time, the setting angle of the small conduit increases with the thickness.When the setting angle is less than 15 • , the beam-arching effect of the small conduit can be effectively utilized; when the angle exceeds 15 • , the small conduit primarily functions as anchors.Additionally, considering the site construction, an excessive angle of the small conduit is not conducive to grouting control.Therefore, the thickness of the reinforced area should not exceed 1.5 m.
As shown in Figure 16a, when the elastic modulus of the reinforced area increases from 1 GPa to 1.5 GPa, the final stress decreases by 0.159 MPa; when it increases from 1.5 GPa to 2 GPa, the final stress only decreases by 0.013 MPa.This is because, as the strength of the reinforced area increases, its bearing capacity correspondingly improves; at the same time, the weight of the rock mass above is also increasing, ultimately acting on the support structure.As can be seen from Figure 16b, the early stress uniformly decreases with the increase in the strength of the reinforced area.When the elastic modulus increases from 1 GPa to 2 GPa, the early stress decreases by a total of 0.337 MPa.Based on the mechanism of the leading conduit, increasing the thickness of the reinforced area enhances the bearing capacity of the surrounding rock over a larger zone.At the same time, the se ing angle of the small conduit increases with the thickness.When the se ing angle is less than 15°, the beam-arching effect of the small conduit can be effectively utilized; when the angle exceeds 15°, the small conduit primarily functions as anchors.Additionally, considering the site construction, an excessive angle of the small conduit is not conducive to grouting control.Therefore, the thickness of the reinforced area should not exceed 1.5 m.
As shown in Figure 16a, when the elastic modulus of the reinforced area increases from 1 GPa to 1.5 GPa, the final stress decreases by 0.159 MPa; when it increases from 1.5 GPa to 2 GPa, the final stress only decreases by 0.013 MPa.This is because, as the strength of the reinforced area increases, its bearing capacity correspondingly improves; at the same time, the weight of the rock mass above is also increasing, ultimately acting on the support structure.As can be seen from Figure 16b, the early stress uniformly decreases with the increase in the strength of the reinforced area.When the elastic modulus increases from 1 GPa to 2 GPa, the early stress decreases by a total of 0.337 MPa.In summary, increasing the length, thickness and strength of the reinforced area has a greater impact on the early stress than on the final stress, which is beneficial for construction stability.Setting the length to 8 m and the thickness to 1.5 m can achieve a good effect, and the strength of the reinforced area should be enhanced as much as possible within the limits allowed by construction.

Comparison between Field Monitoring and Numerical Results
The leading conduit grouting reinforcement measure in this article has been applied in actual construction and has achieved favorable results.The displacement monitors were installed every 2 m within the range of ZK10+700 to ZK10+755 during the construction process.
As shown in Figure 17, a comparison was made between the numerical results and monitoring data for the vault settlement of ZK10+727 section.Due to the idealized conditions in the numerical simulation, there are some discrepancies between the numerical results and monitoring data.The numerical results are smoother than the monitoring data, and the numerical values are generally lower.However, the numerical results and monitoring data exhibit similar trends, and the average error is only 9.12%.This validates the reliability of the numerical model and confirms the feasibility of the leading conduit grouting reinforcement measure proposed in this study.
Buildings 2024, 14, x FOR PEER REVIEW 14 of 16 In summary, increasing the length, thickness and strength of the reinforced area has a greater impact on the early stress than on the final stress, which is beneficial for construction stability.Se ing the length to 8 m and the thickness to 1.5 m can achieve a good effect, and the strength of the reinforced area should be enhanced as much as possible within the limits allowed by construction.

Comparison between Field Monitoring and Numerical Results
The leading conduit grouting reinforcement measure in this article has been applied in actual construction and has achieved favorable results.The displacement monitors were installed every 2 m within the range of ZK10+700 to ZK10+755 during the construction process.
As shown in Figure 17, a comparison was made between the numerical results and monitoring data for the vault se lement of ZK10+727 section.Due to the idealized conditions in the numerical simulation, there are some discrepancies between the numerical results and monitoring data.The numerical results are smoother than the monitoring data, and the numerical values are generally lower.However, the numerical results and monitoring data exhibit similar trends, and the average error is only 9.12%.This validates the reliability of the numerical model and confirms the feasibility of the leading conduit grouting reinforcement measure proposed in this study.

Conclusions
Based on a tunnel project in Jiangsu Province, a numerical model was established to analyze the conversion process from the DSD method to the three-bench method.The mechanical response of the support structure during the excavation method conversion was studied.An analysis was performed on potential adverse effects and a reasonable optimization measure was proposed.The conclusions are as follows: (1) When the excavation method is converted from the DCD method to the three-bench method, the early-stage deformation rate of the surrounding rock accelerates, resulting in the support structure bearing greater pressure.(2) The excavation method conversion leads to a significant increase in stress at the junction of the DCD method and the three-bench method, with the maximum stress occurring at the left arch foot.Moreover, within one excavation step after the

Conclusions
Based on a tunnel project in Jiangsu Province, a numerical model was established to analyze the conversion process from the DSD method to the three-bench method.The mechanical response of the support structure during the excavation method conversion was studied.An analysis was performed on potential adverse effects and a reasonable optimization measure was proposed.The conclusions are as follows: (1) When the excavation method is converted from the DCD method to the three-bench method, the early-stage deformation rate of the surrounding rock accelerates, resulting in the support structure bearing greater pressure.(2) The excavation method conversion leads to a significant increase in stress at the junction of the DCD method and the three-bench method, with the maximum stress occurring at the left arch foot.Moreover, within one excavation step after the conversion, the early stress at this location has already reached 77.95% of the final stress.
Therefore, special attention needs to be paid to the stability at the junction during the conversion process.(3) Setting up a leading conduit during the process of excavation method conversion can help the support structure bear part of the weight of the rock mass, improving the stress state of the support structure.Both early and final circumferential stresses are reduced under the influence of the reinforced area.This suggests that establishing a reinforced area is an effective measure to enhance the stability of the excavation method conversion process.(4) Increasing the length, thickness and strength of the reinforced area can further reduce circumferential stress, and its impact on early stress is greater than that on final stress.
Based on the parameter studies, the length can be set to 8 m, the thickness to 1.5 m, and the strength should be enhanced as much as possible within the limits allowed by construction.For different engineering projects, it is necessary to adjust the parameters of the reinforced area appropriately to achieve better reinforcement effects.
It should be noted that the geological conditions and support structures have been simplified, which may impose some limitations for this study.In future work, it is necessary to establish more accurate models to provide more precise guidance for practical engineering.Additionally, the leading conduit grouting reinforcement used in this paper is just one feasible solution.Further research should explore the effectiveness of other optimization measures to offer more options for different engineering projects.

Figure 1 .
Figure 1.Tunnel longitudinal profile and monitoring section.

Figure 6 .
Figure 6.Evolution process of vault se lement and invert uplift of the surrounding rock at ZK10+727.

Figure 6 .
Figure 6.Evolution process of vault settlement and invert uplift of the surrounding rock at ZK10+727.

Figure 8 .
Figure 8. Distribution of circumferential stress in shotcrete at ZK10+727 with conversion.

Figure 8 .
Figure 8. Distribution of circumferential stress in shotcrete at ZK10+727 with conversion.

Figure 8 .
Figure 8. Distribution of circumferential stress in shotcrete at ZK10+727 with conversion.

Figure 10 .
Figure 10.Reinforcement diagram of the leading conduit.

Figure 10 .
Figure 10.Reinforcement diagram of the leading conduit.

16 Figure 11 .
Figure 11.Finite element model of leading conduit reinforced area.

Figure 11 .
Figure 11.Finite element model of leading conduit reinforced area.

Figure 12 .
Figure 12.Distribution and evolution process of circumferential stress of the shotcrete at ZK10+727 with reinforcement.

Figure 13
Figure13illustrates the process of circumferential stress development at point.The final stress decreases from 8.98 MPa to 8.57 MPa, representing a reduction of 4.57%; early stress decreases from 7.00 MPa to 6.61 MPa, representing a reduction of 5.57%.This indicates that the optimization effect is limited under the current design, and it is necessary to further study the parameters of the reinforced area.

Figure 12 . 16 Figure 13 .
Figure 12.Distribution and evolution process of circumferential stress of the shotcrete at ZK10+727 with reinforcement.Buildings 2024, 14, x FOR PEER REVIEW 11 of 16

Figure 13 .
Figure 13.Evolution process of circumferential stress at point D.

Figure 14 .
Figure 14.The circumferential stress at point D of different lengths: (a) final stress; (b) early stress.

Figure 14 .
Figure 14.The circumferential stress at point D of different lengths: (a) final stress; (b) early stress.

Figure 15 .
Figure 15.The circumferential stress at point D of different thicknesses: (a) final stress; (b) early stress.

Figure 16 .
Figure 16.The circumferential stress at point D of different mechanical parameters: (a) final stress; (b) early stress.

Figure 15 .
Figure 15.The circumferential stress at point D of different thicknesses: (a) final stress; (b) early stress.

Figure 15 .
Figure 15.The circumferential stress at point D of different thicknesses: (a) final stress; (b) early stress.

Figure 16 .
Figure 16.The circumferential stress at point D of different mechanical parameters: (a) final stress; (b) early stress.Figure 16.The circumferential stress at point D of different mechanical parameters: (a) final stress; (b) early stress.

Figure 16 .
Figure 16.The circumferential stress at point D of different mechanical parameters: (a) final stress; (b) early stress.Figure 16.The circumferential stress at point D of different mechanical parameters: (a) final stress; (b) early stress.

Figure 17 .
Figure 17.Comparison of numerical results and monitoring data.

Figure 17 .
Figure 17.Comparison of numerical results and monitoring data.

Table 1 .
Mechanical parameters of surrounding rock and support structure.

Table 2 .
Mechanical parameters of leading conduit and the reinforced area.

Table 2 .
Mechanical parameters of leading conduit and the reinforced area.

Table 3 .
Length of reinforced area.

Table 3 .
Length of reinforced area.

Table 4 .
Thickness of reinforced area.

Table 5 .
Mechanical parameters of reinforced area.

Table 5 .
Mechanical parameters of reinforced area.