CC BY
20
© М.М. Конорев, Г.Ф. Нестеренко, 2002
УДК 622.458
М.М. Конорев, Г.Ф. Нестеренко
ТЕОРЕТИЧЕСКИЕ АСПЕКТЫ ПРОЦЕССОВ ФОРМИРОВАНИЯ И ПОДЪЕМА ПЫЛЕГАЗОВОГО ОБЛАКА ПРИ МАССОВЫХ ВЗРЫВАХ В КАРЬЕРАХ
п
ри расчетах количества вредных газовых выбросов следует учитывать, что в пылегазовое облако (ПГО) попадает лишь часть взрывных газов, поскольку взорванная горная масса является своеобразным задерживающим и фильтрующим слоем. С достаточной для практики степенью точности пропускная эффективность этого «фильтра» может быть определена по формуле:
^пр
(к р - 1)Уб0
о
р
кУ
1
о
(1)
АУ
^ крУ бл
где АУ - изменение объема взорванного блока, м3; У1 - объем взорванной горной массы, м3; Убл - объем блока, подготовленного к взрыву, м3; кр - коэффициент разрыхления.
По классификации В.В. Ржевского [1], значения кртах = 1,4-1,5 для пород I - II категорий, а для пород Ш - V категорий кртах = 1,2-1,3.
Для определения объема газов, выделяющихся в результате вторичного выброса - с момента развала купола взорванной горной массы до выравнивания давления взрывных газов до атмосферного, необходимо учесть объем газов, оставшихся после первичного выброса:
У = — У
ні 7 у г
(2)
С учетом вероятности проникновения газов находим объем газов, выделяющихся в ПГО при вторичном выбросе:
= (1 - тг) ту'- (3)
К р К р
Суммарный объем газов, проникающих в ПГО, составит:
У = У
У н у И1
+Ун2 = (1 -^-У-.
к 2
(4)
Для пород 1-П категорий суммарный объем газов, выброшенных в ПГО, составляет (49-56)%, а для пород Ш-У категорий (31-41)%.
Начальные параметры ПГО: объем газов (Уг), объем газов, проникающих в ПГО (Ун), начальный температурный перегрев 80Т^ - определяются на основе II начала термодинамики - максимального возрастания энтропии системы «взрывные газы - окружающая среда» при ее переходе в состояние равновесия.
В частности, объем газов Уг определяется по формуле [2]:
т/ М з , Н
Уг = —з1п----------
(5)
рз пт
где Мз - масса взрывчатых веществ (ВВ), кг; рз - плотность ВВ, кг/м3;
РзШ ц
Н =-
относительное давле-
рм>ЯГ о
ние в момент взрыва; Я - универсальная газовая постоянная, дж/кг К; Тц
- температура цепной реакции, К; Т0 - температура воздуха в зоне взрыва, К; пт - коэффициент политропы, соответствующий максимальному возрастанию энтропии системы «взрывные газы-окружающая среда» при переходе в состояние равновесия.
Температура ПГО в момент выравнивания давления взрывных газов до атмосферного определяется по формуле [2]:
Т0 = Тцехр 1пН . (6)
После подстановки значения Уг в формулу 4 определяется объем газов, проникающих в ПГО (V,,), и начальный радиус облака (Rio ):
(7)
После выравнивания давления взрывных газов до атмосферного начинается динамический этап - подъем и развитие ПГО, параметры которого рассчитываются на основе I начала термодинамики и уравнений Мещерского (движение тел с переменной массой) с учетом эффекта смешения. Система уравнений включает [2, 3]:
+ аХі = -уа + а(Го - yz)
dz
dWi
2
dz
- + 2аW 2 = 2р8Г1 3
z + z*
Я і = Я іо + 0^*,
z* = 5Я
о
(8)
(9)
(10)
где уа - сухоадиабатический градиент температуры воздуха, К/м (уа~ 0,01К/м); R^ - начальный радиус ПГО, м; а - коэффициент вовлечения окружающего воздуха в ПГО, м-1; р = g/T - параметр плавучести, м/с2 К; g - ускорение свободного падения, м/с2; Т - средняя температура в атмосфере карьера, К; Wi - скорость подъема ПГО, м/с; у - температурный градиент в атмосфере карьера, К/м; 5оТ = (То - То) - начальный температурный перегрев
ПГО, К; Т0 - температура воздуха в месте взрыва, К.
Однородные линейные уравнения с правой частью системы (8-10) решаются путем замены переменной г на у = 1 +г/г*.
Система уравнений (8-10) приобретает вид:
т
а =
ё±+ъ!± = -у* + А То - г*(у-1)] ёу у У
^ + 6^ = 2^Т,
ё'у у
3
«=-----, ^ ^0 у
г*у
(11)
(12)
(13)
Первое уравнение преобразованной системы решается следующим образом. Примем Т=иу и 3
{ - Уагх +-Т0 - Ггх (у-1)] = б(уЖ
у
ёТ) ёу ёи
тогда —- = и-----------+ V—.
ёу ёу ёу
Уравнение 12 приобретает вид:
.ёу 3уч ёи
и(— + —) + V— = б(у). (14)
ёу у ёу
ёу 3у л
Подберем такое решение, что-----------1---= 0, из которого
ёу у
следует у=у . Значение и определяется по уравнению:
(15)
ёи
у~т = <3( у).
ёу
уЪ„_ Г2-у4 - 3 „4
4 4
и =|^_ёу =—^—у --Г*у + Т0у + уг*у . (16)
Т1 = Т0 - уг*(у-1)] + &Т1 = - Г; у -4уг*у + уг* + с^у 3
(17)
При г=0 у=1 Т = Т0 + 5оТ.
10
Из уравнения 17 С1 = 5оТ + ^Га Г)
Тогда изменение Тi температуры ПГО определится по формуле:
Т1 = с1у - у
-3 „(Га + 3Г)
4
г* + уг* = Т0 - уг*(у-1)] + 5Тi. (18)
Изменение перегрева ПГО в процессе его подъема определится по формуле:
(19)
5Тг = С\у-3 - ■(Г^г*у = Су-3 -Ву
С1 =50Т + (Уа -У)г* =50Т + —,
где 1 0 г 4 0 г 4
аг*
а = у а -у; В1 = -4*.
Из уравнения 19 определится уровень выравнивания температуры ПГО с атмосферной (5Т,=0):
С I 4к
у0 = 4 й = 41 + ~ ’
V В1 V г*
(20)
где А =
<5(>Т/
- максимальный подъем ПГО в политропи-
ческой атмосфере, м.
г I 4к
Поскольку у0 = 1 + —,гт = г*(4/1 +---------------------1). (21)
г* "V г*
Второе уравнение системы решается следующим обра-
зом:
Ц;2 = иу,
2
ёW,■ ёу ёи
----— = и — + V—,
ёу ёу ёу
Q( у) = 2/г*8Т1, (22)
ёу 6у ёи
и(~ +—)+у— = Q(у).
ёу у ёу
В последнем уравнении подберем такое значение у, что
ёу 6у -6
----1---= 0, откуда у = у .
ёу у
Значение и определится по формуле:
и=Iёу = 2/г* I(С1у3 - В1у )у 6'у=у4(С1- В у4)
Значение скорости составит:
(23)
(24)
Ц2 = ^.(с - ^ у4) + С2у-6 2 у 2
Постоянная интегрирования С2 определится из 24 при у=1 Wi=0:
„ /г* ^ В1 ч
С2 = —— (С - 2 ). Поскольку из (20)
С1 = В1 у о,то С 2 а (2 у о -1)
16
Общее решение уравнения скорости подъема ПГО следующее:
Ц2 = ^ [у 4 (2 уо4 - у4) - (2 уо4 -1)];
16 у 6
WI =^[у 4(2уо4 - у4) - (2уо - 1)]а/.
4 у
(25)
(26)
При у=у0, на уровне выравнивания температуры, скорость подъема ПГО имеет максимальное значение:
/Зак
уо
2
/ак
2
4 к 3/ (1 + 4к )/2
г*
(27)
Правая часть уравнения 25 имеет два корня, при которых скорость подъема Wi = 0:
у4 = 1(г = 0) и у2 = 2уо -1 = уЦ = 1 +—(г = гц ).
г*
Уровень конвекции (гц), на котором ПГО останавливается и начинается его рассеяние, можно рассчитать по формуле:
( I—1ТГ А
гц = г*(4у2 -1) = г*
и 8к л
4/1 +----------------1
г*
(28)
Время подъема ПГО до любого уровня определится по формуле:
а
t=
Z dz
J —= z*
W,
I
У=1
dy_
Wi
(29)
После подстановки значения Ші из формулы 26 получим:
t=
1
1
(З0)
^ У=^ У 4 (2 У О4 - У 4 ) - (2 У О4 - 1)
Решение уравнения 30 производится путем преобразования подкоренного выражения и с помощью тригонометрической подстановки. В результате вычитания из левой
части и прибавления к правой части значения у* подкоренное выражение приводится к виду:
(2УоУ 4 - У * - У'8) + (1 - 2у4 + у 8) - (у»4 - 1)2 - (У0 - У V -
к у о4 - 1)2в-(У vУ 422
(У О4 -1)2
(З1)
В уравнении 31 производим тригонометрическую подстановку:
\2
( 4 4 ^
УО - У Уо4 -1
= sin2 p
(З2)
В результате подкоренное выражение 31 приобретает
вид:
(у0 -1)2(1 - sin2 p) = (y0 -1)2 cos2p .
Из З2 определим:
У4 = У0 -(У0 - 1)sinP,
dy 4 = -(y4 - 1)d sin p = -(y4 -1) cos pdp.
(ЗЗ)
(34)
(35)
При подстановке 33 и 35 в 30 (У0 _1) и cosф взаимно уничтожаются, и в результате получим:
ґ = -^ [ dф =-^
Значение ф определится из уравнения 32:
p = arcsin
4 4 4 4
Уо - У _ Уо - У
Уо4-1 Уо4-1
(З6)
(З7)
После подстановки значения ф в 36 окончательно полу-
( У 4 - У о4)
t = -
V^P (yо -1)
(З8)
Уравнение 38 можно упростить, если учесть, что
4 4h
у0 = 1 + :
t=
z* 4 4
(У - Уо)
4h-jaP
(З9)
По формуле 39 можно определить время подъема ПГО до любого уровня. В частности, время подъема ПГО до уровня выравнивания температур составит:
1 (40)
tz_ =
Время подъема ПГО до предельного уровня конвекции составляет:
2
№ =^0ф (41)
По приведенным формулам произведены расчеты при следующих условиях:
Тц = 3000 К; Р3 = 750 кг/м3; кр = 1,3; То = 240 К; у =-
0,021 К/м; R = 287 дж/кгК; р = g/T = 9,8/243 = 0,04; пт = 1,46.
В результате расчетов определены: начальный температурный перегрев ПГО 50Т,- = 14 К, время подъема облака до уровня выравнивания температуры ґгт =28,25 с и до уровня
конвекции ґ = 56,5 с.
Результаты расчетов параметров ПГО при массовых взрывах различной мощности приведены в таблице. ПАРАМЕТРЫ ПГО ПРИ МАССОВЫХ ВЗРЫВАХ.
Мз, т 100 400 600 800 1000
V^ мМ0З 72,0 288,0 4З2,0 576,0 720,0
З 0 СП Зм 28,8 115,2 172,8 2З0,4 288,0
R. , м '0 19,0 З0,0 З4,6 З8,2 41,1
z*, м 95,0 150,0 17З,0 191,0 205,5
zx, м 106,0 1З5,0 145,0 152,5 158,0
zw, м 142,4 186,0 201,0 212,5 221,1
Wimax, м/с 1,68 2,З1 2,55 2,7З 2,92
R z м 40,2 57,0 6З,6 68,7 82,2
м z R 47,5 67,2 74,6 80,7 85,З
Аналогичные расчеты можно выполнить при других условиях: у„>у > -0,05 К/м; £р=1,2-1,5; Т0 = (240-300) К; рз = (750-1000) кг/м3, Тц = (2500-3500) К.
К примеру, при прочих равных условиях, принятых при расчетах в таблице, но при изотермии в атмосфере карьера (У = 0), уровень выравнивания температур (гг) и конвекции ^1^) увеличиваются на (55-65)%, время достижения облаком этих уровней возрастает соответственно до 50 с и 100 с, объем ПГО на этих уровнях увеличивается в 2,2 раза, по сравнению с инверсной стратификацией в атмосфере.
------------ СПИСОК ЛИТЕРАТУРЫ
4
1
1. Ржевский В.В. Процессы открытых горных работ.- М.: «Недра», 1980.-520 с.
2. Конорев ММ., Нестеренко Г.Ф. Вентиляция и пылегазоподавление в атмосфере карьеров.- Екатеринбург: ИГД УрО РАН, 2000. - 312 с.
3. Андреев В., Панчев С. Динамика атмосферных термиков. -
Л.:Гидрометиоздат: 1975. - 152 с.
КОРОТКО ОБ АВТОРАХ --------------------------------------------------------------------------------------------------
Конорев Михаил Максимович - доктор технических наук, зав. лабораторией экологии горного производства (ЭГП), Институт горного дела УрО РАН, г.Екатеринбург.
Нестеренко Геннадий Филиппович - кандидат технических наук, ст. научный сотрудник лаборатории ЭГП, Институт горного дела УрО РАН, г. Екатеринбург.
© Б.Ф. Кирин, А.В. Светличный, 2002
УДК 622. 82
Б.Ф. Кирин, А.В. Светличный
К ВОПРОСУ СОЗДАНИЯ АВТОМАТИЧЕСКОЙ СИСТЕМЫ КОНТРОЛЯ КОЛИЧЕСТВА ОСЕВШЕЙ ВЗРЫВООПАСНОЙ УГОЛЬНОЙ ПЫЛИ И ПРИНЯТИЕ РЕШЕНИЙ ПО ПРЕДУПРЕЖДЕНИЮ ВЗРЫВОВ
В
угольных шахтах взрывы газа метана всегда сопровождаются взрывами пыли, охватывающими большие объемы подземного пространства, приводящими к значительным разрушениям и гибели людей.
Угольная пыль взрывается при условии, если ее концентрация колеблется в пределах 5-600 г/м . При начале взрыва газа ударной волной поднимается во взвешенное состояние осевшая пыль, что и создает взрывоопасную концентрацию.
Поэтому очень важно иметь систему, позволяющую своевременно определить количество осевшей пыли и принять меры по ее связыванию или переводу в невзрывоопасное состояние.
На основании наших исследований и исследований других авторов изменение концентрации пыли по длине выработки подчиняется экспоненциальному закону (рис. 1, 2, 3).
Как видно из рисунков концентрация пыли по длине выработки изменяется быстро на первых 50 м от источников пылеобразования. За указанным расстоянием изменение концентрации пыли в сторону уменьшения происходит медленно и на расстоянии 100-120 м она практически остается постоянной.
Закономерность изменения концентрации пыли при изменении скорости движения воздуха в горной выработке остается прежней. Изменяются только количественные значения. При большей скорости движения воздуха абсолютные значения концентраций пыли выше, чем при меньших значениях. Это объясняется тем, что подъемные аэродинами-
ческие силы выше. Они обеспечивают поддержание более крупных частиц пыли во взвешенном состоянии на больших расстояниях от источника пылеобразования. Исходя из экспоненциального закона изменения концентрации пыли по длине выработки предложен ряд формул для определения возможной концентрации пыли на любом расстоянии от источника ее образования при известном первоначальном значении концентрации пыли у источника образования.
Так, по формуле (1) можно определить концентрацию пыли в любом месте выработки, если известны коэффициенты а и в, что в реальных условиях представляет определенную сложность.
v2
2
С = См • а • е ^ , мг/м3 (1)
где См - максимальная концентрация в сечении выработки, мг/м ; а, в - коэффициенты, учитывающие аэродинамические характеристики выработки, концентрацию и свойства пыли соответственно; е - основание натурального логарифма; Y - расстояние от места измерения концентрации до источника пылеобразования, м; R - расстояние от оси выработки до ее бока, м.
Формула (2)
-є
K
1
(2)
K о a + bL
где К — концентрация пыли на расстоянии L от источника пылеобразования, мг/м ; Ко — начальная концентрация пыли у источника пылеобразования, мг/м ; a и b — коэффициенты: a = 1, b = 0,03; L — расстояние от источника пылеобразования до места замера, м.
Формула (3) представлена в виде
Св = Cо • е ti , мг/м3
(3)
где Св — концентрация пыли в выбранном нами месте выработки, мг/м ; Со — первоначальная концентрация пыли в месте пылеобразования, мг/м ; е — основание натурального логарифма; t — коэффициент затухания концентрации, с/м; і - расстояние от места пылеобразова-ния до места, где мы должны знать концентрацию пыли, м.
JS t = m 4 — ,
\Q
(4)
где т - коэффициент, учитывающий оседание частиц пыли; S - площадь поперечного сечения выработки, м2; Q - количество проходящего воздуха, м/с.
Значения т можно принимать для подготовительных выработок 0,052, для очистных забоев 0,064.
Концентрация пыли в сечении выработки значительно изменяется. Поэтому определение концентрации пыли в сечении должно производиться с учетом коэффициента поля концентраций Кс
Кс =
ao • H • е aoh
1 - е
-VoH
(5)
где #0 - —— - коэффициент, учитывающий оседание Кг
частиц пыли в турбулентном потоке, 1/м; V0 - скорость
оседания частиц, м/с; К - коэффициент турбулентной диффузии, с/м2; Н - высота выработки, м; И - высота установки датчика от почвы выработки.
Значения V0 принимаются осредненно в предположении, что на расстоянии 100 м оседают частицы пыли размером > 10 мкм. Тогда для угольной пыли V0=2,56
см/с = 0,0256 м/с. Коэффициент турбулентной диффузии 0 75
Кт - 0,044-V-Re 5 . Где V - кинематическая вяз-
Рис.1. Изменение концентрации по длине выработки а - и=0,3 м/с; Ь
- и=0,4 м/с; с - и=1 м/с; 1 - экспериментальные данные; 2 - теоретические данные
Рис. 2 Отложение взрывоопасного количества пыли происходит в вентиляционном штреке на длине 30-150 м: 1 - шахта “Октябрьская, 6-й вентиляционный штрек; 2 - шахта №1-1бис, 7-й восточный вентиляционный штрек; 3 - шахта №19-20, 8-й вентиляционный штрек участка №12; 4 - шахта №16-17, восточный вентиляционный штрек; 5 - шахта им. Карла Маркса, вентиляционный штрек участка №22; 6
- шахта №13-бис, 1-й восточный вентиляционный штрек; 7 - шахта №17-17-бис, 10-й восточный вентиляционный штрек
Рис. 3. Изменение концентрации пыли по длине подготовительной выработки
кость воздуха, м2/с, V = 14,4 -10 6 м2/с; Re — число Рейнольдса
Re
U ср • Dr
(6)
где иср — средняя скорость движения воздуха в выработке, м/с; Dr — гидравлический диаметр сечения выработки.
Если принять, что скорость оседания частиц в сечении остается постоянной, то концентрация пыли по высоте от почвы выработки может быть определена по формуле
VSY
Св = Со • е Kt , (7)
где Vsy — скорость оседания частиц пыли, м/с; Kt — коэффициент диффузии.
Отношение VSY/Kt для угольной пыли размером 20 мкм составляет 4,2, размером 10 мкм — 0,16, размером 1 мкм — 0,002.
При кажущейся простоте формулы (7) очень трудно подсчитать отношение VsY/Kt для всего ряда дисперсности пыли.
Если принять, что на участке L = 100-120 м оседают фракции пыли более 10 мкм и предположить с определенной долей погрешности, что скорость оседания частиц может быть принята каким-то средним значением, то изменение концентрации пыли по длине выработки (вентиляционный штрек) может быть определена по формуле
-2VSY -і
C = Co • е k'r'Ub , мг/м3 (8)
где Vsy - скорость оседания частиц, м/с; R - гидравлический радиус выработки, м; U - скорость движения воздуха по выработке, м/с; і - расстояние от места пылеобразования до пункта замера, м.
Произведя расчеты по рассмотренным формулам приходим к выводу, что они дают значительный разброс. Расчетные данные от данных измеренных в шахтах отличаются от 15 до 30 %. Наиболее близкие результаты получены при расчете по формуле (8).
Экспериментальные данные показывают, что 80 % пыли оседают на расстоянии 150 м от источника пылеобразования. Это видно из приведенных ниже значений.
V
С учетом изложенного константируем, что датчик измерения концентрации пыли следует устанавливать на расстоянии 150 м от источника пылеобразования.
Распределение пыли, выносимой из лавы на вентиляционный штрек
на участке штрека...............от 10 до 20 м - 25%
на участке штрека...............от 20 до 40 м - 22%
на участке штрека...............от 40 до 60 м - 11%
на участке штрека...............от 60 до 80 м - 7%
на участке штрека...............от 80 до 100 м - 5%
на участке штрека...............от 100 до 120 м - 6%
на участке штрека...............от 150 до 200 м - 3%
на участке штрека...............от 200 до 300 м - 1%
Принимая в дальнейших рассуждениях формулу (8). Определяем значение концентрации пыли на расстоянии 150 м от источника пылеобразования.
Для автоматизированной системы определения количества осевшей пыли поступаем следующим образом.
1. Измеряем прибором концентрацию пыли на
расстоянии 150 м, где установлен стационарный датчик.
2. Обратным счетом в соответствии с формулой (8) по программе определяем концентрацию пыли на расстоянии 10, 50, 100 и 150 м от источника пылеобразования.
3. Разность значений концентраций, измеренных на расстоянии 10, 50, 100 и 150 м, позволяет нам за определенное время и знании количества проходящего воздуха определить массу пыли осевшей на каждом участке за этот промежуток времени.
4. Зная параметры выработки - размеры площади поперечного сечения выработки по, периметру ее рассчитываем толщину слоя осевшей пыли.
5. При известном объеме выработки и количестве осевшей пыли определяется концентрация пыли в данном объеме.
6. Сравниваем полученное значение с критическим с учетом газового фактора и, если полученные значения превосходит критические или близкие к ним, то система выдает соответствующую рекомендацию.
КОРОТКО ОБ АВТОРАХ
Кирин Б.Ф., Светличный А.В. — Московский государственный горный университет.
