Horizontal Subzone Characteristics and Methane Seepage Properties of the Gas Flowing Fracture Zone above the Gob

Gas flowing fracture zone (GFFZ) is an active zone in the gob where pressure-relieved methane can move freely. However, there are very few research findings on the horizontal development characteristics and internal methane seepage properties of GFFZs. In this paper, based on the development height of a GFFZ above the gob, the GFFZ was horizontally divided into the following: a lateral fracture subzone, an O-ring fracture subzone, and a compacted subzone. +e identification criteria for all of these horizontal subzones were given by analyzing the influence of the stress in the coal rock mass on the development of mininginduced fractures. A numerical simulation study was conducted to determine the influences of the mining height and mining depth on the morphological development of all horizontal subzones of the GFFZ, and the simulation results showed that the mining height was the main factor influencing the development characteristics of horizontal subzones of the GFFZ. Both the maximum distance of the lateral fracture subzone beyond the mining boundary and the width of the O-ring fracture subzone increased with the increase of mining height. A physical simulation experiment was performed for extraction of gob methane through a surface borehole, and the experimental results showed that the gas flow state was laminar within the range of the lateral fracture subzone and the compacted subzone but that the gas flow state was turbulent, not following Darcy’s law, within the range of the O-ring fracture subzone.


Introduction
After a coal seam is mined, the equilibrium state of the in situ stress field is destabilized, the overlying rock is deformed and broken, primary and secondary fractures are developed, and so pressure relief occurs and permeability increases; this provides the route and channel for methane movement, creating conditions for methane extraction.According to the pressure relief, desorption, and movement characteristics of methane, overlying coal rock above the gob is divided the into a gas flowing fracture zone (GFFZ), a pressure relief and desorption zone (PRDZ), and an uneasy desorption zone (UDZ) [1,2], as shown in Figure 1.ese zones are described as follows: (1) GFFZ: Vertical penetrative fractures are present.
Relatively full pressure relief and desorption of methane occur within this zone.Methane can rush into the working face along the mining-induced crack-fracture system.Gas movement occurs mainly through free flowing along large cracks and seepage along small fractures.e height of the GFFZ is equal to that of the water flowing fracture zone (WFFZ).
(2) PRDZ: Vertical fractures are not developed.is zone is dominated by tensile fractures along the coal seam.With the transfer of mining-induced stress, methane within the PRDZ is relieved and desorbed to a certain extent and could flow along the coal seam, but little or no gas could rush down to the working face.Gas movement mainly occurs via seepage along the fracture system.e maximum development height of the PRDZ is up to the bottom of the primary key stratum.With compaction of the coal rock in the PRDZ, the desorbed methane is adsorbed again.
(3) UDZ: ere are no mining-induced fractures formed in this zone.e change in the degree of stress in the overlying rock resulting from mining is insu cient to cause desorption of methane, and the occurrence characteristics of methane basically remain unchanged.
At present, research on the horizontal development characteristics of the GFFZ and the methane seepage laws in the GFFZ by scholars from China and other countries remains relatively limited.For example, Palchik described the area that has connection with mine workings as an interconnected fracture zone [3].Computational uid dynamics simulations based on site measurement data from either the tube bundle system or tracer tests can provide detailed information for studying goaf gas ow patterns [4,5].Palchik conducted a direct measurement of methane emission from various sections of a borehole to localize horizontal fractures [6].In the 1990s, Qian and Xu studied the distribution characteristics of mining-induced fractures in the overlying rock strata, proposed an O-ring theory for the distribution of mining-induced fractures, and successfully applied their theory to the layout optimization of boreholes for extraction of pressure-relieved methane in the Taoyuan Coal Mine and the Luling Coal Mine in China [7].
In this paper, identi cation criteria for horizontal subzones above the gob are given by analyzing the in uence of variation in stress in the coal rock mass on the development of mining-induced fractures.A numerical simulation study was conducted to determine how the mining height and mining depth in uence the morphological development of the horizontal subzones based on the above criteria.A physical simulation experimental study was performed with a self-developed simulation experimental system for gob methane extraction to disclose the gas seepage laws in all subzones of the GFFZ.ese research results could provide theoretical guidance for simultaneous extraction of coal and methane.

