Breaking and Instability Movement Characteristics of High-Position Double-Layer Hard Thick Strata due to Longwall Mining

With multilayer hard roof occurring above the mining area, the strata behaviors will be significantly affected by their special characteristics and interaction. In this paper, relying on the mining background of panel 1061 of a coalmine, the breaking instability of double-layer hard thick strata (DLHTS) and the precursory information of disasters are studied by using theoretical analysis, numerical simulation, and physical simulation. Results show that the DLHTS have bearing, blocking, and controlling effects. Prior to DLHTS breaking, the horizontal development of stress concentration zone and plastic zone can be controlled, and the overburden movement is gentle, while during the breaking process, the stress drops abruptly, the plastic zone expands longitudinally rapidly, the overburden moves notably, and the energy release is sharp, which is likely to directly cause strong dynamic pressure. Moreover, the DLHTS can interact with and influence each other; that is to say, the breaking of lower magmatic rock (LMR) can provide sufficient space for that of upper magmatic rock (UMR), while the UMR breaking can promote the LMR periodic breaking in turn and cause its secondary instability, thereby easily causing complex or composite dynamic disasters’ formation. Hence, compared with single hard thick stratum occurrence, the strata behaviors affected by DLHTS breaking and instability are normal and more intense. Study results can provide a scientific basis for the prevention and control of dynamic disaster with hard thick stratum occurrence and important guidance value for promoting safe and efficient production of coalmines.


