Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Рудничная аэрология.-1

.pdf
Скачиваний:
8
Добавлен:
20.11.2023
Размер:
16.52 Mб
Скачать

Пусть имеется запись изменения концентрации газа при стацио­ нарном газодинамическом процессе. Определим для этого процесса все отклонения Ас измеренных (мгновенных) значений концентрации от математического ожидания концентрации и разделим их на группы, или разряды, в зависимости от величин Ас, например Асх = (0 4- 4-0,1)%; Ас2 = (0,1 ч-0,2)% и т. д. Пусть число отклонении Дсх,

цопавших в первый разряд, будет пх, во второй — тг2 и т. д. Их сумма,

т

очевидно, равна числу измерений : N 2 ty(i номер разряда, i=i

т — число разрядов). Теперь можно построить ступенчатый гра­

фик, на оси

абсцисс которого

 

 

п

 

отложим АС/, а на оси ординат —

 

 

NÛC,

 

- ^ - (р и с . 142).

Такой

график

 

 

 

 

 

 

 

называется

г и с т о г р а м ­

 

 

 

 

мо й .

Очевидно, что суммарная

 

 

 

 

площадь

всех

 

ее

прямоуголь­

 

 

 

 

ников равна 1. Чем уже будут

 

 

 

 

прямоугольники,

тем

ближе

 

 

 

 

будет

ступенчатая линия при­

 

 

 

 

ближаться к некоторой кривой,

 

 

 

 

с которой она сольется в пре­

 

 

 

 

деле

при

Дс/

 

0.

Эта кривая

 

 

 

 

называется

 

г р а ф и к о м

 

 

 

 

п л о т н о с т и

 

в е р о я т ­

 

 

 

 

н о с т и

отклонений Ас.

Она

 

-ДС5 -âCi+ 'ACfACfdCj 0 ACf Дсг ACjAC^. Ac5

графически

выражает

 

закон

Рис. 142. Гистограмма

распределения

распределения

Ас.

 

 

 

Основные процессы

метано-

отклонений

измеренных

значений кон­

центрации от ее математического ожи­

выделения в

угольных

шахтах

 

 

дания

 

с достаточной

 

точностью

под­

з а к о н у ,

описываемому к р и *

чиняются

н о р м а л ь н о м у

в о й

Г а у с с а , изображенной на

рис 142.

 

Уравнение этой кривой имеет вид

 

 

 

 

 

 

 

 

/ =

1■ =- е ~ К 4 ) ’

 

(XIV,61)

 

 

 

 

 

 

 

• V 2л.

 

 

 

 

где / — плотность

вероятности. При большом числе наблюдений

 

 

 

 

 

 

 

/ — lim

N Ас

 

 

 

 

 

 

 

 

 

Дс->*0

 

 

 

а — с р е д н е к в а д р а т и ч е с к о е о т к л о н е н и е ,

J

(XIV,62)

Y N

N — общее число измерений; к — номер измерения.

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

максимум случайного процесса практически весьма

трудно.

В по­

добных случаях расчет ведут по с т а т и с т и ч е с к о м у

м а к ­

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

Пусть qcр — средняя величина газовыделения в выработку;

Ятах — его статистическая

максимальная

величина.

Последнюю можно представить в виде

 

 

Îma* =

?cp + Ag,

(X IV ,63)

где Дд — такое превышение величины

газовыделения над средним

ее значением, при котором определенная по формуле

(XIV,63) величина gmax с некоторой

вероятностью будет

больше любого измеренного значения газовыделения.

Коэффициент

 

 

 

 

=

=

Яср

(XIV,64)

 

9ср

 

называется

к о э ф ф и ц и е н т о м

н е р а в н о м е р н о с т и

г а з о в ы д е л е н и я .

 

через среднеквадратическое

Величину

Дq можно определить

отклонение газовыделения

N

 

 

2

(д?*)2

 

 

о =

к=1

N

(XIV,65)

 

 

 

 

 

 

где N — общее

число измерений;

 

 

к — номер

измерения;

 

 

 

ДqK— отклонение измеренной величины газовыделения от его среднего значения дср.

Физически а определяет некоторый средний уровень отклонения фактического газовыделения от его среднего значения (рис. 143). При этом значительная часть фактически замеренных значений q будет превышать уровень qc? -f- а. Однако если этот уровень под­ нять до дср + 2а, число таких превышений станет значительно меньше. В теории вероятностей доказывается, что для нормально распределенной случайной величины (какой является метановыде­

ление) число отклонений, превышающих уровень qcp + 3а, состав­ ляет примерно 0,15% от всех замеренных значений газовыделения. Следовательно, вероятность того, что

