Фильтрационная сушка угольного пласта перед подземной газификацией

УДК 622.27.326

ФИЛЬТРАЦИОННАЯ СУШКА УГОЛЬНОГО ПЛАСТА ПЕРЕД ПОДЗЕМНОЙ ГАЗИФИКАЦИЕЙ

Н.М. Качурин, М.В. Ганин, Г.В. Фридлендер, Е.К. Мосина

Подземная газификация бурых углей в Подмосковном угольном бассейне происходит в фильтрационном канале. Процесс газификации более эффективен при влажности угля, не превышающей 10... 12 %, поэтому необходимо реализовать процесс сушки угля перед началом подземной газификации. Рассматривая элементарный объем угля, который контактирует с фильтрационным воздушным потоком, формально можно предполагать, что - это система «капиллярно-пористое тело - вода» («КПТ-В»). Обоснована физическая модель и представлено математическое описание конвективной сушки угольного пласта в фильтрационном слое. Разработана математическая модель фильтрационной сушки угольного пласта перед подземной газификацией применительно к технологическим схемам Подмосковного угольного бассейна.

Ключевые слова: подземная газификация угля, угольный пласт, фильтрация воздуха, конвективная сушка, влагосодержание, технологическая схема, математическая модель.

Физическая модель и математическое описание конвективной сушки угольного пласта в фильтрационном слое. Подземная газификация бурых углей в Подмосковном угольном бассейне происходит в фильтрационном канале [1 - 3]. Процесс газификации более эффективен при влажности угля, не превышающей 10... 12 %, а начальная влажность угля, как правило, составляет 25...30 %, поэтому необходимо реализовать процесс сушки угля перед началом подземной газификации. Наиболее распространенные технологические схемы подземной газификации реализуют одномерный фильтрационный воздушный поток в угольном пласте от ряда дутьевых скважин к ряду продуктивных скважин [4 - 5].

Рассматривая элементарный объем угля, который контактирует с фильтрационным воздушным потоком, формально можно предполагать, что - это система «капиллярно-пористое тело - вода» («КПТ - В»). Вода, связанная с капиллярно-пористым телом в области положительных температур, находится в жидком состоянии. Капиллярно-пористая структура частично заполнена влагой в жидком состоянии, а остальная часть - это парогазовая смесь [6 - 8]. Количество влаги в том или ином состоянии меняется, поэтому необходимо учитывать количество влаги в капиллярах элементарного объема угля. Обозначим парообразную воду индексом 1, жидкую воду - индексом 2. Следовательно, объемная концентрация связанной воды (с) равна отношению массы воды М к объему угольной частицы V:

где Ы{ и с - соответственно масса и объемная концентрация связанной воды в состоянии /.

Концентрацию связанной воды можно представить в следующем

виде:

M А =~v = i9i

где т - пористость (объем пор и капилляров в единичном объеме древесины); ЬI - степень заполнения пор и капилляров водой /-го фазового состояния; р, - плотность воды в 1-м фазовом состоянии.

Соотношение (2) показывает, что объемная концентрация воды в капиллярно-пористой структуре угольной частицы может изменяться в процессе тепломассопереноса. Поэтому целесообразно использовать также и удельную величину, которая будет представлять собой влагосодержание угля и. Влагосодержание угля определяется по следующей формуле:

М с

и = ^Г = _'

Му Ру

где Му и ру - масса сухой угольной частицы и плотность сухого угля соответственно.

Влагосодержание угля в рассматриваемом элементарном объеме равно сумме значений и и и2 парциального влагосодержания, обусловленного жидким и парообразным состояниями воды. Расчетная схема миграции воды и тепла в произвольном элементарном объеме угля представлена на рис. 1.

Рис. 1. Расчетная схема миграции воды и тепла в угольной частице произвольного элементарного объема V

То есть объем V произвольной формы ограничен с внешней стороны поверхностью контакта с воздухом площадью S. Массоперенос в угольной частице происходит следующим образом. Под действием тепла

277

влага перемещается к поверхности раздела, с которой потом испаряется и уносится воздушным потоком. В общем случае внутри частицы могут действовать источники появления влаги и ее стоки, приводящие к исчезновению влаги. Суммарную производительность источников и стоков обозначим I. Выделим на внешней поверхности угольной частицы малый элемент (настолько малый, что его кривизной можно пренебречь, а миграцию на поверхности выделенного элемента можно считать равномерно распределенной) площадью, равной йБ. Тогда вектор миграционного потока влаги ] через элемент йБ будет перпендикулярен этому элементу и обусловит потерю влаги в единицу времени, равную Следовательно, через всю внешнюю поверхность угольной частицы уйдет количество влаги 3, численно равное сумме потоков через все малые элементы площади, слагающие данную поверхность и в сумме численно равные Б.

