Natural gas hydrates (NGH) are considered as a future clean energy in the era of carbon neutrality. The seepage behaviour of hydrate-bearing sediment (HBS) and its controlling factors are significant for developing effective production strategies. In this study, we aim to elucidate how the heterogeneous spatial distribution of clayey-silty sediment particles, particularly the non-uniform particle size distribution (PSD) affects the permeability of HBS. An algorithm was developed to generate the heterogeneous spatial distribution of sediment particles based on the PSD of hydrate-bearing cores extracted from Japan Nankai Trough (D50 = 122 μm). A series simulation of fluid flow in porous media based on CFD was conducted for the HBS with varying hydrate saturations (SH). The watershed segmentation algorithm was employed for the identification of pore structure, including pore space and pore throat with a detailed analysis on the pore structure evolution with SH. Host sediments with higher uniformity result in initial high permeability but the permeability reduction with increasing SH is more significant compared with that of lower uniformity. Pore throat width reduced to less than 5 μm with SH above 50% compared with its original width of ∼30 μm. Tortuosity exhibits a linear increase with SH, while permeability decreases exponentially with increasing tortuosity. We propose an improved permeability model (KCpro) based on the classical Kozeny-Carman model accounting for the change in specific surface area for non-uniform PSD in HBS with good agreement. This study confirmed the applicability of CFD modelling as a promising tool to study the seepage behaviour of HBS complementary to classical permeability experimental tests. Results of this study provide fundamental understanding on how permeability evolves in HBS with significant PSD variations in natural gas hydrate reservoirs in nature.
Read full abstract