CC BY
48
Г. Я. Полевщиков
д-р техн. наук, проф., заведующий лабораторией Института угля СО РАН
A. А. Рябцев
ведущий инженер Института угля СО РАН
B. П. Титов
начальник отдела Института физики полупроводников СО РАН
УДК 622.121.54
ОПРЕДЕЛЕНИЕ ГАЗОКИНЕТИЧЕСКИХ ХАРАКТЕРИСТИК УГОЛЬНЫХ ПЛАСТОВ
Разработаны научно-технические основы развития метода оценки уровней газодинамической активности угольных пластов от квазистатической метанообильности выработок до собственно внезапных выбросов угля и газа как ряда генетически общих газопроявлений. Приведены результаты натурных наблюдений и их сопоставление на основе показателей потенциальной и кинетической энергии газовой компоненты углеметановых пластов.
Работа выполнена при финансовой поддержке междисциплинарного интеграционного проекта СО РАН № 99 и партнерского интеграционного проекта СО РАН № 100.
Ключевые слова: УГОЛЬНЫЙ ПЛАСТ, ГАЗОНОСНОСТЬ, ВЫДЕЛЕНИЕ МЕТАНА,
УСТОЙЧИВОСТЬ, РАЗРУШЕНИЕ, ТРАНСПОРТИРОВАНИЕ, ВНЕЗАПНЫЙ ВЫБРОС, ЭНЕРГЕТИЧЕСКИЕ ПОКАЗАТЕЛИ
Обеспечение инновационной значимости научных исследований газодинамики угольных пластов требует направленности на создание адаптивных методов регионального и локального прогнозов, текущего контроля, предотвращения и локализации всех видов газопроявлений. Сложность достижения поставленной цели существенно снижается, если представить механизм в виде многостадийного процесса, когда в относительно «слабом» газопроявлении можно обнаружить признаки несостоявшегося в данной горно-технологической ситуации более «сильного» [1]. В результате имеем псевдопоследовательный ряд:
- квазистатическое газовыделение из приконтурной части пласта, приобретающее в условиях интенсивной зональной дезинтеграции приконтурной зоны массива динамический характер [2];
- выдавливания и обрушения угля с повышенным газовыделением (масса разрушенного угля до 20 т), возникающие при нарушении (с участием давления газа) силового равновесия элементов неупругой зоны пласта [1, 3];
- средние динамические газопроявления (масса разрушенного угля 20-100 т), возникающие при физическом или технологическом
научно-технический журнал № 2-2013
78 вестник
вскрытии части пласта, в которой энергетический потенциал газовой компоненты достаточен для ее саморазрушения («волна дробления» [4]), но недостаточен для формирования углегазового потока («волна выброса» [5]) с производительностью, соответствующей скорости фронта дробления;
- собственно внезапные выбросы угля и газа (масса разрушенного угля более 100 т), развивающиеся в зонах пластов с кинетической энергией распада углеметана, достаточной для массопереноса в выработку с производительностью, соответствующей скорости фронта дробления с учетом разрыхления угля и выпадения его основной части на почву полости выброса и выработки. Выброс затухает с ростом необходимого для транспортирования на возрастающее расстояние давления газа в полости и снижения по этой причине диспергирования угля на фронте дробления [1, 4-6].
Представленная последовательность возникновения, развития и затухания газопроявлений достаточно близка известным их классификациям [7] и методам оценки опасности по энергетическим [6, 8-10] и механическим [11] показателям. Но она позволяет рассматривать различные виды газопроявлений как следствия
единого физического процесса, различные стадии развития и затухания которого требуют выполнения соответствующих необходимых и достаточных условий. Акцент на газокинетические характеристики газоносных пластов позволил в предшествующий период описать механизм внезапных выбросов вплоть до создания инженерного метода расчета их параметров не только для угля и метана, но и угля и углекислого газа, соли и углекислого газа, песчаника и метана [1, 12]. А для внезапных выбросов угля и метана разработать методику расчета зон загазования вентиляционной сети и их локализации [13, 14]. Отметим, что без рассмотрения геомеханиче-ского фактора на стадиях средних и собственно внезапных выбросов этот акцент приводит к модели саморегулирующегося процесса, так как сколь угодно малому «живому» сечению углегазового потока может соответствовать столь же малая площадь фронта дробления. Примером этому служат факты внезапных выбросов через шпуры диаметром 42 мм. Но в больших масштабах процесс изменения напряженно-деформированного состояния пласта перед фронтом дробления, возможно, и определяет регистрируемую цикличность собственно внезапных выбросов. Например, установленная физическим моделированием дальность распространения фронта чисто газового дробления угля с коэффициентом крепости, например, 0,3 при давлении газа 5 МПа может составлять всего лишь 23 см [15]. При этом образуются частицы средневзвешенного диаметра 0,16 см, что достаточно для их транспортирования в потоке расширяющегося газа, но при условии, что доля этих частиц в общей массе разрушаемого угля, учи-
тывая коэффициент его разрыхления, должна быть не менее 30% для предотвращения само-подбучивания фронта дробления. Тогда можно принять, что установленные натурными наблюдениями циклические задержки в перемещении фронта дробления связаны с периодом времени на изменение напряжений в приконтурной части пласта, создание макротрещин и насыщение их газом, что близко стадиям высыпания и обрушения угля. Выпавший из потока уголь не является инертной массой, а интенсивно выделяя газ, поддерживает режим газотранспорта раздробленной части угля в выработку.
Развивая эту методологию, на базе новых представлений о существовании метана в угольных пластах [16, 17] был обоснован показатель потенциальной энергии полураспада углемета-на [18, 19]
Е = 0,25 (Х-Х0)/В, кДж/кг , (1)
где Х и Х0 - геологическая газоносность пласта и его остаточная газоносность после разрушения, м3/т; В - коэффициент роста газоносности с увеличением глубины залегания пласта, 1/МПа.
С его применением проанализированы статистические данные шахт Карагандинского, Кузнецкого и Печорского бассейнов о зарегистрированных в период с 1943 г. по настоящее время 197 случаях внезапных выбросов [20, 21]. Из этой выборки исключены выбросы силой менее 20 т. Оставшиеся 114 случаев распределены по группам с градацией по показателю энергии полураспада углеметана через 10 кДж/кг. В каждой группе определена средняя сила выброса (рисунок 1).
научно-технический журнал № 2-2013
вестник 79
Как видим, распределение точек в этих координатах близко к изложенной выше последовательности развития динамических газопроявлений. При Е < 85 кДж/кг явлений не зарегистрировано. Но следует также учесть, что при массе поступившего в выработку угля менее 20 т и отсутствии травм у горнорабочих карточка регистрации газопроявления, как правило, не заполняется. По этой причине границы интервала квазистатической метанообильности выработок с возможными динамическими притоками метана в условиях интенсивного развития зональной дезинтеграции приконтурной части массива следует выделить значениями показателя 0 < Ек < 60 кДж/кг. Тогда выдавливаниям и обрушениям угля с повышенным газовыделением соответствует интервал 60 < Ев < 90 кДж/кг. График силы явлений на интервале 90 < Ес < 150 кДж/кг (20-100 т) имеет максимум, что согласуется с третьей стадией развития - слабые внезапные выбросы. Здесь энергетический потенциал газовой компоненты углеметана уже достаточен для интенсивного саморазрушения геоматериала, но недостаточен для формирования углегазового потока соответствующей производительности, и возникает самоподбучивание фронта разрушения пласта в связи с разрыхлением угля. Но при Есв > 150 кДж/кг вступает в действие режим транспортирования угля в потоке расширяющегося газа и наблюдаются наиболее грозные природно-технологические явления - собственно внезапные выбросы угля и газа. Нижнее значение их силы в 100 т принято с учетом результатов исследований условий загазирования вентиляционных сетей шахт [13, 14]. Установлено, что при выбросах силой более 150 т уже возможны опрокидывания вентиляционных потоков за пределами аварийного выемочного участка с дальнейшим движением газового потока уже с участием общешахтной депрессии, т. е. масштаб аварии становится непредсказуемым.
Обоснованность предложенной оценки уровней газодинамической активности угольных пластов подтверждается анализом результатов ранее выполненных горно-экспериментальных работ по измерению удельной скорости газо-выделения при разрушении угля бурением [1]. Исследования выполнены на семи наиболее выбросоопасных в те годы шахтопластах Кузбасса. Проведено около 150 измерений, часть их них непосредственно в зонах выдавливаний от 5 до 20 т угля со значительным выделением метана. Два замера выполнены за 4 и 9 ч до аварии.
научно-технический журнал № 2-2013
Современная обработка полученных данных с применением показателя энергии полураспада углеметана позволила получить информацию об особенностях газокинетических свойств пластов на интервале 0 < Е < 90 кДж/кг. Ситовой анализ выбуриваемого угля показал, что эта газоотдача соответствует средневзвешенному диаметру частиц 0,8...2,5 мм. В качестве газодинамической характеристики пласта в точке измерений было принято отношение объема газа, выделяющегося за первые 30 с от момента мгновенного отторжения единицы массы угля, к ее начальной газоносности Пг (рисунок 2).
Рисунок 2 - Связь показателя начальной скорости газоотдачи разрушаемого угля и показателя энергии полураспада углеметана
Согласно результатам математического [4, 5] и физического [12] моделирования массопе-реноса угля в потоке расширяющегося газа, для реализации внезапного выброса угля и метана достаточно участия в нем не менее 30 % содержащегося в угле газа. Как видим из графика, это условие выполняется уже при Е = 40 кДж/кг. При этом явная нелинейность связи указывает, что реализация газосодержания угля даже в начальный период его разрушения очень далека от законов расширения свободного газа и может быть описана уравнениями диффузии [22], но с применением эмпирических зависимостей. Их установление является следующим шагом в создании новых инновационных методов ранжирования вида и уровня опасности, контроля эффективности мероприятий по ее снижению. Прежде всего, это оценка способности углеметанового геоматериала в процессе своего распада поддерживать высокое давление газа в стесненных условиях - как развивающейся трещины, так и формирующегося в критической ситуации углегазового потока, т.е. показателей уже кинетической энергии распада. Для этого необходимы измерения не только газоносности пласта
[23], но и скорости газовыделения из разрушенного угля в стесненных условиях, давления газа и температуры.
Выполнение указанных требований к экспериментальному оборудованию обеспечено разработанной авторами конструкцией термобарометрической колбы (рисунок 3), предназначенной для измерения газоносности угля и динамики давления и температуры выделяющегося из него газа в изолированном объеме.
Пробу газоносного угля помещают в колбу и герметизируют крышкой, имеющей соединения с газовой магистралью. Внутри колбы размещены электронные датчики давления и температуры газа и автономная микропроцессорная система для регистрации их показаний во времени. Конструкция позволяет выполнять контролируемые ступенчатые выпуски газа с периодичностью снятия показаний на каждой ступени 5.300 с без замены электропитания в течение 2 сут. Электронная система измерений разработана с учетом взрывоискробезопасности, имеет оперативную память достаточной емкости с возможностью перезаписи файлов на компьютер. Включение системы на начало регистрации магнитное.
Отметим, что в приведенных ниже результатах рассматривается не полная газоносность угля, а ее динамическая составляющая, т. е. основная часть газосодержания материала (1).
Пробы отобраны на пласте с геологической газоносностью 25 м3/т при глубине залеганий около 700 м в процессе бурения из подготовительной выработки барьерной скважины под углом 30 к борту, что позволило отобрать пять проб на близком удалении (0,65-0,9 м) от обнаженной поверхности пласта. Бурение выполнялось с промывкой, по этой причине содержание в пробах свободной воды составляло 17-19 %,
что до настоящего времени практически не исследовалось.
В таблице 1 приведены данные лабораторного анализа проб, а на рисунках 4 и 5 - графики изменения давления газа в колбах, полученные с приведением значений давления и объема выпущенной смеси к условиям: Р =760 мм рт. ст.; t =20 0С. Отметим, что в течение первых 24 ч колбы с пробами доставлялись в лабораторию. Небольшие скачкообразные изменения показаний датчиков связаны с перезарядкой аккумуляторов системы измерений и записью на компьютер содержания ее оперативной памяти. Результаты измерений смотрятся непрерывной линией в связи с их высокой плотностью. Метод обработки данных представлен на примере пробы № 4.
При обработке результатов измерений давления газа в колбах принималось, что рост давления обусловлен поступлением в свободный объем газа, выделяющегося при распаде углеметана. Скорость распада снижается как по мере газоистощения частиц, так и роста давления в колбе с переходом в пределе двухкомпонентного геоматериала в новое метастабильное состояние. Исходя из этого представления, по графику давления определялись эмпирические коэффициенты аппроксимирующей функции на основе уравнения вида
У = Утах [1-Ia exp(-ßX)], (2)
где У и Х - координатные оси; Утах - максимальное значение физической величины; ai -доля ее начального значения в соответствующем члене ряда; ß. - темп снижения в размерности Х для этого члена.
Парциальное давление метана определялось с учетом его концентрации в смеси, замеренной при выпуске. В уравнении (2) экспоненциальный ряд описывает остаточное значение рассматриваемой величины, что позволяет согласно установленным значениям эмпирических коэффициентов рассчитать, например, изменение динамической газоносности при распаде углеметана в стесненных условиях и начальную скорость газовыделения из угля после герметизации пробы, когда давление смеси в свободном объеме колбы близко к атмосферному. В качестве показателя газоотдачи принят объем газа, выделяющегося из единицы массы угля за первые 5 мин после герметизации колбы. В размерности, принятой в рудничной аэрогазодинамике, получим не менее 0,17 м3/т интенсивно увлажненного угля.
научно-технический журнал № 2-2013 А
• о
Рисунок 3 - Термобарометрическая колба для измерения газокинетических характеристик газоносного или искусственно насыщенного газом угля
Таблица 1 - Результаты измерений динамической газоносности проб угля
№ пробы Сухая беззоль-ная масса угля, г Ком- натная темпе- ратура, 0С Атмосферное давление Р мм рт.ст. Давление газа в колбе перед выпуском Р МПа Общий объем выпущ. газа Усш, см3 Концентрация метана, % Динамическая газоносность пробы Хд, см3/г
выпуск
1-й 2-й 1-й 2-й
1 176,8 24 746 0,23 1477 72,5 71,0 6,03 1,42
2 180,7 0,25 1496 74,0 82,0 6,29 1,67
3 182,5 0,21 1482 81,0 86,5 6,68 1,65
4 177,4 0,24 1588 75,0 76,0 6,73 1,60
5 146,3 0,21 1346 73,0 82,0 6,90 1,63
0,3
0,25
0,2
0,15
0,1
0,05
О
а
1
У
&
2000 4000 Т, мин 8000
1-5 номера проб
Рисунок 4 - Изменения давления газовоздушной смеси в колбах
0,2
0,15
0,1
0.05
О
1Па
г—
. •
/У л * * замер
Г ! / — — расчет
/ / и
1000
2000
3000
4000 Ті. мин
Рисунок 5 - Изменения давления газовоздушной смеси в колбе № 4 на первом интервале измерений
Следовательно, например, при отбойке взрывными работами 20 т угля с динамической газоносностью 6,73 м3/т в выработку за 5 мин непосредственно из отбитого угля поступит 3,4 м3 метана, что требует, согласно ПБ, подачи в забой не менее 34 м3/мин воздуха. Знание эмпирических коэффициентов в формуле (2) позволяет определять и интенсивность выделения метана в зоне работы комбайна с учетом его производительности. В то же время, сравнивая эти числа с графиком на рисунке 2, видим, что при сухом угле и выделении из него газа без противодавления уже за первые 30 с получаем Пг = 0,1 или
0,673 м3/т, т. е. в 4 раза больше.
Представленный результат отражает один из основных фрагментов квазистатической ме-танообильности выработок. Но для оценки уровня газодинамической активности пласта необходимы уже показатели кинетической энергии распада углеметана.
Из графика изменения давления газа в свободном (постоянном по величине) объеме колбы (рисунок 5) нетрудно получить данные о работе сжатия смеси газов, совершаемой выделяющимся из пробы метаном (рисунки 6 и 7).
Сопоставляя результаты, представленные на рисунках 1, 6 и 7, видим, что показатели по-
научно-технический журнал № 2-2013
тенциальной и кинетической энергии распада углеметана отличаются на порядок, хотя каждый из них достаточно объективно отражает уровень газодинамической активности пласта. Эта особенность указывает на существенную научную, а в перспективе и прикладную значимость разработки метода оценки кинетической энергии газовой компоненты пласта. Для реализации перспективы необходимо обеспечение оперативного контроля при проведении выработок, например эффективности дегазации пласта в установленных региональным прогнозом по показателю Е опасных зонах, формирующихся под влиянием нелинейных геомеханических процессов [24, 25]. С технической стороны оперативность метода достигается минимизацией свободного объема колбы за счет уменьшения ее объема или увеличения массы пробы (в рассматриваемом эксперименте объем колбы был в четыре раза больше объема пробы). При более полном заполнении выделяющийся из пробы метан к моменту герметизации колбы практически полностью вытеснит находившийся в ней воздух, и парциальное давление будет равно измеряемому, а электронный модуль устройства обеспечит возможность выполнения расчетов с выдачей конечного результата на дисплей колбы непосредственно в шахте.
Рисунок 6 - Изменение энергии сжатия смеси в колбе выделяющимся из пробы метаном
Рисунок 7 - Изменение отношений энергии сжатия смеси в колбе и остаточной газоносности пробы к их замеренным максимальным значениям
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Полевщиков, Г. Я. Динамические газопроявления при проведении подготовительных и вскрывающих выработок в угольных шахтах / Г. Я. Полевщиков. - Кемерово: Институт угля и углехимии СО РАН, 2003. - 317 с.
2. Плаксин, М. С. Анализ технических возможностей систем мониторинга рудничной атмосферы к непрерывному дистанционному контролю газодинамической активности приконтурной части пласта / М.С. Плаксин // Энергетическая безопасность России: труды Международной научно-практической конференции. - Кемерово, 2011. - С. 47-50.
3. Мурашев, В. И. Развитие внезапного выброса угля и газа при проведении подготовительных выработок / В. И. Мурашев, В. Л. Бульбенко // Уголь. - 1978. - № 3. - С .19-22.
4. Христианович, С. А. О волне дробления / С. А. Христианович // Известия АН СССР, ОТН. -1953. - № 12. - С. 1689-1699.
5. Христианович, С. А. О волне выброса / С. А. Христианович // Известия АН СССР, ОТН. - 1953. -№ 12. - С. 1679-1688.
6. Ходот, В. В. Внезапные выбросы угля и газа / В. В. Ходот. - М.: Госгортехиздат, 1961. - 363 с.
7. Зыков, В. С. Внезапные выбросы угля и газа и другие газодинамические явления в шахтах /
В. С. Зыков. - Кемерово: Полиграф, 2010. - 334 с.
8. Полевщиков, Г. Я. Определение параметров внезапных выбросов угля и газа / Г.Я. Полевщиков // Безопасность труда в промышленности. - 1981. - № 2. - С. 41-42.
9. Радченко, П. А. Эволюция волновых процессов и разрушения в многослойных анизотропных пластинах при ударе / П. А. Радченко, А. В. Радченко // Физика экстремальных состояний вещества.-Черноголовка: ИПХФ РАН, 2008. - С.114-116.
10. Федоров, А. В. Разлет сжатой стратифицированной концентрированной смеси / А. В. Федоров, И. В. Леонтьев // Физика горения и взрыва. - 2006. - № 2. - С.70-79.
11. Черданцев, Н. В. Геомеханическое состояние анизотропного по прочности массива в призабойной части подготовительной выработки / Н. В. Черданцев, В. Т. Преслер, В. Е. Ануфриев // Известия вузов. Горный журнал. - 2010. - № 3. - С. 33-39.
12. Gritsko, G.I. Prediction of possible parameters of sudden coal and gas outbursts on the basis of mining-and-geological data and estimation of face-space-zone’s hazardous gas emission into workings / G. I. Gritsko, G.Ja Polevshikov // Proceedings of International high gas emissions and sudden coal and gas outbursts at coal mines. - Australia, 1995. - pp. 171-175.
13. Полевщиков, Г. Я. Зоны загазования шахтной вентиляционной сети при выбросах угля и газа / Г. Я. Полевщиков, Е. С. Розанцев // Безопасность труда в промышленности. - 1979. - № 6. - С. 30-35.
14. Полевщиков, Г. Я. Устойчивость режима проветривания шахты при внезапных выбросах угля и газа / Г. Я. Полевщиков, Е. С. Розанцев // Всесоюзный семинар «Управление вентиляцией и газодинамическими явлениями в шахтах». - Новосибирск, 1979. - С. 8-12.
15. Полевщиков, Г. Я. Изучение процесса разрушения газонасыщенного угля / Г. Я. Полевщиков, Е. С. Розанцев // Физико-технические проблемы разработки полезных ископаемых. - 1979. - № 5. -
научно-технический журнал № 2-2013
ВЕСТНИК 83
С. 87-91.
16. Эттингер, И. Л. Метанонасыщенный угольный пласт как твердый метаноугольный раствор / И. Л. Эттингер // Физико-технические проблемы разработки полезных ископаемых. - 1990. - № 2. - С. 66-72.
17. Малышев, Ю. Н. Фундаментально-прикладные методы решения проблемы угольных пластов / Ю. Н. Малышев, К. Н. Трубецкой, А. Т. Айруни. - М.: ИАГН, 2000. - 519 с.
18. Полевщиков, Г. Я. Влияние параметров устойчивости системы газоносных геоматериалов на газодинамическую активность разрабатываемых пластов / Г. Я. Полевщиков, Т. А. Киряева // Проблемы рудничной аэрологии и безопасной разработки месторождений полезных ископаемых: Научно-техническая конференция к 75-летию со дня рождения профессора И.И. Медведева. - Пермь, 2004. - С. 62-67.
19. Плаксин, М. С. Оценка газодинамической активности углеметановых пластов при ведении горных работ и планирование объемов извлечения попутного метана / М. С. Плаксин, А. А. Рябцев,
B. А. Сухоруков // Вестник Научного центра по безопасности работ в угольной промышленности. - 2010. -№ 1.- С. 43-50.
20. Бирюков, Ю. М. Каталог внезапных выбросов угля и газа (Карагандинский угольный бассейн) / Ю. М. Бирюков, Р. Р Ходжаев, Н. А. Карев. - Калининград: КГТУ, 2006. - 158 с.
21. Систематизированные данные по внезапным выбросам угля и газа на шахтах восточных и северных месторождений страны. - Кемерово, 1974. - 428 с.
22. Назарова, Л. А. Определение коэффициента диффузии и содержания газа в пласте на основе решения обратной задачи по данным измерения давления в герметичной емкости с угольным веществом / Л. А. Назарова, Л. А. Назаров, Г. Я. Полевщиков, Р И. Родин // Физико-технические проблемы разработки полезных ископаемых. - 2012. - № 5. - С. 15-23.
23. Тайлаков, О. В. Разработка методики определения газоносности угольных пластов в процессе ведения горных работ / О. В. Тайлаков, Е. В. Мазаник, А. Н. Кормин // Горный информационно-аналитический бюллетень. Экология, метанобезопасность. - 2011. - № 8. - С. 210-214.
24. Козырева, Е. Н. Газокинетические следствия нелинейных геомеханических процессов в массиве горных пород на шахтах Кузбасса / Е. Н. Козырева, М. В. Шинкевич, Р. И. Родин // Нелинейные гео-механико-геодинамические процессы при отработке месторождений полезных ископаемых на больших глубинах: Труды 2-й Российско-китайской научной конференции. - Новосибирск: ИГД СО РАН, 2012. -
C. 267-272.
25. Козырева, Е. Н. Особенности газогеомеханических процессов на выемочном участке шахты / Е. Н. Козырева, М.В. Шинкевич // Вестник Научного центра по безопасности работ в угольной промышленности. - 2010. - № 2. - С. 28-35.
DETERMINATION OF GAS-KYNETIC Полевщиков Геннадий Яковлевич
PROPERTIES OF COAL SEAMS е-mail: [email protected]
G.Ya. Polevshchikov, A. A. Riabtsev, V. P. Titov
Scientific-technical grounds are developed to Рябцев Андрей Александрович
improve evaluation method of gas-dynamic activity е-mail: [email protected]
levels of coal seams from quasi-static methane
content of the openings to exactly sudden outbursts Титов Виктор Петрович
of coal and gas as a row of genetically common е-mail: [email protected]
gas manifestations. Results of on-site observations
and their comparison on the basis of potential and
kinetic energy data of coal-methane seams gas
component are given.
The work is fulfilled with financial support of
interdisciplinary integration project of SB RASc No
99 and partnership integration project of SB RASc
No 100.
Key words: COAL SEAM, GAS CONTENT,
METHANE EMISSION, DESTRUCTION,
TRANSPORTATION, SUDDEN OUTBURST,
ENERGY DATA
А научно-технический журнал № 2-2013
о4 вестник
