Mechanical Performance of Methane Hydrate–Coal Mixture

: Understanding mechanical behaviors of a methane hydrate–coal mixture are important for its associated application in coal and gas outburst prevention. A triaxial compression apparatus for a specimen of the methane hydrate–coal mixture was developed to measure its strength and deformation properties. Triaxial compression tests were performed on coal briquette samples under different conﬁning pressures and methane hydrate saturation. Strain softening behavior of the methane hydrate–coal mixture can be identiﬁed for all the specimens under various conditions. The larger the methane hydrate saturation or conﬁning pressure, the larger the peak strength, elastic modulus and peak strain. Mathematical relationships for correlating the peak strength with the methane hydrate saturation or conﬁning pressure were also proposed.


Introduction
Coal and gas outburst will eject a large amount of coal and gas very shortly [1]. Disaster will become more serious for deeper coal mining [2,3]. So far, four factors shown in Figure 1 are found to have effects on the gas outbursts [4]. Among these factors, the outburst is very sensitive to gas, considered by most of scholars [5]. Therefore, techniques to prevent this kind of disaster contain: (1) rising up permeability of coal; and (2) decreasing stress in coal seams [6]. Commonly, gas drainage and hydraulic fracturing are adopted to prevent coal and gas outbursts, as well as directional drilling, or water infusion [2,7,8]. Additionally, if a coal seam is strengthened, outbursts can be eliminated [9]. For instance, Wu et al. successfully synthesized artificial methane hydrate (MH), and they used a hydrate method for outburst prevention [10]. Mechanism of the method is illustrated in Figure 2. Coal mine and gas, such as methane, can generate MH with water at a temperature of 0.55 • C and a pressure of 2.77 MPa, and the formed MH is energy-dense. Then, it is possible to reduce gas content and gas pressure significantly, which can decrease or eliminate risks of coal and gas outbursts. Consequently, evaluating the mechanical properties of MH-coal mixtures is very important when applied to outburst prevention.  Schematic diagram for the hydrate-based method for preventing coal and gas outbursts [11].
Up to now, there are only very limited literatures available relevant to mechanical behaviors of MH-coal mixtures [11]. Some investigations have been made related to mechanical behaviors of gas hydrate-bearing sediments. The relationships were obtained between strength parameters, temperature and confining pressure for hydrate-bearing sediments [12,13]. It was found that strength and deformation parameters increase with increasing hydrate saturation for hydrate-bearing sediments [14][15][16]. Low-temperature and high-pressure triaxial apparatus was used to study effects of hydrate saturation, confining pressure and temperature on the strength and deformation parameters for hydrate-bearing sediments [17][18][19][20][21][22][23][24][25]. Recently, direct shear tests were conducted to study the effects of hydrate saturation, temperatures as well as loading rates on mechanical behaviors for hydrate-bearing sediments [26]. In above-mentioned studies, it is found that both hydrate saturation and confining pressures have strong effects on strength and stiffness indices of gas hydrate-bearing sediments.
For MH-coal mixture samples, some preliminary studies were made related to strength characteristics of coal before and after gas hydrate formation under various confining pressures (1, 2 and 3 MPa) and hydrate saturation (25%, 50% and 80%) [11]. The results showed the strength of coal could be enhanced with gas hydrate formation. Additionally, it was found that mechanical properties were strengthened for MH-coal mixture samples, when there were increases in the hydrate  Schematic diagram for the hydrate-based method for preventing coal and gas outbursts [11].
Up to now, there are only very limited literatures available relevant to mechanical behaviors of MH-coal mixtures [11]. Some investigations have been made related to mechanical behaviors of gas hydrate-bearing sediments. The relationships were obtained between strength parameters, temperature and confining pressure for hydrate-bearing sediments [12,13]. It was found that strength and deformation parameters increase with increasing hydrate saturation for hydrate-bearing sediments [14][15][16]. Low-temperature and high-pressure triaxial apparatus was used to study effects of hydrate saturation, confining pressure and temperature on the strength and deformation parameters for hydrate-bearing sediments [17][18][19][20][21][22][23][24][25]. Recently, direct shear tests were conducted to study the effects of hydrate saturation, temperatures as well as loading rates on mechanical behaviors for hydrate-bearing sediments [26]. In above-mentioned studies, it is found that both hydrate saturation and confining pressures have strong effects on strength and stiffness indices of gas hydrate-bearing sediments.
For MH-coal mixture samples, some preliminary studies were made related to strength characteristics of coal before and after gas hydrate formation under various confining pressures (1, 2 and 3 MPa) and hydrate saturation (25%, 50% and 80%) [11]. The results showed the strength of coal could be enhanced with gas hydrate formation. Additionally, it was found that mechanical properties were strengthened for MH-coal mixture samples, when there were increases in the hydrate Up to now, there are only very limited literatures available relevant to mechanical behaviors of MH-coal mixtures [11]. Some investigations have been made related to mechanical behaviors of gas hydrate-bearing sediments. The relationships were obtained between strength parameters, temperature and confining pressure for hydrate-bearing sediments [12,13]. It was found that strength and deformation parameters increase with increasing hydrate saturation for hydrate-bearing sediments [14][15][16]. Low-temperature and high-pressure triaxial apparatus was used to study effects of hydrate saturation, confining pressure and temperature on the strength and deformation parameters for hydrate-bearing sediments [17][18][19][20][21][22][23][24][25]. Recently, direct shear tests were conducted to study the effects of hydrate saturation, temperatures as well as loading rates on mechanical behaviors for hydrate-bearing sediments [26]. In above-mentioned studies, it is found that both hydrate saturation and confining pressures have strong effects on strength and stiffness indices of gas hydrate-bearing sediments.
For MH-coal mixture samples, some preliminary studies were made related to strength characteristics of coal before and after gas hydrate formation under various confining pressures (1,2 Energies 2018, 11, 1562 3 of 14 and 3 MPa) and hydrate saturation (25%, 50% and 80%) [11]. The results showed the strength of coal could be enhanced with gas hydrate formation. Additionally, it was found that mechanical properties were strengthened for MH-coal mixture samples, when there were increases in the hydrate saturation as well as the confining pressure. Therefore, laboratory triaxial apparatus was used to investigate the strength properties of the MH-coal mixtures under various confining pressures (4, 5 and 6 MPa) and high MH saturation (from 50% to 80%). Furthermore, the mathematical expressions were obtained to quantitively predict the strength parameters of the MH-coal mixtures under specified conditions.

