An Analytical Solution and Numerical Modeling Study of Gas Hydrate Saturation Effects on Porosity and Permeability of Porous Media

There are two major types of natural gas hydrate distributions in porous media: pore filling and contact cementing. The difference between these two distribution types is related to hydrate nucleation and growth processes. In the pore filling distribution, hydrate nucleates from a gas-dissolved aqueous phase at the grain boundary and grows away from grain contacts and surfaces into the pore space. In the contact cementing distribution, hydrate nucleates and grows at the gas-water interface and at intergranular contacts. Previous attempts to correlate changes on porosity and permeability during hydrate formation/dissociation were based on the length difference between the pore body and pore throat, and only considered contact cementing hydrate distribution. This work consists of a study of mathematical models of permeability and porosity as a function of gas hydrate saturation during formation and dissociation of gas hydrates in porous media. In this work, first we derive the permeability equation for the pore filling hydrate deposition as a function of hydrate saturation. Then, a more comprehensive model considering both types of gas hydrate deposition is developed to represent changes in permeability and porosity during hydrate formation and dissociation. This resulted in a model that combines pore filling and contact cementing deposition types in the same reservoir. Finally, the TOUGH+Hydrate numerical reservoir simulator was modified to include these models to analyze the response of production and saturation during a depressurization process, considering different combinations of pore filling and contact cementing hydrate distributions. The empirical exponent used in the permeability adjustment factor model influences both production profile and saturation results. This empirical factor describes the permeability dependence to changes in porosity caused by solid phase formation in the porous medium. The use of the permeability exponent decreases the permeability of the system for a given hydrate saturation, which affects the hydraulic performance of the system. However, from published experimental work, there is only a rough estimation of this permeability exponent. This factor could be represented with an empirical equation if more laboratory and field data becomes available.