Аэрогазодинамические процессы на очистных участках шахт и рудников

УДК 622.4/.6

АЭРОГАЗОДИНАМИЧЕСКИЕ ПРОЦЕССЫ НА ОЧИСТНЫХ УЧАСТКАХ ШАХТ И РУДНИКОВ

Н.М. Качурин, И.И. Мохначук, А.А. Поздеев, Г.В. Стась

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

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

Очистные участки шахт и рудников могут конструктивно сильно отличаться, но при этом структурно практически полностью быть идентичными. Структурные элементы любого очистного участка угольной шахты, рудника по добыче руд различных металлов или строительных материалов, а также рудника по добыче калийной руды всегда включают следующие горно-геологические и геотехнологические объекты: очистной забой; выработанное пространство; выработки, примыкающие к очистному забою; поверхности обнажения полезного ископаемого (угольного пласта, рудного тела, калийного пласта и т.п.); отбитую горную массу полезного ископаемого различного гранулометрического состава; воздушные потоки в системе горных выработок очистного участка. Между перечисленными структурными элементами очистного участка существует аэрогазодинамическая связь (рис.1).

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

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

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

Рис. 1. Аэрогазодинамическая структурная схема очистного участка

шахт и рудников:-- свежая струя воздуха;

---- - исходящая вентиляционная струя;

^ - конвективно-турбулентные диффузионные потоки

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

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

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

Если происходит взрывная отбойка полезного ископаемого, процесс конвективно-турбулентной диффузии газов ВВ начинается сразу же после взрыва. Газовые смеси, выделяющиеся из выработанного пространства (это, в первую очередь, С02 и N2 , составляющие основу состава так называемого «мертвого воздуха»), разбавляют кислород в атмосфере очистного участка, и процесс разбавления может быть более значимым по сравнению с процессом поглощения кислорода горным массивом.

Тогда математическое описание предлагаемой физической модели аэрогазодинамики очистного участка шахты или рудника можно представить в следующем виде [1 - 3]:

- фильтрационный перенос газов в пористой сорбирующей среде

еВ 2

а, ч а2

—( шр) + тр) =

ау 7 а п

+

V , °

Л1+

где т - пористость горного массива, пористыми сорбирующими горными породами; р - плотность газа; е - временной масштаб корреляции, имеющий размерность времени; В - норма корреляционного тензора в начальный момент времени, характеризующая газовую проницаемость горного массива; Кр- коэффициент, характеризующий скорость газообмена между твердой фазой и свободным объёмом пор; V- главный вектор скорости фильтрации газа; 0(х1у х2, х3, I) - функция, учитывающая влияние внутренних источников (или стоков) при фильтрации газа в рассматриваемом горном массиве; х1г х2, х3- пространственные координаты; ? - время;

- фильтрационный перенос газов в трещиновато-пористой среде

