Представлены результаты моделирования полей напряжений и анализа прочности породного массива на участке Енисейский (Красноярский край), выбранного для строительства подземной исследовательской лаборатории. Обоснованы варианты граничных условий нагружения моделей по ее границам, результаты моделирования распределения компонент тензора напряжений для четырех вариантов нагружения и оценка стабильности массива с использованием известных моделей прочностного анализа Хука – Брауна, Кулона – Мора, Мизеса и др. Определены закономерности распределения полей напряжений в массиве и отличия, связанные с тектоническими условиями района. Установлено, что локализация зон концентрации интенсивности напряжений зависит от соотношения компонент главных действующих напряжений. Ориентация сжатия в субмеридиональном направлении приводит к росту интенсивности напряжений на 10-15 % относительно иных вариантов моделирования. Зоны аномальных значений интенсивности напряжений расположены внутри блоков, а также в лежачих крыльях тектонических разломов. Модели характеризуются высокими значениями потенциальной энергии формоизменения в зонах разломов (как участков массива, наиболее легко поддающихся деформированию) и их пересечений. Трехмерное моделирование позволяет обнаружить эффекты, которые слабо проявляются в моделях плоского деформированного состояния. Результаты геомеханического моделирования необходимы при планировании экспериментов в подземной исследовательской лаборатории для уточнения изоляционных свойств породного массива при захоронении высокоактивных радиоактивных отходов. Методические подходы трехмерного моделирования применяются геомеханическими и геотехническими службами промышленных предприятий и других опасных производственных объектов (подземные хранилища газа, месторождения полезных ископаемых и др.).