Таким образом, можно записать

Очевидным является тот факт, что миграция влаги и действие источников и стоков приведут к изменению количества влаги в самой частице. Выделим внутри частицы малый объемный элемент (настолько малый, что распределение влаги внутри него можно считать равномерно распределенным) величиной йУ. Тогда изменение количества влаги в этом объемном элементе в единицу времени будет равно (-ди / д1 +1)йУ. Следовательно, во всем объеме V изменение количества влаги будет равно сумме изменений во всех объемных элементах П, слагающих данный объем.

Поэтому можно записать

»-яС-дН".

В соответствии с законом сохранения массы вещества ^П. Тогда из соотношений (4) и (5) получим

Равенство (6) преобразуется по формуле Остроградского - Гаусса к следующему виду:

Ш

У

* (-1

йУ - 0.

При этом равенство (7) справедливо только в том случае, если подынтегральная функция равна нулю, следовательно, справедливо уравнение

Шу

Уравнение (8) представляет собой уравнение неразрывности потока влаги при ее миграции в рассматриваемом элементарном объеме угольной частицы. Аналогичным образом рассмотрим процесс миграции тепла в угольной частице произвольного элементарного объема V (см. рис. 1). Вновь выделим на внешней поверхности частицы малый элемент площадью, равной dS. Тогда вектор теплового потока д, проходящего через элемент dS, будет перпендикулярен к этому элементу и обусловит потерю тепла за счет теплопроводности в единицу времени, равную qdS. Часть потерь тепла будет вызвана миграцией влаги. Этот конвективный вынос тепла через элемент dS будет равен судг](уХ) dS, где суд - удельная изобарная теплоемкость воды; Т - температура системы «КПТ - В»; V - оператор Гамильтона. Следовательно, через всю внешнюю поверхность частицы уйдет количество тепла О, численно равное сумме тепловых потоков через все малые элементы, слагающие данную поверхность. Таким образом, можно записать

0 = +

Теплообмен, обусловленный теплопроводностью и миграцией влаги, приведет к изменению энтальпии самой угольной частицы. Выделим внутри угольной частицы малый объемный элемент величиной dV. Тогда изменение энтальпии этого объемного элемента в единицу времени будет равно -(св урвудТ / Э^У, где св.у и рв.у - теплоемкость и плотность влажного угля соответственно. Следовательно, во всем объеме V изменение энтальпии будет равно сумме изменений во всех объемных элементах Н, слагающих данный объем:

н = -с,.уР,.у №§с/у.

V Э

В соответствии с законом сохранения энергии О=Н. Тогда из соотношений (9) и (10) получим

#

5 V

Равенство (11) преобразуется по формуле Остроградского - Гаусса к следующему виду:

" ЭТ"

у

V

Ш ^(4 + СУд) + СвуРв

dV = 0 .

Равенство (12) аналогично равенству (7) справедливо только в том случае, если подынтегральная функция равна нулю, следовательно, можно записать уравнение

¿IV

Уравнения (8) и (13) образуют незамкнутую систему, так как не определены в явном виде миграционный поток влаги ] и тепловой поток д. Если предположить, что миграционный процесс влаги эквивалентен диффузионному процессу, обусловленному градиентами влагосодержания и температуры, то поток ] может быть определен по первому закону Фика:

. = -

где В - коэффициент диффузии влаги в угле, обусловленный неравномерностью концентрации влаги; - коэффициент термодиффузии влаги в угле.

Рассуждая аналогичным образом относительно теплового потока д, получим, что этот вектор характеризует распространение тепла в древесине за счет механизма теплопроводности. Следовательно, тепловой поток д можно определить по закону Фурье:

q--X УТ,

где X - теплопроводность древесины.

Следовательно, в данном случае будут справедливы и уравнения неразрывности диффузионного и теплового потоков (8) и (13). Однако феноменологические законы для диффузионного и теплового потоков (14) и (15) необходимо дополнить конвективными составляющими. Таким образом, можно записать

.-

Ч--X УТ + своздРвозд\Т,

