Mathematical Modeling of Gas Hydrates Dissociation in Porous Media with Water-Ice Phase Transformations Using Differential Constrains

2D numerical modeling algorithms of multi-component, multi-phase filtration processes of mass transfer in frost-susceptible rocks using nonlinear partial differential equations are a valuable tool for problems of subsurface hydrodynamics considering the presence of free gas, free water, gas hydrates...

Full description

Bibliographic Details
Published in:Mathematics
Main Authors: Natalia Alekseeva, Viktoriia Podryga, Parvin Rahimly, Richard Coffin, Ingo Pecher
Format: Article in Journal/Newspaper
Language:English
Published: MDPI AG 2022
Subjects:
Ice
Online Access:https://doi.org/10.3390/math10193470
https://doaj.org/article/543b5e9c23ed44138cd99418d3964578
Description
Summary:2D numerical modeling algorithms of multi-component, multi-phase filtration processes of mass transfer in frost-susceptible rocks using nonlinear partial differential equations are a valuable tool for problems of subsurface hydrodynamics considering the presence of free gas, free water, gas hydrates, ice formation and phase transitions. In this work, a previously developed one-dimensional numerical modeling approach is modified and 2D algorithms are formulated through means of the support-operators method (SOM) and presented for the entire area of the process extension. The SOM is used to generalize the method of finite difference for spatially irregular grids case. The approach is useful for objects where a lithological heterogeneity of rocks has a big influence on formation and accumulation of gas hydrates and therefore it allows to achieve a sufficiently good spatial approximation for numerical modeling of objects related to gas hydrates dissociation in porous media. The modeling approach presented here consistently applies the method of physical process splitting which allows to split the system into dissipative equation and hyperbolic unit. The governing variables were determined in flow areas of the hydrate equilibrium zone by applying the Gibbs phase rule. The problem of interaction of a vertical fault and horizontal formation containing gas hydrates was investigated and test calculations were done for understanding of influence of thermal effect of the fault on the formation fluid dynamic.