Разработка алгоритма расчета метанообильности очистного забоя с учетом геомеханических процессов в углепородном массиве

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

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

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

Решение обозначенной задачи осуществлялось в два этапа:

1 - численное моделирование напряженно-деформированного состояния (НДС) углепо-родного массива в окрестности очистного забоя;

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

Моделирование НДС геомассива в окрестности очистного забоя осуществлялось с использованием пакета компьютерных программ, разработанных на кафедре геотехнологии в Сибирском государственном индустриальном университете. Комплекс программ адаптирован для решения двумерных задач геомеханики и выявления закономерностей изменения механических, реологических и плотностных свойств углепородного массива [3 - 9].

В качестве объекта исследования принят углепородный массив Байдаевского месторож-

дения в пределах шахтного поля ООО «Шахта «Есаульская», Кузбасс. Исследования проведены в период отработки выемочного столба 26-28 по пласту 26а.

Средняя глубина ведения горных работ составляет 490 м, вынимаемая мощность разрабатываемого пласта 2,1 м, угол падения от 1 до 10°. Горный массив представлен алевролитами и песчаниками, а также включает семь угольных пластов, ближайший из которых пласт 29а залегает в 130 - 150 м выше отрабатываемого пласта 26а. Шахта отнесена к сверхкатегорной по газу метану. Система разработки столбовая, с полным обрушением пород кровли, подвига-нием очистного забоя длиной 300 м по восстанию. Схема выемки угля в лаве - челноковая, снизу вверх по ходу струи свежего воздуха.

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

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

—пч з -^пп + —т

где 1оч.з - метанообильность очистного забоя, м3/мин; 1пл - метанообильность отрабатываемого пласта, м3/мин; /т.г.м. - метанообильность транспортируемой по забойному конвейеру горной массы, м3/мин.

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

3 =

к ёР(х) ц ёх

где 3 - скорость фильтрации газа, м/с; к - коэффициент проницаемости массива; ц - динамическая вязкость газа, Па с; Р - давление газа, МПа; х - расстояние до поверхности очистного забоя, м.

Рис. 1. Вариант результатов моделирования изменения горизонтальных деформаций угольного массива впереди очистного

забоя для заданных условий ведения горных работ (е-10-3): 0 - 18 - изолинии горизонтальных деформаций; I, II, III - характерные линии изменения деформаций

В реальных условиях давление газа и проницаемость угольного массива не постоянны и меняются в процессе ведения очистных работ. С целью учета меняющихся физико-механических свойств массива предлагается использовать формулу, предложенную Б.Г. Тарасовым [11]:

k = 2,08d:

m0 + du/dx 1 + du/dx

где d - эквивалентный диаметр породообразующего зерна; т0 - природная пористость массива; du/dx -деформация массива.

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

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

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

очистных работ предлагается определить по формуле

е0б = j s(x)dx,

где £об - суммарные деформации в краевой части угольного пласта; s - распределение деформации в краевой части пласта по направлению x; l - зона активной фильтрации газа, м.

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

Введем локальную систему координат Olh с началом координат на пересечении линии очистного забоя и почвы отрабатываемого пласта (рис. 1). Ось абсцисс указывает расстояние до поверхности очистного забоя, ось ординат - расстояние до почвы пласта. Согласно результатам проведенных исследований, деформации угольного массива изменяются не только в горизонтальной плоскости, но и в вертикальной (рис. 1). В этой связи определены функции, описывающие зависимость деформаций от расстояния до очистного забоя в следующих характерных горизонтальных линиях: у кровли пласта (на рис. 1 линия обозначена цифрой I); в центральной части мощности пласта (на рис. 1 линия обозначена цифрой II); у почвы отрабатываемого пласта (на рис. 1 линия обозначена цифрой III). Общий вид зависимостей следующий:

Sk = al2 + ü\l + a2;

о

8Ср = Ы2 + ЪХ1 + Ъ2;

8п = С12 + С\1 + С2,

Р, = Ра-

где ек, 8ср и 8п - изменение деформаций по характерным линиям у кровли пласта, в центральной части мощности пласта и у почвы пласта соответственно; I - зона активной фильтрации газа, м; аь Ъь сг- - коэффициенты эмпирических зависимостей по отдельным плоскостям.

