Prediction of gas pressure in thin coal seams in the Qinglong Coal Mine in Guizhou Province, China

Thin coal seams in mines usually lack gas data. Thus, preventing and controlling gas outbursts of thin coal seams are difficult. In this study, a coal structure index, which is used to express the damage degree of coal, was estimated by logging curve. In accordance with the contour line of the floor of the coal seam, structural curvature was calculated to express the complexity of the coal seam structure quantitatively. Subsequently, relationships among the burial depth, thickness, coal structure index, structural curvature were analyzed on the basis of the gas pressure of coal seam. The gas pressure values of the coal seams of Nos. 22, 24, and 27 in the study area were predicted by multiple linear regression (MLR) and were then verified and analyzed. The deviation rate of the MLR method was 6.5%–19.7%, with an average of 13.0%. The average deviation rate between the predicted value and the measured value was 11.6%, except for the measuring point of No. 2, which had a large deviation. Results show that the prediction accuracy of the aforementioned method is acceptable and has practical value in the prediction of gas pressure in thin coal seams without measured data. The results in the gas pressure prediction provide a basis for evaluating the risk of gas outbursts in thin coal seams.


Understanding of coal seam gas occurrence characteristics
In mining areas, the occurrence of gas is generally affected by the coalification degree of coal, the permeability of surrounding rocks, the geological structure, and the burial depth and thickness of coal seams (Lu et al. 2019;Li et al. 2020a).
The spatial distribution of gas content is firstly controlled by geological structure. Regional and mine structures form different sealing conditions, which have a considerable control effect on the distribution of gas content in coal seams (Moore 2012;Gao et al. 2020;Quan et al. 2020). Given a large compressive stress in compressional fault, fault gouge and mylonite with compacted structure are relatively developed, and the permeability of gas is poor (Bustin and Clarkson 1998;Zhang et al. 2016;Tong et al. 2019). Gas migration along the fault or vertical fault plane is relatively difficult, thereby becoming conducive to gas accumulation; meanwhile, the opposite holds for extensional fault. The sealing property is better when the argillaceous content of the roof rock of the coal seam is higher, which is beneficial to gas preservation (Jiang et al. 2010;Zhao et al. 2014b;. The greater the burial depth of the coal seam, the larger the in-situ stress and the greater the distance of gas escape, which is also conducive to gas sealing. The greater the thickness of the coal seam, the more conducive it is to gas preservation (Lu et al. 2018b;Kang et al. 2020;Li et al. 2020b). However, some studies have shown that the depth effect is no longer obvious when the depth of coal seam exceeds the critical point of 800-900 m (Li et al. 2020a).
The geological structure conditions of mining areas considerably influence the present in-situ stress, and the depressions formed by fault structures are high-lying stress areas, that are conducive to gas preservation Li et al. 2020c;Lu et al. 2021). Given the differences in the associated structures and stress fields above and below the neutral layer, the large fold structures have varying control effects on gas occurrence (Li et al. 2017;Lu et al. 2018a;Meng et al. 2019;Zhang et al. 2020a).
Coal structure is an important factor that affects gas occurrence. Various types of tectonic deformed coal have different physical properties, such as porosity and permeability, which affect gas occurrence (Guo et al. 2014;Song et al. 2017;Liu et al. 2020).
In practice, the gas occurrence law in mining areas must be predicted. Through the pre-analysis of the main influencing factors of gas content or pressure, the analytic algorithm, multiple linear regression (MLR), artificial neural networks, and other methods can achieve satisfactory results (Wei et al. 2009;Dai 2016;Li et al. 2019;Wang et al. 2019). Some scholars use geostatistics method to estimate the gas content of coal seam under GIS environment, which also has a certain effect (Vaziri et al. 2015).
The factors that influence gas occurrence have been studied deeply by previous researchers, and various methods have been used to predict the occurrence law of gas in mining areas. However, for the realization of these methods and the discussion of their reliability, further studies remain necessary to improve the accuracy of gas prediction continuously.