где V - вектор скорости нагретого воздуха, фильтрующегося через подготовленный участок угольного пласта; свозд - удельная теплоемкость воздуха; рвозд - плотность воздуха.

Применительно к рассматриваемым физическим условиям можно считать функцию источников равной нулю. Тогда математическое описание процесса конвективной сушки угольного пласта можно представить в следующем виде:

ди + У ди + У ди + У ди - р У 2и + р У Т,

д1 дх у ду дг дТ г ^ ^

— + х ы

У дТ+У дТ+У дТ

х ^ у г

дх ду дг

а У 2Т,

/

где Vx, Vy, Vz - проекции вектора скорости воздуха на оси координат; а -температуропроводность влажного угля; х - отношение объемной теплоемкости воздуха к объемной теплоемкости влажной древеси-

ны;а = х / св.уРв.у; X = (сво3дРво3д Ж. дРв. д).

Математическая модель фильтрационной сушки угольного пласта перед подземной газификацией. Система уравнений (18) представляет собой математическое описание физической модели сушки подготовленного участка угольного пласта фильтрующимся потоком нагрето-

го воздуха. Математическое описание, представленное в виде системы (18), необходимо адаптировать к конкретным технологическим условиям. Это позволит существенно упростить математическое описание процесса за счет физически обоснованных допущений. Рассматривая подготовленный к подземной газификации участок угольного пласта, через который фильтруется нагретый воздух, вновь формально можно предполагать, что -это система «КПТ - В». Структурная модель тепломассопереноса в этом случае может быть представлена в виде схемы, показанной на рис. 2 [9 -11]. Практическая реализация подземной газификации угля в Подмосковном бассейне может быть осуществлена по технологическим схемам, представленным на рис. 3 - 4 [12 - 13].

Тепло

Рис. 2. Структурная схема тепломассопереноса при фильтрационной сушке подготовленного участка угольного пласта перед подземной газификацией

Технологическая схема реализации способа подземной газификации тонких и средней мощности пластов бурого угля содержит вертикальные продуктивные скважины 1, с которыми соединены горизонтальные продуктивные скважины 2, расположенные на границах отрабатываемого участка пласта на расстоянии 50.. .60 м друг от друга. Между горизонтальными продуктивными скважинами 2 по центру отрабатываемого участка пласта расположен ряд нагнетательных скважин 4 с шагом 15.20 м друг от друга. Между вертикальными скважинами 1 на одной оси 6, перпендикулярной линии простирания угольного пласта 7, расположена первая нагнетательная скважина 5. Скважины 1 связаны с дымососом 11.

Рис. 3. Технологическая схема подземной газификации тонких и средней мощности пластов бурого угля для получения электроэнергии

Аналогичным образом можно рассматривать технологическую схему для реализации способа подземной газификации бурого угля и комплексного использования угольного месторождения [13], представленную на рис. 4.

Схема реализации способа комплексного освоения месторождений бурого угля, представленная на рис. 4, содержит ряд дренажных скважин 1, ряд продуктивных скважин 2 газификации угля, ряд нагнетательных скважин 3 для воздуха, ряд нагнетательных скважин 4 для растворителя золошлаковых остатков угля, ряд скважин 5 для откачки продуктивного раствора на поверхность, ряд скважин 6 для заполнения выработанного пространства блока закладочным материалом. Ряды скважин в блоке расположены друг от друга на расстоянии 20.25 м. В каждом ряду располагают 10.12 вертикальных скважин на расстоянии 15.20 м друг от друга. Ряд дренажных скважин 1 подключен к водоводу 7, который соединен с узлом водоподготовки 8. Ряд продуктивных скважин 2 газификации угля подключают к газопроводу 9, соединенному с дымососом 10, который подключен к узлу очистки энергетического газа 11, установленному перед локальной газовой электростанцией 12. Ряд скважин 5 для откачки продуктивного раствора на поверхность подключен к трубопроводу 13, направленному к химико-технологическому узлу 14. Узел аккумулирования 15 диоксида углерода соединен с узлом очистки энергетического газа 11 и локальной газовой электростанцией 12. Узел водоподготовки 8, узел очистки

282

энергетического газа 11, химико-технологический узел 14 связаны с входами узла неутилизируемых отходов 16. Узел неутилизируемых отходов 16 и узел аккумулирования диоксида углерода 15 соединены с закладочным комплексом 17.

18

шякшкшяшш

К потребителю -

