Control Effects of Five Common Solid Waste Backfilling ... - MDPI

0 downloads 0 Views 9MB Size Report
Jan 2, 2019 - In one SBM mine site or digging, however, different backfilling ... Figure 1. The movement of overburden rock layer in Solid Backfill Mining ...
energies Article

Control Effects of Five Common Solid Waste Backfilling Materials on In Situ Strata of Gob Peng Huang 1 , Sam Spearing 2, *, Feng Ju 3, *, Kashi Vishwanath Jessu 2 , Zhongwei Wang 1 and Pai Ning 3 1 2 3

*

School of Mines, Key Laboratory of Deep Coal Resource Mining, China University of Mining and Technology, Xuzhou 221116, China; [email protected] (P.H.); [email protected] (Z.W.) Western Australia School of Mines, Curtin University, Kalgoorlie, WA 6430, Australia; [email protected] State Key Laboratory for Geomechanics and Deep Underground Engineering, China University of Mining and Technology, Xuzhou 221116, China; [email protected] Correspondence: [email protected] (S.S.); [email protected] (F.J.)

Received: 18 September 2018; Accepted: 10 October 2018; Published: 2 January 2019

 

Abstract: Solid backfill mining as a green mining method has already been successfully applied in many mine sites. Higher requirements for the backfilling materials have been put forward in special regions, such as shallow coal seams, ecologically susceptible areas, and sites with building on the surface. The control effects of common backfilling materials on in situ strata of gob need to be studied and compared to ensure the suitable materials are applied in the mine. The meso-structure, stress variation, energy dissipation, and backfilling effects of the five common solid backfilling materials, which are Aeolian sand, gangue, mineral waste residue, coal ash, and loess are analyzed in this paper. The results show that the Aeolian sand and gangue are densely packed and internally hard when compared to other backfill materials. The deformation of the five materials to absorb the same amount of energy in the order of higher to lower was determined as Aeolian sand, gangue, mineral waste residue, coal ash and loess. The Aeolian sand’s strain energy density is 1.67 times larger than the loess, however, the Aeolian’s strain energy density in front of working face is just 32.2% of the loess, which is important to ensure the safety of the working face. The stress changes were monitored in situ with Aeolian sand and gangue as backfill materials. The monitoring results show that the Aeolian sand stress increase rate is quicker than gangue, which can support the roof effectively. Keywords: green mining; solid backfilling material; energy dissipation; compaction property; meso-structure

1. Introduction Solid Backfill Mining (SBM) as a green mining method has already been applied successfully in many coal mines in China. Backfilling materials, such as gangue, coal ash, loess, etc., were delivered to fill the empty area in SBM [1–3]. In one SBM mine site or digging, however, different backfilling material perform different compressive properties. The backfilling material selection not only considers the convenience of obtaining and the cost but also pays attention to its backfilling effect, especially in shallow depth or thin bedrock or high mine height coal mine site. For example, Aeolian sand and loess cost less and are easy to obtain if used as backfilling material in the west of China, due to lots of those material covering the surface. In that region, meanwhile, the thick and shallow coal seam, thin bedrock, and fragile ecological environment all need the high resistance requirement of backfilling material [4–7]. The key to applying SBM successfully is its compaction property, which affects the rock layer deformation, fracture, and energy release [8,9]. During compaction, the backfill material’s porosity,

Energies 2019, 12, 154; doi:10.3390/en12010154

www.mdpi.com/journal/energies

Energies 2019, 12, 154

2 of 14

deformation, stress variation, and energy absorption are important factors affecting the backfill effect [10]. In the last decades, many theoretical studies evolved on the gangue deformation during the compaction test. Wu [11] studied solid backfilling materials compaction characteristics from stress and strain relationship. Huang et al. [12] researched the strata movement control with backfilling the gangue and coal ash mixture into gob. Liu [13] studied the gangue compaction mechanical property in SBM. Actually, the horizontal distortion is almost unconfined in a mine site [14,15]. The deformation being considered solely as an influence factor is insufficient. Therefore, the material stress and strain energy during compaction are put forward as another index to indicate the backfilling effect. Little research has been conducted on energy release of backfilling materials. Zhou et al. [16] tested energy dissipation of gangue in the loading process under conditions of different particle sizes, loading rates, and first-time stress loads. Less research focuses on or compares with the five common solid backfilling materials’ stress. Thus, this paper compares five different types of solid backfill in regard to meso-structure by using scanning electron microscopy, compaction properties, stress distribution, and energy dissipation. Results from monitoring of stress variation with two different backfill materials from the mine site are discussed. 2. Basic Principles of SBM and Influence Factor 2.1. Basic Principles of SBM In SBM, solid waste materials, such as Aeolian sand, gangue, coal ash, mineral waste residue, and loess, are transported through a vertical pipe and then delivered to the backfilling area with the belt conveyor and the self-moving loading mechanism. The unload hole is well-organized to unload the backfilling materials to the empty area and then using the punning structure to make the backfilling materials denser. When compared with the conventional face layout, a belt conveyor in the tailgate in an SBM face delivers the filling materials to a conveyor in the gob side face behind the shield supports [17]. Therefore, in SBM, the face layout allows simultaneous mining and backfilling operations. 2.2. Key Influence Factor in SBM With the advancing working face, the roof begins to sink and the backfilling materials are compacted progressively. The backfilling material plays a very important role that supports the roof in the gob, which requires the backfilling material to undergo small deformation during compaction and absorb more energy from the roof [18]. The ideal result with the backfill that needs to be achieved is to form a fissure and bending zones in the overburden rock layers to the limit where the formation of a caving zone is initiated [19], as shown in Figure 1. Therefore, different materials have a different effect on the cracks and deformation of the overburden rock layer in the gob.

