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.

This work considers the isothermal process of incompressible viscous fluid filtration in an oil-saturated, fractured-porous reservoir. A study of the pressure and water saturation distribution process is carried out for a case in which a production well is put into operation. For this problem, i.e., a mathematical model in a two-dimensional formulation, a numerical method and a parallel algorithm are proposed. The mathematical model of two-phase filtration is written in accordance with the classical laws of continuum mechanics and Darcy’s law and also includes a function of fluid exchange between low-permeability pores and high-permeability natural fractures within the framework of the Warren–Root model. The numerical solution is based on the finite-difference method and a splitting scheme of physical processes and spatial coordinates. For a split system with respect to piezoconductivity, an implicit finite-difference scheme with fixed saturations is constructed, and with respect to saturation transfer, explicit and implicit difference schemes are constructed. For parallel implementation of the developed numerical approach, a method based on geometric parallelism is selected. Testing of the developed method is performed using the example of calculating liquid mass transfer for a wide range of parameters. To verify the model, the obtained calculated pressure curves are compared with field data recorded by a deep-well measuring device. The results allow for estimation of the distribution of reservoir pressure and water saturation depending on the permeability of the fracture set and the pore part. The obtained results allow for monitoring of well operations, reducing unexpected accident risks and optimizing the development system in order to increase oil production in fractured-porous reservoirs. Computational experiments confirm the efficiency of the developed numerical algorithm and its parallel implementation.

The process of two-phase filtration in a carbonate formation of fractured-pore type is considered. A mathematical model in a spatially two-dimensional formulation is proposed, a numerical method for solution and a parallel algorithm for its implementation are developed. The mathematical model is based on the Buckley-Leverett approach. The reservoir takes into account the exchange of fluids between low-permeability pores and natural fracturing, specified within the framework of the dual porosity model. The numerical algorithm is based on the use of the finite difference method and the splitting scheme by physical processes. To speed up calculations, a parallel algorithm based on two-dimensional domain decomposition is used. Numerical experiments were carried out, which showed that the developed algorithm is highly efficient and allows one to calculate the necessary characteristics of the modeled physical process.