Numerical method for calculating heat and mass transfer of two-phase fluid in fractured-porous reservoir

The work proposes an approach to numerical modeling based on splitting by physical processes for a non-isothermal problem of filtration in a fractured-porous medium. The task is complicated by the presence of a two-phase fluid and dual porosity of the reservoir. The system of equations defining the model is complex and is described by a system of strongly non-linear partial differential equations. The use of the splitting method according to physical processes makes it possible to simplify the solution algorithm while maintaining the equivalence to the conservative difference approximation of the original equations and ensuring the stability of the problem solution. In the numerical solution, the approximations of differential operators obtained in the framework of the finite difference method are used. The implementation of the numerical algorithm is based on the matrix sweep method. To test the method, a series of computational experiments were carried out. Calculations have shown that the developed methodology is correct and allows one to simulate the dynamic operating conditions of wells.