With the advancing working face, the roof begins to sink and the backfilling materials are compacted progressively. The backfilling material plays a very important role that supports the roof in the gob, which requires the backfilling material to undergo small deformation during compaction and absorb more energy from the roof [18]. The ideal result with the backfill that needs to be achieved is to form a fissure and bending zones in the overburden rock layers to the limit where the formation Energies 2019, 12, 154 3 of 14 of a caving zone is initiated [19], as shown in Figure 1. Therefore, different materials have a different effect on the cracks and deformation of the overburden rock layer in the gob.

Figure 1. The movement of overburden rock layer in Solid Backfill Mining (SBM).

The five common solid backfilling materials, which are Aeolian sand, gangue, mineral waste residue, coal ash, and loess, are distributed in different ways. In Eastern China, lots of gangue is produced from mine sites; in the north of Chian, the loess and Aeolian sand are covered in the surface, which means the materials are cheap; in the western area of China, lots of mineral waste residue are produced from the open pit. Considering the cost of backfilling materials and the convenience of obtaining them, the common backfilling materials compaction property needs to be studied. 3. Backfilling Material Meso-Structure 3.1. Source of Backfill Materials The Aeolian sand was selected from the Yuheng diggings, Shaanxi province. The washing gangue size was filtered under 30 mm from Dongpin coal mine, Shanxi province. The coal ash was selected from Yangzhuang coal mine power plant, Huaibei city, and then mixed with the water to make its moisture content reach 21%. The loess and mineral waste residue are also from Dong Ping coal mine. The mineral waste residue size was sieved under 50 mm. 3.2. Scanning Electron Microscope Test Five backfilling materials’ mesostructure were analyzed by using an FEI Quanta TM 250 scanning electron microscope, as shown in Figure 2.

make its moisture content reach 21%. The loess and mineral waste residue are also from Dong Ping coal mine. The mineral waste residue size was sieved under 50 mm. 3.2. Scanning Electron Microscope Test Five Energies 2019,backfilling 12, 154

materials’ mesostructure were analyzed by using an FEI Quanta TM4 of 250 14 scanning electron microscope, as shown in Figure 2.

Figure 2. SEM micrographs of five materials.

As shown in Figure 2, Aeolian sand particles are very small, uniform, dense, and hard, which means efficient pressure transfer among particles. Coal ash is fluffy internally and the particles have an irregularly flaky structure. The gangue particles are dense, small particles filled with few large particles. The mineral waste residue has a large portation of bulky mineral. The loess has developed intergranular mesopores and micropores, with good connectivity between them, and loose structure. 4. Backfilling Material Compaction Property and Energy Dissipation Test 4.1. Test Equipment and Scheme The YAS-5000 electro-hydraulic servo-controlled rock mechanic test system, manufactured by Changchun Kexin Test Instrument Company, was employed for the backfill material compaction tests. The literature [20] gives detailed information about steel circular cylinder for testing. The uniaxial compaction tests were performed on five selected backfilling materials by using the steel circular cylinder, for 0–10 MPa stress in the axial direction, with a loading rate of 1 kN/s. 4.2. Test Result The stress-strain relationships of the five backfilling materials under compaction test were derived, as shown in Figure 3.

uniaxial compaction tests were performed on five selected backfilling materials by using the steel circular cylinder, for 0–10 MPa stress in the axial direction, with a loading rate of 1 kN/s. 4.2. Test Result Energies 2019, 12, 154

5 of 14

The stress-strain relationships of the five backfilling materials under compaction test were derived, as shown in Figure 3. Test datas: Aeolian Sand Fitting curves:

Gangue

Mineral Waste Residue

Coal Ash

Loess

10 σ e = 0.143e(18.280ε

e − 2.920)

2

R = 0.999

σ d = 0.296e(16.825ε

8

d

−1.808)

σ c = 0.417e(17.713ε

6

c −1.642)

2

R = 0.999

σ b = 0.5221e(17.304ε

b − 0.939)

Initial stage Developing stage

Stress/MPa

R 2 = 0.999

2

R = 0.998

4

σ a = 0.668e(47.787ε

a − 0.817)

2

R = 0.986

2

0 0.0

0.1

0.2

0.3

0.4

Strain Figure 3. 3. Stress-strain Stress-strain relationship relationship of of the the five five backfilling backfilling materials. Figure materials.