Current development status of thin coal seams
In China, the thin coal seams were mined early to meet the energy demand of economic construction. The thickness of thin and extremely thin coal seam in China is specified as 0.8-1.3 m and less than 0.8 m, respectively. Other countries generally set the thickness of thin coal seam as 0.6-0.8 m and 1.5-1.6 m at the lower and upper limits, respectively.
In Southern China, including Guizhou Province, the coal bearing strata are mainly in Triassic. The greatest feature of coal seam development is the existence of many coal seams and thin coal seams. At present, literature reports on thin coal seam are scarce, among which Chinese scholars mostly focus on mining simulation experiment (Ma et al. 2008) and fully mechanized mining technology in thin coal seam (Wang et al. 2011;Wang et al. 2012;Zhao et al. 2014a;Liu et al. 2017) and geological or engineering conditions related to thin coal seam mining. The literature related to gas includes the pressure relief mining of protective layer to prevent gas outburst (Li 2014;Cao et al. 2018) and gas drainage technology in thin coal seam (Li 2011;Wang et al. 2017). Other countries in the world that have carried out research on thin coal seams include Poland and Ukraine, mainly Poland. Relevant reports include geological resources and distribution characteristics of thin coal seams (Krowiak 2011), the analysis of reserve proportion of thin coal seams (Dyczko 2007), economic rationality evaluation of thin coal seam mining (Piwniak et al. 2007;Saluga 2008), and discussion on unmanned mining of thin coal seams with advanced consciousness (Litvinsky 2007).
The Guizhou Province is located in the copulae of the circum-Pacific tectonic tract and the Western Tethyan tectonic domain. Its tectonic evolution controls coal-forming environments and gas occurrence. During the Yanshan movement, a series of folds, overthrust fault, and nappe tectonics were formed and resulted the development of tectonic deformed coal. These factors led to serious coal and gas outburst in Guizhou Province (Lu et al. 2020;Zhang et al. 2020b). The Qinglong Coal Mine in the Western Guizhou Province has developed many thin coal seams. The Chinese law stipulates that gas outburst risk assessment must be conducted for thin coal seams with a thickness of more than 0.3 m (State Administration of Coal Mine Safety 2019). However, obtaining the measured gas data of thin coal seam is difficult in general. So is Qinglong Coal Mine. Under this background, the gas pressure prediction of thin coal seams was conducted by using relevant methods in this study, thereby providing a basis for evaluating the risk of gas outburst in thin coal seams.
Perhaps other countries outside China do not have the urgency of mining thin coal seams at present, but coal resources are non-renewable, and the possibility to continue mining thin coal seams in future is great because of energy demand. Therefore, the research experience and results in the gas prediction of "thin coal seams" are of great significance to the development and utilization of thin coal seams in future.

Geological conditions of coal seam gas
The Qinglong Coal Mine is approximately located 110 km Northwest of Guiyang City, Guizhou Province. The length of the mine field from NE to SW is approximately 7.89-9.93 km, the width from NW to SE is approximately 2.31-3.20 km, the area is 20.65 km 2 , and the mining elevation is between + 1300 and + 700 m (Fig. 1). The strata exposed in the study area from old to new are the Maokou Formation (P 2 m) of the Middle Permian, the Emeishan basalt formation (P 3 β), the Longtan (P 3 l) and Changxing Formations (P 3 c) of the Upper Permian, and the Yelang (T 1 y) and Maocaopu Formations (T 1 m) of the Lower Triassic, and Quaternary (Q). The coal bearing strata in the study area include the Longtan Formation (P 3 l) of the Upper Permian, which can be divided into P 3 l 1 and P 3 l 2 in accordance with lithologic association. The thickness of P 3 l is 158.50-188.30 m, and the average is 172.05 m.
The structural features of the study area mainly include a series of wide and gentle folds with NE or NNE trends. The study area is located in the Northwest wing of the Gelaozhai anticline, and its structural characteristics are consistent with the regional tectonic framework. The mine field is dominated by NE trending structures with normal and reverse faults, large-scale and long extensions, and large dip angles. The secondary faults are mainly reverse faults with smallscale and short extensions (Fig. 1). Mining and excavation data show that small faults and folds are well developed in the mine field. Most of them are faults with a throw of 0.5-10 m, and the dip angle of the small fold strata varies from 3 to 90°.
The entire coal bearing strata has 15-26 coal seams, including 9 minable and partially minable coal seams (Nos. 2,3,16,17,18,22,24,27,and 30). The coal seams of Nos. 16 and 18 are the main minable coal seams. The average thickness of all coal seams is 16.26 m, and the coal bearing coefficient is 9.5%. Although the distribution of these thin coal seams is unstable, partially minable or not, the gas content is high and has the risk of gas outburst. For example, the average thickness of the coal seam of No. 22 is 1.23 m, and gas outburst occurs many times when the coal seam has been exposed.

