Study on the Permeability Evolution and Its Formation Mechanism of Xiaojihan Aquifer Coal Seam under Plastic Flow

Study on permeability evolution of an aquifer coal seam in Western China is of great significance for preventing water inrush disaster and realizing water-conserving coal mining. The permeability evolution of an aquifer coal seam is related to a loading path closely under plastic flow. In this work, permeability variations of the Xiaojihan water-bearing coal seam and Longde nonwater coal seam are researched using a transient method under plastic flow. The experiment results indicated the following: (1) Under the same axial strain, the permeability, relative residual strain, and confining pressure influence coefficient of Xiaojihan coal specimens all decrease in plastic flow with the increase of loading-unloading times and confining pressure, while the permeability recovery coefficient increases during this process. (2) The permeability of Xiaojihan water-bearing coal specimens decreases with the growth of axial strain in plastic flow, resulting in the increase of relative residual strain and reinforcement of plasticity. Besides, the confining pressure influence coefficient decreases and the permeability recovery coefficient decreases slightly with the axial strain. (3) Finally, the permeability of Xiaojihan coal specimens is greater than that of Longde coal specimens, while the confining pressure influence coefficient and permeability recovery coefficient of Longde coal specimens are greater than those of Xiaojihan coal specimens. The closure rate of internal cracks of the water-bearing coal specimen is lower than that of the nonwater coal specimen, which is beneficial for water storage and transport.


