Numerical Investigation into Evolution of Crack and Stress in Residual Coal Pillars under the Influence of Longwall Mining of the Adjacent Underlying Coal Seam

Longwall mining of the adjacent coal seam with the presence of residual coal pillars overlying the seam can result in abnormal strata pressure and severe overburden failure, which poses a significant threat to mining safety./e threat is mainly manifested in the form of intense coal or rock burst and hazardous interconnection between gobs. /is study employed the universal distinct element code (UDEC) to investigate themicroscopic failuremechanism of the overlying residual coal pillars under the influence of longwall mining of an adjacent underlying coal seam in Yuanbaowan coal mine, China. Using the Voronoi method, we innovatively visualized the evolution of cracks in residual pillars, revealed the mechanism behind the failure of pillars, and explored the evolution and distribution of abutment stress. Also, strata movement characteristics during underlying panel extraction have been surveyed. Based on the modeling results, effective measures are proposed to ensure safe mining under residual coal pillars. /is study might provide a certain reference for safe extraction of multiple seams in Datong Coalfield, China, and also in the central and western Appalachian Basin, United States, where many mining activities are carried out under residual pillars.


Introduction
Coal seams in Datong coalfield show distinctive characteristics.Most of them occur close to each other and are gently inclined shallow thick coal seams with fewer faults and hard roof.In the past, small-scale lane mining and room-and-pillar mining methods were employed to extract shallow coal seams [1].As a result, a significant number of coal pillars are left in gobs.Residual coal pillars along with irregular gobs bring a series of difficulties to longwall mining of underlying adjoining coal seams [2][3][4][5][6][7][8].Effects of multiple seam interactions and pillar instability include roof falls, rib spalling, floor heave, and bumps which can seriously disrupt mining operations and threaten the safety of miners [7][8][9][10][11][12][13][14][15][16][17].In early 2006, a West Virginia coal miner was killed by a rib roll which occurred in a high-stress zone beneath a remnant structure in an overlying mine [18].
Critical issues standing out in safe extraction of the underlying abutting coal seams under residual pillars are the effective control measures of roofs [4,5,[7][8][9][10].Many scholars have made a detailed study on the relation between overburden rock structures and support-surrounding rock systems using theoretical analysis, numerical modeling, and in situ experiments [4-8, 10-17, 19-25].Based on the catastrophe theory, Wang et al. established a cusp catastrophe model for the study of the stability of residual coal pillars left by room mining in shallow seams and explored the catastrophic failure characteristics of coal pillars [3][4][5].Using physical simulation and 3DEC numerical modeling, both movement characteristics of the overlying strata of shallow coal seams and the hydraulic support jamming mechanism in longwall mining of seams under residual pillars are investigated.Yang et al. regarded two abutting seams as a composite structure system and proposed a simplified model to investigate the stress and movement characteristics of the system [6].With the proposed model, the critical condition for roofs under residual pillars to maintain stability is derived; also, the effects of the location of immediate roof fractures and the dimension of residual pillars on roof stability are investigated.According to the loading characteristics of hydraulic support, Wang divided the strata overlying the fully mechanized coalface of a seam under a gob into four types, analyzed the caving characteristics of the strata of various types, and explored the corresponding bearing characteristics of hydraulic support [19].Ju et al. examined the failure characteristics of overburden rock structures as the coalface under pillars is progressing [20][21][22][23].
ey concluded that hydraulic support jamming is mainly attributable to the structural instability of critical rock blocks above residual pillars.Based on this, they proposed some measures against hydraulic support jamming.It can be seen that the above studies have made bold attempts to study issues concerning mining under residual pillars.To be more specific, the issues cover the failure characteristics of residual coal pillars, the stability of roof under residual coal pillars, and the weighting characteristics of roofs [2,[26][27][28][29][30][31][32].However, studies to date have mainly investigated pillar and roof failures from a macro perspective, whereas few efforts have been devoted to the microscopic mechanism behind these failures.A good microunderstanding of the evolution of stress and crack in residual coal pillars and the progressive collapse process during longwall mining of underlying seam is of great importance to make clear the relation among hydraulic supports, roofs, and pillars.erefore, to achieve safe mining under residual pillars, it is of great necessity to carry out microstudies on the progressive failure of pillars.
is study took the #6 coal seam in the Yuanbaowan Mine, China, as an example and employed the universal distinct element code (UDEC) to investigate the failure mechanism of residual coal pillars under the influence of longwall mining of the underlying adjoining seam.Above the #6 coal seam lays the #4 coal seam where residual pillars are located.Using the Voronoi method [33][34][35][36][37][38][39][40][41], we discretized meshes for the #4 coal seam and interburden between #4 seam and #6 seam, to better visualize the evolution of cracks.e blocks that can be either rigid or deformable are separated by interfaces which are viewed as contacts.With the custom-developed Fish language, we quantified the number of the tensile and shear cracks.As a result, the progressive failure of residual pillars was satisfactorily modeled.Also, we modeled the effect of residual pillars on the scope of influence and distribution of the front abutment stress in the underlying abutting seam during panel extraction.Based on the modeling results, some measures are also proposed to ensure safe mining under residual coal pillars. is study might have a certain reference value for the safe extraction of coal seams in Datong or Shuozhou Coalfield, China, where many mining activities are implemented under residual pillars.
Affiliated to China National Coal Group Corporation, the Shanxi Huayu Yuanbaowan coal mine is a consolidated mine.Mining starts from the #6 coal seam, with 6104 fully mechanized face as the first coalface.
e 6105 fully mechanized face is the successive coalface.Above 6104 and 6105 coalfaces lay in the #4 coal seam where pillars are left in place as a result of small-scale lane mining and room-andpillar mining (Figure 1).e underground site survey shows that, though apparent cracks are observed on the surface of some coal pillars, the roof of mined rooms in the #4 coal seam and most residual pillars in gobs are still in fairly good condition (Figure 2), demonstrating favorable bearing capacity. Figure 2(a) is a site photo that shows the exposed roof of mined rooms in the #4 coal seam.According to the mining practice of the nearby mines, the #4 seam could pose a significant threat to longwall mining of the #6 coal seam.e threat mainly manifests itself in the form of intense roof weighting and hazardous interconnection between the gobs, as illustrated in Figure 3.
Resulting from the extraction of the #6 seam and the failure of the residual coal pillars in the #4 seam, the intense roof weighting could exert dynamic impacts on the coalface in production, leading to rock burst and hydraulic support jamming or damage.e interconnection between the gobs results from the unfavorable caving of the roof above the coalface.Once the mined-out areas are interconnected, hazardous or poisonous gas could enter into the #6 seam coalface, posing a significant threat to personal safety.

Mining and Geological Conditions
Yuanbaowan mine is about 4.2 km long in the East-West direction and 2.1 km wide in the North-South direction.e #4 seam and #6 seam belonging to the carboniferous Taiyuan Formation are the two primary minable coal seams in the mine.Mining starts from the #6 seam, above which lays the #4 seam where residual pillars are left in place.e #6 seam is 150-160 m deep, with a thickness of 2.7-5.2m (on average 3.4 m). e #4 seam is 7.2 m thick on average, overlying the #6 seam with a vertical distance of 8.56-19.72 m (on average 16 m).Between the #6 seam and #4 seam are the coarse sandstone, fine sandstone, and sandy mudstone.Figure 4 shows the stratigraphic column of the borehole YZK2102 in 6105 longwall working face.
e 6104 face is the first coalface in the #6 seam. is coalface is 2.8 high and 180 m long, with an advance length of 1080 m. e longwall fully mechanized top caving method is employed in the coalface to extract the coal resources.e 6105 face is the successive coalface.Figure 5 is the projection plane that depicts the spatial relationship between the gob in the #4 seam (the area delineated by the red dotted line) and the coalfaces in the #6 seam.

Model Construction.
e residual coal pillars left in the #4 seam act as the carriers that connect the seam roof and floor.erefore, the stability of these coal pillars exerts a the coal pillar is 6 m high and 15 m wide.To model the most undesirable situation, we design the width of the room as 30 m in the model, which is the maximum width of the room according to our site survey.Note the longwall panel in the lower seam is excavated along the direction perpendicular to the axis of the pillars in the upper seam.In this case, the residual coal pillars could su er the maximum loading from the roof strata and exert more severe e ects on the oor.
In this study, we choose the universal distinct element code software (UDEC) as it can realistically model the crack initiation and propagation in rocks.In UDEC, the computational domain is discretized into blocks using a nite number of intersecting discontinuities.A physical discontinuity is created when the stress level at the interface between blocks exceeds a threshold value either in tension or shear.Using the polygonal block, rock failure is captured either in the form of plastic yielding of the rock matrix or displacements of the discontinuities due to Voronoi tessellation.e objective of this study is mainly to investigate the evolution of crack and stress within the residual coal pillars and failure mechanism of the pillars during longwall mining of the underlying seam.Meanwhile, considering the e ciency of model convergence, which largely depends on the size and number of blocks, we incorporate two pillars and three rooms in the upper #4 coal seam in the model.Figure 6 shows the con guration of the longwall model using the polygonal command in UDEC.e model is 250 m in length and 120 m in height.In the vertical direction, the model is constructed in strict accordance with the stratigraphic column shown in Figure 4.
As the emphasis of this study is to investigate the microscopic mechanism of failure within the residual coal pillars in the upper seam and within the interburden, we discretize both the pillars and interburden using the Voronoi method [33][34][35][36][37][38][39][40][41].Given the sizes of the whole model and the pillars and rooms within the #4 seam, we set the average edge length of the polygons within the pillars as 0.5 m, which is approximately one-tenth of the shortest side length of the pillar's cross section (6 m, the height of the coal pillar).To facilitate the study of the roof failure during the panel advance in #6 seam, we discretize the interburden between seams into polygons as well.e average edge length of these polygons is 1.5 m, which is also about tenth of the thickness of the interburden between the #4 seam and #6 seam (16.6 m). is setting can well represent crack generation  Shock and Vibration and propagation in the coal mass and can reliably reveal the failure mode of the pillar.At the same time, it could not be time consuming.e other strata, including the #6 seam, are discretized into coarse rectangular blocks with various dimensions.e bedding planes between the layers are simulated by horizontal persistent joints.Besides, to deliver more realistic modeling results, we use the CRACK and JSET command offered by UDEC to incise the polygons and blocks.In this way, we incorporate the preexisting discontinuities including bedding planes and cross joints into the pillars.e preexisting discontinuities are also assigned with suitable parameters.

Parameters Selection.
In UDEC, calibration plays a critical role, owing to inconsistency in parameters such as Young's modulus (E) and uniaxial compressive strength (UCS).To achieve this calibration, the uniaxial compressive experiment was carried out to derive the complete stressstrain curve of the specimens, using a CMT5303 apparatus (Shenzhen San Si Co., China).en, a similarly sized sample was built in UDEC with a loading rate of 0.1 m/s to calibrate the simulation parameters.With the obtained Young's modulus, the bulk modulus K and shear modulus G for the blocks in the numerical model can be calculated through the following equations [33,34,37]: where E and μ are Young's modulus and Poisson's ratio, respectively.e normal stiffness and shear stiffness for the contacts in the numerical model can be derived from the following equations [33,34,37]: where ΔZ min is the smallest width of the zone adjoining the contact in the normal direction.Previous scholars used polygonal blocks in UDEC to perform sensitivity analysis on the cohesion, internal friction angle, and tensile strength of the joints [34,42,43].ey concluded that (1) the model tensile strength is only affected by the joint tensile strength, (2) the model internal friction angle is only affected by the joint internal friction angle, and (3) the model cohesion is affected by both the joint cohesion and internal friction angle, with the former playing a dominant role.
e values of the joint cohesion, internal friction angle, and tensile strength used in this study were based on the laboratorytested specimens and adjusted according to references [34,35].
e simulation parameters for the rock mass are summarized in Table 1.e vertical stress exerted on the top of the model was calculated from the total weight of the overlying strata.e basal boundary is fixed in the vertical direction, and the lateral boundaries are fixed in the horizontal direction (Figure 6).Governed by the Mohr-Coulomb criterion, the model runs to an initial equilibrium loaded by gravity.

Excavation
Sequence.In the model, three rooms are to be excavated firstly in the #4 seam, leaving two coal pillars in place.e three rooms are excavated simultaneously due to the following two reasons: (1) the #4 seam was mined out about two years ago and (2) this study focuses on the effect of longwall mining of the underlying seam on the stability of the residual coal pillars.en, the #6 seam is to be excavated.e open-off cut of the longwall panel in the #6 seam is located right below the left edge of the Room I. e advance of the panel, from the left side to the right side of the model, is simulated using a stepwise excavation.Each stage involves a 10 m advance of the panel.For each stage, 10 5 time steps have been run to relieve stress and to allow the roof to cave.e total advance of the longwall panel is 120 m.To obtain a thorough understanding of the progressive failure of the coal pillars, we set up monitoring regions in the two pillars.A FISH function was custom developed to monitor and number the tensile and shear cracks generated within the pillars.To gain a thorough understanding of the progressive failure of the roof strata, we embed the measurement line in the immediate roof of the #6 seam to record the vertical stress.

Failure Mechanism of Residual Coal Pillar
Subject to the Influence of Longwall Mining of Underlying Seam  Shock and Vibration generation and propagation in rocks surrounding a roadway [35].ey defined debonding between contacts along the normal direction as the tensile failure and debonding between contacts along the shear direction as the shear failure.en the threshold was determined for tensile failure and shear failure, respectively.To capture the progressive failure of the residual coal pillars, we developed a FISH program to differentiate and count different kinds of failures.We set up monitoring regions in Pillar A and Pillar B to compute the changes of the tensile and shear cracks.Statistics for the number of the tensile and shear cracks versus model time steps are plotted in Figure 7.
As shown in Figure 7, following the excavation of the three rooms in the #4 seam, tensile and shear cracks in Pillar A and Pillar B firstly mushroom and then keep steady.It can be seen from Figure 7 that tensile cracks significantly outnumber the shear cracks, with the former about two times the number of the latter.As the panel in the #6 seam advances from A to B, increases of the tensile and shear cracks are not perceptible, meaning that Pillar A is mildly affected.At this point, the coalface is right below the left edge of Pillar A. As the panel progresses from B to C underneath Pillar A, cracks in both pillars start to grow.Unlike Pillar B, Pillar A is heavily affected over this course as evidenced by the significant increase in the number of tensile cracks.As the panel advances from C to D, cracks in both pillars grow much faster, especially in Pillar A. e fast growth of the cracks in Pillar A is attributable to dual influences of mining and roof caving.At this point, the panel advances past Pillar A by 5 m.en the panel continues to progress and advances from D to E. Over this course, tensile cracks keep the momentum to proliferate, while the shear cracks remain dormant, exhibiting slight growth.Different from Pillar A, Pillar B is less affected as evidenced by the insignificant increase of the tensile and shear cracks.en the panel progresses from E to F. As the panel advances under Pillar B, an interesting phenomenon appears in Figure 7. e number of tensile cracks in Pillar A fluctuates. is is caused by the closing of cracks.As for Pillar B, cracks start to grow, with the tensile cracks in particular.en the panel continues to advance.Over the course of panel advance from F to G, in Pillar A, tensile cracks still fluctuate in number, while shear cracks gain rapid growth.Pillar B over this course is severely affected as evidenced by the rapid growth of the tensile and shear crack.As the panel advances from G to H, tensile and shear cracks in Pillar A stabilize and remain steady.e reason is that mining-induced influence on Pillar A weakens.However, tensile and shear cracks in Pillar B keep growing, indicating that Pillar B is still subject to the intense mininginduced influence.

Visualization of Crack Evolution in the Residual Coal
Pillar.To visualize the development of the cracks, we handpicked some screenshots of Pillar B at various time steps, as shown in Figure 8.
To better reveal the failure mechanism of the residual pillars under the influence of longwall mining of the underlying seam, we use the FISH function offered by UDEC to visualize the evolution of tensile and shear cracks in the pillar.As mentioned above, we defined debonding between contacts along the normal direction as the tensile failure and debonding between contacts along the shear direction as the shear failure.en the threshold was determined for tensile failure and shear failure, respectively.Another FISH program was developed to identify the polygons where tensile failure or shear failure occurred and group the polygons according to different kinds of failures.en the two groups were highlighted in different colors.With this custom-developed FISH program, we visualize the tensile and shear cracks in the pillar.Figure 9 presents a more graphic illustration that sheds more light on crack evolution in Pillar B.
As shown in Figure 9(a), following the extraction of the three rooms in the #4 seam, cracks develop in the two ribs of the residual pillar.e cracks extend inward in a V-shape pattern in the left rib and in an arc-shape pattern in the right rib. ese two failure patterns are consistent with the site observation.As the panel in the #6 seam advances, cracks in the pillar gradually propagate inward.e failure depth in the rib that is close to the advancing panel (left rib) is larger than that in the other rib.Different from the right rib, tensile failure in the left rib is the dominant failure mechanism, as shown in Figures 9(b) and 9(c).As the panel

Evolution of Vertical Stress in the Residual Pillar.
To understand the evolution of vertical stress in Pillar B after the extraction of the #4 seam, we selected some stress contour graphs at various time steps, as shown in Figure 10. Figure 10(a) is the distribution of vertical stress after the extraction of the #4 seam.It can be seen that, after the removal of the #4 seam, the two ribs of Pillar B are in the plastic state.e plastic zone is about 0.7-0.8m deep.Vertical stress in the two ribs decreases and is lower than the primitive stress (3.25 MPa).e central part of the pillar is in the elastic state, where the stress ranges from 8.5-10.0MPa.Between the plastic zone and elastic zone lies the unloading zone, where the stress is highly concentrated, above 12.0 MPa.In some areas in the unloading zone, the vertical stress is up to 12.8 MPa.
erefore, the maximum stress concentration coefficient is 3.94.Overall, the vertical stress in the pillar is distributed in a saddle-shaped pattern.
As the panel in the #6 seam progresses, the unloading zone in the two sides of the pillar gradually propagates towards the center and interconnects with each other at the center of the pillar, as shown in Figures 10(b) and 10(c).It is the same with the stress concentration zones in the two ribs of the pillar.As the panel advances, the stress concentration zones gradually transfer to the center.As a result, the original elastic zone disappears.Instead, a single high-stress zone forms up, as shown in Figure 10(c).
As the advancing panel reaches the underneath of Pillar B, subject to more violent mining-induced influence, the left side of the pillar suffers more severe damage, as shown in Figure 10(c).e unloading zone transfers inward, and the high-stress zone narrows.e peak stress at the right side of the high-stress zone is 21.1 MPa (stress concentration coefficient: 6.49).As the panel continues to progress, the left side of the pillar fails completely and the high-stress zone in the pillar gradually vanishes.Stress only concentrates on the right side of the pillar, which is 12.3 MPa (stress concentration coefficient: 3.78).As shown in Figure 10(e), at this point, the unloading zone has covered most area of the pillar.Later, as the panel continues to advance, the pillar fails entirely, as shown in Figure 10(f).Shock and Vibration induced influence, the vertical stress in Pillar B increases significantly and distributes in a ladder-shaped pattern as shown in Figure 11(b4).In this stage, the peak vertical stress at two sides of the pillar moves inward and overlaps, resulting in an extensive stress concentration zone where the stress ranges from the primitive level to the maximum level 17.4 MPa.In this stage, the vertical stress in the elastic zone rises to a high level, with an average value of 15.1 MPa (stress concentration coefficient: 4.64).Different from the "saddle shape" stage, the fracture zones at two sides of the pillar moves inward, and the elastic zone narrows.less than that at the right side of the pillar.In the loadbearing zone, the maximum vertical stress reaches 19.2 MPa (stress concentration coefficient: 5.9).In this stage, the pillar has reached its maximum bearing capacity.In this stage, the pillar loses the bearing capacity completely.e vertical stress along the horizontal center line of the pillar distributes in a wavy pattern as shown in Figure 11(d4).e residual stress decreases to a level that is less than the primitive stress.
In conclusion, longwall mining of the #6 seam results in stress redistribution.As the panel advances, the elastic zone in the pillar gradually narrows and the plastic zone in the pillar acts as a carrier that transfers stress from roof strata to floor strata.As the fracture zones at two sides of the pillar move inward, the effective load-bearing zone narrows, leading to significant stress increase.In this stage, the stress in the pillar can exert an adverse impact on the safe extraction of the underlying adjoining seam under the residual pillars.Ultimately, as the panel advances past the pillar, the pillars along with the roof strata above the #6 seam fail.e pillars are only left with the residual bearing capacity.

Overlying Strata Movement of a Longwall
Mining Face in an Abutting Seam beneath a Room Mining Gob e distance between the seams is a crucial factor that affects the safe removal of the seam underlying the mined-out area or gob.For example, if the underlying seam is adjacent to the overlying mined-out area, it is highly possible that the caving zone resulting from longwall mining of the underlying seam could interconnect with the overlying mined-out area.e interconnection could lead to the rapid ingress of the gob water and gas into the coalface in production, posing a significant threat to mining safety.Also, the stability of the irregular coal pillar left in the mined-out area can have a substantial effect on the behavior of the mining pressure in the underlying panel.erefore, it is of great engineering significance to investigate the characteristics of overlying strata movement in the longwall mining of an adjoining seam underlying a gob left by room mining.12 presents the progressive failure of the strata overlying the #6 seam as the panel advances.In the initial stage, room-andpillar mining is modeled in the 3 m thick #4 seam, which results in the exposing of the immediate roof.As shown in Figure 12(a), in this stage, the roof strata maintain good integrity and only a small fraction of cracks appear.is modeling result is consistent with the site survey.

Progressive Failure of Overlying Strata. Figure
After the panel advances by 40 m, roof strata above the #4 seam still maintain good integrity, with no occurrence of caving.ough the larger area of roof strata above the #6 seam is exposed after panel advance, the roof strata are in good condition.Only a small portion of rocks caves to the gob, as shown in Figure 12(b).is is due to the reason that the immediate roof of the #6 seam is made of 7.5 m thick coarse sandstone, which possesses a higher compressive strength.erefore, it can maintain its integrity and stability over a larger span.
As shown in Figure 12(c), with the further advance of the panel, a massive amount of traverse cracks appears in the immediate roof of the #6 seam and the beam-like roof bends and fractures.As a result, a significant portion of rocks caves to the gob.28 m high above the #4 seam lays an 8 m thick siltstone, which is the key stratum [44][45][46].In this stage, bed separation occurs between this rock layer and the 3-1 seam below.
After the panel advances by 80 m, the immediate roof of the #6 seam fractures entirely and cracks begin to develop in the overlying 3.5 m thick fine sandstone and 5 m thick sandy mudstone due to deformation.At this stage, the immediate roof of the #4 seam fractures and cracks start to develop in the main roof and proliferate all the way up to the rock layer beneath the key stratum.As the roof strata above the #6 seam sag or subside, separation below the key stratum exacerbates (Figure 12(d)).
As shown in Figure 12(e), after the panel progresses by 100 m, a great number of cracks appear in the fine sandstone and sandy mudstone above the #6 seam and the overlying strata deform and subside further.Figure 12(f ) illustrates the state of the overlying strata after the panel advances by 120 m.It can be seen that the roof strata above the #6 seam cave entirely and the roof strata above the #4 seam subside as a result of the caving.It can also be seen that the separation below the key stratum extends over a larger area.

Evolution of Vertical Stress in #6
Seam Roof.In the mining of seams in proximity, residual coal pillars act as carriers that transfer stress from the overlying strata to the underlying strata, resulting in stress redistribution in the underlying seam.
erefore, the residual pillar-induced stress concentration exerts a profound influence on safe mining under pillars.To probe into the evolution of vertical stress in the roof of the #6 seam during panel extraction, we set up 16 measurement points in the immediate roof of the #6 seam to record the stress.As shown in Figure 13, at an interval of 10 m, the measurement points are distributed over a straight line that is 0.8 m away from the #6 seam.
According to the location of the stress measurement points, we plot the stress versus model time steps in three graphs as presented in Figure 14.
Figure 14(a) presents the evolution of the vertical stress recorded by the measurement points under solid coal seams (P1, P2, P15, and P16).Measurement points P1 and P2 record the rear abutment stress.e stress shows the following characteristics: (1) as the panel advances, the stress gradually increases; (2) when the panel is advancing beneath the coal pillar, the stress rises significantly; and (3) during the first roof weighting, the stress also rises rapidly.After the   Figure 14(c) presents the evolution of the vertical stress recorded by the measurement points under mined-out area (P3, P4, P8, P9, and P13).ese measurement points record the stress within the stress relief zone.After the extraction of the #4 seam, the stresses are low, ranging from 0.2 MPa to 2.2 MPa.
When the advancing panel has not reached the measurement point, the stress is slightly influenced by the front abutment stress, and the maximum stress is only about 1.9 MPa.After the panel has advanced past the measurement point, the stress reduces to about zero and undergoes slight fluctuation.14 Shock and Vibration

Measures to Guarantee Safe Longwall Mining beneath Room Mining Gob
From the above study, we obtain the following findings: (1) the scope of influence of the front abutment stress during longwall mining beneath the residual pillars is 40-50 m; (2) the front abutment stress reaches its peak about 10 m ahead of the advancing panel; and (3) under the residual pillar, the maximum front abutment stress could reach 13.2 MPa (stress concentration coefficient: 4.1).erefore, as the panel approaches and advances beneath the pillars left in the #4 seam, it is very likely that the mining pressure behaves violently.As a result, a host of problems can occur, such as hydraulic support jamming, coal pillar spalling, serve roadway deformation, and uncontrolled roof caving, posing a significant threat to mining safety.With reference of previous research achievement, we propose the following workable measures to guarantee safe longwall mining in panel 6104 and panel 6105 in Yuanbaowan mine (Figure 15).
(1) Handle the mined-out area in the #4 seam.As the #4 seam was extracted using unconventional mining method, the mined-out area and the residual pillars are unmatched in dimension.erefore, the residual pillars cannot sustain long-term stability and the presence of these pillars can cause catastrophic stress concentration.Two measures can be taken to address this problem.e first one is to force the roof strata to cave artificially.In places where a large area of roof strata are exposing, the loosening blasting method can be employed to fracture the roof and residual pillars, forcing the roof strata to cave.e second one is to backfill the mined-out area with sand or quick-setting cement.
(2) Optimize longwall mining in the #6 seam.Measures such as selecting a suitable mining height, quickening panel advance, strengthening support management, and using proper hydraulic support can bring the mining pressure under control.
e site survey shows that mined-area in the #4 seam is in a complicated situation.Besides, loosening blasting and backfilling are not only difficult to implement but also cost ineffective.erefore, optimizing the longwall mining in the #6 seam is a preferred choice.

Conclusions
Sustainable development of coal industry is critically significant for the national economy and social development.Hence, it is very important to increase the coal recovery rate, to extract coal seams with complex conditions, and to reduce mine accidents.Recent years have witnessed the increase of longwall mining in adjoining seams beneath room mining gob.Longwall mining of the abutting seam with the presence of the residual coal pillars overlying the seam can result in abnormal strata pressure and severe overburden failure, which poses a significant threat to mining safety.e threat mainly manifests itself in the form of intense roof weighting and hazardous interconnection between the gobs.erefore, it has become a significant challenge to achieve safe mining under residual pillars.To address this challenge, we took the #6 coal seam in the Yuanbaowan Mine, China, as an example and numerically investigated the failure mechanism of the residual coal pillars under the influence of longwall mining in a seam beneath room mining gob and revealed the dynamic evolution of stress and crack in residual coal pillars.We also proposed effective measures to ensure safe extraction under the residual coal pillars.e main conclusions from this study are summarized as follows: (1) A FISH function was custom developed to monitor and number the tensile and shear cracks generated within the residual pillars.It is found that, after the extraction of the #4 seam, tensile cracks significantly outnumber the shear cracks in the pillars, with the former about two times the number of the latter.As the panel advances in the #6 seam, cracks in the pillar ribs gradually propagate inward.e rib that is more adjacent to the advancing panel suffers more severe failure, with the tensile failure playing a dominant role.As the panel continues to progress, the cracks in the center of the pillar gradually coalesce with the cracks in the two ribs.e shear failure band is intermittently distributed in the residual pillar, and tensile failure is the dominant failure mechanism in the central part of the pillar.
(2) Based on the integrated analysis, we divided the vertical stress distribution along the horizontal center line of the pillar into four stages: (1) "saddle shape" stage; (2) "ladder shape" stage; (3) "oblique triangle" stage; and (4) "waveform" stage.(3) e characteristics of overlying strata movement in the longwall mining of an adjoining seam beneath a room mining gob were revealed.In the initial stage, after the extraction of the #4 seam, the roof strata above the #4 seam maintain good integrity and only a small fraction of cracks appear.e modeling results show good agreement with the site survey.After the panel advances by 40 m, roof strata above the #4 seam still maintain good integrity.ough a larger area of roof strata above the #6 seam is exposed after panel advance, the roof strata are in the favorable condition.Only a small portion of rocks caves to the gob.With the further progress of the panel, a massive amount of traverse cracks appears in the immediate roof of the #6 seam and the beam-like roof bends and fractures.After the panel advances by 80 m, the immediate roof of the #6 seam fractures entirely and cracks begin to develop in the overlying 3.5 m thick fine sandstone and 5 m thick sandy mudstone.At this stage, the immediate roof of the #4 seam fractures and cracks start to develop in the main roof and propagate all the way up to the rock layer beneath the key stratum.After the panel advances by 120 m, the roof strata above the #6 seam cave entirely and the Shock and Vibration roof strata above the #4 seam subside as a result of the caving.At this stage, separation below the key stratum extends over a larger area.(4) As the panel progresses, the rear abutment stress gradually increases.When the panel is advancing beneath the coal pillar, the stress rises significantly.
After the panel advances by over 100 m, the stress levels off.Under the influence of the pillars left in the #4 seam, the scope of influence of the front abutment stress during longwall mining is 40-50 m, and the front abutment stress reaches its peak about 10 m ahead of the advancing panel.(5) e overlap of the front abutment stress induced by mining in the #6 seam with the concentrated stress in the pillar floor can reach a maximum stress level up to 13.2 MPa (stress concentration coefficient: 4.1).erefore, as the panel approaches and advances beneath the residual pillars, it is very likely that the mining pressure behaves violently.As a result, a host of problems can occur, such as hydraulic support jamming, coal pillar spalling, serve roadway deformation, and uncontrolled roof caving, posing a significant threat to mining safety.To guarantee safe longwall mining in Yuanbaowan mine, we propose some workable measures.Shock and Vibration

Figure 2 (
b) shows the roof caving in one laneway.Subject to long-term loading, weathering, and water intrusion, the residual coal pillars dotted around the gob exhibit varying degrees of erosion and spalling (Figures 2(b)-2(d)).

Figure 6 :
Figure 6: Model geometry and boundary conditions of the model built using the polygonal command in UDEC [33].

Figure 7 :
Figure 7: Number of tensile and shear cracks versus model time steps.

4. 4 . 4 .
Fourth Stage: Waveform.In this stage, the panel in the #6 seam advances by 110 m and advances past Pillar B by 20 m.Under the combined influence of mining and roof caving, the pillar fails ultimately, with the fracture zone covering nearly the whole pillar as shown in Figure 11(d3).

12
Shock and Vibration panel progresses by over 100 m, the stress levels off, indicating that the stress has reached a new equilibrium.As for measurement points P15 and P16, they record the front abutment stress.e recorded stress shows distinctive characteristics: (1) when the panel reaches the point that is right below Pillar A, the stress slowly rises and soon stabilizes; (2) after the panel advances past Pillar B, the stress rises significantly; and (3) when the advancing panel is 20 m away from measurement point P15, the stress rises rapidly, indicating that the area now is within the scope of influence of the front abutment stress.

Figure 14 (
b) presents the evolution of the vertical stress recorded by the measurement points located in the floor of #4 seam with concentrated stress caused by the residual pillars (P5, P6, P7, P10, P11, and P12) and the boundary of the mining area (P14).estresses recorded by these measurement points can be divided into three stages based on their characteristics.In the first stage, the stresses rise due to the extraction of the #4 seam; in the second stage, the stresses increase as the advancing panel in the #6 seam is approaching; in the third stage, the stresses fall to a very low level due to the complete failure of the residual coal pillars.Just take the measurement point P11 as an example.After the extraction of the #4 seam, this point locates within the stress concentration zone caused by the presence of the residual coal pillar.e vertical stress is about 6 MPa; when the panel in the #6 seam advances by 30 m, the recorded stress starts to rise; when the panel advances by 40 m, the recorded stress rises significantly; when the panel advances by 70 m, the recorded stress rises rapidly and reaches the maximum level of 13.2 MPa (stress concentration coefficient: 4.1); as the panel advances further, the recorded stress plummets to about zero.e stress recorded by measurement point P7 shares similar characteristics with measurement point P11.erefore, according to the recorded stress, we conclude that the scope of influence of the front abutment stress during longwall mining beneath the residual pillars is 40-50 m, and the front abutment stress reaches its peak about 10 m ahead of the advancing panel.

Figure 13 :Figure 14 :
Figure 13: Stress measurement points embedded in the immediate roof of the #6 seam.

Figure 15 :
Figure 15: Block diagram of proposed measures to guarantee safe longwall mining beneath room mining gob.

Table 1 :
Simulation parameters for the rock mass and contact.