<7тах = <7ср + 3а

(XIV,66)

будет больше любой замеренной величины газовыделения, равна 0,9985, или 99,85%. Определение qmax со столь высокой вероятностью достаточно для надежного расчета количества воздуха. Формула (XIV,66) выражает п р а в и л о т р е х с и г м а .

Из уравнений (XIV,64) и (XIV,66) имеем

ки= 1 + 3 — .

(XIV,67)

Рис. 143. График га­ зовыделения в выра­ ботку

По данной формуле определяются коэффициенты неравномер­ ности газовыделения в шахтах, которые в зависимости от условий могут изменяться в основном от 1,2 до 2,5.

Ввиду того что пульсации газовыделения Aqyвызываемые отдель­ ными факторами, примерно одинаковы при различных средних его значениях дср, значение Кнуменьшается с ростом дср (см. рис. 8). Зависимость Ки (qcp) различна для различных условий.

§ 89. Методика определения характеристик газодинамических процессов

Определение числа наблюдений. Из математической статистики известно, что среднеквадратическая погрешность оа при определе­ нии среднеквадратического значения некоторой величины (напри­ мер, концентрации сср) подсчитывается через среднеквадратическую ошибку одного измерения о и число измерений N по формуле

о

оА (XIV,68) 7 Г *

Величина а обычно дается в паспорте прибора. Например, для современных анализаторов метана о = ±(0,1 —0,3)% СН4. По­ грешность оа определяется из условий требуемой точности наблюде­ ния. Так, если требуется концентрацию сср= 1% измерить с точ-

ностыо до 0,1 ее величины, то OA = -JQ сср = 0,1%. Таким образом,

зная or и аА, по уравнению (XIV,68) легко определить требуемое число наблюдений N.

На рис. 144 изображена зависимость N (сср) для погрешностей прибора о = ± 0 ,1 ; 0,2 и 0,3% при аА = ± 10% от сср. Из рисунка следует, что для измерения малых абсолютных значений величины

Рис. 144. Определение числа наблюдений в зависимости от абсолютной величины из­ меряемой концентрации газа

большой

погрешностью

тре­

буется

производить

 

значи­

тельное

число

 

измерений.

Для сокращения

числа наб­

людений

 

следует

пользо­

ваться более точными прибо­

рами. Так,

при

 

сср =

 

0,25%

для

определения

этой

вели­

чины с

погрешностью

оА =

=

0,1сср при а =

± 0,3%

по­

требное

число

 

измерений

N

=

60,

а при

а =

±0,1%

N

=

14.

 

 

 

 

 

 

 

 

Определение средней вели­

чины газовыделения.

 

Сред­

няя величина газовыделения

в

выработке

определяется

как

разность средних

деби­

тов

газа

в

вентиляционном

потоке в ее конце qcp%к и в

начале

дср#н:

9ср=

ЯсР. к“ Яср. н* (XIV,69)

В свою очередь, для определения дср; к ц Яср.н замеряются дебит воздуха в начале QH и в конце QK выработки, средние концентрации газа сср н и сср к, а затем вычисления ведутся по формулам:

сср. к Яср. к = <?к “ Ш Г

S (XIV.70)

сср. н 9ср. н = <?н 100

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

Однако в шахтных условиях они обычно изменяются во времени, что вносит погрешность в измерение дср. Например, если при одно­ временных измерениях в начале и конце выработки в ее начале появилось газовое облако (от взрывных работ, работы комбайна и т. п.), это увеличит сср н, что уменьшит qcp; если в момент замера

облако будет в конце выработки, это вызовет увеличение сср к и qcp. Чтобы избежать влияния нестационарностн газовыделения, замеры сср. н 11 Сср. к надо производить в одном и том же объеме воздуха, т. е.

со с д в и г о м

в о в р е м е н и на величину

 

 

t =

(XIV,71)

 

иср

 

где L — расстояние между замерными пунктами;

 

иср — средняя

скорость движения воздуха между ними.

то, что

Сдвигать время замеров целесообразно несмотря на

нестационарность

Q требует производства одновременных

замеров.

Объясняется это тем, что степень нестационарностн расхода воздуха меньше, чем газовыделения.

Определение коэффициента неравномерности газовыделения. Для

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

Согласно (XIV,67) Ки зависит от средней величины газовыделе- ] пя qcр и среднеквадратической величины пульсаций газовыделения а. Определение qcp было рассмотрено выше. Остановимся на опре­ делении а.