К потребителю -

К потребителю

В закладочный массив выработанного пространства

Рис. 4. Технологическая схема подземной газификации бурого угля и комплексного использования угольного месторождения

Для технологических схем, показанных на рис. 3 - 4, можно фильтрационный поток считать одномерным. Рассматривая одномерный процесс, можно упростить систему уравнений

ды ТТды ^ д u _ д Т — + V— = D —т + DT —-dt дх дх дх

T T — + а — дt дх

д2Т дх2:

a

где V - средняя скорость фильтрационного движения нагретого воздуха вдоль слоя дробленой древесины; а - параметр конвективного тепломассообмена; V=Vx; а = %Vx.

Решение второго уравнения системы (19) для краевых условий Т(х,0) = Т = const, Т(0,t) = Т = const, lim Т(х,t) ф да получено с использо-

х^да

ванием преобразования Лапласа и имеет вид

7

8

9

>

T(х,t) = T + 0,5 Tc

erfc

С I—Л

vr= - 0,5 J1

2v at M a

+

+ exp

v a у

erfc

WT= - 0,5 J1

2 Vat \ a

Первую и вторую производные температурного поля (20) можно представить следующим образом:

— = -0,5 Tc

дх

1

V^at

exp

г I—л

- 0,5 J -

2 V at V a

+

+exp

v a у

exp

Л i—V

VT= + 0,5jt

2V at V a

у

+ — exp

a

v a у

x erfc

'-f- + 0,5у I? 2 V at V a

