Research Article  Open Access
Guangzhe Deng, Rui Zheng, "Reconstruction of 3D Micro Pore Structure of Coal and Simulation of Its Mechanical Properties", Advances in Materials Science and Engineering, vol. 2017, Article ID 5658742, 9 pages, 2017. https://doi.org/10.1155/2017/5658742
Reconstruction of 3D Micro Pore Structure of Coal and Simulation of Its Mechanical Properties
Abstract
This article takes the low permeability coal seam in the coalfield of South Judger Basin in Xinjiang, as a research object. The pore structure characteristics of coal rock mass in low permeability coal seam were analyzed quantitatively using scanning electron microscopy (SEM) through the methods of statistics and digital image analysis. Based on the pore structure parameters and the distribution function of the coal rock mass, a threedimensional porous cylinder model with different porosity was reconstructed by FLAC3D. The numerical simulation study of reconstructed pore model shows that (1) the porosity and the compressive strength have obvious nonlinear relation and satisfy the negative exponential relation; (2) the porosity significantly affects the stress distribution; with the increase of micro porosity, the stress distribution becomes nonuniform; (3) the compressive failures of different models are mainly shear failures, and the shape of fracture section is related to porosity; (4) the variation of seepage coefficient of the pore reconstruction model is consistent with the development of micro cracks. The micro mechanism of the deformation and failure of coal and the interaction of multiphase flow with porosity are revealed, which provides a theoretical reference for the clean development of the low permeability coal seam.
1. Introduction
With the adjustment of China’s energy structure and the increasing demand for clean mining, exploring the intrinsic relationship between the coal pore structure and macro and micro mechanical properties has significant theoretical meaning for solving mining and other fields of engineering problems, such as prefracturing of coal seams, coal bed methane mining, underground coal gasification, and unconventional oil and gas resources [1–4].
Domestic and foreign scholars pay much attention to the study of the relationship between the micro structure of coal seam and its macro and micro mechanical properties. Firstly, the pore structure of the coal seam is the reservoir space and also the spatial of the gasliquid flow [5, 6], which is the main factor that affects the coal gas storage capacity and the seepage characteristics. The pore structure of coal has a direct influence on the gas adsorption and desorption characteristics [7]. Secondly, the study of the characteristics of the micro pore structure is helpful for analyzing the influence of pore structure on the crack propagation, coal bed methane, and gas outburst [8, 9], which are helpful for the prediction of coal seam gas resources and coal and gas outburst risk.
At present, the acquisition and analysis of the pore structure characteristics, the macroscopic mechanical properties, and the deformation and failure behavior of coal and rock mass mainly depend on the methods of field sampling and indoor experiment and numerical simulation. The relationship between porosity and compressive strength and shear strength of sandstone is studied [10]. Hatzor and Palchik [11] probed into the effect of the porosity of dolomite on the crack stress intensity of rock uniaxial compressive strength and fracture and gave an empirical relationship between porosity and stress intensity. The discrete element model is used to study the relationship between stress and strain, acoustic emission characteristics, strain localization, and shear reinforcement [12]. The threedimensional statistical model of rock pore structure is constructed, and the influence of the pore on the stress distribution, the failure mode, and the fracture connectivity of rock are analyzed [13]. The threedimensional pore model is established and the influence of pore structure parameters on the mechanical properties of porous sandstone is discussed [14]. Based on the model of micro pore structure of rock, the inner mechanism of the physical and mechanical behavior of rock is analyzed [15]. The relationship between pore structure parameters and physical and mechanical properties of coal and rock mass is studied, and an empirical relationship model is established for evaluation, classification, and prediction of physical parameters [16]. Most scholars study the change of macroscopic properties of the coal and rock mass, by experimental means, which can only indirectly reflect the effect of pore structure on the physical and mechanical properties. But the quantitative analysis of the internal mechanism between the coal pore structure and macroscopic physical and mechanical properties is less. On the basis of statistical analysis, the threedimensional pore reconstruction model provides a new approach for the quantitative study of the mechanism of pore structure parameters of coal rock.
Based on the above analysis, selecting the typical low permeability coal seam in the coalfield of South Judger Basin in Xinjiang as the research object, we constructed a threedimensional cylinder pore model of low permeability. We analyze the influence of the distribution law of micro pores on the physical and mechanical properties of the low permeability coal rock using the numerical simulation method.
2. Pore Structure Parameter Test
The experiment was made with the low permeability coal seam in the coalfield of South Judger Basin in Xinjiang, and the pore structure characteristics of low permeability coal seam were obtained by scanning electron microscope. Using Otsu threshold banalization processing method [17], SEM image of coal samples under different magnification is processed (Figure 1) and pore structure parameters, such as pore diameter, number, location, spacing, and statistical properties, are extracted. A large number of CT and SEM experiments show that the number of pores satisfies an approximate uniform distribution along the circumference, and the pore space and pore size are in normal distribution and exponential distribution [18–20].
(a) Coal sample SEM diagram
(b) Banalization processing graph
According to the SEM data of coal samples, the pore statistical model with different porosity is established. The pore distribution in pore statistical model is uniform and the pore size distribution is exponential distribution. The pore size distribution function is expressed aswhere , , are constants and is the aperture. The diameter of the image is calculated, and the ratio of the number of apertures to the total aperture is calculated (Figure 2).
The circle is divided into 25 equal intervals. If the number of holes in each interval is in uniform distribution, the following should be satisfied:
The number of pores in the processed image is calculated, and the numerical value of the ratio of the number of pores distributed in different circumferential intervals to the total number of pores is obtained. The specific distribution of pore number is shown in Figure 3.
For a given specimen height section of the pore, according to the plane polar coordinates method to determine the location of its pores, the statistical function of pore coordinate is expressed as
The partial coordinate information of pores on the circular surface in the = 50 mm direction is shown in Table 1.

