Актуальность исследования обусловлена необходимостью обеспечения безопасности на Соколовском подземном руднике, а также повышения полноты и качества выемки руды в сложных горно-геологических условиях – при наличии покрывающей толщи обводненных рыхлых отложений. Цель: определить оптимальные параметры устойчивых камер и целиков для предотвращения обрушения горных пород и проникновения песчано-глинистых масс в подземные выработки. Объект: массив горных пород на Северном участке Соколовского месторождения. Методы. Многовариантное численное моделирование методом конечных элементов в двумерной постановке выполнялось в программе Rocscience RS2, что позволило учесть большое количество факторов, влияющих на состояние массива. При расчетах приняты во внимание не только физико-механические свойства пород и действующие в массиве напряжения, но и структурные характеристики массива, а также степень техногенного воздействия. В качестве модели поведения массива использовалась модель Хоека–Брауна, отличительным преимуществом которой является ее нелинейность. Так как реальный массив характеризуется неоднородностью прочностных свойств и напряженно-деформированного состояния, был выполнен предусмотренный в программном комплексе Rocscience вероятностный анализ, учитывающий изменение граничных условий в установленных интервалах для создаваемой модели. В результате численного моделирования определено изменение напряженно-деформированного состояния массива, на основании чего вычислены коэффициент устойчивости окружающего массива, величины и направления его смещений, а также установлены зоны распределения вероятности разрушения массива в заданных интервалах граничных условий. Определены параметры камер и целиков, обеспечивающие общую устойчивость, минимизирующие возможность прорыва водогрязевых масс в выработки. Выводы. Сохранение устойчивости кровли подземных камер при отработке очистных рудных блоков на шахте «Соколовская» позволит снизить риск проникновения песчано-глинистых масс в подземные выработки, обеспечивая безопасное и эффективное ведение горных работ, а также повысить полноту и качество выемки руды. Выполненное методом конечных элементов многовариантное численное моделирование влияния очистных выработок на массив позволило определить оптимальные геометрические параметры камер и целиков, обеспечивающие устойчивость окружающего массива. Применение программного комплекса Rocscience RS2 в результате учета большого количества факторов, влияющих на состояние массива, значительно повышает надежность оценки устойчивости массива.