РОЗРОБКА ЕНЕРГОЗБЕРіГАЮЧОї ТЕХНОЛОГії ФУНКЦіОНУВАННЯ БіОГАЗОВОї УСТАНОВКИ У СКЛАДі КОГЕНЕРАЦіЙНОї СИСТЕМИ Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

В роботi запропонована технологiя функщонування бюгазовог установки, що дозволяв встановлювати температуру теплоноыя, що грie, на входi в теплообмт-ник, вбудований в метантенк, при вимiрюваннi темпе-ратури теплоноЫя, що грie, на виходi з теплообмтника. Використання ттегрованог системи ощнки змти темпе-ратури зброджування, здобутог на основi математич-ного та логiчного моделювання у складi когенеращйног системи, забезпечув постшний вихи) бюгазу, свовчас-не видвантаження збродженого сусла та завантаження свiжого матерiалу. Ощнена ефективтсть розробленог технологИ

Ключовi слова: технологiя, бюгазова установка, температура зброджування, математичне та логiчне моделювання, когенеращйна система

□-□

В работе предложена технология функционирования биогазовой установки, которая позволяет устанавливать температуру греющего теплоносителя на входе в теплообменник, встроенный в метантенк, при измерении температуры греющего теплоносителя на выходе из теплообменника. Использование интегрированной системы оценки изменения температуры сбраживания, полученной на основе математического и логического моделирования в составе когенерационной системы обеспечивает постоянный выход биогаза, своевременную выгрузку сброженного сусла и загрузку свежего материала. Оценена эффективность разработанной технологии Ключевые слова: технология, биогазовая установка, температура сбраживания, математическое и логическое моделирование, когенерационная система -□ □-

УДК 621. 31