3. ThreeDimensional Model of Pore Structure
According to the statistics derived from the pore position, number, size parameters, and the distribution function, using the Monte Carlo method, simulated annealing algorithm, and Matlab program, two groups of number series consisting of the parameters and distribution function are generated [21, 22]. Each row of the 3column random sequences in accordance with the characteristics of uniform distribution represents the threedimensional coordinates of each pore, where two groups of onerow and multiplecolumns random number sequence, accordingly, represent the random distribution of pore diameter. At the same time, the number of pores can be determined according to the measured porosity and pore size distribution, and the distribution of the number of pores can meet the regularity of the uniform distribution.
Because the space position of the pore, pore size distribution, and the number of changes are random, the difficulty of the finite element program exists in the grid, element size, and quantity control, so FLAC3D is chosen to construct the pore model. Considering the limitation of minimum size of the unit of FLAC program, scaleup of pore radius obtained in the experiment was made, at the same time keeping the porosity and porosity distribution unchanged and thus making the model and the actual coal samples have good geometric similarity.
Importing the generated random number series into the FLAC3D through the selfcompiled program and, then, using FISH language, a threedimensional cylinder is built by pore model with a diameter 50 mm, height 100 mm, and maximum pore radius 3.6 mm (Figure 4).
4. Numerical Simulation of Uniaxial Compression of ThreeDimensional Pore Reconstructed Models
The numerical simulation of the uniaxial compression of the 3D reconstruction model with different pores was carried out in FLAC3D. In the experiment, the upper and lower boundaries are simultaneously moved at a rate of 5 × 10^{−7} m/s (Figure 5), and the FISH language is used to monitor the stress and strain parameters in the model. Take MohrCoulomb model as the material constitutive model, and the physical and mechanical parameters of the material are given in Table 2.