Из функций (5) - (7), описывающих изменение деформаций по отдельным плоскостям, получена обобщенная зависимость изменения деформаций в пределах всей зоны активной фильтрации газа в краевой части пласта:

8об = (ук2 + ухк + у 1) /2 + (вк2 + в2к + в2) I +

+

где у-, вг-, - коэффициенты эмпирической зависимости.

С учетом полученных зависимостей функция (4) примет следующий вид:

I к

8об = Ц (ук2 + у к + у]) I2 + (вк2 + в2к + в2) I +

0 0

+ + g2к + g2)dldк.

Для определения градиента газового давления, используемого в формуле (2), предлагается применить зависимость, предложенную В.А. Колмаковым [10]:

Для того, чтобы определить среднее значение градиента газового давления в зоне повышенной фильтрации пласта предлагается следующая формула:

Рср =

I

где Рср - среднее значение градиента газового давления, МПа/м; Р'х - первая производная от газового давления, МПа; I - зона активной фильтрации пласта, м.

Интенсивность метановыделения с поверхности очистного забоя зависит от времени его обнажения. Через определенный промежуток времени, после выемки полосы угля, глубина зоны активной фильтрации пласта увеличивается, а давление газа в краевой части отрабатываемого пласта снижается, что приводит к уменьшению дебита выделяющегося метана. Таким образом, формула (10) справедлива для давления газа на вновь обнаженном участке очистного забоя, который соответствует длине полосы угля, вынимаемой очистным комбайном за одну минуту. С целью учета изменения градиента газового давления с течением времени формулу (10) предлагается модернизировать следующим образом:

Рх = Р0 + (Р -Р0) • е

212

где Рх и Р0 - соответственно давление газа на любом расстоянии от обнаженной поверхности пласта и начальное давление в пласте, МПа; Р1 - давление газа на обнаженной поверхности пласта, МПа; х - расстояние до поверхности пласта; Щ) - глубина зоны разгрузки пласта.

Начальное давление газа в пласте определяется при помощи формулы Г.Д. Лидина [12]

Р0 = 0,01(Яг - Я0)13+2,

в которой Р0 - давление газа в пласте на расчетной глубине, МПа; Нг и Н0 - соответственно глубина разработки и зоны газового выветривания, м.

Условимся считать, что давление газа на поверхности очистного забоя в течение процесса фильтрации метана с краевой части отрабатываемого пласта равно атмосферному давлению:

Рх = Р0 + (Р -Р0) • е

2(1 2(г )Ф(< ))2

где - функция ползучести угля.

Изменение реологических свойств угольного массива с течением времени предлагается определить при помощи формулы Ж.С. Ержа-нова [13]:

Ф,

1 +

ы

1-а

1 - а

где 5 и а - параметры ползучести, характеризующие реологические свойства угля; ^ - время цикла выемки угля, с.

Параметр 5 рассчитывается так:

5 = 4,73 10-3 - 6,9-10-9- Е0;

здесь Е0 - модуль упругости угля, т/м [14], определяемый как

0

2

х

х

Ео = 200 осж,

где осж - предел прочности угля при сжатии, т/м2.

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

Ч = — Рр,

Ц

где ц - удельный поток метана в угольном пласте в направлении к свободной поверхности забоя, м3/(м2 мин); к(В) - коэффициент проницаемости с учетом изменения деформаций угольного массива; ц - динамическая вязкость газа, Пас.

Модернизация формулы (2) состоит в следующем:

- для определения коэффициента проницаемости массива применяется формула Б.Г. Тарасова (3), при вычислении по которой используются результаты численного моделирования и функции изменения деформаций угольного массива впереди очистного забоя (5) - (7), применение которых позволяет учитывать изменение фильтрационных свойств массива под влиянием геомеханических процессов;

- при вычислении градиента газового давления за основу принята формула В.А. Колма-кова (10), в которой учитывается изменение давления в зависимости от расстояния до поверхности очистного забоя, однако она была модернизирована (14) за счет того, что глубина зоны активной фильтрации газа принимается по результатам проведенного моделирования, а учет влияния фактора времени на градиент газового давления осуществляется по функции ползучести Ж.С. Ержанова (15), что позволило частично автоматизировать вычисление градиента давления и снизить трудоемкость расчета;