Height of the Gas Flowing Fracture Zone
With vertical broken fractures developed, the GFFZ has a relatively high permeability, and the coal seam methane within the GFFZ is fully relieved and desorbed.e GFFZ has characteristics identical to those of a traditional water owing fracture zone, and both the gas owing fracture zone and the water owing fracture zone have vertical penetrative fractures allowing uid to ow upward and downward.
erefore, identi cation of the height of the GFFZ can be obtained by identifying the height of the water owing fracture zone.e hard rock strata, which have controlling e ects on the movement of overlying rock strata, are called key strata.A hard rock stratum that has a controlling e ect on a group of rock strata is called a subkey stratum, and a hard rock stratum that has a controlling e ect on its upper layer through surface strata is called a primary key stratum [8,9].Zhu et al. indicated that the position of a primary key stratum has a signi cant in uence on the height of the water owing fracture zone [10]; when the distance between the primary key stratum and the coal seam is <(7-10)M (where M is the mining height), the broken cracks in the primary key stratum will run through to become water owing fractures, and the water owing fracture zone would be developed at the bedrock top.When the distance between the subkey stratum and the coal seam is <(7-10)M, the broken cracks in the subkey stratum will run through to become water owing fractures, the broken cracks in the overlying rock strata, which were broken synchronously under the control of this subkey stratum, would also run through to become water owing fractures, and the water owing fractures would be developed to the nearest bottom of the subkey stratum above the critical height of (7-10) M.
Figure 2 shows how the key stratum position can be used to identify the height of the water owing fracture zone.e steps are as follows: Step 1. Collect borehole column data from the working face.
Step 2. Identify the position of the key stratum of overlying rock under the speci c borehole column conditions using the key stratum identi cation software KSPB.
Step 3. Calculate the height from the position of the key stratum to the mined coal seam, and identify whether the broken cracks in the key stratum are penetrative.If the height from the position of the key stratum to the mined coal seam is >(7-10)M, the broken cracks in this key stratum are not penetrative; if it is <(7-10)M, the broken cracks in this key stratum are penetrative, and the broken cracks in the overlying rock stratum controlled by the key stratum are also penetrative.Advances in Civil Engineering Step 4. Determine the height of the water flowing fracture zone.When the overlying primary key stratum is located within the critical height (7-10)M, water flowing fractures will be developed at the bedrock top, and the height of the water flowing fracture zone will be greater than or equal to the bedrock height; when the overlying primary key stratum is located beyond the critical height (7-10)M, water flowing fractures will be developed at the bottom of the nearest subkey stratum above the critical height (7-10)M, and the height of the water flowing fracture zone will be equal to the height from this key stratum to the mined coal seam.In this paper, the height of the GFFZ above the gob is identified using the above identification method for the height of the water flowing fracture zone based on the key stratum position.

Horizontal Subzones of the GFFZ above the Gob
3.1.Horizontal Subzones.Based on the development laws and morphological characteristics of gas flowing fractures in different areas, the GFFZ above the gob was horizontally divided into the following subzones: (I) lateral fracture subzone, (II) O-ring fracture subzone, and (III) compacted subzone (Figure 3).ese subzones are described as follows.
3.1.1.Subzone I. Lateral fracture subzone: After the coal seam is mined, the overlying coal rock stratum gets broken along the breaking line; according to the theory of an elastic foundation beam, the overlying rock stratum breaking line is located at a certain distance from the mining boundary.In addition, the coal rock masses beyond the breaking line are destroyed because of the decrease of horizontal stress, and large quantities of secondary fractures are generated in the coal rock masses under the effect of concentrated stress, greatly increasing the permeability of the coal rock masses; therefore, there exists a lateral fracture subzone with relatively developed fractures in a certain range beyond the mining boundary.erefore, mastering the morphological development characteristics in this subzone has great significance for gob methane extraction.

Subzone III.
Compacted subzone: With the constant advance of the mining working face, the overlying rock above the middle of the gob subsides, the broken and caved coal rock masses get compacted, and the mining-induced fractures inside the compacted subzone reclose; thus, the coal rock masses in this subzone have relatively low permeability.

Influence of the Variation in Stress in the Coal Rock Mass on Development of Mining-Induced Fractures.
During mining of a coal seam, the variation in stress in a coal rock mass significantly influences the deformation and failure characteristics of the coal rock mass itself.Figure 4 shows how the variation of the stress field influences the distribution of mining-induced fractures after mining of the coal seam.For the lateral fracture subzone, after the overlying coal rock above the gob is broken, the horizontal stress in the coal rock masses beyond the breaking line decreases significantly; meanwhile, under the effect of lateral abutment pressure, plastic failure occurs in the coal rock masses, significantly raising the permeability of the coal masses.For the O-ring fracture subzone, the vertical stress in the coal rock mass near the coal wall decreases under the supporting effect of the coal wall, delamination fractures can still be maintained continuously, and the permeability is relatively high.Owing to the subsidence of the overlying coal rock, the vertical stress above the middle of the gob reverts to its in situ stress state, the broken and caved coal rock is compacted, mining-induced fractures are closed, and permeability decreases.It can be discerned from the above analysis that the variation of the stress field in the coal rock mass was the intrinsic cause for the formation and development of fractures in the coal rock mass; therefore, the  stress distribution characteristics in coal rock masses after mining of a coal seam can be used as identi cation criteria for the horizontal subzones of the GFFZ used in this paper.

