Разработаны вычислительные алгоритмы для решения задач переноса. Построены численные схемы для одного уравнения переноса с запаздыванием для двумерных задач, где учитывают формирование концентраций химических веществ флюида и горных пород. Проведенные сравнение исходных численных вычислений соответствует с аналитическими расчетами. Проведен анализ физической и численной устойчивости решений систем уравнений переноса в одномерном и двумерном случаях подземного выщелачивания с реакционным запаздыванием