Apparatus
The schematic diagram of the apparatus used is shown in Figure 3. The system mainly involved two parts. One was for hydrate formation and the other one was for triaxial compression. The detailed information can be referred to Gao et al. [11].
Energies 2018, 11, x FOR PEER REVIEW 3 of 13 saturation as well as the confining pressure. Therefore, laboratory triaxial apparatus was used to investigate the strength properties of the MH-coal mixtures under various confining pressures (4, 5 and 6 MPa) and high MH saturation (from 50% to 80%). Furthermore, the mathematical expressions were obtained to quantitively predict the strength parameters of the MH-coal mixtures under specified conditions.

Apparatus
The schematic diagram of the apparatus used is shown in Figure 3. The system mainly involved two parts. One was for hydrate formation and the other one was for triaxial compression. The detailed information can be referred to Gao et al. [11].

Specimen Preparation
The briquette specimens were prepared using the outburst coal from Taoshan Coal Mine, Heilongjiang Longmay Mining Holding Group Co., Ltd. (Harbin, China). Samples with 50 mm in diameter and 100 mm in length were finally obtained after compression. We have prepared six samples for the mercury intrusion porosimetry (MIP) tests. Figure 4 shows the variations of the accumulated intruded mercury volume with the pressure. Based on the pore classification in literature [27], we calculated the pore volume distribution and Table 1 showed the detailed pore distribution condition. It is clear that for tested coal samples, Macro-pores and Meso-pores accounted for most part of the pore volume. For six MIP samples, the average total pore volume was 0.1322 mL/g. As for mechanical tests, samples used were assumed to have the same porosity in this work, due to the same briquetting pressure and time for every sample prepared with a coal particle range from 60-80 mesh [28]. Purity of methane was 99.99%. Water was distilled. The detailed preparation procedure can be found in Reference [11].

Specimen Preparation
The briquette specimens were prepared using the outburst coal from Taoshan Coal Mine, Heilongjiang Longmay Mining Holding Group Co., Ltd. (Harbin, China). Samples with 50 mm in diameter and 100 mm in length were finally obtained after compression. We have prepared six samples for the mercury intrusion porosimetry (MIP) tests. Figure 4 shows the variations of the accumulated intruded mercury volume with the pressure. Based on the pore classification in literature [27], we calculated the pore volume distribution and Table 1 showed the detailed pore distribution condition. It is clear that for tested coal samples, Macro-pores and Meso-pores accounted for most part of the pore volume. For six MIP samples, the average total pore volume was 0.1322 mL/g. As for mechanical tests, samples used were assumed to have the same porosity in this work, due to the same briquetting pressure and time for every sample prepared with a coal particle range from 60-80 mesh [28]. Purity of methane was 99.99%. Water was distilled. The detailed preparation procedure can be found in Reference [11].