As shown shown in in Figure Figure 3, 3, the of five As the stress-strain stress-strain curves curves of five backfilling backfilling materials materials show show an an exponential exponential distribution trend. In the initial stage (0–2 MPa), the strain increase is larger and the deformation rate distribution trend. In the initial stage (0–2 MPa), the strain increase is larger and the deformation rate is faster, which accounts for 70–80% of the total strain. With the gradual increase of pressure, the filling is faster, which accounts for 70–80% of the total strain. With the gradual increase of pressure, the material is gradually compacted and the strain increment decreases gradually. The compaction stress filling material is gradually compacted and the strain increment decreases gradually. The compaction state of the five materials from high to low is Aeolian sand, gangue, mineral waste residue, coal ash, stress state of the five materials from high to low is Aeolian sand, gangue, mineral waste residue, coal and loess. ash, and loess. The strain ofof material in in thethe process of compaction is given as [16]: The strain energy energydensity densityofofunit unitvolume volume material process of compaction is given as

[16]: υε = (1 −

µλh ) r

Z ε 1 0

σdε

(1)

where µ is the lateral pressure coefficient, λ is the friction coefficient of the filling material and the compact cylinder, σ is the compaction stress (MPa), r is the steel circular cylinder’s inner diameter (mm), and h is the filling height (m). The strain energy density-strain curve for the five backfilling materials have been calculated and shown in Figure 4.



r

0

where μ is the lateral pressure coefficient, λ is the friction coefficient of the filling material and the compact cylinder, σ is the compaction stress (MPa), r is the steel circular cylinder’s inner diameter (mm), and h is the filling height (m). Energies 2019, 12, 154 of 14 The strain energy density-strain curve for the five backfilling materials have been calculated6 and shown in Figure 4. 500

Test datas: Aeolian Sand Fitting curves:

Gangue

Mineral Waste Residue

ν e = 6.632e(18.316ε

Coal Ash

Loess

e − 2.923)

ν d = 13.921e(17.014ε

400

d

−1.791)

2

R = 0.999

ν c = 18.292e(18.114ε

c −1.639)

2

R = 0.999

300

ν b = 20.198e(18.265ε

b − 0.921)

2

R = 0.999

200

ν a = 9.829e(51.342ε

Initial stage Developing stage

Strain energy density/kJ•m

-3

R2 = 1

a − 0.904)

2

R = 0.999

100

0 0.0

0.1

0.2

0.3

0.4

Strain Figure4. 4.Strain Strainenergy energydensity-strain density-strain relationship relationship of of the the five five backfilling backfilling materials. materials. Figure

As shown in Figure 4, the strain energy density-strain curve of the five materials are also As shown in Figure 4, the strain energy density-strain curve of the five materials are also exponentially distributed. In the initial stage (0–2 MPa), the strain energy density accounts for exponentially distributed. In the initial stage (0–2 MPa), the strain energy density accounts for only only to 10% to 30% ofentire the entire strain energy density. a constant strain energydensity, density,each eachmaterial’s material’s 10% 30% of the strain energy density. At At a constant strain energy 3 deformation is is different. different. For For example, example, at the Aeolian Aeolian sand sand deformation at aa strain strain energy energy density density value value of of 100 100 kJ/m kJ/m3,,the strain is 0.063, meanwhile, the loess’s strain is 0.221. Therefore, to absorb the same amount of energy, strain is 0.063, meanwhile, the loess’s strain is 0.221. Therefore, to absorb the same amount of energy, the leoss leoss deformation deformation is is larger. larger. the Solid backfilling materials are not not elastic elastic materials. materials. There There is is no no accurate accurate parameter parameter to to express express Solid backfilling materials are the elastic modulus. The modulus in the compression process is defined as the slope of the trend line the elastic modulus. The modulus in the compression process is defined as the slope of the trend line within aa range range of of the the stress-strain stress-strain curve. curve. The The expression expression is is given given as: as: within σ − σ0 EEs == σ ii − σ 0 s εε ii −−εε0 0

(2) (2)

where σ (MPa) and ε are the arbitrary point stress and strain, respectively, and σ (MPa) and ε0 are the where σii (MPa) and εii are the arbitrary point stress and strain, respectively, and0 σ0 (MPa) and ε0 are zero point stress and strain, respectively. The relationship of the five backfilling materials’ modulus the zero point stress and strain, respectively. The relationship of the five backfilling materials’ and stress is shown in Figure 5. modulus and stress is shown in Figure 5.

Energies 2019, 12, 154

7 of 14

Energies 2018, 11, x FOR PEER REVIEW

6 of 13

140

Test datas:

Aeolian Sand Fitting curves

Coal Ash

Mineral Waste Residue

Loess

Ec = 3.37σ c + 3.82

120 Eb = 7.81σ b + 4.70

R 2 = 0.99

R 2 = 0.99

100

Modulus/MPa

Gangue

Ed = 2.92σ d + 3.07

Ea = −0.66σ a 2 + 18.5σ a + 11.08

R 2 = 0.99

R 2 = 0.99

80

60

Ee = 2.48σ e + 1.32 R 2 = 0.99

