A Novel Permeability Prediction Model for Deep Coal via NMR and Fractal Theory

: A quantitative description of the permeability of deep coal is of great signiﬁcance for improving coalbed methane recovery and preventing gas disasters. The Schlumberger–Doll research (SDR) model is often used to predict rock permeability, but it has inherent defects in characterizing the pore structure of deep coal. In this study, a permeability model with fractal characteristics (FCP model) is established for deep coal based on nuclear magnetic resonance (NMR). The constants in the SDR model are theoretically explained by the relevant parameters in the FCP model. Centrifugation and NMR experiments were performed to determine the optimal centrifugal force and dual T 2 cutoff values. The results show that the coal samples are mainly composed of micrometer and nanometer pores. The adsorption pores account for the largest proportion, followed by the percolation pores and migration pores. In addition, the prediction accuracy of the FCP model is signiﬁcantly higher than that of the other three models, which provides a fast and effective method for the evaluation of deep coal permeability.


Introduction
Due to the decrease in shallow coal resources in China, deep mining is gradually becoming normal. However, deep coal is in a high-stress and high geothermal environment, and the permeability of the coal seam decreases sharply [1,2]. This situation will lead to frequent accidents and difficulty in achieving safe and efficient mining. The magnitude and spatial distribution of permeability represent the ability of the coal seam to allow fluid flow and output [3]. In the mining of deep coal resources, accurate acquisition of coal seam permeability plays an important role in preventing water inrush and preventing gas outbursts.
NMR technology can be effectively applied to the rock strata with complex lithology to accurately measure the permeability of rock mass [4,5]. At present, a variety of NMR permeability prediction models have been established and applied. The SDR model is a classical NMR permeability prediction model. The SDR model mainly estimates permeability based on the T 2g [6,7]. However, as an empirical formula, the SDR model still has many shortcomings. To improve the SDR model, Zhou found that the geometric mean value of T 2 could better characterize the permeability of shale and obtained a single-parameter NMR permeability model [8]. To further reflect the influence of saturated and bound fluids on permeability, a permeability model is proposed based on the dual T 2 cutoff values to better predict shale permeability [9].
The Timur-Coates (Coates) model uses the T 2 cutoff value to determine the distribution of fluid in the pores and porosity to predict permeability [10]. Many scholars have improved the Coates model because the T 2 cutoff values of different rocks are complex and difficult to accurately determine. A permeability model based on dual T 2 cutoff values is established. The model divides the fluid in tight sandstone into three parts to further describe the seepage characteristics of the fluid [11]. The pore connectivity factor was obtained by displacing the fully saturated carbonate core with heavy water. The Coates model was modified according to pore connectivity [12].
In addition to improvements or innovations based on the SDR model or Coates model, some scholars have proposed composite methods to predict permeability. The permeability is predicted by transverse surface relaxivity, logarithmic mean of transverse relaxation time, and formation factor. The obtained results are more accurate than the classical permeability model [13]. The partial least square method is used to determine the relationship between transverse surface relaxivity segments and permeability, and permeability can be predicted [14]. Smith proposed the use of NMR responses to provide equations without any external input, with good results in shale and carbonate reservoirs [15]. The pore network geometry and permeability of the fractured rock are analyzed using OpenFOAM as the main flow solver [16][17][18].
At present, the SDR model is widely used in the prediction of rock permeability, but it still faces the following key problems.
(1) SDR model lacks theoretical support. Whether the constants in the SDR model have practical significance needs to be studied. (2) The SDR model is fitted by measuring the permeability of a large number of rocks, which is uneconomical and time-consuming. (3) The pore structure is assumed to be homogeneous in the SDR model, but it is complex and heterogeneous in deep coal. The application of the SDR model is limited. The equation of the SDR model is as follows [6]: (1) where m and n are constants in the model.
The research on the theoretical improvement of the SDR model is not enough [19,20], and these problems need further study. In this study, a permeability prediction model for deep coal was established to solve the above problems. The FCP model provides theoretical support for the SDR model. The physical meaning of the constants in the SDR model is explained by the relevant parameters in the FCP model. Meanwhile, the FCP model focuses on predicting the permeability of tight rocks with low porosity. According to the experimental results, the pore size distribution of the samples was determined, and the evolution law of pore water content was analyzed. The predicted permeability value of the FCP model is closer to the measured value. Compared with other permeability prediction models, the accuracy of the FCP model is verified.

