CC BY
13
УДК 622.831
A.А.КОЗЫРЕВ д-р техн. наук, зам. директора по научной работе, [email protected] И.Э.СЕМЕНОВА, канд. техн. наук, старший научный сотрудник
B.В.РЫБИН, канд. техн. наук, старший научный сотрудник И.М.АВЕТИСЯН. м.иьПиий научный сотрудник
Горный институт Кольского научного центра РАН. Апатиты
A.A.KOZYREV, Dr. in eng. sc., deputy director on research work, innas(a>goi.kolasc.net.ru I.E.SEMIONOVA PhD in eng. sc., senior research assistant V.V.RYBIN PhD in eng. sc., senior research assistant I.M.A VETISYAN, junior research
Mining institute of the Kola research center o/RAS, Apatity
ОЦЕНКА ЧИСЛЕННЫМ МЕТОДОМ ПОЛЯ НАПРЯЖЕНИЙ И ДЕФОРМАЦИЙ В ОКРЕСТНОСТИ КРУПНОЙ КАРЬЕРНОЙ ВЫЕМКИ ПРИ ПОСЛЕДОВАТЕЛЬНОМ УЧЕТЕ ГЕОМЕХАНИЧЕСКИХ, ГЕОЛОГИЧЕСКИХ И ГОРНО-ТЕХНИЧЕСКИХ ФАКТОРОВ
Представлены результаты моделирования напряженного состояния массива пород методом конечных элементов в окрестности крупной карьерной выемки. Модель позволяет учитывать влияние крупных разломов. Граничные условия заданы, исходя из результатов измерений параметров напряжений in situ. Определено положение опасных зон в приборто-вом массиве.
Ключевые словах математическое моделирование, поля напряжений и деформаций, учет геологических факторов, определение опасных зон.
NUMERICAL ASSESSMENT OF STRESS-STRAIN FIELD IN VICINITY OF LARGE OPEN PIT MINE TAKING INTO ACCOUNT GEOMECHANICAL, GEOLOGICAL AND MINING FACTORS
The results of modeling of rock mass stress state in the vicinity of large open pit mine using finite-elements method are represented. The model allows to take into account the influence of large-scale fractures. The initiation of boundary conditions is based on the results of determination of stress parameters in situ. The locations of risk areas near the slope are defined.
Key words', mathematical modeling, stress-strain fields, account of geological factors, determination of hazardous zones.
Одним из эффективных путей минимизации затрат на разработку месторождений полезных ископаемых открытым способом является увеличение угла наклона бортов карьеров. Однако в этом случае остро встает проблема управления геомеханическим состоянием горно-рудной прпродно-техничес-кой системы, которое должно сочетать прогноз напряжений и деформаций в окрестности карьерной выемки и разработку профилактических мероприятий, снижающих риск
потери устойчивости горно-технических конструкций.
Для прогнозирования напряженно-деформированного состояния (НДС) в отечественной и зарубежной практике успешно применяются численные модели, учитывающие геологические и горно-технические факторы и дающие возможность оценки геомеханического состояния массива пород при различных вариантах развития работ в пространстве и времени [1,2, 4-5].
Параметры напряженного состояния массива пород карьера рудника «Железный»
Абсолютная отметка измерительной станции, время проведения измерений, место измерений Яср, м Результаты измерений методом разгрузки Результаты расчета исходя только из собственного веса вышележащих пород
МПа СТгаш, МПа Джюх, град а41 , МПа о1* , МПа ^гор 5
+ 10 м, 2008, северный борт 247 27 8 -17 7,4 2,4
+ 10 м, 2007, северо-восточный борт 247 20 3 -24° 7,4 2,4
+ 25 м, 2007, северо-восточный борг 232 24 14 -15° 7 2,3
+ 40 м, 2005, северо-восточный борт 217 13 4 18* 6,5 2,2
+ 94 м, 2000, северо-восточный борг 163 20 4 18° 5 1,6
+ 94 м, 2000, северо-восточный борг 163 22 9 -14° 5 1,6
+ 118 м, 2005, северо-западный борт 90 5 -1 6° 2,7 0,9
+ 142 м,2001, северо-западный борт 65 8 -2 43° 2 0,6
+ 154 м, 2005, юго-восточный борт 52 11 3 88* 1.6 0,5
+ 166 м, 2003, юго-восточный борт 40 3 -1 13° U 0,4
Условные обозначения: Нч, - средняя глубина положения измерительной станции от первоначального рельефа, м; с„ш - максимальная компонента напряжений, МПа; гт.щц - минимальная компонента напряжений, МПа; -угол наклона ст|Ш>: к горизонту', град.; о4"' -уН - вертикальная компонента. МПа; ст':р - ХуН - горизонтальная
ь верт гор
компонента, МПа.
Рис.1. Схемы задания граничных условий д ля расчетных вариантов: а - гравитационная постановка; б- с учетом тектонических сил < в - с учетом тектонических сил Рь >
Объектом исследований является карь-ер рудника «Железный» ОАО «Ковдор-ский ГОК», проектная глубина которого составляет 1000 м, а угол наклона борта 60°. Натурными исследованиями выявлено действие в массиве пород, включающем месторождение, повышенных горизонтальных напряжений (см. таблицу).
Для проведения модельных исследований НДС был выбран метод конечных элементов (МКЭ) в упругой постановке, что позволило: задать адекватные исследуемому объекту граничные условия по напряженному состоянию; учесть неоднородности различного иерархического уровня и с достаточной детальностью промоделировать геометрию очистных пространств. Расчеты выполнены в среде программного комплекса Sigma3D, разработанного в Горном институте КНЦ РАН [3].
Границы модели выбраны в соответствии с принципом Сен-Венана, согласно которому зона влияния особенности распространяется на расстояние не более трех ее характерных размеров (в нашем случае радиусов карьера по длинной R « 1150 м и по короткой оси г« 830 м, а также его максимальной проектной глубины Нтз> = 1000 м). В итоге размер области моделирования -8000 х 7500 х 3050 м3.
Рис.2. Коэффициенты концентрации стпш в окрестности карьерной выемки с глубиной 1000 м в вертикальных сечениях: а, б - при действии максимальной тектонической нагрузки ^ по короткой оси карьера;
в, г - по длинной оси карьера
С использованием разработанной модели рассчитано несколько вариантов задачи распределения параметров НДС в массиве, которые можно разделить на две группы по граничным условиям: чисто гравитационное напряженное состояние (а\=уН, 02 — 03 = = h/H, где у, Н и X - соответственно объемный вес, глубина и коэффициент бокового распора) и гравитационно-тектоническое напряженное состояние. Причем варианты с учетом тектонических сил были просчитаны как при условии действия их максимальной составляющей по короткой оси карьерной выемки (совпадающей с направлением оси У модели), так и по ее длинной оси (совпадающей с направлением оси X модели). Схемы задания граничных условий представлены на рис. 1.
Для каждого вида граничных условий рассчитывали НДС: для исходного состояния массива пород; при выемке карьерных запасов, соответствующей его реальной глубине, до абсолютной высотной отметки +10 м (глубина карьера 340 м); при проектной глубине карьерной выемки до отметки -650 м (глубина карьера 1000 м).
Заданная максимальная горизонтальная пригрузка составила на уровне дневной поверхности 1 МПа, а на уровне дна модели -65 МПа. Коэффициент бокового отпора с учетом действия тектонической нагрузки Л"б = 0,5. Учитывая общие закономерности распределения тектонических напряжений с глубиной, предположили, что на глубинах более 400-500 м зависимость отах =/(Я) будет нелинейной.
Приведем результаты расчета НДС в окрестности проектной карьерной выемки с учетом действия тектоники. Отмечено образование зон концентрации напряжений у сопряжений дна и бортов карьера. Причем при действии максимальных горизонтальных напряжений по короткой оси карьера значения коэффициентов концентрации (АГатах) достигают 6,2. При действии максимальных горизонтальных напряжений по длинной оси карьера Котах < 2,2. Уровень напряжений в окрестности верхних и средних участков соответствующих бортов, перпендикулярных заданным горизонтальным напряжениям, ниже исходного (рис.2, а, г; 3, а, в).
Рис.3. Коэффициенты концентрации ст^ в окрестности карьерной выемки с глубиной 1000 м в горизонтальных сечениях: а, б - при действии максимальной тектонической нагрузки Р\ по короткой оси карьера; в, г - по длинной оси карьера; а, в- сечение по гор. -20 м: б, г - сечение под дном карьера, гор. -660 м
-500 -
Рис.4. Коэффициенты концентрации сттах в массиве пород, включающем ослабленную зону (выделена штриховкой), в вертикальных сечениях при глубине карьера 340 м
В бортах карьера, согласных с направлением действия формируются зоны концентраций напряжений: при действии по короткой оси карьера (рис.2, б и 3,а) А"атах находится в пределах 1,2-2,6; при действии ,РТ по длинной оси карьера (рис.2, б и 3, в) Котях = 1,1-1,5. Причем более высокий уровень концентрации сжимающих напряжений отмечен в окрестности части борта с вертикальными уступами.
Следует заметить, что существует предположение, подкрепляемое результатами исследований искривления стволов длинных геолого-разведочных скважин, об ориентации максимальной компоненты главных напряжений в изучаемом массиве пород в направлении по длинной оси карьера, что является более благоприятным с точки зрения его устойчивости.
Следующий вариант был просчитан с учетом мощной ослабленной зоны в юго-восточном борту карьера при текущей геометрии карьерной выемки с глубиной 340 м и заданной тектонической пригрузкой по длинной оси карьера Модуль упругости элементов ослабленной зоны задавали на порядок ниже модуля упругости вмещающего массива. В этом случае образуется зона концентрации сттах по участку борта с промоделированной зоной (рис.4), особенно на сопряжении ослабленной зоны с карьерной выемкой и во вмещающем массиве ниже по борту. Значения коэффициентов концентрации А"атах в этом случае достигают 2,6. В аналогичном варианте без учета ослабленной зоны этот участок борта был разгружен, а локальная область концентраций наблюдалась лишь на сопряжении борта с дном карьерной выемки.
Анализ распределения напряжений под дном карьерной выемки показывает увеличение как размеров области концентраций напряжений, так и абсолютных значений напряжений в массиве вмещающих пород в непосредственной близости от ослабленной зоны
В результате выполненной работы:
1. Изучены особенности распределения напряжений в окрестности глубокой карьерной выемки, находящейся в массиве по-
род с гравитационно-тектоническим типом НДС, которые связаны прежде всего с тем, что борт карьера пересекает несколько зон с различным соотношением гравитационных и тектонических сил.
2. Определено местоположение зон концентраций сжимающих напряжений в зависимости от направления действующих в массиве пород тектонических сил при наличии и при отсутствии ослабленной зоны, расположенной в пределах юго-восточного борта карьера. Выявленные области концентраций максимальной компоненты напряжений ата>; позволяют сделать вывод о меньшем уровне действующих напряжений при направлении ,ГТ по длинной оси и о большем уровне действующих напряжений при направлении /ч по короткой оси выемки.
Показано, что наличие ослабленной зоны приводит к росту сжимающих напряжений в 1,7 раза и увеличению размеров области их концентрации.
Дальнейшими направлениями исследований являются: постановка и расчет вариантов с учетом наиболее крупных разломных зон; построение локальных моделей отдельных уступов борта карьера с заданием граничных условий из предыдущего этапа моделирования и уточненный анализ НДС.
Работа выполнена при финансовой поддержке Российского фовда фундаментальных исследований (проекты №07-05-13579; 06-05-64681 и 06-05-03006).
ЛИТЕРАТУРА
1. Зубков A.B. Геомеханика и криотехнология / УрО РАН. Екатеринбург, 2001. 335 с.
2. Козырев A.A. К проблеме устойчивости бортов глубоких карьеров в высоконапряженных массивах / А.А.Козырев, В.В.Павлов, В.А.Мальцев // Проблемы разработки глубоких карьеров: Материалы между народного симпозиума по открытым горным работам «Мир-ный-91» / Удачный, 1991. С. 129-134.
3.КозыревA.A. Трехмерное моделирование напряженно-деформированного состояния массива горных пород как основа прогноза удароогасносш на рудниках ОАО «Апатит» / А.А.Козырев, И.Э.Семенова, А.А.Шестов // Сборник докладов международной конференции «Геодинамика и напряженное состояние недр Земли». Новосибирск, 2007.
4. Stacey T.R. Slope Stability in High Stress and Hard Rock Conditions. Slope Stability 2007. P. 187-200.
5.Stead D.. CogganJ.S., ElmoD., YanM. Modelling Brittle Fracture in Rock Slopes - Experience Gained and Lessons Learned. Slope Stability 2007. P.239-252.
REFERENCES
1. Zubkov A. V. Geomechanics and Creotechnology / UrO RAS. Yekaterinburg, 2001. 335 p.
2. Kozyrev A.A. On slope stability of deep open-pit mines in highly-stressed rock mass / A.A.Kozyrev, V.V.Pavlov, V.A.Maltsev // Problems in mining of deep
open-pit mines: Proc. Int. Symp. on open-pit mining, «Mirny-91» /Udachny, 1991.P.129-134.
3. Kozyrev A.A. 3D-modeling of stress-strain state of rock mass as the basis for prediction of rockburst-hazard at ore mines of the Apatit Co. / A.A.Kozyrev, I.E.Semionova, A.A.Sliestov // Proc. Int. Conf. «Geodynamics and Stress State of the Earth's Interior». Novosibirsk, 2007.
4. Stacey T.R. Slope Stability in High Stress and Hard Rock Conditions. Slope Stability 2007. P. 187-200.
5.iSteadD., CogganJ.S., ElmoD., YanM. Modelling Brittle Fracture in Rock Slopes - Experience Gained and Lessons Learned. Slope Stability 2007. P.239-252.
