Розглядається нестаціонарна тривимірна модель потоку грунтових вод та транспорту забруднень у насичено-ненасичених середовищах. Основою моделі є зв"язані рівняння вологопереносу і масопереносу в ненасиченій зоні та зв"язані рівняння Бусінеска і масопереносу в грунтових водах. Для розв"язування відповідних крайових задач використовується неявний метод скінченних різниць на нерівномірних сітках. Створено пакет програм, що складається з розрахункового модуля та модуля візуалізації на основі ПС MapInfo. Наведено результати моделювання для будівельного майданчика об"єкту "Вектор" (30-км зона Чорнобильської АЕС).