Principles of NMR
The basic principle of NMR technology is to measure the relaxation properties of hydrogen by arranging hydrogen nuclei under an applied magnetic field. Relaxation responses of various hydrogen-containing components (water, methane, etc.) can be obtained by specific RF scans of hydrogen nuclei. Equation (2) is the expression for the transverse relaxation time T 2 [21]: where T 2 is the transverse relaxation time (ms) and ρ 2 denotes transverse surface relaxivity (µm/s); S and V refer to the surface area (µm 2 ) and volume (µm 3 ), respectively. F s refers to a constant value dependent on geometric pore shape. F s is 2 for columnar and 3 for spherical pores, respectively. From Equation (2), there is a proportional relationship between the T 2 value and porosity, as shown in Figure 1. The larger the value of T 2 , the larger the corresponding pore radius. The pore structure of porous media can be quantitatively characterized by the T 2 spectrum [22][23][24]. where T2 is the transverse relaxation time (ms) and ρ2 denotes transverse surface relaxivity (μm/s); S and V refer to the surface area (μm 2 ) and volume (μm 3 ), respectively. Fs refers to a constant value dependent on geometric pore shape. Fs is 2 for columnar and 3 for spherical pores, respectively. From Equation (2), there is a proportional relationship between the T2 value and porosity, as shown in Figure 1. The larger the value of T2, the larger the corresponding pore radius. The pore structure of porous media can be quantitatively characterized by the T2 spectrum [22][23][24].

Permeability Expression
Permeability reflects the capacity of a rock to allow fluid to pass; the connected pores are the prerequisite for seepage. qi is the flow rate of the capillary with radius ri by the Hagen-Poiseulle equation [25]: where vi is the volume of capillary i, vi = πri 2 L0, μ is the viscosity of the fluid, Δp is the pressure difference, and L0 is the actual length of the capillary.
The capillary with different pore sizes in Figure 2 reflects the pore structure of the rock. According to the fractal geometry theory, the cumulative size of pores in porous media is expressed as follows [26]: where N is the cumulative number of pores, ε is the length scale, Df is the fractal dimension, and rmax is the maximum pore radius.

Permeability Expression
Permeability reflects the capacity of a rock to allow fluid to pass; the connected pores are the prerequisite for seepage. q i is the flow rate of the capillary with radius r i by the Hagen-Poiseulle equation [25]: where v i is the volume of capillary i, v i = πr i 2 L 0 , µ is the viscosity of the fluid, ∆p is the pressure difference, and L 0 is the actual length of the capillary.
The capillary with different pore sizes in Figure 2 reflects the pore structure of the rock. According to the fractal geometry theory, the cumulative size of pores in porous media is expressed as follows [26]: where N is the cumulative number of pores, ε is the length scale, D f is the fractal dimension, and r max is the maximum pore radius.
where T2 is the transverse relaxation time (ms) and ρ2 denotes transverse surface relaxivity (μm/s); S and V refer to the surface area (μm 2 ) and volume (μm 3 ), respectively. Fs refers to a constant value dependent on geometric pore shape. Fs is 2 for columnar and 3 for spherical pores, respectively. From Equation (2), there is a proportional relationship between the T2 value and porosity, as shown in Figure 1. The larger the value of T2, the larger the corresponding pore radius. The pore structure of porous media can be quantitatively characterized by the T2 spectrum [22][23][24].

