Space-Time Evolution Laws of Zonal Disintegration in Deep UndergroundCavesBased onCoordinationDeformation between the Bolt Body and Surrounding Rock

Zonal disintegration refers to the special phenomenon whereby fractured zones and intact zones appear alternately in deep-buried surrounding rocks under high stress conditions, which are different from that of the shallow rock mass. Because the divisional rupture law in engineering practice is closely related to the force characteristics of the bolt body, this paper analyzed stress distribution rules of the same bolt body at different times and that of different bolt bodies at the same time in the case of zonal disintegration based on coordination deformation between the bolt body and surrounding rock. .e nonlinear rheological mechanics model of rock mass on the elastic-plastic interface under the maximum support pressure was established. It puts forward the theoretical calculation formula about the mechanics criterion and breakdown moment of the zonal disintegration. Using the mechanics model of interaction between a bolt and the surrounding rock, the distribution locations along the bolt body of the anchor neutral points and its maximum axial forces were discussed with the multiple theoretical neutral points. Furthermore, the location and width of each fracture zone were back analyzed. .e results show that the rock mass on the elasticplastic interface of the surrounding rock has a significant creep effect after the excavation of the deep underground cave. While maximum deviator stress of the rock mass is more than its long-term strength, the rock mass will fracture along a radial direction and come into the fractured zone. .e multiple redistribution of the surrounding rock stress will generate alternate distribution phenomena of the fractured zone and intact zone. Meanwhile, the distribution regularity of the peaks and troughs interval of the displacement of surrounding rock leads to multiple neutral points along the anchor length direction. .e computed results of zonal disintegration through the back analysis can reflect the actual space-time evolution laws of zonal disintegration in deep underground caves.