[dOI: 10.15587/1729-4061.2015.44252|

РОЗРОБКА ЕНЕРГО-ЗБЕР1ГАЮЧО1 ТЕХНОЛОГИ ФУНКЦ1ОНУВАННЯ Б1ОГАЗОВО1 УСТАНОВКИ У СКЛАД1 КОГЕНЕРАЦ1ЙНО1 СИСТЕМИ

£. £. Чайковська

Кандидат техшчних наук, старший науковий ствроб^ник, доцент Кафедра теоретично!', загальноТ та нетрадицтноТ енергетики Одеський нацюнальний пол^ехычний уыверситет пр. Шевченка, 1, м. Одеса, УкраТна, 65044 E-mail: eechaikovskaya@gmail.com

1. Вступ

В умовах енергозбереження когенерацшт технологи, що виробляють як електричну енерпю, так i теплоту ввд одного джерела енергп мають додатковий iмпульс для розвитку щодо використання бюгазу як альтернативного джерела енергп [1-4]. Але функщонування бю-газових установок супроводжуеться непостшним вихо-дом бюгазу iз-за складносп визначення точного термшу змши потоюв свiжоi та зброджено сировини. В умовах недостовiрного використання ощнки змши темпера-тури зброджування при вимiрюваннi, що обумовлено значною тепловою акумулюючою емнiстю сусла, витра-ти на тдтримку процесу зброджування складають до 20-30 % виробленоi енергii. Б^ьш того, цiлодобове функцiонування бiогазових установок при непостш-ному виходу бюгазу та ввдключення когенерацiйних систем в години найменшого споживання потребуе до-даткового обладнання для збертння бiогазу та збродженого матерiалу. Цим обгрунтовуеться актуальшсть даноi роботи.

2. Аналiз лкературних даних i постановка проблеми

Статичш методи удосконалення когенерацiйних технологiй, що мають у своему складi основш скла-

довi: первинний двигун, електрогенератор, систему утилiзацii теплоти, систему контролю й управлшня дозволяють, наприклад, на основi економiчного ана-лiзу встановити дощльшсть доповнення тепловоi схеми за рахунок пда^ву повiтря пiсля компре-сора, але без можливостi ощнки И ефективностi в реальних умовах функщонування когенерацiйноi системи [5]. Не зб^ виробництва та споживання електричноi енергii та теплоти також не врахований, наприклад, в робой [6] щодо запропонованоi елек-тронноi системи подачi бiогазу з цiллю зменшення шкiдливих викидiв. В роботi [7] встановлено зв'язок виробництва бюгазу зi змшуванням рiзних сиро-винних матерiалiв, але без узгодження iз спожи-ванням. Iнтенсифiкацiя теплообмiну в метантенках розглянута, наприклад, в роботах [8, 9] щодо удосконалення конструктивних параметрiв бiогазовоi установки та ощнки змши температури зброджування за часом, без урахування и змши в об'ему бю-реактора, що ускладнюе шдтримку функщонування метантенка при змт якост сировини, ii температури, т. ш. в узгодженш iз споживанням. А в робой [10] представлено методику визначення теплообмшу зовнiшньоi поверхнi метантенка з навколишшм се-редовищем щодо визначення витрати теплоноия, що грiе та ощнки юлькоси здобутого бiогазу також без узгодження iз споживанням.

© е. е.

Запропоноваш засоби тдтримки функцiонування бiогазових установок не враховують значну теплову аку-мулюючу eмнiсть сусла та базуються на вимiрюваннi температури зброджування щодо змiни витрати сусла для тдтримки процесу зброджування, що може пору-шити обов'язковий баланс потокiв свiжоi та зброджено! сировини для здобуття постшного виходу бiогазу. Вико-ристання ж змши витрати теплоносiя, що ^е, може негативно вплинути на актившсть процесу зброджування. У зв'язку iз залежнiстю температури свiжого сусла, що завантажують, вiд температури навколишнього середо-вища запропонованi засоби не дають також вiдповiдь на встановлення точного термшу вiдвантаження збродже-ного сусла та завантаження е^валентно! витрати свь жого матерiалу, що виконують через 4-6 годин на добу. Тому тдтримка функцiонування бюгазово! установки повинна вiдбуватись на основi iнтегрованоi аналогично! оцiнки змши температури зброджування на рiвнi прийняття ршень для пiдтримки головних складових процесу здобуття бюгазу: виходу бюгазу, вiдвантаження збродженого сусла та завантаження свiжого матерiалу [11]. Для своечасного прийняття ршень необхiдно з щею цiллю прогнозувати змiну температури зброджування, вимiрюючи температуру теплоносiя, що ^е, на виходi з теплообмiнника, вбудованого в метантенк, змша яко! ввдбуваеться ранiше, нiж змiна температури зброджування та встановлювати температуру теплоноая, що грiе, на входi в теплообмшник [12-14].

3. Мета та задачi дослщження

Мета роботи - розробити енергозбер^аючу техно-логiю функщонування бiогазовоi установки у складi когенерацшно! системи.

Поставлена мета може бути досягнена при вико-наннi таких задач:

- обгрунтувати необхiднiсть розробки штегровано! аналiтичноi оцiнки змiни температури зброджування на рiвнi прийняття ршень для пiдтримки функщо-нування бюгазово! установки у складi когенерацiйноi системи;

- розробити структурну схему комплексного ма-тематичного моделювання динамiки бiогазовоi установки та теплообмшника, вбудованого в метантенк, щодо здобуття еталонно! ощнки змiни температури зброджування;

- розробити структурну схему лопчного моделювання контролю працездатност бiогазовоi установки для здобуття тдсумково! функцiональноi оцiнки змi-ни температури зброджування;

- розробити структурну схему лопчного моделювання прийняття ршень на тдтримку функщону-вання бiогазовоi установки у складi когенерацiйноi системи;

- розробити структурну схему лопчного моделювання вдентифжацп стану бiогазовоi установки для тдтвердження прийнятих рiшень;

- запропонувати штегровану систему оцiнки змiни температури зброджування щодо енергозбертючо! технологи функцiонування бюгазово! установки;

- оцiнити практичну значушдсть здобутих резуль-татiв.

4. Енергозбери'аюча технологiя функцiонування бюгазово! установки у складi когенерацшно! системи

Основою для розробки енергозбер^аючо! техно-логii функщонування бюгазово! установки е здобуття як еталонно!, так i функщонально! шформа-цii щодо прийняття ршень на пiдтримку процесу зброджування. Розроблена математична модель ди-намiки бюгазово! установки для аналiтичноi оцiнки змши температури зброджування при зм^ температури теплоноия, що грiе у зв'язку iз недосто-вiрним використанням вимiрювання температури зброджування, що обумовлено значною тепловою акумулюючою емнiстю сусла. Система нелшшних диференцiальних рiвнянь включае рiвняння стану щодо оцiнки фiзичноi моделi бiогазовоi установки, рiвняння енергп теплоносiя, що грiе, для тдтримки процесу зброджування, рiвняння енергii сусла, що зброджуеться та рiвняння теплового балансу для стшки теплообмiнника, вбудованого в метантенк. Вщмшною рисою математично! моделi е рiвняння енергii сприймаючого середовища - сусла, що опи-суе змшу температури зброджування як у чаи, так i вздовж просторово! координати осi теплообмiнника, що ствпадае з напрямком потоку руху середовища. В результатi реалiзацii математично! моделi динамiки бiогазовоi установки отримана передатна функщя за каналом: «температура зброджування -температура теплоноия, що гр1е» та виконана математична обробка передатно! функцii щодо можли-востi переходу з обласи частоти до реального часу функщонування бюгазово! установки [12-14].

У зв'язку зi змiною температури сусла, що завантажують, ввд змiни температури оточуючого середовища та не можливютю змшювати витрату сусла впро-довж процесу зброджування, щоб не порушити баланс потокiв збродженого та свiжого матерiалу, необхiдно встановити таю рiвнi функцiонування бiогазовоi установки щодо змши температури теплоноая, що ^е, на входi та на виходi iз теплообмiнника, при яких би цей баланс потоюв збертвся. Так, з цiею цiллю за-пропоновано структурну схему (рис. 1) комплексного математичного моделювання динамжи бюгазово! установки та теплообмiнника, вбудованого в метантенк. Використано, наприклад, вихщш даш конструк-тивно-режимно! реалiзацii метантенка, що виробляе 352,5 м3/добу бiогазу. В межах запропоновано! циклiч-но! структури при зберiганнi постшними витрати по-токiв сировини та теплоноая, що грiе, встановлено наступш рiвнi функцiонування бiогазовоi установки щодо змши температури теплоноая, що ^е, на входi в теплообмiнник та на виходi з теплообмiнника: перший рiвень: 55 оС-44,15 оС; другий рiвень: 52,1 оС-43,06 оС; третш рiвень: 49,9 оС-42,3 оС;четвертий рiвень: 47,5 оС-41,35 оС;п'ятий рiвень: 45 оС-40,3 оС, що вщ-поввдають змiнi температури сусла на входi в бiогазову установку: 5 оС, 10 оС, 14 оС, 18 оС, 22 оС в залежност вiд температури оточуючого середовища, вщповщно (рис. 1). Запропонована структурна схема дозволяе, здобувши постiйнi часу та коефвденти математично! моделi динамiки температури зброджування, визнача-ти гранично припустиму змшу температури зброджування для встановлених рiвнiв функщонування [12-14].

Для розробки енергозберiгаючоi технологи функ-цюнування бiогазовоi установки на основi метода графа причинно-наслщкових зв'язкiв [11] представлено лопчну модель контролю працездатност бiогазовоi установки, що надае можливють отримати вiд блоку контролю СТс (рис. 2) пiдсумкову iнформацiю: (1-4).

При безперервному вимiрюваннi температури те-плоносiя, що грiе, на виходi з теплообмiнника, вбу-дованого в метантенк, що порiвнюеться з еталонним значенням рiвня функцiонування, з використанням логiчноi структури в межах циклу запропоновано тд-тримувати розряд чи заряд бiогазовоi установки, що супроводжуеться виходом бюгазу, на основi здобутоi iнформацiйноi оцiнки (1), (2), вщповщно. Здобуття ж пiдсумковоi iнформацii (3) при досягненш температури теплоноия, що грiе, на виходi з теплообмшника гранично низького рiвня - 37,04 оС, потребуе прийнят-тя ршення на вiдвантаження збродженого сусла та завантаження свiжого матерiалу щодо встановлення температури теплоноая, що ^е, на входi в теплооб-мшник, на гранично верхньому рiвнi -55 оС [12-14].

Пiдтвердження правильностi прийнятих ршень виконано на основi запропонованоi структурноi схеми iдентифiкацii стану бiогазовоi установки у складi ко-генерацiйноi системи (рис. 4), що розроблена на основi графа причинно-наслiдкових зв'язюв [11].

На основi запропонованих структурних схем (рис. 1-4), що апробоваш для тдтримки функцюну-

вання бiогазовоi установки для фiксованого штер-валу часу [12-14], розроблено штегровану систему тдтримки функцюнування бiогазовоi установки (табл. 1, рис. 5), що надае можливють на основi аналь тичноi оцiнки змiни температури зброджування при-ймати рiшення на змiну температури теплоноия, що грiе, на входi в теплообмiнник, вбудований в метантенк, при безперервному вимiрюваннi температури теплоноия, що грiе на виходi з теплообмiнника. Запро-понована система дозволяе виконувати заряд когене-рацiйноi системи впродовж термшу розряду та заряду бiогазовоi установки щодо виходу бюгазу та встанов-лювати точний термш вiдвантаження збродженого матерiалу та завантаження свiжого сусла щодо забез-печення тдтримки розряду когенерацiйноi установки.

Так, наприклад, через 2600 с тсля завантаження свiжого матерiалу температура теплоноая, що грiе, на виходi з теплообмшника, вбудованого в метантенк, зменшилась до 43,6 оС, то для тдтримки виходу бюгазу необхщно встановити температуру теплоноая, що ^е, на входi в теплообмiнник 52,1 оС. А якщо ж температура теплоноия, що грiе, на виходi з теплообмiнника зменшилась до 37,04 оС , то у зв'язку iз неможливiстю пiдтримувати температуру зброджування ввд 34 оС до 36 оС необхiдно виконати ввдвантаження збродженого сусла та завантажити свiжий матерiал, встановивши температуру теплоносiя, що ^е, на входi в теплооб-мiнник 55 оС.

Рис. 1. Структурна схема комплексного математичного моделювання динамки бюгазовоТ установки та теплообмшника, вбудованого в метантенк: dв , — дiаметр теплообмi нника, внутршнш, зовнiшнiй, м, вщповщно; Gв, Gз витрата теплоносiя, що грie, сусла, кг/с, вiдповiдно; ^ , tз, t5 — температура сусла на входi в бюгазову установку в залежностi вщ температури оточуючого середовища, К, вщповщно; ^их. — температура сусла на виходi iз бiогазовоТ установки, К; а в — коефщieнт тепловiддaчi вщ теплоносiя, що грie, до стшки теплообмiнника, вбудованого в метантенк, Вт/(м2-К); аз — коефiцieнт тепловiддaчi вщ стiнки теплообмiнникa, вбудованого в метантенк, до сусла, Вт/(м2-К); к — коефщieнт теплопередaчi, Вт/(м2^К)

Блок контролю ЕПЛИЕ1Е. СТ\ -*- (СВД, (х0 (т).-Т; (т)(СЕ (т)(0):С3(т)(0):

Блок ошнки зшни температури спнки теплообмшника,, СТ2 — (СТ2 (т)С/(т)(е(т)С-)(+)(г_ (т)(о)(о), ^ (!)(-)(-))))), гс2(т).

Блок ошнки зшни коефйцента Кв. СТз

— (СГ3 (т)(А- (т)(А'Е (т)(+)(-)))), (т).

Блок ошнки зшни температуры зброджування, СТд

(СТч(т)О{тХ^Х0)}}),^4(т},

Блок ошнки зшни дина шч них

параметр ¡в, СТз *

(СТ5 (тХ^(тХЛг(г) / дги (т) (-)(+)))), (т).

Блок контролю працездагносп,

СТс

(СТ (гхадг) 1. ш (г)< CT.pDip.S5pX V /

СГРО!

(СГс(гХА/(г)/ (г) >

грч»^ СО < !))■

3.(СГДг)(Дг(г) ,.„(*)< 0)).

л. (сгс(тхмгу < го-

Рис. 2. Стуктурна схема лопчного моделювання контролю працездатностi бюгазовоТ установки: СТ — контроль поди; Z — логiчнi вiдносини; d — динамiчнi параметри; t — температура зброджування, К; ^ — температура теплоноая, що грie, К; х — впливи; f — параметри, що дiагностуються; у — вихщж параметри; К— коефiцieнти математичного опису;

1(9с - tв.

К =:

G„

- — коефщieнт математичноТ моделi динамiки бiогазовоТ установки, де 0 — температура стшки

теплообмшника, К, вiдповiдно; G — витрата теплоноая, що грie, кг/с; 1ндекси: с — контроль прaцездaтностi; вх., вих.— вхiд i вихiд теплоноая, що грie; верх.— верхнiй рiвень функцiонувaння; в — теплоносiй, що грie; з — сусло; розр., ст. — розрахункове, стале значення параметра; 0, 1, 2 — початковий стацюнарний режим, зовжшш, внутршж параметри; 3 — коефiцieнти рiвнянь динaмiки; 4 — суттeвi параметри, що дiaгностуються; 5 — динaмiчнi параметри

Таблиця 1

1нтегрована система шдтримки функцiонувaння бiогaзовоТ установки

Змiна температури зброджування Час, т, 100 с

13 26 39 52 65

Завантаження; ^ =55 оС з.вх ^ 1 - - - -

Розряд-заряд; ^З.вих=43,6 °С - 0,8874 - - -

Прийняття ршення; ^з.вх 52,1 С - - 0,8866 - -

Розряд-заряд; 1 =42 6 °С - - - 0,8130 -

Прийняття рiшення; ^з.вх 49,9 С - - - - 0,8119

Змiна температури зброджування Час, т, 100с

78 91 104 117 130

Розряд-заряд; 1 =415 оС из.вих ^ 0,6871 - - - -

Прийняття рiшення; 1з.вх=47,5 оС - 0,6823 - - -

Розряд-заряд; ^з.вих-40 С - - 0,4872 - -

Прийняття рiшення; 1 =45 оС з.вх ^ - - - 0,4870 -

£з.вих=37,04 оС. Вiдвантаження - - - - 0

Рис. 3. Структурна схема шдтримки процесу здобуття бюгазу на рiвнi прийняття ршень, де СТ — контроль поди; t — температура зброджування, К; tз —температура

теплоносiя, що грie, К; tзе — еталонне значення температури теплоноая, що грie, К; т — час, с. 1ндекси:

с — контроль прaцездaтностi; i — число рiвнiв функцiонувaння; верх. — перший рiвень функцiонувaння; розр. ст. — розрахункове, стале значення параметра

Рис. 4. Структурна схема лопчного моделювання щентифкацм стану бюгазовоТ установки: ST — щентифшащя стану; Р — властивост бiогазовоТ установки; Z — лопчж вiдносини; d — динамiчнi параметри; t — температура зброджування, К; f — параметри, що дiагностуються; у — вихiднi параметри; К— коефщieнти

т (б0 - '

математичного опису; K =-

G„

— коефiцieнт математично'1 моделi динамiки бюгазовоТ' установки,

де 0 — температура стшки теплообмшника, К, вщповщно; G — витрата теплоноая, що rpie, кг/с; 1ндекси: s — стан; верх. — верхнш рiвень функцiонування; розр., ст. — розрахункове, стале значення параметра; 0, 2 — початковий стацюнарний режим, внутршш параметри; 3 — коефщieнти рiвнянь динамiки; 4 — суттeвi параметри,що дiагностуються;

5 — динамiчнi параметри

процесу здобуття 6iora3y: виходу бюгазу, ввдвантаження збродженого сусла та завантаження свiжого матерiалy.

2. Запропоновано структурну схему комплексного математичного моделювання динамжи бюгазовоТ установки та теплообмшника, вбудованого в метантенк, що надае можливють визначати гранично припу-стиму змшу температури зброджування для встанов-лених рiвнiв функцюнування бюгазовоТ установки.

3. Запропоновано структурну схему лопчного моделювання контролю працездатност бюгазовоТ установки для здобуття тдсумково'Т функцюнально'Т оцш-ки змши температури зброджування щодо порiвняння iз гранично припустимою.

4. Запропоновано структурну схему лопчного моделювання вибору здобутоТ тдсумково'Т функционально! шформацп щодо пiдтримки головних складових процесу здобуття бюгазу: виходу бюгазу, ввдвантаження збродженого сусла та завантаження свiжого матерiалy.

5. Запропоновано структурну схему лопчного моделювання вдентифжацп стану бюгазовоТ установки щодо тдтвердження правильност прийнятих рiшень пiсля входження в новi умови фyнкцiонyвання.

6. Запропонована штегрована система пiдтримки функцюнування бюгазовоТ установки, що на основi штегровано'Т аналiтичноï оцiнки змiни температури зброджування дозволяе встановлювати температуру теплоноая, що ^ie на входi в теплообмiнник, вбу-дований в метантенк при вимiрюваннi температури теплоноая, що грie на виходi iз теплообмiнника та виконувати своечасне вщвантаження збродженого сусла та завантаження свiжого матерiалy.

7. Запропонована технологiя функцюнування бь огазово'Т установки у складi когенерацiйноï системи

Рис. 5. 1нтегрована система пiдтримки функцiонування бюгазовоТ' установки: ^-вх — температура теплоносiя,

що грie на входi в теплообмшник, вбудований в метантенк, що встановлена на рiвнi прийняття рiшень;

t.3

температура теплонос1я, що гр1е на виход|

з теплообмiнника, вбудованого в метантенк, що вимiрюeться при функцюнуванш бюгазовоТ' установки

5. Висновки

В результат проведених дослiджень встановлено, що: 1. Когенерацшш технологи з використанням бюгазу потребують iнтегрованоï системи оцшки змiни температури зброджування щодо тдтримки головних складових

надае можлившть, використовуючи, наприклад, те-пловий насос при виробництвi 352,5 м3/добу бiогазу здобути економiю бiогазу 25,4 тис. M3/piK, що при тд-вищеннi товарностi бiогазовоi установки на 13,94 % знижуе собiвартiсть виробництва електроенергп та

теплоти в межах 20-30 %. Р1чна економiя енергп в оди-ницях умовного палива складае 19,5 т. у. п., а грошовий еквiвалент додатково виробленоi енергii - близько 100 тис. гри./рж [10-12].

Лiтература

1. Железная, Т. А. Перспективы производства электрической энергии из биомассы в Украине [Текст] / Т. А. Железная, Е. Н. Олейник, А. И. Гелетуха // Промышленная теплотехника. - 2013. - Т. 35, № 6. - С. 67-75.

2. Rade, M. C. Techno-Economic Analysis of Biogas Powered Cogeneration [Text] / C. M. Rade, Z. Kuzmanovic // Journal of Automation and Control Engineering. - 2014. - Vol. 2, Issue 1. - P. 89-93. - doi: 10.12720/joace.2.1.89-93

3. Doseva, N. Advanced exergatic analysis of cogeneration system with a biogas engine [Text] / N. Doseva // 14th SGEM GeoConference on Energy and Clean Technologies Conference Proceedings, June 19-25. - 2014. -Vol. 1. - P. 11-18. doi: 10.5593/sgem2014/b41/s17.002

4. Moedinger, F. Innovate biogas Multi-Stage Biogas Plant and Novel Analytial System [Text] / F. Moedinger, F. Ragazzi, M. Ast, P. Foladori, E. C. Rada, R. Binnig // Energy Procedia. - 2012. - Vol. 18. - P. 672-680. doi: 10.1016/j.egypro.2012.05.082

5. Тодорцев, Ю. К. Выбор экономической схеми регенерации тепла когенерационной энергетической установки [Текст] / Ю. К. Тодорцев, О. С. Тарахтий, А. Н. Бундюк // Восточно-Европейский журнал передовых технологий. -2015. - Т. 2, № 8(74). -

C. 17-22. doi: 10.15587/1729-4061.2015.40401

6. Daingade, P. S. Electronically operated fuel supply system to control air fuel ratio of biogas engine [Text] / P. S. Daingade, S. D. Yadav // 2013 International Conference on Energy Efficient Technologies for Sustainability, 2013. - P. 740-743. doi: 10.1109/iceets.2013.6533476

7. Talukder, N. Technical and economic assessment of biogas based electricity generation plant [Text] / N. Talukder, A. Talukder,

D. Barua, Anindya // 2013 International Conference on Electrical Information and Communication Technology (EICT), 2014. - P. 1-5. doi: 10.1109/eict.2014.6777854

8. Ратушняк, Г. С. 1нтенсифшащя теплообмшу та термостабшзащя бiореакторiв [Текст] / Г. С. Ратушняк, В. В. Джеджула // Вюник Вшницького пол^ехшчного шституту. - 2006. - № 2. - С. 26-31.

9. Ратушняк, Г. С. Моделювання нестацюнарних режимiв теплообмшу в бюгазових реакторах [Текст] / Г. С. Ратушняк, В. В. Джеджула, К. В. Анохша // Вюник Хмельницького нацюнального ушверситету. - 2010. - № 2. - С. 142-145.

10. Остапенко, Д. В. Тепловые процессы в метантенке при сбраживании биомассы [Текст] / Д. В. Остапенко, О. В. Чеботарева, О. В. Сербин // Восточно-Европейский журнал передовых технологий. - 2007. - Т. 6, № 5(30). - С. 18-20.

11. Чайковская, Е. Е. Оптимизация энергетических систем на уровне принятия решений [Текст] / Е. Е. Чайковская // Промышленная теплотехника. - 2013. - Т. 35, № 7. - С. 169-173.

12. Чайковська, 6. 6. Шдтримка функцюнування бюгазово'1' установки у склад1 когенерацшно'1' системи [Текст] / 6. 6. Чай-ковська, Б. I. Молодковець // Вюник Нацюнального техшчного ушверситету «ХП1». Збiрник наукових праць. Серiя: Ме-хашко-технолопчш системи та комплекси. - 2014. - № 60(1102). - С. 31-36.

13. Чайковська, 6. 6. Розробка методу тдтримки функцюнування бюгазово'1' установки у склад1 когенерацшно'1' системи [Текст] / 6. 6. Чайковська, Б. I. Молодковець // Технологический аудит и резервы производства. - 2015. - Т. 1, № 1 (21). - С. 41-46. doi: 10.15587/2312-8372.2015.37190

14. Чайковська, 6. 6. Комплексне моделювання бюгазово'1' установки у склад1 когенерацшно'1' системи [Текст] / 6. 6. Чайковська, Б. I. Молодковець // Вюник Нацюнального техшчного ушверситету «ХП1». Збiрник наукових праць. Серiя: Енергетичш та теплотехшчш процеси й устаткування. - 2015. - № 17(1126). - С. 135-143.