Permeability Expression
Permeability reflects the capacity of a rock to allow fluid to pass; the connected pores are the prerequisite for seepage. qi is the flow rate of the capillary with radius ri by the Hagen-Poiseulle equation [25]: where vi is the volume of capillary i, vi = πri 2 L0, μ is the viscosity of the fluid, Δp is the pressure difference, and L0 is the actual length of the capillary.
The capillary with different pore sizes in Figure 2 reflects the pore structure of the rock. According to the fractal geometry theory, the cumulative size of pores in porous media is expressed as follows [26]: where N is the cumulative number of pores, ε is the length scale, Df is the fractal dimension, and rmax is the maximum pore radius.  From Equation (4), the number of pores distributed between r and r + dr is [26]: The total flow rate Q can be expressed as: The expression for v i is: where η i is the fraction of capillary i in the pore volume, A is the cross-sectional area, and L is the representative length. The total flow Q is expressed as: The following equation is the expression of dη i .
The relation between porosity and D f is as follows [26]: where ϕ is porosity, D is the Euclidean dimension of the space, and D = 2 in two dimensions. By combining Equations (9) and (10), the total flow Q is obtained.
The average pore radius r ave can be obtained with the aid of Equation (4) [27,28].
Combined with Equation (10), the maximum pore radius r max can be expressed as: Substituting Equation (13) into Equation (11), the expression for Q is obtained.
The permeability equation is obtained: where τ is tortuosity.
T is given by the following expression [29][30][31]: The relaxation time of rock is in good agreement with its internal aperture distribution [32,33]. Based on Equation (2), r ave is represented by the geometric mean of T 2 .
where T 2g is the geometric mean value of T 2 , and F s is 2 for columnar. Substituting Equation (17) into Equation (15), the permeability model is: The FCP model in Equation (18) provides a theoretical explanation for the SDR model. We found that the constant m in the SDR model is closely related to the fractal dimension D f , tortuosity τ, porosity ϕ, and transverse surface relaxivity ρ, and the value of the constant n is affected by the fractal dimension D f .

Coal Samples
The raw coal, from the Pingdingshan Coal Mine in Henan Province, China, is between 1006 m and 1137 m deep. By official standards, the raw coal was processed into 8 cylindrical coal samples with a diameter of 25 mm and a length of 50 mm, as shown in Figure 3. According to the composition analysis of the X-ray fluorescence spectrum, the coal samples do not contain substances that affect the NMR results.
where τ is tortuosity. Τ is given by the following expression [29][30][31]: The relaxation time of rock is in good agreement with its internal aperture distribution [32,33]. Based on Equation (2), rave is represented by the geometric mean of T2.
where T2g is the geometric mean value of T2, and Fs is 2 for columnar. Substituting Equation (17) into Equation (15), the permeability model is: The FCP model in Equation (18) provides a theoretical explanation for the SDR model. We found that the constant m in the SDR model is closely related to the fractal dimension Df, tortuosity τ, porosity φ, and transverse surface relaxivity ρ, and the value of the constant n is affected by the fractal dimension Df.

Coal Samples
The raw coal, from the Pingdingshan Coal Mine in Henan Province, China, is between 1006 m and 1137 m deep. By official standards, the raw coal was processed into 8 cylindrical coal samples with a diameter of 25 mm and a length of 50 mm, as shown in Figure 3. According to the composition analysis of the X-ray fluorescence spectrum, the coal samples do not contain substances that affect the NMR results.

Experimental Procedures and Facilities
(1) Coal samples were weighed, and their length and diameter were measured. The permeability of the coal samples was measured by the gas injection porosimetry method. (2) Coal samples were dried at 80 • C for 24 h until the weight no longer changed. After cooling to room temperature, coal samples were placed into the vacuum saturation device. Dry pumping was carried out for 2 h with negative pressure of 0.1 MPa, then soaking in water with a pressure of 15 MPa for 24 h to make coal samples fully saturated with water. (3) An NMR experiment was conducted on the coal sample. As shown in Figure 4, the instrument model used is MacroMR12-150 H-I. The magnets used are permanent magnets with a magnetic field intensity of (0.3 ± 0.05) T, main frequency of 12.42 MHz, and temperature of (32.00 ± 0.02) • C. The echo-spacing TE = 0.1 ms, the 90 • pulse width P1 = 10 µs, the 180 • pulse width P2 = 20 µs, and the echo number NECH = 15,000. (4) To determine the optimal centrifugal force of the coal sample, coal sample F was selected to conduct centrifugal experiments with different centrifugal forces.  (5) The remaining coal samples were centrifuged by optimal centrifugal force and tested by NMR. The experimental parameters were the same as in Step (3) to obtain the NMR porosity and T 2 cutoff value. The NMR data were processed via Numai NMR data analysis software V1.0, and the T 2 spectrum was inverted by the SIRT algorithm.
(3) An NMR experiment was conducted on the coal sample. As shown in Figure 4, the instrument model used is MacroMR12-150 H-I. The magnets used are permanent magnets with a magnetic field intensity of (0.3 ± 0.05) T, main frequency of 12.42 MHz, and temperature of (32.00 ± 0.02) °C. The echo-spacing TE = 0.1 ms, the 90° pulse width P1 = 10 μs, the 180° pulse width P2 = 20 μs, and the echo number NECH = 15,000. (4) To determine the optimal centrifugal force of the coal sample, coal sample F was selected to conduct centrifugal experiments with different centrifugal forces. (5) The remaining coal samples were centrifuged by optimal centrifugal force and tested by NMR. The experimental parameters were the same as in Step (3) to obtain the NMR porosity and T2 cutoff value. The NMR data were processed via Numai NMR data analysis software V1.0, and the T2 spectrum was inverted by the SIRT algorithm.

