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.