Introduction
According to data statistics, about 30% of the Chinese coalmines are covered with a hard thick roof, such as hard thick magmatic rock in Yangliu coalmine in Huaibei, thick red sandstone in Baodian coalmine in Jining, and huge thick conglomerate in Yima coalmine in Henan [1][2][3]. e highposition hard thick roof has the characteristics of high strength and large distance from the coal seam, which easily induced the overburden movement and structure, and strata behaviors' manifestation under mining disturbance has obvious variation. For instance, the large bed separation is prone to form beneath the hard thick roof, the roof will lose its stability once it is broken, thereby inducing disasters such as roof dynamic loading, hurricane, and rock bursts [4][5][6][7].
Moreover, if under the condition of coal seam group mining, it is easy to communicate the abandoned goaf at the upper part, so that the harmful gas, ponding, and hidden fire area can be connected to the working face, which can form potential safety hazards such as gas burst, water permeability, and fire [8][9][10].
Considerable studies have been conducted to study the breaking law and dynamic response of hard and thick strata with the working face mining. For instance, in terms of the breaking law of hard thick strata, a mechanical model of Winkler orthogonal beam for a hard thick stratum according to its spatial structure and boundary foundation characteristics was established, and then its breaking span was reasonably predicted [11]. In consideration of the crack propagation speed of hard roof under the feature scale of a coalmine, the evolution of the roof breaking was analyzed, then the focal model during the hard roof breaking process was built, and the far-field vibration effect was quantitatively described [12]. e sudden collapse characteristics of hard roof in steep seam were analyzed by using physical simulation and numerical simulation [13]. e breaking-movement characteristics of a hard thick stratum and the rock burst-induced mode and mechanism were studied in the vicinity of fault [1,2]. Based on the field microseismic monitoring, the fracture of the overlying double-layer hard thick roof was studied, and the long-hole presplit blasting technology was proposed to weaken the strength of the hard thick roof [3]. Based on the Hoek Brown criterion, rock mechanics test of hard roof in repeated mining was carried out, and the relationship between failure strength and confining pressure of hard roof with different damage was analyzed [14].
On the response of hard thick rock breaking, the evolution process and morphological characteristics of bed separation under the key layer were studied, and then a mechanical analysis model and prediction method of maximum bed separation were proposed [8,15]. e prediction model of overburden movement and deformation is established, and the strata behaviors of surface subsidence and crack were studied [4,16]. Based on the formation mechanism of abutment stress, a new theoretical calculation method of abutment stress with overlying multilayer hard thick strata occurrence was proposed [17,18]. In view of the mining advancing towards faults and passing through faults, the displacement-stress characteristics of hard thick stratum were studied, and then the inducing-mechanism of rockburst was revealed [1,2]. e evolution of mining fissures and abutment stress affected by the high-position main key stratum and these effects on gas migration were studied [16,19]. Based on the engineering background of panel 1411 in Huafeng coalmine, the mechanism of rockburst caused by hard roof was studied, and the characteristics of acoustic emission and other pieces of precursor information were analyzed [20]. e strength index of strata appearance was defined for the first time, a new method for grading and predicting the strength of pressure coming from the hard roof was proposed, and four quantitative indexes of different factors and levels were given [21].
Previous studies normally have analyzed the characteristics of the breaking and movement of the lower position main roof or the single-layer hard thick stratum and the causes of the dynamic disasters. In fact, when the occurrence position of hard thick stratum is high which is far away from the coal seam, the suspension above the goaf is large after mining activities, and the characteristics of stress displacement and energy release will vary significantly; moreover, when multilayer hard thick strata are overlaid, the movement of key layers affects each other and plays an important role in disturbance pressure or abutment stress. Hence, it is necessary to carry out a detailed study and analysis on the characteristics of strata behavior and its disaster induction mechanism under the condition of highposition multilayer hard thick strata.
In this paper, based on the mining background of overlying double-layer hard thick strata (DLHTS) above working face, the strata behaviors of DLHTS breaking and instability, such as overlying structure, plastic zone expansion, stress, and displacement distribution, and the precursor information-induced disasters are studied and analyzed, thereby revealing the interaction effect among the key layers and their breaking motion on the appearance of strata behaviors. Study results will provide a scientific basis for the prevention and control of dynamic disasters with hard thick strata conditions and provide an important guiding value for promoting the safe and efficient production of coalmines. During the mining period of panel 10414 of Yangliu coalmine, which has double layers (three layers in local position) of magmatic rock above coal seam marked as 10#, from December 2010 to August 2011, the statistical analyses of mine pressure observation and support resistance are carried out in order to obtain the variation curve of working resistance of the typical support (marked as 11#), as shown in Figure 1. It can be seen that panel 10414 advances less than 160 m before April (corresponding to the gray area in Figure 1), the periodic weighting interval is monitored to be about 12 m, and the average working resistance is lower than 30 MPa as a whole of which the maximum is not more than 34 MPa; then, after April (corresponding to the blue area in Figure 1), the support resistance is significantly higher than that of before April, the average value exceeds 32 MPa, and the maximum value reaches 44 MPa. e primary reason is that the overlying strata maintain stable and have little effect on the supports due to small range mining before April, causing working resistance to slightly vary; however, after April, the goaf range expands obviously, causing the overlying magmatic rock to tend to its ultimate suspension span, thereby leading to the support significantly loaded and increase of support resistance. On April 17, when the face advances to about 220 m, the large-scale suspended lower magmatic rock (LMR) breaks and moves sharply due to not bearing its own weight and overlying load, resulting in dynamic loading and large-area damage of support according to the field observations [11]. According to the analysis of the field observations, due to the mining disturbance, the fracture and bed separation develop and expand below the magmatic rock, resulting in a large amount of free gas which is rich in the 10 coal seams previously moving upwards along these cracks and gathering in the bed separation. en, the main key stratum (i.e., the upper magmatic rock (UMR)) will settle and move with the face advancing to 525 m, and the bed separation will close rapidly with squeezing the gas and water in the bed separation, leading to water-gas blowout phenomenon through 2# gas drainage hole [22].

Problems Arising and
From the aforementioned case studies of dynamic disasters of panel 10414, it can be seen that the breaking and instability movement of double layers of hard thick magmatic rock are the main cause of various dynamic disasters. Hence, it is necessary to carry out an in-depth and detailed study on the characteristics of strata behaviors and the precursor information of induced disasters with the  occurrence of hard thick strata, especially with multilayer hard thick strata, so as to explore the interaction effect among the key layers and their breaking motion on the appearance of strata behaviors. e comprehensive mechanized mining method is adopted for mining, and the free caving mode is used for roof management. is mining area is close to no. 104 mining area with similar geological conditions, so the weak burst liability identified by 104 mining area is applicable to 106 mining area.