Из выражения (XIV,65) следует, что абсолютная величина а рассчитывается по абсолютной величине отклонений отдельных измеренных значений газовыделения от его средней величины. В свою очередь, эти отклонения зависят от естественных пульсаций газо­ выделения и ошибок измерения (погрешность прибора, недостаточ­ ная квалификация наблюдателя, несовершенство методики заме­ ров); последние увеличивают отклонения и их желательно макси­ мально снизить.

Если измеренное значение газовыделения в выработке больше среднего, установить, насколько оно ошибочно, без специального анализа невозможно. Однако если газовыделение q меньше сред­ него, можно утверждать, что оно ошибочно, по крайней мере тогда, когда оно меньше нуля (отрицательно), ибо газовыделение по физи­ ческому смыслу существенно положительно. Иными словами, Aq — = 7ср — 7ср- В то же время в области газовыделений больших, чем qcр, часто А? = q —qcp> qc? вследствие больших ошибок измерения.

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

существенно ошибочных

значений q

qcp. Исследования показы­

вают, что такая выборка

не уменьшает надежности о вследствие

уменьшения общего числа замеров, связанного с отбрасыванием значений q > gcp.

 

Исключить ошибку измерения можно также, пользуясь формулой

где' а — среднеквадратическое

отклонение

измеряемой

величины

 

газовыделения от ее

среднего значения дср, включающее

 

и ошибку измерения;

 

 

 

 

°и ~~ среднеквадратическая ошибка, зависящая от ошибки изме­

 

рения концентрации газа и количества воздуха.

 

Г л а в а XV

 

 

 

ДИНАМИКА АЭРОЗОЛЕЙ В ГОРНЫХ ВЫРАБОТКАХ

 

§ 90. Общие сведения о динамике аэрозолей

 

 

Под динамикой аэрозолей понимают движение частиц аэрозоля,

их

гравитационное, электростатическое и

инерционное

осаждение,

а

также процесс перехода порошкообразных тел в

аэрозольное

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

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

§ 91. Уравнение движения для частиц аэрозоля

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

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

Пусть в неподвижном воздухе со скоростью v перемещается

частица с массой т . Если она движется с ускорением

то испыты­

вает лобовое сопротивление R, вызванное инерцией среды, приводи­ мой в движение, и пропорциональное ускорению. Сила, приложен­ ная к частице, может быть представлена в виде

F =

+

(XV.1)

где т г — масса среды, захваченной движением (с ускорением);

т ! ~~ ло®овое сопротивление.

Чтобы найти величину т \ рассмотрим элемент среды размером dxdydz в системе координат, связанной с телом. Масса этого объема рQdx dy dz, а кинетическая энергия всей среды, приведенной в дви­ жение,

 

 

 

Т =

dxdydz,

(XV,2)

где

Ро — плотность

среды (р0 = const);

 

 

и — скорость движения элемента среды.

 

 

Интегрирование проводится по всему объему среды. Энергию

выразим через скорость движущегося тела:

 

И Л И

 

Т = -£? Ч

( т ) 2 dxdydz

(XV.3)

 

 

 

 

 

 

 

Т =

0 - 2 ^ - ,

(XV,4)

где

Q — интеграл,

имеющий размерность объема.

 

 

Кинетическую энергию приведенной в движение среды можно

представить в

виде

 

 

 

 

 

 

T = - ^ - v K

(XV,5)

Приравнивая

выражение (XV,5) к (XV,4), получаем

 

 

 

 

m' = Qp0.

(XV,6)

Масса т ! зависит от формы и размеров частицы, направления ее движения и от вязкости среды. Тогда лобовое сопротивление тела выразим формулой

Д = а р о ^ .

(XV,7)

Присоединенную массу среды можно рассматривать как кажу­ щееся увеличение массы тела, если оно движется с ускорением в идеальной среде (без трения). Для движущегося шара радиусом г

д==Т * Т яг3’

(xv,8)

т. e. объем присоединенной массы равен половине объема шара. Для цилиндра длиной /, движущегося в направлении нормали к его боковой поверхности,

Q = пгЧ,

(XV,9)

т. е. объем присоединенной массы равен объему цилиндра. В случае эллипсоида, движущегося перпендикулярно к его большой оси вращения, Q имеет промежуточное значение между половиной и всем объемом фигуры. Если частица полностью увлекается потоком (со скоростью и), то сила, действующая на нее, равна силе, приложен­ ной к объему среды V вытесняемой частицей (в отсутствие послед­ ней). Количество движения среды в этом случае будет

m'u = p0Vu,

(XV,10)

а сила

 

/Г1 = 4 ( Р о ^ ) = Ро^-^-.

(XV,И)

Если частица не полностью увлекается средой, т. е. получает неполное количество движения, то среда приобретает дополнитель­ ное количество движения, равное т ' (и — и), а на частицу действует дополнительная сила, равная

