Оценка изменения термодинамических свойств подземных вод при подземном сжигании угля

Розроблена математична модель теплоперенесення в пі-дземних водах покриваючих порід вугільного пласта при його спа-люванні. Для гірничо-геологічних умов ділянки "Ольхова Нижня" за допомогою програми чисельного моделювання ModFlow 2009 отримані параметри зміни розмірів і форми теплових аном...

Повний опис

Збережено в:

Бібліографічні деталі
Дата:	2014
Автори:	Садовенко, И.А., Инкин, А.В.
Формат:	Стаття
Мова:	Russian
Опубліковано:	Український науково-дослідницький і проектно-конструкторський інститут гірничої геології, геомеханіки і маркшейдерської справи НАН України 2014
Назва видання:	Наукові праці УкрНДМІ НАН України
Онлайн доступ:	http://dspace.nbuv.gov.ua/handle/123456789/99532
Теги:	Додати тег Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:	Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:	Оценка изменения термодинамических свойств подземных вод при подземном сжигании угля / И.А. Садовенко, А.В. Инкин // Наукові праці УкрНДМІ НАН України. — 2014. — № 14. — С. 48-59. — Бібліогр.: 7 назв. — рос.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine

id	irk-123456789-99532
record_format	dspace
spelling	irk-123456789-995322016-05-01T03:02:02Z Оценка изменения термодинамических свойств подземных вод при подземном сжигании угля Садовенко, И.А. Инкин, А.В. Розроблена математична модель теплоперенесення в пі-дземних водах покриваючих порід вугільного пласта при його спа-люванні. Для гірничо-геологічних умов ділянки "Ольхова Нижня" за допомогою програми чисельного моделювання ModFlow 2009 отримані параметри зміни розмірів і форми теплових аномалій, що формуються у водопроникних породах. A mathematical model of heat transport in groundwater which saturate the overlying rocks of the coal seam by burning it. For the geological site conditions "Olhova Nizhnyaya" using simulation ModFlow 2009 estimates of changes in the size and shape of thermal anomaly formed in permeable rocks above the roof of the coal seam, depending on the angle of his fall, and the stage of combustion. 2014 Article Оценка изменения термодинамических свойств подземных вод при подземном сжигании угля / И.А. Садовенко, А.В. Инкин // Наукові праці УкрНДМІ НАН України. — 2014. — № 14. — С. 48-59. — Бібліогр.: 7 назв. — рос. 1996-885X http://dspace.nbuv.gov.ua/handle/123456789/99532 622.278 ru Наукові праці УкрНДМІ НАН України Український науково-дослідницький і проектно-конструкторський інститут гірничої геології, геомеханіки і маркшейдерської справи НАН України
institution	Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection	DSpace DC
language	Russian
description	Розроблена математична модель теплоперенесення в пі-дземних водах покриваючих порід вугільного пласта при його спа-люванні. Для гірничо-геологічних умов ділянки "Ольхова Нижня" за допомогою програми чисельного моделювання ModFlow 2009 отримані параметри зміни розмірів і форми теплових аномалій, що формуються у водопроникних породах.
format	Article
author	Садовенко, И.А. Инкин, А.В.
spellingShingle	Садовенко, И.А. Инкин, А.В. Оценка изменения термодинамических свойств подземных вод при подземном сжигании угля Наукові праці УкрНДМІ НАН України
author_facet	Садовенко, И.А. Инкин, А.В.
author_sort	Садовенко, И.А.
title	Оценка изменения термодинамических свойств подземных вод при подземном сжигании угля
title_short	Оценка изменения термодинамических свойств подземных вод при подземном сжигании угля
title_full	Оценка изменения термодинамических свойств подземных вод при подземном сжигании угля
title_fullStr	Оценка изменения термодинамических свойств подземных вод при подземном сжигании угля
title_full_unstemmed	Оценка изменения термодинамических свойств подземных вод при подземном сжигании угля
title_sort	оценка изменения термодинамических свойств подземных вод при подземном сжигании угля
publisher	Український науково-дослідницький і проектно-конструкторський інститут гірничої геології, геомеханіки і маркшейдерської справи НАН України
publishDate	2014
url	http://dspace.nbuv.gov.ua/handle/123456789/99532
citation_txt	Оценка изменения термодинамических свойств подземных вод при подземном сжигании угля / И.А. Садовенко, А.В. Инкин // Наукові праці УкрНДМІ НАН України. — 2014. — № 14. — С. 48-59. — Бібліогр.: 7 назв. — рос.
series	Наукові праці УкрНДМІ НАН України
work_keys_str_mv	AT sadovenkoia ocenkaizmeneniâtermodinamičeskihsvojstvpodzemnyhvodpripodzemnomsžiganiiuglâ AT inkinav ocenkaizmeneniâtermodinamičeskihsvojstvpodzemnyhvodpripodzemnomsžiganiiuglâ
first_indexed	2025-07-07T08:06:29Z
last_indexed	2025-07-07T08:06:29Z
_version_	1836974699522818048
fulltext	Наукові праці УкрНДМІ НАН України, № 14, 2014 Transactions of UkrNDMI NAN Ukraine, № 14, 2014 48 УДК 622.278 ОЦЕНКА ИЗМЕНЕНИЯ ТЕРМОДИНАМИЧЕСКИХ СВОЙСТВ ПОДЗЕМНЫХ ВОД ПРИ ПОДЗЕМНОМ СЖИГАНИИ УГЛЯ Садовенко И. А., Инкин А. В. (Государственное ВУЗ "НГУ", г. Днепропетровск, Украина) Розроблена математична модель теплоперенесення в пі- дземних водах покриваючих порід вугільного пласта при його спа- люванні. Для гірничо-геологічних умов ділянки "Ольхова Нижня" за допомогою програми чисельного моделювання ModFlow 2009 отримані параметри зміни розмірів і форми теплових аномалій, що формуються у водопроникних породах. A mathematical model of heat transport in groundwater which saturate the overlying rocks of the coal seam by burning it. For the geological site conditions "Olhova Nizhnyaya" using simulation ModFlow 2009 estimates of changes in the size and shape of thermal anomaly formed in permeable rocks above the roof of the coal seam, depending on the angle of his fall, and the stage of combustion. Введение. По оценкам ведущих нефтегазовых компаний, начало ХХI века ознаменовалось значительным увеличением до- ли угля в топливно-энергетическом балансе ряда стран мира [1]. Рост его потребления был вызван, с одной стороны, ограничен- ностью запасов нефти и природного газа, а также перемещением разрабатываемых месторождений в труднодоступные районы, с другой стороны, проблемами безопасности атомных электро- станций и высокими затратами на захоронение радиоактивных отходов. Вместе с тем, увеличение добычи и переработки угля тради- ционными способами привело к превращению угледобывающих регионов в зоны экологического бедствия. Так, работа небольшой Наукові праці УкрНДМІ НАН України, № 14, 2014 Transactions of UkrNDMI NAN Ukraine, № 14, 2014 49 угольной электростанции мощностью 200 МВт сопровождается ежегодным выбросом в атмосферу 100 тыс. т твердых частиц, 15 тыс. т сернистых соединений и 2 тыс. т окислов азота [2]. Расширение сферы использования угольных ресурсов воз- можно путем их преобразования на месте залегания в экологиче- ски более чистый газообразный энергоноситель. Наиболее техни- чески освоенным процессом такой переработки является подзем- ное сжигание угля (ПСУ). Однако, несмотря на почти вековой научно-промышленный опыт развития, оно до сих пор не полу- чило широкого распространения ввиду значительного рассеива- ния вырабатываемой тепловой энергии в породном массиве [3]. Постановка задачи исследований. Для повышения КПД данной технологии в работе [4] уже была установлена возмож- ность отбора и количественная характеристика теплового потока поступающего из реакционного канала в выше залегающий водо- носный горизонт. Целью данной работы является количествен- ный анализ пространственно-временной динамики формирования гидрогеотермических полей, образующихся в водоносном пласте при сжигании угля. Это позволит обосновать рациональные тех- нологические параметры извлечения тепла из откачиваемых нагретых вод. Материал и результаты исследований. Двумерная в плане миграция тепла в подземных водах описывается уравнениями, приведенными в работе [5]. Принципиально важным этапом раз- работки модели является корректная формулировка граничного условия по температуре на подошве водоносного горизонта над реакционным каналом. Для определения температуры подземных вод в данной части водоносного пласта выделяется объемный блок малой толщины в виде параллелепипеда, расположенный непосредственно над нагретым водоупором, где происходит теп- лообмен (рис. 1). Баланс тепла в блоке определяется на основа- нии равенства потока тепла снизу от реакционного канала (Q0) трем составляющим: конвективному выносу тепла фильтрацион- ным потоком (Q2), убыли тепла в результате замещения теплой воды холодной (Q1), кондуктивного потока тепла вверх (Q3), а также тепла, затрачиваемого на изменение температуры подзем- ных вод и горных пород непосредственно в блоке. Наукові праці УкрНДМІ НАН України, № 14, 2014 Transactions of UkrNDMI NAN Ukraine, № 14, 2014 50 При учёте баланса тепла в блоке приняты следующие допу- щения: 1) баланс тепла определяется для последовательности пе- риодов (например, длительностью 1 сут), в течение каждого из которых теплофизические параметры осредняются; 2) боковая теплоотдача из блока не учитывается, поскольку она на порядок меньше теплоотдачи через его верхнюю границу; 3) коэффициенты граничного условия задачи теплопереноса определяются на основе предварительно рассчитанного поля ско- ростей фильтрации; 4) температура в блоке осредняется по ее начальным и конечным значениям в каждый период; 5) температура воды и породы в блоке одинакова и выравнивает- ся условно мгновенно при втекании и вытекании воды с разной температурой. Рис. 1. Схема баланса тепла в блоке водоносного пласта над кровлей реакционного канала при сжигании угля При сделанных допущениях изменение температуры воды и пород в блоке определяется уравнением баланса тепла: BTTQQQQ ⋅−=−−+ )( 013210 , (1) где пппwww VсVcB ρρ += , (2) yxqQ ∆∆=0 , (3) cwATQ =1 , cwww TcVzyA ρτ ⋅⋅⋅∆⋅∆= , (4) Tcw, Q1 T0, T1 Q2 Q3 Δ Δ x Δ y Водоупор Δ z Q0 Наукові праці УкрНДМІ НАН України, № 14, 2014 Transactions of UkrNDMI NAN Ukraine, № 14, 2014 51 2 01 2 TT AQ + ⋅= , (5) ) 2 () 2 ( 0101 3 cwcw T TT DT TT z yxQ − + ⋅=− + ∆ ∆⋅∆= λ . (6) Здесь Тcw – фоновая температура воды в водоносном гори- зонте (температура холодной воды); Т0, Т1 – соответственно температура воды и пород в блоке сетки размерами ΔхΔуΔz в начале и в конце периода времени осреднения τ ; q – мощность теплового потока от реакционного канала; V – скорость фильтрации; λ – теплопроводность водоносного пласта; www Vc ,,ρ – соответственно плотность, теплоемкость и объем воды в блоке; ппп Vс ,,ρ – плотность, теплоемкость и объем пород в блоке. Подставляя выражения (2)–(6) в (1), получим 2/2/ )2/()( 00 1 DAB ВТTTDAQ T cw ++ +−⋅+− = τ . Прибавим и вычтем из числителя ( ) 2 0TDA+ , тогда 2/2/ )()( 0 01 DAB TTDAQ ТT cw ++ −⋅+− += τ , (7) если рассчитывается временной ряд температур то 2/2/ )()( 1 1 DAB TTDAQ ТT cwii ii ++ −⋅+− += − − τ , (8) где Ti – температура в блоке сетки в течение i-го периода осреднения. Изменение теплоемкости С и теплопроводности λ водонос- ного пласта в зависимости от температуры Т выражается следу- ющими уравнениями [6] )20(1,020 −+= TnCCп (9) Наукові праці УкрНДМІ НАН України, № 14, 2014 Transactions of UkrNDMI NAN Ukraine, № 14, 2014 52 )1/( 2020 Тαλλλ += (10) где С20, 20λ – соответственно теплоемкость и теплопровод- ность водоносных пород в нормальных условиях (при Т = 20 °С); n,α – постоянная (зависящая от типа пород) и поправочный ко- эффициент соответственно. Изменение коэффициента фильтрации в зависимости от температуры воды учитывается изменением вязкости воды μ (Па·с) )000221,00337,01(/000183,0 2ТТ ++=µ . (11) В результате последовательных расчетов по формуле (8) можно определить переменную во времени температуру, которая задается в блоке конечно-разностной сетки при численном моде- лировании. В остальных блоках над водоупором вне реакционно- го канала задаются нулевые потоки тепла и воды. Коэффициенты граничных условий были определены по формуле (8) в программе Mathcad для горно-геологических усло- вий участка «Ольхово нижнее», расположенного в пределах Чи- стяково-Снежнянского угленосного района Донецкого бассейна. Каменоугольные отложения данной территории (рис. 2) пред- ставлены смоляниновской свитой (С2 3), в разрезе которой между отрабатываемыми угольными пластами h8 и h11 развиты несколь- ко нерабочих угольных пластов (h9, h10 н, h10 1) с углами падения от 10 до 30° [7]. Данные пласты имеют простое строение и повсе- местное распространение в пределах выделенного участка. Угли относятся к антрацитам (Wа = 2,5 %, Ас = 17 %, Sоб с = 1,8 %). Для свиты С2 3 характерно развитие обводненной толщи «бабаков- ских» песчаников h10Sh11 мощностью 50–60 м, залегающих выше угольного пласта h10 1. В пределах этой толщи подземные воды имеют гидрокарбонатно-кальциевый состав с минерализацией до 3 г/л и не используются для питьевого водоснабжения. Учитывая требования, предъявляемые к угольным пластам, отрабатывае- мым способом ПСУ [3] и фактические данные о строение участ- ка, можно считать, что на данной территории оптимальным будет разработка угольного пласта h10 1 с использованием в качестве теплоносителя воды песчаников h10Sh11. Наукові праці УкрНДМІ НАН України, № 14, 2014 Transactions of UkrNDMI NAN Ukraine, № 14, 2014 53 Рис. 2. Литолого-стратиграфический разрез Смоляниновс- кой свиты: 1, 2 – соответственно рабочие и нерабо- чие угольные пласты; 3 – водоносный горизонт, нагреваемый при сжигании угля На рисунке 3 приведены результаты расчета температуры блока обводненных песчаников, залегающих над сжигаемым угольным пластом при различных скоростях фильтрации подзем- h11 h10 1 h10 н h9 h8 h7 1 h7 h6 2 h6 1 h6 h5 C2 3 H6 1 H6 H5 1 H5 0 H4 1 3 2 Наукові праці УкрНДМІ НАН України, № 14, 2014 Transactions of UkrNDMI NAN Ukraine, № 14, 2014 54 ных вод, соответствующих периодам отбора, простоя и закачки воды в условиях участка «Ольхово нижнее»: Тcw =12 °C; Δх = Δу = 25 м; Δz = 1 м; пρ = 1900 кг/м3; λ20 = 1,07 Вт/м·°C; τ = 1 сут; n = 2,33;α = 10-3 м/Вт; С20 = 900 Дж/кг·°C; Vw = 125 м3; Сw = 4182 Дж/кг·°C; Vп = 500 м3; ρw= 1000 кг/м3. Величина теплового потока q, поступающего из реакцион- ного канала в водоносный горизонт, рассчитывалась по методике [4] при мощности водоупора 5 м. Значение скоростей фильтрации подземных вод V было определено в результате решения гидро- динамической задачи методом численного моделирования в про- грамме ModFlow 2009. Рис. 3. Изменение температуры пород водоносного горизон- та непосредственно над кровлей сжигаемого уголь- ного пласта с углом падения 10°: 1 – V=0,92 м/сут (природный режим фильтрации); 2 – V=1,13 м/сут (закачка воды в пласт); 3 – V=1,45 м/сут (отбор под- земных вод) Моделирование теплопереноса было выполнено для «Баба- ковской» водоносной толщи при закачке и отборе теплоносите- лей в течение одного года. Для решения поставленной задачи ис- пользовалась разработанная и протестированная в [5] численная модель с нестационарными источниками воды и тепла, позволя- ющая описывать переходные режимы. 1 t, сут Т, °C 2 3 Наукові праці УкрНДМІ НАН України, № 14, 2014 Transactions of UkrNDMI NAN Ukraine, № 14, 2014 55 Предложенная модель позволяет определять количество от- бираемого тепла в зависимости от скорости фильтрации воды в пласте, его наклона, мощности, расположения скважины относи- тельно реакционного канала, нестационарного характера тепло- выделения и т.д. Дебит скважины, пробуренной до подошвы во- доносного горизонта, принимался в соответствии с годовой ди- намикой потребления тепловой энергии жилищно-коммунальным сектором: 0–30 сут (сжигание угольного пласта); 30–180 сут (от- бор теплоносителя в холодное время года с дебитом 200 м3/сут); 180–245 сут (простой); 245–335 сут (закачка воды в пласт летом с расходом 333 м3/сут); 335–365 сут (подготовительные работы). На рисунках 4, 5 приведены размеры и формы гидрогеотермиче- ских полей, возникающих в «бабаковских» водоносных песчани- ках при сжигании угольного пласта h10 1 в различные периоды времени. Анализируя полученные результаты можно отметить тес- ную взаимосвязь между температурой водоносного горизонта, скоростью фильтрации подземных вод и углом наклона пласта. Конфигурация изотерм не совпадает с изогипсами водоносного горизонта – они располагаются асимметрично по отношению к реакционному каналу и вытягиваются по движению потока (рис. 4, а). Ореол нагретых вод постепенно растет вверх, прибли- жаясь к кровле водоносного горизонта, и одновременно растека- ется по ней. При этом по мере удаления от водоупора термиче- ские зоны в пласте снижают свою температуру и интенсивно пе- ремещаются от теплового источника с разрывом тепловых ано- малий (рис. 4, б). Расширение во времени области теплового воздействия происходит вследствие конвективного переноса нагретых вод под влиянием природного потока, а также под воздействием молеку- лярной (кондуктивной) теплопроводности. При этом ковективная составляющая теплопереноса имеет преобладающие значение в горизонтальном направлении, а кондуктивная – в вертикальном. Наукові праці УкрНДМІ НАН України, № 14, 2014 Transactions of UkrNDMI NAN Ukraine, № 14, 2014 56 Рис. 4. Распределение уровня подземных вод (изогипсы) и их температуры (оттенки серого цвета) в водонос- ном горизонте, залегающем над угольным пластом при его отработке способом ПСУ: а), б) – по окон- чании сжигания угля (t = 30 сут) соответственно на высоте 1 и 3 м от водоупора; в) при завершении от- бора нагретых вод (t = 180 сут); г) после последую- щего простоя (t = 245 сут); 1 – эксплуатационная скважина над реакционным каналом. Стрелками по- казаны направления течения подземных вод После ликвидации теплового источника (окончания сжига- ния угольного пласта, рис. 4 в, г) распространение теплового воз- действия в подземных водах обусловлено растеканием сформи- рованного ранее ореола нагретых вод и рассеиванием тепла в горном массиве. 50 100 150 200 250 300 150 200 250 300 350 400 30 30 40 40 50 50 60 50 100 150 200 250 300 150 200 250 300 350 400 30 30 40 40 50 50 60 а) г) IV IV х, м у, м х, м у, м Т,°С 50 100 150 200 250 300 150 200 250 300 350 400 30 30 40 40 50 50 55 55 в) III III у, м 50 100 150 200 250 300 150 200 250 300 350 400 30 30 40 40 50 50 55 55 б) х, м у, м Т,°С 1 1 х, м Т,°С Т,°С I I Наукові праці УкрНДМІ НАН України, № 14, 2014 Transactions of UkrNDMI NAN Ukraine, № 14, 2014 57 Рис. 5. Превышение температуры подземных вод в обвод- ненных породах над кровлей угольного пласта при его сжигании по различным профилям (I – I, II – II, III – III и IV – IV) спустя t=30 сут (а), t=100 сут (б), t=180 сут (в) и t=245 сут (г), при угле падения водо- носного горизонта: 1–10°; 2–20°; 3–30° х, м х, м х, м Т - Т cw , ° C Т - Т cw , ° C Т - Т cw , ° C Т - Т cw , ° C 1 1 1 1 2 2 2 2 3 3 3 3 а) б) в) г) I I II II III III IV IV х, м Наукові праці УкрНДМІ НАН України, № 14, 2014 Transactions of UkrNDMI NAN Ukraine, № 14, 2014 58 После завершения отработки угольного пласта способом ПСУ (рис. 5, а) максимум температурной аномалии (Т–Тcw) рас- полагается в пределах блока водоносных пород над реакционным каналом (х = 20 м). По мере увеличения угла наклона водоносной толщи он медленно перемещается по потоку подземных вод, при этом его максимальное значение уменьшается с 70 до 45 °C. Наибольшие температуры подземных вод характерны для наибо- лее пологого водоносного пласта. Разницы температур нагретых и естественных подземных вод по профилю I – I при х = 70 м в пластах с углами падения 10 и 20° соответственно равны 20 и 40 °C, а при х = 100 м разницы температур по тому же профилю в пластах с углами падения 20 и 30° равны 5 и 10 °C. Таким обра- зом, возникающие температурные аномалии не постоянны во времени и в пространстве. Продвижение теплового фронта навстречу движению потока не происходит на всех этапах экс- плуатации скважины, даже при отборе воды, что свидетельствует о преобладании конвективного теплопереноса над кондуктивным и указывает на необходимость увеличения дебита скважины, от- бирающей теплоноситель. Выводы. Разработана математическая модель для расчета теплопереноса в подземных водах, отражающая изменение их температуры в процессе сжигания угольного пласта и отбора теплой воды. Идентификация модели выполнена при гидрогеоло- гических и теплофизических параметрах, соответствующих гор- но-геологическим условиям смоляниновской свиты Чистяково- Снежнянского угленосного района, пригодной для разработки способом ПСУ. Расчеты теплового баланса в блоке водоносной толщи, зале- гающей над реакционным каналом, позволили оценить размеры и формы тепловых аномалий, формирующихся в обводненных по- родах кровли угольного пласта в зависимости от угла его падения и стадии отработки. Размеры аномалий в течение одного года не превышают 250 м в длину и 50 м в ширину. Превышение темпе- ратуры над фоновой изменяется от 20° до 70 °С. Полученные конфигурации термальных зон отражают течение воды в обвод- ненных песчаниках и могут использоваться для определения мест отбора нагретых вод в разные периоды времени. Наукові праці УкрНДМІ НАН України, № 14, 2014 Transactions of UkrNDMI NAN Ukraine, № 14, 2014 59 Результаты исследований позволяют, обосновано подойти к определению эксплуатационных запасов нагретых вод и возмож- ности их отбора при конкретизации геотехнологической схемы, которая предполагает закачку воды через систему нескольких скважин, с синхронизацией годового графика температур и пери- одов сжигания угольного пласта. СПИСОК ССЫЛОК 1. На рынке энергетического угля ожидается увеличение спроса и цен [Электронный ресурс] / И. Марьяненко // Электронный журнал энергосервисной компании «Экологические систе- мы». –– 2011. –– № 11. –– Режим доступа к журн. : http://www.esco.co.ua/journal/2011_11/art164.htm. –– Название с экрана. 2. Батурина И. М. Горное дело и окружающая среда. Геодина- мика недр : Учебное пособие для вузов / Батурина И. М., Пе- тухов И. М., Батурин А. С. –– М : Издательство Московского государственного горного университета, 2009. –– 120 с. 3. Теория и практика термохимической технологии добычи и переработки уля : Монография / [под. ред. О. В. Колоколова]. –– Днепропетровск : НГА Украины, 2000. –– 281 с. 4. Садовенко И. А. Численное исследование теплового поля во- круг подземного газогенератора / И. А. Садовенко, Д. В. Руда- ков, А. В. Инкин // Збірник наукових праць НГУ. –– 2012. –– № 39. –– С. 11––20. 5. Рудаков Д. В. Моделирование теплопереноса в водоносном горизонте при аккумуляции и отборе тепловой энергии / Д. В. Рудаков, И. А. Садовенко, А. В. Инкин // Науковий віс- ник НГУ. –– 2012. –– № 1. –– С. 40––45. 6. Гончаров С. А. Термодинамика: Учебник / Гончаров С. А. –– М : Издательство Московского государственного горного университета, 2002. –– 440 с. 7. Геология месторождений угля и горючих сланцев СССР. Т. 1. Угольные бассейны и месторождения юга Европейской части СССР / [Лагутина В. В., Левенштейн М. Л., Попов В. С. и др.]; под ред. И. А. Кузнецова. –– М. : Госгеолтехиздат, 1963. –– 210 с.

Оценка изменения термодинамических свойств подземных вод при подземном сжигании угля

Репозитарії

Схожі ресурси