Numerical Determination of Anisotropic Permeability for Unconsolidated Hydrate Reservoir: A DEM–CFD Coupling Method

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.