Identi cation Index for the Horizontal Subzone.
Figure 5 shows the distribution of vertical stress and mininginduced fractures in the overlying rock strata of the mining site after the mining of a coal seam [11].As can be seen from Figure 5, there are basically no mining-induced fractures formed outside the peak of the lateral abutment pressure, whereas there are large quantities of mining-induced fractures distributed inside the peak of the lateral abutment pressure.is can be explained as follows: inside the peak of the lateral abutment pressure, plastic failure occurs in the coal rock masses, and the fractures in the coal rock masses developed and ran through.erefore, the position of the peak of the lateral abutment pressure can be used as the identi cation criterion for the outer boundary of the lateral fracture subzone.
Qu indicated that the degree of pressure relief r 0.3 (the ratio of the vertical stress in coal after pressure relief to the original value) was used as the identi cation criterion for the boundary of the O-ring fracture subzone, based on the constitutive model of rock with primary crack deformation during triaxial pressure relief [12].It can be seen from   Advances in Civil Engineering Figure 5 that, in the middle of the gob, when the vertical stress rose to 30% of the in situ stress, the mining-induced fractures basically disappeared.

Factors Influencing Morphological Development of Horizontal Subzones
At present, scholars from China and other countries have made great strides in quantitatively studying the height of the GFFZ, but studies on the factors in uencing the morphological development of the horizontal subzones of the GFFZ are lacking.Consequently, in this paper, a numerical calculation model for dip direction was established based on the production geological conditions at the K8206 fully mechanized working face (FMWF K8206) in Yangquan Coal Mine no. 3 using the Universal Distinct Element Code (UDEC) numerical simulation software, and a numerical simulation study was conducted to determine the in uences of the mining height and mining depth of the working face on the morphological development of the horizontal subzones of the GFFZ.

Brief Introduction to UDEC.
UDEC is a twodimensional numerical software that simulates the quasistatic or dynamic response to loading of media containing multiple, intersecting joint structures.e discontinuous medium is represented as an assemblage of discrete blocks, while the discontinuities are treated as boundary conditions between blocks.Large displacements along discontinuities and rotations of blocks can occur.UDEC utilizes an explicit solution scheme that can model complex, nonlinear behaviours [13].In the coal seam, the original methane content is 7.13 m 3 /t, and the original methane pressure is 0.20-0.30MPa. e strike longwall retreat mining method and fully mechanized top-coal caving technology have been used for FMWF K8206, and the roof is managed using the fully caving method.

Production Geological Conditions at the Working
Based on the geological borehole column at FMWF K8206, the positions of the key strata at the fully mechanized working face were identi ed by using the key stratum identi cation software KSPB; the results are shown in Figure 6.As can be discerned from the gure, there are ve subkey strata and one primary key stratum existing in the overlying rock strata above FMWF K8206.
e distance from the primary key stratum to the mined coal seam is 172.6 m.According to the identi cation method for the height of the water owing fracture zone based on the position of the key stratum, the height from the position of the key stratum to the mined coal seam is >(7-10)M, so water owing fractures in FMWF K8206 would be developed to the nearest key stratum at heights >7-10 greater than the mining height; that is, in subkey stratum 3, it is 125.21 m from the mined coal seam.
Based on the production geological conditions at FMWF K8206, a 2D numerical model for dip direction was established.e model had a length of 500 m and a height of 186 m, the coal seam thickness was 7 m, and the working face length was 250 m, with a 125 m coal pillar left on each side.
e upper rock layer with a thickness of 314 m in the model was converted into a uniform load of 7.85 MPa applied to the top interface of the model.e left and right boundaries of the model were xed in the horizontal direction, and the bottom boundary was xed in the vertical direction.
Coal is usually assumed to be an isotropic homogeneous medium, and the commonly used yield criteria mainly include the Mohr-Coulomb yield criterion, the Drucker-Prager yield criterion, the Gri th yield criterion, the Hoek-Brown yield criterion, and twin shear strength criterion.Among them, the Mohr-Coulomb yield criterion can well re ect the strength characteristics of geotechnical materials, and its linear characteristics can make the procedures and equations of analytical calculation more concise.erefore, the material constitutive model used was an ideal elastoplastic constitutive model with a Mohr-Coulomb yield criterion.
is model can simulate the in uence of mining height and depth on the morphological development of horizontal subzones by enabling the thickness of the coal

