The work is devoted to modeling the phase transformations of gas hydrate inclusions in porous media. Studying these problems makes it possible to elaborate various technologies for the development of...

The purpose of this paper is to present a one-dimensional model that simulates the thermo-physical processes for methane hydrate decomposition in porous media. The mathematical model consists of equations for the conservation of energy, gas, and liquid as well as the thermodynamic equilibrium equation for temperature and pressure $\left(P-T\right)$ in the hydrate stability region. The developed model is solved numerically by using the implicit finite difference technique on the grid system, which correctly describes the appearance of phase, latency, and boundary conditions. The Newton–Raphson method was employed to solve a system of nonlinear algebraic equations after defining and preparing the Jacobean matrix. Additionally, the proposed model describes the decomposition of methane hydrate by thermal catalysis of the components that make up the medium through multiple phases in porous media. In addition, the effect of thermodynamic processes during the hydrate decomposition on the pore saturation rate with hydrates a7nd water during different time periods was studied in a one-dimensional model. Finally, in a one-dimensional model over various time intervals, $t=1,\hspace{0.17em}10,\hspace{0.17em}50\hspace{0.17em}\mathrm{s}$, the pressure and temperature distributions during the decomposition of methane hydrates are introduced and investigated. The obtained results include more accurate solutions and are consistent with previous models based on the analysis of simulations and system stability.

The thermodynamically equilibrium model for splitting by the physical processes of a two-component three-phase filtration fluid dynamics with gas hydrate inclusions is considered in the paper, for which a family of two-layer completely conservative difference schemes of the support operators method with time weights profiled in space is constructed. On the irregular grids of the theory of the support-operators method applied to the specifics of the processes of transfer of saturations and internal energies of water and gas in a medium with gas hydrate inclusions, methods of directwind approximation of these processes are considered. These approximations preserve the continual properties of divergence–gradient operations in their difference form and are related to the velocity field providing saturations transfer and internal energies of fluids. Fluid dynamics with gas hydrate inclusions are also calculated on the basis of the proposed approach, in particular, in areas of severe pressure depression in the collector space.

The work is devoted to numerical modeling of the processes of heat and mass transfer of a two-phase fluid in the environment of a production well, which is necessary for monitoring the development of fractured-porous reservoirs. This work proposes an efficient approach to constructing a solution to the problem. To solve the problem, a model of the “double medium” type is used, where the pore part of the reservoir is considered as the first medium, and the system of natural fractures is considered as the second medium. For the resulting mathematical model, the difference schemes with time weights are constructed based on the algorithm of splitting by physical processes, which ensure the correctness and consistency of fluxes in the fracture system and the pore reservoir. In the numerical solution, the approximations of differential operators obtained in the framework of the finite difference method are used. For the parallel implementation of the developed numerical approach, the domain decomposition method and the matrix sweep algorithm are chosen. The program implementation is made using the MPI standard. Computational experiments are carried out, the results of which confirm the effectiveness of the developed numerical algorithm and its parallel implementation. In numerical experiments, the distributions of pressure and temperature near an operating production well are obtained, on the basis of which it is possible to adjust the operation of wells in order to increase production.

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.