Formation of Methane Hydrates in Coal Samples and Experimental Procedure
Proper thermodynamic conditions are needed for MH formation, which are low temperature or high pressure [29]. Calculation of MH saturation can be referred to Reference [11]. Note that the volume of samples was assumed constant during hydrate formation process, and water in the coal sample was completely transformed to hydrate with methane under a specified temperature and pressure [13,26,[30][31][32]. The hydrate formation and triaxial compression process are schematically shown in Figure 5. We chose an appropriate thermodynamic condition (0.5 °C, 4 MPa) to ensure hydrate formation in coal samples. In Figure 5, the left part in blue above the MH phase boundary represents a stable region for the hydrate formation.

Formation of Methane Hydrates in Coal Samples and Experimental Procedure
Proper thermodynamic conditions are needed for MH formation, which are low temperature or high pressure [29]. Calculation of MH saturation can be referred to Reference [11]. Note that the volume of samples was assumed constant during hydrate formation process, and water in the coal sample was completely transformed to hydrate with methane under a specified temperature and pressure [13,26,[30][31][32]. The hydrate formation and triaxial compression process are schematically shown in Figure 5. We chose an appropriate thermodynamic condition (0.5 • C, 4 MPa) to ensure hydrate formation in coal samples. In Figure 5, the left part in blue above the MH phase boundary represents a stable region for the hydrate formation.

Formation of Methane Hydrates in Coal Samples and Experimental Procedure
Proper thermodynamic conditions are needed for MH formation, which are low temperature or high pressure [29]. Calculation of MH saturation can be referred to Reference [11]. Note that the volume of samples was assumed constant during hydrate formation process, and water in the coal sample was completely transformed to hydrate with methane under a specified temperature and pressure [13,26,[30][31][32]. The hydrate formation and triaxial compression process are schematically shown in Figure 5. We chose an appropriate thermodynamic condition (0.5 °C, 4 MPa) to ensure hydrate formation in coal samples. In Figure 5, the left part in blue above the MH phase boundary represents a stable region for the hydrate formation.   (1) Tightness of experimental system checking. After placing moist coal sealed with thermal contractible tubes on pedestal, the experimental system was injected with N 2 using a pressure pump; then some solutions of foaming agents were applied on the pipeline surface to observe whether bubbles appear. If no bubbles were identified, a gas intake valve was closed to ensure the sealing of the system, and tightness condition was regarded as fine if the gas pressure remained constant for approximately 12 h. Finally, the confining pressure and the methane gas pressure were increased to 0.5 MPa and 0.3 MPa, respectively. Then, methane was expelled repeatedly to ensure no air exited in the coal specimen and the pipeline.
(2) MH formation in the coal specimen. After the tightness of the system was checked, the moist sealed coal specimen was placed inside the cylinder chamber, shown as point (a) in Figure 5, and was injected with methane gas. For each test, methane gas was injected with a pressure around 4 MPa for approximately 24 h to ensure gas adsorption, shown as point (b) in Figure 5. After the methane was injected, the pressure and the temperature were adjusted to the designed condition (a temperature of 0.5 • C, a gas pressure of 4 MPa) for gas hydrate formation, indicated as the point c in Figure 5. Hydrate formation was considered to complete when no obvious gas pressure varied for approximately 6 h. All the hydrate formation process continued approximately 48 h after 24 h methane adsorption.
(3) Triaxial compression tests on the MH-coal mixture. When MHhydrate formation was completed, the axial loading was conducted manually. The gas pressure was maintained at around 4.0 MPa, the confining pressures were kept at 4.0, 5.0, and 6.0 MPa, and the temperature remained at 0.5 • C. Variables of triaxial compression tests are summarized in Table 2.   Figure 6 shows a stress-strain curve of the tested coal sample. This curve can be divided into four stages. OA represents stage 1, namely a linear elastic region. In this stage, the particles inside the coal pores were squeezed without generation of cracks. AB is the post-yielding and pre-failure stage. In this region, a shear movement occurred between the particles, leading to the formation of cracks. However, the MH-coal mixture could still sustain some further deformation without losing its capacity. BC shows the post-peak regime, with the propagation and the extension of cracks. Therefore, the stress of the coal sample gradually decreased. CD illustrates the residual state of the mixture. Figure 7 illustrates the results from triaxial compression tests for specimens with different hydrate saturation (50%, 60%, 70% and 80%), subjected to various confining pressures (4, 5 and 6 MPa). It is clear that all stress-strain curves exhibit a similar trend. Deviator stress reached a maximum value at Energies 2018, 11, 1562 6 of 14 a certain strain between 2% and 5%. Then, strain softening happened. The residual state happened till the strain exceeded 10%. Afterwards, deviator stress reached the residual state. In Figure 7a,b, for S h of 50% and 60%, the confining pressure of 6 MPa resulted in higher stiffness in the linear elastic region. For the specimen with S h of 70% and 80% (Figure 7c,d), the strain softening behavior was more obvious under the high confining pressure of 6 MPa. This implies that the confining pressure acts as a vital part in the mechanical behavior of the tested coal samples.

Triaxial Compression Test
Energies 2018, 11, x FOR PEER REVIEW 6 of 13 maximum value at a certain strain between 2% and 5%. Then, strain softening happened. The residual state happened till the strain exceeded 10%. Afterwards, deviator stress reached the residual state. In Figure 7a,b, for Sh of 50% and 60%, the confining pressure of 6 MPa resulted in higher stiffness in the linear elastic region. For the specimen with Sh of 70% and 80% (Figure 7c,d), the strain softening behavior was more obvious under the high confining pressure of 6 MPa. This implies that the confining pressure acts as a vital part in the mechanical behavior of the tested coal samples.      maximum value at a certain strain between 2% and 5%. Then, strain softening happened. The residual state happened till the strain exceeded 10%. Afterwards, deviator stress reached the residual state. In Figure 7a,b, for Sh of 50% and 60%, the confining pressure of 6 MPa resulted in higher stiffness in the linear elastic region. For the specimen with Sh of 70% and 80% (Figure 7c,d), the strain softening behavior was more obvious under the high confining pressure of 6 MPa. This implies that the confining pressure acts as a vital part in the mechanical behavior of the tested coal samples.       Moreover, a remarkable increase in the strength was observed with the increasing saturation and confining pressure. This is due to that MH filled in the voids within the coal sample enhanced its structure and cementation, while higher confining pressure generated higher densification inside the specimen.

Correlations between Peak Strength and Saturation, and between Peak Strength and Confining Pressure
Relationships between the peak strength (failure strength) (σ d(peak) ) and the hydrate saturation (S h ), along with relationships between the peak strength and the confining pressure, are displayed in Figures 8 and 9. Table 3 summarizes the typical strength and the deformation indices from tests. It was observed in Figure 8 that saturation had a significant impact on the peak strength of specimen. For each confining pressure, the peak strength increased with the increased saturation. At the confining pressure of 4 MPa, when S h was increased from 50% to 80%, the peak strength increased from 5.86 MPa to 8.68 MPa. In addition, the relationship between the peak strength and the hydration saturation can be captured by a quadratic polynomial function. This trend is similar to the finding by Miyazaki et al. [16], who discovered the peak strength increases at higher saturation. For high-saturation samples, hydrate particles assemble collectively to increase the cementation of the specimen [26]. Waite et al. [29] revealed that bulk density of hydrate-bearing sediments also increases at higher saturation, indicating that the bonding structure inside the specimen is enhanced. Moreover, a remarkable increase in the strength was observed with the increasing saturation and confining pressure. This is due to that MH filled in the voids within the coal sample enhanced its structure and cementation, while higher confining pressure generated higher densification inside the specimen.

Correlations between Peak Strength and Saturation, and between Peak Strength and Confining Pressure
Relationships between the peak strength (failure strength) (σd(peak)) and the hydrate saturation (Sh), along with relationships between the peak strength and the confining pressure, are displayed in Figures 8 and 9. Table 3 summarizes the typical strength and the deformation indices from tests. It was observed in Figure 8 that saturation had a significant impact on the peak strength of specimen. For each confining pressure, the peak strength increased with the increased saturation. At the confining pressure of 4 MPa, when Sh was increased from 50% to 80%, the peak strength increased from 5.86 MPa to 8.68 MPa. In addition, the relationship between the peak strength and the hydration saturation can be captured by a quadratic polynomial function. This trend is similar to the finding by Miyazaki et al. [16], who discovered the peak strength increases at higher saturation. For highsaturation samples, hydrate particles assemble collectively to increase the cementation of the specimen [26]. Waite et al. [29] revealed that bulk density of hydrate-bearing sediments also increases at higher saturation, indicating that the bonding structure inside the specimen is enhanced.    Moreover, a remarkable increase in the strength was observed with the increasing saturation and confining pressure. This is due to that MH filled in the voids within the coal sample enhanced its structure and cementation, while higher confining pressure generated higher densification inside the specimen.

Correlations between Peak Strength and Saturation, and between Peak Strength and Confining Pressure
Relationships between the peak strength (failure strength) (σd(peak)) and the hydrate saturation (Sh), along with relationships between the peak strength and the confining pressure, are displayed in Figures 8 and 9. Table 3 summarizes the typical strength and the deformation indices from tests. It was observed in Figure 8 that saturation had a significant impact on the peak strength of specimen. For each confining pressure, the peak strength increased with the increased saturation. At the confining pressure of 4 MPa, when Sh was increased from 50% to 80%, the peak strength increased from 5.86 MPa to 8.68 MPa. In addition, the relationship between the peak strength and the hydration saturation can be captured by a quadratic polynomial function. This trend is similar to the finding by Miyazaki et al. [16], who discovered the peak strength increases at higher saturation. For highsaturation samples, hydrate particles assemble collectively to increase the cementation of the specimen [26]. Waite et al. [29] revealed that bulk density of hydrate-bearing sediments also increases at higher saturation, indicating that the bonding structure inside the specimen is enhanced.       Figure 9 shows a larger confining pressure leads to a higher peak strength. Taken S h of 50% as an example. σ d(peak) could increase from 5.86 MPa to 9.38 MPa. The peak strength was increased because the higher confining pressure limited the movement of coal particles, resulting in more energy to break the cementation between particles. Furthermore, a linear relationship can be built between the peak strength and the confining pressure.
As shown in Figure 8, a relationship between σ d(peak) and S h was established: where c, b, a are the fitting parameters, which should be related to the confining pressure, based on the observation in Figure 9. By regression analysis, correlations between c, b, a and confining pressures σ 3 were obtained as follows: Combining Equations (1)-(4), the following formula was achieved for characterizing the impact of the saturation and the confining pressure upon the peak strength:

Correlation of Residual Strength with Hydrate Saturation and Confining Pressure
In underground mine engineering, the post-peak mechanical properties, such a residual strength of rock or coal mass, are essential for design of the project. Figure 10 plots the relationship between the residual strength (σ d(res) ) and the hydration saturation (S h ) for the MH-coal mixture. It was obviously identified that σ d(res) increased with S h . σ d(res) increased from 4.43 MPa to 7.67 MPa for the specimen under 4 MPa. Moreover, the quadratic polynomial function could be adopted to capture the relationship between σ d(res) and S h . Figure 11 presents that σ d(res) increased linearly with σ 3 . σ d(res) increased from 4.43 MPa to 7.40 MPa for the MH-coal mixture when S h increased from 50% to 80%. As shown in Figure 12, σ d(res) increased linearly with σ d(peak): Equation (6)

Correlations of Peak Strain with Hydrate Saturation and Confining Pressure
A peak strain (failure strain) is the axial strain corresponding to the peak strength. In general, ε p increased with S h (Figure 13). For the specimen under 6 MPa, when S h was increased from 50% to 80%, ε p climbed rapidly from 2.838% to 5.124%. It indicated that higher hydrate saturation led to the failure of specimen at a larger strain level. It can be seen from Figure 14 that ε p increased with the increased σ 3 obviously. For instance, when σ 3 was increased from 4 MPa to 6 MPa, ε p increased from 3.176% to 5.124% for the specimen with S h of 80%.

Correlations of Peak Strain with Hydrate Saturation and Confining Pressure
A peak strain (failure strain) is the axial strain corresponding to the peak strength. In general, εp increased with Sh ( Figure 13). For the specimen under 6 MPa, when Sh was increased from 50% to 80%, εp climbed rapidly from 2.838% to 5.124%. It indicated that higher hydrate saturation led to the failure of specimen at a larger strain level.
It can be seen from Figure 14 that εp increased with the increased σ3 obviously. For instance, when σ3 was increased from 4 MPa to 6 MPa, εp increased from 3.176% to 5.124% for the specimen with Sh of 80%.

Correlations of Elastic Modulus with Hydrate Saturation and Confining Pressure
Elastic modulus (Et) was determined as a ratio of 50% of the peak stress to the corresponding axial strain. The variation of Et with σ3 is plotted in Figure 15. It states clearly that Et increased with σ3. For Sh of 50%, when σ3 ranged from 4 MPa to 6 MPa, Et climbed from 343.65 MPa up to 514.02 MPa. Et of the specimen with Sh of 80% was about twice of that of specimen with Sh of 50% under 6 MPa. A similar phenomenon can also be found in the effect of the confining pressure on the strength. Therefore, it indicates that the mechanical performance (both strength and stiffness) of the MH-coal mixture was significantly influenced by the hydrate saturation and the confining pressure.

Correlations of Elastic Modulus with Hydrate Saturation and Confining Pressure
Elastic modulus (E t ) was determined as a ratio of 50% of the peak stress to the corresponding axial strain. The variation of E t with σ 3 is plotted in Figure 15. It states clearly that E t increased with σ 3 . For S h of 50%, when σ 3 ranged from 4 MPa to 6 MPa, E t climbed from 343.65 MPa up to 514.02 MPa. E t of the specimen with S h of 80% was about twice of that of specimen with S h of 50% under 6 MPa. A similar phenomenon can also be found in the effect of the confining pressure on the strength. Therefore, it indicates that the mechanical performance (both strength and stiffness) of the MH-coal mixture was significantly influenced by the hydrate saturation and the confining pressure.

Discussion
The presented experimental studies are limited to coal briquette samples due to difficulty in obtaining hydrate saturation in coal samples. For each coal sample, we assumed the same porosity due to the same molding method, and water was completely converted to MH. Hence, the theoretical hydrate saturation was higher than actual saturation, which could not be measured due to limitations of the apparatus. Additionally, axial and confining pressures were loaded manually, which inevitably resulted in errors of mechanical properties. The manual loading mode led to the fact that experiments could only be done at a low confining pressure. Moreover, it can be seen from Equations (5) and (6) that the peak strength was closely related to the saturation and the confining pressure, and the residual strength increased linearly with the increased peak strength. Therefore, a value of hydrate saturation is very important for reliably predicting the peak strength and the residual strength.

Discussion
The presented experimental studies are limited to coal briquette samples due to difficulty in obtaining hydrate saturation in coal samples. For each coal sample, we assumed the same porosity due to the same molding method, and water was completely converted to MH. Hence, the theoretical hydrate saturation was higher than actual saturation, which could not be measured due to limitations of the apparatus. Additionally, axial and confining pressures were loaded manually, which inevitably resulted in errors of mechanical properties. The manual loading mode led to the fact that experiments could only be done at a low confining pressure. Moreover, it can be seen from Equations (5) and (6) that the peak strength was closely related to the saturation and the confining pressure, and the residual strength increased linearly with the increased peak strength. Therefore, a value of hydrate saturation is very important for reliably predicting the peak strength and the residual strength. However, the predicted peak strength and residual strength may be higher than their actual strengths, due to higher theoretical saturation. This phenomenon will be more obvious under 6 MPa. Furthermore, the predicted mathematical model was only limited to high hydrate saturation and low confining pressures. In future studies, a more rational porosity measurement method will be investigated and thus corresponding hydrate saturation can be calculated properly. The loading system will be loaded at a constant loading rate, and experiments can be conducted under a relatively high confining pressure.

Conclusions
This study carried out laboratory simulation of MH formation in coal samples and corresponding mechanical tests at various confining pressures, as well as saturation. Strength and deformation characteristics were discussed through the obtained stress-strain curves. Strain softening was identified for all tested specimens. For some cases, like hydrate saturation of 70% and 80%, the strain softening behavior was more obvious under a high confining pressure of 6 MPa. Mechanical properties of the MH-coal mixture were enhanced due to the increasing saturation, together with confining pressures. Higher hydrate saturation or confining pressures also led to a larger failure strain of the specimen. A linear relationship was established between the residual strength and the peak strength of the MH-coal mixture. Moreover, both the peak and residual strengths showed clear linear relationships with the confining pressure, and the quadratic polynomial relationships between the peak strength and the hydrate saturation, as well as the quadratic relationships between the residual strength and the hydrate saturation were identified. Finally, comprehensive mathematical expressions for the failure strength considering both the saturation and the confining pressure were also proposed.