4.1. Effect of Porosity on Compressive Strength
For the porosity of 0%, 1.86%, 3.7%, 7.4%, and 11.2%, a total of models were subjected to uniaxial numerical tests. The stressstrain curve of model with porosity 0% is shown in Figure 6. According to the stressstrain curves of different models, the variation of compressive strength with porosity is obtained (Figure 7).
From test result analysis, with the increase of porosity, the compressive strength of different pore model decreases and porosity and compressive strength shows negative exponent relationship. The model with porosity of 11.2% has the minimum compressive strength of 31 Mpa, which is 8.58% lower than the maximum compressive strength 33.91 Mpa of the model with porosity of 0%. Because the pore size distributions of different porosity models in the experiment have the same exponential distribution, the number of pores and the distribution of pores have the greatest influence on the compressive strength of the model. With the increase of porosity in the model, the damage of the model increases and the compressive strength decreases.
The experimental results show that there is a nonlinear relationship between porosity and compressive strength.
4.2. Mechanical Properties of Different Porosity Models
Stress slice of the cylinder model has a porosity of 0%~11.2% at its height of 50 mm (Figure 8).
The experimental results show that the stress distribution of different porosity models is obviously different. The stress distribution of the porosity model with porosity 0% is symmetrical, and the stress distribution becomes more inhomogeneous with the increase of porosity.
The plastic zone of different models at the peak value of stressstrain curve shows that the failure zone is in the plastic deformation stage (Figure 9). According to the stressstrain curve of the rock under uniaxial compression, the cylinder is in the stage of unstable fracture development, and the rupture of this stage is developing until the specimen is completely destroyed. The results show that the failure zone is in the plastic deformation stage.
(a) 0%
(b) 1.86%
(c) 3.7%
(d) 7.4%
(e) 11.2%
In this paper, the uniaxial compression of the cylinder model with different porosities was carried out in the same loading rate and time step, and the plastic zone of the failure process of each model was obtained (Figure 10).
(a) 0%
(b) 1.86%
(c) 3.7%
(d) 7.4%
(e) 11.2%
The results of numerical simulation are analyzed as follows:(1)At the porosity of 0% damage shows symmetric form. In the initial state of model loading symmetrical shear crack appears at the upper and lower ends of the model. With the increase of loading steps, shaped crack extends gradually to the middle position of the specimen. At the peak of loading the cylinder rupture and the cracks at the upper and lower ends convergent in the middle, cylindrical damage is mainly shear failure.(2)The destruction of the model with porosity 1.86%~3.7% basically shows symmetric form; at the initial loading stage, crack shear failure model is irregular. With the increase of the number of steps of loading, crack failure gradually extends to the middle position; the crack failure is shape, influenced by some of the cracks in the pores of slight deformation; when the model is destroyed, cracks at both ends converge to the middle and extend to both sides, and the failure mode is mainly shear failure.(3)The destruction of the model with porosity 7.4%~11.2% shows asymmetrical form; at the initial loading complex failure cracks appear. With the increase of loading steps, destruction crack expands rapidly; crack of the right and left sides gradually merged and extended to the middle position; at the upper end and the lower end “” shaped failure formed. When the model is damaged, the cracks at the two ends converge in the middle and extend to both sides, with damage in the form of shear failure.(4)For different models, the shaped conjugate slant shear damage is the most common; that is, the angle between the normal vector of damage and the axis of specimen surface is 60 degrees.
Based on the analysis of shear strain increment cloud diagram of different porosity models, the expansion trend of the failure surface is obtained (Figure 11). The shear strain area of different porosity model meets the breaking angle law of 45° + . The shear strain of the model with 0% porosity is symmetrical, with increasing porosity of shear strain area having the characteristics of uneven distribution; the main reason is the uneven stress distribution caused by the different pore diameter and distribution in the cylinder.(1)When the porosity is smaller, the distribution of the pore has less cumulative effect.(2)When the porosity is larger, the increase of the pore number and the pore size have a significantly increasing effect on the cumulative effect of the stress.
5. Numerical Simulation of Seepage Flow in Pore Reconstruction Model
Seepage numerical simulation was taken on a 3D pore reconstruction model with porosity of 7.4%. In the loading process, certain axial pressure , lateral pressure , and pore pressures , were added initially; then displacement control loading was used in the whole loading process, with loading displacement increment = 0.01 mm, = 4 MPa, = 2.3 MPa, and = 3.8 MPa.(1)The evolution law of the coal sample loading stepvertical loadpermeability coefficient is described; the consistency between the change of permeability coefficient and the crack growth is verified (Figure 12). At the same time, the initiation and propagation of crack play a very important role in the selection of flow paths.(2)Compared with the pore waterfree pressure test, the pore water pressure test specimen is easier to produce micro crack and to be destroyed, and the failure surface is more complex, the peak value of vertical load is also reduced by 3%, and the number of steps is also two steps ahead (Figure 13).(3)When the models have the same porosity, their permeability variations changed considerably, which indicates that the pore structure characteristics have important influence on the permeability of the model.
6. Conclusions
(1)Based on SEM test, we have obtained the pore structure parameters of low permeability coal seam by MCMC and digital image processing methods. And on this basis, threedimensional cylindrical pore models with different porosity (1.86%, 3.7%, 7.4%, and 11.2%) have been constructed.(2)Uniaxial compression tests of different pore models reveal the influence of pore structure parameters on the physical and mechanical properties of coal and rock mass:①porosity and compressive strength curve have negatively exponential relationship, and the quantity and distribution of pore have the greatest influence on the compressive strength of the model;②under different porosity, the stress distribution is different. With the increase of porosity, the stress distribution gradually becomes nonuniform;③the compressive failures of different pore models are mainly shear failures, and the shape of fracture section is related to porosity;④the pore structure parameters have different influence on the stress concentration of the model. When the porosity is smaller, the distribution of the pore has less cumulative effect; when the porosity is larger, the increase of the pore number and the pore size have an obvious increasing influence on the cumulative effect of the stress.(3)The evolution law of the coal sample loading stepvertical loadpermeability coefficient is described. The consistency of the crack growth and the variation of the permeability coefficient are verified. At the same time, the influential mechanism of the pore structure parameters on the permeability coefficient of coal and rock is revealed.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work is supported by major scientific and technological projects of “13115” Science and Technology Innovation Project in Shaanxi (2009ZDKG57), science and technology to support major projects in Xinjiang (2014AB025), and major projects of natural science and education in Shaanxi (2015JS060).
References
 G. Z. Deng, S. B. Wang, and B. X. Huang, “Research on behavior character of crack development induced by hydraulic fracturing in coalrockmass,” Chinese Journal of Rock Mechanics and Engineering, vol. 23, no. 20, pp. 3489–3493, 2004 (Chinese). View at: Google Scholar
 W. J. Li, J. J. Luo, X. X. Liang et al., “Rule of porous structure of semicoke in the underground coal gasification,” Coal Conversion, vol. 32, no. 2, pp. 10–13, 2009. View at: Google Scholar
 F. J. Dong, W. Sun, W. W. Chen et al., “Influence of micropore structure on water flooding effectiveness in the low permeability sandstone reservoir: a case of San Jianfang reservoir,” Journal of Northwest University (Natural Science Edition), vol. 40, no. 6, pp. 1041–1045, 2010 (Chinese). View at: Google Scholar
 J.J. Fan, Y.W. Ju, S.B. Liu, and X.S. Li, “Micropore structure of coals under different reservoir conditions and its implication for coalbed methane development,” Journal of the China Coal Society, vol. 38, no. 3, pp. 441–447, 2013 (Chinese). View at: Google Scholar
 G. Z. Deng, “Study on seepage features of coal bed cracks under controlling stress,” Journal of China Coal Society, vol. 25, no. 6, pp. 593–598, 2000 (Chinese). View at: Google Scholar
 X. D. Zhang, Y. Qin, and S. X. Sang, “Status quo and prospect of studies on CBM reservoir adsorptive characteristics,” Coal Geology of China, vol. 17, no. 01, pp. 16–21 + 29, 2005 (Chinese). View at: Google Scholar
 X. M. Meng, Study on the Pore Structure of Coals and Characteristics of Gases Adsorption on Coals, Shandong University of Science and Technology, 2007.
 G. Z. Deng and W. S. Zhu, “An experiment research on the crack propagation in rock mass,” Journal of Experimental Mechanics, vol. 17, no. 2, pp. 177–183, 2002 (Chinese). View at: Google Scholar
 Y.G. Zhang, Z.M. Zhang, and Y.X. Cao, “Deformedcoal structure and control to coalgas outburst,” Journal of the China Coal Society, vol. 32, no. 3, pp. 281–284, 2007 (Chinese). View at: Google Scholar
 L. Vernik, M. Bruno, and C. Bovberg, “Empirical relations between compressive strength and porosity of siliciclastic rocks,” International Journal of Rock Mechanics and Mining Sciences and, vol. 30, no. 7, pp. 677–680, 1993. View at: Publisher Site  Google Scholar
 Y. H. Hatzor and V. Palchik, “A microstructurebased failure criterion for Aminadav dolomites,” International Journal of Rock Mechanics and Mining Sciences, vol. 35, no. 6, pp. 797–805, 1998. View at: Publisher Site  Google Scholar
 B.S. Wang, Y. Chen, H.K. Ge, L.L. Song, and T.F. Wong, “A discrete element model of porous rock deformation,” Chinese Journal of Geophysics, vol. 48, no. 6, pp. 1336–1342, 2005 (Chinese). View at: Publisher Site  Google Scholar
 Y. Ju, Y. M. Yang, Z. D. Song et al., “A statistical model for porous structure of rocks,” Science in China (Ser E), vol. 38, no. 7, pp. 1026–1041, 2008 (Chinese). View at: Google Scholar
 Y. M. Yang, Y. Ju, H. B. Liu, and H. J. Wang, “Influence of porous structure properties on mechanical performances of rock,” Chinese Journal of Rock Mechanics and Engineering, vol. 28, no. 10, pp. 2031–2038, 2009. View at: Google Scholar
 Y. Ju, J. T. Zheng, M. Epstein, L. Sudak, J. Wang, and X. Zhao, “3D numerical reconstruction of wellconnected porous structure of rock using fractal algorithms,” Computer Methods in Applied Mechanics and Engineering, vol. 279, pp. 212–226, 2014. View at: Publisher Site  Google Scholar
 B. M. Yu, “Advances of fractal analysis of transport properties for porous media,” Advances in Mechanics, vol. 33, no. 3, pp. 333–346, 2003. View at: Google Scholar
 N. Otsu, “A threshold selection method from GrayLevel histograms,” Automatica, vol. 11, no. 285296, pp. 23–27, 1975. View at: Google Scholar
 G. Desbois, J. L. Urai, P. A. Kukla, J. Konstanty, and C. Baerle, “Highresolution 3D fabric and porosity model in a tight gas sandstone reservoir: a new approach to investigate microstructures from mm to nmscale combining argon beam crosssectioning and SEM imaging,” Journal of Petroleum Science and Engineering, vol. 78, no. 2, pp. 243–257, 2011. View at: Publisher Site  Google Scholar
 L. M. Keller, L. Holzer, R. Wepf, and P. Gasser, “3D geometry and topology of pore pathways in Opalinus clay: Implications for mass transport,” Applied Clay Science, vol. 52, no. 12, pp. 85–95, 2011. View at: Publisher Site  Google Scholar
 L. M. Keller, P. Schuetz, R. Erni et al., “Characterization of multiscale microstructural features in opalinus clay,” Microporous and Mesoporous Materials, vol. 170, pp. 83–94, 2013. View at: Publisher Site  Google Scholar
 C. L. Yeong and S. Torquato, “Reconstructing random media,” Physical Review E. Statistical, Nonlinear, and Soft Matter Physics, vol. 57, no. 1, pp. 495–506, 1998. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Jiao, F. H. Stillinger, and S. Torquato, “A superior descriptor of random textures and its predictive capacity,” Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 42, pp. 17634–17639, 2009. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 Guangzhe Deng and Rui Zheng. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.