Численное моделирование и оценка закономерностей развития зон концентрации напряжений в массиве горных пород при влиянии горных работ

УДК 622.8+539.3

Аннотация

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

На основе полученных результатов обеспечивается возможность оценки геомеханического состояния массива горных пород, заблаговременного выявления потенциально опасных зон концентрации напряжений и областей релаксации, а также прогнозирования устойчивости подземных сооружений. Результаты исследования могут быть использованы в качестве научно-практической основы для выбора рациональной последовательности ведения горных работ на рудниках Ангренского рудоуправления, обоснования параметров крепления горных выработок и разработки технических мероприятий, направленных на повышение устойчивости горных выработок и безопасности ведения подземных горных работ.

Abstract

The article presents the results of research aimed at assessing the geomechanical state of the rock mass under mining conditions at the mines operated by the Angren Mining Administration. The relevance of the study is determined by significant changes in the initial stress state of the rock mass resulting from increasing mining depth, growing mining volumes and intensity, as well as the expansion of mined-out areas. Under these conditions, mining operations induce a redistribution of the natural stress field, which may lead to the formation of localized zones of stress concentration and areas of stress relief within the rock mass. These geomechanical processes can cause deformation of the contours of underground structures, rock failure, and a reduction in the stability of mine excavation roofs, thereby adversely affecting the safety of mining operations. The main objective of the study is to establish the patterns of changes in the stress–strain state of the surrounding rock mass resulting from mining operations at the mines of the Angren Mining Administration, to assess the mechanisms of stress redistribution, and to determine the conditions governing the formation of zones of stress concentration and relaxation.

Based on the results obtained, it is possible to assess the geomechanical state of the rock mass, identify potentially hazardous zones of stress concentration and areas of stress relaxation in advance, and predict the stability of underground structures. The results of the study can be used as a scientific and practical basis for selecting a rational sequence of mining operations at the mines of the Angren Mining Administration, substantiating the parameters of mine excavation support, and developing technical measures aimed at improving the stability of mine excavations and ensuring the safety of underground mining operations.

Ключевые слова: численное моделирование; концентрация напряжений; горные выработки; геомеханика; перераспределение напряжений; напряжённое состояние; горным давления; релаксационные зоны; тектонических нарушений.

Keywords: numerical modeling; stress concentration; mine workings; geomechanics; stress redistribution; stress state; rock pressure; relaxation zones; tectonic disturbances.

Введение

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

На рудном поле месторождения «Кызыл-Олма» участки Самарчук являются основными участками ведения горных работ. При отработке промышленных запасов на данных участках деформационные процессы в окружающем массиве горных пород распространяются на значительные расстояния, вплоть до приповерхностной зоны. Кроме того, для обеспечения сохранности подземных сооружений возникает необходимость оставления охранных целиков в определённом количестве и в соответствии с установленным порядком их расположения.

Развитие фронта горных работ на месторождении требует реализации комплекса определённых мероприятий, направленных на обеспечение устойчивости и сохранности подземных сооружений. Одновременное развитие горных работ на различных участках месторождения существенно усложняет процесс обеспечения устойчивости подземных сооружений, поскольку реализация необходимых защитных мероприятий требует значительных финансовых инвестиций и продолжительного времени. Разлом Гушсай, расположенный в центральной части месторождения, имеет мощность до 400–550 м и протяжённость до 850–1100 м. Рудное тело, расположенное в центральной части месторождения. Его мощность изменяется от 30–40 м до 180–200 м, а протяжённость достигает 1000 м.

Материалы и методы

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

Результаты и обсуждение

В связи с этим построение реальных расчётных схем без применения специальных методов расчёта становится практически невозможным. Реализация таких расчётов осуществляется исключительно с использованием компьютерных технологий. Среди указанных методов расчёта в первую очередь следует выделить метод конечных элементов (МКЭ) [2].

В процессе завершения горных работ на подэтажах горизонта +925, Характер распределения тангенциальных напряжений τm в массиве горных пород после проведения горных работ на площади месторождения «Кызыл-Aлма» представлены на рисунке 1.

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

Для количественной оценки характера деформирования и возможного разрушения породного массива выполнено численное геомеханическое моделирование с использованием программного комплекса ROCSCIENCE (Examine 2D). Результаты расчетов показали, что сформировавшаяся вокруг очистного блока зона разрушения имеет сравнительно небольшие размеры и не оказывает существенного влияния на устойчивость окружающего массива при рассматриваемых условиях.

При построении расчетной модели первоначально формируется геометрическая схема исследуемого участка, после чего определяются соответствующие начальные и граничные условия. Особое внимание при моделировании напряжённо-деформированного состояния трещиноватого массива уделяется учёту его природной структурной неоднородности. В зависимости от характера пространственной ориентировки и взаимного расположения трещин породный массив может быть представлен тремя основными структурными типами: массивом с системой протяжённых параллельно ориентированных трещин; массивом с выраженной блочной структурой; массивом с системой трещин, характеризующихся неупорядоченной, или хаотической, пространственной ориентировкой.

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

