Numerical studies of gas hydrate systems: sensitivity to porosity-permeability models

The formation of sub-seafloor gas hydrates in marine environments can be described as a coupled transport and thermodynamic process inside a host sediment matrix undergoing structural evolution. The transport processes are driven by the sedimentary load and induced overpressure gradients, controlled by sediment permeability. In order to accurately model the resulting fluid flow profile, the decrease of sediment permeability during hydrate precipitation has to be taken into account, which affects both the transport of solutes and sediment compaction. In this paper, we investigate how total hydrate abundance is affected by regions of low permeability which deflect the flow field in their vicinity. For this purpose, a two-dimensional numerical hydrate system model was set up which permits to quantify this effect in scenarios where changes in water depth cause lateral variations of the thickness of the hydrate stability field, as well as of hydrate saturation and sediment permeability. The microscopic structure of gas hydrate crystals in the host sediment matrix defines the evolution of the permeability reduction during hydrate formation. Grain-coating precipitates have a stronger tendency to clog flow paths through pore throats than do pore-filling precipitates. Our results clearly show that these pore-scale processes affect the large-scale flow field and hydrate abundance. The sensitivity depends on the model geometry and, for a 5° slope of the seafloor, 4.1% relative difference is predicted for the hydrate saturation according to different porosity-permeability relationships.