- синтез уравнений Дарси, Б.Г. Тарасова, В.А. Колмакова, Ж.С. Ержанова и результатов численного моделирования обеспечит вычисление дебита метана, выделившегося с поверхности очистного забоя, с использованием вместо постоянных параметров, которые не в полной мере соответствуют реальным горнотехнологическим условиям, переменных значений деформаций и пористости массива. Предложенный подход позволит учитывать

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

На основе модернизированных уравнений Дарси, Б.Г. Тарасова, В.А. Колмакова общее количество метана, выделившегося с поверхности отрабатываемого пласта в исходящую струю, предлагается определять по формуле

Qпл = Яв.об + Яюб,

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

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

Яв.об = Цв.об ^в.об m,

где цв.об - удельный дебит метана, выделившегося с вновь обнаженной поверхности очистного забоя м3/(м2мин); £в.об - длина вновь обнаженного участка очистного забоя, м; т -мощность отрабатываемого пласта, м.

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

Яюб = Цюб Ьюб т,

где ц1об - удельный дебит метана, выделившегося с поверхности очистного забоя до начала нового цикла выемки угля, м3/(м2 мин); Ь1об -длина очистного забоя за вычетом вновь обнаженного участка, м.

С целью учета метановыделения из транспортируемого по забойному скребковому конвейеру угля предлагается использовать поправочный коэффициент Кту, который определяется следующим образом [1]:

К = а 0,6 4 /от + /птк + /кш ;

т у 2 4 60&т.оч 60&т.п^ 60&т.кш

здесь Кту - коэффициент, учитывающий степень дегазации отбитого угля при его транспортировании по выработкам участка; а2 - коэффициент, характеризующий газоотдачу отбитого угля; /оч - длина забойного скребкового конвейера, м; 5т.оч - скорость транспортирования угля по очистному забою, м/с; /птк - длина подлавного перегружателя, м; 5птк - скорость транспортирования угля по подлавному пере-

Рис. 2. Схема вентиляции выемочного участка 26-28 ООО «Шахта «Есаульская», Кузбасс: 1 - направление движения очистного забоя; 2 - опережающая сбойка; 3 - ходок 26-53; 4 - вентиляционный штрек 26-28;

5 - газодренажный ходок 26-53; 6 - вентиляционный штрек 26-31; 7 - магистральный путевой штрек 26-5; 8 - изолирующая перемычка; 9 - свежая струя воздуха; 10 - исходящая струя воздуха; 11 - выработанное пространство; 12 - датчик метана системы АГК, с которого считывались показания газа метана при проведении эксперимента;

13 - место замера скорости воздуха

гружателю, м/с; 1к.ш - длина конвейера, расположенного в конвейерном штреке, м; 3к.ш -скорость транспортирования угля по конвейерному штреку, м/с.

Коэффициент, характеризующий газоотдачу отбитого угля а2, равен

а2 = 0,25 а3,

где а3 - коэффициент, характеризующий газоотдачу угля в массиве.

Таким образом, формула (1) примет следующий вид:

предназначенного для непрерывного контроля концентрации метана на угольных предприятиях.

Таким образом, за период проведения эксперимента сделано более 100 замеров скорости воздуха и проведен газовый анализ 90 забое-дней.

Согласно «Инструкции по контролю состава рудничного воздуха, определению газообильности....» [15] средний расход газа метана, проходившего в пункте измерения в исходящей струе очистного забоя за определенный период времени, определяется по формуле

^оч.заб

- 1пл ( 1 + К.у),

L, = 0,01

X 0уч X см

где /оч.заб - метанообильность очистного забоя, м3/мин; 1пл - метанообильнось отрабатываемого пласта м3/мин.

С целью оценки адекватности предложенного алгоритма расчета ожидаемой метано-обильности очистного забоя проведен шахтный эксперимент. Для исследования аэрогазовой обстановки за основу принят инструментальный метод. В ремонтные и добычные смены в вентиляционном штреке 26-28, в 20 м выше сопряжения с очистным забоем, при помощи переносного анемометра АПР-2, предназначенного для измерения средней скорости движения воздушных потоков в шахтах, были проведены замеры скорости воздуха (рис. 2). Кроме того, использовались показания скорости воздуха, полученные горными мастерами участка вентиляции и техники безопасности (ВТБ). Учет концентрации метана в исходящей струе очистного забоя осуществлялся при помощи стационарного датчика метана (ДМС01),