Introduction
It is necessary to set up safe coal pillars with different supporting and isolating functions during coal seam excavation. In the excavation process, plastic flow occurs in some areas of a coal pillar after shear yield occurring in these areas, and the deformation of coal pillars causes the change of permeability correspondingly. Multiple redistribution of the stress field also leads to the coal pillar being continuously in the abutment pressure area or pressure relief area; that is to say, the coal pillar is under a cyclic loading and unloading process. Actually, cyclic loading and unloading loads are common in the underground engineering [1][2][3][4]. For example, the front coal body is affected by the periodic breaking of the overlying basic rock stratum structure in the mining process of the coal mining face, and the moving support pressure on the work-ing face shows obvious periodicity. Coal outburst and coal caving support lifting and crushing of top coal caused by unloading confining pressure generated when Shimen is connected to coal. In fact, it is a process of repeated loading and unloading of coal, and unloading is more likely to cause instability of coal rock.
The Xiaojihan coal mine is located in the northeastern part of the nationally planned Jurassic coalfield in the northern Shanxi Province, located in Xiaojihan Township, Yuyang District, Yulin City. It is the first 10-million-ton modernized mine in the Yanheng mining area. Its main recoverable coal seam is No. 2 and No. 4-2 coal seams. With the development of the Xiaojihan coal mine, China's coal mines were first discovered, especially in the western mining areas. The mining coal seam is a special geological phenomenon of the main aquifer of the mine. According to the actual measurement, the water inflow from the minefield is mainly the aquifer water of the No. 2 coal seam, accounting for more than 80% of the total water inflow of the mine. Since the No. 2 coal seam developed is a high-pressure aquifer [5], the original geological exploration report is far from the actual situation, causing water accumulation or flood disaster in roadways easily.
The brittle failure occurs in the ordinary coal pillars due to its low water content, while plastic behavior is observed obviously in the Xiaojihan coal mine with coal seam as the main aquifer [6][7][8]. Under the repeated action of mining disturbance, the plastic zone of the coal pillar expands continuously, and plastic flow occurs at various points in the plastic zone. The deformation of the coal pillar induces the variation of porosity and permeability. It is worth noting that the permeability of the coal pillar in the plastic zone is dependent on its stress-strain state, but the single value corresponding relationship between permeability and strain is no longer existing. Meanwhile, the position of phase points in the plastic zone is determined by the flow path. Thus, the permeability of the coal pillar after failure depends not only on the current stress-strain state but also on the plastic flow path of points in the coal pillar plastic zone. It is necessary to study the permeability of coal by understanding the evolution of permeability under various hysteresis loops.
The coal measured rock is under different stress states and levels due to mining excavation. In order to study the evolution of permeability of coal measured rock under different stress states and levels, lots of experiments are conducted to study the mechanism of permeability, especially focused on the relation between permeability and stress-strain curve, as well as the relation between permeability and confining pressure. Permeability tests of coal measured rock during the whole stress-strain process were conducted by scholars [9][10][11][12][13][14][15] indicating that permeability decreased slightly in the initial stage with the increase of axial strain, which was caused by the closure of primary cracks. Then, permeability increased slowly with the increase of axial strain due to initiation of new cracks and expansion of original cracks. And then, the permeability increased in magnitude since the main cracks were formed in the elastoplastic stage. Finally, the permeability decreased to be sϒ table due to the compaction of induced fractures. Some experts [16][17][18] studied the permeability of coal measured rock under cyclic loading and unloading, concluding that the unloading curve of rock did not coincide with the loading curve, forming a closed plastic hysteresis loop, and also obtained the evolution of residual strain with the cyclic times and the damage evolution of coal measured rock. Besides, the relation between confining pressure and permeability was researched [19][20][21][22][23][24][25][26]. The experiments and numerical results indicated that the relationship of confining pressure and coal permeability could be described by an exponential function. In all, most of the studies of coal permeability were conducted on an intact coal specimen since the flow characteristics are more regular and the experimental methods are relatively simple.
Based on the above, the permeability of coal is determined by its deformation and stress state. However, the following aspects are not considered in the previous studies.
(1) At present, most of the testing methods are to increase the axial strain or stress gradually under constant confining pressure, which is a simple loading mode since the influence of different stress states on coal permeability and the irrecoverability of deformation is not considered. (2) It is assumed that a specific strain/stress state corresponds to a specific permeability; that is, the permeability of a rock specimen has a fixed expression, without considering the effect of plastic flow. According to Reference [27], elastoplastic behavior has a significant influence on the permeability of coal, especially on the coal with massive water content and fractures. Furthermore, the research of the plastic flow of coal in the main aquifer of a coal seam is seldom reported. (3) For the analysis of the coal seam as the main aquifer, scholars [28,29] mainly focus on the geological environment of coal seam and the comparison of its strength and permeability with those of upper and lower strata. However, the comparison of mechanical properties and permeability between regular coal seam and main aquifer coal seam is not well understood under the same geological environment.
In this work, the permeability tests are conducted on the Xiaojihan water-bearing coal seam and Longde nonwater coal seam with different cyclic loading and unloading under plastic flow. Then, the reason of water storage in an aquifer coal seam is analyzed through the comparison of permeability between the aquifer coal seam and the nonaquifer coal seam. The study of permeability evolution of the aquifer coal seam under plastic flow is of great significance to the study of the mining stress field, mining fracture field, and water seepage field, as well as the study of the seepage instability behavior of a waterbearing coal seam and prevention of water inrush disaster.