Advances in Civil Engineering
seam and the top load to be changed.e structural schematic diagram of the model is shown in Figure 7.
Table 1 lists the mechanical parameters of the overlying rock strata above the gob based on the lithologic measurement results of the roof at FMWF K8206.

Influence of Mining Height on Morphological Development of Horizontal Subzones of the GFFZ.
To study the influence of the mining height on the morphological development of horizontal subzones of the GFFZ, numerical calculation models were established for simulation with mining heights of 1, 2, 3, 4, 5, 6, and 7 m, respectively, based on the production geological conditions at FMWF K8206. Figure 8 shows the laws of variation of the morphological development of the horizontal subzones of the GFFZ under different mining heights according to the simulation results.

Advances in Civil Engineering
As can be discerned from Figure 8, variations in the morphology of the lateral fracture subzone were mainly re ected by the variation of the outer boundary of the lateral fracture subzone, and the maximum distance of the lateral fracture subzone beyond the mining boundary increased with the increase of mining height.As a whole, with the increase of mining height, the average width of the O-ring fracture subzone increased gradually; when the mining height increased from 1 to 7 m, the average width of the O-ring fracture subzone increased gradually from a minimum of 22.8 m to a maximum 57.2 m; the increase gradient in the average width of the O-ring fracture subzone decreased gradually with the increase of mining height.With the increase of mining height, the average width of the compacted subzone decreased gradually; when the mining height increased from 1 to 7 m, the average width of the compacted subzone decreased gradually from a maximum of 184.6 m to a minimum of 100.1 m.

In uence of Mining Depth on Morphological Development of Horizontal Subzones of the GFFZ.
To study the in uence of mining depth on the morphological development of the horizontal subzones of the GFFZ, the morphological development characteristics of the horizontal subzones of the GFFZ were simulated at mining depths of 300, 500, and 700 m by changing the boundary load based on the production geological conditions at FMWF K8206.e vertical stress on the horizontal measuring line was extracted under the conditions in various schemes, and the laws of variation of the morphological development of the horizontal subzones of the GFFZ were obtained using the boundary identi cation criteria for the horizontal subzones of the GFFZ, as shown in Figure 9.
As can be discerned from Figure 9, with the increase of mining depth, the outer boundary line of the lateral fracture subzone tended to move away from the mining boundary, which, however, was not very conspicuous.When the mining depth increased from 300 to 700 m, the minimum and maximum distances of the lateral fracture subzone from the mining boundary increased by 1.2 and 2 m, respectively.With the increase of mining depth, the average width of the O-ring fracture subzone decreased gradually.In the process during which the mining depth increased from 300 to 700 m, the average width of the O-ring fracture subzone decreased by 10.6 m.With the increase of mining depth, the average width of the compacted subzone increased gradually.In the process during which the mining depth increased from 300 to 700 m, the average width of the compacted subzone increased gradually by 22.8 m.

Experimental Apparatus.
A self-developed "simulation experimental system for gob methane extraction" (Figure 10)

Advances in Civil Engineering
was used in this experiment, and the main components of this system included five parts: the main frame and support base, the fan power system, the methane injection system, the mining simulation system, and the monitoring and control system.e function of all parts in this system is as follows: (i) e main frame and support base: e size of the main frame is 3.2 m × 1.8 m × 1.8 m, enclosed within acrylic glass and using a sealing washer to ensure that the experimental subject is isolated from the outside air.e size of the support base is 3.3 m × 1.9 m × 1.9 m. (ii) e fan power system: is provides different suction pressures for the vertical well through the use of small vacuum pumps.(iii) e methane injection system: is supplies methane at different flow rates, from various sources such as residual coal and adjacent coal seams.(iv) e mining simulation system: A lifting screw is used to simulate coal seam mining.e caving zone at the bottom of the main frame descends by downward motion of the screw.e mining simulation system can simulate a mining height of 2-10 m.
(v) e monitoring and control system: Data from sensors in the monitoring system are converted into electrical signals by the data acquisition module and then sent to the engine control console of the control system, in which software processes and stores the experimental data.

