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