CC BY
14
Оригинальная статья
УДК 533.6:622.4 © Г.А. Колегов, А.Ю. Крайнов, 2021
Метод учета аэродинамических параметров
выработанных пространств в моделях
*
шахтных вентиляционных систем
DOI: http://dx.doi.org/10.18796/0041-5790-2021-7-33-38
КОЛЕГОВ Г.А.
Аспирант кафедры математической физики ФГАОУ ВО «Национальный исследовательский Томский государственный университет» (ТГУ), 634050, г. Томск, Россия, e-mail: zergferr@gmail.com
КРАЙНОВ А.Ю.
Заведующий кафедрой математической физики ФГАОУ ВО «Национальный исследовательский Томский государственный университет» (ТГУ), 634050, г. Томск, Россия, e-mail: akrainov@ftf.tsu.ru
Выработанные пространства угольных шахт являются важным элементом рудничныхвентиляционных сетей. В статье предложен метод включения выработанных пространств угольных шахт в математические модели шахтных вентиляционных систем посредством представления пористой среды сетью фиктивных ветвей, учитывающих свойства элементарных репрезентативных объемов. Аэродинамические параметры фиктивных ветвей предлагается вычислять с использованием квадратичных сопротивлений, соблюдая законы Кирхгофа. Обоснован и верифицирован способрасче-та аэродинамических сопротивлений ветвей выработанных пространств согласно формуле Эргуна. Разработанный метод задания аэродинамического сопротивления выработанного пространства был применен для оценки градиентов давления в выбранном выработанном пространстве шахты «Распадская», для которого имеются измеренные значения перепадов давления. Приведены результаты расчетов. Определена практическая значимость полученных результатов, обозначены пути дальнейшего развития методологии. Ключевые слова: вентиляция, вентиляционная сеть, математическая модель, выработанное пространство, аэродинамическое сопротивление. Для цитирования: Колегов Г.А., Крайнов А.Ю. Метод учета аэродинамических параметров выработанных пространств в моделях шахтных вентиляционных систем // Уголь. 2021. № 7. С. 33-38. 001: 10.18796/0041-5790-20217-33-38.
ВВЕДЕНИЕ
Математическое моделирование является неотъемлемой частью современных производственных процессов. Моделирование шахтовой вентиляции решает задачи как организации обеспечения объектов проветривания требуемыми расходами воздуха, так и разработки планов ликвидации аварий. От качества моделей зависит эффективность мероприятий, разрабатываемых на их основе. Качество определяется условиями и объектами, учитываемыми моделями, а также соответствием вводных данных и получаемых результатов фактически наблюдаемой ситуации в шахте. Помимо этого, инструменты и методы моделирования, широко применяемые на горнодобывающих предприятиях, должны обладать достаточной простотой при сохранении определенного уровня точности и не быть требовательными к вычислительным мощностям. Программный комплекс «Вентиляция 2» [1], наиболее распространенный на Российских предприятиях подземной угледобычи, соответствует этим требованиям, позволяя пользователю решать широкий спектр задач аэродинамики в масштабах шахтной вентиляционной сети (ШВС). В программном обеспечении (ПО) [1] вентиляционные сети представлены связным ориентированным графом, в котором соблюдаются законы сетей (нулевая сумма входящих в узел и исходящих из него объемных расходов и нулевая сумма депрессий замкнутого контура), а также закон сохранения энергии. При этом депрессия 1-й ветви в контуре определяется в соответствии с выражением:
(1)
где Q. - объемный расход в ветви /', м3/с; Я. - аэродинамическое сопротивление ветви /', кмюрг; И. - перепад давления (депрессия) в ветви I, даПа.
Выражение (1) называют квадратичным законом сопротивления, описывающим движение воздушных потоков по горным выработкам при развитом турбулентном режиме [1]. В таком подходе величину сопротивления определяют геометрические формы каналов, и не учитываются вязкие трения слоев потока. Аэродинамические со-
* Работа выполнена при финансовой поддержке Минобрнауки России в рамках выполнения государственного задания № 0721-2020-0036.
противления определяются прямыми замерами во время воздушно-депрессионных съемок, а также расчетным путем на основании экстраполяции коэффициентов аэродинамических сопротивлений, полученных эмпирическим путем для характерных типов выработок.
Выработанные пространства угольных шахт, применяющих столбовую систему разработки с полным обрушением, представляют собой пористые среды, в которых преобладают фильтрационные режимы движения газов [2]. Определение аэродинамических параметров выработанных пространств является задачей, решавшейся с момента возникновения рудничной вентиляции как научного направления [3]. На современном этапе преобладает подход к их расчету с использованием средств вычислительной гидродинамики [4, 5, 6, 7, 8]. Данный подход обеспечивает точность вычислений, учитывает взаимодействие выработанного пространства с ШВС в виде данных, задаваемых на основании одноразовых измерений и не позволяет оперативно учитывать изменяющиеся условия в шахте [6]. Существуют исследования, позволяющие учитывать аэродинамику в выработанных пространствах, рассчитываемую методами конечных объемов, в модели ШВС [9, 10]. Однако такие подходы требуют значительных вычислительных мощностей. В настоящей работе представлен метод расчета аэродинамики в выработанном пространстве совместно с расчетом аэродинамики проветривания сети выработок шахты.
ОПИСАНИЕ МЕТОДА
Для возможности включения выработанного пространства в существующие средства моделирования ШВС необходимо определить метод описания пористой среды выработанных пространств ветвями связных графов. Для перехода от трехмерной задачи фильтрации потоков через каналы переменного сечения и пространственной ориентации обрушенной среды используем подход репрезентативного элементарного объема, описанный Теруелом и Аддином [11], представленный на рис. 1.
Полагается, что зона Ь формируется из участков В, характеризующихся вариациями размеров структурных элементов выработанного пространства и расстояний между ними р. Предположим, что один элементарный участок В с усредненными по нему характеристиками пористости и эквивалентного диаметра структурного элемента может быть представлен ветвью в модели ШВС.
В работе [4] сопротивление движению потока в остаточных каналах выработанного пространства задано формулой Эргуна:
В
Рис. 1. Формирование зоны выработанного пространства из элементарных объемов
Ар _ 150г|и(1-£) 1.75ри2(1-е)
Т~ Ц? + 7? ,
где п - динамическая вязкость, Па-с; р - плотность, кг/м3; и - скорость фильтрации, м/с; е- пористость; - эквивалентный диаметр структурного элемента, м.
Первое слагаемое в выражении обозначим линейным удельным падением давления, а второе - квадратичным удельным падением давления. Предположим, что выражение (2) возможно использовать для расчета интегральных аэродинамических сопротивлений фиктивных ветвей выработанного пространства в квазистационарной постановке задачи с целью включения этих ветвей в общую модель вентиляционной сети. Все средства моделирования шахтной вентиляции используют в настоящий момент квадратичные сопротивления. Предполагая типичные условия эксплуатации шахт, оценим влияние каждой переменной в уравнении Эргуна на результирующие удельные потери давления при фиксированных значениях других величин. Динамическая вязкость и плотность изменяются в зависимости от температуры при постоянном давлении. Динамическая вязкость определяется выражением: з
Л = Ло
т0+с
т+с
(3)
'о у
где п0 - величина вязкости при температуре Т0, Па-с; Т -температура, К; Т0 - контрольная температура, Т0= 291,15 К; С - постоянная Сазерленда, для воздуха С = 120 К.
Плотность изменяется в зависимости от температуры в соответствии с уравнением состояния идеального газа:
Р =
рМ ЯТ'
(4)
гдер - давление, Па;М- молярная масса, кг/моль; Я - универсальная газовая постоянная, Дж/(моль-К).
Проведем оценку вклада первого слагаемого в выражении (2) в общую величину удельного падения давления Запишем отношение первого слагаемого ко второму, получим выражение:
150т](1 — е)
А =
(5)
1,75^эр и
Оценим максимальную величину (5) в интервалах изменения переменных, входящих в (5): 10 < Т < 290С; 0,1 < и < 1 м/с; 0,1 < е < 0,9; 0,2 < < 5,0 м. Такие интервалы выбраны исходя из характерных величин для обрушенного пространства [12, 13]. При изменении температуры в выбранном интервале вязкость газа (3) изменяется в интервале 1,787 • 10"5 < т] < 1,881 ■ 10"5 Па-с , плотность газа (4) 1,235 > р > 1,157 кг/м3. Выбирая крайние значения величин из выбранных интервалов их изменения, которые приводят к получению максимального значения ве-
личины А, получим: А =
150-1,881-Ю"5 (1-0,1)
= 0,0627.
1,75-0,2-1,157-0,1
Первое слагаемое выражения (2) в интервале выбранных значений параметров обрушенной среды не превышает 6,27%. Изменение температуры через изменение динамической вязкости и плотности влияет на изменение удельных перепадов давления в наименьшей степени.
Ь
Поэтому принимаемый в рамках модели всей ШВС изотермический режим позволяет без значительной для конечного результата погрешности рассматривать эти величины постоянными. Следовательно, в рамках общей ШВС, в которой используется квадратичный закон сопротивления, формулу Эргуна можно использовать в следующем виде:
Др_1,75ри2(1-е)
Т" ~сйг .
(6)
Предполагая постоянными плотность газовой смеси, а также скорость в фиктивных ветвях, получим удельное падение давления, зависящее от параметров выработанного пространства: приведенной пористости и эквивалентного диаметра структурных элементов. На рис. 2 приведена зависимость удельных квадратичных падений давления от эквивалентного диаметра (ДР) при постоянной пористости и падений давления от пористости при постоянном эквивалентном диаметре (ДРЕ ).
Из рис. 2 следует, что в максимально зажатой среде больший вклад в изменение давления вносит пористость, при увеличении пористости и эквивалентного диаметра эквивалентный диаметр вносит более значимый вклад.
Разделив выражение (6) на квадрат расхода воздуха и переведя перепад давлений в декапаскали (даПа), получим формулу для вычисления удельного сопро-
0 2
й (м)
4 6
10
200
150
< 100
Д
50
0
320
240
160
Д
80
0
0,6
6
тивления фиктивных ветвей выработанного пространства, зависящего от распределения пористости вдоль выемочного столба:
Д.,
2,8р(1-е)
(7)
71 "эе
Такая форма выражения сопротивления пригодна для использования в программных продуктах, использующих квадратичный закон перепадов давлений, в частности, в ПО «Вентиляция». Предполагается задание сетки с шагом, меньшим характерного изменения величины пористости вдоль выемочного столба. Пористость возможно вычислять на основании экспорта из специализированного геомеханического ПО наподобие РЬАСЭР либо вычислять на основании изменения коэффициента проницаемости, определяемого по соответствующим методикам [4, 11, 12, 13].
АПРОБАЦИЯ И ВЕРИФИКАЦИЯ МЕТОДИКИ
Разработанный метод задания аэродинамического сопротивления выработанного пространства был применен для оценки градиентов давления в выработанных пространствах выемочных столбов, включенных в общую модель ШВС. Для выработанного пространства лавы № 5а-7-34 шахты «Распадская» было сделано предположение о равномерном распределении пористости по ширине столба и переменным по его длине. Удельные сопротивления рассчитывались для каждого значения пористости, присвоенного фиктивным ветвям расчетной сети на основании (7) и приведены в табл. 1.
Шаг сетки принимался равным около 30 м. На рис. 3 приведено принятое распределение пористости выработанного пространства, рис. 4 иллюстрирует результирующее
Таблица 1
Значения удельных сопротивлений, используемые в расчетной сети выработанного пространства лавы № 5а-7-34
Рис. 2. Изменение квадратичного удельного падения давления
Пористость Е Эквивалентный диаметр й, м , КМ
0,6 3,0 0,002636
0,5 3,0 0,005694
0,2 2,3 0,543007
Рис. 3. Моделируемое распределение пористости вдоль выемочного столба 5а-7-34
8
Рис. 4. Распределение избыточных давлений вдоль выемочного столба 5а-7-34
распределение избыточного давления в выработанном пространстве.
Как следует из результатов, вслед за пористостью, падение давления распределено равномерно вдоль отработанной части выемочного столба. При этом избыточное давление в действующих выработках варьируется от 2371 до 1172 Па, а давление внутри столба падает от 2317,6 до 1963,1 Па.
Единственным способом контроля расчетных характеристик выработанного пространства являться сопоставление расчетных перепадов давления на изоляционных перемычках в сбойках с фактическими значениями, предполагая, что распределение перепадов давлений в действующих выработках также соответствует действительности. Изменения перепадов давлений на перемычках также можно использовать для оценки реакции моделируемой пористой обрушенной среды на изменения во внешнем контуре, а также для моделирования связей выработанного пространства с атмосферой или выработками, пройденными по пластам-спутникам. Полученные перепады давлений на перемычках и расходы, а также результаты расчета депрессий на перемычках в модели и показания датчиков СКПА, усредненные за неделю, приведены в табл. 2.
Таблица 2
Распределение перепадов давлений и расходов на контролируемых перемычках
Название Расход, м3/с Депрессия из расчета по модели, Па Депрессия на датчиках СКПА, Па
Сбойка №4 0,24 130 100
Сбойка №6 0,50 558 600
Сбойка №8 0,44 435 400
Одновременно расходы служат способом контроля направления утечки - при его изменении расходы станут отрицательными. Для верификации полученных в модели ШВС значений использовались данные системы станций контроля параметров атмосферы (СКПА), применяющейся на шахте «Распадская». На рис. 5 приведены изменения перепада давления на перемычках в сбойках № 4, № 6 и № 8 между вентиляционным штреком № 5а-7-36 и конвейерным штреком № 5а-7-34.
В рамках погрешностей датчиков СКПА можно говорить о совпадении расчетных значений с усредненными по времени данными о фактических перепадах давлений на контролируемых перемычках.
Р, Па
СБ,
Сб. №4
Сб. №8
1600
иоо
ноо
1000
800
600
400
200
лг V
/ ТГ 41
/ \
/ к
^ _____т
1 Л "Т ■
иг^
Рис. 5. Перепады давлений на контролируемых СКПА перемычках, Па
01.03.2021
02.03.2021
03.03.2021
04.03.2021
05.03,2021
Дата
ВЫВОДЫ
Предложен метод включения выработанных пространств угольных шахт в математические модели шахтных вентиляционных систем посредством представления пористой среды сетью фиктивных ветвей, учитывающих свойства элементарных репрезентативных объемов. Обоснован и верифицирован способ расчета аэродинамического сопротивления ветвей выработанных пространств.
Использование подхода разделения пористой среды выработанного пространства на репрезентативные элементарные объемы позволяет моделировать аэродинамику в выработанном обрушенном пространстве с использованием ПО «Вентиляция».
Численный анализ формулы Эргуна демонстрирует вклад сил вязкости, не превышающий 6,5% в суммарные удельные потери давления при скоростях фильтрации, составляющих десятки сантиметров в секунду, что позволяет использовать лишь квадратичное слагаемое для оценки градиентов давления в сети фиктивных ветвей выработанного пространства. Факторами, определяющими сопротивление выработанного пространства, являются распределения пористостей и эквивалентных диаметров структурных элементов.
При включении выработанных пространств в модели шахтных вентиляционных систем верифицируемыми индикаторами адекватности расчетов являются перепады давлений на перемычках, отделяющих выработанное пространство от действующих выработок.
Список литературы
1. Палеев Д.Ю., Лукашов О.Ю. Программа расчета вентиляционных режимов в шахтах и рудниках // Горная промышленность. 2007. № 6. С. 20-23.
2. Ушаков К.З. Аэромеханика вентиляционных потоков в горных выработках. М.: Недра, 1975. 167 с.
3. Скочинский А.А., Комаров В.Б. Рудничная вентиляция. М.-Л.: Углетехиздат, 1949. 442 с.
4. Кобылкин С.С. Обоснование метода расчета параметров вентиляции шахт на основе объемного моделирования аэрогазодинамических процессов: дис. ... канд. техн. наук: спец. 05.26.03: защищена 2011 / С.С. Кобылкин. М., 2011. 161 с.
5. Говорухин Ю.М., Фрянов В.Н., Палеев Д.Ю. Численное моделирование взаимодействующих геомеханических и аэродинамических процессов на выемочном участке // Известия Тульского государственного университета. Науки о Земле. 2019. Вып. 2. С. 288-299.
6. Gilmore R.C. Computational fluid dynamics modeling of underground coal longwall gob ventilation systems using a developed meshing aproach. Dr. Mechanical Engineering diss. Colorado School of Mines, 2015. 280 p. URL: http://hdl.handle. net/11124/20176 (дата обращения: 15.06.2021).
7. Lolon S. Computational Modeling of Barometric Pressure Fluctuation Effects on Explosive Methane-Air Mixtures in a Longwall Mine Gob. PhD (Mining and Earth Systems Engineering) diss., 2017. 190 p. DOI: http://dx.doi.org/10.13140/ RG.2.2.22535.19360.
8. Jing Shen, Mingran Chang. Effect of nitrogen ingection on the variation of oxidation in coal mine gob based on numerical simulation // The Open Fuels and Energy Science Journal.
2016. Vol. 9. P. 47-54.
9. Wedding W.C. Multiscale modeling of mine vintilation system and flow through the gob. Dr. of Philosophy diss. Lexington, Kentucky, 2014. 249 p. URL: https://uknowledge.uky. edu/cgi/viewcontent.cgi?article=1014&context=mng_etds (дата обращения: 15.06.2021).
10. Wu F-I, Luo Y, Chang X-t. Coupling simulation model between mine ventilation network and gob flow field. The Southern African Institute of Mining and Metallurgy, 2019. URL: https://www.researchgate.net/publication/337580509 (дата обращения: 15.06.2021).
11. Teruel E.F., Uddin R. Macroscopic modeling of turbulence in porous media flows, Mecánica Computacional Vol. XXVII, págs. San Luis, Argentina, 2008. P. 265-280.
12. Long Fan, Shimin Liu. A conceptual model to characterize and model compaction behavior and permeability evolution in coal mine gobs // International Journal of Coal Geology.
2017. Vol. 172. P. 60-70.
13. Петушкеев Б.Л. Решение трехмерной задачи газовой динамики и переноса метана в угольной шахте с использованием параллельных вычислений. Дис. ...канд. физ.-мат. наук: спец. 01.02.05: защищена 2010 / Б.Л. Петушкеев. Томск, 2010. 150 с.
SAFETY
Original Paper
UDC 533.6:622.4 © G.A. Kolegov, A.Yu. Krainov, 2021
ISSN 0041-5790 (Print) • ISSN 2412-8333 (Online) • Ugol' - Russian Coal Journal, 2021, № 7, pp. 33-38 DOI: http://dx.doi.org/10.18796/0041-5790-2021-7-33-38
Title
METHOD TO ACCOUNT FOR AERODYNAMIC PARAMETERS OF MINED-OUT SPACES IN MODELS OF MINE VENTILATION SYSTEMS
Authors
Kolegov G.A.1, Krainov A.Yu.1
'Tomsk State University, Tomsk, 634050, Russian Federation Authors Information
Kolegov G.A., Postgraduate student of Mathematical physics department, e-mail: zergferr@gmail.com
Krainov A.Yu., Head of Mathematical physics department, e-mail: akrainov@ftf.tsu.ru
Abstract
Underground coal mines' gob regions are an integral part of the ventilation process. The paper proposes an approach to incorporating isolated gob regions within longwall coal mines into mine ventilation network models by ascribing elementary representative volumes' properties to virtual edges
to simulate porous media. Aerodynamic properties of the virtual edges are determined by their quadratic resistances in accordance with Krchhoff's laws. A way to calculate the aerodynamic resistances of the virtual edges using Ergun's formula is determined and verified by comparing the results with factual pressure differences on stoppings in the «Raspadskaya» underground mine. The practical significance of the results is considered and the ways to improve the methodology in the future are proposed.
Keywords
Ventilation, Mine ventilation network, Mathematical model, Gob, Aerodynamic resistance.
References
1. Paleev D.Yu. & Lukashov O.Yu. Software for calculation of ventilation modes in soft-rock and hard-rock underground mines. Gornaya promyshlennost', 2007, (6), pp. 20-23. (In Russ.).
2. Ushakov K.Z. Aeromechanics of ventilation air flows in mine excavations. Moscow, Nedra Publ., 1975, 167 p. (In Russ.).
3. Skochinsky A.A. & Komarov V.B. Mine Ventilation. Moscow-Leningrad, Ug-lechtekhizdat Publ., 1949, 442 p. (In Russ.).
4. Kobylkin S.S. Justification of a method to calculate mine ventilation parameters based on 3D modeling of aerogasdynamic processes. PhD (Engineering) diss. Moscow, 2011, 161 p. (In Russ.).
5. Govorukhin Yu.M., Fryanov V.N. & Paleev D.Yu. Numerical modeling of interacting geomechanical and aerodynamic processes in the extraction zone. Izvestiya Tul'skogo gosudarstvennogo universiteta, Nauki o Zemle, 2019, (2), pp. 288-299. (In Russ.).
6. Gilmore R.C. Computational fluid dynamics modeling of underground coal longwall gob ventilation systems using a developed meshing aproach. Dr. Mechanical Engineering diss. Colorado School of Mines, 2015, 280 p. Available at: http://hdl.handle.net/11124/20176 (accessed 15.06.2021).
7. Lolon S. Computational Modeling of Barometric Pressure Fluctuation Effects on Explosive Methane-Air Mixtures in a Longwall Mine Gob. PhD (Engineering) diss., 2017, 190 p. DOI: http://dx.doi.org/10.13140/RG.2.2.22535.19360 .
8. Jing Shen & Mingran Chang. Effect of nitrogen ingection on the variation of oxidation in coal mine gob based on numerical simulation. The Open Fuels and Energy Science Journal, 2016, (9), pp. 47-54.
9. Wedding W.C. Multiscale modeling of mine vintilation system and flow through the gob. Dr. of Philosophy diss. Lexington, Kentucky, 2014, 249 p. Available at: https://uknowledge.uky.edu/cgi/viewcontent. cgi?article=1014&context=mng_etds (accessed 15.06.2021).
10. Wu F-I., Luo Y. & Chang X-t. Coupling simulation model between mine ventilation network and gob flow field. The Southern African Institute of Mining and Metallurgy, 2019. Available at: https://www.researchgate.net/ publication/337580509 (accessed 15.06.2021).
11. Teruel E.F. & Uddin R. Macroscopic modeling of turbulence in porous media flows, Mecánica Computacional, Vol. XXVII, págs. San Luis, Argentina, 2008, pp. 265-280.
12. Long Fan & Shimin Liu. A conceptual model to characterize and model compaction behavior and permeability evolution in coal mine gobs. International Journal of Coal Geology, 2017, (172), pp. 60-70.
13. Petushkeev B.L. Solving a three-dimensional gas dynamics and methane transfer problem in a coal mine using parallel computing. PhD (Engineering) diss. Tomsk, 2010, 150 p. (In Russ.).
Acknowledgments
This work was financially supported by the Ministry of Education and Science of the Russian Federation as part of State Assignment No. 0721-2020-0036.
For citation
Kolegov G.A. & Krainov A.Yu. Method to account for aerodynamic parameters of mined-out spaces in models of mine ventilation systems. Ugol, 2021, (7), pp. 33-38 (In Russ.). DOI: 10.18796/0041-5790-2021-7-33-38.
Paper info
Received March 24,2021 Reviewed June 06,2021 Accepted June 15,2021
Холдинг «СДС-Уголь» выпустил в реки Кузбасса мальков редких пород рыб
УГОЛЬ
Полугодовалых мальков хариуса сибирского, тайменя, ленка тупорылого в рамках программы восполнения водных биологических ресурсов выпустили в приток реки Томь экологи предприятий АО ХК «СДС-Уголь».
В мероприятии, которое прошло 20 июня, приняли участие губернатор Кузбасса Сергей Цивилёв и президент Холдинга «Сибирский Деловой Союз», председатель Совета директоров АО ХК «СДС-Уголь» Михаил Федяев. Организаторами выступили АО ХК «СДС-Уголь» и КАО «Азот».
Работы по искусственному воспроизводству популяции исчезающих видов рыб проводятся угольным Холдингом в теплое время года совместно с инспекторами Федерального агентства по рыболовству. Начиная с 2018 г. АО ХК «СДС-Уголь» выпустило в кузбасские реки около 600 тыс. мальков.
300^1
КУЗБАСС
«Стоит отметить, что мы не только выпускаем в реки мальков рыб, которые занесены в Красную Книгу, -рассказывает Любовь Тургенева, начальник управления экологической безопасности и охраны окружающей среды АО ХК «СДС-Уголь». - Мы принимаем необходимые меры для того, чтобы наши реки стали чище, регулярно проводим экологические акции по очистке водоемов. Это системная работа по снижению негативного воздействия на природу Кузбасса, проводимая в рамках единой ESG-политики Холдинга «СДС-Уголь» и концепции «Чистый уголь - зеленый Кузбасс».
Сохранение экологического баланса, выполнение компенсационных мер, минимизация негативного воздействия на природу являются приоритетными направлениями для АО ХК «СДС-Уголь». С 2018 г. АО ХК «СДС-Уголь» занимает одну из ведущих позиций в рейтинге экологической ответственности горнодобывающих и металлургических компаний Всемирного фонда дикой природы (WWF).
Чистый уголь - ЗЕЛЕНЫЙ КУЗБАСС