Materials and Testing Method
2.1. Area of Research. Yuheng Mining Area is located in the border area between Maowusu Desert and Loess Plateau, which is a typical arid and water-deficient area. The water content of coal seams in mining areas of Hongshixia, Xihongdun, Boluo, Yuandatan, Kekegai, Dahaize, Shiliutai, Wusuhaize, Hongshiqiao, and Longde is extremely low (Figure 1). In the early 21st century, the special geological structure of a coal seam as the main aquifer was first found in the Xiaojihan coal mine of the Yuheng Mining Area, which  Figure 1: Minefield distribution map of the Yulin mining area. 2 Geofluids is rare even in south China. Therefore, "the study of the cause of coal seam becoming the main aquifer" has become a hot research topic in the field of coal science and technology. The coal seam of the Xiaojihan coal mine in Yanheng Mining Area is an aquifer. The water seepage source during the tunneling process of the coal seam roadway is mainly the No. 2 coal seam fissure water and the coal roof sandstone fissure water. Field observations and studies have shown that the No. 2 coal seam is stable in layer structure and simple in structure, but the coal seam fissures are developed, and the direct bottom plate has a certain thickness of the aquifuge, which has a good water storage space [5]. The Xiaojihan No. 2 coal seam is the main aquifer due to its roof rock of water-rich characteristics. The Quaternary unconsolidated groundwater is mainly supplied by precipitation and lateral recharge, partly by desert condensate recharge, irrigation return water, and canal water leakage recharge. Fractured diving in the bedrock weathering zone is recharged directly by atmospheric precipitation infiltration in the exposed area, while the rest is recharged by infiltration of overlying unconsolidated layer diving. The recharge of deep confined water mainly comes from the infiltration of the overlying aquifer, and the overflow channel is mainly the tortuous path of fracture in the aquifer. Thus, adequate source and recharge of fissure water provide favorable conditions for the formation of the main aquifer in a coal seam. The coal seam groundwater runoff roadmap is shown in Figure 2. from the Xiaojihan water-bearing coal seam and adjacent mining area Longde no-water coal seam with a size of 50 mm × 100 mm. By measuring the wave velocity of coal samples, a more uniform coal sample is screened out. In order to study the basic properties of coal related to its permeability, such as initial fracture characteristics and initial porosity, scanning electron microscopy (SEM) tests and mercury intrusion tests are conducted on coal specimens.
A scanning electron microscope Quanta 250 is used in the SEM tests. The scanning image of the coal specimen is obtained with a magnification of 400x. Then, the optimal threshold of the image is determined by selecting the threshold algorithm using a self-developed program. Finally, the fractures in coal specimens are obtained through binary processing of the SEM image as shown in Figure 3.
It can be seen from Figure 3 that a large number of fissures are observed in the Xiaojihan coal sample, and the fissures are curved and have no obvious trend. Moreover, there are many cracks, and the coal structure is loose. Meanwhile, less cracks are found in the Longde coal mine, and the cracks are rarely staggered, and the opening of the fissures is basically unchanged.
An Autopore 9500 mercury porosimeter is used in mercury intrusion tests, and the experimental results of Xiaojihan coal specimens and Longde coal specimens are shown in Table 1.
In order to compare the initial pore structure of these two different coal specimens directly, Figure 4 is drawn correspondingly.
As shown in Table 1 and Figure 4, the initial porosity and permeability of the water-bearing coal seam are 9.89% and  3 Geofluids 8:51 × 10 −15 m 2 , respectively. The initial porosity, total mercury amount, initial permeability, total pore area, and average pore size of Xiaojihan coal specimens are all greater than those of Longde coal specimens, indicating that the water storage capacity of the Xiaojihan coal seam is stronger than that of the Longde coal seam.

Testing
Method. The methods of conventional triaxial tests and permeation tests under plastic flow are introduced based on the test principle and equipment.

Testing
System. The schematic diagram of the penetration test system is shown in Figure 5. The penetration test system is composed of an axial loading system (MTS816), confining pressure system (oil pump), permeation loop (hydraulic pump), and permeation meter (self-made cylinder). Its main functions included the following: (a) applying axial load to the rock sample through the MTS816 rock mechanic control system; (b) applying confining pressure to the rock sample through the oil pump; besides, A is the inlet end of the oil; and (c) applying pore pressure (osmotic pressure) to the upper and lower end faces of the rock sample. B and C are the inlet ends of the bottom water and top water, respectively. The permeation circuit consists of a per-meameter, a manual pump, a pressure sensor, a paperless recorder, a pressure gauge, and a storage tank.
A standard cylindrical coal sample is placed in the cylinder. The porous disk is arranged above and below the coal sample as shown in Figure 5. The role of the porous disk is to distribute fluid pressure and equalize the percolation velocity. The upper steel plate is designed as a spherical surface to match the convex upper surface of the contact. These two parts can be rotated relative to each other to adjust the parallelism between the coal sample, the porous disk, the steel plate, and the piston. The design above avoids eccentric loads.
In order to study the mechanism of coal permeability evolution under plastic flow, it is necessary to measure the radial residual strain. And a bow sensor is introduced, which is connected to a static strain guage at the central hole of the pin shaft as shown in Figure 6. The data acquisition system consists of a static resistance strain indicator, a pressure sensor, a paperless recorder, and a notebook computer. The system can collect radial strain and pore pressure values in real time.

