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