Pore-Scale Modellings on the Impacts of Hydrate Distribution Morphology on Gas and Water Transport in Hydrate-Bearing Sediments

Gas and water transport behavior, which is controlled by the pore characteristics and capillarity in hydrate-bearing sediments (HBS), is one of the key factors affecting the gas production. Hydrate distribution morphology (HDM) can significantly influence the pore structures of HBS, affecting the relative permeabilities of gas and water. To elucidate the impacts of HDM in microscopic scale, a phase-field lattice Boltzmann (LB) model is developed to describe the gas and water transport in HBS.To simulate the transport of immiscible fluids, which exist obvious density and viscosity contrasts, a phase-field LB model with the conservative form of interface-tracking equation is developed to suppress the spurious currents at phase interfaces. To describe the fluid-solid interactions, the bounce-back condition is applied for both solid phases (hydrate and grains) to achieve the non-slip condition and the wettability condition is applied for grains and hydrate to describe the wettability behavior. To improve the numerical stability, the multi-relaxation-time (MRT) collision operator is applied and the discretization schemes with 8th order accuracy for the gradient operator are selected. In this work, we first validated our model by applying several benchmark cases aiming at fluids with obvious density contrasts such as the layered Couette/Poiseuille flows, Rayleigh–Taylor instability. Then the synthetic geometries of the pore-filling and grain-coating HBS with several hydrate saturation (Shyd) were constructed by guaranteeing the same extent of connectivity. Then the steady-state relative permeability measurement and drainage capillary pressure measurement processes were simulated by the LB model for two HDM cases under several Shyd. The results showed that in the hydrophilic HBS, the relative permeability of gas in the pore-filling case is obviously larger than that in the grain-coating case at the same Shyd, and larger capillary pressure can be obtained in the pore-filling case. In addition, as the Shyd increased, it would notably enhance these differences of fluids relative permeability and capillary pressure between two HDM cases. Because the HDM can not only influence the pore space structures but also the wettability of the porous medium by creating solid surfaces of varying wettability, the distribution and transport of fluid phases in different HDM cases can be obviously affected. The phase-filed LB model applied in this study is capable to handle and suppress the spurious currents at phase interfaces, ensuring a satisfactory numerical stability and accuracy. Thus, the real density and viscosity contrasts between the water and gas under the in-situ thermodynamic conditions can be considered in the simulation. The impacts of HDM on the gas and water transport were quantitively analyzed by simulating multiphase flow processes in HBS.