CC BY
22
УДК 622.831:622.016.6
А.Ю. Ермаков
ВЫДЕЛЕНИЯ МЕТАНА С ПОВЕРХНОСТИ
ОБНАЖЕНИЯ УГОЛЬНОГО ПЛАСТА ПРИ ВЫСОКОЙ СКОРОСТИ ПОДВИГАНИЯ ОЧИСТНЫХ И ПОДГОТОВИТЕЛЬНЫХ
ЗАБОЕВ
Теоретические исследования показали, что для прогнозирования процесса формирования метановой опасности в подготовительных выработках необходимо учитывать конечную скорость распространения давления метана в угольном пласте. Закон сопротивления при фильтрационном движении метана в угольном пласте представляет собой функциональную связь газового потока с градиентом давления газа и локальной скоростью изменения газового потока. Процессы фильтрации метана в угольном пласте описываются уравнением гиперболического типа, а применение уравнений параболического типа является физически обоснованным для длительных периодов времени. Интенсивность метановы-деления связана с планограммой работ в очистном забое и в зависимости от скорости подвигания забоя в цикле метановыделение нарастает за счет увеличения площади газо-отдающей поверхности. Газоотдающая поверхность находится под перепадом давления, а площадь частично дегазированной поверхности обнажения уменьшается с той же скоростью, с которой увеличивается площадь не дегазированной газоотдающей поверхности, поэтому при работе выемочного комбайна формируется стационарное состояние.
Ключевые слова: выделение, метан, поверхность обнажения, угольный пласт, скорость подвигания, очистной, подготовительный, забой.
Современные технологии и технические средства отработки угольных пластов обеспечивают увеличение производительности очистных участков на порядок.
Сравнительный анализ аэрогазодинамического состояния шахт Кузбасса показывает, что, во-первых, шахты характеризуются исключительно высокой газообильностью, во-вторых, возможностью возникновения аварийных загази-рований горных выработок и появлением взрывоопасных концентраций метана. Теоретические исследования [2—12]
DOI: 10.25018/0236-1493-2018-4-0-98-105
показали, что для прогнозирования процесса формирования метановой опасности в подготовительных выработках необходимо учитывать конечную скорость распространения давления метана в угольном пласте.
При этом закон сопротивления при фильтрационном движении метана в угольном пласте представляет собой функциональную связь газового потока с градиентом давления газа и локальной скоростью изменения газового потока для произвольно выбранной точки в рассматриваемой области угольного пласта.
ISSN 0236-1493. Горный информационно-аналитический бюллетень. 2018. № 4. С. 98-105. © А.Ю. Ермаков. 2018.
А газовая проницаемость угольного пласта представляет собой обобщенную характеристику свойств метана, коллек-торских свойств угольного и релаксации процесса фильтрационного движения метана, при этом численное значение газовой проницаемости пропорционально третьей степени эффективной пористости горного массива.
Процессы фильтрации метана в угольном пласте описываются уравнением гиперболического типа, а применение уравнений параболического типа является физически обоснованным для длительных периодов времени. Современные технологии и технические средства отработки угольных пластов обеспечивают увеличение производительности
очистных участков на порядок. Разумеется, что это приводит к тому, что в несколько раз возрастает скорость под-вигания очистных и подготовительных забоев. Рассматривая угольные пласты, которые отрабатывают по технологии шахта — лава, процесс фильтрационного движения метана считать ламинарным и одномерным. Расчетные схемы выделения метана с поверхности обнажения разрабатываемого угольного пласта в подготовительные выработки и очистные забои для таких условий представлены на рис. 1.
Физически обоснованно следующее математическое описание поля давлений метана в угольном пласте, имеющем поверхность обнажения
Рис. 1. Расчетные схемы выделения метана с поверхности обнажения разрабатываемого угольного пласта: подготовительная выработка (а); очистной забой (б)
dp2 . d2 p2 _ dp^ . cfp^ _ дУ
dt + r dt2 "x dx2 dt + r dt2 "x ax2
p (z,0) = po = const, —p (z,0) = 0, p(0, t) = p1 = const, limp *a>,
при z
(2)
где р — давление свободного метана в трещиновато-пористой структуре угля; г — пространственная координата; t — время; tr — период релаксации процесса ламинарной фильтрации метана в угольном пласте; х — пьезопроводность угольного пласта, зависящая от фильтрационных и сорбционных свойств угля; р0, р1 — давление свободного метана в угольном пласте в начальный момент времени и на поверхности обнажения пласта соответственно.
Решение уравнения (1) для условий (2) имеет следующий вид
2 2 1 v^=Ц
Pl - Po I
0,5z(x)-05 ] + 0,5z(x)-05 J(e2 - z2trx-1 )-05 exp(-0,59t-1 )x
(3)
xi
0,5t-1 (e2 - z2trx-1 )0
dQ\a
t - z
( x-1 )0
где I.
o,5t-1 (e2 - z2trx-1 )0
— модифицированная функция Бесселя первого по-
рядка; ст0 [г - z(tr / х)0,5 ] — единичная функция Хэвисайда.
Из соотношения (3) следует аналитическая зависимость для расчета метановы-деления с единичной площади поверхности обнажения угольного пласта [1]
J1 = J0exp (-0,51 / 1Г) (0,5* / 1Г), (4)
где ]1 — метановыделение с единичной площади поверхности обнажения угольного пласта; J0 — начальная скорость газовыделения; /0(0,5^г) — модифицированная функция Бесселя нулевого порядка для аргумента, записанного в круглых скобках.
В период проведения выработки дебит метана в подготовительную выработку J2 с элементарной поверхности обнажения пласта с/Э можно определить как
Ш2 = птХехр (-0,5t / гг) х 10 (0,5t / гг),
где п — количество поверхностей обнажения угольного пласта, контактирующих с атмосферой подготовительной выработки; тп — мощность разрабатываемого угольного пласта; — скорость подвигания подготовительного забоя. Интегрируя уравнение,получим t
J2 () = птХJ01ехр (-0,5£ / г, ) (5)
0
где Т — срок проведения подготовительной выработки.
Следовательно, максимальное значение метановыделения Jm будет при t = Т, т.е. Jm = J0(T). В период обособленного проветривания выработки дебит метана в подготовительную выработку J3 с элементарной поверхности обнажения угольного пласта с/Э можно определить как
Ш3 = птп11У1 ]0 ехр(-0,5* / *г) х /0 (0,5* / *г),
где — проектная длина подготовительной выработки. Интегрируя уравнение, получим
J3 () = Jm ехр[-0,5((-Т)/ I,] 10 [0,5((-Т)/ I,] при t > Т. (6)
Модифицированная функция Бесселя нулевого порядка может быть представлена в интегральном виде, тогда, вводя обозначение т = 0,5^г, можно записать
J2 (т) = 0,318птп1гУ±J001 (т), при т < т1 (7)
Jm (х± ) = а^пт^А (х±); (8)
->з (т) = -402 (т-т1) при т > т1 (9)
где т1 = т1 — безразмерный срок проведения подготовительной выработки;
т п
0± (т)= ]ехр(-£) |[ехр(cos£) + ехр( cos
(10)
Рис. 2. Графики зависимостей &/т) и &2(т) для различных интервалов значений переменной т; период проведения подготовительной выработки (а); период обособленного проветривания подготовительной выработки (б)
п
02 (т) = 0,159 ехр[-(т -т1 )] Цехр[{т-т1 )cos + ехр[-(т-т± )cos С]} (11)
0
Функции 01(т) и 02(т) представляют собой безразмерные значения метано-выделения с поверхности обнажения угольного пласта в соответствующие моменты времени, которые присваиваются переменной т. Графики функций 01(т) и 02(т) представлены на рис. 2.
Результаты вычислительных экспериментов, проведенных для возможных интервалов изменения исходных данных, отражающих условия проведения подготовительных выработок в шахтах Кузбасса, позволили получить инженерные формулы для расчета этих функций. Аппроксимации функций 01(т) и 02(т) приведены в табл. 1—2.
Аппроксимации, заданные соотношениями табл. 1, характеризуются значениями корреляционного отношения от 0,995 до 0,999, а заданные соотношениями табл. 2 — корреляционными отношениями от 0,967 до 0,996. Пог-
Таблица 1
Аппроксимации зависимости 01 = 01(т)
решность аппроксимации не превышает 1%. Из формул (7)—(8) следует, что интенсивность метановыделения непосредственно связана со скоростью подвигания подготовительного забоя и проектной длиной подготовительной выработки.
Анализ зависимостей показывает, что в процессе проведения выработки мета-новыделение нарастает за счет увеличения площади обнаженной поверхности угольного пласта. Поверхность обнажения пласта вначале находится под перепадом давления приблизительно равным значению Р°; - Р2. Затем происходит дегазация краевой части угольного пласта, контактирующей с атмосферой подготовительной выработки, в соответствии с закономерностью (6). Аппроксимация кривой газового истощения приведена в табл. 2. Метановыделение с поверхности обнажения разрабатываемого угольного
Интервал безразмерного периода времени т Теоретически обоснованная аппроксимирующая формула
т е [0, 10] 01 (т) = 0,3787т3 -2,4678т2 + 6,7909т-0,367
т е [10, 100] 01 (т) = 0,0003т3 - 0,0329т2 +1,5865т +1,177
Таблица 2 Аппроксимации зависимости 02 = 02(т)
Интервал безразмерного периода времени т Теоретически обоснованная аппроксимирующая формула
т е [0, 3] е2 (т) = -0,0467т3 + 0,3314т2 - 0,8283т + 0,9957
т е [3, 10] 02 (т) = -0,0003т3 + 0,0074т2 - 0,0738т + 0,4012
т е [10, 50] е2 (т) = 5 • 10-3 т2 - 0,047т + 0,0908
т е [50, 100] е2 (т) = 2 • 10-6 т2 - 0,0008т + 0,0908
пласта в очистной забой также описывается уравнением (3) и условиями (4). Следовательно, соотношение (5) и закономерность (6) отражают динамику метанообиль-ности очистного забоя, обусловленную газовыделением с поверхности обнажения пласта. Поэтому с учетом принятых обозначений можно записать, что
l-2/(2trV2) п
Js = 0,318mntrV2Jo J exp (-x)j[exp (т cos Z) + exp (-т cos Z)]CZCT (12)
0 0
где Js — метановыделение с поверхности обнажения угольного пласта в очистном забое; V2 — скорость подачи очистного комбайна; L2 — длина очистного забоя.
Js = 0,318nmntrV2J001(Т2), (13)
где т2 = L2/2trV2); т2 — безразмерное значение длительности прохождения очистным комбайном всей лавы.
Из формул (12)—(13) следует, что интенсивность метановыделения непосредственно связана с планограммой работ в очистном забое. Анализ зависимости свидетельствует о том, что при выемочном цикле метановыделение нарастает за счет увеличения площади газоотдающей поверхности. Газоотдающая поверхность также находится под перепадом давления, приблизительно равным значению Р0 - p2, а площадь частично дегазированной поверхности обнажения уменьшается с той же скоростью, с которой увеличивается площадь не дегазированной газоотдающей поверхности, поэтому при работе выемочного комбайна формируется стационарное состояние, определяемое формулой (12).
СПИСОК ЛИТЕРАТУРЫ
1. Качурин Н. М., Борщевич A. M., Качурина О. Н., Бухтияров A. A. Безопасность геотехнологий добычи угля по газовому фактору // Безопасность жизнедеятельности. — 2010. — № 5. — С. 24—28.
2. Socolov E. M. et al. System of imitation for forecasting the 137Cs migration in the radioactive trace zone at the Chernobyl Atomic Power Station failure / International Symposium on Radiation Safety. — Moscow, 1994. — Pp. 101—103.
3. Kachurin N. M. Conceptual rules of the monitoring of the «Environment — Human Health» system in the Russian Federation / The 2-nd International Symposium «Mining and Environmental Protection». — Belgrade, 1998. — Pp. 21—26.
4. Kachurin N. M., Babovnikov A. L. Gassing during the break and transport of coal in a re-treatlongwall / Development of new technologies and equipment for mine haulage and hoisting. — Budva, 2005. — С. 245—249.
5. Siemek J., Rajtar J. Simulation of gas ouflow from porousfissured media // Arch. Mining. Sci. — 1989. — 34, № 1. — Pp. 119—128.
6. Яновская М.Ф. О скорости десорбции метана из разрушенного угля / Проблемы рудничной аэрологии: сборник статей. — М.: Госгортехиздат, 1959. — С. 32—37.
7. Васючков Ю. Ф. Диффузия метана в ископаемых углях // Химия твердого топлива. — 1976. — № 4. — С. 76—79.
8. Качурин Н. М. Фильтрация газа в угольных пластах при конечной скорости распространения давления / Подземная разработка тонких и средней мощности угольных пластов: сборник статей. — Тула: ТулПИ, 1983. — С. 56—62.
9. Качурин Н.М. Прогноз газовыделений в подготовительные выработки с использованием гиперболического уравнения фильтрации / Физико-технические проблемы управления воздухообменом в горных выработках больших объемов: сборник статей. — Л., 1983. — С. 83.
10. Качурин Н. М. Газовыделение на очистных участках шахт Подмосковного бассейна при изменении давления воздуха / Подземная разработка тонких и средней мощности угольных пластов: сборник статей. — Тула: ТулПИ, 1983. — С. 74—80.
11. Качурин Н. М., Ермаков А. Ю., Сенкус Вал. В. Аэрогазодинамика очистных и подготовительных забоев при отработке мощных пологих пластов. — Кемерово: АИ «Кузбассвузиздат», 2017. — 287 с.
12. Ярунин С.А., Бухны Д.И. Расчет газодинамического состояния призабойной части угольного пласта / Вентиляция и газодинамические явления в шахтах: сборник статей. — Новосибирск, 1981. — С. 6—12. [¡233
КОРОТКО ОБ АВТОРE
Ермаков Анатолий Юрьевич — кандидат технических наук, управляющий филиалом ООО «Сибнииуглеобогащение», г. Прокопьевск.
ISSN 0236-1493. Gornyy informatsionno-analiticheskiy byulleten'. 2018. No. 4, pp. 98-105.
A.Yu. Ermakov
METHANE RELEASE FROM EXPOSED COAL AT HIGH RATE OF ADVANCE OF DEVELOPMENT AND PRODUCTION FACES
The theoretical research shows that prediction of methane hazard in development headings should account for the final rate of methane pressure propagation in coal seams. The law of resistance during methane seepage in coal seams is the functional relation between the gas flow, the gas pressure gradient and the gas flow velocity variation. Gas permeability of coal seams is a generalized characteristic of properties of methane, reservoir characteristics of coal, and relaxation of methane flow; in this case, the numerical value of gas permeability is proportional to the third degree of effective porosity of rock mass. Methane flow processes in a coal seam are described by the hyperbolic equation, and the application of the hyperbolic equations is physically valid for the long-term periods of time. The rate of methane release is connected with the operating schedule of production faces and, depending on the face advance rate per methane release cycle, is increasing due to the growth of the gas-emitting exposed surface area. The gas-emitting area is subjected to differential pressure, and the area of partly degassed exposed surface is diminishing at a rate of enlargement in the area of non-degassed emitting surface; for this reason, steady-state conditions are formed during operation of cutter-loaders.
Key words: release, methane, exposed coal surface, coal seam, advance rate, production face, development face.
DOI: 10.25018/0236-1493-2018-4-0-98-105
AUTHOR
Ermakov A.Yu., Candidate of Technical Sciences, Manager of the Branch, LLC «Sibniiugleobogaschenie», Prokopyevsk, Russia.
REFERENCES
1. Kachurin N. M., Borshchevich A. M., Kachurina O. N., Bukhtiyarov A. A. Bezopasnost' zhizne-deyatel'nosti. 2010, no 5, pp. 24-28.
2. Socolov E. M. et al. System of imitation for forecasting the 137Cs migration in the radioactive trace zone at the Chernobyl Atomic Power Station failure. International Symposium on Radiation Safety. Moscow, 1994, pp. 101-103.
3. Kachurin N. M. Conceptual rules of the monitoring of the «Environment Human Health» system in the Russian Federation. The 2-nd International Symposium «Mining and Environmental Protection». Belgrade, 1998, pp. 21-26.
4. Kachurin N. M., Babovnikov A. L. Gassing during the break and transport of coal in a retreatlong-wall. Development of new technologies and equipment for mine haulage and hoisting. Budva, 2005, pp. 245-249.
5. Siemek J., Rajtar J. Simulation of gas ouflow from porousfissured media. Arch. Mining. Sci. 1989. 34, no 1, pp. 119-128.
6. Yanovskaya M. F. Problemy rudnichnoy aerologii: sbornik statey (Problems of mine aerology: collection of articles), Moscow, Gosgortekhizdat, 1959, pp. 32—37.
7. Vasyuchkov Yu. F. Khimiya tverdogo topliva. 1976, no 4, pp. 76—79.
8. Kachurin N. M. Podzemnaya razrabotka tonkikh i sredney moshchnosti ugol'nykh plastov: sbornik statey (Underground development of thin and average power of coal seams: collection of articles), Tula, TulPI, 1983, pp. 56—62.
9. Kachurin N. M. Fiziko-tekhnicheskie problemy upravleniya vozdukhoobmenom v gornykh vy-rabotkakh bol'shikh ob"emov: sbornik statey (Underground development of thin and average power of coal seams: collection of articles), Leningrad, 1983, pp. 83.
10. Kachurin N. M. Podzemnaya razrabotka tonkikh i sredney moshchnosti ugol'nykh plastov: sbornik statey (Physical-technical problems of control of air in mines large volumes: collection of articles), Tula, TulPI, 1983, pp. 74—80.
11. Kachurin N. M., Ermakov A. Yu., Senkus Val. V. Aerogazodinamika ochistnykh i podgotovitel'nykh zaboev pri otrabotke moshchnykh pologikh plastov (Aerogasdynamics cleaning and preparatory faces during the mining of thick flat seams), Kemerovo, AI «Kuzbassvuzizdat», 2017, 287 p.
12. Yarunin S. A., Bukhny D. I. Ventilyatsiya igazodinamicheskieyavleniya vshakhtakh: sbornik statey (Ventilation and gas — dynamic phenomena in mines: collection of articles), Novosibirsk, 1981, pp. 6—12.
FIGURES
Fig. 1. Models of methane release from exposed surface of coal seam under mining: development heading (a); production heading (b).
Fig. 2. Graphs of 01(x) and ©2(t) for different ranges of the variable t: during driving development heading (a); in the period of isolated ventilation of development heading (b).
TABLES
Table 1. Approximation of the relation 01 = 01(t). Table 2. Approximation of the relation ©1 = OJt).
A
НОВИНКИ ИЗДАТЕЛЬСТВА «ГОРНАЯ КНИГА»
Кузнецов Ю.Н., Курцев Б.В., Стадник Д.А., Стадник Н.М. Научные основы формирования геоинформационной базы прогнозирования и оценки запасов угольных месторождений: Монография
Год: 2017 Страниц: 126 ISBN: 978-5-98672-475-1 UDK: 622.1; 622.28
Осуществлен анализ результатов научных исследований, тенденций и направлений совершенствования теории и практики создания систем автоматизированного проектирования горных предприятий с учетом позитивной динамики развития трехмерного моделирования угольных месторождений. Установлено, что при наличии в настоящее время мощных и высокоразвитых систем автоматизированного проектирования на первый план выходят вопросы совершенствования методов трехмерного моделирования угольных месторождений, позволяющих эффективно прогнозировать характеристики горно-геологических условий их залегания и зонировать их в автоматизированном режиме. Обоснована необходимость разработки и сформирована структура интегрированной геоинформационной базы прогнозной оценки запасов угольных месторождений, освоение запасов которых должно быть реализовано на базе технологических систем, адаптивных к специфике их природной аккумуляции. Предложены нейросетевая модель распознавания природных геоструктур, способная «работать» в условиях неполной и нечеткой информации о месторождении, и алгоритм ее практической реализации.
Ю.Н. Кузнецов Б.В. Курцев Д. А. Стадник ИЖ Стадник
НАУЧНЫЕ ОСНОВЫ ФОРМИРОВАНИЯ ГЕОИНФОРМАЦИОННОЙ БАЗЫ
ПРОГНОЗИРОВАНИЯ И ОЦЕНКИ ЗАПАСОВ УГОЛЬНЫХ МЕСТОРОЖДЕНИЙ
