Program package for two-dimensional modeling offluidodynamic processes in hydrate-containing environment of sedimentary basins

Currently, the investigation of gas hydrate deposits is very relevant. The study of sedimentary basins containing gas hydrates is complicated by the presence of complex geometric regions in them (interlayers, shear and low-permeability zones, etc.). Areal and profile calculations have spatial stratigraphic unstructured areas, their arbitrary location and variable effective thickness of layers. This work focuses on the study of the problem of two-dimensional filtration in a zone of three-phase hydrate equilibrium (gas, hydrate, water), taking into account the possible low permeability of rocks. To solve this problem, an original filtration model with a non-classical form of the law of motion is used. Based on the support operator method, the effective computational algorithms that distinguish the hyperbolic (for the transfer of water and hydrate saturation) and piezoconductive-dissipative parts of the problem have been developed. Their program implementation on grids of irregular structure, which make it possible to model the processes of dissociation of gas hydrates in sedimentary basins of complex geometric and lithological structure were realized. Numerical calculations for a model domain with parameters corresponding to a real sedimentary basin in the South China Sea were realized.