40

20

0

0

2

4

6

Stress/MPa

8

10

Figure Figure 5. 5. Relationship Relationship cave cave of of five five materials materials modulus modulus and and stress. stress.

As As shown shown in in Figure Figure 5, 5, the the modulus modulus and and stress stress curve curve of of the the five five materials materials is is aa linearly linearly distributed distributed trend. However, the modulus change rate is different. Under the same stress, the Aeolian sand has has the trend. However, the modulus change rate is different. Under the same stress, the Aeolian sand highest modulus, which means this material has a high non-deformability. The fitting function can be the highest modulus, which means this material has a high non-deformability. The fitting function used to simulate the filling body deformation. can be used to simulate the filling body deformation. 5. Filling Effect of Five Solid Backfilling Materials Using Simulation Method 5.1. Engineering Background Dong Ping Coal Mine is located west of Mt. Taihang with a typical rolling terrain. The first backfill panel was 15,061, the average length along the strike direction was 286 m and the average width along the dip direction was 84 m. The panel is mined in the #15 coal seam, a lower segment of the Taiyuan group which has an average height of 6.8 m. The average depth of the coal seam is 160 m. The mining was conducted in the lower part of the coal seam with 3.0 m height and the remaining 3.8 m served as the immediate roof. The backfill material was selected as gangue. The bore hole columnar section and the rock stratum’s basic mechanic properties tested in laboratory are shown in Figure 6.

Figure 6. Simulation model.

5. Filling Effect of Five Solid Backfilling Materials Using Simulation Method 5.1. Engineering Background Dong Ping Coal Mine is located west of Mt. Taihang with a typical rolling terrain. The first backfill panel was 15,061, the average length along the strike direction was 286 m and the average width along the dip direction was 84 m. The panel is mined in the #15 coal seam, a lower segment of the Taiyuan group which has an average height of 6.8 m. The average depth of the coal seam is 160 m. The mining was conducted in the lower part of the coal seam with 3.0 m height and the remaining 3.8 m served as the immediate roof. The backfill material was selected as gangue. The bore hole

Energies 2018, 11, x FOR PEER REVIEW

7 of 13

columnar Energies 2019,section 12, 154 and the rock stratum’s basic mechanic properties tested in laboratory are shown 8 of in 14 Figure 7.

Energies 2018, 11, x FOR PEER REVIEW

140

Test datas:

6 of 13

Aeolian Sand Fitting curves

Coal Ash

Gangue

Mineral Waste Residue

Loess

Ec = 3.37σ c + 3.82

120 Eb = 7.81σ b + 4.70

R 2 = 0.99

2

R = 0.99

Modulus/MPa

100

Ed = 2.92σ d + 3.07

Ea = −0.66σ a 2 + 18.5σ a + 11.08

80

R 2 = 0.99

R 2 = 0.99

60

Ee = 2.48σ e + 1.32 R 2 = 0.99

Figure 7. 40 Figure 6. Rock Rock stratum stratum columnar columnar condition condition and and basic basic mechanical mechanical properties. properties.

Simulation Model 5.2. Simulation 20

filling effect effect of of five five different different FLAC3D numerical simulation software was used to investigate the filling backfilling materials. A numerical model was built based on the size of panel 15,061, with backfilling materials. A numerical model was built based on the size of panel 15,061, with a total 0 width in y-direction 0 of 188 m, the 2 length of one 4 zone is 2.56 m total length 8 in x-direction 10 of 380 m, the Stress/MPa length of one zone is 4 m, height in z-direction z-direction of length 70 m, the length of one zone is 4 m, and a barrier pillar width of 50 m, as shown in Figure 7. The geostress in mine site is main affected by the vertical Figure 5. Relationship6.cave of five materials modulus and stress. by monitoring monitoringthe thecrustal crustalstress. stress.Thus, Thus, uniform vertical stress, MPa, applied at stress by thethe uniform vertical stress, 2.5 2.5 MPa, waswas applied at the the of top of the model. The Mohr–Coulomb model was adopted for the rock strata and strain soften top the model. The Mohr–Coulomb model was adopted for the rock strata and strain soften model As shown in Figure 5, the modulus and stress curve of the five materials is a linearly distributed model thethe coal seam. mining coal seam, filling with thebackfilling backfilling materials which used the in coal seam. AfterAfter mining the the coal seam, filling with the materials trend.inused However, modulus change rate is different. Under the same stress, the Aeolian sand has such as dynamic deformation compaction property, modulus, operated by writing the fish language the highest modulus, which means this material has a high non-deformability. The fitting function into theused backfilling model. can be to simulate the filling body deformation. Due to the large number of irregular joints and cracks in the rock mass, a series of rock strengths and equations was created to calibration model [21–23]. The strength of the in situ rocks was determined from the laboratory scale strength using the relationship suggested by Bai [24]:

σr = (

Er q ) σl El

(3)