где 1м - средний расход газа метана, проходившего в пункте измерения, м3/мин; 0уч -расход рудничного воздуха в исходящей струе выемочного участка, м3/мин; См - концентрация метана в исходящей струе выемочного участка по данным системы аэрогазового контроля (АГК), %; пв - количество измерений расхода рудничного воздуха; пм - количество учитываемых значений метана за исследуемый период по данным системы АГК.

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

п

п

в

м

12108 -6 -4 -2 -

х Ж ХЙЙЙВЯОЗ* хх

• M

• • •

^ | 0 1 2 3 4 5 6 Скорость движения очистного комбайна, м/мин

Рис. 3. Зависимости фактической (•) и расчетной (х) абсолютной метанообильности очистного забоя от скорости движения очистного комбайна

му алгоритму рассчитана ожидаемая абсолютная метанообильность для тех же скоростей.

На рис. 3 изображены результаты расчета фактической и расчетной абсолютной метанообильности очистного забоя 26-28, откуда следует, что разработанный алгоритм обеспечивает прогноз максимально возможной метано-обильности очистного забоя.

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

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Инструкция по применению схем проветри-

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

2. Руководство по проектированию вентиляции угольных шахт. - Макеевка-Донбасс, 1989.

3. Риб С.В., Фрянов В.Н. Разработка комплекса

проблемно-ориентированных программ для численного моделирования напряженно-деформированного состояния неоднородных угольных целиков // Горный информационно-аналитический бюллетень. 2015. № 3. С. 367 - 371.

4. Риб С.В., Волошин В.А., Максимов А.А.,

Борзых Д.М., Никитина А.М., Фрянов В.Н. Численное моделирование методом конечных элементов напряженно-деформированного состояния углепород-ного массива при переходе очистным забо-

ем передовой выработки // Горный информационно-аналитический бюллетень. 2015. № 3. С. 414 - 422.

5. Петрова О.А., Васильев П.В., Фрянова О.В.,

Фрянов В.Н. Модель формирования напряжений, деформаций и повреждений в углепородном массиве при интеграции гравитационного и геотектонического полей напряжений. Материали за VIII международна научна практична конференция «Новината за напреднали наука - 2012». Т. 26. Технологии. - Бял ГРаД-БГ ООД. София: 2012. С. 10 - 19.

6. Фрянов В Н., Петрова О.А., Петрова Т В.

Комплекс проблемно-ориентированных программ для моделирования формирования и распределения опасных зон в газоносном геомассиве: свидетельство о регистрации электронного ресурса. № 21123, дата регистрации 03 августа 2015 года // Хроники объединенного фонда электронных ресурсов «Наука и образование» № 08-09 (75-76) август-сентябрь 2015. С. 2. Режим доступа: http://ofernio.ru/portal/ newspaper05.php.

7. Fryanov V.N., Pavlova L.D. Prospects for ap-

plication of robotic mechano-hydraulical excavation of gas-bearing coal layers at great depths // IOP Conference Series: Earth and Environmental Science. 2016. Vol. 45. (012002). Режим доступа:

http://iopscience.iop.org/issue/1755-1315/45/1.

8. Tsvetkov A.B., Pavlova L.D., Fryanov V.N.

Construction of the approximant of complete diagram for rock deformation // IOP Conference Series: Earth and Environmental Science. 2016. Vol. 45. (012009). Режим доступа: http://iopscience. iop.org/issue/1755-1315/45/1.

9. Pristupa Yu.D., Fryanov V.N., Pavlova L.D.

Construction of a conceptual model of transport system for a coal mining region // IOP Conference Series: Earth and Environmental Science. 2016. Vol. 45. (012014). Режим доступа: http://iopscience.iop. org/issue/1755-1315/45/1.

10. Колмаков В.А. Метановыделение и борьба с ним в шахтах. - М.: Недра, 1981. - 134 с.

11. Тарасов Б.Г., Колмаков В.А. Газовый барьер угольных шахт. - М.: Недра, 1978. - 200 с.

