Under the action of dynamic loadings such as earthquakes and volcanic activities, the mechanical properties of gas-hydrate-bearing sediments will deteriorate, leading to a decrease in the stability of hydrate reservoirs and even inducing geological disasters such as submarine landslides. In order to study the effect of dynamic loading on the mechanical properties of hydrate sediments, triaxial compression tests of numerical specimens were carried out by using particle flow code (PFC2D), and the macro-meso mechanical behaviors of specimens were investigated. The results show that the loading frequency has a small effect on the stiffness of the hydrate sediment, while it has a large effect on the peak strength. The peak strength increases and then decreases with the increase in loading frequency. Under the same loading frequency, the peak strength of the hydrate sediment increases with the increase in loading amplitude, and the stiffness of the specimen decreases with the increase in loading amplitude. The maximum shear expansion of the specimen changes with the movement of the phase change point and the rearrangement of the particles. The maximum shear expansion of the specimen changes with the movement of the phase change point and the change of the bearing capacity of the particles after the rearrangement, and the more forward the phase change point is, the stronger the bearing capacity of the specimen in the plastic stage. The shear dilatancy angle and the shear dilatancy amount both increase linearly with the increase in loading amplitude. The influence of loading frequency and amplitude on the contact force chain, displacement, crack expansion, and the number of cementation damage inside the sediment is mainly related to the average axial stress to which the specimen is subjected, and the number of cracks and cementation damage of the sediment specimen increases with the increase in the average axial stress to which the sediment specimen is subjected. As the rate of cementation damage increases, the distribution of shear zones becomes more obvious.

Studying the mechanical characteristics of hydrate-bearing sediments (HBS) contributes to the comprehensive understanding of the mechanical behavior in environments with natural gas hydrate (NGH) occurrences. Simultaneously, the distribution patterns of hydrates significantly influence the strength, deformation, and stability of HBS. Therefore, this paper employs particle flow code (PFC) to conduct biaxial discrete element simulations on specimens of HBS with different hydrate distribution patterns, revealing the macroscale–mesoscale mechanical properties, evolution patterns, and destructive mechanisms. The results indicate that the strain-softening behavior of HBS specimens strengthens with the increase in hydrate layer thickness, leading to higher peak strength and E₅₀ values. During the gradual movement of the hydrate layer position (A_y) from both ends to the center of the specimen (A_y = 0.40 mm → A_y = 20 mm), the strain-softening behavior weakens. However, when A_y = 20 mm, the specimen exhibits evident strain-softening behavior again. Moreover, with an increase in the angle between the hydrate layer and the horizontal direction (α) greater than 20°, the peak strength of the specimen increases, while E₅₀ shows an overall decreasing trend. The influence of axial loads on the hydrate layer in specimens varies with α, with larger contact forces and fewer cracks observed for higher α values.

In order to study the macro-meso shear mechanical characteristics of natural gas hydrate-bearing sediments, the direct shear simulations of natural gas hydrate-bearing sediment specimens with different saturations under different normal stress boundary conditions were carried out using the discrete element simulation program of particle flow, and the macro-meso shear mechanical characteristics of the specimens and their evolution laws were obtained, and their shear damage mechanisms were revealed. The results show that the peak intensity of natural gas hydrate-bearing sediments increases with the increase in normal stress and hydrate saturation. Hydrate particles and sand particles jointly participate in the formation and evolution of the force chain, and sand particles account for the majority of the force chain particles and take the main shear resistance role. The number of cracks produced by shear increases with hydrate saturation and normal stress. The average porosity in the shear zone shows an evolutionary pattern of decreasing and then increasing during the shear process.