where σl (MPa) and Er (GPa) are the compressive strength and elastic modulus of the laboratory intact rock, respectively; σr (MPa) and El (GPa) are the compressive strength and elastic modulus of the rock mass, respectively; q is an empirical parameter taking the value of 0.7 in this paper. Similarly, the elastic modulus of the in situ rocks was determined from the laboratory scale elastic modulus using the relationship suggested by Bai [24]:

Er = 100.0186 RQD−1.19 El

(4)

where the RQD was evaluated by the integrity of rock mass by the borehole peep. The laboratory scale uniaxial compressive strength and the elastic modulus are shown in Table Figure 6.(4) Simulation model. Simulation model. 1, which were reduced by Equations Figure (3) and7. to simulate the models. The density, friction angle, cohesion, and Poisson's ratio of rock were not altered. The tensile strength was set to 10% of the large5. Filling of Fivestrength. Solid Backfilling Materials Using Simulation Method scale rock Effect compressive 5.1. Engineering Background Dong Ping Coal Mine is located west of Mt. Taihang with a typical rolling terrain. The first backfill panel was 15,061, the average length along the strike direction was 286 m and the average width along the dip direction was 84 m. The panel is mined in the #15 coal seam, a lower segment of the Taiyuan group which has an average height of 6.8 m. The average depth of the coal seam is 160 m. The mining was conducted in the lower part of the coal seam with 3.0 m height and the remaining

Energies 2019, 12, 154

9 of 14

Due to the large number of irregular joints and cracks in the rock mass, a series of rock strengths and equations was created to calibration model [21–23]. The strength of the in situ rocks was determined from the laboratory scale strength using the relationship suggested by Bai [24]: σr = (

Er q ) σl El

(3)

where σl (MPa) and Er (GPa) are the compressive strength and elastic modulus of the laboratory intact rock, respectively; σr (MPa) and El (GPa) are the compressive strength and elastic modulus of the rock mass, respectively; q is an empirical parameter taking the value of 0.7 in this paper. Similarly, the elastic modulus of the in situ rocks was determined from the laboratory scale elastic modulus using the relationship suggested by Bai [24]: Er = 100.0186RQD−1.19 El

(4)

where the RQD was evaluated by the integrity of rock mass by the borehole peep. The laboratory scale uniaxial compressive strength and the elastic modulus are shown in Table 1, which were reduced by Equations (3) and (4) to simulate the models. The density, friction angle, cohesion, and Poisson’s ratio of rock were not altered. The tensile strength was set to 10% of the large-scale rock compressive strength. Table 1. The mechanical parameters of rock and coal. Intact Rock

Lithology Packsand Sandy mudstone Gritstone Sandy mudstone Packsand Mudstone Limestone Coal Sandy mudstone Packsand

Rock Mass

RQD

EI /GPa

σ c /MPa

32.3 15.8 27.1 14.5 34.8 20.3 15.2 7.4 18.2 34.1

47.6 35.5 50.7 33.6 47.4 35.2 21.5 9.6 36.8 45.6

89 91 93 90 95 89 83 62 81 89

ER /GPa

σ cR /MPa

σ tR /MPa

18.0 9.6 17.9 8.4 25.0 11.3 6.5 1.3 7.2 19.0

31.6 25.0 37.9 23.0 37.6 23.4 11.9 2.8 19.2 30.3

3.2 2.5 3.8 2.3 3.8 2.3 1.2 0.3 1.9 3.0