Experiment Scheme.
In order to simulate the evolution of coal seam deformation and water seepage caused by unloading of underground roadway excavation and the evolution of permeability versus confining pressure of coal specimens from the main aquifer coal seam under plastic flow, the conventional triaxial tests are conducted to make the coal specimens reach the state of plastic zone; then, the confining pressure is loaded and unloaded to measure the permeability during the process. The specific procedures are as follows:   Geofluids (a) The confining pressure of 2 MPa is applied to coal specimens at a loading rate of 0.5 MPa/s; then, the specimen is loaded in the axial direction to postpeak stress state A * with a loading rate of 0.003 mm/s (σ 1 * = 90%σ max ). And the axial strain is ε 1 * at this state and remains constant (Figure 7). If the pressure of the hysteresis' starting point, which means σ 1 * , closes to the peak, the pressure's hysteresis would be more plump and the change of permeability of rock samples under plastic flow would be much more comprehensive. The change of permeability under plastic flow can be reflected much more comprehensively. TestStar controls axial pressure of state A * under 90 percent of peak pressure, which is MTS816 experiment systems' software. However, if exceeding 90%, it may lead to out of control, and the coal sample would be crushed. Because of the large discreteness of coal, it is not easy to determine the state of each coal sample reaching the plastic zone. Therefore, 90% of the peak stress of each coal sample in the postpeak stage is taken as the initial value of the plastic zone (b) Xiaojihan coal seam buried depth is 350 m [30], its covered rock formation's average density is 2600 kg/m 3 , and vertical ground stress is 8.92 MPa.
According to the definition of the pressure coefficient [31], lateral pressure vector is 1.11, so we obtain horizontal ground stress 9.81 MPa and the pore pressure difference 1. Penetration tests were carried out separately at the above 21 stress state points ( Figure 8). Axial strain ε * 1 is divided into five grades, and plastic zone state ε * 1 is determined according to the actual situation. The initial confining pressure of each water-bearing specimen is 2 MPa. In order to compare the permeability of two kinds of coal specimens under the same stress state, if the strain ε * 1 of the nonaquifer coal sample is different from that of the aquifer coal sample, the axial strain is applied to make it reach the same level (ε * 1 ). Pore pressure difference at upper and lower ends of the permeability test under each confining pressure is 1.5 MPa (ΔP = 1:5 MPa) as shown in Table 2.

Experiment
Principle. Due to the low permeability of coal [32], 90% of peak stress of the coal specimen in the   5 Geofluids postpeak stage is taken as the initial value of the plastic zone. Although macrocracks could be observed in the coal specimen, the seepage velocity is still at a low level, and water could not be collected at the outlet. Therefore, the transient method is applied to determine the permeability of the coal specimen, and the principle of this method is shown in Figure 9.
The height and cross-sectional area of the coal specimens are H and A, respectively. At the beginning of the experiment, the initial pressures at the two ends of the coal specimen are P 10 and P 20 ðP 10 > P 20 Þ. The axial pressure gradient in the coal specimen is ðP 20 − P 10 Þ/H. During the process of infiltration, the liquid in tank 1 moves to tank 2 by passing through the coal specimen. As a result, the pressure of tank 1 continuously decreases, while that of tank 2 continuously increases, and the pressure gradient gradually decreases until the pressures of two tanks reach an equilibrium state (see Figure 10). We assume that the mass flow rate of the liquid entering the coal specimen from tank 1 is q. If the coal specimen is water-saturated, the mass flow rate of the liquid entering water tank 2 from the coal specimen is q . The seepage velocity in the rock sample is V = q/ρA. Water is considered as compressible, and its compression factor is defined as follows: where c f is the compressibility of water and ρ is the density of water. Equations (2) and (3) can be derived from dρ = −qdt/B and q = ρAV: Equation (4) can be induced from equations (2) and (3): According to Darcy's law, the seepage velocity is defined as Substituting equation (5) into (4), the following equation can be obtained: During the test, pressure data are collected at equal intervals τ, with a total of n acquisitions. At the end of data collection, the pressures in the water tanks at t f = n π are P 1f and P 2f , respectively. The integral of equation (6) is obtained: A typical confining pressure-permeability hysteresis curve is shown in Figure 11. In order to describe the influence of axial strain on the confining pressure-permeability hysteresis curve, two parameters named confining pressure influence coefficient and permeability recovery coefficient are introduced here. And the confining pressure influence coefficient is defined in  Figure 9: Principle of the transient method. The permeability recovery coefficient is defined as where k ðiÞ e is the permeability at the end of cycle times i, i = 1, 2, 3, 4.