Breaking Analysis of Overlying Magmatic Rock in Panel 1061
According to drill data of boreholes 2010-4, 9-5, and 1061-1 in panel 1061 or its vicinity, the thickness of magmatic rock and the distance between magmatic rock and coal seam vary greatly. Hence, in order to improve the accuracy and rationality of the calculation, this section takes the average value of the exposed magmatic rock thickness and the distance between magmatic rock and coal seam, as listed in Table 1. e analysis of the breaking span of magmatic rock can be divided into two aspects: the advancing direction and the strike direction of the working face. e former is to calculate the initial weighting step of the working face and the latter is to calculate the ultimate stable span when magmatic rock breaking; that is, when the span of overlying magmatic rock strata in the strike direction of the working face is large, unstable breaking may occur.

Calculation of Ultimate Stable Span of Magmatic Rock in
Strike Direction. With the panel 1061 mining, the overlying magmatic rock will subside and move with the goaf increase.
At this time, the magmatic rock can be regarded as the fixed support beam in the strike direction of the face. When it reaches the ultimate suspension span and breaks, considering the influence of overburden fracture angle [25,26], the ultimate stable span in strike direction b 0 can be expressed as where b is the ultimate span of magmatic rock in strike direction, m; k is the cracking coefficient of rock which reflects the development of fractures in rocks, k � 0.25-0.75 [25]; θ is the fracture angle of strata, θ � 75° [25,26]. LMR is an inferior key stratum; it needs to bear its selfweight and the overburden load between the two magmatic rocks, as expressed by q l � 25 kN m − 3 × 94 m � 2.35 MPa, while the UMR is the main key stratum, resulting in bearing its self-weight and all the strata load on it, as calculated by q u � 25 kN m − 3 × 384 m � 9.60 MPa. Hence, based on a double-layer magmatic rock being regarded as fixed support beams without considering the foundation effect of underlying strata, the ultimate span of two layers of magmatic rock in strike direction is obtained by substituting the above parameters into equation (1)  Due to the length of panel 1061 of 189 m, it is larger than the ultimate span of LMR and smaller than that of UMR   which indicates that the LMR will break with a certain mining distance, but the UMR may not. However, if the adjacent working face continues to be mined, the UMR will also break once the suspension length of UMR caused by two adjacent working faces mined exceeds its ultimate span. It should be noted that the parameter q in equation (1) only considers the self-weight and its follower layers when calculating LMR breaking span and does not consider the disturbance load caused by the UMR suspension [18]. Hence, the actual ultimate span should be smaller than the calculated value; that is, it is more likely to break.

Calculation of Initial Weighting
Step of Magmatic Rock. According to the comprehensive histogram of 106 mining area and thin plate theory, the magmatic rock can be regarded as thin plates because the ratio of its thickness to the strike span meets certain conditions, and then the initial weighting step of magmatic rock is obtained [25]: where a is the LMR ultimate span in advancing direction, m; L is the initial weighting step of LMR, m; H is the thickness of LMR and its follower layer, m; c is the strata density, 2500 kN/m 3 ; and α is the dip angle of the coal seam, 8°.
Combining equations (1)-(3), the initial weighting step of LMR can be obtained as L � 188.42-233.39 m. Results show that the advancing distance L is close to the length of panel 1061 (i.e., 189 m), resulting in forming the phenomenon of "one-time square" on the project site. Hence, the strata behaviors are obvious and release a large amount of elastic strain energy, which may pose a serious threat to the equipment and safety production and should strengthen the monitoring and prevention of the face.