In Flac3D software, the modulus of the material is defined by the bulk modulus and shear modulus. The equation of the bulk modulus and shear modulus used by the elastic modulus is as follows: (

G= K=

Er 2(1+ µ ) Er 3(1−2µ)

(5)

where the G is the shear modulus, GPa; K is the bulk modulus, GPa; µ is Poisson’s ratio. The constitutive behaviour of the five backfilling materials can use the fitting formula in the Figure 3. Similarly, the relationship of elastic modulus and the stress can be adopted from the equations in Figure 5. For example, which relationship of elastic modulus and the stress for gangue is: Er = 3.984σb + 6.142

(6)

where the σb is the stress applied to the backfilling materials, MPa. The bulk modulus and shear modulus of the gangue can be written as: (

G= K=

3.984σb +6.142r 2(1+ µ ) 3.984σb +6.142 3(1−2µ)

(7)

Energies 2019, 12, 154

10 of 14

By using the FISH language code embedded in the FLAC3D software, the stress in gangue can be monitored to recalculate the bulk modulus and shear modulus of the gangue by using the equation (7). Then reassign the bulk modulus and shear modulus in every 20 steps to simulate the gangue deformation property. Similarly, for other backfilling materials, deformation characteristics were simulated with their respective equations from Figure 5. 5.3. Stress Distribution of Five Backfilling Materials The backfilling material’s modulus and stress follow the path shown in Figure 5. The backfilling material’s vertical stresses were obtained as shown in Figure 8. Energies 2018, 11, x FOR PEER REVIEW

9 of 13

Figure8.8.Five Fivebackfilling backfillingmaterials materialscompaction compactionstress stressininSBM. SBM. Figure

The Thepressure pressurefrom fromthe theoverlying overlyingroof roofwill willbe beapplied appliedto tothe thefront frontof ofthe theworking workingface faceand andbackfill backfill materials materialsin inSBM SBM[25]. [25].As Asshown shownin inFigure Figure8,8,in inSBM, SBM,the thebackfilling backfillingmaterial’s material’sstress stressisishigher higherwhile while the theabutment abutmentpressure pressureisisreduced. reduced. This Thisindicates indicates that thatthe thebackfilling backfilling materials materials can canshare sharepart partof ofthe the pressure pressurewith withthe theboundary boundarycoal coalpillar. pillar.The Thehigher highercompaction compactionstress stressof ofthe thebackfilling backfillingmaterials materialsalso also signifies Aeolian sand and coal ash. Meanwhile, thethe abutment pressure in signifiesthe thebetter betterresistance, resistance,such suchasas Aeolian sand and coal ash. Meanwhile, abutment pressure front of working face is lower with Aeolian sand and coal ash than other backfilling materials, which in front of working face is lower with Aeolian sand and coal ash than other backfilling materials, which is important to ensure the safety of the working face, because the Aeolian sand and coal ash backfilling material bear more pressure from the overlying strata when compared to other materials.

5.4. Strain Energy Density Distribution of Five Backfilling Materials According to elastic mechanics, the strain energy density can be expressed as follows [26]:

Energies 2019, 12, 154

11 of 14

is important to ensure the safety of the working face, because the Aeolian sand and coal ash backfilling material bear more pressure from the overlying strata when compared to other materials. 5.4. Strain Energy Density Distribution of Five Backfilling Materials According to elastic mechanics, the strain energy density can be expressed as follows [26]: υε =

1 2 [σ + σ22 + σ32 − 2µ(σ1 σ2 + σ1 σ3 + σ2 σ3 )]. 2E 1

(8)

Using the FISH language and Equation (5) to configure a new variate: strain energy density, the five backfilling materials’ strain energy density are shown in Figure 9. Energies 2018, 11, x FOR PEER REVIEW

10 of 13

Figure9.9.Five Fivebackfilling backfillingmaterials materialscompaction compactionstrain strainenergy energydensity densityininSBM. SBM. Figure

As Asshown shownininFigure Figure9,9,backfilling backfillingmaterials materialscan canabsorb absorbthe theenergy energythat thatisisinduced inducedfrom fromthe the overlying overlyingroof roofdeformation deformationand andcracks. cracks.The Themore moreenergy energythe thefilling fillingmaterials materialsabsorb, absorb,the theless lessenergy energy release working face, which is a benefit to the safety a mine. example, sand backfilling releaseinin working face, which is a benefit to theofsafety of For a mine. For Aeolian example, Aeolian sand material hasmaterial a higherhas compaction resistance compared other materials, energy in backfilling a higher compaction resistanceto compared to othertherefore, materials,large therefore, large energy in stope is absorbed by the Aeolian. The Aeolian sand strain energy density is 1.67 times larger than the loess, however, the Aeolian strain energy density in front of the working face is just 32.2% of the loess. In special regions, such as shallow coal seams, ecologically susceptible areas, rock burst prone and sites with building on the surface, the backfilling materials should select the higher compaction resistance materials, such as Aeolian sand and gangue.

Energies 2019, 12, 154

12 of 14

stope is absorbed by the Aeolian. The Aeolian sand strain energy density is 1.67 times larger than the loess, however, the Aeolian strain energy density in front of the working face is just 32.2% of the loess. In special regions, such as shallow coal seams, ecologically susceptible areas, rock burst prone and sites with building on the surface, the backfilling materials should select the higher compaction resistance materials, such as Aeolian sand and gangue. 6. Engineering Application 6.1. Monitoring Scheme In order to compare the different backfilling materials stress change and filling effect, two mine sites were selected. One of the mine sites is Dongpin coal mine and the other is the Taiyuan coal mine, which is located at Yulin, Shanxi province. The Aeolian sand is used as backfilling material in Taiyuan Energies 2018, 11, x FOR PEER REVIEW 11 of 13 coal mine and gangue in the Dongpin coal mine. A series of rock stress sensors are installed in the backfilling materials to the stress in the two mine sites, the schematic of which is shown in backfilling materials to monitor monitor in Energies 2018, 11, x FOR PEER REVIEW the stress in the two mine sites, the schematic of which 11 is of shown 13 Figure Figure10. 10. backfilling materials to monitor the stress in the two mine sites, the schematic of which is shown in Figure 10.

Figure 10. Backfilling material stress monitor scheme. Figure 10. Backfilling material stress monitor scheme. Figure 10. Backfilling material stress monitor scheme.

6.2. Backfilling Material Stress 6.2.6.2. Backfilling Material Stress Backfilling Stress sand dynamic stress distribution during compaction in each mine site is The gangueMaterial and Aeolian The gangue and sanddynamic dynamic stress distribution during compaction eachsite mine site shown in Figure 11. The gangue andAeolian Aeolian sand stress distribution during compaction in eachinmine is shown ininFigure is shown Figure11. 11. 66

Initial Initialstage stage

Developing stage Developing stage

Gangue Gangue Aeolian sand

5

5

Aeolian sand

4

4

Stress/MPa Stress/MPa

StableStable stage stage

Roof begin to sink

Roof stoped to sink

Roof begin to sink

Roof stoped to sink

3

3 2

2 1

1 0

0

20

40

60

Working face advancing distance/m

80

20 40 60 80 Figure 11. Gangue Gangue dynamic stress distribution in SBM. Figure 11. dynamic stress distribution in SBM. Working face advancing distance/m

100

100

The stress in the backfill started to increase afterstress the face advancedinby 42 m (Figure 11) which Figure 11. Gangue dynamic distribution SBM. indicates the roof started to sink in the gob. The stress increase rate at its peak in the Aeolian sand is higher theingangue indicating that to theincrease roof can be stabilized at aadvanced quicker rate Aeolian sand. The than stress the backfill started after the face bywith 42 m (Figure 11) which It also shows that the peak stress where the roof is stabilized is lower for Aeolian sand (i.e., 2.7 MPa) indicates the roof started to sink in the gob. The stress increase rate at its peak in the Aeolian sand is where as for gangue the peak stress is 5 MPa. The peak stress in the Aeolian Sand is reached when

Energies 2019, 12, 154

13 of 14

The stress in the backfill started to increase after the face advanced by 42 m (Figure 11) which indicates the roof started to sink in the gob. The stress increase rate at its peak in the Aeolian sand is higher than the gangue indicating that the roof can be stabilized at a quicker rate with Aeolian sand. It also shows that the peak stress where the roof is stabilized is lower for Aeolian sand (i.e., 2.7 MPa) where as for gangue the peak stress is 5 MPa. The peak stress in the Aeolian Sand is reached when the advancing face has reached 55 m, while the peak stress for the gangue is not reached until 63 m. This indicates the Aeolian sand is a better backfill material than the gangue. 7. Conclusions 1.

2.

3.

4.

5.

The basic principles of SBM were introduced and the key factor that influences the strata control is the backfilling material. In some special areas, the requirement of strata control is higher, which does not only need considering of the cost of backfilling material but also needs investigation its compaction property. In microscopic view, five common backfilling materials meso-structure were tested, the Aeolian sand and gangue are dense, hard, and the interior structure is more stable than the other three materials, which means pressure among particles is transferred effectively. In macroscopic view, five common backfilling materials compaction laboratory test shows that in the initial stage (0–2 MPa), the strain increase is larger and the deformation rate is faster, which accounts for 70–80% of the total strain. At a constant strain energy density, each material’s deformation is different. At a strain energy density value of 100 kJ/m3 , the Aeolian sand strain is 0.063; meanwhile, the loess’s strain is 0.221. Under the same stress, the Aeolian sand has the highest deformation modulus which means this material has a high non-deformability. In macroscopic view, the simulation of SBM with five common backfilling materials were built, the result shows that higher compaction stress of the backfilling materials signify better resistance, which is found in Aeolian sand and coal ash. Meanwhile, the abutment pressure in front of working face is lower with Aeolian sand and coal ash than others backfilling materials, which is important to ensure the safety of the working face. The more energy the filling materials absorb, the less energy release onworking face, therefore increasing the safety of a mine. For example, the Aeolian sand strain energy density is 1.67 times larger than the loess, therefore, the Aeolian strain energy density in front of working face is just 32.2% of the loess, which results in lower abutment pressure. In mine site, gangue and Aeolian Sand dynamic stress distribution during compaction were monitored. In situ measurements determined that the stress increase rate is faster in Aeolian Sand, which can provide support for the roof more effectively and quickly.

Author Contributions: For this paper, P.H. completed experimental tests and wrote the paper; S.S. designed the article structure and proposed many useful recommendations; F.J. put forward study ideas; K.V.J. revised the English writing style; Z.W. collected and analyzed the data from mine site; P.N. conducted the numerical simulation and analyzed the data. Acknowledgments: This work was supported by the National Natural Science Foundation of China (51674241), and the National Basic Research Program of China (2013CB227905), and the Postgraduate Research & Practice Innovation Program of Jiangsu Province (KYCX18_1972). Conflicts of Interest: The authors declare no conflict of interest.

References 1. 2.

Zhang, J.X.; Jiang, H.Q.; Deng, X.J.; Ju, F. Prediction of the Height of the Water-Conducting Zone Above the Mined Panel in Solid Backfill Mining. Mine Water Environ. 2014, 33, 317–326. [CrossRef] Zhang, J.X.; Sun, Q.; Zhou, N.; Jiang, H.Q.; Germain, D.; Abro, S. Research and application of roadway backfill coal mining technology in western coal mining area. Arab. J. Geosci. 2016, 9, 58. [CrossRef]

Energies 2019, 12, 154

3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13. 14. 15. 16. 17. 18. 19. 20. 21. 22.

23.

24. 25. 26.

14 of 14

Miao, X.X. Progress of fully mechanized mining with solid backfilling technology. J. China Coal Soc. 2012, 8, 1247–1255. Xu, Y.N.; Li, Z.P.; Chen, S.B.; Chen, H.Q.; Yuan, H.C. Effect of coal mining collapses of the Daliuta coal mine on land desertification. Geol. China 2008, 1, 157–162. Zhang, F.; Liu, W.S. A numerical simulation on the influence of underground water flow regime caused by coal mining—A case study in daliuta, shenfu mining area. J. Saf. Environ. 2002, 4, 30–33. Wang, H.L.; Li, W.J.; Chen, Y.J. Impact on ground water by coal mining in daliuta, shenmu area. Geol. Shaanxi 2002, 20, 89–96. Lu, C.P.; Dou, L.M.; Wang, Y.F.; Du, T.F. Microseismic effect of coal materials rockburst failure induced by hard roof. Chin. J. Geophys. 2010, 53, 450–456. Zhang, J.X.; Zhou, Y.J.; Huang, Y.L. Integrated Technology of Fully Mechanized Solid Backfill Mining. J. China Coal Soc. 2012, 11, 10–13, 27. Miao, X.X.; Zhang, J.X.; Guo, G.L. Study on waste-filling method and technology in fully- mechanized coal mining. J. China Coal Soc. 2010, 35, 1–6. Miao, X.X.; JU, F.; Huang, Y.L.; Guo, G.L. New development and prospect of backfilling mining theory and technology. J. China Univ. Min. Technol. 2015, 44, 391–399, 429. Wu, X.G. Mechanical Properties of Solid Filling Materials Research and Application; China University of Mining & Technology: Xuzhou, China, 2014. Huang, P.; Li, B.Y.; Xiao, M.; Chen, Z.W.; Pei, Y.L. The Design of Critical Filling ratio in close distance coal seams by upward backfill mining technology. J. Min. Saf. Eng. 2016, 4, 597–603. Liu, Z. Compaction Properties of Gangue and Its Application in Backfilling Coal Mining; China University of Mining & Technology: Xuzhou, China, 2014. Wang, L.; Zhang, X.N.; Guo, G.L.; Zha, J.F. Quality control system framework for fully mechanized mining. J. China Coal Soc. 2013, 38, 1568–1575. Wang, L. Study on Strata Movement Mechanism and Deformation Prediction of Coal Mining with Solid Waste Compacted Filling. Ph.D. Thesis, China University of Mining & Technology, Xuzhou, China, 2014. Zhou, N.; Han, X.L.; Zhang, J.X.; Li, M. Compressive deformation and energy dissipation of crushed coal gangue. Powder Technol. 2016, 297, 220–228. [CrossRef] Deng, X.; Zhang, J.X.; Wit, B.W.; Ju, F. Pressure Propagation Characteristics of Solid Waste Backfilling Material During Compaction and Its Applications In Situ. Geotech. Geol. Eng. 2016, 34, 1631–1642. [CrossRef] Zhang, J.X.; Deng, X.J.; Zhao, X.; Ju, F.; Li, B.Y. Effective control and performance measurement of solid waste backfill in coal mining. Int. J. Min. Reclam. Environ. 2017, 31, 91–104. [CrossRef] Huang, Y.L.; Zhang, J.X.; Zhang, Q.; Nie, S.J.; An, B.F. Strata Movement Control due to Bulk Factor of Backfilling Body in Fully Mechanized Backfilling Mining Face. J. Min. Saf. Eng. 2012, 29, 162–167. Ju, F.; Huang, P.; Guo, S.; Xiao, M.; Lan, L.X. A roof model and its application in solid backfilling mining. Int. J. Min. Sci. Technol. 2017, 1, 139–143. [CrossRef] Mo, S.; Canbulat, I.; Zhang, C.; Oh, J.; Shen, B.; Hagan, P. Numerical investigation into the effect of backfilling on coal pillar strength in highwall mining. Int. J. Min. Sci. Technol. 2018, 28, 281–286. [CrossRef] Wu, S.C.; Zhang, S.; Zhang, G. Three-dimensional strength estimation of intact rocks using a modified Hoek-Brown criterion based on a new deviatoric function. Int. J. Rock Mech. Min. Sci. 2018, 107, 181–190. [CrossRef] Zhang, Y.; Ren, F.Y.; Yang, T.H.; Wang, S.Y.; Zhang, W.F.; Yu, M.X. An Improved Rock Mass Characterization Method Using a Quantified Geological Strength Index and Synthetic Rock Mass Model. Rock Mech. Rock Eng. 2018, 1–16. [CrossRef] Bai, Q.S.; Tu, S.H.; Chen, M.; Zhang, C. Numerical modeling of coal wall spall in a longwall face. Int. J. Rock Mech. Min. Sci. 2016, 88, 242–253. [CrossRef] Zhang, J.X.; Zhou, N.; Huang, Y.L.; Zhang, Q. Impact law of the bulk ratio of backfilling body to overlying strata movement in fully mechanized backfilling mining. J. Mining Sci. 2011, 47, 73–84. [CrossRef] Xu, Z.L. Mechanics of Materals; Higher Education Press: Beijing, China, 2004. © 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).