Introduction
Since the twenty-first century, the development and utilization of deep rock mass engineering has been the most active research direction in the field of civil engineering. A large number of deep rock mass engineering projects, with depths of more than 1000 m, such as mine roadways, mountain traffic tunnels, water diversion tunnels, nuclear waste disposal wells, and oil combat readiness storage projects, are being built successfully. According to incomplete statistics, there have been 100 metal mines with depths of more than 1000 m built in foreign countries, predominantly in South Africa and Russia. In China, the mining depth of many mines has also exceeded 1000 m and it is expected that many mines will reach depths between 1000 and 1500 m in the next 20 years. e maximum excavation depth of Zhongnanshan highway tunnel in Qinling is approximately 1600 m and that of the right line of Liutongzhai railway tunnel, which is under construction, is approximately 1400 m. In addition, the water diversion tunnel in Jinping II hydropower station lies in 2600 m underground, and the depth of Jintan gas underground storage in Jiangsu Province is also over 1000 m depth. e occurrence environment of deep rock mass engineering is generally characterized by "three highs," namely, high ground stress, high ground temperature, and high-pore water pressure as well as strong artificial disturbance. e surrounding rock shows a significant time effect and leads to many engineering disasters such as zonal disintegration, rock burst, and extruding deformation [1,2]. Among these disasters, the zonal disintegration of surrounding rock is typical of deep rock mass engineering and has attracted great attention from experts and scholars in the field of rock mechanics engineering, both in China and abroad. In the 1980s, Shemyakin found the phenomenon of "zonal rupture" of surrounding rock by resistivity meter in deep Na>l mine, as shown in Figure 1(a); Fang et al. [3] also observed the objective natural phenomenon of zonal disintegration, as shown in Figure 1(b), through field monitoring of a deep roadway surrounding rock in Chinese Jinchuan nickel mine.
In order to understand the formation mechanism and distribution law of zonal disintegration, domestic and foreign experts have conducted a considerable amount of research in this field. In foreign countries, the authors in [4] thought that the mechanism of fracture of deep roadway surrounding rock under unloading conditions is similar to the formation of the penetrating crack of the test specimen under lateral pressure conditions; the authors in [5] studied the deformation failure of the defective medium and the zonal disintegration of the surrounding rock in a deep roadway using non-Euclidean geometry; Metlov et al. [6] analyzed the zonal rock failure of the surrounding rock based on a nonequilibrium thermodynamics view; Mirenkov [7] pointed out that the shear stress around the hole increased at the position of o times excavation radius; and Qi et al. [8] established the strain gradient equation for the zonal disintegration of surrounding rock in deep underground caves. In China, Bai [9] first discovered a very interesting engineering problem, the zonal disintegration of surrounding rock, by field monitoring the stress state of the bolt body. Subsequently, the authors in [10,11] carried out a considerable amount of pioneering research on zonal disintegration. In recent years, mechanics and instability criterion [12][13][14][15] of the zonal disintegration of surrounding rock in deep roadways have been studied and numerical simulation tests of the zonal disintegration under different working conditions have been carried out [16][17][18]. e evolution law of the zonal disintegration under the condition of different hole shapes, loading methods, and characteristic parameters of surrounding rock has also been researched through similar material model tests [11,[19][20][21][22]. In addition, the time effect problem of the zonal disintegration has been thoroughly discussed in [23][24][25].
In conclusion, despite a series of research achievements on zonal disintegration obtained using on-site monitoring, model tests, theoretical analysis, and numerical simulation, this extremely complex rock engineering problem, which is closely related to space and time, remains the subject of many research questions regarding the mechanism and spatiotemporal evolution law of zonal disintegration. e full-length anchoring bolt is widely used as an effective active support technology for the surrounding rock of a deep cave. Although the stress and strain of surrounding rock are distributed alternately in waves and troughs, as pointed out in the literature [3,20,21], and the actual stress state of the pull-pressure alternation was emphasized in [8] on the condition of zonal disintegration, research still neglects the idea that the basic law of zonal disintegration may be back analyzed according to the stress state of the bolt body. In this paper, the authors first establish the nonlinear rheological mechanical model of the rock mass on the elastoplastic interface of the surrounding rock, put forward the mechanics fracture criterion for the rock mass, and determine the fracture zone quantities as well as its formation moment. Second, the mechanical model to describe interaction between the surrounding rock and the bolt body on the condition of zonal disintegration will be established. A new approach for determining the location and width of each fracture zone is given in line with the principle and the mechanical model of coordinated deformation of bolt and surrounding rock. Finally, one engineering case study is presented to demonstrate the effectiveness, accuracy, and applicability of the proposed method for the zonal disintegration.

Analysis of the Fracture Process of Surrounding Rock in a Deep Cave
In the initial stage of the cave excavation, the second redistribution of surrounding rock stress occurs. Meanwhile, the maximum deviatoric stress originates from the difference between the tangential stress and the radial stress of the cave, which makes the deformation of surrounding rock develop into an elastic-plastic state and controls the steady state of the deep cave [19]. On account of the interior wall of the cave, a free surface, the surrounding rock produces lateral tensile expansion only towards the inside cave under the action of tangential maximum support pressure. On condition that the tensile deformation of surrounding rock reaches its ultimate strain, the first fractured zone, that is, pseudo-free surface, will arise in the surface area of the cave wall [26]. e first fractured zone appears impossibly after the surrounding rock stress is released because the level of ground stress is usually low. For a deep rock mass under high ground stress, the outer boundary of the first fractured zone generated after the stress is released is equivalent of a new excavation boundary. us, the surrounding rock stress is redistributed again. e surrounding rock stress is released again, and the second fractured zone comes into being in case the new ground stress meets the failure criteria of the rock mass [27]. us, in this way, the phenomenon of zonal disintegration will disappear until the maximum radial of the surrounding rock strain under the action of tangential maximum support pressure is lower than its ultimate tensile strain [10]. Ultimately, zonal disintegration in a deep underground cave is generated, as depicted in Figure 2.

Coordination Deformation Laws between Surrounding Rock and the Bolt Body under Conditions of Continuous
Deformation. In general, the surrounding rock is subjected to continuous elastic-plastic deformation and even failure after the excavation of the cave, and it develops into the elastoplastic zone and fractured zone accordingly. e rock mass in each deformation zone has a different corresponding radial displacement rate, which leads to the objective phenomenon that the closer to the wall of the cave, the bigger the radial displacement rate of the surrounding rock. A section of the bolt body near the excavation face of the cave tends to prevent rock mass in the fractured zone from moving to the cave, where the surface produces negative friction pointing to the cave. Because the rock mass deformation rate in the elastic-plastic zone is lower than that in the fractured zone, some positive friction pointing to the deep surrounding rock is generated on the surface of the other section of bolt body under the pullout action of the previous bolt body. e interface between the positive friction and negative friction is called the neutral point of the bolt body, on which the friction and displacement of the bolt body relative to the surrounding rock are zero and the axial tension is maximized. e displacement of surrounding rock in different depths and the internal force of the bolt body regularities of distribution are shown in Figure 3 [28].

Law of the Coordinated Deformation of the Surrounding
Rock Bolt Body in the Case of Zonal Disintegration. As mentioned previously, zonal disintegration will take place in the surrounding rock of the just-excavated cave under certain mechanical conditions. In accordance with the coordinated deformation principle of surrounding rock and the bolt body, the distribution of the friction resistance on the anchor surface varying with time and space is shown in Figure 4. e procession of zonal disintegration is a progressive repeated destruction process, which is closely associated with time and space, as can be seen in Figures 2 and 4. e coordinated deformation process between surrounding rock and the bolt body includes the following stages: (1) e first fractured zone comes into being after the rock mass near the cave surface has fractured at time t 0 . e secondary distribution of surrounding rock stress occurs afterwards. A set of reverse friction resistance is distributed on the surface of the l 0 section of the bolt body, which is in a tensile state. (2) e second fractured zone takes shape after the rock mass on the elastic-plastic interface of the secondary stress distribution has fractured at time t 1 . e third distribution of surrounding rock stress occurs soon afterwards. A set of reverse frictional resistance is distributed on the surface of the l 1 section of the bolt body, which is in a tensile state. e bolt body of this section has the tendency of elongation which extrudes the l 0 section of the bolt body and moves it towards the cave. Meanwhile, the surrounding rock around the l 0 section of the bolt body will prevent that from moving to interior of the cave. us, some positive frictional resistance is distributed on the surface of the l 0 section of the bolt body, which is in a pressed status. (3) Similarly, the third fractured zone takes shape after the rock mass on the elastic-plastic interface of the   third stress distribution has fractured at time t 2 . e fourth distribution of surrounding rock stress occurs soon afterwards. e l 2 section of the bolt body, acted on by a set of reverse frictional resistance, is in a tensile state. It is known from the above analysis that the l 2 section of the bolt body is acted on by some frictional resistance pointing to the exterior of the cave. erefore, it generates compressive deformation because it is in a pressed state. Due to the continuity of the bolt body, the l 0 section of the bolt body tends to displace towards the exterior of the cave which is prevented by the surrounding rock around the l 0 section of the bolt body and some negative frictional resistance is distributed on the surface of that bolt body, which is in a tensile status.
It is worth mentioning that Figure 4 also obviously indicates that there is always a stress state of tensionpressure alternation for the same bolt body at different times or for the different bolt bodies at the same time under the condition of zonal disintegration, which has been strongly verified by field measurements and laboratory tests. e tension-pressure alternating phenomenon of a certain section of the bolt body in surrounding rock with zonal disintegration was monitored in Quantai mine by Qi et al. [8], as shown in Figure 5. In addition, the radial strain of the surrounding rock in the deep cave is distributed with alternate peaks and troughs, varying in wave, which is also shown in Figure 6 [25].
Based on the geological mechanics model test of deep rock zonal disintegration, the authors in [20] pointed out that the radial displacement curve of the surrounding rock distributes with alternating peaks and troughs; the peak part of the surrounding rock with displacement maximum is the fractured zone, and the trough part of the surrounding rock with displacement minimum is the nonfractured zone, as shown in Figure 7.
e stress state of the full-column anchor bolt laid out in the surrounding rock presents a distribution law of the alternating tension and pressure due to the fluctuant change in the surrounding rock radial displacement curve with peaks and troughs, as is consistent with the stress distribution law of the bolt body at different moments or different sections of the bolt body, as illustrated in Figure 4.

Deviatoric Stress Maximum Calculation Acted on the Rock Mass of the Elastoplastic Interface in Surrounding
Rock. e surrounding rock of the cave is regarded as an ideal elastoplastic medium with continuous, homogeneous, and isotropic characteristics as well as remarkable creep behavior.
e deformation of the surrounding rock can be simplified as an engineering problem of plane strain, while the axial length of the cave is much larger than its transverse width.
When the (k + 1)-th fractured zone takes shape, the location under the support maximum pressure of surrounding rock is a new elastoplastic interface on which the radial horizontal stress σ p rk , tangential vertical stress σ p θk , and axial horizontal stress σ p zk correspond to the minimum principal stress σ 3 , maximum principal stress σ 1 , and intermediate principal stress σ 2 , which can be written as where k belongs to the natural number such as 0, 1, and 2. e deviator stress maximum on the elastic-plastic interface between the tangential direction and radial direction can be calculated as follows [29,30]: where R e s(k−1) is the outer diameter of the k-th fractured zone, which is equal to the (k + 1)-th excavation radius. In case of k � 0, R e s(−1) symbolizes the corresponding design excavation   Advances in Civil Engineering radius whose value is R 0 and R e pk is the (k + 1)-th plastic zone outer diameter that occurs after the (k + 1)-th fracture or excavation of the surrounding rock.
Considering the dilation and expansion of the surrounding rock cracks during the stress redistribution, the mechanical properties of the surrounding rock will deteriorate continuously, and the effective cohesion and internal friction angle can be approximated as where c 0 and φ 0 are the initial cohesive force and initial internal friction angle of original rock, respectively, and m and n are strength reduction factors, with ranges from 0.2 to 0.8. e support reaction force acted on the interior surface of the cave loaded by bolt or on the interior surface elastoplastic interface due to the expansion deformation of the fractured zone can be computed with the following formula [31]: where R i s(k−1) is the inner diameter of the k-th fractured zone of the surrounding rock and τ s is the residual shear strength of surrounding rock, which can be determined by the following formula [32]: where τ p is the maximum of the shear strength.
It is necessary to point out that the support counterforce on the interior surface of the cave or elastoplastic interface of the surrounding rock is considered to be 0 kPa for the full grouted non-prestressed anchor.

Attenuation Characteristics of Long-Term Strength of Rock
Mass on the Elastic-Plastic Interface. After the tunnel excavation, there are two cases where the rock mass failure occurs under the effect of the support pressure maximum [33]: (1) rock mass quickly breaks down in a short time, while support pressure maximum is more than the instantaneous compressive strength of the rock mass; (2) although the support pressure maximum is lower than the instantaneous strength of the rock mass, the long-term strength of the rock mass decreases continuously along with time passing and the failure of the rock mass occurs until the long-term strength is lower than the support pressure maximum. Figure 8 shows the attenuation characteristic curve of the long-term strength of rock mass with time.
According to the Maxwell rheological model and the test results [29], the relaxation equation of surrounding rock on the elastic-plastic interface can be expressed as follows: where σ c 0 and σ c t are instantaneous compressive strength and the compressive strength at t moment, respectively, and a and b are undetermined constants, whose value ranges are a � 0 ∼ 3 and b � 50 ∼ 150. According to the relationship between the principal stress maximum and the principal stress minimum on the condition of rock triaxial compression, the triaxial instantaneous limit strength can be written as follows [33]:

Advances in Civil Engineering
e compressive strength of the rock mass at any moment in case of triaxial compression could be inferred from equations (6) and (7):

Mechanical Criterion of Zonal Disintegration.
As previously mentioned, the transverse tensile deformation of the rock mass on the elastoplastic interface is only towards to the cave because of the free face on the cave wall or the exterior interface of the fractured zone. According to Griffith's strength theory, the rock mass deformation on the elastoplastic interface changes from a plastic state into a brittle radial crack state in the case that the tensile strain of the rock mass under the support pressure maximum reaches its ultimate value. e mechanics criterion of zonal disintegration can be established as follows [31]: where σ P rk is computed as where σ t represents the tensile strength of surrounding rock, which can be obtained as follows: In this section, m is a coefficient related to the type and integrity of rock mass with 0.001 ≤ m ≤ 25 and s represents complete coefficient of rock mass with 0 ≤ s ≤ 1. σ p θk and σ p zk are obtained as follows: It can be seen from the above analysis that while the radial stress, tangential stress, and axial stress on the (k + 1)th elastic-plastic interface meet equation (9), the rock mass will produce cracking damage and develop into the (k + 1)-th fractured zone under the effect of support pressure maximum.

Moment Analysis of the Zonal Disintegration.
According to the intensity theory of deviatoric stress, the rock mass on the elastic-plastic interface will transform from a plastic state into a fractured state in the case that the deviatoric stress maximum between the tangential direction and radial direction is more than the long-term strength of the rock mass. e following equation can be obtained in line with equations (2) and (8): On the basis of equation (13), the formation moment of the (k + 1)-th fractured zone on the elastic-plastic interface can be derived, as shown in the following equation: It is stressed that the compressive stress acting on the elastic-plastic interface is greater than its triaxial instantaneous limit strength, which leads to the break of the rock mass instantaneously in the case of t k ≤ 0.

Calculation Model of Bolt Neutral
Points. For the convenience of analysis, this paper assumes the following: (1) non-prestressed full-length anchor bolts are used for the initial support of surrounding rock, and their length meets computational requirements; (2) there is no relative sliding between the bolt body and the surrounding rock that is homogeneous, continuous, and isotropic; and (3) the tensile strength of the surrounding rock is far below that of the bolt body. It takes the rock mass within the scope of 1 m along the axis of the cave as a research object. e peak stress of surrounding rock transmits to the depth of surrounding rock continuously by multiple redistributions.
e new neutral points appear constantly along the bolt body during the process of zonal disintegration of surrounding rock, the fractured zone with large deformation rate, and the nonfractured zone with low deformation rate on either side of it, as shown in Figure 9.

Analysis of Anchor Neutral Position.
As can be seen from the mechanical model of the bolt body in Figure 9, under the condition of zonal disintegration, a set of frictional resistances of opposite direction and equal size acted on the bolt body within any zonal disintegration range, which is in static equilibrium. erefore, the mechanical model of the Advances in Civil Engineering (k + 1)-th bolt body can be established, as demonstrated in Figure 10.
A microunit AB of the bolt body in Figure 10 is taken as an analytic target which is at a distance of x from the free face of the cave. e following equation can be derived: According to the theory of material mechanics [34], equation (16) can be derived: where u x and σ x as well as ε x are the displacement of the bolt body relative to surrounding rock and the stress as well as the strain of the rock mass and E b is the elasticity modulus of the bolt body, respectively. e following can be obtained from equation (16): where A is the cross section area of the bolt body. erefore, by combining equations (15) and (17), the differential equation of the displacement can be derived as Based on Hookeʼs law, the relationship between shear stress and relative displacement may be expressed as follows: where K is the shear stiffness coefficient on the interface between the bolt body and surrounding rock. According to equation (19), equation (18) can be further rewritten as if Hence, equation (20) can be further rewritten as e characteristic equation of equation (22) can be shown as where λ is the characteristic root of equation (22), which can be written as It can be seen from equation (24) that the real part and imaginary part of the characteristic root are α � 0 and β � α 1 , respectively. us, the general solution of equation (22) can be expressed as where c 1 and c 2 are both undetermined coefficients. According to Figure 9, the surrounding rock displacement and the axial force on the interior end of the bolt body can be determined as Combining equations (17), (25), and (26), the undetermined coefficients c 1 and c 2 can be written as Substituting equation (27) into (25), then On the neutral point of the bolt body, the relative displacement between the surrounding rock and bolt body is zero, namely, e solution to equation (29) can be written as where k � 0, 1, and 2.

Advances in Civil Engineering
Equation (30) clearly indicates that there are a few neutral points along the bolt body in the case of zonal disintegration.

Range Determination of the Fractured Zone and Nonfractured Zone
According to Figures 9 and 10, the mechanical model of the bolt body within any partition range of the surrounding rock can be established, as shown in Figure 11. It can be seen from equation (30) and Figures 10 and 11 that the distance from any neutral point along the bolt body to the cave center can be expressed as

Calculation of Surrounding Rock Zone Width.
In engineering practice, the bolt body moves along with the radial displacement of the surrounding rock. e surrounding rock and the bolt body deform in coordination within any partition range. At any distance x from the free face of the cave, the relative displacement between the bolt and surrounding rock is defined as where u r(x) is the radial displacement of surrounding rock at any distance x from the free face of the cave and u b(x�O k M k ) is the radial displacement of the bolt body at any neutral point along the bolt body. According to the theory of rock mass mechanics, the displacement of the surrounding rock at the neutral point can be calculated as where B is a coefficient related to surrounding rock deformation, which is determined with the following formulation: where E r and μ are the elasticity modulus and Poissonʼs ratio, respectively; P 0 is the stress of the primary rock before the excavation; and R p 0 is the external diameter of the plastic zone while the secondary stress state of the surrounding rock has an elastic-plastic distribution that can be written as e surrounding rock and bolt body deform coordinately after the surrounding rock is supported by the bolt. According to the properties of the neutral point of the bolt, it can be seen that the displacement of the bolt body is equal to that of surrounding rock, namely, e displacement of the surrounding rock at optional position is expressed as where r is the distance from the optional position to the cave center. e shear force on the one side of the neutral point is equal to that on the other side of the neutral point; that is to say, the sum of the shear along the length l k of the partition bolt is zero, and therefore, Ignoring the elongation of the bolt itself and substituting equations (19) and (33) as well as (9)−(37) into equation (38) yields Equation (40) about the radius of neutral point and surrounding rock partition width can be deduced from equation (39): erefore, the parameter l k , that is, the (k + 1)-th partition range of the surrounding rock, can be inferred from equations (31) and (40). Figure 12 shows the width calculation model of the surrounding rock under the condition of zonal disintegration. ereinto, R e sk and R e pk are the relative external diameters of the (k + 1)-th fractured zone and nonfractured zone, respectively. Figure 11: Computational model of the zonal thickness of surrounding rock. 8 Advances in Civil Engineering

Analysis for Determining the Relative External Diameter of the Nonfractured Zone.
When the surrounding rock is reinforced by full-length anchor bolts, the stress acted on the nonfractured zone meets the Molar-Coulomb plastic equation, that is, where P is the component of the nonfractured zone. en, equation (42) can be inferred as where C 0 is an integral constant which is computed by equation (43): According to equations (41)-(43), the surrounding rock stress in the plastic zone can be deduced as where p k is zero for the non-prestressed full-length anchor bolt. Due to the coordinated deformation between the bolt body and surrounding rock, the displacement of the any unit along the full-length anchor bolt can be written as follows [30,35]: where G is the shear modulus of surrounding rock that is computed as e axial force of the bolt body at the free face of the cave is zero and, in the case of R 0 ≤ r ≤ ρ k , is expressed as where C is an integration constant which is determined by Substituting equations (48) into (47), we get In a similar way, under the condition of ρ k < r < r bk , the axial force of the bolt body is ascertained as follows: where r bk is the distance between the right end of the bolt and the center of the tunnel within the range of each surrounding rock partition and is calculated by the following formula: In the case of r � ρ k , equation (49) is equal to equation (50), from which the radius of the neutral point can be derived as follows: At the neutral point, the axial force of the bolt body reaches its maximum, that is, where b is a coefficient about bolt spacing and lithology and b < 1, u b 0 is the initial displacement of the free face of the cave, R pk is the relative outer diameter of the nonfractured zone, and r bk is the distance from the external end of the bolt body to the cave center.

Advances in Civil Engineering
While the surrounding rock stress redistributed (k + 2)-th, the surrounding rock turns from a nonfractured into a fractured state on the condition that the transverse tensile strain of the rock mass on the elastoplastic interface exceeds its ultimate tensile strain. en, the supporting pressure is transferred to the depth of surrounding rock. On the basis of equations (2) and (53), the relative external diameter of the (k + 1)-th nonfractured zone can be derived as follows [30]: where A k and t are determined as follows: e axial force maximum at the neutral point of the bolt body can be computed by the following equation [28]: erefore, the relative external diameter of the (k + 1)-th nonfractured zone can be obtained from equations (54)−(57).

Analysis for Determining the Relative External Diameter of the Fractured Zone.
If the stress acted on the elastoplastic interface of the surrounding rock meets the condition illustrated in equation (9), the rock mass will turn from a plastic state into a brittle fractured state and develop into the (k + 1)-th fractured zone, in which the relative external diameter is calculated with the following formulation: R e sk � R e pk 1 1 + sin φ (1/t) . (58) Substituting equations (54) into (58) leads to (59) In conclusion, the width of the (k + 1)-th fractured zone and the (k + 1)-th nonfractured zone may be evaluated as where d sk and d nsk are the width of the (k + 1)-th fractured zone and the (k + 1)-th nonfractured zone, respectively.

Case Description.
e proposed approach is applied in a water diversion tunnel, which has a total length of 81.78 km, located in the Qinling mountain area of Shaanxi Province. e depth range of the tunnel is from 500 to 2000 m, and its average buried depth is up to 912 m. Within the mileage of K45 + 180～K45 + 370, the surrounding rock is very weak and liable to large plastic deformation, rock spalling, and constriction which lead to difficulties in the tunnel construction. e design parameters, physical and mechanical parameters of surrounding rock, and anchoring parameters of the tunnel are shown in Tables 1 and 2. According to the data presented in Tables 1 and 2, each zonal width of surrounding rock in this engineering can be calculated by equations (31), (40), and (60). e radial stress, tangential stress, and axial stress on the elastic-plastic interface obtained from equations (10) and (12) are substituted into equation (9), and the stability of the elastoplastic interfacial rock mass may be discriminated. e quantity of zonal disintegration is demonstrated in Table 3. Table 3 shows that equation (9) is false in the case of k � 4. erefore, there are four fractured zones existing in the surrounding rock of the water diversion tunnel.
Based on equation (35), the external diameter of the plastic zone on the condition of elastoplastic state is 2.12 m before the surrounding rock is anchored. According to equations (54), (57), (58), and (60), the parameters such as R e p , N max , R e s , d s , and d ns can be derived. For this project, the full grouted non-prestressed anchor is adopted, for which the support counterforce is considered to be 0 kPa, and then, the parameter M, that is, the deflection stress maximum, can be calculated by equation (2). Furthermore, the creep rupture moment of the rock mass on each elastoplastic interface can be derived by equation (14), as shown in Table 4.
Combining the distribution situation of zonal disintegration in Table 3 and fractured moment listed in Table 4, the time-space distribution law of zonal disintegration about this cave is shown in Table 5 and Figure 13.

Verification.
In order to verify the reliability of the above theoretical calculation results, the TYGD10 television imager for rock drilling was used to monitor the fracture development of surrounding rock on-site. e distribution of partial typical fractures in deep surrounding rock is demonstrated in Figure 14.
e above analysis found that there are four fractured zones generated in the surrounding rock of this tunnel, with further details as follows. e first serious fractured zone can be considered as the traditional loose zone that ranges from 0.0 to 0.42 m and outside of which is the first nonfractured zone. As for the second fractured zone, the width of which is from 5.18 to 5.53 m, the damage degree of the rock mass within this range is relatively serious, and the next is the second nonfractured. e failure mode of the third fractured zone with a width range from 6.95 to 7.41 m mainly performs cracks. Last, next to the third nonfractured zone is the fissured fourth fractured zone, with a width range from 10.83 to 11.05 m. e comparison between theoretical calculation results and field measurement values of zonal disintegration is shown in Table 6.  (9)      Based on the determination of the location, thickness, and fractured moment of the zonal disintegration in this tunnel, reasonable optimization about the original support scheme as well as its relevant parameters will be discussed in another paper due to the space limitation of the present paper.

Conclusions
is study presents the applicability of revealing the spacetime evolution laws based on the mechanical model and principle of coordinated deformation of interaction between anchor and surrounding rock. Both the stress distribution rules of the same bolt body at different moments and that of different bolt bodies at the same time in the case of zonal disintegration are discussed in detail. e results of this study indicate that, with the nonlinear rheological mechanics model of rock mass on the elastic-plastic interface and the stress distribution laws of the bolt body, the rupture moment and location, the width, and the quantities of each fractured zone can be accurately predicted. Specific conclusions are listed as follows: (1) e surrounding rock stress will be redistributed repeatedly after excavation of a deep cave. Due to the remarkable rheological properties of rock mass, the long-term strength of the rock mass constantly decreases with time passing. While the tensile stress of the rock mass on the elastic-plastic interface reaches its ultimate tensile stress, the surrounding rock may come into zonal disintegration. (2) e stress and strain of surrounding rock at different depths have the distribution characteristics of tension and pressure alternation after zonal disintegration occurs. Because of the coordinated deformation between the bolt and surrounding rock, multiple groups of opposing frictional resistance and neutral points without relative displacements exist along the bolt body. According to the mechanical model of the bolt in the case of zonal disintegration and neutral point theory of the bolt body, the location and the width of each fractured zone may be analyzed inversely. (3) e zonal disintegration of surrounding rock not only has the objective natural law of fractured zone and nonfractured zone alternate distribution but also is a progressive destruction process closely related to time. e mechanics criterion and the specific time required for the rock mass on the elastic-plastic interface to take the shape of the fractured zone can be ascertained by the nonlinear rheological mechanics model under the maximum support pressure. (4) In fact, it is easy to determine the excavation scheme and its surrounding rock support parameters of a deep-buried rock mass engineering. e special zonal disintegration behavior of surrounding rock can be determined by the value of the anchor rod pretightening force. erefore, it is possible to prevent the occurrence of zonal disintegration of surrounding rock by determining reasonable bolt parameters, which is in order to provide reliable technical support for the rapid construction and safe operation for deep-buried rock mass engineering [36][37][38][39][40][41].