In-depth research on the mechanical properties and constitutive models of gas hydrate-bearing sediments (GHBSs) is fundamental for achieving efficient hydrate exploration and geological disaster prevention. In the current study, a bounding surface model for GHBSs is developed based on the principle of thermodynamics. By choosing an appropriate dissipation function and free energy function, a yield surface function containing three shape parameters can be obtained. Considering the filling and bonding effects of hydrates, and introducing the hydrate strength evolution parameter, a thermodynamics-based bounding surface model for GHBSs is established using a non-associated flow rule. Then, the explicit substeping scheme with error control is implemented to develop a UMAT subroutine for the proposed model and integrated into the ABAQUS. Compared with the drained monotonic triaxial shear data indicates that the proposed model can adequately capture the shear behaviors of sandy, silty sandy, and clay-silty GHBSs under different stress levels and saturations. In addition, the model demonstrates good applicability and feasibility in undrained cyclic triaxial shear tests and boundary value problem analysis.

Shear strength of hydrate-bearing sediment is an essential parameter for assessing landslide potential of hydrate reservoirs under exploration conditions. However, the characteristics and simulation of this shear strength under varying dissociation conditions have not been thoroughly investigated. To this end, a series of triaxial compression tests were first carried out on sediments with varying initial hydrate saturations along dissociation pathways. Combining measured data with microscale analysis, the underlying mechanism for the evolution of shear strength in hydrate-bearing sediment was studied under varying partial dissociation pathways. Moreover, a shear strength model for hydrate-bearing sediment was proposed, taking into account the hydrate saturation and the unhydrated water content. Apart from the parameters derived from the hydrate characteristic curve, only one additional model parameter is required. The proposed model was validated using measured data on hydrate sediments. The results indicate that the proposed model can effectively capture the shear strength behavior of hydrate-bearing sediment under varying dissociation paths. Finally, a sensitivity analysis of the model parameters was conducted to characterize the proposed model. (c) 2025 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/ 4.0/).

The mechanical behavior of Methane Hydrate-Bearing Sediment (MHBS) is essential for the safe exploitation of Methane Hydrate (MH). In particular, the pore size and physicochemical characteristics of MHBS significantly influence its mechanical behavior, especially in clayey grain-cementing type MHBS. This study employs the Distinct Element Method (DEM) to investigate both the macroscopic and microscopic mechanical behavior of clayey grain-cementing type MHBS, focusing on variations in pore size and physicochemical characteristics. To accomplish this, we propose a Thermo-Hydro-Mechanical-Chemical-Soil Characteristics (THMCS) DEM contact model that incorporates the effects of pore size and physicochemical characteristics on the strength and modulus of MH. This THMCS model is validated using experimental data available in the literature. Using the proposed contact model, we conducted a series of investigations to explore the mechanical behavior of MHBS under conventional loading paths, including isotropic and drained triaxial tests using the DEM. The numerical results indicate that smaller pore sizes and lower water content-key physicochemical characteristics resulting from variations in electrochemical properties and the intensity of the electric field-can lead to reduced shear strength and stiffness due to the increased breakage of aggregates and weakened cementation. Additionally, heating was found to further accelerate the process of structural damage in MHBS.

The mechanical properties of hydrate-bearing sediment (HBS) are significantly influenced by dynamic loads, which is an engineering safety issue and deserves serious attention in the exploitation of hydrates. However, relevant research is scarce. Considering the safety of mining and the purpose of establishing an optimal mining strategy, studying the dynamic mechanical properties of HBS under cyclic loading is fundamental for the safe exploitation. In this study, the effects of deviator stress amplitude on the mechanical properties of the HBS with different hydrate saturation are investigated from the particle-scale. The results show that the number of cycles decreases with increased deviator stress amplitude and hydrate saturation, and these factors also affect the failure mode of HBS. The non-cooperative displacement of the particles inside the specimen causes the deformation of the HBS specimen. In addition, from the perspective of fabric evolution, high hydrate saturation corresponds to a larger initial coordination number. The deviator stress amplitude governs the variation of the coordination number, and the contact force distribution reveals significant anisotropy in HBS. The energy analysis shows that during the deformation of HBS, the input energy is mainly dissipated unidirectionally through the dissipated energy, which indicates that HBS is a strongly dissipative system.

The mechanics of methane hydrate-bearing sediments (MHBS) have been broadly investigated over recent years in the context of methane-gas production or climate-change effects. Their mechanical investigation has mainly been carried out using models constructed from experimental data obtained for laboratory-formed MHBS. Along with the dominant effects of hydrate saturation and morphology within the host soil pores, this study recognizes the effective pressure at which the hydrate is formed as a key factor in the MHBS mechanics. A state-of-the-art experimental study has been conducted in order to isolate the effect of the hydrate formation pressure, for use as a model parameter. Two generalized mechanical prediction models that incorporate the effect of the hydrate formation pressure are developed in this work: (a) an analytical shear strength prediction, and (b) an empiric graphical model for predicting volumetric changes along a given stress path. The models are related to a novel data representation which enables the analysis of a few individual test outcomes as a whole, through a volume-change mapping that describes the complex influence of the volumetric effect of hydrate in MHBS, under combined hydrostatic and deviatoric loading scenarios. In this study, we delve into a specific configuration of hydrate morphology, hydrate saturation, and host soil type, enabling a distinctive fundamental geotechnical investigation and the development of a conceptual modeling approach. The paper describes the approaches by which the MHBS properties can be extracted for other MHBS samples (than those examined in this work) having different host soils and hydrate pore-space morphologies.