Determination of the Optimal Centrifugal Force
The purpose of the centrifugation experiment is to obtain the centrifugal force suitable for coal samples. The optimal centrifugal force can achieve accurate calibration of the T2 cutoff value of NMR. The accurate selection of centrifugal force is very important. If the centrifugal force is small, the T2 cutoff value will be larger, the bound water saturation will be larger, and the seepage capacity of coal will also be underestimated. To obtain the best centrifugal force of deep coal, centrifugal tests were carried out on the selected sample F under the centrifugal force of 0.53 MPa, 0.72 MPa, 0.94 MPa, 1.19 MPa, and 1.46 MPa. NMR was used to measure the coal sample after each centrifugation. Figure 5 shows the NMR T2 spectra of sample F with the characteristics of three peaks. With the increase in centrifugal forces, the cumulative porosity of sample F decreases gradually. According to the pore classification method [34][35][36][37], the peak S1 (T2 < 2.5 ms) of the T2 spectrum was defined as adsorption pores. The peak S2 (2.5~100 ms) was defined as the percolation pores. The peak S3 (T2 > 100 ms) was defined as migration pores. In Figure 6, with the increase in centrifugal force, the water content of each pore decreases gradually. The water content of adsorption pores and total pores has the same trend and similar value. It is higher than percolation pores and migration pores.