Experimental Schemes.
In this simulation experiment, FMWF K8206 in Yangquan Coal Mine no. 3 was taken as a prototype.Loose masses and prefabricated blocks were used to simulate the overlying rocks of the caving and fracture zones, respectively; the loose masses simulating the caving zone were formed by mixing pebbles with different grain sizes and yellow sand, and the dimensions and laying sequence of the prefabricated blocks simulating the fracture zone depended on the breaking dimensions of the overlying hard rock strata (Figure 11).e surface borehole was simulated by a metal tube 3 mm in diameter.In this experiment, the surface borehole was placed 65 m from the open-off cut.In scheme 1, the surface borehole was located 10 m beyond the boundary of the gob; based on scheme 1, the borehole locations in the other schemes were moved 30 m inwards towards the gob in turn (Figure 12).Under the conditions of scheme 1, the

Analysis of Experimental Results
. Figure 13 shows the curves of the mixture flow rate at initial extraction versus suction pressure for different surface borehole locations.As can be seen from Figure 13, the initial methane flow rate from the borehole increased linearly with suction pressure under the conditions in schemes 1, 4, and 5, indicating that the lateral fracture subzone and compacted subzone had relatively low permeability and that the methane-air mixture gas had a relatively slow flow speed; the gas flow speed and the gas pressure difference had a directly proportional relationship, basically following Darcy's law; and the methane flow state in this zone was laminar.Under the conditions in schemes 2 and 3, with the increase of suction pressure, the initial methane flow rate no longer increased linearly; because the delamination fractures in the O-ring fracture subzone were developed and the gas flow speed was relatively fast, the methane flow state in the O-ring fracture subzone was turbulent under the surface borehole extraction conditions and Darcy's law was no longer applicable.

Conclusions
Based on the height of a GFFZ above the gob, the GFFZ was divided into the following: a lateral fracture subzone, an O-ring fracture subzone, and a compacted subzone.Identification criteria for the horizontal subzones were given by analyzing the influence of the variation in stress in the coal rock mass on the development of mining-induced fractures: the position of the peak of the lateral abutment pressure was used as the identification criterion for the outer boundary of the lateral fracture subzone, and the degree of stress pressure  A numerical simulation study was conducted to ascertain how the mining height and mining depth in uenced the development characteristics of the horizontal subzones of the GFFZ.e results showed that the mining height was the main factor in uencing the morphological development of horizontal subzones of the GFFZ and that the maximum distance of the lateral fracture subzone from the mining boundary increased stepwise with the increase of mining height.e average width of the O-ring fracture subzone increased gradually with the increase of mining height.In the case of the xed working face length, the width of the compacted subzone decreased constantly with the increase of mining height.
A physical simulation experiment study on gob methane extraction from surface boreholes was performed with a selfdeveloped simulation experimental system for gob methane extraction, and the simulation results showed that the extraction mixture ow rate from all surface boreholes exhibited a constantly attenuating trend; under the same suction pressure, the methane extraction mixture ow rate from a surface borehole in the O-ring fracture subzone was relatively high, whereas that from a surface borehole in the compacted subzone was relatively low.e methane ow states in di erent horizontal subzones were analyzed, and the results showed that the methane ow states in the lateral fracture subzone and compacted subzone were laminar but that in the O-ring fracture subzone did not follow Darcy's law and was turbulent.

Figure 1 :
Figure 1: e overlying three zones for pressure relief, movement, and extraction of methane in mining of a coal seam group.

7
Height of WFFZ = height from the nearest key stratum above (7-10)M to the coal seam Collect the borehole column data of the working face Is the primary key stratum located within (7-10)of the overlying key stratum Calculate the height from the position of the key stratum to the mined coal seam (

Figure 2 :
Figure 2: Prediction method flow chart for determining the height of the water flowing fracture zone based on the position of the key stratum.

Figure 4 :Figure 3 :
Figure 4: In uence of the variation in stress in coal rock masses on the development of mining-induced fractures.

4
Face.FMWF K8206 is the working face for mining #15 coal seam in Yangquan Coal Mine no. 3. Its elevation is 503.6-596.3m; the elevation of the ground surface is 1025-1168 m; and the working face has a strike length of 1579 m, a dip length of 252.2 m, and an average advance speed of 3 m/d.For the working face, the coal seam has a total thickness of 6.8 m and a net coal thickness of 6.8 m, and it dips 1 °-7 °, averaging 5 °.

Figure 5 :
Figure 5: Distribution of vertical stress and mining-induced fractures after mining of a coal seam.

FIGURE 7 :
FIGURE 7: Structural schematic diagram of the model.

Figure 8 :
Figure 8: Laws of variation of the horizontal subzones of the GFFZ under different mining heights: (a) lateral fracture subzone; (b) O-ring fracture subzone; (c) compacted subzone.

Figure 11 :
Figure 11: Schematic diagram of the trapping foam board and yellow sand.

Table 1 :
Physical and mechanical properties of rock strata.