Results and Discussion
In this section, the coal sample is first brought to the state of the plastic zone by a conventional triaxial compression test. Secondly, the hysteresis curve of permeability is obtained by the permeability test under plastic flow, and the variation law of permeability is analyzed.

Results of Conventional Triaxial Compression Tests.
A conventional triaxial compression test was carried out on the coal sample to make it reach the state of the plastic zone. The confining pressure is kept as 2 MPa. The stress-strain curves of two kinds of coal specimens are shown in Figure 12, and the damage load, peak stress, and peak strain are shown in Table 3.
In order to obtain the same ε 1 * with the water-bearing coal specimen, the non-water-bearing coal specimen is in the plastic zone state of 70%-90% peak stress by controlling the MTS816 program. If ε 1 * of the non-water-bearing coal specimen is in a large difference with that of the waterbearing coal specimen, the non-water-bearing coal specimen is regarded as unqualified. The results are shown in Table 4.

Variation of Relative Residual Strain under Plastic Flow.
The evolution of permeability is closely related to the deformation of coal specimens. In order to study the mechanism of permeability evolution, plastic flow tests are conducted on water-bearing coal specimens under constant strains of 1.61%, 1.72%, 1.79%, 1.81%, and 2.01%, respectively. The hysteretic curves of circumferential strain-circumferential stress (confining pressure) are shown in Figure 13.
It can be seen from Figure 13 that the hysteresis curve of radial strain-radial stress of the coal sample under constant axial strain has the following characteristics. When confining pressure is unloaded to the initial value at each stage of axial strain, the circumferential strain lags behind the circumferential stress, and the circumferential strain could not return to the initial value, resulting in irreversible deformation. Therefore, the loading-unloading curve of circumferential strain-circumferential stress does not form a closed hysteresis loop.
Relative residual strain is introduced to describe the plastic deformation [33]. The transverse strain corresponding to the lower limit circumferential stress point of each cycle in the circumferential strain-circumferential stress hysteresis curve is extracted, and the relative residual strain of the cycle times i is calculated according to where Δε 3i is the relative residual transverse strain at cycle times i and ε 3i is the residual transverse strain at cycle times i. The relative residual transverse strain of each cycle under each stage of axial strain is shown in Table 5.
The increase of the cyclic coefficient corresponds to the increase of confining pressure. In order to observe the influence of axial strain and confining pressure on relative residual strain, the relative residual strain versus confining pressure under different strain states is drawn according to Table 5 (see Figure 14).
The following is concluded from Table 5 and Figure 14. (a) The relative residual strain is relatively large at the initial increase of confining pressure under the same axial strain. The relative residual strain decreases slowly with the increase of confining pressure, which indicates that plastic deformation is produced mainly in the first cyclic loading and unloading with the increase of confining pressure, resulting in a large irrecoverable deformation. And the coal specimens are mainly under a neutral load during the last three loading and unloading cycles; that is to say, the elastic deformation of the coal specimens is produced after plastic deformation. (b) The relative residual strain increases with the increase of axial strain under the same confining pressure, which shows the irrecoverable deformation and plasticity increase with the increase of axial strain.

