При разработке угольных месторождений подземным способом активизируются опасные геодинамические явления и процессы на локальном и региональном уровнях. Для определения риска их развития необходимо выполнять оценку напряженного состояния района размещения угольных шахт. В статье представлена методика численного моделирования полей напряжений в районах разработки угольных месторождений подземным способом. Было смоделировано напряженно-деформированное состояние фрагмента геологической среды в районе г. Ленинска-Кузнецкого, где эксплуатируются несколько крупных угольных шахт (им. Кирова, им. Рубана, им. 7-го ноября и др.). Расчеты напряжений были выполнены методом конечных элементов в программном продукте «COMSOL». В качестве расчетной схемы была использована карта разломной тектоники района центрального Кузбасса. Граничные условия и физико-механические свойства модели задавались на основе геолого-геофизических исследований территории. В результате были получены и проанализированы компоненты тензора напряжений в ситуации до и после выработки угольного месторождения.