Determination of the Optimal Centrifugal Force
The purpose of the centrifugation experiment is to obtain the centrifugal force suitable for coal samples. The optimal centrifugal force can achieve accurate calibration of the T 2 cutoff value of NMR. The accurate selection of centrifugal force is very important. If the centrifugal force is small, the T 2 cutoff value will be larger, the bound water saturation will be larger, and the seepage capacity of coal will also be underestimated. To obtain the best centrifugal force of deep coal, centrifugal tests were carried out on the selected sample F under the centrifugal force of 0.53 MPa, 0.72 MPa, 0.94 MPa, 1.19 MPa, and 1.46 MPa. NMR was used to measure the coal sample after each centrifugation. Figure 5 shows the NMR T 2 spectra of sample F with the characteristics of three peaks. With the increase in centrifugal forces, the cumulative porosity of sample F decreases gradually. According to the pore classification method [34][35][36][37], the peak S1 (T 2 < 2.5 ms) of the T 2 spectrum was defined as adsorption pores. The peak S2 (2.5~100 ms) was defined as the percolation pores. The peak S3 (T 2 > 100 ms) was defined as migration pores. In Figure 6, with the increase in centrifugal force, the water content of each pore decreases gradually. The water content of adsorption pores and total pores has the same trend and similar value. It is higher than percolation pores and migration pores.  In Figure 6, the water content of adsorption pores decreases greatly when the centrifugal force increases from 0.53 MPa to 1.19 MPa. When the centrifugal force increases from 1.19 MPa to 1.46 MPa, the water content of adsorption pores is less than 3%. The water content of percolation pores and migration pores does not change obviously when the centrifugal force is greater than 0.53 MPa. The expected experimental results were achieved at a centrifugal force of 1.46 MPa. In addition, Yang et al. [37] conducted centrifugal experiments on coal samples and selected 1.46 MPa as the best centrifugal force. This result is consistent with the centrifugal force we selected, indicating that the centrifugal experiment is effective. Based on 1.46 MPa centrifugal force, NMR experiments were performed on the remaining seven coal samples before and after centrifugation.  In Figure 6, the water content of adsorption pores decreases greatly when the centrifugal force increases from 0.53 MPa to 1.19 MPa. When the centrifugal force increases from 1.19 MPa to 1.46 MPa, the water content of adsorption pores is less than 3%. The water content of percolation pores and migration pores does not change obviously when the centrifugal force is greater than 0.53 MPa. The expected experimental results were achieved at a centrifugal force of 1.46 MPa. In addition, Yang et al. [37] conducted centrifugal experiments on coal samples and selected 1.46 MPa as the best centrifugal force. This result is consistent with the centrifugal force we selected, indicating that the centrifugal experiment is effective. Based on 1.46 MPa centrifugal force, NMR experiments were performed on the remaining seven coal samples before and after centrifugation. Figure 7 shows the comparison of NMR responses of samples under the saturated state and the bound water state after centrifugation. Coal samples before and after the centrifugal T2 spectra have the characteristics of the three peaks. In Figure 7, S2 and S3 belong to a whole, and their corresponding pores are well connected. S1 is independent of S2 and S3, and the connectivity between adsorption pores and other pores is poor. After centrifugation, the peak of S1 did not decrease significantly, while the peak of S2 and S3 decreased significantly. The corresponding peak values of S2 and S3 are accompanied by a trend of rightward migration. This is the performance of water flowing from partial  Figure 7 shows the comparison of NMR responses of samples under the saturated state and the bound water state after centrifugation. Coal samples before and after the centrifugal T 2 spectra have the characteristics of the three peaks. In Figure 7, S2 and S3 belong to a whole, and their corresponding pores are well connected. S1 is independent of S2 and S3, and the connectivity between adsorption pores and other pores is poor. After centrifugation, the peak of S1 did not decrease significantly, while the peak of S2 and S3 decreased significantly. The corresponding peak values of S2 and S3 are accompanied by a trend of rightward migration. This is the performance of water flowing from partial percolation pores to migration pores. This indicates that the connectivity between the percolation pores and the transport pores is higher than that of the adsorption pores.

Calibration of T 2 Cutoff Values
After centrifugation, the water content of the adsorption pores of the samples did not change obviously. Yao [38] gave the T 2 boundary between capillary water and film water, T 2 = 0.25 ms. This indicates that the film water mainly exists in the adsorption pores. The dual T 2 cutoff values method was applied to remove the effect of water film thickness on pore size. In Figure 7, the adsorption pores are divided into two parts by the first cutoff value T 2c1 . The water in the adsorption pores between 0 ms and T 2c1 is called immovable film water, and the water between T 2c1 and 0.25 ms is called partially movable film water. Adsorption pores filled with immovable film water are not considered to be involved in water seepage [38,39]. Therefore, we removed this part of adsorption pores to minimize the influence of water film thickness on pore diameter. Table 1 presents the porosity of the eight coal samples ranging from 1.67% to 3.83% with an average of 2.51%, and the permeability ranges from 0.021 × 10 −3 µm 2 to 0.272 × 10 −3 µm 2 with an average of 0.113 × 10 −3 µm 2 . The porosity and permeability of the coal samples are complicated and belong to low porosity and low permeability coal. dual T2 cutoff values method was applied to remove the effect of water film thickness o pore size. In Figure 7, the adsorption pores are divided into two parts by the first cuto value T2c1. The water in the adsorption pores between 0 ms and T2c1 is called immovab film water, and the water between T2c1 and 0.25 ms is called partially movable film wate Adsorption pores filled with immovable film water are not considered to be involved water seepage [38,39]. Therefore, we removed this part of adsorption pores to minimi the influence of water film thickness on pore diameter.  Table 1 presents the porosity of the eight coal samples ranging from 1.67% to 3.8 with an average of 2.51%, and the permeability ranges from 0.021 × 10 −3 μm 2 to 0.272 × 1 μm 2 with an average of 0.113 × 10 −3 μm 2 . The porosity and permeability of the coal samp are complicated and belong to low porosity and low permeability coal.