12. Лидин Г.Д. К вопросу о закономерности выделения метана из угля, отторгнутого из массива. - В кн.: Управление газовыделением и пылеподавлением в шахтах. - М.: Наука, 1972. С. 50 - 52.

СибГИУ, 1988. - 139 с.

15. Инструкция по контролю состава рудничного воздуха, определению газообильности и установлению категорий шахт по метану и/или диоксиду углерода. Утверждена приказом № 704 Федеральной службы по экологическому, технологическому и атомному надзору 06.12.2012.

© 2017 г. И.А. Поздеев Поступила 1 февраля 2017 г.

УДК 622.831

А.А. Исаченко, А.А. Петров ООО «Распадская угольная компания»

ОБОСНОВАНИЕ ПО РЕЗУЛЬТАТАМ ЧИСЛЕННОГО МОДЕЛИРОВАНИЯ ПАРАМЕТРОВ КРЕПИ СБОЕК КАПИТАЛЬНЫХ ВЫРАБОТОК, ПРОЙДЕННЫХ В НЕОДНОРОДНОМ УГЛЕПОРОДНОМ МАССИВЕ

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

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

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

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

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

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

Вентиляционный ^УКтнпл ЫйТранспортный уклон пл. 1 / — уклон пл. 1

Транспортный уклон пл. 2

20

§ с <4 ¡4 10

^ а

;т о 0

с? «3

£ «

о ^ г ч -10

40

о «к -20

10 0 10 20 30 40 50 60 70 80 90 100110 Расстояние от западного бока вентиляционного уклона пласта 1, м

Рис. 1. Схема расположения уклонов

Цель настоящей работы - обосновать рекомендации для обеспечения устойчивости горных выработок при взаимном влиянии элементов сложной системы «весьма сближенные пласты 1 и 2, четыре охраняемые угольными целиками уклона, сбойки между уклонами пластов».

Для достижения поставленной цели решены следующие задачи:

- разработана методика исследований и проведена адаптация пакета компьютерных программ к горно-геологическим и горнотехническим условиям весьма сближенных пластов 1 и 2;

- проведено численное моделирование напряженно-деформированного состояния уг-лепородного массива в окрестности сбоек уклонов;

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

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

Методы исследований: экспертная оценка, математическое моделирование численными методами механики горных пород, обобщение результатов исследований.

Для моделирования напряженно-деформированного состояния геомассива использованы следующие исходные данные:

- план горных выработок по пласту 2;

- геологические колонки по геологоразведочным скважинам (рис. 2);

- рабочие чертежи пространственного положения следующих уклонов: вентиляционно-

го уклона пласта 1, конвейерного уклона пласта 2, транспортно-го уклона пласта 1, транспортного уклона пласта 2;

- прогнозный геологический разрез по транспортному уклону пласта 1.

Мощность пород между пластами 1 и 2, представленная переслаиванием аргиллитов и алевролитов, по трассе уклонов составляет 7,09 м. В кровле пла-ста 2 залегает алевролит мелкозернистый сухой без расслоения, а в почве пла-ста 1 - алевролит мелкозернистый (рис. 2).

Учитывая отсутствие результатов натурных измерений деформаций по-род в сбойках уклонов в качестве основного при исследовании принят метод конечных элементов [1 - 3] с использованием авторских пакетов программ KLATRAT, разработанных на кафедре геотехнологии Сибирского государственного индустриального университета [4 - 6]. В качестве исходных приняты данные геологической колонки. Геометрическая модель представлена в виде вертикального разреза, перпендикулярного осям уклонов. Начало условной системы координат принято на пересечении кровли пласта 2 и западного бока вентиляционного уклона пласта 1 (рис. 1). Размеры модели по простиранию приняты с учетом возможного влияния очистного выработанного пространства выемочных столбов пласта 2 на состояние пород в окрестности уклонов. Размеры геометрической модели по простиранию - 480 м. Глубина залегания пласта 2 - 270 м.

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

Первый этап проведения сбойки с присеч-кой пород кровли пласта 1 (рис. 3). Выявлено начало активного влияния пласта 2 на устойчивость пород кровли сбойки на расстоянии 4,5 м от почвы пласта 2.