Изложены алгоритмы прогноза напряженно-деформированного состояния неоднородного слоистого физически нелинейного породного массива. Моделирование базируется на применении комплекса вычислительных методов: вариационного, дискретного продолжения по числовому параметру, квазилинеаризации нелинейных краевых задач, конечных разностей, матричной прогонки и общего итерационного процесса.