Experimental Development of Coal-Like Material with Solid-Gas Coupling for Quantitative Simulation Tests of Coal and Gas Outburst Occurred in Soft Coal Seams

Solid-gas coupling coal-like materials are essential for simulating coal and gas outbursts and the long-term safety study of CO2 sequestration in coal. However, reported materials still differ substantially from natural coal in mechanical, deformation and gaseous properties; the latter two aspects are common not considered. There is a lack of a definite and quantitative preparation method of coal-like materials with high similarity for future reference. Here, 25 groups of raw material ratios were designed in the orthogonal experiment using uniaxial compression, shearing and adsorption/desorption tests. Experiment results indicated that the coal-like materials were highly similar to soft coals in properties mentioned above. And range analysis revealed the key influencing factors of each mechanical index. The gypsum/petrolatum ratio controls the density, compressive strength, elastic modulus, cohesion and deformation characteristic. The coarse/fine coal powder (1–2 and 0–0.5 mm) controls the internal friction angle and is the secondary controlling factor for compressive strength and elastic modulus. The effect of coal particle size on the sample strength was studied using scanning electron microscope (SEM). When the gypsum/petrolatum ratio increased, the deformation characteristics changed from ductile to brittle. The different failure modes in the samples were revealed. The coal powder content is a key in the gas adsorption/desorption properties and an empirical formula for estimating the adsorption capacity was established. Based on the range analysis of experimental results, a multiple linear regression model of the mechanical parameters and their key influencing factors was obtained. Finally, a composition closely resembling the natural coal was determined, which differs by only 0.47–7.41% in all parameters except porosity (11.76%). Possible improvements and extension to similar materials are discussed. The findings of this study can help for better understanding of coal and gas outburst mechanism and stability of CO2 sequestration in soft coal seams.