Mining Responses before and after Magmatic Rock Breaking
According to the aforementioned analysis, the structural evolution of overlying strata and the strata behaviors may be different from those of conventional strata affected by the occurrence of hard thick strata. Hence, during the mining process of panel 1061, mining responses before and after the magmatic rock breaking are analyzed by using numerical simulation and physical simulation, in order to directly describe and quantify the effect of hard thick strata on the mining responses.

Evolution Analysis of Mining Stress and Plastic Zone.
According to the actual geological condition of panel 1061, a UDEC numerical model with the size of 700 m (length) × 450 m (height) is established to simplify several strata in order to simulate the mining state of 550 m, and then the mining stress evolution and plastic zone expansion affected by the DLHTS breaking are studied. Herein, the model boundary, loading method, and conditions are all applied, referring to the actual conditions and previous references [27][28][29][30]. Limited by the length of the paper, this section will directly show and describe the nephogram of mining stress and plastic zone before and after the DLHTS breaking.

Evolution Analysis of Plastic Zone.
Affected by the mining activities, the fractures in the overlying rock above goaf develop and expand, which shows that the plastic zone grows and expands until the overlying strata completely break and move. Figure 4 shows the distribution and evolution of plastic zone during the mining process with several typical mining states. From Figure 4, it can be seen that, with the goaf range increasing, the plastic zone is constantly expanding, and tensile, shear, or volume failures in different locations and regions occur. However, due to the block and control of the structure of magmatic rock, the development of plastic zone before and after the magmatic rock breaking shows a certain difference.
In Figures 4(a) and 4(b), when the working face advances 160 m, the fracture mainly occurs below the LMR and only shows horizontal expansion due to block effect of LMR, leading to plastic zone gathering LMR beneath and presenting much tensile, shear, or volume failure in the model; at this moment, LMR is in the state of energy accumulation and stress rises due to large-scale suspension, which provides sufficient energy conditions for the subsequent dynamic disasters. As the face advancing to 190 m shown in Figure 4(c), fractures develop in LMR caused by its bending and deformation because its suspension reaches its ultimate span and drives the upper strata to move harmoniously, resulting in the plastic zone obviously extending to the bottom of UMR, but the LMR can still maintain relative stability. As shown in Figure 4(d), the fracture in LMR continues to develop and the obvious tensile fracture through LMR occurs with the face advancing to 220 m, and the range of plastic zone expands obviously to the bottom of UMR, which indicates that the LMR has been completely broken and drives the upper layer to move harmoniously. At this time, attention should be paid to the energy released by the breaking and violent movement of magmatic rock and its induced harm. Meanwhile, the simulation results of LMR breaking (220 m) are basically consistent with the theoretical analysis results (188.42-233.39 m).
As shown in Figure 4(e), we can see that a new plastic failure zone through the whole LMR is formed which shows that LMR breaks periodically and drives its upper layer to move synchronously which lead to an obvious vertical tensile fracture in UMR, indicating that the UMR begins to bend and break but still has enough bearing capacity. In  Shock and Vibration initial and periodical break (reflected by multiple vertical plastic zones running through the UMR); it indicates that the UMR is basically completely breaking and we should pay attention to the impact of UMR breaking and its secondary instability effect on LMR.