F 2 = m' - ^ ( u - v ) ( т ' ^ т ) .

(XV,12)

Полная сила, действующая на тело, движущееся в среде,

F = P o V ^ + m '- lr {u - v ).

(XV,13)

Выразив F через т получим уравнение движения частицы в среде

<x v ’14>

Обозначив плотность частицы через р', уравнение (XV,14) можем переписать в следующем виде:

Ро-£+ Р'|=(Ро + Р')1Г.

(XV,15)

Умножим обе части этого равенства на dt и проинтегрируем по времени, тогда

 

(Ро+ Р > =

(Р о + Р > + *.

(XV, 16)

где

к — некоторая постоянная

величина.

0. При р' =

=

Если скорость и = 0, то v =

0 и, следовательно, к =

р0 имеем и = и. Это условие выполнимо при одинаковых удель­

ных весах частиц и среды. Выражение (XV,16), так же как и (XV,14),

2Г>8

является уравнением^движения частицы в потоке. Если среда дви­

жется с постоянной скоростью, то член р0У ^

= 0. Тогда при нали­

чии внешней силы F

(t)

можно написать

 

т

чь

= F (*>+ т ' 4 <“ - * ) •

(XV,17)

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

т'-^ -(и — и) = 6лгг) (ц— у),

(XV, 18)

где г| — вязкость среды.

Для частиц другой формы в последний член вводят поправочный

коэффициент ф'. Тогда уравнение (XV,17) принимает вид

 

m -^ = F (О + блГПф' (u —v),

 

(XV/19)

где г — линейный размер

частицы.

 

 

 

§ 92. Исследования

характера распространения пыли

 

в выработках по направлениям х, у, z

 

 

Вначале рассмотрим закон изменения концентрации пыли в атмо­

сфере по высоте. Количество пыли, проходящее в

1 сек

под

дей­

ствием силы тяжести сверху

вниз через горизонтальную

площадку

в 1 см2, равно usc2 (где с2 — концентрация пыли

вдоль

оси

Oz).

Число частиц, пролетающих через эту площадку

в обратном

на­

правлении благодаря турбулентной диффузии, равно DT^

(где DT

коэффициент турбулентной диффузии). Приравнивая эти выражения и разделяя переменные, получаем

dc2

(XV,20)

С2

Так как зависимость Dr от у вблизи почвы сильная, а вдали -- слабая, то в центральной области выработки считаем коэффициент турбулентной диффузии постоянным и, интегрируя уравнение (XV,20), получаем

lnc2= - ^ y + k .

(XV,21)

Если у =

0, то k = In cZo (где cz%— концентрация пыли в атмосфере

у почвы

выработки). Тогда выражение (XV,21) имеет вид

cÿ= c2#e x p ( —

( XV, 22)

Наблюдения показали, что для частиц с г =

1 мк концентрация

по вертикали практически не меняется. Для полидисперсного аэро­ золя (2—60 мк) наблюдается неравномерность распределения кон­ центрации по вертикали.

с,%

80

60

40

20

О

1

2

3

Ч

5

6

х,м

Рис. 145. Относительное падение концентрации частиц по длине при различных значениях vsi R , ис?:

Кривая

 

VS,

R, Л1

иср ,

Кривая

 

Vs,

R,

 

м/сек

 

м/сек

 

м/сек

 

1

7 10-*

0,8

0,1

7

7 - 10-3

0,6

2

7

10-4

1,0

10,0

8

4

10—7

1,5

3

7

10“ 4

1,5

/|,0

0

7

10-4

1,5

4

 

0

1,5

0,1

10

7-10-4

1,3

5

7 - 10-4

0,6

/i,0

11

7 * 10-4

1,0

6

7«10-*

0,6

10,0

12

4 - 10-7

1,5

1Хср ,

м/сек

1,0

1,0

0,1

0,1 о оо

Рассмотрим теперь изменение концентрации пыли по длине выработки х, характеризующее отложение пыли в последней. Урав­ нения, описывающие движение частиц, дают возможность в прин­ ципе вычислить концентрации в различных сечениях выработки. Однако получить из этих уравнений распределение концентраций по сечению выработки на расстоянии х от забоя и среднюю концен­ трацию в сечении довольно сложно. Значения сср в сечениях можно вычислять для различных параметров потока — его средней ско­ рости иср, гидравлического радиуса R и различных однородных по составу аэрозолей со скоростями осаждения vB с помощью траекто­ рий частиц, полученных на счетных машинах. Графики изменения концентрации частиц по длине выработки при различных значениях vs, Л , иср изображены на рис. 145. Из графиков следует, что относи-