Introduction
Coal and gas outburst accidents are among the most serious disasters affecting coal mining [1][2][3][4] and a problem in almost all major coal-producing countries [2][3][4].Coal and gas outburst is an extremely complex gas dynamic phenomenon, in which large amounts of coal and gas are often ejected at a very fast rate from the coal rock seam to the mining space in a very short period of time (a few seconds to a few minutes).Such outbursts could destroy underground facilities, damage the ventilation system, cause a large number of casualties and even induce secondary accidents such as gas burning or explosion [2,[4][5][6][7].Today, more than half of the disasters due to coal and gas outbursts occur in China [8], causing major economic losses and casualties.
Many factors could cause coal and gas outbursts, including in-situ stress, gas pressure, geological structure, physical and mechanical properties of coal, mining methods and so forth.[6][7][8][9][10].At present, there is no coherent model that can fully reveal their internal mechanism.In addition, the sudden, transient and dangerous nature makes it almost impossible to observe or study these processes on site.Therefore, simulation tests have been rapidly developed and used to explore the mechanism and process of coal and gas outbursts [5,9,[11][12][13][14][15].However, coal seams that are prone to coal and gas outbursts often have high gas content and low strength and so it is difficult to retrieve samples from the field for large-scale simulation tests.Hence, the development of similar solid-gas coupling coal-like materials in the lab is indispensable for studying the mechanism of coal and gas outbursts.
In addition, there is growing realization in recent years among researchers that coal seams with strong gas adsorption/desorption capacity may be able to capture and store CO 2 [16][17][18][19][20][21].However, the key challenge here is the uncertain impact of CO 2 on the mechanical properties of coal and how this will affect the long-term safety and stability of storage.On the other hand, natural coal has extremely complex composition and physical structure, as well as strong heterogeneity-even two samples close to each other in the same coal seam could have significantly different mechanical properties.This seriously hinders the analysis of laboratory test results and the understanding of the above uncertainty.So, homogeneous and reproducible coal-like materials (also called reconstituted coal sample) will provide significant advantages for studies related to CO 2 sequestration in coal.
The key to solve the above two problems is to develop solid-gas coupling coal-like materials that closely resemble the target natural coal in their mechanical and gaseous properties and deformation features.Table 1 lists representative papers on preparing coal-like materials.There are mainly three types of methods.The first is to directly press the coal powder without any added substances [5,12,22].For example, in 1953, a one-dimensional coal and gas outburst simulation test was completed using briquettes that were cold pressure formed [22].Skoczylas [12] used fine coal powder to produce a series of briquettes (40 mm in diameter, 110 mm in height and with porosities of 11.2-32.0%)and the test results indicated that they had similar gaseous characteristics to natural coal.However, the coal samples prepared by only molding pressure had very low strength.The second type is pressing coal powders with water, oil or diesel but without any binder [11,20,[23][24][25][26][27].Jasinge et al. [20] prepared reconstituted brown coal sample by compacting 0-1 mm coal particles, in order to study the effect of coal swelling on its permeability in the laboratory.While there was a marked similarity in permeability between natural brown coal and reconstituted specimens, the strength of this type of briquette sample was still lower than natural coal and the gas adsorption/desorption characteristics were also quite different.The third type is pressing a powder mixture consisted of coal powder, cement, sand, lime and other substances [6,10,13,18,19,21,[28][29][30].Hu et al. [13] made coal-like materials using coal powder, cement, water, sand and activated carbon.They studied the effect of the proportion of each component on the density and mechanical properties of the samples.In 2017 and 2018, Wang et al. [10] and Zhao et al. [6] studied the gas and CO 2 adsorption characteristics of such coal-like materials.Zhang et al. [30] analyzed the influence of different material ratios on the elastic modulus.This last type of coal sample and its preparation method represent outstanding progress in the development of coal-like materials.Nevertheless, the deformation characteristics and failure modes of similar materials, which play an important role in the preparation, forming and developing stages during coal and gas outburst [13], are often not considered and rarely compared with those of natural coal.A clear discussion of related researches can be seen in Figure 1.
The abovementioned studies made important contributions to laboratory research on coal and gas outbursts and CO 2 sequestration.However, at present, there are still the following problems in the development of coal-like materials: (1) The coal powders used in existing literature have different particle size and molding pressure, molding time and additives are also different.In future simulation test, repeated adjustments are usually needed to determine a suitable proportion of materials.There is still a lack of definite and quantitative preparation method of coal-like materials for future reference.
(2) The similarity remains low between coal-like materials and natural coal in mechanical and physicochemical properties, especially in that the former still have lower strength and very few indexes were measured.Their deformation features and adsorption and desorption characteristics are often not considered and still significantly different from those of natural coal [6,10,13,19,29].As a result, most simulation studies only reported qualitative rather than quantitative results [13,31].Today, the understanding of coal and gas outburst mechanism remained poor, the predictions were inaccurate and coal and gas outburst accidents still account for a large proportion of accidents in coal mines [14,32,33].
In this study, we propose a systematic and quantitative development method to determine a solid-gas coupling coal-like material that could closely resemble given natural coal in the mechanical properties, adsorption/desorption features and deformation characteristics.First, according to the physical and mechanical parameters of natural coal samples, raw materials were selected, that is, coarse and fine coal powder, gypsum, petrolatum and light calcium carbonate.Second, the orthogonal method was used to design 25 groups of schemes with different raw material ratios.Meanwhile, the sample preparation method was optimized to improve the homogeneity and reproducibility.Then, density, uniaxial compressive strength (UCS), elastic modulus (EM), cohesion, internal friction angle (IFA), deformation and failure modes, porosity and adsorption/desorption characteristics of samples were tested.Third, range analysis revealed the key influencing factors of each mechanical and physicochemical index and the levels of influence of each factor on a given index.A multivariate linear model was proposed to predict the mechanical parameters of the materials.Fourth, the effects of different particle sizes on the sample strength were examined by scanning electron microscope and the effects of gypsum/petrolatum ratio on the deformation evolution of samples were also analyzed.Finally, a material ratio scheme closely resembling the natural coal was identified.Related improvements and remaining problems in the preparation of similar material were also discussed.This study would contribute to the prevention and control of coal and gas outbursts, as well as other physical simulation tests related to soft coal seams.

Preparation and Testing of Natural Coal Samples
Studies have shown that most of the coal and gas outburst occur in tectonically deformed coal (TDC) [34][35][36].This is mainly because TDC is formed under mono-or multi-phase strong tectonic movements.Its original structure has been deformed or damaged, resulting in larger adsorption surface area and higher gas content than the primary coal [37].In this paper, the natural coal samples

Preparation and Testing of Natural Coal Samples
Studies have shown that most of the coal and gas outburst occur in tectonically deformed coal (TDC) [34][35][36].This is mainly because TDC is formed under mono-or multi-phase strong tectonic movements.Its original structure has been deformed or damaged, resulting in larger adsorption surface area and higher gas content than the primary coal [37].In this paper, the natural coal samples were selected from a typical tectonically deformed coal seam from No. B-1 Coalbed, Zhengzhou, Western Henan, China, where there have been frequent coal and gas outburst disasters over the last 60 years [7].According to the new structural-genetic classification system [38], the TDC samples belong to medium-rank and mylonitic structure coal and are characterized by low strength and strong gas adsorption/desorption [32].
Because this coal seam is soft and fragile, it is difficult to obtain large samples using the traditional drill core sampling method.The natural coal standard samples were obtained successfully only after many attempts and improvements.Specifically, a square iron sampler was used to cut the underground coal seam to obtain a sample of approximately 100 mm × 100 mm × 150 mm and the sample was quickly placed in a sealed bag packaging for preservation.The sample was then carefully polished in the lab using a grinding machine.Samples of two target sizes were prepared: 50 mm × 50 mm × 100 mm for uniaxial compression test and 50 mm × 50 mm × 50 mm for variable angle compression shear test.The sample standard requirements and test procedures are in accordance with methods for determining the physical and mechanical properties of coal and rock (Part 7 [39] and Part 11 [40]).The sample loading rate was 0.5 mm/min.The basic physical and mechanical parameters of coal sample are as follows: apparent density =1.280 g/cm 3 , UCS = 1.72 MPa, elastic modulus = 126.35MPa, cohesion = 0.17 MPa and internal friction angle = 27 • .However, the success rate of natural coal sample preparation remained only approximately 5%.The proximate analysis results of this TDC are as follows: moisture content Mad = 0.95%, ash content Ad = 10.42%,volatile matter content Vad = 13.95%,porosity = 15.3%, gas adsorption constants a = 39.789 and b = 1.113 and index of initial velocity of gas diffusion ∆p = 26-29.
The stress-strain curve in Figure 2 contains the typical stages of initial compaction, linear elastic deformation, strength hardening and softening.Notably, it includes a longer compaction phase than primary coal seam, as well as a plastic or shearing deformation stage.The peak-failure strain of TDC reached 1.73%, which was significantly larger than those of brittle coal or rock (approximately 0.5%) [41].In addition, microstructure test (Figure 3) showed that the coal body was very broken.The original structure was almost completely destroyed and was composed of particles of different sizes and the texture was loose.This was the reason why the density and strength of TDC were both usually smaller than those of the primary coal seam.
Processes 2019, 7, x FOR PEER REVIEW 6 of 25 with methods for determining the physical and mechanical properties of coal and rock (Part 7 [39] and Part 11 [40]).The sample loading rate was 0.5 mm/min.The basic physical and mechanical parameters of coal sample are as follows: apparent density =1.280 g/cm 3 , UCS = 1.72 MPa, elastic modulus = 126.35MPa, cohesion = 0.17 MPa and internal friction angle = 27°.However, the success rate of natural coal sample preparation remained only approximately 5%.The proximate analysis results of this TDC are as follows: moisture content Mad = 0.95%, ash content Ad = 10.42%,volatile matter content Vad = 13.95%,porosity = 15.3%, gas adsorption constants a = 39.789 and b = 1.113 and index of initial velocity of gas diffusion Δp = 26-29.
The stress-strain curve in Figure 2 contains the typical stages of initial compaction, linear elastic deformation, strength hardening and softening.Notably, it includes a longer compaction phase than primary coal seam, as well as a plastic or shearing deformation stage.The peak-failure strain of TDC reached 1.73%, which was significantly larger than those of brittle coal or rock (approximately 0.5%) [41].In addition, microstructure test (Figure 3) showed that the coal body was very broken.The original structure was almost completely destroyed and was composed of particles of different sizes and the texture was loose.This was the reason why the density and strength of TDC were both usually smaller than those of the primary coal seam.

Similarity Criteria and Index
Coal and gas outburst is a complex solid-gas coupling process.The currently accepted view is that it is mainly caused by the combined action of in-situ stress, gas pressure and physical and mechanical properties of coal [7][8][9][10].The coal and gas outburst consists of three successive stages [13], among which the static deformation and failure of coal occur during the important outburst preparation stage, while the fracturing of gaseous coal and the movement of pulverized coal and gas

Similarity Criteria and Index
Coal and gas outburst is a complex solid-gas coupling process.The currently accepted view is that it is mainly caused by the combined action of in-situ stress, gas pressure and physical and mechanical properties of coal [7][8][9][10].The coal and gas outburst consists of three successive stages [13], among which the static deformation and failure of coal occur during the important outburst preparation stage, while the fracturing of gaseous coal and the movement of pulverized coal and gas occur during the forming and developing stage.Therefore, several similarity criteria (C i in Equation (1) below) have been proposed to describe the whole process of coal and gas outburst, by using mechanical, deformation and energy models and considering the experimental conditions [6]: where γ, n, ϕ, p, σ, E, c and l, are the volumetric weight, porosity, internal friction angle, gas pressure, compressive strength, elastic modulus, cohesion and length, respectively.Judging from the above, the similarity index of the coal-like materials should include the physical and mechanical parameters (density, uniaxial compressive strength, elastic modulus, internal friction angle, etc.) and deformation properties and failure mode, as well as porosity, absorption constants (a, b) and desorption index (∆p) of the initial gas diffusion velocity.

Composition Selection
Coal-like materials generally consist of aggregates, binders and additives.The raw materials should be selected according to the following principles: (1) similar to the natural coal material, (2) abundant and low-cost and (3) safe, non-toxic and environment-friendly.In the coal and gas outburst tests [6], the density of coal-like material should be the same as that of natural coal.So, the natural coal powder was selected as the aggregate.Another important reason for using coal powder was that it had good adsorption characteristics similar to the natural coal.However, the particle sizes of coal powders reported in Table 1 were different.Some studies [30,42] have shown that larger coal particles (greater than 3 mm) would undergo secondary crushing during the molding and pressing process, which would have an adverse effect on molding.Therefore, two types of coal powders with particle sizes of 1-2 mm and 0-0.5 mm were used respectively as the coarse and fine aggregates.The fine coal powder was the main aggregate to ensure the molding quality.Second, the binder was selected.Normally, the binder has the greatest influence on the strength and deformation parameters of the material.As the strength of sample obtained by pressing only coal powder was too low, gypsum was selected to adjust the strength parameters of coal-like materials.Portland cement was not used here, because the cement hardening time was too long (more than 28 days).Its strength would also change with time, which could lead to strength instability of coal-like materials [19] and poor homogeneity and reproducibility.On the other hand, petrolatum was selected to adjust the deformation characteristics of the coal-like materials.
The molding pressure has an important effect on the density and porosity of briquette [5,12].Excessive molding pressure would damage aggregates with large particle size and affect the molding quality.Therefore, the molding pressure was designed based on the ground stress level of the coal seam.The vertical stress level of the coal seam was approximately 7.2 MPa.In addition, in order to allow the gypsum to play a bonding role, it was necessary to add an appropriate amount of water.After repeated tests, it was found that when the moisture exceeded 11%, the material was too wet and would affect the molding, while the raw material to water ratio of 10:1 was relatively reasonable.Because the density of gypsum (2.3 g/cm 3 ) is higher than that of natural coal, it was also necessary to use minuteness super white light calcium carbonate with very stable mechanical properties to match the target density while using different material ratios.The selected material compositions are shown in Table 2 and Figure 4.  Ultra-fine, ultra-white, light calcium carbonate, 0.1 μm < particle size ≤ 1 μm, stable performance, density 0.54 g/cm 3  Water Tap water

Design Scheme
The components and their ratios both have an impact on the physical and mechanical properties of coal-like materials.The orthogonal test method was used to design 25 groups of ratios, as shown in Tables 3 and 4. The advantage of this method was that the schemes included a combination of any two factors at all levels and reduced the amount of experimental work [43].

Design Scheme
The components and their ratios both have an impact on the physical and mechanical properties of coal-like materials.The orthogonal test method was used to design 25 groups of ratios, as shown in Tables 3 and 4. The advantage of this method was that the schemes included a combination of any two factors at all levels and reduced the amount of experimental work [43].

Sample Preparation Method
(1) The raw materials were weighed accurately using an electronic scale according to the given ratios based on a total amount of 1000 g. ( 2) The dry components (coal powder, gypsum powder and calcium carbonate powder) were mixed and stirred uniformly.(3) A mixture of water and petrolatum was heated to 45 to 50 • C to melt it into a liquid and then quickly mixed with the materials.(4) The well-mixed coal-like materials were divided into five equal parts by weight and put into the mold.Each layer was pressed by a servo press machine at the speed of 120 mm/min.The target pressure was 14.137 kN (the calculated stress was 7.2 MPa) and the load was maintained for 10 s.The actual load-time curve during the layer-by-layer compaction press is shown in Figure 5a.Finally, the whole sample was compacted by the press at the molding pressure of 7.2 MPa and maintained for 10 min, mainly to make the overall force on the whole sample consistent.The loading path on the whole sample can be seen in Figure 5b.It should be noted that the surface between two adjacent layers should be roughened before filling in the next layer of material to increase the adhesion between layers.The steel mold had an inner diameter of 50 mm, a height of 120 mm and a wall thickness of 10 mm.The mold was approximately 20 mm higher than the sample, making it convenient for the fifth material to be added to the mold at one time.The inner surface of the mold had to be very smooth to reduce the friction.( 5) With a self-designed demolding device, as shown in Figure 5d, the molded sample was pushed out at a speed of 5 mm/min by a press.The sample slowly entered the hollow steel pipe (wall thickness 10 mm, height 200 mm, inner diameter 100 mm, with soft material placed on the bottom to prevent sample damage).Then, the sample was cured at room temperature for 10 to 15 days.The water content was tested continuously in this process.When the water content was the same as that of the natural coal, the sample was sealed and stored in a plastic bag in time.In particular, it should be pointed out that many reconstituted coal samples were destroyed during the extrusion process in the past [18], while by using the designed demolding method here, the success rate of sample preparation was 100%.The use of a press machine with quantitative loads for compacting the materials was also an improvement over the manual pressing used in previous studies, which was difficult to quantify and lowered the homogeneity and reproducibility of samples.
Processes 2019, 7, x FOR PEER REVIEW 10 of 25 extrusion process in the past [18], while by using the designed demolding method here, the success rate of sample preparation was 100%.The use of a press machine with quantitative loads for compacting the materials was also an improvement over the manual pressing used in previous studies, which was difficult to quantify and lowered the homogeneity and reproducibility of samples.

Sample Tests
The sample preparation met the requirements suggested by International Society for Rock Mechanics [44].The height of the sample was more than twice the diameter.For each group of formulation scheme, 6 samples of Φ 50 mm × 102 mm were made (a total of 150 samples).Two samples in each group were used for the uniaxial compression test (Figure 6a), three were cut into 6 pieces Φ 50 mm × 50 mm in size to perform the variable angle shear tests (Figure 6b).The spare sample would be used in the isothermal adsorption and desorption test.
The CSS-44300 universal testing machine (Changchun testing machine research institute, Changchun, China) with electro-hydraulic servo was used to perform the uniaxial compression and variable angle shear tests.The test steps referred to references [39,40].The sample loading rate was 0.5 mm/min and the surface unevenness did not exceed 0.2 mm.

Sample Tests
The sample preparation met the requirements suggested by International Society for Rock Mechanics [44].The height of the sample was more than twice the diameter.For each group of formulation scheme, 6 samples of Φ 50 mm × 102 mm were made (a total of 150 samples).Two samples in each group were used for the uniaxial compression test (Figure 6a), three were cut into 6 pieces Φ 50 mm × 50 mm in size to perform the variable angle shear tests (Figure 6b).The spare sample would be used in the isothermal adsorption and desorption test.
The CSS-44300 universal testing machine (Changchun testing machine research institute, Changchun, China) with electro-hydraulic servo was used to perform the uniaxial compression and variable angle shear tests.The test steps referred to references [39,40].The sample loading rate was 0.5 mm/min and the surface unevenness did not exceed 0.2 mm.

Results of Orthogonal Test
The mechanical property test results are shown in Table 5.The property parameters of the 25 groups of coal-like materials covered a wide range: density of 1.171-1.305g/cm 3 , uniaxial compressive strength of 0.55-2.20 MPa, elastic modulus of 32.06-196.18MPa, cohesion of 0.068-0.248MPa and the internal friction angle of 16.8-35.6°.These samples could be used for physical simulation of most soft coal seams [13,19].

Results of Orthogonal Test
The mechanical property test results are shown in Table 5.The property parameters of the 25 groups of coal-like materials covered a wide range: density of 1.171-1.305g/cm 3 , uniaxial compressive strength of 0.55-2.20 MPa, elastic modulus of 32.06-196.18MPa, cohesion of 0.068-0.248MPa and the internal friction angle of 16.8-35.6• .These samples could be used for physical simulation of most soft coal seams [13,19].In addition, according to the orthogonal experimental theory, range analysis (also called visual analysis method) was used to reveal the influence of various factors on results from the orthogonal experiment [43].That is, the results corresponding to the same level of each factor were averaged and the range was obtained by subtracting the minimum average value from the maximum average value for each level.The range value reflected the influence of different levels of a given factor on the index of interest.A large range indicated that different levels of this factor had strong and different influences on the test results.From the test results in Table 5, the average and range of each level of all factors affecting the density, compressive strength, elastic modulus, cohesion and internal friction angle could be obtained, as shown in Table 6.The effect of different levels of various influencing factors on each index could also be obtained, as shown in Figure 7.
The following observations can be found from Table 6 and Figure 7. (1) For density, the gypsum/petrolatum played the main controlling role.The calcium carbonate content was the secondary factor determining the material density.In Figure 7a, the density of coal-like material increased with increasing gypsum/petrolatum ratio but decreased linearly with increasing calcium carbonate content.(2) For UCS, the gypsum/petrolatum ratio had the highest degree of influence on the compressive strength.The coarse/fine aggregate ratio played a secondary controlling role, followed by the aggregate/binder ratio.Factor D (calcium carbonate content) had almost no influence.In addition, Figure 7b shows that (i) the compressive strength tends to be positively correlated with the gypsum/petrolatum ratio, indicating that increasing the gypsum content could improve the material strength significantly.(ii) The compressive strength first increased then decreased with increasing coarse/fine aggregate ratio, indicating that a reasonable ratio of aggregate particles of different sizes had an important influence on the material strength.When the fine aggregate/coarse aggregate content was approximately 70%, the aggregate had the highest strength.(iii) The compressive strength of coal-like material and aggregate/binder showed a negative correlation.( 3) For E, the influence of each factor on the elastic modulus was similar to that on the compressive strength, in the same order of C > B > A > D. (4) For cohesion, the influence of each factor was in the order of C > A > B > D. In Figure 7d, the cohesion of material increased with the increase of gypsum/petrolatum ratio but decreased with increasing aggregate/binder ratios.(5) Regarding IFA, the influence of each factor on internal friction angle was in the order of B > C > A > D. The coarse/fine aggregate had the most significant influence.As shown in Figure 7e, as the coarse aggregate content increased, the internal friction angle increased linearly.This indicated that the friction and bite force between coarse aggregates were significantly larger than those between fine aggregates.Under the condition of ensuring the molding quality, increasing the coarse aggregate content is beneficial to enhancing the internal friction angle of the material.

Deformation Characteristics and Failure Modes of the Coal-Like Samples
The range analysis indicated that the gypsum/petrolatum ratio is the main factor affecting the strength indexes of the specimens.Since samples in the 25 groups used in orthogonal uniaxial compression tests showed different deformation processes, it was necessary to study the effect of gypsum/petrolatum ratio on the deformation features.To this end, coal-like materials with aggregate/binder = 6:1, coarse/fine aggregate = 3:7 and gypsum/petrolatum = 5:5, 6:4, 7:3, 8:2 and 9:1 were prepared and used in deformation tests.Figure 8 shows the stress-strain curves of the five types specimens with different levels of gypsum/petrolatum ratios in the uniaxial compression test.All samples contained the typical stages of initial compaction, linear elastic deformation, strength hardening and softening and a residual phase.Samples with low gypsum/petrolatum ratio had lower strength, longer plastic or ductile deformation phase before peak stress and a gentler drop after the peak value.While the sample with gypsum/petrolatum = 9:1 had the greatest strength and shortest plastic deformation phase, the stress declined linearly after the peak to reach a very small residual value, which corresponds to a brittle characteristic.
More importantly, the damage modes of the specimens with different gypsum/petrolatum ratios under uniaxial compression test could be divided into five main types, named as A, B, C, D and E. Figure 9 demonstrates the physical damage modes and corresponding schematic diagrams of the specimens.Type A is squeezing failure, which belongs to ductile damage and occurred in the sample with gypsum/petrolatum = 5:5, that is, a high content of petrolatum.The major features include apparent lateral convex deformation and the lack of obvious large fracture surfaces, as shown in Figure 9(A-1,A-2).Type B is wedge failure, which mainly occurred in the sample with gypsum/petrolatum = 6:4.As shown in Figure 9(B-1,B-2), large wedge-shaped coal blocks slipped from the sample, so that large fissure surfaces could be found in the specimens.Type C is plastic shear failure and mainly occurred in the sample with gypsum/petrolatum = 7:3.This type is characterized by shear cracks and damage everywhere, the absence of obvious large fissure surface, the coal body collapsing in the form of scales and noticeable dilative shear deformation (Figure 9(C-1,C-2)).Type D is brittle shear failure and occurred in the sample with gypsum/petrolatum = 8:2.In this case, first a major fissure appeared and then it propagated and extended through the whole sample.When the sample broke down, there was a large final fracture surface about 45-60 • from the specimen axial direction (Figure 9(D-1,D-2)).Type E is brittle splitting failure and occurred in the sample with gypsum/petrolatum = 9:1, that is, high gypsum content.Such samples always split quickly after the peak stress.Figure 9(E-1,E-2) present the ultimate destruction characteristic of type E: one or several vertical cracks ran through the top and bottom of the sample and their directions were approximately parallel to the axis of the specimen.
In addition, the occurrence of coal and gas outbursts must be affected by the deformation and failure mechanism of the coal [13,37].According to the different deformation mechanisms, different types of TDC are divided into three series of deformation and ten classes, that is, the brittle, the ductile and the brittle-ductile deformation [38].In this study, upon decreasing the gypsum content or increasing the petroleum content, the deformation of the coal-like materials transitioned from ductile to brittle.The failure modes included typical squeezing damage, wedge splitting, plastic shear failure, brittle shear and brittle fracturing.These varied deformation features are useful in different coal and gas outburst simulation tests.A

Physicochemical Properties
Samples in the four groups of 4, 8, 12 and 17 had similar mechanical parameters (compressive strength, elastic modulus, cohesion, internal friction angle, etc.) and deformation features to the natural coal.Thus, they were selected for further testing of porosity and gas adsorption/desorption.A

Physicochemical Properties
Samples in the four groups of 4, 8, 12 and 17 had similar mechanical parameters (compressive strength, elastic modulus, cohesion, internal friction angle, etc.) and deformation features to the natural coal.Thus, they were selected for further testing of porosity and gas adsorption/desorption.

Physicochemical Properties
Samples in the four groups of 4, 8, 12 and 17 had similar mechanical parameters (compressive strength, elastic modulus, cohesion, internal friction angle, etc.) and deformation features to the natural coal.Thus, they were selected for further testing of porosity and gas adsorption/desorption.

Porosity
Porosity refers to the proportion of the volume of pores to the total volume.It can be calculated from the apparent density ρ a and true density ρ t of material according to the following equation: After the sample was dried, its apparent density could be obtained by weighing using an electronic scale and size measurement using a Vernier caliper.The true density was measured by the automatic industrial analyzer.From the porosity of the three groups of samples (Table 7), it can be seen that their density and porosity were highly consistent with those of the natural coal.So, these coal-like materials could satisfy the porosity and volumetric weight similarity ratio shown in Equation ( 1) [6]:

Adsorption and Desorption Parameters
The gas adsorption constants a and b, used to quantify the adsorption capacity of coal, were determined by the isotherm adsorption test results using methane of 99.9% purity at 30 • C in combination with the Langmuir equation [6].The parameter ∆p represents the initial velocity of coal gas diffusion [45].In this paper, the TDC and the samples 4, 8, 12 and 17 were crushed and grains with the particle size of 0.2 to 0.25 mm were sifted out.For each sample group, 6 portions of 3.5 g each were weighed.Each portion was placed in the WFC-2 initial gas diffusion velocity analyzer (Zhengzhou huazhi electronic technology co.LTD, Zhengzhou, China).Air tightness check and vacuum treatment were performed on the test instrument.The gas adsorption and desorption of samples were tested under the conditions of 99.9% pure methane, 0.1 MPa and 20 • C. The averages of the test results are shown in Table 8.It can be seen that the gas adsorption and desorption properties of coal-like materials were consistent with those of the natural coal and could meet the test requirements [6].Coal has a strong gas adsorption capacity mainly because of its good porosity and relatively large specific surface area [36,37].Therefore, under the same molding pressure conditions, the coal mass ratio had the most impact on these characteristics.Tables 7 and 8 also verify this feature: the higher the coal content, the larger the porosity of the coal-like materials and the stronger the gas adsorption and desorption capacities.On the other hand, higher contents of gypsum, petrolatum and calcium carbonate reduce the porosity and the methane adsorption/desorption, which is consistent with existing research results [13,46].Besides, previous results [47] presented that coal particles below 3.35 mm in size have almost no influence on the gas adsorption capacity.Therefore, the total amount of fine and coarse coal powder can be used to quantitatively estimate the adsorption constants a, as presented in Equation ( 3) and Figure 10 (with R 2 = 0.9967).
Y a = 0.4841X − 4.25 (3) where Y a is the value of gas adsorption constant a and X is the content of coal powder.
Processes 2019, 7, x FOR PEER REVIEW 18 of 25 of fine and coarse coal powder can be used to quantitatively estimate the adsorption constants a, as presented in Equation ( 3) and Figure 10 (with R 2 = 0.9967).
0.4841 4.25 where Ya is the value of gas adsorption constant a and X is the content of coal powder.Combining Tables 4 and 7, it can be seen that the four coal-like material groups of 4, 8, 12 and 17 had a coal powder ratio of 81.43 to 87.11%, which almost coincides the overall range for groups 1 to 25 (80.57to 89.09%).This indicated that the samples in group 1 to 25 all had similar porosity and gas adsorption/desorption characteristics as the natural coal.

Multivariate Linear Model for Predicting the Mechanical Parameters of Coal-Like Material
According to our range analysis (in Section 3.1.2),each mechanical index has multiple influencing factors.To determine the suitable proportion of materials used in physical simulation, repeated adjustments are usually needed.Thus, results from previous studies can hardly be used directly in future research and it is difficult to promote the rapid development of similar materials.Therefore, it is important to establish a multivariate linear model between the mechanical property indexes and their main influencing factors (as shown in Table 6 and Figure 7).The model can be described by the equations in (4): , , , , where Ci, ai, bi, ci, di and ei (I = 1, 2, 3, 4, 5) are the fitting constants.The contents of coarse aggregate, fine aggregate, gypsum, petrolatum and calcium carbonate were set as X1(%), X2(%), X3(%), X4(%) and X5(%); and the indexes of compressive strength, elastic modulus, cohesion, internal friction angle and density were c Y σ , YE, YC, YF and YD, respectively.
It should be noted that in Figure 11, the elastic modulus of the material increased exponentially (with an excellent correlation of R 2 = 0.9144) against the compressive strength, which is consistent Combining Tables 4 and 7, it can be seen that the four coal-like material groups of 4, 8, 12 and 17 had a coal powder ratio of 81.43 to 87.11%, which almost coincides the overall range for groups 1 to 25 (80.57to 89.09%).This indicated that the samples in group 1 to 25 all had similar porosity and gas adsorption/desorption characteristics as the natural coal.

Multivariate Linear Model for Predicting the Mechanical Parameters of Coal-Like Material
According to our range analysis (in Section 3.1.2),each mechanical index has multiple influencing factors.To determine the suitable proportion of materials used in physical simulation, repeated adjustments are usually needed.Thus, results from previous studies can hardly be used directly in future research and it is difficult to promote the rapid development of similar materials.Therefore, it is important to establish a multivariate linear model between the mechanical property indexes and their main influencing factors (as shown in Table 6 and Figure 7).The model can be described by the equations in (4): where C i , a i , b i , c i , d i and e i (I = 1, 2, 3, 4, 5) are the fitting constants.The contents of coarse aggregate, fine aggregate, gypsum, petrolatum and calcium carbonate were set as X 1 (%), X 2 (%), X 3 (%), X 4 (%) and X 5 (%); and the indexes of compressive strength, elastic modulus, cohesion, internal friction angle and density were It should be noted that in Figure 11, the elastic modulus of the material increased exponentially (with an excellent correlation of R 2 = 0.9144) against the compressive strength, which is consistent with previous studies [13].More importantly, their similarity ratios are the same in the similarity criteria (as shown in Equation ( 1)).Therefore, Equation (4) could be simplified to (5), Then, according to the orthogonal test result in Table 5, the multivariate Equation (5) were fitted using Origin program to identify the regression Formula (6).The R 2 values of Y σ c , Y C , Y F and Y D were 0.838, 0.965, 0.945 and 0.937, respectively.
Processes 2019, 7, x FOR PEER REVIEW 19 of 25 with previous studies [13].More importantly, their similarity ratios are the same in the similarity criteria (as shown in Equation ( 1)).Therefore, Equation ( 4) could be simplified to (5), , , , , Then, according to the orthogonal test result in Table 5, the multivariate Equation ( 5) were fitted using Origin program to identify the regression Formula (6).The R 2 values of c Y σ , YC, YF and YD were 0.838, 0.965, 0.945 and 0.937, respectively.
According to the physical and mechanical parameters of the raw coal, the five-variable primary equations can be solved using the Mmult and Minverse functions in the Excel program to obtain the corresponding composition (26.6% coarse aggregate, 62.2% fine aggregate, 8.5% gypsum, 1.7% petrolatum and 1% calcium carbonate).The gas absorption capacity of this material was estimated by Equation (3) to be a = 38.738.
Then, samples were made according to the determined composition and their indexes were tested in the same conditions.Table 9 shows the differences of each index of mechanical, deformation and gaseous properties between the natural coal and developed coal-like material.The difference in mechanical and gas adsorption/desorption parameters was very small (approximately 0.47 to 11.76%).Figure 12 is the stress-strain curve of the coal-like material, showing evolution characteristics of strength and deformation (comprehensive damage mode of brittle shearing and plasticity) very consistent to that of the natural coal.
According to the physical and mechanical parameters of the raw coal, the five-variable primary equations can be solved using the Mmult and Minverse functions in the Excel program to obtain the corresponding composition (26.6% coarse aggregate, 62.2% fine aggregate, 8.5% gypsum, 1.7% petrolatum and 1% calcium carbonate).The gas absorption capacity of this material was estimated by Equation (3) to be a = 38.738.
Then, samples were made according to the determined composition and their indexes were tested in the same conditions.Table 9 shows the differences of each index of mechanical, deformation and gaseous properties between the natural coal and developed coal-like material.The difference in mechanical and gas adsorption/desorption parameters was very small (approximately 0.47 to 11.76%).Figure 12 is the stress-strain curve of the coal-like material, showing evolution characteristics of strength and deformation (comprehensive damage mode of brittle shearing and plasticity) very consistent to that of the natural coal.As stated in Section 3.1.2,a reasonable ratio of aggregate particles of different sizes has an important influence on the material strength and elastic modulus.Meanwhile, the particle sizes of coal powders reported in Table 1 were different.To study the effect of particle size on the preparation and strength of samples, coal-like materials with aggregate/gypsum = 10:1, coarse (1-2 mm)/fine (0-0.5 mm) coal particles = 1:9, 3:7 and 5:5 were prepared.Their microstructures were examined by SEM (Figure 13).
On the one hand, the morphology diagrams (Figure 13a) of the three cut specimens show that, when the proportion of fine aggregate was high (coarse/fine coal particles = 1:9 and 3:7), the sample was relatively flat in morphology and relatively dense in structure, while a high proportion of coarse aggregate (coarse/fine coal particles = 5:5) resulted in an uneven sample morphology, a looser structure and a low strength (shown in Figure 7b).On the other hand, in specimen No. 1 shown in Figure 13b, the gypsum (needle shapes) were stuck together when 90% of the coal had a particle size of 0-0.5 mm.This indicates that too much fine coal particles makes it difficult to mix the sample evenly with the cementing agent, resulting in poor strength stability and low homogeneity of the sample, which will have an important impact on the analysis of test results.However, in specimen No. 3 the gypsum was uniformly attached to the surface of coal particles and the very small coal particles were also firmly bonded to the larger ones.Therefore, a large proportion of coarse coal particles is beneficial to the uniform distribution of cementing agent.When the ratio of coarse/fine coal particles was about 3:7, the cement and aggregate were well mixed, and the specimen also had a compact structure.As stated in Section 3.1.2,a reasonable ratio of aggregate particles of different sizes has an important influence on the material strength and elastic modulus.Meanwhile, the particle sizes of coal powders reported in Table 1 were different.To study the effect of particle size on the preparation and strength of samples, coal-like materials with aggregate/gypsum = 10:1, coarse (1-2 mm)/fine (0-0.5 mm) coal particles = 1:9, 3:7 and 5:5 were prepared.Their microstructures were examined by SEM (Figure 13).
On the one hand, the morphology diagrams (Figure 13a) of the three cut specimens show that, when the proportion of fine aggregate was high (coarse/fine coal particles = 1:9 and 3:7), the sample was relatively flat in morphology and relatively dense in structure, while a high proportion of coarse aggregate (coarse/fine coal particles = 5:5) resulted in an uneven sample morphology, a looser structure and a low strength (shown in Figure 7b).On the other hand, in specimen No. 1 shown in Figure 13b, the gypsum (needle shapes) were stuck together when 90% of the coal had a particle size of 0-0.5 mm.This indicates that too much fine coal particles makes it difficult to mix the sample evenly with the cementing agent, resulting in poor strength stability and low homogeneity of the sample, which will have an important impact on the analysis of test results.However, in specimen No. 3 the gypsum was uniformly attached to the surface of coal particles and the very small coal particles were also firmly bonded to the larger ones.Therefore, a large proportion of coarse coal particles is beneficial to the uniform distribution of cementing agent.When the ratio of coarse/fine coal particles was about 3:7, the cement and aggregate were well mixed, and the specimen also had a compact structure.
In addition, it can be seen that gypsum with its extremely small particle size may reduce the gap between particles and the porosity would be reduced to a certain extent, which helps to reduce the permeability of reconstituted coal.In addition, it can be seen that gypsum with its extremely small particle size may reduce the gap between particles and the porosity would be reduced to a certain extent, which helps to reduce the permeability of reconstituted coal.

Optimization of Sample Preparation Method and Deformation Characteristic Controlling Binder
Samples made by directly pressing the coal powder have lower strengths larger deformation than the natural coal.Consequently, studies in the past [10,13,18,19,21,29] usually used Portland cement to improve the strength of coal-like materials.However, the mechanical properties of cement keep changing over time.The molding and air curing time (usually 28 days) also have significant influence on the mechanical strengths of the sample and should be considered as additional influencing factors in the preparation process.As a result, the sample preparation process becomes more complex, which makes it hard to ensure the homogeneity and reproducibility.In order to solve this problem, this study used gypsum instead of cement and improved the sample preparation method as discussed in Section 2.3.3 (i.e., using quantitative compaction samples and self-designed molding and demolding equipment).However, gypsum as a binder is less effective than cement in improving the mechanical parameters of materials and so it is more suitable for the preparation of low-strength (0 to 2.5 MPa) coal-like materials.More importantly, gypsum and petrolatum were selected as the compound binder to adjust the deformation characteristics of coal-like materials.Similar materials with ductile to brittle deformation features could be obtained by adjusting the mass ratio of gypsum/petrolatum.

Porosity Adjustment Method
Among the properties summarized in Table 9, the porosity difference was large (11.73%) while the difference in adsorption and desorption parameters was smaller (3.63-7.41%).This indicated that the specific surface area was not very different between the natural coal and coal-like material.Rather, the gap between particles was larger in the coal-like materials.Studies have shown that the porosity of coal-like materials decreases with increasing molding pressure [12,48].Thus, the porosity could be adjusted by changing the molding pressure slightly but it could not be reduced by adding

Optimization of Sample Preparation Method and Deformation Characteristic Controlling Binder
Samples made by directly pressing the coal powder have lower strengths larger deformation than the natural coal.Consequently, studies in the past [10,13,18,19,21,29] usually used Portland cement to improve the strength of coal-like materials.However, the mechanical properties of cement keep changing over time.The molding and air curing time (usually 28 days) also have significant influence on the mechanical strengths of the sample and should be considered as additional influencing factors in the preparation process.As a result, the sample preparation process becomes more complex, which makes it hard to ensure the homogeneity and reproducibility.In order to solve this problem, this study used gypsum instead of cement and improved the sample preparation method as discussed in Section 2.3.3 (i.e., using quantitative compaction samples and self-designed molding and demolding equipment).However, gypsum as a binder is less effective than cement in improving the mechanical parameters of materials and so it is more suitable for the preparation of low-strength (0 to 2.5 MPa) coal-like materials.More importantly, gypsum and petrolatum were selected as the compound binder to adjust the deformation characteristics of coal-like materials.Similar materials with ductile to brittle deformation features could be obtained by adjusting the mass ratio of gypsum/petrolatum.

Porosity Adjustment Method
Among the properties summarized in Table 9, the porosity difference was large (11.73%) while the difference in adsorption and desorption parameters was smaller (3.63-7.41%).This indicated that the specific surface area was not very different between the natural coal and coal-like material.Rather, the gap between particles was larger in the coal-like materials.Studies have shown that the porosity of coal-like materials decreases with increasing molding pressure [12,48].Thus, the porosity could be adjusted by changing the molding pressure slightly but it could not be reduced by adding inorganic matter.For example, sand will decrease the coal powder content and gas adsorption/desorption capacity and the latter is the most important index in coal and gas outburst simulation tests [6].
3.5.4.Limitation and Future Direction to Improve the Similarity of Solid-Gas Coupling Coal-Like Materials It can be seen from Table 9 that the difference in gaseous parameters between coal-like material and natural coal was larger than that that in the mechanical parameters.This was mainly because the coal-like material contained inorganic matters (gypsum and calcium carbonate) and petrolatum that do not have gas adsorption properties.Future new binders that could simultaneously increase the strength and adsorption characteristics of coal-like materials will further improve the similarity and broaden the application prospects.

Application Prospects of Solid-Gas Coupling Coal-Like Materials
Quantitative simulation tests of coal and gas outbursts and CO 2 sequestration need to ensure that the coal-like material has highly consistent mechanical, deformation and gas adsorption/desorption characteristics with the natural coal [13,[19][20][21].In the past, such similarity was difficult to achieve [6,10,13,18,19,29], mainly because the deformation characteristics and gas adsorption/desorption characteristics were less considered.This study tested multiple indexes including the density, compressive strength, elastic modulus, cohesion, internal friction angle, porosity and gas adsorption and desorption parameters.The deformation features and damage modes of coal-like materials were also considered.Finally, samples with highly consistent stress-strain curves to the natural coal were obtained.
Since the tectonically deformed coal seam is extremely soft and fragile, it remains difficult to obtain natural coal samples that meet the test requirements and the related triaxial/long-term mechanical properties are still unclear.Existing studies have shown that the strength and deformation characteristics of reconstituted coal-like materials and natural coal were similar [20,23,27].In view of the highly consistent mechanical properties, the homogeneity and the reproducibility of the coal-like material developed in this study; the synthetic material could be used in place of natural coal sample to study the triaxial and rheological mechanical properties of soft tectonically deformed coal.The relevant experimental results will help guide the design of roadway support for soft coal seams, the stability control of the gas drainage borehole [49,50] and the prediction of coal and gas delayed outburst.

Conclusions
In this study, solid-gas coupling coal-like materials were developed according to the characteristics of low strength, strong plastic deformation and large gas adsorption capacities of natural TDC.Coal powder was used as the aggregate, gypsum and petrolatum as the composite binder, calcium carbonate as the additive and vertical ground stress as the molding pressure.The orthogonal method was used to design 25 groups of material ratio schemes.The physical and mechanical parameters, deformation features and gaseous constants were obtained by carrying out uniaxial compression, variable angle shearing and gas adsorption/desorption tests.The following conclusions were obtained.
(1) The range analysis indicates that the gypsum/petrolatum ratio played the main role in controlling the sample density, compressive strength, elastic modulus, cohesion and deformation characteristic.The coarse coal powder (1-2 mm)/fine coal powder (0-0.5 mm) ratio was mainly in control of the sample internal friction angle and also had an important influence on the compressive strength and elastic modulus.When the fine coal powder/aggregate ratio was about 70%, the sample had the highest strength.With the gypsum/petrolatum ratio increasing from 5:5, 6:4, 7:3, 8:2 to 9:1, the deformation characteristics of coal-like materials transitioned from ductile to brittle; and the failure modes included typical squeezing damage, wedge splitting, plastic shear failure, brittle shear and brittle fracturing.The content of coal powder played a decisive role in the gas adsorption and desorption characteristics of coal-like materials.An empirical formula between coal powder content and gas adsorption index was also established.
(2) The microstructures of specimens indicate that an even mixing with the cementing agent is difficult when there is a large ratio of fine coal particles, resulting in poor strength stability and low homogeneity of the sample.When the ratio of coarse/fine coal particles was about 3:7, the binder and aggregate were relatively well mixed and the specimen had a compact structure.
(3) A multivariate linear model was established between the mechanical characteristic indexes and their main influencing factors.Then, the desired material ratio was determined using this model and the developed coal-like material showed high similarity in its mechanical and physicochemical properties to the natural coal sample.Meanwhile, developing new binders that simultaneously improve sample strength and gas adsorption/desorption performance would help further reduce the difference between coal-like material and natural coal.
(4) The method for preparing solid-gas coupling coal-like material was optimized and the developed samples had good reproducibility and homogeneity.This study lays the foundation for future quantitative work in the physical simulation of coal and gas outbursts, the long-term safety study of CO 2 sequestration with coal in the lab and other experiments on the mechanical properties of soft coal seams.

Figure 1 .
Figure 1.The difference between this work and the previous studies.

Figure 2 .
Figure 2. The stress-strain curve of target natural coal.

Figure 3 .
Figure 3.The structure of natural coal under scanning electron microscope.

Figure 3 .
Figure 3.The structure of natural coal under scanning electron microscope.

Figure 5 .
Figure 5. Sample molding and demolding process.(a) The load-time curve during the compaction process of each layer of material.(b) The load-time curve during the compaction process of the whole sample.(c) Sample compression molding process.(d) Sample demolding process.

Figure 5 .
Figure 5. Sample molding and demolding process.(a) The load-time curve during the compaction process of each layer of material.(b) The load-time curve during the compaction process of the whole sample.(c) Sample compression molding process.(d) Sample demolding process.

Figure 7 .
Figure 7.The effect of different levels of key influencing factors on the physical and mechanical indexes: (a) density; (b) UCS; (c) EM; (d) cohesion; (e) IFA.

Figure 7 .
Figure 7.The effect of different levels of key influencing factors on the physical and mechanical indexes: (a) density; (b) UCS; (c) EM; (d) cohesion; (e) IFA.

Figure 8 .
Figure 8. Stress-strain curves of coal-like material specimens with different gypsum/petrolatum ratios.

Figure 9 .
Figure 9. (a) The different physical damage modes and (b) the corresponding schematic diagrams of the specimens.

Figure 9 .
Figure 9. (a) The different physical damage modes and (b) the corresponding schematic diagrams of the specimens.

1 Figure 9 .
Figure 9. (a) The different physical damage modes and (b) the corresponding schematic diagrams of the specimens.

Figure 10 .
Figure 10.The quantitative relation between coal powder content and the gas absorption constant, a and the desorption index, Δp.

Figure 10 .
Figure 10.The quantitative relation between coal powder content and the gas absorption constant, a and the desorption index, ∆p.

Figure 11 .
Figure 11.Relationship between modulus of elasticity and uniaxial compressive strength of samples.

Figure 11 .
Figure 11.Relationship between modulus of elasticity and uniaxial compressive strength of samples.

Figure 12 .
Figure 12.The stress-strain curves of coal-like material and natural coal.

3. 5 .
Discussion 3.5.1.Microscopic Analysis of the Effect of Pulverized Coal Particle Size on Sample Strength

Figure 12 .
Figure 12.The stress-strain curves of coal-like material and natural coal.

3. 5 .
Discussion 3.5.1.Microscopic Analysis of the Effect of Pulverized Coal Particle Size on Sample Strength

Figure 13 .
Figure 13.SEM images of three specimens composed with different coarse/fine coal particle ratios: (a) the morphology diagrams and (b) internal structures of the three cut specimens.

Figure 13 .
Figure 13.SEM images of three specimens composed with different coarse/fine coal particle ratios: (a) the morphology diagrams and (b) internal structures of the three cut specimens.

Table 1 .
Representative literature methods for preparing coal-like materials.
Processes 2019, 7, x FOR PEER REVIEW 5 of 25 Figure 1.The difference between this work and the previous studies.

Table 3 .
The 4 factors and 5 levels of the orthogonal test design.

Table 4 .
The 25 groups of material composition ratio schemes.

Table 5 .
The physical and mechanical parameters of each group.

Table 6 .
Range analysis of the physical and mechanical indexes.

Table 6 .
Range analysis of the physical and mechanical indexes.

Table 7 .
Porosity of the samples.

Table 8 .
Adsorption and desorption parameters of samples.

Table 9 .
Difference analysis between natural coal and determined coal-like material.

Table 9 .
Difference analysis between natural coal and determined coal-like material.