книги / Рудничная аэрология.-1
.pdfПусть имеется запись изменения концентрации газа при стацио нарном газодинамическом процессе. Определим для этого процесса все отклонения Ас измеренных (мгновенных) значений концентрации от математического ожидания концентрации и разделим их на группы, или разряды, в зависимости от величин Ас, например Асх = (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. Из графиков следует, что относи-