diT = 0,5 Tc { 1 дх2 ' c 1 2^

х

У

I 3.3 ГТ la t -sja t

exp

i I—\2

ТГ- - 0'5у#

2 V at V a

+

+ ^ exp a

/ухЛ

erfc

v a у

+ 0,5 Jt

2Vat V a

1

yjnat

exp

Гухл

v a у

x

x

3У х

2a 2at

exp

+ 0,5Jt

2V at V a

С учетом соотношения (22) первое уравнение системы (19) можно представить как

ди тгди ^ д2ы ri ч — + V— = D -г- + /(х,t), дt дх дх2 v ;

где / ( х, t ) = D

T äi2

Краевые условия для уравнения (23) имеют следующий вид: и (х, 0) = u = const, и (0, t) = ис (t), lim u(x,t) ф да,

где ис (^) - влагосодержание передней границы слоя угля, через который

нагретый воздух поступает в осушаемый участок угольного пласта.

Динамика влагосодержания передней границы слоя угля моделируется следующим дифференциальным уравнением:

^ = К

где К - константа скорости сушки вещества угля; ки - константа равновесия содержания влаги в воздухе и в угле; d - влагосодержание воздуха.

Полагая, что влагосодержание древесины в начальный момент времени было равно ц0, получим решение уравнения (25) в виде

ис

Конвективная составляющая в уравнении системы (23) существенно больше диффузионного потока, поэтому диффузионным членом данного уравнения можно пренебречь. В этом случае общее решение системы уравнений (23) имеет вид

и(х, t) = [ис (t - х / V) + / (х, t - х / V)]ст0(г - х /V )-

и

[1 -а0 (t - х / V)] ,

0 при ^ / V,

где а0(г - х /V) - единичная функция Хевисайда; а0(г - х /V ) = <!

[1 при t > х / V.

Устойчивый технологический режим наступает за время, в несколько раз превышающее величину Ьу/У, где Ьу - длина подготовленного к газификации участка угольного пласта, равная расстоянию между рядами дутьевых и продуктивных скважин. Поэтому окончательно функцию (27) с учетом зависимости (22) можно записать следующим образом:

Г 1

{1 х У

о~Г ~ГГ7—ГГ 2л/п{у1 а х3 Уа х

ш Зх у

х ехр

Г I—л2

ТП -

2У ах V а

У

+ ехр а2

Гухл

V а у

ег&

тг=+0,5Ур

2 ах V а

1 (ух^

па х

ехр

3 у

х

V а

2а 2ах

ехр

С I—л2

+ 0,5^Р

2у/ ах \ а

где х = 1 - х/У.

Среднее значение влагосодержания угля на участке угольного пла-стаи определим по формуле

Цу

и

ср

1 7

(г) = — I и (х, г) dх.

Ц 0

^ (t ) = Кd - и0 - ^

ср

Цу

ехр

-к I г -

V

- ехр

> +

+0,5 оттсР

где Г (Ц, г) - функция, учитывающая динамику градиента температуры в

участке угольного пласта, подготовленном к газификации.

Таким образом, разработана замкнутая математическая модель, позволяющая провести вычислительные эксперименты по динамике энер-

гоемкости сушки угля на участке угольного пласта, подготовленного к подземной газификации, при различных технологических режимах. Для автоматизации вычислений разработан алгоритм моделирования основных технологических параметров сушки угольного пласта перед подземной газификацией.

Список литературы

1. Грязев М.В., Качурин Н.М., Захаров Е.И. Горнодобывающая отрасль в экономике Тульской области. Состояние и перспективы // Известия Тульского государственного университета. Науки о Земле. 2015. Вып. 2. С. 57-66.

2. Методические принципы и системный подход к обращению с отходами производства и потребления на территориях угледобывающих регионов / Н.М. Качурин, В.В. Факторович, Е.К. Мосина, Л.Л. Рыбак // Сб. науч. тр. междунар. науч.-техн. конф. «Проблемы безопасности и эффективности освоения георесурсов в современных условиях, посвященной 25-летию Горного института Уро РАН и 75-летию основателя и первого директора института, члена-корреспондента РАН Аркадия Евгеньевича Красноштейна» / Горный институт УрО РАН . 2014. С. 123-127.

3. Экологически безопасная геотехнология комплексного освоения месторождений бурого угля / Н.М. Качурин, В.И. Ефимов, В.В. Факторович, Е.К. Мосина // Безопасность труда в промышленности. 2014. № 10. С. 65-70.

4. Перспективы экологически безопасного использования отходов производства на территориях горнодобывающих регионов / Н.М. Качурин, В.И. Ефимов, Е.К. Мосина, В.В. Факторович // Безопасность труда в промышленности. 2014. № 9. С. 81-84.

5. Перспективы восстановления и комплексного развития Подмосковного буроугольного бассейна / Н.М. Качурин, С.А. Воробьев, П.В. Васильев, С.М. Богданов // Горный журнал. 2016. №2. С. 30-35.

6. Лыков А.В. Теория тепло- и массопереноса. М.; Л.: Госэнергоиз-дат, 1968. 535 с.

7. Лыков А.В. Теория теплопроводности. М.: Энергия, 1952. 392 с.

8. Лыков А.В. Теория сушки. М.: Энергия, 1968. 471 с.

9. Martin Grebner, Bernd Meyer. Performance and exergy analysis of the current developments in coal gasification technology // Fuel. No 116. 2014. P. 910-920.

10. Ganesh R. Kale *, Bhaskar D. Kulkarni, Ranjit N. Chavan. Combined gasification of lignite coal: Thermodynamic and application study // Journal of the Taiwan Institute of Chemical Engineers. No 45. 2014. P. 163-173.

11. Sankar Bhattacharya*, Kazi Bayzid Kabir, Klaus Hein. Dimethyl ether synthesis from Victorian brown coal through gasification e Current status,

and research and development needs // Progress in Energy and Combustion Science. 2013. No 39. P. 577-605.

12. Способ подземной газификации тонких и средней мощности пластов бурого угля: пат. РФ №2522785. Опубл. 20.07.2014. Бюл. №20.

13. Способ комплексного освоения месторождения бурого угля: пат. РФ №2526953. Опубл. 27.08.2014. Бюл. №24.

Качурин Николай Михайлович, д-р техн. наук, проф., зав. кафедрой, ecology_tsu_ tula@, mail.ru, Россия, Тула, Тульский государственный университет,

Ганин Михаил Павлович, инженер, ecology_tsu_ tula@, mail.ru, Россия, Тула, ПАО «Тулачермет»,

Фридлендер Григорий Владимирович, асп., ecology_tsu_ tula@, mail.ru, Россия, Тула, Тульский государственный университет,

Мосина Екатерина Константиновна, асп., ecology_tsu_ tula@, mail.ru, Россия, Тула, Тульский государственный университет

FILTRATION DRYING COAL SEAM BEFORE UNDERGROUND GASIFICATION N.M. Kachurin, M. V. Ganin, G. V. Fridlender, E.K. Mosina

Underground gasification of brown coals in Moscow Coal Basin is realized in filtration channel. Gasification process is more efficiency by coal humidity not exceeding 10 - 12 %. That is why the process of drying must be realized before starting underground gasification. If we consider elementary volume which contacting with filtration air flow we can assume that it is the system of "capillary-porous material - water" (CPM - W). The physical model and mathematical description of convective drying coal seam in filtration channel are substantiated. The mathematical model of filtration drying coal seam before starting underground gasification for technological schemes of the Moscow Cola Basin was created.

Key words: cola underground gasification, coal seam, air filtration, convective drying, humidity, technological scheme, mathematical model.

Kachurin Nikolai Michailovich, doctor of technical sciences, professor, head of a chaicr, ecology_tsu_tula@ mail.ru, Russia, Tula, Tula State University,

Ganin Michail Pavlovich, engineer, ecology_tsu_ tulaa mail.ru, Russia, Tula, OPC "Tulachermet",

Fridlender Grigoryi Vladimirovich, postgraduate, ecology tsu tulaa, mail.ru, Russia, Tula, Tula State University,

Mosina Ekaterina Konstantinovna, postgraduate, ecology tsu tulaa, mail.ru, Russia, Tula, Tula State University

Reference

1. Grjazev M.V., Kachurin N.M., Zaharov E.I. Gornodobyvajushhaja otrasl' v jeko-nomike Tul'skoj oblasti. Sostojanie i perspektivy. Izvestija Tul'skogo gosudar-stvennogo uni-versiteta. Nauki o Zemle. 2015. Vyp. 2. S. 57-66.

2. Metodicheskie principy i sistemnyj podhod k obrashheniju s othodami pro-izvodstva i potreblenija na territorijah ugledobyvajushhih regionov/ N.M. Kachurin, V.V. Faktorovich, E.K. Mosina, L.L. Rybak// Sb. nauch. tr. mezhdunar. nauch.-tehn. konf. «Prob-lemy bezopasnosti i jeffektivnosti osvoenija georesursov v sovremennyh uslo-vijah. Posvjashhennoj 25-letiju Gornogo instituta Uro RAN i 75-letiju osnovatelja i pervogo direktora instituta chlena-korrespondenta RAN Arkadija Evgen'evicha Kras-noshtejna»/ Gornyj institut UrO RAN . 2014. S. 123-127.

3. Jekologicheski bezopasnaja geotehnologija kompleksnogo osvoenija mestorozh-denij burogo uglja/ N.M. Kachurin, V.I. Efimov, V.V. Faktorovich, E.K. Mosina// Bez-opasnost' truda v promyshlennosti, 2014. № 10. S. 65-70.

4. Perspektivy jekologicheski bezopasnogo ispol'zovanija othodov proizvod-stva na territorijah gornodobyvajushhih regionov/ N.M. Kachurin, V.I. Efimov, E.K. Mosina, V.V. Faktorovich // Bezopasnost' truda v promyshlennosti, 2014. № 9. S. 81-84.

5. Perspektivy vosstanovlenija i kompleksnogo razvitija Podmoskovnogo bu-rougol'nogo bassejna / N.M. Kachurin, S.A. Vorob'ev, P.V. Vasil'ev, S.M. Bogdanov // Gornyj zhurnal, 2016. №2. S. 30-35.

6. Lykov A.V. Teorija teplo- i massoperenosa/ M.-L. Gosjenergoizdat. 1968. 535s.

7. Lykov A.V. Teorija teploprovodnosti. M.: Jenergija. 1952. 392 s.

8. Lykov A.V. Teorija sushki. M.: Jenergija. 1968. 471 s.

9. Martin Grebner, Bernd Meyer. Performance and exergy analysis of the current developments in coal gasification technology // Fuel. No 116. 2014. P. 910-920.

10. Ganesh R. Kale *, Bhaskar D. Kulkarni, Ranjit N. Chavan. Combined gasification of lignite coal: Thermodynamic and application study // Journal of the Taiwan Institute of Chemical Engineers. No 45. 2014. P. 163-173.

11. Sankar Bhattacharya*, Kazi Bayzid Kabir, Klaus Hein. Dimethyl ether synthesis from Victorian brown coal through gasification e Current status, and research and development needs // Progress in Energy and Combustion Science. No 39. 2013. P. 577-605.

12. Sposob podzemnoj gazifikacii tonkih i srednej moshhnosti plastov buro-go uglja: pat. RF №2522785. Opubl. 20.07.2014. Bjul. №20.

13. Sposob kompleksnogo osvoenija mestorozhdenija burogo uglja: pat. RF №2526953. Opubl. 27.08.2014. Bjul. №24.