The Permeability versus Confining Pressure and Axial
Strain under Plastic Flow. Plastic flow permeation tests of two kinds of coal specimens under different constant strains are conducted, and the number of loading and unloading cycles was four. The experiment results are shown in Figures 15-19, and permeability of water-bearing and non-water-bearing coal specimens from the beginning of each cycle to the reverse point of each cycle is shown in Tables 6-10.

Geofluids
In order to study the influence of axial strain and confining pressure on permeability under plastic flow, the confining pressure influence coefficient and permeability recovery coefficient of water-bearing and non-water-bearing coal specimens under different axial strains are obtained as shown in Tables 11-14. 3.3.1. The Evolution of Permeability versus Confining Pressure. As shown in Figures 15-19, the permeability of the coal specimen decreases exponentially with the increase of confining pressure under the same axial strain, which is in good agreement with previous studies [18][19][20]. At the initial stage of confining pressure loading, the permeability of water-bearing coal decreases rapidly. Then, the decrease of permeability becomes slower with the increase of confining pressure. According to the pore fracture theory of coal rock mass [34], water flow passage in coal measured media contains pores and cracks. Since the stress surface of the pore is arch and the stress is compressive, the deformation of pore volume is relatively small. While coal crack lacks the support of a skeleton matrix, its deformation is relatively large under compression. Therefore, the main sequence of the deformation of the fluid passage caused by the increase of the confining pressure is the following: the deformation of the fracture under pressure and the closure ⟶ the deformation of the partial pores. During the whole process, the compression deformation and closure of coal cracks and pores become more and more difficult, resulting in the stability of permeability with the increase of confining pressure.
As shown in Tables 11 and 12, permeability of waterbearing coal specimens is particularly large at the low confining pressure stage. Since most of the cracks and microchannels   8 Geofluids are closed after the plastic deformation of the coal specimen due to confining pressure loading, a few of microcracks that are compressible remained. Only large pore volume could be compressed, and compressive deformation is formed, resulting in the change of permeability. Therefore, the influence of confining pressure on permeability of the second The results indicate that the permeability reduction rate decreases with the increase of confining pressure.
The relationship between permeability and confining pressure is several-for-one under the same axial strain as shown in Figures 15-19. The permeability of the waterbearing coal specimen in the unloading stage is less than that in the loading stage, which indicates that the irreversible deformation of the coal specimen occurs during the confining pressure loading stage. And there is an obvious stress lag effect in the recovery of permeability of water-bearing coal specimens, which corresponds to the stress-strain lag in the section.
The permeability recovery coefficient grows with the increase of confining pressure as shown in Tables 13 and  14. During four unloading processes, the permeability increases to a certain extent, but it could not recover to the initial value. And the recovery of permeability gets slower with the increase of confining pressure. The permeability loss rates of coal specimens in the last three loading and unloading processes are much smaller than that in the first loading and unloading process. The results show that the waterbearing coal specimens are mainly subjected to neutral loading in the last three loading and unloading processes; that is to say, the elastic deformation after plastic deformation can be restored to a greater extent after unloading. However, the plastic deformation of coal specimens is dominant due to the existence of opening cracks and throats during the first loading and unloading process. For coal specimen X-2.01%, the permeability of the coal specimen decreases from 2:48 × 10 −14 m 2 to 1:98 × 10 −14 m 2 in the first unloading process, and the permeability recovery coefficient is 0.80. The permeability recovery coefficients are 0.80, 0.84, 0.85, and 0.86 during the four unloading processes, respectively.

The Evolution of Permeability versus Axial Strain.
The confining pressure influence coefficient of the waterbearing coal specimen decreases with the increase of axial strain under the same confining pressure shown in Tables 11 and 12. According to Table 5, the plastic deformation becomes greater with the increase of axial strain due to the closure of lots of cracks and pores during this process. Thus, the rest of microcracks and pores not compressed becomes less, decreasing the influence of confining pressure on permeability.
The permeability recovery coefficient of the waterbearing coal specimen decreases with the increase of axial strain under the same confining pressure shown in Tables 13 and 14, indicating that the plastic deformation increases with the growth of axial strain. According to Table 5 and Figure 13, when the coal specimen is compressed more compactly, the irrecoverable deformation of the waterbearing coal specimen becomes greater, resulting in the decrease of the permeability recovery coefficient.