Distribution of Coal Pore Size
The pores are assumed to be capillaries in Section 2.2. The coal samples in this paper belong to tight cores, and the empirical value of ρ 2 is 10 µm/s. The T 2 spectra of watersaturated coal samples in Figure 8 correspond to the pore structure of samples. According to Equation (2), the T 2 spectra of samples are converted to pore radius. In Figure 8a, there are three peaks (S1, S2, S3) in the T2 spectra of coal samples. The S1 peak distribution of coal samples is similar, and the peak value is around 0.1 ms. However, the peak intensity is slightly different, and the distribution of S2 and S3 peaks is significantly different. In Figure 8b, the S1 peaks of samples B, C, and D almost coincide, indicating that their adsorption pore structures are similar. We find that the porosity, T2g, and T2c1 of the samples are indeed similar in Table 1. The difference between the samples is mainly reflected in the S2 peak and S3 peak, which represent the percolation pores and migration pores. The pore structures of samples E and H in Figure 8c also have similar characteristics.
The pore distribution frequency of coal samples and the corresponding pore structure diagram are shown in Figure 9. The spherical pores in the pore distribution diagram represent pores with different pore sizes. The pore size distribution and pore connectivity are characterized by the number and arrangement of these spherical pores. Comparing spherical pores between different samples, the complexity of the pore structure can be reflected, and the relationship between pore structure and permeability can be under- In Figure 8a, there are three peaks (S1, S2, S3) in the T 2 spectra of coal samples. The S1 peak distribution of coal samples is similar, and the peak value is around 0.1 ms. However, the peak intensity is slightly different, and the distribution of S2 and S3 peaks is significantly different. In Figure 8b, the S1 peaks of samples B, C, and D almost coincide, indicating that their adsorption pore structures are similar. We find that the porosity, T 2g, and T 2c1 of the samples are indeed similar in Table 1. The difference between the samples is mainly reflected in the S2 peak and S3 peak, which represent the percolation pores and migration pores. The pore structures of samples E and H in Figure 8c also have similar characteristics.
The pore distribution frequency of coal samples and the corresponding pore structure diagram are shown in Figure 9. The spherical pores in the pore distribution diagram represent pores with different pore sizes. The pore size distribution and pore connectivity are characterized by the number and arrangement of these spherical pores. Comparing spherical pores between different samples, the complexity of the pore structure can be reflected, and the relationship between pore structure and permeability can be understood. The samples with complex pore structures have lower permeability. From the pore network composition ratio in Figure 9, adsorption pores account for the largest proportion of total pores in the coal samples. Adsorption pores are specially developed in samples B, C, D, and F. Percolation pores are a small part of the pore system. The development of migration pores in sample G is higher than that in other coal samples. The connectivity between adsorption pores and percolation pores is lower than that between percolation pores and migration pores.
The distribution of adsorption pores in coal samples is similar, and the pore volume proportion is high based on Figures 8 and 9. The coal samples are mainly composed of a relatively uniform distribution of micrometer and nanometer pores. The distribution of percolation pores and migration pores is relatively complex. Figure 10 shows the pore water content of samples before and after centrifugation at 1.46 MPa. The results in Table 2 show that the adsorption pores volume accounts for 67.60~85.82% of the total pore volume, the percolation pores volume accounts for 3.20~19.86%, and the migration pores volume accounts for 10.98~19.59%. The proportion of pore water content shows that the adsorption pores account for most of the total volume of coal samples in Pingdingshan, followed by percolation pores and migration pores. After centrifugation, the water content distribution of the eight samples changed obviously. Residual water mainly exists in the adsorption pores. The water content in the adsorption pores of the samples increased by 6.94% at least and by 19.07% at most. The water content of percolation pores does not change. The decrease in water content in migration pores is most obvious, and the decrease range of water content is between 6.57% and 18.75%. The change in water content in pores is consistent with the conclusion by Shen [40]. The developed migration pores and percolation pores increase the drainage, while the adsorption pores mainly contribute to the residual water content.

