Abstract
<jats:p>This article presents a mathematical model and algorithm for numerically solving a problem for studying the hydrodynamic process of in-situ leaching in a heterogeneous porous medium, taking into account mass transfer kinetics, the law of layer deformation dependent on the elastic coefficient, changes in porosity depending on pressure, and the permeability of the ore reservoir for the purpose of developing ore deposits. The developed mathematical apparatus allows for a comprehensive study of the properties of the ore seam, optimization of the location of production and injection wells, determination of changes in porosity under pressure, and analysis and provision of a factor ensuring the protection of groundwater from pollution. The model under consideration is described using a system of second-order parabolic partial differential equations with initial, boundary, and internal conditions. Since the analytical solution to this problem was too complex or impossible for numerical integration, a second-order finite-difference scheme was used.</jats:p>