Evolution Analysis of Mining Stress.
Mining activities break the original stress balance and produce stress concentration or stress reduction in surrounding rocks, which has a certain influence on mining face, roadway layout, coal pillar reservation, and so on. In particular, high-stress concentration is likely formed with the occurrence of hard thick stratum because of the characteristics of mechanical and structure. Figure 5 shows the nephogram of stress distribution and evolution with typical mining states. It can be seen that, due to the occurrence of two layers of magmatic rock, the influence range and value of front abutment stress before breaking will continue to increase, while after breaking the energy will be released at the same time. Combined with the distribution of the plastic zone in Section 4.1.1, the variation of tensile stress in magmatic rock is analyzed as an example. During the mining process of 190 m, the LMR continue to suspend, affected by mining activities and the abutment stress rises, and the LMR above the goaf gradually runs into the stress reduction zone; then, the LMR is bent and deformed with the face advancing to 190 m, resulting in a large range of tensile failure zone, and the maximum tensile stress reaches 9 MPa, as shown in Figure 4(c) and 5(c). When it continues to advance to 220 m, the LMR appears a tensile fracture throughout the whole layer (Figure 4(d)) which presents that it has been broken, the maximum tensile stress can reach 13 MPa ( Figure 5(d)), and it drives the upper layer to move synchronously. en, the periodic breaking of LMR drives the synchronous movement of its upper strata with the face mining 250 m so that the UMR suspension increases and the tensile fracture appears in UMR (Figure 4(e)); at present, the maximum tensile stress is only 8 MPa ( Figure 5(e)) and the UMR still has enough bearing capacity. When the face advances to 370 m, the plastic failure zone in UMR runs through the whole layer and presents the initial breaking, and the maximum tensile stress rises to 12 MPa (Figure 5(f )); however, the value becomes 1.1 MPa (Figure 5(g)) when the face advances to 550 m.
As shown in Figure 6, the variation of tensile stress before and after double layers of magmatic rock breaking can be divided into three stages. Stage a corresponds to the mining process of the first 220 m; the LMR is continuously suspended, bent, and subsided, and the tensile stress consequently rises with the maximum of 13 MPa; this process shows that stress concentration and energy accumulation are affected by LMR large-scale suspension, until the tensile stress of LMR exceeds its ultimate strength and it breaks completely with advancing to 220 m (as shown in Figure 4(d)). Subsequently, the tensile stress decreases sharply, and the high concentration stress decreases and the energy is released; however, the UMR begins to span, bend, and settle continuously because of LMR periodic breaking, and the tensile stress increases and the maximum recovers to 12 MPa; this process shows stress concentration and energy accumulation caused by UMR suspension and corresponds to the stage b (advancing to 370 m). As the face continues to advance, the UMR experiences the initial and periodic break, leading the tensile stress to drop suddenly, and the elastic strain is released, corresponding to the stage c in Figure 6.

Characteristics of Structure and Movement of Double
Layers of Magmatic Rock. In order to show the structure and movement characteristics of hard thick strata more intuitively, taking the double layers of magmatic rock of Yangliu coalmine as the prototype, then a physical model test is carried out based on similarity theory with the model size of 3.0 m × 0.4 m × 1.5 m (length × width × height), and the similarity coefficients are 1/200 (geometry), 1/300(strength), and 1/1.5 (density). In the model test, river sand, calcium carbonate, and gypsum are chosen as equivalent materials, and the actual strata strength is needed to be converted into the simulated strength of equivalent material according to the similarity coefficients; then, the optimum material ratio of each stratum in model is determined through the literature reference and multiple comparative tests [31][32][33]. In addition, the designed mining height and depth in model are 8 m and 600 m and can simulate to advance to 500 m; moreover, some iron needs to be imposed on top of the model to compensate for the upper failed strata to simulate [34]. Record the structure and movement of overburden by photography.

