Computational modeling of methane hydrate dissociation in a sandstone core

Hydrate dissociation in a porous sandstone core was studied using a computer modeling approach. It was assumed that the hydrate was dispersed in the pores of the core. Using FLUENTTM code, an axisymmetric model of the core was developed and solved for multiphase flows during the hydrate dissociation. The core model contained three separate phases: methane hydrate, methane gas, and liquid water. At the start of simulation, the valve at one end of the core was opened exposing the core to low pressure; hydrate began to dissociate and methane gas and water began to flow. The depressurization was controlled by adjusting the pressure of the outlet valve. A comprehensive Users’ Defined Subroutine (UDS) for analysis of hydrate dissociation process into the FLUENT code was developed. The new UDS uses the kinetic model introduced by Kim et al. [Kim. H.C., Bishnoi, P.R., Heidemann, R.A., Rizvi, S.S.H., 1987. Kinetics of methane hydrate decomposition. Chemical Engineering Science 42, 1645–1653.] and can model multiple zones dissociation and multiphase flows. Variations of relative permeability of the core were included using Corey's model. The new model allows for variation of the porosity with hydrate saturation. For different core temperatures and various outlet valve pressures, the spatial and temporal variations of temperature, pressure, and flow fields in the core were simulated. The time evolutions of methane gas and water flow rate at the outlet were also evaluated. It was shown that the rate of hydrate dissociation in a core was a sensitive function of surrounding environment temperature, outlet pressure condition, and permeability.