Main influencing factors and indexes of gas pressure
Gas pressure is an important index for the risk assessment of gas outbursts in coal seam. However, measured data of gas pressure in thin coal seam at the initial stage of mining are lacking due to the unstable distribution and partial minability. Therefore, the gas pressure of thin coal seams must be predicted to ensure production safety.
The coal seam of Nos. 16 and 18 in the study area have substantial gas data. Data on burial depth, coal thickness, coal structure index, and structural curvature of coal seam were analyzed by single-factor analysis, and gas pressure have a certain correlation with the aforementioned factors (Fig. 2). Among these four indexes, data on burial depth and thickness of coal seam could be directly obtained, whereas those on coal structure index and tectonic curvature must be indirectly obtained through calculation.

Calculation of coal structure index
Among the commonly used logging curves in coal field geological exploration, the apparent resistivity curve (DLW) and artificial gamma curve (HGG) are the obvious indexes that reflect the difference in the degree of coal damage. Given that pore and fracture are well developed in damaged coal, the water content of coal increases under the immersion of drilling fluid, which leads to the decrease in the amplitude of DLW and the increase in the amplitude of HGG in coal seams with a high degree of damage. The changes in these two curves are the main bases for judging the degree of damage in coal (Fu et al. 2003). The variation of well diameter can also be used as a reference. The coal with a higher damage degree usually causes the enlargement of the well diameter ( Fig. 3). At this point, the degree of coal damage can be categorized into three levels: the value of 0-1 represent no to small damage, the value of 1-2 represent moderate damage, and the value of 2-≥ 3 represent relatively large damage to large damage.
In accordance with the above principles, the values of DLW, HGG, and borehole diameter change were calculated to determine the degree of coal failure comprehensively.
Apparent resistivity (Unit: Ω·M): Artificial gamma (Unit: CPS) and change in borehole diameter (Unit: mm): where X i is the reading value of the logging curve of the borehole, X min is the minimum reading value of the logging curve, and X max is the maximum reading value of the logging curve reading.

Calculation of construction curvature
Structural curvature is a dimensionless value calculated in accordance with the shape, density, and bending degree of the coal seam floor's contour line, which can reflect the deformation degree of the coal seam floor well. The curvature of the contour line of the coal seam floor can reflect the local change in the morphology of the coal seam floor; and this curvature is equivalent to the reflection of the structural form of the coal seam. The calculation formula is expressed as follows (Shen et al. 2010): where K is the curvature value of a point in the coal seam floor; z is the coal seam floor elevation, z = f (x, y) , z � = f x , and z �� = f � x . The contour lines of the coal seam floor were divided into many square cells with 500 m × 500 m. The curvature of the floor contour lines in the North, East, South, and West directions were calculated at the center point of the cells, and the maximum value was taken as the structural curvature of the calculation point. Through this method, the structural curvature of the floor of the coal seam with No. 22 in the study area was calculated (Fig. 4).

Prediction by MLR
With Y taken as variable, and burial depth, coal thickness, coal structure index, and structural curvature taken as explanatory variables X, X = { X 1 , X 2 , X 3 , X 4 . X i ∈ (0, + ∞)}, on the basis of 62 groups of data (hereinafter referred to as "early measured data") measured in the early stage of the coal seams of Nos. 16 and 18, MLR analysis was conducted in accordance with the following formula: where Y is the gas pressure, MPa; X is the explanatory variable; X 1 is the burial depth, m; X 2 is the coal thickness, m; X 3 is the coal structure index (dimensionless); X 4 is the structural curvature (dimensionless); b is the regression coefficient, b = (b 0 , b 1 , b 2 , b 3 , b 4 ); and n is the sample number of MLR, n = 62.
where Y is the gas pressure, MPa.
The predicted gas pressure values of the coal seams of Nos. 22, 24, and 27 were accordingly plotted. The plots show that the gas pressure of each coal seam from shallow to deep is approximately 1.2-3.0 MPa, thereby gradually increasing from NE to SW overall (Fig. 5).

