Abstract
Gas occurrence in coal seams shows zonation, and the outburst risk prediction index is significantly affected by coal seam gas. For different coal mine, the sensitivity of prediction index of outburst risk is inconsistent, and the different area of the same coal mine can present different sensitivities to the same index. In this study, we measured the industrial analysis of coals and determined the degree of outburst risk of each coal seam using coal samples with the particle size of 1 ∼ 3 mm. The relationships between the equilibrium gas pressure and gas desorption index of drill cuttings, between the desorption index of drill cuttings and the adsorption pressure, and between the critical values of the indexes of drill cuttings were also analyzed using the fuzzy clustering method. Based on the analysis of the particle size distribution of coal samples, the influence of grain size on the critical value of the drill cuttings desorption index was studied. The results showed that the critical values of the desorption index of drill cuttings were different under different adsorption equilibrium gas pressures, and those critical values increased with increasing adsorption equilibrium gas pressure. The desorption index of drill cuttings and adsorption equilibrium gas pressure had a linear relationship, and the higher the degree of outburst risk was, the greater the slope of the fitting function and the smaller the intercept; under the same gas pressure of adsorption equilibrium, there was a linear relationship between Δh2 and
Keywords
Introduction
Coal is the main nonrenewable energy resource consumed in China. Due to the advancement of coal mining in recent years, mining depths have reached 1,300 m and are estimated to reach 1,500 m within the next 20 years (Fan et al., 2019a; Hasiah et al., 2013; Levine, 1986). As mining depth increases, the
A number of scholars have developed sensitivity indexes to predict outburst-prone coal seams through research on the desorption properties of coal. For instance, the variation
Gas desorption indexes.
These sensitivity indexes play a very important role in outburst risk prevention. In China, after years of application and verification, it is found that the outburst risk prediction index and its critical value are inconsistent with the actual outburst risk of the mine, which is the phenomenon of ‘outburst occurrence with low index value’ (Liu et al., 2018; Nandi and Walker, 1975; Yang et al., 2018). Cheng et al. (2016) established the relationship between the gas desorption index drill cuttings Δh2 and gas pressure in coal seams; Sun and Sun (2000) conducted a similar simulation study on the
The above studies are outburst risk predictions with single-index. Due to the uncertainty of individual indexes, it is difficult to obtain accurate outburst risk judgments. To avoid the occurrence of ‘outburst occurrence with low index value’, multi-index prediction should be adopted as much as possible (Wang et al., 2018). Practice shows that different geology predictions may be different for different geological conditions. Even for the same working face under different conditions of mining technology, there may be different prediction sensitivity indexes; conversely, different prediction indexes result from different geological and mining process conditions. The working face has different sensitivities whether the predictors used are sensitive to determining the accuracy of the forecasting risk. Therefore, research on the sensitivity of outburst risk prediction indexes has become a very important technical problem in the process of outburst prediction.
In this paper, the sensitivities and critical values of
Experimental research
Experimental principle
Generally, the outburst risk prediction for a mining face is performed before mining, so the regional outburst risk prediction method (including gas pressure
The determining principle of the outburst risk prediction index of the working face is as follows: Coal-containing gas was instantly exposed to the atmosphere or instrument similar to atmospheric environmental conditions, according to the volume or variable volume pressure desorption principle determined at different times of gas desorption quantity or different times of gas desorption rate. According to the relationship between the measured data and the exposure time of the coal samples, the corresponding mathematical treatment was carried out to obtain the gas desorption index of drill cuttings.
When using the constant volume gas desorption instrument or the variable volume gas desorption instrument, indications on the instrument are the pressure values, which need to be converted into the desorption amount according to the coefficient provided by the instrument. The conversion formula (The Chinese Academy of Sciences Institute of Mathematical Statistics, 2008) is:
The gas desorption law of the coal sample obeys the following formula:
Because
When using a variable volume and variable pressure gas desorption instrument for measurement, the gas desorption values at the times
Experimental equipment
The definition of a gas desorption instrument is as follows: an electronic or physical measuring instrument that is based on an equal volume and variable capacity pressure desorption principle.
In this paper, according to the Chinese standard “determination method of gas desorption index of drill cuttings” (AQ/T 1065–2008) (The Chinese Academy of Sciences Institute of Mathematical Statistics, 2008), a device for measuring the gas desorption index of drill cuttings was designed.
A schematic diagram of the test device is shown in Figure 1, which mainly includes a charging system, a high-pressure vessel, a vacuum system, and a desorption system. The high-pressure vessel was connected to the desorption system through check valves. The compression release valve was installed on the charging system to achieve the specified gas pressure.

Schematic diagram of the test system for determining the gas desorption index of drill cuttings.
The application of both an equal volume and variable volume gas desorption data acquisition system is the most significant feature of the experimental equipment. To obtain enough pressure data during adsorption/desorption, the time interval of desorption pressure recording was reduced to approximately 10 seconds. In general, the number of data measurements was unrestricted, and the final gas pressure was almost reduced to the level of atmospheric pressure. To ensure the accuracy of the experimental data, the high-pressure vessel was placed in the constant temperature water bath, and the temperature setting was the same as the temperature of the coal seam.
Test instrument requirements:
Constant volumetric gas desorption instrument (for determination of “ Variable capacitance type gas desorption instrument (for determination of Δh2) with pressure ranges of 0 kPa–3 kPa, error ±2%. Vacuum pump with a flow rate of 1000 ml/s, and the vacuum limit was not greater than 0.1 Pa. Coal sample tank with a volume of 300 ml, coal sample cup with a volume of 8.6 ml (wall thickness should be no more than 2 mm), connection pipe space volume of 30 ml, and the volumetric error should not exceed 5%. When the coal sample tank should bear a gas pressure of 6 MPa or 4 Pa, the pressure drop in 30 h should be less than 1% under the condition that the coal sample tank did not hold the test sample. The pressure gauge range was 6 MPa, and the accuracy was 0.4. The tested gas methane concentration should not be less than 99.9%. The electronic balance sensibility was 0.1 g.
According to the basic requirements of the above design, the experimental system for
Coal sample preparation
The coal samples were obtained from the Fuxiang mine (FX) of Shandong Energy Guizhou Co., Ltd., which is located in Bijie City, Guizhou Province. The geographical location of the FX is shown in Figure 2. The to-be-mined coal seams in this coal mine are abundant, and the geological occurrence is common. The possibilities of outburst occurrence are different; thus, the possibility of outburst proneness in this coal mine cannot be accurately predicted by a single method. Therefore, a combination of laboratory measurements and field tests were used in this study. As shown in Figure 3, the coal samples were collected from coal seams #6, #8 and #13 and marked as N-6, N-8 and N-13, respectively. Each coal sample was taken from the same site, and the temperatures of the above three coal seams were also measured.

Geographical location of the FX mine.

Coal seam histogram.
Sampling and sample preparation methods (flow chart shown in Figure 4) were as follows:

Sampling and sample preparation process.
Sampling method-relevant regulations were carried out from “coal and rock physical and mechanical properties determination method first parts: sampling General Provisions” (GB/T 23561.1-2009) (The People's Republic of China Ministry of Machine Building, 2009a); on-site sampling of #6, #8, and #13 coal seams produced samples of more than 1 kg.
According to the relevant provisions of the “preparation method of coal sample” (GB/T474-2001) (The People's Republic of China Ministry of Machine Building, 2008), the sample preparation method first broke the coal sample, and then the broken coal samples were screened between 1 mm and 3 mm; the mass of the test coal sample was not less than 20 g, and the particle size was 1 mm-3 mm.
To analyze the relationship between the size of the coal sample and the critical value of the desorption index of drill cuttings, the particle size distribution of the 1-3 mm coal sample was analyzed using a particle size analyzer.
Experimental methods
According to the experimental objective, the experiment could be divided into three groups: basic coal property characterization (group I), explosive index determination (group II) and coal adsorption performance characterization (group III).
Experimental results
Proximate analysis
The results of the proximate analysis are shown in Table 2. Table 2 shows that all the coal samples had the same coal rank, and the true and apparent relative densities presented minor differences. Moreover,
Results of the proximate analysis.
Outburst risk degree classification of coal seams
The results of gas pressure determination, coal strength coefficient measurement and initial gas diffusion velocity measurement are displayed in Table 3.
Results of gas pressure determination (
Note: where
Table 3 shows that different coal seams (N-6, N-8 and N-13) had various outburst risk degrees, namely, having different values of gas pressure, coal strength coefficient and initial gas diffusion velocity of coal. In addition, the values of gas pressure could be ranked as N-6 > N-8 > N-13, the values of the coal strength coefficient could be ranked as N-13 > N-8 > N-6, and the values of the initial gas diffusion velocity of coal could be ranked as N-6 > N-8 > N-13. Generally, the gas pressure represented the amount of gas stored in the coal seam, the coal strength coefficient represented the strength of the coal, and the initial gas diffusion velocity of coal represented the gas diffusion capacity of coal during the initial desorption period. Therefore, the coal seam N-6 had the largest amount of stored gas, lowest coal strength and strongest gas diffusion capacity at the initial desorption period. Coal seam N-13 had the lowest amount of stored gas, largest coal strength and weakest gas diffusion capacity at the initial desorption period of.
To classify the outburst risk degree of the coal seams according to the “Provision of prevention and control of coal and gas outburst” issued by China, the outburst risk degree classification criteria are summarized and listed in Table 4. In accordance with the abovementioned classification criteria, the outburst risk degrees of N-6, N-8 and N-13 were obtained and are displayed in Table 3. Table 3 shows that the outburst degree of the three coal seams could be ranked as N-6 > N-8 > N-13. In addition, the selection of the target coal mine was reasonable because the coal samples from the same coal mine had three different outburst risk degrees.
Outburst degree classification criteria of coal seams.
The threshold value of drill cuttings desorption
Judging from the outburst risk division of a coal seam, the greater the gas pressure was, the greater the outburst risk dgree of a coal seam; the gas pressure of a coal seam with 0.74 MPa represented the critical pressure for the coal seam with outburst risk, while the drill cuttings desorption index could reflect the residual gas pressure in the coal seam. When the adsorption equilibrium gas pressure was 0.74 MPa, the desorption index of drill cuttings was the critical value for the evaluation of the mining face.
The critical values of the desorption index of laboratory cuttings in three locations (N-6, N-8 and N-13) in FX were determined, and the results were obtained from the average values of those tests. Table 5 is the critical value of the desorption index of drill cuttings under different adsorption equilibrium gas pressures.
Under different gas pressure values, the critical values of
Table 5 shows that the measured values of N-6 were different under different adsorption equilibrium gas pressures, and the values increased with increasing adsorption equilibrium gas pressure. The larger that the gas pressure of adsorption equilibrium was, the larger the value would be and vice versa. Under different adsorption equilibrium gas pressures, the Δh2 value was different, and the Δh2 value increased with increasing adsorption equilibrium gas pressure. The larger the adsorption equilibrium gas pressure was, the greater the value of Δh2 and vice versa. N-8 and N-13 follow the same rule as N-6, which means that, through the test data, the critical values of desorption indexes
Meanwhile, by comparing the test data of the desorption index of drill cuttings in different coal seams, it was found that the critical values were different under the same adsorption gas pressure; the desorption indexes of drill cuttings
Discussion
Relationship between critical value and gas pressure
From the comprehensive mechanism of coal and gas outbursts, it is known that the power of coal and gas outbursts mainly be produced by gas pressure. For the conditions of
Thresholds of four indicators for the identification of outburst-prone coal seams.
The critical value of a single index for assessing the outburst risk of coal seams is given in Table 6, in which the failure type of coals was an objective evaluation index. The initial velocity of gas emission and the coefficient of coal strength were determined by experiment; accordingly, the accuracy was high. The gas pressure in coal seams needed to be measured in the field, but there are many related problems, such as the duration of the process, the facilities affected by the quality of the borehole sealing and data inaccuracy. In contrast, the desorption index of drill cuttings had the characteristics of short measuring time and accurate results, and it was the key parameter for predicting coal and gas outburst risk in the mining face. This could reflect the gas pressure of the mining face to a certain extent. Therefore, the residual gas pressure in the working face could be deduced accurately by finding the relationship between the critical value of the desorption index of drill cuttings and the gas pressure. At the same time, the risk degree of the mining face could be accurately determined.
In addition, the “coal and gas outburst prevention regulations” (2009) (State Administration of Work Safety of China and State Administration of Coal Mine Safety of China, 2009) give the critical value of drill cuttings desorption index prediction, as shown in Table 7.
Reference critical values of the outburst risk prediction index.
The critical reference values of the outburst risk prediction index given in Table 7 showed that when the coal seam gas pressure was 0.74 MPa, the desorption index of drill cuttings in the corresponding mining face was 200 Pa and 0.5 mL/(g·min0.5). In view of the complexity of coal seam gas occurrence conditions, the metamorphic degree and adsorption/desorption performance of coal were quite different. “Prevention and control regulations on coal and gas outburst” (2009) (State Administration of Work Safety of China and State Administration of Coal Mine Safety of China, 2009) stipulates that the critical risk value of each index should be determined according to the measured data from each mine. Therefore, the critical value of the coal seam gas pressure outburst index and the critical value of the gas desorption index of drill cuttings outburst prediction might not have a one-to-one correspondence. When the gas pressure of the coal seam was 0.74 MPa, the measured values of desorption indexes Δh2 and

The critical values of indexes
Figure 5 shows the measured points
From Figure 5, it can be seen that the critical values of
According to the Langmuir monolayer adsorption theory, the greater the gas pressure was, the greater the amount of gas adsorbed by coal; similarly, the amount of gas will be larger during desorption. The gas desorption index of drill cuttings characterizes the gas desorption and release ability of coal in unit time. Therefore, the greater the adsorption equilibrium gas pressure was, the greater the desorption index of drill cuttings.
The trend values of
Regression analysis of the gas desorption index of drill cuttings under different adsorption equilibrium pressures.
From Table 8, it can be seen that the gas desorption index of drill cuttings was a linear function with the gas pressure of adsorption equilibrium, and the slopes of the fitted curves were all different. Concerning the slopes of
First, this showed that there was a linear function between the gas desorption index and the equilibrium gas pressure in the coal seams, which was inconsistent with the power function relationship between
According to the theory of gas adsorption and desorption, the gas in coal exists in the adsorbed state and free state. The amount of adsorbed gas is affected by pore characteristics and pore pressure, and the amount of free gas is affected by the pore volume of the coal body, the coal seam temperature and the pore pressure (Xue et al., 2020). Adsorbed gas is deposited on the micropore surface and particle structure, and most of the gas will be adsorbed in the coal body. The process of coal adsorption/desorption is also a percolation diffusion process (Black, 2011). The measurement process of the gas desorption index of drill cuttings is similar to the adsorption equilibrium gas desorption diffusion process, which includes not only the seepage of pore free gas but also the desorption and diffusion of adsorbed gas. The slope of the fitting curve of the desorption index of drill cuttings reflects the desorption rate of coal, which is the direct embodiment of gas occurrence characteristics of coal; the intercept reflects the pore feature of coal seams and the pore volume of coal bodies and is an inherent attribute of coal seams. Compared with nonoutburst coal, the characteristic of outburst coal in gas desorption is that the gas desorption velocity is high at the beginning and decreases rapidly with time. Therefore, the critical value of the desorption index of N-6 drill cuttings is lower than that of nonoutburst N-13; that is, the greater the slope of the trend fitting function is, the higher the risk of coal seam outburst. This is consistent with the assessment results of outburst risk of N-6, N-8 and N-13.
The analysis of the critical values of

Relationship between the critical values of indexes
According to Figure 6(a) to (c), the critical values of the Δh2 index of each seam increased linearly with increasing critical value of the
The trends of the
Regression analysis between
As shown in Table 9, K1 of each coal seam and the trend of Δh2 showed a linear function. The slope of the fitted curve was greater than 1, and the slope of the fitted curve was as follows: N-6 < N-8 < N-13; the intercept values were as follows: N-6 > N-8 > N-13. The correlation coefficients of each fitting function were all greater than 0.98, and the trends of the fitting curves were in good agreement with the measured data. First, the critical value of Δh2 increased more with increasing rate of gas pressure than the critical value of
Judging from the outstanding risk assessment results of N-6, N-8 and N-13, it could be concluded that the higher the risk of coal seam outburst was, the smaller the slope of the fitting trend of the desorption index of drill cuttings, and therefore, the greater the slope of the fitting curve of the desorption index of coal seam cuttings was, the lower the risk of coal seam outburst. The higher the risk of coal seam outburst was, the greater the intercept value of the fitting trend line of the desorption index of drill cuttings. Therefore, the smaller the intercept curve of the desorption index of coal seam drill cuttings was, the lower the risk of coal seam outburst. Thus, coal seam cuttings desorption index values could be fitted with trend line slopes and intercepts to determine the degree of risk of coal outbursts. The specific coal seam could be determined according to the measured value, which provided a new technical means for the risk assessment of coal seam outbursts.
Index sensitivity and critical value
Based on the measured data, the sensitivities and critical values of the desorption index of drill cuttings of coal seams N-6, N-8 and N-13 were further determined by the fuzzy clustering method.
The fuzzy clustering method is a method based on the measured data, through the establishment of a mathematical model of predictive sensitive indicators, to determine the abstract sensitivity index with specific numerical values (Beamish and Crosdale, 1998; Ruppel et al., 1972). Its steps are as follows:
(1) All measured data are nondimensional
(2) All data are standardized
(3) Sensitivity function of the prediction index for coal and gas outburst
Taking
Establish the sensitivity function of the forecast index. Differences in forecasting index analysis are different from those in mathematical statistics. The deviation in mathematical statistics reflects the degree of deviation between each discrete index and its mathematical expectation. In the forecasting index analysis, the discrepancy was reflected between each pair of forecasting indexes and the critical index. Because the deviation has positive and negative points, the average of the sum of squares of each index deviation can reflect the discrete degree of the index as a whole. If the descriptive index is defined as
The sensitivity function of the index is defined as follows:
The sensitivity functions
The sensitivity of the gas desorption index of drill cuttings in each coal seam in Table 4 is calculated using (5)–(9). The results are shown in Table 10.
Sensitivity of gas desorption index of drill cuttings.
From Table 10, the sensitivity calculation results of the gas desorption index of drill cuttings showed that the calculation results for the coal seams indicate
The sensitivity histogram of the critical values of the desorption index of drill cuttings in coal seams N-6, N-8 and N-13 are shown in Figure 7.

Sensitivity histogram of the desorption index of drill cuttings in each coal seam.
It can be seen from Figure 7 that
When the N-6 adsorption equilibrium gas pressures were 0.741 MPa, 0.742 MPa and 0.739 MPa, the values of Δh2 were less than the reference critical value of 200 Pa given in Table 4; they were 180 Pa, 187 Pa and 179 Pa, respectively. Similarly, the values of
If N-6 used the critical value of the desorption index given in Table 6 as the prediction index of the outburst risk of the mining face, the true coal seam gas pressure range would be 0.79 MPa-0.81 MPa > critical gas pressure 0.74 MPa. That is, when the measured cuttings desorption index was 188 Pa < Δh2 reference critical value 200 Pa, 0.47 ml/(g
As shown in Figure 8, the critical values of Δh2 and

Critical value of the desorption index of drill cuttings in each coal seam.
In summary, the maximum critical value of Δh2 in N-6, N-8 and N-13 was 188 Pa, and the maximum critical value of
Particle size distribution of sensitivity indicators
The desorption index of coal samples reflects the law of releasing gas after the exposure of underground coal (Zhou et al., 2020). The gas desorption index of drill cuttings reflects the gas pressure of coal seams, while different indexes show different sensitivities and critical values due to the zonation characteristics of gas occurrence in coal seams. The desorption index of drill cuttings value of

Analysis of the particle size distribution of coal particles with 1–3 mm diameter, where (a) is N-6, (b) is N-8, and (c) is N-13.
As shown in Figure 9(a), it can be seen that for N-6, particle sizes of 2092 μm–2637 μm accounted for 79.9% of the total, while minerals with particle sizes of 2334 μm–2637 μm accounted for 46.49% of the total. From Figure 9(b), it can be seen for N-8 that particle sizes of 1934 μm–2635 μm accounted for 71.54% of the total, while minerals with a grain size of 2170 μm–2635 μm accounted for 53.86% of the total. From Figure 9(c), it can be seen for N-13 that particle sizes of 2034 μm–2834 μm accounted for 78.25% of the total, while minerals with particle sizes of 2297 μm–2834 μm accounted for 58.47% of the total. From the overall trend, the mineral particle distributions of N-6 and N-13 are relatively similar. The mineral particles of N-8 with dimensions less than 2000 μm are slightly larger than those of N-6 and N-13. Under the same sieving condition, more coal particles would be easier to crush, or in other words, the coefficient of solidity of the coal would be smaller. According to the single index of outburst risk evaluation, the smaller the coefficient of solidity of coal was, the higher the risk of coal seam outburst.
The results show that under the same pressure conditions, the smaller the particle size of coal samples were, the greater the desorption index value of drill cuttings (Beamish and Crosdale, 1998; Zhou et al., 2020). According to the Langmuir monolayer adsorption model, the greater the specific surface area of coal was, the larger the amount of gas adsorbed under the same temperature and pressure conditions (Liu et al., 2019; Toda and Toyoda, 1972; Fan et al., 2019d) and the smaller the amount of desorption gas. Under the same adsorption equilibrium gas pressure condition, the critical values of
Field verification results of the desorption index of drill cuttings
When using the desorption index of drill cuttings to investigate the effect of the regional measurements on the 12619 (N-6 coal seam) mining face of FX, 3 holes with a diameter of 42 mm and a hole depth of 8–10 m were drilled into the front of the working face. First, the borehole should be drilled in the soft coal seam as much as possible. Second, at least one of the boreholes should be located in the middle of the cross-section of the excavation alley, and it also should be parallel to the direction of the heading face. Third, the final position of other boreholes should be outside the contour line on both sides of the alley section at 2–4 m. Moreover, the indexes

Determined drilling arrangement diagram.
The gas desorption index Δh2 of N-6 was Δh2 = 88.15 + 127.39

Desorption index of actual drill cuttings.
As shown in Figure 11, a total of 70 Δh2 values were determined during the preparation of the mining face, and the range of the Δh2 values of the measured gas desorption index of drill cuttings was 152∼221 Pa. In the actual determination of Δh2, there were 8 data values greater than 182 Pa, and the rest were less than 182 Pa. According to the relevant data of N-6, there were 7 gas dynamic phenomena occurring in the 8 boreholes of the actual Δh2 greater than 182 Pa. The reliability of Δh2 was high, and it was reasonable to predict the gas outburst risk of each coal seam with the desorption index Δh2 of drill cuttings. From the determined requirements and the process of finding the gas desorption index Δh2 of drill cuttings, the result was influenced by the exposure time of the coal sample, the type of coal destruction, the size of the coal sample and the sealing property of the instrument. In terms of the exposure time of coal samples, the gas desorption indexes of drill cuttings were usually measured using the orifice sampling method in the field. Therefore, it cannot be guaranteed that all the drill cuttings are taken from a predetermined depth, causing different gas desorption times for the coal samples and leading to errors when measuring the gas desorption index of drill cuttings. In addition, the operational skills of coal mine technicians also affect the exposure time.
Conclusions
In this paper, through on-site sampling, a laboratory-designed experimental device using the isobaric pressure desorption principle is used to measure the gas desorption index of coal samples in different time periods, and the test results have been theoretically analyzed and discussed. Finally, the following conclusions have been drawn:
Under the conditions of different adsorption equilibrium gas pressures, the critical values of the desorption index of drill cuttings were different, and the critical value increased with increasing adsorption equilibrium gas pressure. The trend slope of Δh2 was greater than the trend slope of The fuzzy clustering method was used to analyze the measured data. From the numerical calculation, it was concluded that Δh2 was more sensitive than Through repeated tests, the mathematical expectation was used to determine the critical value of the desorption index of drill cuttings of N-6, N-8 and N-13. The sensitivity indexes Δh2 were 182 Pa, 181 Pa, and 188 Pa, and the auxiliary predictive indexes Coal seams with different outburst risk levels had different sensitivities to the desorption index of drill cuttings; the sensitivity of different coal seams to the desorption indexes of the same drill cuttings was different for the same mine; the sensitivity and critical value of the desorption index of drill cuttings in different coal seams were also different.
Footnotes
Acknowledgements
The authors also thank the editor and anonymous reviewers very much for their valuable advice.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work is financially supported by the National Science and Technology Major Project of China (Grant No. 2016ZX0504004), the State Key Research Development Program of China (Grant No. 2016YFC0801404 and 2016YFC0801402), the National Natural Science Foundation of China (51674050), the mechanism of coal bed gas desorption and induced dynamic disasters in Guizhou Karst area (KY[2015]353) and the experimental study on coal bed adsorption capacity and improvement of gas recovery rate (LH[2015]7107), which are gratefully acknowledged.
