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