Brief Analysis on the Cause of Water Containing in
Xiaojihan Coal Seam. In order to compare the confining pressure influence coefficient and permeability recovery coefficient of two kinds of coal specimens under the same stress level, Figure 20 is drawn according to  The confining pressure influence coefficient and permeability recovery coefficient of nonwater coal specimens are greater than those of water-bearing coal specimens under the same stress level, which shows that the permeability of the water-bearing coal specimen is more insensitive to the change of stress. In the process of periodic displacement of the coal seam, the closure rate of internal cracks of the water-bearing coal specimen is lower than that of the nonwater coal specimen, which is beneficial for water storage and transport.
The order of magnitude of permeability of the waterbearing coal seam is greater than that of the nonwater coal seam under the same stress level. Water from roof strata can easily penetrate into the aquifer coal seam, which is the main reason for the coal seam in the Xiaojihan coal mine as the main aquifer.
There are sufficient fissure water source and recharge in the Xiaojihan No. 2 coal seam as stated in Section 2.1. From      Figure 19: The permeability evolution during the loading-unloading process with the axial strain 2.01%. Table 6: The permeability evolution during the loading-unloading process with the axial strain 1.61%.  Table 7: The permeability evolution during the loading-unloading process with the axial strain 1.72%.   Table 9: The permeability evolution during the loading-unloading process with the axial strain 1.81%.

Conclusions
In view of the special geological structure of the "coal seam in the western mining area as the main aquifer," the permeability tests of the Xiaojihan coal specimens and Longde coal specimens under plastic flow are conducted. The main conclusions are obtained as follows.
(1) The relation between confining pressure and coal permeability is negative exponential. The permeability, relative residual strain, and confining pressure influence coefficient of Xiaojihan coal specimens all decrease with the increase of loading-unloading times and confining pressure. However, the influence of confining pressure on coal permeability during the second loading process is much weaker than that during the first loading process. The plastic deformation is mainly produced during the first loading and unloading process, while the water-bearing coal specimens are mainly subjected to neutral loading in the last three loading and unloading processes, which is the elastic deformation after plastic deformation. The plasticity produced with an increase in the number of times of confining pressure loading and unloading is reduced, and the influence on the permeability is reduced (2) The permeability, confining pressure influence coefficient, and permeability recovery coefficient of the water-bearing coal specimen decrease with the axial strain, but the relative residual strain increases with the axial strain. It is indicated that the plastic deformation of the coal sample increases, and the influence of the confining pressure on the permeability decreases with the axial strain (3) The reasons for the main aquifer of the Xiaojihan coal seam are mainly explained in the following aspects. First, the Xiaojihan coal seam contains plenty of primary cracks with wider width and strong capacity of water storage. Second, the permeability of the Xiaojihan coal specimen is greater than that of the Longde coal specimen; thus, water from roof strata is more likely to penetrate into Xiaojihan coal seams. Besides, the confining pressure influence coefficient and permeability recovery coefficient of Xiaojihan coal specimens are smaller than those of Longde specimens, which indicates that the permeability of the Xiaojihan coal specimen is more insensitive to the change of stress. In the process of periodic displacement of the coal seam, the closure rate of internal cracks of the Xiaojihan coal specimen is lower than that of the Longde coal specimen, which is beneficial for water storage and transport. Those are the main reasons of Xiaojihan coal seams as the main aquifer The plastic behavior of coal specimens from the coal seam as the main aquifer is extremely obvious. This paper mainly studies the variation of permeability of the aquifer with confining pressure and axial strain under plastic flow. However, the theoretical relationship between stress (strain) and permeability under different paths in plastic flow is subject to study in the future.

Data Availability
The data used in the manuscript is from experiment. you can contact Jingna Guo (tb18030001b4@cumt.edu.cn) if you want learn more about our work.