CC BY
10
I. АКТУАЛЬНО I. URGENT
| С.Б. Романченко // S. B. Romanchenko romanchenkosb@mail.ru
д-р техн. наук, доцент, ведущий научный сотрудник ФГБУ ВНИИПО МЧС России, Россия, 143903, Московская область, г. Балашиха, мкр. ВНИИПО, д. 12 doctor of technical sciences, assistant professor, leading researcher of FGBU VNIIPO MChS of Russia, microdistrict 12, VNIIPO, Balashikha, Moscow Region, 143903, Russia
УДК 622.81
АСПЕКТЫ ИССЛЕДОВАНИЯ ДИНАМИКИ РУДНИЧНЫХ АЭРОЗОЛЕЙ В ПРОЦЕССЕ ИХ ПЕРЕМЕЩЕНИЯ ПО СЕТИ ГОРНЫХ ВЫРАБОТОК
MINE AEROSOL DYNAMICS RESEARCH ASPECTS IN THE PROCESS OF THEIR MOVEMENT ALONG THE MINE OPENING NETWORK
Перемещение аэрозолей в пределах одной горной выработки сопровождается динамическими процессами агрегации и седиментации частиц. В различных сечениях горных выработок аэрозоли отличаются по вещественному составу, массе витающих частиц (выражаемых объемной концентрацией С, мг/м3) и имеют существенные изменения дисперсного состава аэрозоля. В статье рассмотрены методы исследований процессов пылевой динамики, определение численных показателей и функций дисперсного состава на современной лабораторной базе - электронная микроскопия, лазерный анализ. Показана аппаратная зависимость результатов определения функций плотности распределения частиц от числа контролируемых диапазонов эквивалентных диаметров или каналов лазерного анализатора. Предложены аппаратно-независимые способы аппроксимации функций. Рассмотрено применение фракционно-массовых функций распределения при исследовании динамики рудничных аэрозолей.
The movement of aerosols within the same mine is accompanied by dynamic processes of particles' aggregation and sedimentation. In different sections of mine workings, aerosols differ in material composition, the mass of floating particles (expressed by the volume concentration C, mg /m3) and have significant changes in aerosol dispersed composition. The article reviews the dust dynamics process research methods, determining the dispersed composition numerical indicators and functions on a modern laboratory base - electron microscopy, laser analysis.
The particle distribution density functions' determining results hardware dependence on the number of the laser analyzer equivalent diameters' controlled ranges or channels is shown. Hardware-independent methods of function approximation are proposed. The application of fractional-mass distribution functions in the study of mine aerosol dynamics is considered.
Ключевые слова: УГОЛЬНЫЙ АЭРОЗОЛЬ, ТВЕРДАЯ ДИСПЕРСНАЯ ФАЗА, ДИСПЕРСНЫЙ СОСТАВ, ЭКВИВАЛЕНТНЫЙ ДИАМЕТР, МОДА ДИАМЕТРА, СРЕДНЕЕ КВАДРАТИЧЕСКОЕ ОТКЛОНЕНИЕ, РАСПРЕДЕЛЕНИЕ ГАУССА, МИКРОСКОПИЧЕСКИЕ ИССЛЕДОВАНИЯ,
СКАНИРУЮЩАЯ МИКРОСКОПИЯ, ЛАЗЕРНЫЙ АНАЛИЗ, СТРУКТУРА ЧАСТИЦ, ПЛОТНОСТЬ ПЫЛИ
Key words: COAL AEROSOL, SOLID DISPERSED PHASE, DISPERSED COMPOSITION, EQUIVALENT DIAMETER, DIAMETER MODE, MEAN QUADRATIC DEVIATION, GAUSSIAN DISTRIBUTION, MICROSCOPIC STUDIES, SCANNING MICROSCOPY, LASER ANALYSIS, PARTICLE STRUCTURE, DUST DENSITY
Введение
К аэрозолям относятся многофазные системы, в которых мелкодиспергиро-ванные жидкие частицы и/или твердые частицы (соответственно жидкая и твердая дисперсные фазы) поддерживаются определенное время во взвешенном состоянии внутри газо-но-сителя. Рудничный воздух (газ-носитель) именуется также дисперсионной средой. Применительно к угольным шахтам в месте разрушения горного массива обязательным является применение систем орошения, результатом работы которых является создание водных аэрозолей с размером капель 40-60 мкм (мелкодисперсные системы) или ~15 мкм в системах эжекторного орошения. Водные аэрозоли в месте разрушения массива смешиваются с твердыми угольными, породными или угле-породными частицами, и в шахтах происходит образование разнородных 3-х фазных аэрозолей с одной жидкой фазой (вода) и не менее 2-х твердых фаз: угольные и породные частицы.
Твердая дисперсная фаза рудничных аэрозолей содержит различимые в микроскоп угольные и породные частицы (рисунок 1).
Твердые частицы в основе имеют неправильную изометрическую форму, доля продолговатых частиц не превышает 15%. Из-за количественно-качественных изменений на микроуровне происходят изменения в структуре и угольные пылевые частицы имеют существенно
большую пористость чем частицы угля фракциями более 1 мм (рисунок 2). Это находит отражение в величине фактической плотности частиц (отношение массы к видимому объему, на который воздействуют аэродинамические силы).
Перемещение аэрозолей в пределах одной горной выработки сопровождается динамическими процессами агрегации и седиментации частиц. В различных сечениях горных выработок (например, начальном и конечном) аэрозоли отличаются по вещественному составу, массе витающих частиц, выражаемых объемной концентрацией С, мг/м3, и имеют существенные изменения дисперсного состава. Методы исследований указанных процессов, численных показателей и функций, ориентированных на процессы пылевой динамики, рассмотрены в рамках настоящей статьи.
1. ПРОБЛЕМНЫЕ ВОПРОСЫ ПРИ ОПРЕДЕЛЕНИИ ФУНКЦИЙ РАСПРЕДЕЛЕНИЯ ПЫЛЕВЫХ ЧАСТИЦ
Дисперсный состав (или дисперсность) является существенным фактором, влияющим на показатели пожаровзрывобезопасности горючих пылей - максимальное давление взрыва - Pmx и нижний концентрационный предел распространения пламени (НКПР или НПР). На основе работ МакНИИ [5], KD Barbara (Польша) [3,4], работ ВНИИПО [8,9,13 и других исследователей [14,15] экспериментально установлено,
Рисунок 1. Угольные и породные частицы рудничного аэрозоля («П» - порода) Figure 1. Mine aerosol coal and rock particles ("P" - rock)
*> '
1 яь i ■ Л
Л v-ус ш - ■*
RaI ^ i' г
Щ ■ .А
' -т-v 5 !
Л
г 1
* л Л
У i ?
SEC 2GkV VUniOmm
Рисунок 2. Угольные частицы с эквивалентным диаметром 17-55 мкм Figure 2. Coal particles with an equivalent diameter of 17-55 microns
что величина Ртах существенно возрастает при уменьшении среднего диаметра частиц (ё). При уменьшении ё с 850 мкм до 150 мкм рост Ртх незначителен, однако при дальнейшем снижении размеров частиц происходит резкий рост давления взрыва: для углей с У" >40% происходит почти десятикратный рост давления взрыва при уменьшении размеров эквивалентных диаметров частиц с 150 до 10 мкм. Это особо подчеркивает актуальность постоянного контроля дисперсного состава, который выражается весовыми или объемными функциями распределения ^(й) и F(d)), а также весовыми или объемными функциями плотности распределения (&(й) иЯй)) [1,2].
В многочисленных исследованиях последнего периода [7, 16] микроскопическими исследованиями выявлена повышенная трещинова-тость частиц малых эквивалентных диаметров, численно выражающаяся в росте удельной поверхности при уменьшении диаметра частиц. Одновременно с этим 3^ сканирующие микроскопы позволяют явно выделить образующуюся на поверхности угольных частиц поверхностную пленку толщиной 120-150 нм [16]. Пленка содержит окислы сложного угольного вещества, физически отделена от основной части частиц, окислы имеют плотность ниже плотности не-окисленного вещества. Процентная доля пленок в видимом объеме частиц резко возрастает при увеличении дисперсности: для частиц с й>100 мкм она составляет 0,4 %, для 10 микронных частиц она возрастает до 4,4%, а в частицах в один микрон доля объема окислительной пленки в объеме частиц составляет 39%. Суммарно указанные факторы делают различными плотность частиц разных диаметров, даже для однородного аэрозоля. Физически возможно измерение
только насыпной плотности частиц, которая для мелких фракций угольной пыли (0<20 мкм) снижается до 0,25-0,30 г/см3 при плотности угля ~1,3 г/см3. При пересчете от насыпной к фактической плотности из мерного объёма (в рассмотренном случае 1 см3) вычитается объем пустот между частицами. Как показывают элементарные математические преобразования, объем пустот между частицами практически не зависит от степени дробления сыпучего вещества, при этом фактическая плотность для мелких частиц остаётся в 2-3 раза менее плотности угля и примерно в 2 раза менее плотности воды.
Существенный разброс плотности частиц как для однородных, так и для разнородных рудничных аэрозолей делает проблемным переход от объемных функций распределения F(d) и функций плотности распределения Дй) частиц (по классификации Н.А.Фукса [1,2 ) к весовым ^ и g соответсвенно). В работах Н.А.Фукса плотность частиц для однородного аэрозоля у принимается равной константе. В этом случае весовое распределение пропорционально объемному. При разнородном составе аэрозоля, при разнородных или агрегированных аэрозолях (уФсот() микроскопические или лазерно-оптические измерения позволяют получать только функции распределения по объему и их сложно адаптировать к нормативным требованиям, ориентированным на массовые показатели (мг/м3 и т.д.).
Необходимо отметить, что для функции Дй) в классической теории вероятности выделено её основное свойство: функция плотности распределения Дй) не универсальна, она существует только для непрерывных случайных величин. Число витающих частиц N в рудничных аэрозолях крайне велико (#>>109), в пределах нескольких метров горной выработки речь может
8
Рисунок 3. Функции f(d) и F(d) при их численном определении Figure 3. Functions f (d) and F (d) in their numerical determination
идти о числе частиц порядка N ~'0'5. Несмотря на большое число частиц это число N конечно, и построение функции плотности распределения возможно только её аналитическим построением или аппроксимацией. На практике происходит попытка экспериментального определения Дй) по ограниченному числу экспериментальных точек, что делает погрешность такого эксперимента неопределенной.
Вторым важным свойством функции плотности распределения Дй) является значение величины интеграла Дй) в пределах от 0 (формально от -да) до максимально возможной величины частиц йтш=84'^850 мкм1, относимых к пыли (формально до +да):
Г= 1 }0 (1) При вычислении Дй) не в долях, а в процентах значение интеграла (1) равно 100%.
При экспериментальном определении Дй) эксперимент сводится к установлению интервалов частиц от сС. до сС. + Дй, в пределах которых функция предполагается равной константе и установлении величины/(й). Число указанных интервалов для наиболее распространенных в настоящее время моделей лазерных анализаторов составляет 62^66 ( принимает значения от 1 до 62^66) и экспериментально определенное значениеДй) имеет вид гистограммы (рисунок 3).
Увеличение числа экспериментальных точек на уровне теоретических рассуждений должно приводить к уточнению получаемых значений Дй). Именно это предположение описано в фундаментальных работах Н.А.Фукса. Рисунок 4 (3 левых фрагмента) является воспроизведением рисунка №5 из работ [1,2]. В этом случае не
1 В РФ это частицы до 850 мкм, в США, Австралии и др.
угледобывающих странах к взрывоопасным относятся частицы, проходящие через сито № 20 (841 мкм). К витающей пыли, аналогично нормам РФ, в большинстве стран относятся частицы, проходящие через сито № 200 (74 мкм).
учтены особенности функций с фиксированной суммой (100%). Сумма 3-х слагаемых значений Дй) для 3-х интервалов (нижний левый фрагмент рисунка 4) должны дать 100%, так же как и сумма 8 и 16 слагаемых на вышерасположенных фрагментах. Очевидно, что это уменьшает абсолютное значение каждого слагаемого и графически гистограмма должна перемещаться не снизу вверх, а сверху вниз при увеличении числа интервалов диаметров частиц. Результаты перехода от 62 интервалов к 124 интервалам лазерного анализатора наглядно демонстрируют указанный эффект на правом фрагменте рисунка 4. Вместо перемещения вверх (левые фрагменты рисунка 4) рост числа экспериментальных интервалов с 62 до 124 приводит к резкому снижению значений Дй) - максимальное значение функции плотности распределения достигавшееся на величине диаметра частиц в 75 мкм при 62 диаметрах /62(75)=8,87%0 снижается практически в 2 раза при увеличении числа интервалов до 124 /■ Ш(75М,43%0).
Из проведенного анализа видно, что получение функции плотности распределения частиц является неустойчивым и аппаратно-зави-симым процессом. Однако независимо от числа контролируемых интервалов частиц постоянными величинами при проведенных экспериментах является модальное значение М=75 мкм (соответствует максимуму функции Дй)). Также постоянным является интервал размеров частиц от сС0=0,2 мкм до йтш='45 мкм, экспериментально определенный в рассмотренной пробе пыли. Это позволяет сделать вывод о возможности получения аналитического описания (математической зависимости) по достоверным аппаратно-независимым параметрам.
В качестве одной из экспериментально подтвержденных версий рассмотрено получе-
ние нормально распределенной функции на основе аппаратно-независимых данных от лазерных анализаторов дисперсного состава пыли. С учетом центральной предельной теоремы из классической теории вероятности - случайная величина (в нашем случае - сС - диаметр пылевых частиц), образующаяся в результате суммирования независимых либо слабо зависимых случайных процессов, распределена по нормальному закону Гаусса с параметрами Мс и а:
(2)
Одним из неоспоримых преимуществ распределения (2) является безусловное выполнение условия (1). Доказательство проводится методом замены переменных и приведению выражений к функции Лапласа. Одновременно необходимо отметить, что другие отмеченные в литературе законы распределения пылевых частиц: «нормально-логарифмический», Роллера, Разина-Рамблера и др. не приводят математических обоснований и граничных условий выполнения (1). Это особенно актуально с учетом сложности интегрирования функций, содержащих логарифмы с основанием 10 для так называемого «нормально-логарифмического» закона распределения. Формально это не вполне корректное название - оно не содержит нормальных логарифмов, а представляет собой квазинормальную функцию асимметричного «колоколоо-
бразного» вида, содержащую вместо диаметра й десятичный логарифм диаметра
В значительном числе экспериментов 6,7] показана близость возможного представления плотности распределения витающих аэрозольных частиц (с размерами до 74-100 мкм) как в виде нормального, так и нормально-логарифмического законов.
На основе этого можно предположить, что а - среднеквадратическое отклонение случайной величины с1 может быть определено:
= ¿тах-бд (3)
6
Таким образом, аппроксимация функций плотности распределения частиц (в аппаратно-независимых условиях) может достигаться путем экспериментального получения моды Мс и а - среднеквадратического отклонения диаметров частиц с последующим представлением функции плотности распределения частиц в виде уравнения Гаусса.
2. ФУНКЦИИ ДЛЯ ИССЛЕДОВАНИЯ ДИНАМИЧЕСКИХ СОСТОЯНИЙ АЭРОЗОЛЯ
Стендовые и шахтные исследования динамики аэрозолей показывают, что в процессе перемещения аэрозоля как в пределах одной выработки, так и разветвленной сети выработок, происходит одновременное изменение массы витающих частиц и изменяется дисперсный состав. Для отображения указанных изменений в
Рисунок 4. Фактический (справа) и «традиционный» (слева) [1] вид дифференциальной функции распределения
(объемная функция по классификации Н.А.Фукса) для различного числа контролируемых интервалов частиц Figure 4. Actual (right) and "traditional" (left) [1] view of the differential distribution function (volume function according to the classification of N.A. Fuchs) for a different number of controlled particle intervals
10
Рисунок 5. Дисперсный состав пыли в различных сечениях выработки Figure 5. Dispersed composition of dust in various mine sections
составе витающей пыли применим прием, при котором функция плотности распределения частиц в каждом из исследуемых сечений выработки умножается на концентрацию пыли в данном сечении. Полученные функции именуются как фракционно-массовые [7].
На рисунке 5 представлены функции плотности распределения пыли, определенные соответственно в сечениях на расстоянии 3 м, 6 м, 9 м и 12 м от начала выработки (по ходу вентиляционной струи).
Поскольку функцияf(d) на рисунке 5 выражается в процентах, то независимо от массы витающих частиц сумма всех значений f(d) во всех сечениях одинакова или 100%.
Даже при наличии изменений в дисперсном составе аэрозоля функции f(d) для различных сечений налагаются и трудно различимы (рисунок 5). Переход к фракционно-массовым функциям (рисунок 6) позволяет упростить процесс их исследования, наиболее понятной является их визуальное отличие: одновременно видна масса витающих частиц (площадь фигуры) и дисперсный состав пыли в каждом сечении горной выработки.
Преимущества фракционно-массовых функций для исследования динамики аэрозолей очевидны из сопоставления рисунков 5 и 6. Дополнительно к указанным свойствам фракци-онно-массовые функции позволяют выполнять дополнительный контроль достоверности измерений. Преимущества фракционно-массовых функций изложены также в работе [6,7]. При отсутствии внутренних источников пылевыделе-ния достигается условие непревышения каждой из фракций пыли на последующих этапах или визуально (рисунок 6) «эффект вложенности» кривых.
ВЫВОДЫ
Рисунок 6. Динамика массовой концентрации и дисперсного состава пыли Figure 6. Dynamics of dust mass concentration and dispersed composition
1. Определение функции плотности распределения пылевых частиц f(d) при помощи лазерного анализа (независимо от моделей и фирм производителей) определяется величиной i - количеством контролируемых интервалов размеров частиц (ё.+Дё), а также шириной этих интервалов (Дё). Изменение числа интервалов приводит к изменению значений f(d), и экспериментальное получение f(d) является аппа-ратно-зависимым процессом. Существенным недостатком современного лазерного анализа является отсутствие встроенного контроля по уравнению (1) для экспериментально определенной f(d), что делает неопределенной полученную погрешность.
2. Аппаратно-независимым для мономодальных аэрозолей является получение моды размеров диаметра частиц Md и а - среднеква-дратического отклонения.
3. Исходя из центральной предельной теоремы теории вероятности размеры угольных частиц после их образования могут быть распределены по нормальному закону Гаусса, использующего два указанных выше параметра Мё и а. Практическая возможность описания плотности распределения частиц в рудничных аэрозолях законом Гаусса подтверждена экспериментально для диапазона размеров частиц 1-74 мкм [6,7]. Это открывает существенные возможности контроля и уточнения результатов экспериментального определения f(d) на применяемых физических принципах (ситовый или лазерный анализ, седиментационные весы и т.д.).
4. Для исследований динамики аэрозолей по мере их перемещения по сети горных выработок предложена и практически апробирована фракционно-массовая функция, одновременно фиксирующая изменения в концентрации аэрозоля и в его дисперсном составе.
СПИСОК ЛИТЕРАТУРЫ
1. Фукс Н. А., Механика аэрозолей.- М.: Издательство Академии наук СССР, 1955.- 352 с.
2. Fuchs N.A. The mechanics of aerosols. Oxford: Pergamon Press, 1964. -360 p.
3. Cybulski W. Wybuchy pylu weglowego i ich zwalczanie. Katowice, Wydaw. "Slask"., 1973.- 451 с.
4. Lebecki.K.Zagrozenia pyiowe w gornictwie. Katowice: Giowny Instytut Gornictwa, 2004, -399 s.
5. Нецепляев М.И., Любимова А.И., Петрухин П.М.. Борьба со взрывами угольной пыли в шахтах. М.: Недра, 1992. - 298 с.
6. Лебецки К.А., Романченко С.Б. Пылевая взрывоопасность горного производства. -М.: Горное дело, 2012.-464 с.
7. Романченко С.Б., Руденко Ю.Ф., Костеренко В.Н. Пылевая динамика в угольных шахтах.-М.: Горное дело, 2011.-256 с.
8. Романченко С.Б., Девликанов М.О. Влияние дисперсного состава угольной пыли на показатели взрывоопас-ности.// Вестник научного центра по безопасности работ в угольной промышленности. - № 2-2019, с.16-23.
9. Романченко С.Б., Костеренко В.Н. Полномасштабные исследования взрывов угольной пыли и критерии эффективности средств локализации.// Вестник научного центра по безопасности работ в угольной промышленности. - 2018. №4. С.6-20.
10. Грин Х., Лейн В. Аэрозоли-пыли, дымы и туманы / пер. с англ. Л.: Химия, 1969. - 428 с.
11. Райст П. Аэрозоли. Введение в теорию. М.: Мир, 1987. 280 с.
12. Фрич. Измерение размера частиц.- Обзор продукции фирмы Фрич. Идар-Оберсен, 2010, - 28 c.
13. Шебеко Ю.Н., Навценя В.Ю., Копылов С.Н., Замышевский Э.Д., Шебеко Д.Ю. Экспериментальные исследования искробезопасности материалов в различных взрывоопасных средах.// Пожарная безопасность.-2000, №4.-С.122-126.
14. Шоль Э.В. Возникновение взрыва метана и угольной пыли и их предотвращение. Глюкауф, 21/22, 1989.-С.9-11.
15. Кэшдоллар К., Херцберг М. Камера объемом 20 л для исследования взрываемости порошков и газов.Ежеме-сячный журнал американского института физики. т.56, 1985, №4, с.108-115.
16. Романченко С.Б., Трубицин А.А., Кубрин С.С. Проблемы определения фактической плотности угольных частиц в процессах витания и седиментации // Вестник научного центра по безопасности работ в угольной промышленности. - 2020. №1. С.6-14.
REFERENCES
1. Fuchs, N.A. (1955). Mekhanika aerozoley [The mechanics of aerosols]. Moscow: Izdatel'stvo Akademii nauk SSSR [in Russian].
2. Fuchs, N.A. (1964). The mechanics of aerosols. Oxford: Pergamon Press [in English].
3. Cybulski, W. (1973). Wybuchy pylu weglowego i ich zwalczanie. Katowice, Wydaw. "Slask". [in Polish].
4. Lebecki, K. (2004). Zagrozenia pyiowe w gornictwie. ^towice: Giowny Instytut Gornictwa [in Polish].
5. Netsepliaev, M.I., Liubimova, A.I., & Petrukhin, P.M. (1992). Bor'ba so vzryvami ugol'noy pyli v shakhtakh [Coal dust explosion suppression in mines]. Moscow: Nedra [in Russian].
6. Lebecki, K.A., & Romanchenko, S.B. (2012). Pylevaya vzryvoopasnost' gornogo proizvodstva [Mining dust explosion hazard]. Moscow: Gornoie delo [in Russian].
7. Romanchenko, S.B., Rudenko, Yu.F., & Kosterenko, V.N. (2011). Pylevaya dinamika v ugol'nykh shakhtakh [Dust dynamics in coal mines]. Moscow: Gornoie delo [in Russian].
8. Romanchenko, S.B., & Devlikanov, M.O. (2019). Vliyaniye dispersnogo sostava ugol'noy pyli na pokazateli vzryvoopasnosti [Coal dust dispersed composition influence on explosion hazard indicators]. Vestnik nauchnogo tsentra po bezopasnosti rabot v ugolnoi promyshlennosti - Herald of Safety in Mining Industry Scientific Center,2, 16-23 [in Russian].
9. Romanchenko, S.B., & Kosterenko, V.N. (2018). Polnomasshtabnyye issledovaniya vzryvov ugol'noy pyli i kriterii effektivnosti sredstv lokalizatsii [Full-scale studies of coal dust explosions and criteria for the locating means effectiveness]. Vestnik nauchnogo tsentra po bezopasnosti rabot v ugolnoi promyshlennosti - Herald of Safety in Mining Industry Scientific Center, 4, 6-20 [in Russian].
10. Green, H., & Lein, V, (1969). Aerozoli - pyli, dymy i tumany [Aerosols - dusts, fumes and mists]. Leningrad: Khimia [in Russian].
11. Raist, P. (1987). Aerozoli. Vvedeniye v teoriyu [Aerosols. Introduction to theory]. Moscow: Mir [in Russian].
12. Frich. (2010). Izmereniye razmera chastits.- Obzor produktsii firmy Frich. Idar-Obersen [Fritsch. Particle size measurement. - Fritsch product overview. Idar-Obersen [in Russian].
13. Shebeko, Yu.N., Navtsenia, V.Yu., Kopylov, S.N., Zamyshevsky, E.D., & Shebeko, D.Yu. (2000). Eksperimental'nyye issledovaniya iskrobezopasnosti materialov v razlichnykh vzryvoopasnykh sredakh [Experimental studies of materials' spark-proofness in various explosive environments]. Pozharnaia bezopasnost - Fire safety, 4, 122-126 [in Russian].
14. Shol', E.V. (1989). Vozniknoveniye vzryva metana i ugol'noy pyli i ikh predotvrashcheniye [Methane and coal dust explosion initiation and their prevention]. Gluckauf, 21/22, 9-11 [in Russian].
15. Keshdollar, K., & Khertsberg, M. (1985). Kamera ob"yemom 20 l dlya issledovaniya vzryvayemosti poroshkov i gazov [20 liters volume chamber for investigating the explosiveness of powders and gases]. Yezhemesyachnyy zhurnal amerikanskogo instituta fiziki - American Institute of Physics Monthly Journal, v. 56, 4, 108-115 [in Russian].
16. Romanchenko, S.B., Trubitsyn, A.A., & Kubrin, S.S. (2020). Problemy opredeleniya fakticheskoy plotnosti ugol'nykh chastits v protsessakh vitaniya i sedimentatsii [Problems of determining the actual density of coal particles in the processes of floating and sedimentation], Vestnik nauchnogo tsentra po bezopasnosti rabot v ugolnoi promyshlennosti - Herald of Safety in Mining Industry Scientific Center, 1,6-14 [in Russian],