ар а ^ а2р Л а2р ( л

' ах2 1=1 с^2

2

а1 а11=1 ах2 1=1 ах2

1 12

где ф=кТ(ртф)~ ;ц=к6~ кТ! ; р - давление газа в трещинах; кт , кб- газовая проницаемость трещин и породных блоков соответственно; т0 - пористость породных блоков, слагающих трещиновато-пористую среду; ¡и - динамическая вязкость газа; ¡5- коэффициент сжимаемости метана; I - среднее значение характерного размера пористых породных блоков; О!(х1, х2, х3, I) - функция, учитывающая влияние внутренних источников (или стоков) при фильтрации газа в рассматриваемом горном массиве нарушенной структуры;

- диффузионный перенос газов в пористой сорбирующей среде

| + ^У (оУ) = ё!У [(+ + ) вгаё с ] + д, да = к д [а р (с)" а ],

а

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

- дегазация отбитого полезного ископаемого

л

- КРх,

- = оэ

а э

д х 2 дх

—г +—

чдг г дгу

где .х - газоносность рассматриваемых кусков горной массы полезного ископаемого; ОЭ - эффективный коэффициент диффузии газа в моделируемых кусках горной массы; Кр - константа скорости распада диффундирующего газа (процесс имеет место при диффузии радиоактивных газов, в остальных случаях Кр= 0); г - текущий радиус сферической системы координат для изотропной однородной сферы, моделирующей кусок отбитой горной массы;

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

— + ё1у

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

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

Адаптированные уравнения фильтрационно-диффузионного переноса газов в горном массиве и вентиляционных струях воздуха, как прави-

>

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

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

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

к

Jf =— ,

у и

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

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

3в = — DgradCi,

где - вектор диффузионного потока газа; Э - коэффициент диффузии, характеризующий диффузионную проницаемость рассматриваемого горного массива; ¡и - динамическая вязкость газа; Сг концентрация /-го газа.

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

Таблица 1

Базовые математические модели фильтрационно-диффузионного переноса газов в горном массиве и отбитой горной массе

№ п/п Вид аэрогазодинамического процесса Адаптированное уравнение, описывающее аэрогазодинамический процесс Начальные и граничные условия

1 Выделение метана с поверхности обнажения угольного пласт ар2 а у а у а г а2 ^ ах2' где ^ - период релаксации; ху.п. - пьезопроводность угольного пласта p(x, 0) = p0 =const, а / at [p(x,0)] = о, p(0,t) = pc = const , lim p ф ro x—ro

2 Газовыделение с поверхности обнажения калийных пластов ар2 а2р2 . _ , к = к Л +юexp( рд 1), аг ах2 д где к - пьезопроводность калийной соли; ю, Р - параметры десорбции газа p2(x, 0) = pi = const , p2(0,t) = p2 = const, limp2(x,t) Фго x—ro

3 Метановыделение из подработанных горных пород ар а3р а2р Л 9 =Ф 9 • p(z, 0) = p0 = const p(+0,t) = p0-(p0-p( xexp(-rf *<t), lim p(z,t)=const z—0

4 Метановыделение из куска отбитого угля ах _ (а2х 2 ах ^ — = °Э —г+ — а ^аг г аг) x(r, 0) = x3 = const, x(R,t) = xR(t) ,где R - радиус эквивалентной сферы.

5 Метановыделение из отбитого угля в очистной забой ахо ахо —о+ут —- = -!Д, аг т а^ Д' где хо - газоносность угля на конвейере; ут- скорость транспортирования угля xo&0) = 0, Xo (0,t) = X3, где Хз - остаточная газоносность краевой части угольного пласта

6 Выделение радона с поверхности обнажения горного массива Н2А ё . ХА + ] = 0, Э ёх2 где А - удельная активность газовой смеси; X - постоянная радиоактивного распада; I -интенсивность образования радона в угольном пласте A п=0, lx=0 ' lim A фго x—^ro

7 Поглощение кислорода поверхностью обнажения горного массива аск ^ а2ск и = ^гч о кпС„ ' а Э ах2 п к где ск - концентрация кислорода; к п - константа скорости поглощения кислорода углем. ck (x,0) = 0, cK (0,t ) = c0 = cons lim c„ ф ro к x—ro

)х

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

Адаптация уравнения (5) с учетом таких допущений позволила получить базовые математические модели конвективно-турбулентной диффузии газовых примесей в выработках очистного участка, представленные в табл. 2. Решение уравнений, приведенных в табл. 2, для конкретных начальных и граничных условий позволяет получить теоретические зависимости для оценки газовой опасности очистных участков шахт и рудников. При этом следует отметить тот факт, что имеет место единый системный подход к решению задачи на базе общих математических моделей фильтрационно-диффузионного переноса газов в пористых средах и вентиляционных струях воздуха.

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

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

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

Таблица 2

Базовые математические модели конвективно-турбулентной диффузии газовых примесей в выработках очистного участка

№ п/п Вид аэрогазодинамического процесса Адаптированное уравнение, описывающее аэрогазодинамический процесс Начальные и граничные условия

1 Перенос метана в очистном забое зС зС ^ —+и— = У I, (г), а зх ^ где и - средняя скорость воздуха; I! -интенсивность газовыделения из ь го источника C(x, 0) = C0 = const, C(0,t) = CH = const

2 Перенос метана в выработках очистного участка зС з , ^ ^ з2С , ч а + зХ ( иС) = 0 ^с+11. (г), где и = и(х) C(x, 0) = C0 = const, C(0,t) = CH = const, lim C(x,t) ф ж

3 Перенос радона в пределах очистного участка дА дА --+ и-= -ХА + 21, сЧ дх где 21 - суммарные выделения радона в воздух в выработок очистного участка из различных источников A(x, 0) = A0 = const, A(0,t) = A = const, lim A(x,t) ф ж

4 Динамика концентрации кислорода в пределах очистного участка буро-угольной шахты +;* (^ >=^ - кск, где СК - концентрация кислорода в атмосфере горных выработок; К -константа скорости поглощения кислорода из атмосферы очистного участка CK(x,0) = C0 = const, CK(0,t) = CH = const, lim CK (x,t) ф ж

5 Перенос газовой примеси на очистном участке калийного рудника зС з / ч з2С — +—(иС) = БЭ—- а зх э зх2 C(x, 0) = C0 = const, C(0,t) = CH = const, lim C(x,t) ф ж

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

Рис. 3. Подсистема контроля и управления процессом обеспечения

воздухом очистных участков

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

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

Список литературы

1. Качурин Н.М. Перенос газа в породоугольном массиве // Известия вузов. Горный журнал. 1991. Вып. 1. С. 43 - 47.

2. Аэрогазодинамика углекислотообильных шахт/ Качурин Н.М. [и др.]. М.: МГГУ. 2005. 345 с.

3. Качурин Н.М., Стась Г.В., Качурина О.Н. Математическая модель выделения радона с поверхности обнажения разрабатываемого угольного пласта // Известия ТулГУ. Экология и безопасность жизнедеятельности. Вып. 7. Тула: Изд-во ТулГУ. 2004. С.187 - 190.

Качурин Николай Михайлович, д-р техн. наук, проф., зав. кафедрой, есо1о£\П1Ш1. tula.ru, Россия, Тула, Тульский государственный университет,

Мохначук Иван Иванович, канд. экон. наук, председатель, ,ecology@tsu. tula.ru, Россия, Москва,Российскийнезависимыйпрофсоюз работников угольной промышленности,

Поздеев Александр Александрович, ген. директор, ecology@tsu. tula.ru, Россия, Пермь, УК Западно-Уральского машиностроительного концерна»),

Стась Галина Викторовна канд. техн. наук, доц., galina_stas@,mail.ru, Россия, Тула, Тульский государственный университет

AEROGASDYNAMICS PROCESSES A T DIFFERENT MINES PRODUCTION FACES

N.M. Kachurin, I.I. Mochnachuk, A.A. Pozdeev, G.V.Stas

The system approach to solving dynamics tasks of gas emission at different mines production faces is proposed. It's shown that structural elements of coal mine and mine by mining different metals or constructional materials and mine by mining potassium ore production faces always consists of identical mining and geological and geotechnological objects. Aerogasdynamics connection exists between the structural elements of any production face. Mathematical models of aerogasdynamics processes in production face workings make possible to get reliable evaluation of gas situation and to be theoretical base for dynamical method of calculating quantity of air.

Key words: aerogasdynamics, gas emission, filtration, diffusion, gas admixture, air, ventilation jets, production face, mathematical model.

Kachurin Nikolay Mikhailovich, doctor of technical sciences, professor, the head of chair, ecology@tsu.tula.ru, Russia, Tula, Tula State University,

Mochnachuk Ivan Ivanovich, candidate of econ. sciences, , chairman of Russian Independent Union of mining workers, ecology@tsu.tula.ru, Russia, Moscow,

Pozdeev Alexander A., General Director, ecology@tsu.tula.ru, Russia, Perm, Management company of West-Ural machine-building concern,

Stas Galina Viktorovna, candidate of technical sciences, docent,galina_stas@mail.ru, Russia, Tula, Tula State University