CC BY
44
ГОРНОЕ ДЕЛО
УДК 669.187
ПРОГНОЗ ГАЗОВЫХ СИТУАЦИЙ В УГОЛЬНЫХ ШАХТАХ В ПЕРИОДЫ ПАДЕНИЯ АТМОСФЕРНОГО ДАВЛЕНИЯ
Н.М. Качурин, С.А. Воробьев, О.А. Афанасьев, Д.Н. Шкуратский
Рассмотрены математические модели диффузионных процессов в угольных пластах и выработанных пространствах в периоды экстренных газовыделениий. Получены значения дебита газовой смеси для периодов уменьшения атмосферного давления.
Ключевые слова: диффузионные процессы, газовыделение, выработанное пространство, фильтрация газа.
В связи с незначительней избыточным давлением газовой смеси в выработанных пространств, газообильности очистных и подготовительных участков зависит от колебаний статического давления воздуха в горных выработках[5]. При трехсторонней примыкании выработанного пространства к очистным выработкам (вентиляционный штрек - лава - конвейерный штрек) дебит углекислого газа в периоды уменьшения атмосферного давления достигает 1,5 м /мин. Уменьшение статического давления воздуха создает избыточное давление газовой смеси, находящейся в порах и пустотах обрушенных пород. В результате изотермического расширения газовая смесь фильтруется к поверхности контакта выработанного пространства с вентиляционной струей [4], [9].
В качестве основного модельного структурного элемента, поглощающего газы из атмосферы и выделяющего газы в рудничную атмосферу шахт Подмосковного бассейна [1], [2], является пористая сорбирующая среда, которая должна отражать физико-химические свойства, присущие как породоугольному массиву нарушенной структуры (отбитый уголь, выработанное пространство), так и породоугольному массиву ненарушенной структуры (обнажение поверхности разрабатываемого пласта и вмещающих пород.) [6], [7].
Рассмотрена модель динамики концентрации для кислорода [8]. Если принять, что в призабойном пространстве по всей его длине концентрация кислорода одинакова и зависит только от времени [3], то уравнение баланса количества кислорода в проветриваемом объеме можно записать в
152
следующем виде:
— = — (I + Q(c - с) ),
dt W г o ''
(1)
I = -n FdC г dx
+1м.в.'с,
x=0
c(0) = сн = const.
Рассмотрим частный случай, когда 1м в (t) = I = const, а для по-
dc
верхности обнажения угольного пласта
dx
x=0
Ik = c
k¥ в
С1
4DTo,
тогда имеем:
) = с0 - Вк^ИкКоIсрс + — е ^
и . (2)
На рис. 1 изображен график зависимости динамики концентрации
кислорода в призабойном пространстве, в период экстренных газовыделе-
ний.
Рис. 1. График концентрации кислорода в атмосфере призабойного пространства
Дебит углекислого газа из выработанного пространства определяется в соответствии с законом Дарси по формуле:
I’’ (t):
св.пhoLkn f j(t) dt
J Ґ4. _\U ’
(3)
гатфСО ^ -т)1,5
где - концентрация углекислого газа в газовой смеси выработанного пространства; ^, L - высота и длина поверхности контакта выработанного пространства с вентиляционной струей.
Математическая модель динамики концентрации кислорода в атмо-
153
t
сфере подготовительной выработки в период экстренных газовыделений учитывает, что диффузионный поток кислорода в угольном пласте определяется по закону Фика и для одномерного процесса разбавления кислорода «мертвым» воздухом с учетом конвективно-турбулентной диффузии можно записать, что на исходящей струе концентрация кислорода будет изменяться по следующему закону:
с(Ьп.в. ^ = сн |ехр(-М) -1ехр[-кг Ц - т)] • 0(ЬП В., т)ёт| + 0,5с0О(Ьпв., ^ +
г 1 1
1 - ехр(-кг1)] -1 о(Ьп.в^ т)ат+1 ехр[-кг(1 - т)] • о(Ьп.в^ т)ат I, (4)
+л
где
0(х,1) =
ехр (-хБт 0,5) ег&
х
л/(0,25О,ю2 + р + К „)!
хегій
2^
0,5х(Б т І)-0,5 +уІ (0,25Бтю2 +Р + К п )1
ехр (хЭт0,5 )х
1.
1;
кг = (0 + 1ср )П-1 + Кп, Л = (ОСо + 1срС1)^1[(д + ^ )ОТ1 + Кп Г1; Ьпв. - протяженность подготовительной выработки; сн, с0 - начальная концентрация кислорода и концентрация кислорода в свежей струе воздуха соответственно; 0 - среднее количество воздуха; Бт - коэффициент турбулентной диффузии кислорода в вентиляционной струе.
Тогда условие аэрогазодинамической безопасности по фактору обескислороживания воздуха имеет следующий вид: с(ЬПВ ,1) £ 0,2. Анализ показал, что, в конечном счете, устанавливается некое стационарное распределение концентрации кислорода. Поэтому пренебрегая турбулентной диффузией целесообразно использовать следующее уравнение:
Озп^С
= -(Р + К п )
А
Ь + К п
(5)
8-л(х)ёх
где Qзп - количество воздуха, подаваемое в подготовительный забой; Б - площадь поперечного сечения подготовительной выработки; ^(х) - коэффициент доставки воздуха.
Решение уравнения (5) с учетом того, что с(Ьпв ) = ПДК, получено алгебраическое уравнение для расчета количества воздуха, необходимого для проветривания подготовительной выработки в периоды экстренных газовыделений:
0з.пС0 + ІС1
0з.п.+ I + К п8п.в.Ьп,
■ +
0з.пС0 + ІС1
0,.+ I + К п8п.в.Ьп,
х
в. У
х ехр
(Р+К п )3П.,
0з.п. 0
154
'п.в.
| л(х)ёх
= ПДК.
(6)
0
0
С
0
Вычислительные эксперименты с использованием уравнения (6) показали, что в период экстренных газовыделений потребное количество воздуха для условий Подмосковного бассейна составляет 230 - 420 м /мин.
Математическая модель динамики концентрации кислорода в атмосфере очистного участка в период экстренных газовыделений с учетом, того, что диффузионный поток кислорода в угольном пласте определяется по закону Фика, а для одномерного процесса разбавления кислорода «мертвым» воздухом за счет конвективной диффузии позволила получить следующую теоретическую закономерность:
с(х,0 = (сн - А) ехр[-(Р + Кп1 Со - (сн - А) ехр
-(Р+К п)
і
X
и у
-X
х ехр
Р+К п
X
и
(7)
Зависимость (7) показывает, что среднее распределение концентрации кислорода по ходу движения воздуха стремиться к некому стационарному продольному профилю концентрации. Если учесть, что на исходящей концентрация кислорода должна быть не менее ПДК, то можно определить количество воздуха, необходимое для проветривания очистного забоя в период экстренных газовыделений по кислородному фактору Ок из решения следующего уравнения:
Ок Со + 1с1
ОК со + Ісі
ОК +1 + К п8(ів.п + 10Ч)
(Р + К п )8
ОК + I + К пБ(1в.п + 1оч )
X
х ехр
О
(1в.п + 1оч )
ПДК.
(8)
Вычислительные эксперименты показывают, что дополнительное количество воздуха составляет 45 - 87 % от расчетной величины.
Решение алгебраических уравнений (6) и (8) реализованы на языке Си в виде отдельных модулей, изображенных на рис. 2.
Расчеты количества воздуха для проветривания подготовительных выработок по фактору обескислороживания воздуха с использованием данной программы.
На рис. 3 представлены результаты расчета количества воздуха для проветривания подготовительной выработки в виде сводной диаграммы. Исходные данные брались по недействующим шахтам Тульской области.
Таким образом, на основании выполненных теоретических и практических исследований уточнены закономерности колебаний атмосферного давления, уточнена математическая модель процесса фильтрационного переноса и динамики газовыделений через выработанное пространство в периоды снижения атмосферного давления, а также разработана математическая модель разбавления и переноса газовых примесей в вентиляцион-
о
ном потоке горных выработок очистных и подготовительных участков, что позволяет уточнять объем воздуха, необходимого для проветривания.
й* Расчет кол-ва воздуха
- для проветривания подготовительной выработки в периоды экстренны« газовыделений
0.2007
0,1
СО-начальная концентрация кислорода в шахтном воздухе С1 -концентрация кислорода в ВП ' аКа
5- площадь поперечного сечения подготовительной выработки ^-количество "мертвого воздуха", выделяющегося в единицу времени коэф.доставки воздуха
К_р-константа скорости поглощения кислорода |__п.в.-длина воздухоподающей выработки ПДК
Q, і 1с,
Q,
Qs„ + I + KnS„L„. ^ Q^+I + K^L,,
хехр
J Ч(х)*.
= ПДК.
177 мЗ/мин, Т очность=З.Б853457892Є161Е-ОЄ
Расчет кол-во воздуха
для проветривания очистного забоя в период экстренных газовыделений СО-начальная концентрация кислорода в шактном воздухе С1 -концентрация кислорода в ВП
0,2007
0,001
1000
alfa
Б - площадь поперечного сечения подготовительной 1_5г-количество "мертвого воздуха", выделяющегося в единицу времени ^ 1__оч-длина лавы
0,0001 К_р-константа скорости поглощения кислорода
1__п.в.-длина воздухоподающей выработки
ПДК
0,2
QÍCO + IC!
<ї+і+кпа(і^,+ій) (P + KJS
q;+i*k„s,i,„ L)
хехр
= гщк.
197 мЗ/мин, Т очность=3,63312500128998Е-0Є Расчет І
Рис. 2. Форма для расчета количества воздуха
Рис. 3. Результаты расчет количества воздуха для подготовительной выработки
Анализируя результаты, можно сделать вывод о том, что в период экстренных газовыделений потребное количество воздуха в атмосфере подготовительной выработки составляет 130... 420 м/мин. Дополнительное количество воздуха в атмосфере очистного участка составляет 45. 87 % от расчетной величины.
Список литературы
1. Быков Л.Н., Соколов Э.М., Левин Е.М. Состав рудничной атмосферы шахт восточного Донбасса и методы оценки уровня газовыделений
и эффективности проветривания//Уголь Украины. 1967. № 5. С. 45-47.
2. Воронин В.Н. Основы рудничной аэрогазодинамики. М.: Угле-техиздат. 1961. 365 с.
3. Качурин Н.М., Гусев Н.Д. Газовыделение из угольного пласта в подготовительные выработки при разработке глубоких горизонтов // Известия вузов. Горный журнал. №8. 1984. С. 46-50.
4. Касимов О.И., Капиев Р.Э. О точности определения фактического газовыделения на выемочных участках // Вопросы проветривания шахт Донецкого бассейна. 1969. С. 113-122.
5. Пыхачев Г.Б., Исаев Р.Г. Подземная гидравлика. М.: Недра, 1973.
537 с.
6. Скочинский А.А., Лидин Г.Д., Гердов М.А. О явлениях быстрого кислородного обеднения в подземных выработках // Известия АН СССР. ОТН. №11. 1943. С. 251-273.
7. Руководство по проектированию вентиляции угольных шахт. М.: Недра, 1975. 238 с.
8. Соколов Э.М., Качурин Н.М. Углекислый газ в угольных шахтах М.: Недра. 1987. 142 с.
9. Соколов Э.М., Качурин Н.М., Кузнецов А.А. Газовыделение в тупиковые выработки шахт Подмосковного бассейна // Вентиляция шахт и рудников. 1979. С. 72-77.
Качурин Николай Михайлович, д-р техн. наук, проф., зав. кафедрой, ecology@tsu. tula. ru, Россия, Тула, Тульский государственный университет,
Воробьев Сергей Александрович, канд. техн. наук, ген. директор, [email protected], Россия, Москва, Издательский дом «Руда и металлы»,
Афанасьев Олег Александрович, соискатель, [email protected], Россия, Тула, Тульский государственный университет,
Шкуратский Дмитрий Николаевич, ген. директор, director@gallurgy. ru, Россия, Пермь, ОАО «Галлургия»
FORECASTING GAS SITUATION IN COAL MINES DURING PERIOD OF DECREASING A TMOSPHERIC PRESSURE
N.M. Kachurin, S.A. Vorobiev, O.A. Afanasiev, D.N. Shkuratskiy
Mathematical models of diffusion processes in coal layers and the developed spaces in times of emergency gas filtering were considered. The values of flow rate of the gas mixture for the periods of reduction of atmospheric pressure are provided.
Key words: diffusion processes, outgassing, open area, gas filtration.
Nikolay Mikhaylovich Kachurin, doctor of technical sciences, professor, Headof a chair, ecology@tsu. tula. ru, Russia, Tula, Tula State University,
Sergey Alexandrovich Vorobiev, candidate of technical sciences, general director, vorobjov@rudmet. ru Russia, Moscow, The Publishing House “Ore and Metals”,
Oleg Alexandrovich Afanasiev, postgraduate, leader-express@,tula. net, Russia, Tula, Tula State University,
Dmitriy Nikolaevich Shkuratskiy, general director, director@gallurgy. ru, Russia, Perm, Company “Gallurgiy”
УДК 669.187
МАТЕМАТИЧЕСКИЕ МОДЕЛИ МЕТАНОВЫДЕЛЕНИЯ В ПОДГОТОВИТЕЛЬНЫЕ И ОЧИСТНЫЕ ЗАБОИ ИЗ ОТБИТОГО УГЛЯ
Н.М. Качурин, С. А. Воробьев, А.Н. Качурин, И.В. Сарычева
Обоснованы математические модели изменения газоносности отбитого угля на очистных и подготовительных участках угольных шахт. Получены уравнения динамики дегазации отбитого угля в подготовительных и очистных забоях. Приведены результаты вычислительных экспериментов и даны расчетные формулы для инженерной практики прогноза метановыделения при проектировании вентиляции очистных и подготовительных забоев.
Ключевые слова: метан, газоносность, метановыделение, диффузия, отбитый уголь, подготовительный забой, очистной забой, математическая модель.
В подготовительном забое в процессе проведения подготовительных выработок, проводимых частично или полностью по углю, происходит разрушение угля исполнительными органами проходческих комбайнов. Отбитый уголь дробится на блоки, которые можно заменить эквивалентными сферами, радиус которых соответствует некоторой эффективной величине, определяемой гранулометрической кривой, представляющей собой закон распределения разрушенного угля по фракциям различного размера. При очистных работах происходит нарушение структуры угольных пластов [1]. В процессе выемки угольный пласт разрушается исполнительными органами очистных комбайнов. Уголь также дробится на блоки, которые тоже целесообразно заменить эквивалентными сферами. Тогда можно ввести следующие допущения: кусок отбитого угля заменяется эквивалентной сферой, которая дегазируется в диффузионном режиме; движущей силой диффузионного переноса является градиент остаточной газоносно-