The Permeability Predicted Results
To verify the accuracy of the FCP model for deep coal permeability prediction, it is compared with the predicted values of the other three models based on the measured permeability values of samples. of percolation pores does not change. The decrease in water content in migration pores is most obvious, and the decrease range of water content is between 6.57% and 18.75%. The change in water content in pores is consistent with the conclusion by Shen [40]. The developed migration pores and percolation pores increase the drainage, while the adsorption pores mainly contribute to the residual water content.   [10].
where s is the model constant.
The equation of the extended SDR model in coal is as follows [41,42]: where c is the model constant. Based on measured permeability, the extended Coates model is fitted, and the equation is obtained: The extended SDR model is given by the following expression: Liu obtained the permeability equation by the pore size ratio λ [43]: Table 3 shows the parameters needed to calculate the SDR model and Coates model. The permeability calculation results of the three models are shown by comparing them with the measured permeability in Figure 11a-c. The results show that the correlation coefficients with measured permeability are 0.61, 0.83, and 0.89, respectively. Liu's model has the lowest prediction accuracy, its confidence interval is larger, and the number of samples in the confidence interval is the least.
The calculation effect of the FCP model is shown in Figure 11d, and its correlation coefficient reaches 0.98. The correlation coefficient between measured permeability and predicted permeability of the FCP model is the highest, which indicates that the prediction accuracy of the FCP model is better than that of other models. The confidence interval of the FCP model is the narrowest, and the number of samples in the confidence interval is also the largest. The results show that the stability and reliability of the FCP model are higher than those of the three models. The accuracy and validity of the FCP model for deep coal permeability prediction are verified on this basis. samples in the confidence interval is the least. The calculation effect of the FCP model is shown in Figure 11d, and its correlatio coefficient reaches 0.98. The correlation coefficient between measured permeability an predicted permeability of the FCP model is the highest, which indicates that the predictio accuracy of the FCP model is better than that of other models. The confidence interval the FCP model is the narrowest, and the number of samples in the confidence interval also the largest. The results show that the stability and reliability of the FCP model a higher than those of the three models. The accuracy and validity of the FCP model f deep coal permeability prediction are verified on this basis.

Conclusions
In this study, we propose a permeability model with fractal characteristics to improv the SDR model. The parameters in the FCP model were obtained through centrifugatio and NMR experiments, and the FCP model was validated through comparison with oth models. The following conclusions were obtained: (1) A permeability prediction model suitable for deep coal is established in this pape The FCP model is an improvement of the SDR model. The constants in the SDR mod are given theoretical explanations based on the fractal dimension Df, tortuosity τ, p rosity φ, and transverse surface relaxivity ρ in the FCP model. (2) The optimum centrifugal force of the samples was determined to be 1.46 Mpa by th evolution of pore water content. After centrifugation, the adsorption pores volum accounts for 67.60~85.82% of the total pore volume, the percolation pores volum accounts for 3.20~19.86%, and the migration pores volume accounts for 10.98~19.59% This indicates that the adsorption pores account for most of the total volume of co

Conclusions
In this study, we propose a permeability model with fractal characteristics to improve the SDR model. The parameters in the FCP model were obtained through centrifugation and NMR experiments, and the FCP model was validated through comparison with other models. The following conclusions were obtained: (1) A permeability prediction model suitable for deep coal is established in this paper.
The FCP model is an improvement of the SDR model. The constants in the SDR model are given theoretical explanations based on the fractal dimension D f , tortuosity τ, porosity ϕ, and transverse surface relaxivity ρ in the FCP model.
(2) The optimum centrifugal force of the samples was determined to be 1.46 MPa by the evolution of pore water content. After centrifugation, the adsorption pores volume accounts for 67.60~85.82% of the total pore volume, the percolation pores volume accounts for 3.20~19.86%, and the migration pores volume accounts for 10.98~19.59%. This indicates that the adsorption pores account for most of the total volume of coal samples, followed by percolation pores and migration pores.