В рамках метода конечных элементов (МКЭ) определение механических параметров массива осуществляется с использованием усовершенствованного варианта методики, основанной на критерии прочности Хука–Брауна (Hoek–Brown). Данный подход позволяет учитывать влияние структурных особенностей и степени нарушенности массива на его прочностное состояние и деформационные свойства. Для массива, характеризующегося системой параллельно ориентированных трещин, расчёт упругих характеристик выполняется на основе соответствующих зависимостей с дополнительным учётом состояния трещин. В частности, отдельно рассматривается случай открытых (несомкнутых) трещин, поскольку их наличие существенно влияет на эффективные деформационные характеристики и перераспределение напряжений в окружающем массиве.

Модуль упругости массива в вертикальном направлении:

Модуль упругости массива в горизонтальном направлении:

Коэффициент Пуассона:

где: μᵢ – геометрическая характеристика i-й системы трещин, отражающая ориентацию; плотность и особенности пространственного расположения трещин; -угол внутреннего трения сыпучего материала; E и ν – соответственно модуль упругости и коэффициент Пуассона пород, слагающих массив горных пород.

Результаты численного моделирования напряжённо-деформированного состояния массива горных пород по компоненте горизонтальных напряжений после отработки двух маломощных рудных тел, расположенных на горизонте +925 м центрального участка рудника «Кызыл-Алма», представлены на рисунке 2. Полученные результаты характеризуют пространственное перераспределение горизонтальных напряжений в массиве горных пород под влиянием проведённых горных выработок и позволяют оценить изменение напряжённого состояния в окрестности отработанных рудных тел. На основе результатов моделирования установлены области повышенной и пониженной концентрации горизонтальных напряжений, формирующиеся вследствие изменения первоначального равновесного состояния массива при ведении горных работ.

В процессе численного моделирования были приняты следующие механические характеристики рудного тела и вмещающих горных пород, определяющие их напряжённо-деформированное состояние и механическое поведение при воздействии горных работ: объёмный вес горных пород γ = 0,025 МН/м³; модуль Юнга для рудного тела E = 3850 МПа, для вмещающих горных пород E = 6100 МПа; коэффициент Пуассона руды – 0,27; удельный вес горных пород принят равным 0,022 МН/м³.

На рисунке 3 представлен график распределения тангенциальных напряжений τm вдоль контура рудного тела 10 Б. Согласно результатам моделирования, значения тангенциальных напряжений по контуру рудного тела изменяются в диапазоне от 3 МПа до 44 МПа. Такой диапазон изменения напряжений характеризует неоднородность их распределения по контуру выработанного пространства, сформированного в результате отработки рудного тела.

По горизонтальной оси OX графика отложена длина контура выработанного пространства, образовавшегося в результате отработки 10 Б рудного тела. По мере перемещения вдоль контура наблюдаются изменения величины тангенциальных напряжений τ_m, обусловленные перераспределением напряжений в массиве горных пород после удаления рудного тела. Полученная зависимость позволяет оценить характер распределения тангенциальных напряжений по контуру выработанного пространства и выделить участки с повышенными значениями напряжений.

На данном графике представлена динамика неравномерного распределения тангенциальных напряжений τ_m в массиве горных пород по периметру рудного тела, на основании чего можно сделать следующие выводы: наблюдается зона высоких напряжений в начальном участке; затем зона релаксации; зона сильной локальной концентрации в интервале 105–112 м; после нее резкое снижение напряжений на участке 113–128 м и, наконец, зоны повторно возрастающего опорного давления в заключительной части. Подобное состояние обычно объясняется взаимным влиянием фронта очистных работ, тектонических нарушений и зон опорного давления

Заключение

Проведённое моделирование напряжённого состояния массива горных пород в зоне формирующегося выработанного пространства при разработке месторождения на горизонте +925 м центрального участка рудника «Кызыл-Алма», а также в пределах транспортного штрека этого горизонта позволило установить характерные особенности перераспределения горизонтальных напряжений в массиве:

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

2. В угловой части, примыкающей к висячему боку рудного тела выработанного пространства, фиксируются максимальные значения тангенциальных напряжений (рисунок 2). Данный участок является локальным очагом концентрации напряжений и оценивается как зона с высокой вероятностью развития потенциальных деформационных процессов.

3. Согласно результатам моделирования, величина растягивающих и сжимающих напряжений, формирующихся вблизи боковых стенок камеры, ниже соответствующих пределов прочности горных пород, что свидетельствует об относительно ограниченном риске резкого нарушения устойчивости на данных участках.

Использование установленных закономерностей формирования зон повышенных напряжений в дальнейшем позволит научно обоснованно разрабатывать геотехнологические схемы и оптимизировать параметры ведения горных работ при освоении участков, расположенных в пределах зон концентрации динамических проявлений.

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

На основе полученных закономерностей могут быть обоснованы рациональные параметры систем разработки, последовательность и направление ведения очистных работ, размеры очистных пространств, параметры оставляемых целиков и мероприятия по управлению горным давлением. Кроме того, появляется возможность корректировки параметров крепления, схем поддержания горных выработок и мер по обеспечению устойчивости массива в условиях повышенной концентрации напряжений.

Применение данного подхода позволит повысить безопасность ведения горных работ, снизить вероятность возникновения динамических проявлений горного давления, обеспечить устойчивое состояние горных выработок и повысить эффективность освоения месторождения в сложных геомеханических условиях.

Список литературы

References