The mechanical characteristics of gas hydrate-bearing sediments (HBS) are important for evaluating reservoir stability. The interbedded formation of HBS is common in target mining reservoirs. Existing studies on the triaxial mechanical properties of HBS are primarily based on homogeneous and isotropic samples. Therefore, the stress–strain law of the target mining reservoirs cannot be predicted accurately. In this study, a series of sediment models with interlayers of coarse and fine mineral grains were established based on the PFC3D code, and the influence of the layered distribution characteristics of sediment particles and hydrates on the macroscopic mechanical behaviour of the reservoir was comprehensively analysed. The triaxial compression simulation results indicate that the peak strength, deformation modulus, and cohesion of the layered HBS are significantly lower than those of the homogeneous model. The deformation modulus of the reservoir is mainly affected by the fine-grained layer without hydrates. When the coarse and fine grains correspond to different mineral components, the two minerals are heterogeneous in terms of their micromechanical parameters, which can further reduce the macroscopic mechanical parameters of the HBS. In addition, the layered distribution of hydrate results in significant anisotropy of the reservoir. This study constitutes a reference regarding the control mechanism of gas hydrate reservoir strength.

Natural gas hydrate (NGH) is considered as a type of clean energy to replace coal and oil. During exploitation, permeability is one of the key parameters controlling production efficiency, reservoir stability, and greenhouse gas sequestration. Limited by experimental and numerical simulation tools, in current research, the directionality of permeability is usually ignored. In this work, a DEM–CFD coupling simulation method is developed to compute the anisotropic permeability. The sedimentary process of reservoir sediments is reconstructed, enabling the acquisition of numerical models that possess pore structures consistent with the actual fabric characteristics. The fluid transport process in various directions can be simulated with a finite element method. Taking the natural gas hydrate reservoir in the Shenhu area of the South China Sea as an example, the proposed method is validated and applied to explore the effect of compaction stress on permeability anisotropy. With the increase in compaction stress, the permeability anisotropy exhibited a rapid initial increase, followed by a sustained stabilization. The primary cause is the rearrangement of sediment particles. The non-spherical particles are driven to align in a predominantly horizontal orientation, thereby enhancing anisotropy. The proposed method provides a tool for the efficient exploitation of hydrate resources.

The permeability of natural gas hydrate (NGH) turbidite reservoirs typically exhibits significant anisotropy, with anisotropy being a crucial basis for evaluating reservoir production. The presence of hydrates, as a crucial constituent of the solid framework, not only impacts the overall permeability but also influences the permeability anisotropy. To investigate the saturation sensitivity of permeability anisotropy, a series of simulations are performed by integrating particle flow and computational fluid dynamics methods to construct the homogeneous and layered numerical samples and compute the evolution of permeability anisotropy. It is shown that the permeability is isotropic for homogeneous sediments and the isotropy remains unchanged regardless of variations in hydrate saturation. The permeability of layered sediments, in contrast, exhibits significant anisotropy due to the presence of dominant channels within the coarse layer. For uniformly distributed hydrates, the more effective blockage in coarse layers results in a reduction in anisotropy. While for preferentially distributed hydrates, the excess blocking of coarse layers makes the dominant channels transfer to the fine layers, the further blocking causes a U-shaped anisotropy–saturation curve characterized by a decrease–increase transformation. During the reservoir production process, the preponderance channels blocked by hydrates will be cleared and the horizontal permeability will significantly increase. As a result, the production efficiency of horizontal wells may exceed expectations. The findings offer a parameter support for production estimation and environmental assessment.

To explore the influence of microfracture development caused by chemical dissolution on the mechanical properties of limestone, this paper presents a new numerical simulation and quantitative analysis method. First, the dissolution rate was determined by the theory of chemical kinetics, and a differential equation that can be solved for results of the fracture evolution process by COMSOL Multiphysics was established to describe the microfracture's expansion. The fractal dimensions of the microfractures were found to have a linear relationship with the damage variables at different time periods through analysis of the simulation results with the fractal geometry method using fracture width as the index, which proves that the evolution of damage has a fractal nature. After that, a damage evolution equation was fitted to predict the deterioration in rock mechanical properties under hydrochemical actions and the predictive uniaxial compressive strength of limestone is seen to be in agreement with experimental test results. The application of the fractal geometry method has important engineering significance as it relates the development of microscopic fractures to changes in the macroscopic mechanical properties and predicts the mechanical properties of the rock under chemical damage.