Evolution of Overburden Structure.
During the mining process, overlying strata will collapse and move, but it is significantly affected by the structure and mechanical characteristics of hard thick strata. Hence, several typical mining states are chosen to study the structural shape and evolution of overburden under the effect of hard and thick strata, as shown in Figure 7.
It can be seen that double layers of magmatic rock play a decisive role in the development of bed separation and structural evolution of the overburden, and they interact with and influence each other. In Figure 7(a) and 7(b), during the mining process of 220 m, the overlying strata are moved, and the bed separation develops to the bottom of LMR and expands laterally, thus resulting in LMR beginning to bending and subsidence with obvious tensile failure occurring with advancing to 220 m; the results have good agreement with the development law of plastic zone in Figure 4(a)-4(c), and the bed separation below provides enough space for LMR moving and breaking. When the face advances to 236 m, as shown in Figure 7(c), the LMR makes initial breaking and instability, driving its upper strata to move synchronously because of the fractures development at both ends of LMR and then running through the whole layer; the LMR and its follower layer move significantly, and the bed separation develops rapidly to the bottom of UMR, which is similar to the description results of Figure 4(d). As the face continues to advance to 330 m (in Figure 7(d)), the periodic breaking of LMR increases the suspension span of UMR and provides enough space for UMR breaking; and then the UMR and its follower layer subside and move with the face advancing to 410 m, resulting in the top of the model presenting obvious displacement which verifies that the UMR is the main key stratum and can control all strata above it; the results are consistent with the development law of plastic zone in Figure 4(e). When the face continues to advance to 470 m, several times periodic breaking of LMR appears and UMR suspension span increases, causing the UMR initial breaks and all follower layers move synchronously and significantly, which is similar to the description in Figure 4(

Movement Analysis of Double Layers of Magmatic
Rock. In order to better quantitatively describe the movement of double layers of magmatic rock, two displacement monitoring lines are set at the bottom of two layers of magmatic rock in the physical model, and these specific locations are calibrated with the red-dotted line in Figure 8(a) and 8(b); then, a reflector is fixed with 10 cm interval, and the total station without prism is used for displacement measurement. Figure 8 shows the movement variations of the two layers of magmatic rock before and after their breaking. From Figures 7(a)-7(c) and 8(a), we can see that when the working face advances to 220 m, the LMR reaches ultimate span and begins to bend and settle with small displacement; then, it advances to 236 m, and LMR suddenly breaks and loses stability, causing its displacement to sharply increase with the maximum value of 5.50 m which is likely to release a large amount of elastic strain energy suddenly, thus inducing dynamic disasters, such as rockburst and shock bump.
As shown in Figures 7(d)-7(f ) and 8(b), when the working face continues to advance to 330 m, the UMR shows slight deformation and settlement affected by the initial and periodic breaking of LMR; while it advances to 340 m, due to the heterogeneity characteristic and weak plane occurrence in magmatic rock, the suspended part of UMR is sheared and destroyed in advance, which presents direct collapse and settlement with the maximum displacement of 5.0 m; it is easy to cause such dynamic problems as water-gas inrush in the bed separation and necessary to strengthen the monitoring and protection brought by the weak plane of rock mass. en, the UMR shows obvious bending and settlement with the face advancing of 410 m, while it breaks and loses stability while advancing to 470 m and its maximum displacement rises to 5.6 m, thus easily resulting in the release of elastic strain energy and inducing dynamic disasters. Hence, with the hard thick strata occurrence above the working face, it is easy to cause a sudden and sharp move of hard thick rock, and the displacement presents mutation characteristic. Figure 9 shows the movement states of LMR before and after UMR breaking. It can be seen that, before UMR completely breaks (corresponding to advancing 410 m), the LMR has experienced multiple periodic breaking and after a long time of settlement and compaction, its displacement is basically stable; however, the LMR is still in a cantilever state before next periodic breaking at this time. After UMR complete breaking (corresponding to advancing to 470 m), the movement range and displacement of LMR increase significantly as depicted in Figure 9 with blue-dotted box, and the displacement of LMR that has moved previously also increases slightly; it shows that the breaking and movement of UMR not only can promote the periodic breaking of LMR but also may cause the secondary impact on its stability.
In conclusion, the overlying high-position magmatic rocks play a role of bearing, blocking, and controlling, and they have a significant influence on the mining responses (i.e., overlying rock structure, plastic zone expansion, stress, and displacement distribution) before and after they break. Prior to DLHTS breaking, the horizontal development of stress concentration zone and the plastic zone is controlled, and the movement of overburden is gentle, while during the breaking process, the stress drops abruptly, the plastic area expands longitudinally rapidly, the structure of the magmatic rock is unstable, the displacement of overburden increases obviously, and the energy release is sharp. Meanwhile, the DLHTS can interact and influence each other, that is, the LMR breaking can provide sufficient space for UMR movement, while the UMR breaking can promote the LMR periodic breaking in turn and cause its secondary instability. Hence, when DLHTS occurs above the working face, the direct impact caused by DLHTS breaking and movement should be paid more attention to and also the LMR periodic breaking and the secondary instability caused by the UMR breaking. Generally speaking, compared with the influence of single-layer hard thick strata breaking, the  strata behaviors produced by DLHTS breaking and instability are much stronger, and it is easier to induce strong dynamic pressure phenomena, such as strong support dynamic load, rockburst, shock bump, and coal and gas outburst.

Engineering Case
In this section, the support pressure in panel 1061 will be monitored in order to further explain the influence of breaking and instability of hard thick strata on the appearance of strata behaviors. In panel 1061, 12 sets of pressure monitors for fully mechanized support are installed, which are, respectively, arranged on hydraulic supports marked as 5#, 15#, 25#, 55#, 65#, 75#, 85#, 95#, 105#, and 115# from the head to the tail. Figure 10 shows the variation curve of typical support pressure monitored in June 2014. According to Figure 10 depiction and the information provided by the coalmine, the pressure of most supports increased after June 18, among which the pressure of 5, 15, 25, 55, 65, 75, and 115# supports increased significantly, especially after June 23, the pressure of some supports has exceeded 40 MPa, and the maximum of 65# support can reach 48.4 MPa. By referring to the field monitoring and analysis, it is found that initial breaking happens to the highposition hard thick magmatic rock after June 23, and the magmatic rock lost its stability and moved violently, causing the hydraulic support to press in a large range and dynamic loading phenomenon obvious; at present, the dynamic pressure coefficient of the supports is more than 2.0 in 15, 25, 55, 95, and 105# supports, and maximum value has reached 2.58 in support 95#, as shown in Table 2. e results show that the high-position hard thick magmatic rock loses stability and moves rapidly after breaking, resulting in producing strong dynamic impact and releasing a large amount of elastic strain energy; this conclusion has good agreement with the results of Sections 3 and 4. en, the energy transmits to the surrounding rock in the form of shock wave, which causes the strata behavior of strong dynamic pressure of the supports.
From the case studies of panel 10414 in Section 2.1, the breaking and instability of double layers of hard thick magmatic rock are the main cause of various dynamic disasters, which seriously threatens the safe production during coal mining. Meanwhile, the analysis of Section 4 shows that when DLHTS occur above the working face, the hard thick stratum can interact and influence each other, which is likely to cause complex or composite dynamic disasters, such as  rock bursts, water inrush from bed separation, surface subsidence, and crack. Hence, compared with the influence of single-layer hard thick strata breaking, the strata behaviors produced by DLHTS breaking and instability are much stronger, and it is easier to induce strong dynamic pressure phenomena.

Conclusions
In this paper, based on the mechanical and structural characteristics of double-layer hard thick strata and their interaction effect, the studies for breaking instability of DLHTS on strata behaviors and the precursory information characteristics of the disasters are carried out by multiscale and multimethods. We can obtain several conclusions as follows: (i) e high-position DLHTS have bearing, blocking, and controlling effects. Prior to DLHTS breaking, the horizontal development of stress concentration zone and plastic zone can be controlled, and the overburden movement is gentle, while during the breaking process, the stress drops abruptly, the plastic zone expands longitudinally rapidly, the overburden moves obviously, and the energy release is sharp, which is easy to directly cause the appearance of strong dynamic pressure. (ii) e MLHTS can interact with and influence each other. e LMR breaking can provide sufficient space for UMR movement, while the UMR breaking can promote the LMR periodic breaking in turn and cause its secondary instability, which is easy to form complex or composite dynamic disasters.
Hence, when DLHTS occurs above the working face, the direct impact caused by DLHTS breaking and movement should be paid more attention to and also the LMR periodic breaking and the secondary instability caused by the UMR breaking. Compared with the influence of single-layer hard thick strata breaking, the strata behaviors produced by  multilayer hard thick strata breaking and instability are much stronger. Study results can provide a scientific basis for the prevention and control of dynamic disaster with multilayer hard thick strata occurrence and provide important guidance value for promoting the safe and efficient production of coalmines.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare no conflicts of interest.