Systematic deviation revision of different batches of data
Six months after the completion of this study, three gas test holes were constructed in a mining area, and some test data were obtained (Table 3). These newly measured data are called "current measured data". The "early measured data" mentioned above are only the data of the thicker coal seams of Nos. 16 and 18; whereas other thin coal seams lack data. The purpose of this test is to obtain the gas pressure of the thin coal seam. Therefore, the current measured data include not only the gas pressure data of the thin coal seam but also the gas pressure data of the coal seams of Nos. 16 and 18. In accordance with these test data, the prediction results can be verified and compared.
In accordance with the "early measured data" of gas pressure in the coal seams Nos. 16 and 18, a contour map of gas pressure was drawn, and test holes Z1-Z3 are superimposed in Fig. 6. The early measured gas pressure values of holes Z1-Z3 could be obtained from the figure. However, a certain deviation exists between the values  (Table 4).
As shown in Table 4, the values between the previous and later measurements are inconsistent, and the deviation rate is discrete. Determining which of the data are more reliable is impossible because they are all measured values at different times; thus, the maximum and minimum values of the deviation rate are eliminated, and then the average value of the middle values of the deviation rate is regarded as the system error. In this manner, the system error of the current measured value is considered 0.22 MPa. Therefore, the current measured value needs to be corrected. In other words, 0.22 MPa should be subtracted from the current measured value to match the early measured value.

Discussion on prediction results
In accordance with the above adjustment principles, the predicted gas pressure values of the coal seams of Nos. 22, 24, and 27 were verified and compared (Table 5). As shown in Table 5, the deviation rate of the prediction results of MLR ranges from 6.5 to 19.7%, with an average of 13.0%, with the correction values of the measured values taken as the verification standard.
In addition, among the six measuring points, the prediction effect of measuring point No. 2 (coal seam No.22 of borehole Z2) is poor, and the reasons need to be discussed further. The predicted values of other measuring points are close to the measured values. The average deviation rate between the predicted and measured values is  11.6%, except for measuring point No. 2; this result is acceptable in actual production. The research in gas prediction has two trends, one is to select different prediction methods, and the other is to choose different influencing factor indexes in the research. The main method different from this paper is artificial neural network (ANN) (Wei et al. 2009;Li et al. 2019). Although this method have its advantages, it first needs to learn and train the model and ensure the convergence of the training process, which requires enough data for learning, and even needs to optimize the parameters of the model. The use of ANN has also been attempted in the research of this case, but its effect is not as good as that of MLRs, so it is not adopted.
As for the influencing factor indicators, the buried depth of coal seam is generally recognized as one of the main influencing factors, and other influencing factor indicators vary from researchers to researchers, such as the roof and floor lithology of the seam, coal physical indicators, coal quality indicators, and geophysical detection signal data. Of course, these data are more or less related to coal seam gas content. However, these indicators are not necessarily easy to obtain, especially the physical or chemical indicators of coal, which need to be tested. In addition, geophysical signal data are actually the indirect response of the object of relevant influencing factors. The indicators used in this case covered the main influencing factors of coal seam gas, which can be easily obtained directly or through calculation. Hence, it is conducive to the popularization and application of this method.
In summary, the MLR method has practical application value in the prediction of gas pressure in thin coal seams of coal mines. Given that the appropriate and quantitative influencing factors are selected, a satisfactory prediction effect can be achieved.

Conclusions
(1) Gas occurrence in coal seam is influenced by many factors, such as burial depth, coal thickness, coal structure, and geological structure, and the quantitative indexes of which can be used to predict gas pressure.
(2) The structural condition of coal can be expressed by the structural index of coal, which can reflect the degree of coal damage. The structural index of coal can be evaluated synthetically by using DLW, HGG, and drilling diameter change.
(3) Geological structure can be expressed by the development degree of faults and folds. For the structural deformation of multiple coal seams, the structural curvature can be calculated in accordance with the contour line of the coal seam floor, which can describe the structural development degree of different coal seams objectively and quantitatively. (4) In the study of gas pressure prediction of several thin coal seams in the Qinglong Coal Mine in Guizhou Province, the deviation rate of the prediction results of MLR ranges from 6.5 to 19.7%, with an average of 13.0%. Except for that of measuring point No. 2, the average deviation rate between the predicted and measured values of the five other measuring points is 11.6%, which is acceptable in actual production. This result shows that the MLR method has practical value in gas pressure prediction in thin coal seams without measured data.

Conflict of interest
No conflict of interest exits in the submission of this manuscript, and manuscript is approved by all authors for publication.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.