CC BY
22
УДК И.А. Поздеев, И.М. Поздеева, П.В. Васильев
622.834.2:
51®.б2: ИССЛЕДОВАНИЕ ЗАВИСИМОСТИ
519.254 ^
ПАРАМЕТРОВ НАПРЯЖЕННО-
ДЕФОРМИРОВАННОГО
СОСТОЯНИЯ
УГЛЕПОРОДНОГО МАССИВА ОТ РАСПОРА СЕКЦИИ КРЕПИ МЕХАНИЗИРОВАННОГО ЗАБОЯ
Обозначена актуальная научно-практическая задача и предложен вариант ее решения. Проведены исследования зависимости параметров напряженно-деформированного состояния углепородного массива в окрестности очистного забоя от распора механизированной крепи. Выявлены виды функций, наиболее адекватно описывающих зависимости исследуемых параметров. Ключевые слова: комплексно-механизированный забой, параметры напряженно-деформированного состояния геомассива, распор секции механизированной крепи, численное моделирование, метод наименьших квадратов, эмпирическая функция, коэффициент корреляции.
В соответствии с результатами анализа причин внеплановых остановок очистного забоя, проведенного авторами, следует, что в процессе ведения горных работ в комплексно-механизированном очистном забое (КМЗ) наблюдаются такие явления как вывалы негабаритных кусков пород кровли на забойный конвейер, отжим угля с поверхности очистного забоя, высыпание пород кровли между секциями крепи, при их передвижке, пучение пород почвы. Данные явления являются причинами остановки процесса добычи угля из-за необходимости разрушения негабаритных кусков пород для транспортирования по забойному конвейеру, неполной передвижки секции механизированной крепи и снижения поперечного сечения ходового отделения и лавы, что затрудняет проход людей по ходовому отделению механизированной крепи и создает дополнительное сопротивление для движения воздуха, а также повышается вероятность травмирования рабочих очистного забоя. Дезинтеграция пород кровли и почвы на блоки способствует миграции
ISSN 0236-1493. Горный информационно-аналитический бюллетень. 2016. № 7. С. 313-327. © 2016. И.А. Поздеев, И.М. Поздеева, П.В. Васильев.
метана по образовавшимся системам трещин, что приводит к загазированию очистного забоя и выработок выемочного столба.
В связи с этим актуальной научно-практической задачей является создание системы оперативного мониторинга и прогноза геомеханических и газодинамических явлений, протекающих в очистном забое при отработке выемочного столба, с целью обоснования технологических и технических решений, реализация которых обеспечит безопасную эффективную работу КМЗ.
Для решения поставленной задачи проведено численное моделирование НДС углепородного массива в окрестности очистного забоя с использованием пакета компьютерных программ, разработанных в Сибирском государственном индустриальном университете. Комплекс программ адаптирован для решения двумерных задач и выявления закономерностей изменения механических, реологических и плотностных свойств углепород-ного массива, а также прогноза дебита и давления метана [1—4].
По результатам моделирования выявлялась зависимость параметров НДС массива горных пород в окрестности очистного забоя от горно-геологических и горнотехнических факторов. В общем виде зависимость между исследуемыми параметрами можно выразить в виде следующей функции:
и = Дт, Н, а, Сж, Я, 3), (1)
где и — параметры НДС геомассива; т — мощность пласта; Н — глубина ведения горных работ; а — угол залегания; ссж — предел прочности пород/угля при сжатии; Я — распор механизированной крепи, 3 — скорость подвигания очистного забоя.
Объектом исследования принят углепородный массив шахтного поля ООО «Шахта «Есаульская» расположенный в Байда-евском месторождении Кузбасса в период отработки выемочного столба 26—28 по пласту 26а.
Глубина ведения горных работ изменяется от 410 до 580 м, вынимаемая мощность разрабатываемого пласта — 2,1 м, угол падения — от 1 до 100. Горный массив представлен алевролитами и песчаниками, а также включает 7 угольных пластов, ближайший из которых пласт 29а залегает в 130—150 м выше отрабатываемого пласта 26а. Шахта отнесена к сверхкатегорным по газу.
Схема подготовки шахтного поля — панельная. Система разработки столбовая, с полным обрушением пород кровли, подвига-нием длинного очистного забоя длиной 300 м по восстанию и выемкой угольного пласта с помощью комплекса 2КМ-138, включа-
ющего комбайн KSW-460NE, забойный конвейер RYBNIK-850, секции М-138, перегружатель ПСП-308.94. Схема выемки угля в лаве — односторонняя, снизу вверх по ходу струи свежего воздуха.
Для изучения геомеханических процессов, протекающих в окрестности очистного забоя, выбрана геометрическая модель в виде вертикального разреза по разведочной геологической скважине № 1689. Размер модели по горизонтальной оси OX составляет 410 м, по вертикальной оси OY — 690 м, где зона над-работки ограничена координатой Y = -115 м, а поверхность — Y = 575 м. Геометрическая модель объекта исследования делится на элементарные области (конечные элементы), в каждой из которых, искомая непрерывная величина описывается дифференциальным уравнением, в результате чего задача сводится к решению системы линейных уравнений [1—4]. Результаты деления модели представлены на рис. 1.
В статье рассматривается влияние величины распора механизированной крепи на изменение параметров НДС массива.
Моделирование геомеханических параметров проводилось при величине распора механизированной крепи от 600 до 1200 кН/м2 с шагом изменения параметра 200 кН/м2.
Для детального анализа и выявления зависимости НДС геомассива от распора механизированной крепи выбраны харак-
-200 -100 0 100 200 300 Расстояние от монтажной камеры, м
4776
Рис. 2. Схема КМЗ: 1—7 — характерные точки
терные точки, изменение НДС пород в которых наиболее полно характеризует состояние очистного забоя и вмещающих пород отрабатываемого пласта (рис. 2).
В процессе исследования определялись следующие параметры:
• главные вертикальные и горизонтальные составляющие тензора напряжений массива пород;
• коэффициенты концентраций вертикальных и горизонтальных напряжений;
• горизонтальные и вертикальные смещения;
• отношение остаточной прочности пород к исходной.
Для выявления вида зависимости параметров НДС геомассива от распора механизированной крепи при наличии экспериментальных данных, в том числе полученных при моделировании, выбран метод наименьших квадратов. Сущность метода заключается в определении коэффициентов эмпирической зависимости такой функции от аргумента, при которой сумма квадратов отклонений экспериментальных данных от расчетных будет наименьшей. Определение коэффициентов осуществляется посредством решения системы уравнений методом Гаусса либо матричным методом [5].
Для выявления функции, наиболее адекватно описывающую зависимость исследуемых параметров, значения, полученные в результате моделирования параметров НДС массива, были аппроксимированы линейными и нелинейными уравнениями. Виды используемых уравнений приведены в таблице.
На нижеследующих рисунках приведены результаты моделирования и функции, которые наиболее адекватно описывают зависимость исследуемых параметров, из которых следует, что:
1) в точке 1 (см. рис. 2):
• изменение величины распора механизированной крепи в пределах 600—1200 кН/м2 несущественно влияет на вертикаль-
Вид функций для выявления зависимостей параметров НДС геомассива от распора механизированной крепи
Название эмпирической зависимости Функция
Линейная y = ax + b
Квадратичная y = ax2 + bx + c
Логарифмическая y = a ■ ln(x) + b
Степенная y = axm
Показательная (экспоненциальная) y = a ■ emx
Дробно-рациональная y = x/ax + b
ные напряжения в обнаженной части кровли пласта у очистного забоя (рис. 3, а);
• распор механизированной крепи в большей степени влияет на изменение горизонтальных напряжений (рис. 3, б);
• увеличение распора механизированной крепи способствует снижению вертикальных и горизонтальных смещений пород кровли пласта в незакрепленной механизированной крепью части (рис. 4);
• изменение распора механизированной крепи в пределах 600—1200 кН/м2 не в значительной степени влияет на состояние пород кровли у призабойного пространства (рис. 5);
величина распора механизированной крепи К, hH/mZ
| SOD 700 900 1100 1300
• экспериментальные значения
Нид функции -Дробно рациональная (r=o,S98i]
и,.
б)* и<" X
§- &0О /И) 900 ПШ 1300
Г л: С:: механизированной «рея* R иН/М?
а)
Распор механизированной крепи И, кН/м2
500
700
поо
13оо
4, И
1 ? I ^
| п *
® 5 -29
-30
•
1 —2,5 ■ 1(Г?-Яг + 0,006383-
Экспериментальные значения Вид функции: -Квадратичная (п=1)
к 24,4
X X г\2
ь 24
Щ 5 £ Ч « ,4 23,8
1 = 23,6
16 23,4
га ей 23,2
| 3:
X 23
= 1,875 • ю-*-5?- 0,00136. В + 25,442
• Экспериментальные значения Вид функции: -Квадратичная {г~1|
8-
500
700
900
1100
1зоо
Распор механизированной крепи И, кН/м2
Рис. 4. Закон распределения вертикальных (а) и горизонтальных (б) смещений геомассива в зависимости от распора механизированной крепи в точке 1: г — коэффициент корреляции
2) в точках 2 и 3 (см. рис. 2):
• повышение величины распора механизированной крепи способствует увеличению сжимающих вертикальных и горизонтальных напряжений в породах кровли над верхним перекрытием крепи (рис. 6);
• увеличение распора механизированной крепи способствует уменьшению вертикальных смещений пород кровли пласта над верхним перекрытием секции крепи (рис. 7, а). Как видно из рис. 7, б, на изменение горизонтальных смещений величина распора механизированной крепи влияет в меньшей степени;
• изменение распора механизированной крепи в пределах 600—1200 кН/м2 в незначительной мере влияет на состояние
3,7484' 1(Г*-Яг-0,0П5-Д+ 4,5909
Рис. 6. Закон распределения вертикальных (а) и горизонтальных (б) напряжений геомассива в зависимости от распора механизированной крепи в точках 2, 3: г — коэффициент корреляции
Рис. 8. Закон распределения отношения остаточной прочности горных пород к исходной в зависимости от распора механизированной крепи в точках 2, 3: г — коэффициент корреляции
горных пород над верхним перекрытием крепи в районе передней гидростойки (рис. 8);
3) в точке 4 (см. рис. 2):
• при увеличении распора механизированной крепи сжимающие вертикальные напряжения и растягивающие горизонтальные напряжения в верхней части очистного забоя снижаются, но снижение напряжений незначительно (рис. 9);
• увеличение распора механизированной крепи способствует снижению смещений угольного массива в верхней краевой части выемочного столба (рис. 10);
• на горизонтальные смещения распор механизированной крепи влияет в меньшей степени (рис. 10);
а) 3
Величина распора механизированной крепи В, кН/м2
500 700 900 1100 1300
х *
з ь « | 5 | -9,45
' -9,5
9,55
-9,6
Е £
о, = 5,637- 10"5- Л +■ 0,00023
Экспериментальные значения Вид функции --Квадратичная 999989)
б)
в щ (1,6
2 к 9
| 5 | 0,58
I I 5 |1|
0,56 0,54
0,5>
аг - -6,4 -Н8 - 10 я- & + 2,379 ■ 10 5 - Я +
Я - 9,721 0,601
Экспериментальные Вид функции: -КвлдрйГИЧН,|Г| 995935)
500 /00 да «00 3 300
Распор глехзнижированной чрепн Я, кН/ы2
Рис. 10. Закон распределения вертикальных (а) и горизонтальных (б) смещений геомассива в зависимости от распора механизированной крепи в точке 4: г — коэффициент корреляции
• в точке 4 формируются положительные горизонтальные смещения, это свидетельствует о том, что смещения направлены вглубь угольного массива (рис. 10, б);
• увеличение распора механизированной крепи способствует увеличению остаточной прочности угольного массива в 0,5 м ниже кровли пласта, но увеличение незначительно (рис. 11);
4) в точке 5 (см. рис. 2):
• увеличение распора механизированной крепи способствует увеличению вертикальных сжимающих напряжений (рис. 12, а), но изменения вертикальных напряжений незначительны;
• увеличение распора механизированной крепи способствует снижению горизонтальных сжимающих напряжений в породах почвы пласта у очистного забоя (рис. 12, б);
• при увеличении распора механизированной крепи увеличиваются вертикальные смещения, т.е. происходит процесс
0,29
I ЪЪ £
I о. 5. £
11!
II £ I
□,25
Вт - Ь2' 19 ■ 1<Га ■ Л 2 н- Л,989 ■
* Экспериментальные значения ВИД функции:
— Квадратичная ((=0,99998)
1<Г5-Я +0,21918
500 700 900 1100
Распор механизированной крепи Я, кН/м2
1300
а) 5
Распор механизированной крепи Я, КН/м7
/ТО 900 1100 1Э(Ю
-1,315
-1,317
-1,319
%" 11 1,3 Л
-1,3 73
1\ -Ш
1,377
1,379
© 1,331
\
\
а„ = 4.0812-10 0 ■ Л - 9,2577 ■
• ЭКСПерИМРНТЗЛ^НМР
значения вид функции:
—■Колдратичная
276
б)
500
Распор механизированной крепи Я, нН/м2
700 900 1100
к Е 1Л 1,36
1 £ -1,37
От С В 1.38 -
% 1,39 -
1,4
сг = -6,387-1 о-8 ¡е- + 0.00018-
Экспериментальные
данные
ВИД функции!
- Киадрат ичнля (г=ЩЧ99/)
Рис. 12. Закон распределения вертикальных (а) и горизонтальных (б) напряжений геомассива в зависимости от распора механизированной крепи в точке 5: г — коэффициент корреляции
вдавливания почвы, но смещение пород незначительны, в результате чего можно утверждать, что влияние величины сопротивления крепи в исследуемых интервалах на почву несущественно (рис. 13, а);
• горизонтальные смещения при увеличении распора механизированной крепи уменьшаются (рис. 13, б), но уменьшение горизонтальных смещений так же, как и вертикальных, несу-
0,68
| и °'67
И "8 8 П И
ь а. х. г
1 О 3 а 0,65
л 5 * ь
О а * 4
Е О 0,64
<тос1.= 1,358-Ю"' -йг- 0,000212-« +
0.727е)
• Экспериментальные значении Вид функции:
-^Коадра!ичнан (г=0,9б37)
500 700 900 1100
Распор механизированной крепи Я, кН/мг
1300
Рис. 14. Закон распределения дезинтеграции горных пород в зависимости от распора механизированной крепи в точке 5: г — коэффициент корреляции
щественны (менее 1 мм) в результате чего можно утверждать, что влияние величины сопротивления крепи в исследуемых интервалах на почву незначительно;
• увеличение распора механизированной крепи способствует улучшению состояния пород почвы пласта (рис. 14). Как видно из графика, увеличение остаточной прочности пород наблюдается при распоре крепи равного величине 1200 кН/м2, но изменение незначительно (менее 0,05);
5) в точках 5 и 6 (см. рис. 2):
• увеличение распора механизированной крепи способствует повышению вертикальных напряжений пород почвы под основанием крепи (рис. 15, а).
а)
31
3Q5
s зс II ъ
¡9 5
I I 7S,5 1 |
* - -
и; <=-i,2 МО"*-Я* - 0.00124- R + 31.58-1
Б-10-"-^ - 0,00067 ■ Л + 28,147
--
70D 900 1100
Распор менэнизированнои ирепи R, кИ/гл2
б)
20,7
0 | 70.5 3 *
1 * № С с
х х
I I мл
& 5
i6,25 ■ 10 J0 /?г -0,000849- R + 21.O^N.
5кс перине н га льны е ы+ачепия '. trm.iii) вид функции:
-Квадратичная ((=0,99999]
Экспериментальные значении (тоикэ7) Рид функции:
Эксперимента ьныр значения (точ«Э6) Аид функции:
(г=Л,99999] Экспе ри мечт а л ьиые значения (точка71 Вид функции:
1300" — Квадратичная (in),9999]
500 700 900 1100
Распор механизированной крепи R, кН/м2
Рис. 16. Закон распределения вертикальных (а) и горизонтальных (б) смещений геомассива в зависимости от распора механизированной крепи в точках 6, 7: r — коэффициент корреляции
• горизонтальные напряжения при увеличении распора крепи переходят из растягивающих в сжимающие (рис. 15, б);
• при увеличении распора механизированной крепи снижаются вертикальные и горизонтальные смещения (рис. 16);
• увеличение распора механизированной крепи способствует повышению остаточной прочности пород почвы в районе передней гибростойки (рис. 17);
• увеличение распора механизированной крепи влечет повышение остаточной прочности пород почвы пласта в районе задней гидростойки крепи. Изменении прочности пород наб-
людается при увеличении распора крепи с 600 до 800 кН/м2. При дальнейшем увеличении распора крепи остаточная прочность практически не меняется. Изменение прочности колеблется в пределах 0,5%.
По результатам анализа полученных при моделировании данных сформулированы следующие выводы:
• распор механизированной крепи незначительно влияет на состояние пород незакрепленной крепью части кровли у очистного забоя;
• распор секции механизированной крепи, превышающий контактную прочность пород кровли, приводит к их разрушению, в результате чего возможны высыпания разрушенных пород между перекрытиями секций крепи в процессе их передвижки. С целью исключения возможных неблагоприятных и травмоопасных явлений во время рабочего процесса распор механизированной крепи рекомендуется устанавливать на таком уровне, чтобы напряжения, возникающие в массиве, не превышали контактную прочность пород;
• изменение распора механизированной крепи в пределах 600—1200 кН/м2 существенно не влияет на вероятность отжима угольной пачки с поверхности очистного забоя;
• при недостаточном распоре механизированной крепи может наблюдаться эффект заклинивания пласта породами кровли и почвы, что может привести к разрушению угля при формировании напряжений, превышающих предел прочности угля;
• увеличение распора механизированной крепи способствует уменьшению вероятности пучения пород почвы пласта и выдавливанию забойного скребкового конвейера у очистного забоя;
• увеличение распора крепи способствует уменьшению заклинивания пласта породами кровли и почвы и снижает вероятность перемещения секций мехкрепи в сторону очистного забоя;
• увеличение распора крепи в пределах 600—1200 кН/м2 несущественно повлияет на состояние пород почвы пласта под скребковым конвейером очистного забоя;
• увеличение распора механизированной крепи способствует уменьшению пучения почвы под основанием секции крепи, но, при изменении величины распора крепи в пределах 600— 1200 кН/м2 влияние на почву отрабатываемого пласта незначительно;
• при величине распора мехкрепи 600—1200 кН/м2 изменение прочности пород незначительно и находится в пределах
0,63—0,70. Однако при отношении остаточной прочности пород к исходной равного и более 0,7 деформации горных пород можно отнести к упругим. В результате чего можно предположить, что увеличение распора крепи способствует улучшению состояния пород почвы и, как следствие, исключению образования систем трещин, которые могут стать путями миграции метана из пород почвы и вероятных суфлярных выделений.
СПИСОК ЛИТЕРАТУРЫ
1. Риб С.В., Фрянов В.Н. Разработка комплекса проблемно-ориентированных программ для численного моделирования напряженно-деформированного состояния неоднородных угольных целиков // Горный информационно-аналитический бюллетень. — 2015. — № 3. — С. 367-371.
2. Риб С. В., Волошин В. А., Максимов А. А., Борзых Д. М., Никитина А. М., Фрянов В. Н. Численное моделирование методом конечных элементов напряженно-деформированного состояния углепородно-го массива при переходе очистным забоем передовой выработки // Горный информационно-аналитический бюллетень. — 2015. — № 3. — С. 414—422.
3. Петрова О. А., Васильев П. В., Фрянова О. В., Фрянов В. Н. Модель формирования напряжений, деформаций и повреждений в углепород-ном массиве при интеграции гравитационного и геотектонического полей напряжений / Материали за VIII международна научна практична конференция «Новината за напреднали наука — 2012». Т. 26. Технологии. — Бял ГРАД-БГ ООД. София: 2012. — С. 10—19.
4. Фрянов В. Н., Петрова О. А., Петрова Т. В. Комплекс проблемно-ориентированных программ для моделирования формирования и распределения опасных зон в газоносном геомассиве: свидетельство о регистрации электронного ресурса. № 21123, дата регистрации 03 августа 2015 года // Хроники объединенного фонда электронных ресурсов «Наука и образование» № 08—09 (75—76) август-сентябрь 2015. — С. 2. Режим доступа: http://ofernio.ru/portal/newspaper05.php
5. Кремер Н.Ш. Теория вероятностей и математическая статистика: учебник для студентов вузов. — 3-е изд., перераб. и доп. — М.: ЮНИ-ТИ-ДАНА, 2010. — 551 с. EES
КОРОТКО ОБ АВТОРАХ
Поздеев Игорь Андреевич — аспирант, помощник начальника участка ВТБ, e-mail: igor-pozdeev0910@mail.ru,
Сибирский государственный индустриальный университет, Поздеева Ирина Михайловна — ведущий инженер, e-mail: irina.pozdeeva4387@yandex.ru,
ООО «Новокузнецкий институт геофизики и проектирования», Васильев Павел Валентинович — кандидат технических наук, генеральный директор, e-mail: vasilyevpv@bk.ru, ООО «Сибирская экспертная организация».
Gornyy informatsionno-analiticheskiy byulleten'. 2016. No. 7, pp. 313-327.
I.A. Pozdeev, I.M. Pozdeeva, P.V. Vasil'ev RESEARCH OF DEPENDENCE OF PARAMETERS OF THE INTENSE DEFORMED CONDITION OF THE ROCK MASS ON RESISTANCE OF SECTION FIX OF THE MECHANIZED FACE
The actual scientific and practical task has been designated and the version of her decision has been offered. Researches of dependence of parameters of the condition of the rock mass in a neighborhood of a clearing face from a resistance mechanized are conducted fix. Types of the functions which are most adequately describing dependences of the studied parameters have revealed.
Key words: the complex mechanized face, parameters of the condition of the rock mass, resistance of the section mechanized fix, numerical modeling, least square method, empirical function, correlation coefficient.
AUTHORS
Pozdeev I.A., Graduate Students, Assistant Chief of Ventilation and Occupational Safety Department, e-mail: igor-pozdeev0910@mail.ru,
Siberian State Industrial University, 654007, Novokuznetsk, Russia, Pozdeeva I.M., Leading Engineer, e-mail: irina.pozdeeva4387@yandex.ru, Novokuznetsk Institute of Geophysics and Design, 654066, Novokuznetsk, Russia,
Vasil'ev P.V., Candidate of Technical Sciences, General Director, e-mail: vasilyevpv@bk.ru,
Siberian Expert Organization LTD, 653000, Prokop'evsk, Russia. REFERENCES
1. Rib S. V., Fryanov V. N. Gornyy informatsionno-analiticheskiy byulleten'. 2015, no 3, pp. 367—371.
2. Rib S. V., Voloshin V. A., Maksimov A. A., Borzykh D. M., Nikitina A. M., Fryanov V. N. Gornyy informatsionno-analiticheskiy byulleten'. 2015, no 3, pp. 414—422.
3. Petrova O. A., Vasil'ev P. V., Fryanova O. V., Fryanov V. N. Materiali z,a VIII mezh-dunarodna nauchna praktichna konferentsiya «Novinata z,a naprednali nauka 2012». T. 26. Tekhnologii (Материали за VIII международна научна практична конференция «Новината за напреднали наука 2012», vol. 26. Технологии), Byal GRAD-BG OOD. Sofiya, 2012, pp. 10—19.
4. Fryanov V. N., Petrova O. A., Petrova T. V. Kompleks problemno-orientirovannykh programm dlya modelirovaniya formirovaniya i raspredeleniya opasnykh zon v gazonos-nom geomassive: svidetel'stvo o registratsii elektronnogo resursa, no 21123, data regis-tratsii 03 avgusta 2015 goda. Khroniki ob"edinennogo fonda elektronnykh resursov «Nauka i obrazovanie» № 08—09 (75—76) August-September 2015, pp. 2, available at: http://ofern-io.ru/portal/newspaper05.php
5. Kremer N.Sh. Teoriya veroyatnostey i matematicheskaya statistika: uchebnik dlya studentov vuzov. 3-e izd. (Probability theory and mathematical statistics: Textbook for high schools, 3rd edition), Moscow, YuNITI-DANA, 2010, 551 p.
UDC 622.834.2: 519.62: 519.254
