generacja rozproszona
Optymalizacja doboru mocy elektrociepłowni biomasowej z modułem ORC w projekcie modernizacji ciepłowni komunalnej z kotłami węglowymi typu WR
Wprowadzenie
Stale zaostrzane przepisy Unii Europejskiej w zakresie redukcji emisji gazów cieplarnianych i innych zanieczyszczeń czy wyczerpywanie się zasobów naturalnych węgla to główne przyczyny transformacji energetycznej która obecnie dzieje się na naszych oczach. Sektorem, w którym obecnie obserwuje się znaczący przyrost działań modernizacyjnych jest ciepłownictwo systemowe, które jak powszechnie wiadomo oparte jest głównie na węglu. W efekcie działań związanych z redukcją emisji gazów będących produktami procesu spalania, poprawą efektywności energetycznej oraz obniżaniem kosztów związanych w emisją CO2 wiele przedsiębiorstw ciepłowniczych podejmuje działania w zakresie modernizacji majątku wytwórczego. W większości przypadków wiąże się to ze zmianą konfiguracji układu technologicznego, wdrożeniem nowych technologii w tym kogeneracji i trójgeneracji. Obecnie w Polsce obserwuje się dwa główne kierunki takich działań. W ramach pierwszego ciepłownie węglowe przekształcane są do postaci elektrociepłowni przez wdrożenie modułu kogeneracji gazowej, zwykle opartej na silniku tłokowym. Zwiększa to efektywność i elastyczność pracy układu a przede wszystkim pozwala wytworzyć własną energię elektryczną, spełniając przy tym wszystkie normy emisyjne. Kierunkiem drugim jest budowa układu technologicznego opalanego biomasą, bądź w postaci kotła na biomasę bądź też w postaci bloku elektrociepłowniczego. Z coraz większym zainteresowaniem spotykają się również rozwiązania wykorzystujące promieniowanie słoneczne, pompy ciepła i zasobniki ciepła.
W zakresie skojarzonego wytwarzania ciepła i energii elektrycznej z biomasy dostępnych jest obecnie kilka rozwiązań technologicznych, a wybór właściwego dla danego systemu ciepłowniczego nie jest zadaniem łatwym. Szczególnie dotyczy to układów w małej skali w rozproszonych systemach ciepłowniczych, w których planuje się wykorzystanie lokalnie dostępnych zasobów surowca. Jednym z możliwych rozwiązań jest tu technologia organicznego obiegu Rankine’a (ORC), którą omówiono w [1]. W praktyce moduły ORC zintegrowane z kotłami na biomasę stanowią w Europie najbardziej rozpowszechnioną technologię w małych elektrociepłowniach komunalnych. W pracach [2, 3] Tartière i współpracownicy stwierdzili, że łączna moc zainstalowana elektrowni ORC w 2016 r. przekroczyła 2,7 GW w ponad 705 projektach, w których rozmieszczono 1754 jednostki ORC. Udział elektrowni opalanych biomasą wynosi 11%. W 332 biomasowych blokach ORC zainstalowano 301 MW mocy elektrycznej. Zainstalowana moc elektryczna pojedynczego układu mieści się w przedziale od 0,2 do 13,0 MW. Większość systemów ma jednak moc poniżej 1,5 MW. W Polsce zrealizowano 11 projektów tego typu, przy czym wszystkie moduły ORC pochodzą od jednego wytwórcy, włoskiej firmy Turboden s.r.l. Przykładem jest tu projekt zrealizowany w latach 2012 – 2013 w Miejskim Przedsiębiorstwie Gospodarki Komunalnej w Krośnie, gdzie do istniejącej kotłowni węglowej dobudowano blok kogeneracyjny z kotłem biomasowym firmy VAS oraz modułem ORC Turboden 14 CHP SPLIT. Inwestycja przyczyniła się do redukcji zużycia węgla o około 15 tys. Mg/rok, wytwarzając przy tym ponad 7 tys MWh energii elektrycznej powstałej ze spalania lokalnych zasobów biomasy [4]. Jednakże aspekty ekonomiczne w tego typu rozwiązaniach silnie uzależnione są od systemu wsparcia OŹE, systemu handlu uprawnieniami do emisji CO2 (EUA) czy rynkowych cen biomasy i cen sprzedawanej energii elektrycznej na rynku bilansującym [5].
W niniejszym opracowaniu przedstawiono efekty zabudowy biomasowej jednostki kogeneracyjnej w układ kotłów węglowych, bazując na doświadczeniach eksploatacyjnych ex post rzeczywistego układu w zakładzie MPGK Krosno. Przedstawiono również jak modernizacja wpływa na osiągane wskaźniki energetyczne, ekonomiczne oraz ekologiczne. Głównym celem opracowania jest optymalizacja doboru wielkości elektrociepłowni biomasowej z modułem ORC dla stanu referencyjnego ciepłowni MPKG Krosno (jeszcze przed modernizacją). Dostęp do zarchiwizowanych danych pomiarowych systemu SCADA Elektrociepłowni Krosno pozwolił na zidentyfikowanie oraz opracowanie kluczowych charakterystyk pracy modułu kogeneracji biomasowej w rzeczywistych warunkach. Ponadto, wykorzystano dane parametrów technicznych podobnych jednostek kogeneracyjnych w różnych zakresach mocy, na podstawie dostępnych danych od producentów podobnych modułów ORC. W analizie założono obecne wartości cen paliw, energii elektrycznej oraz uprawnień do handlu emisją CO2 jak i również uwzględniono trend zmian ich wartości w kolejnych latach. Zebrane dane eksploatacyjne umożliwiły opracowanie modelu matematycznego na potrzeby zadania optymalizacji doboru mocy elektrociepłowni biomasowej ORC z uwagi na kryterium finansowe. Przeprowadzone obliczenia pozwoliły na weryfikację czy powstały układ w MPGK Krosno po zmianach cen, w tym świadectw pochodzenia energii elektrycznej, jakie zaszły w latach 2012 – 1019 znajduje się obecnie w zakresie parametrów optymalnych.
Opis analizowanego układu referencyjnego
Układ technologiczny elektrociepłowni w Krośnie obejmuje zintegrowane podsystemy składające się z zespołu kotłów wodnych rusztowych typu WR oraz kogeneracji opalanej biomasą. Technologia kogeneracyjna składa się z komory spalania biomasy wraz z wymiennikami ciepła oleju termalnego firmy VAS oraz modułu ORC firmy Turboden o mocy elektrycznej znamionowej 1317 kW oraz mocy cieplnej 5350 kW. Sumaryczna moc cieplna zainstalowana elektrociepłowni wynosi 41,5 MW, w tym 34,8 MW stanowi moc zainstalowana w czterech kotłach węglowych. Uproszczony schemat zintegrowanego układu kogeneracyjnego z kotłami węglowymi przedstawia na rys. 1, natomiast parametry nominalne układu przedstawiono w tab. 1.
Moc modułu ORC została dobrana tak, aby w okresie letnim blok elektrociepłowniczy pracował z około 50% swojej mocy znamionowej, będąc przy tym jedynym źródłem ciepła w układzie. Temperatura wody w układzie kogeneracyjnym nie powinna przekraczać 95°C, dlatego w momencie rozpoczęcia sezonu grzewczego, elektrociepłownia pracuje szeregowo z kotłami węglowymi, które dogrzewają wodę do wartości wymaganej zgodnie z charakterystyką sieci ciepłowniczej. Szczytowe obciążenie grzewcze systemu wynosi około 30 MW, a maksymalna temperatura wody przesyłanej do sieci wynosi około 125 °C. Uporządkowany wykresy zapotrzebowania na ciepło oraz temperatury sieci ciepłowniczej przestawiono na rys. 3. Dane te stanowią podstawowy zbiór danych wejściowych zadaniu optymalizacji doboru mocy źródła ciepła
Optymalizacja doboru mocy kogeneracji biomasowej
Główne zadanie optymalizacji w projekcie nadbudowy węglowej ciepłowni komunalnej blokiem kogeneracyjnym polega na doborze optymalnej mocy bloku kogeneracyjnego ORC. Obecne doświadczenia eksploatacyjne z obiektu w Krośnie dowodzą, że taki system kogeneracyjny może efektywnie pracować przez długi czas nawet przy zmiennych warunkach obciążenia sieci i właściwości biomasy. Średnia roczna dyspozycyjność układu w latach 2013 – 2018 przekroczyła 90%. Przy projektowaniu układu należy jednak wziąć pod uwagę fakt, iż projekt inwestycyjny jest silnie uzależniony od sytuacji rynkowej w zakresie mechanizmów wsparcia. Każda jednostka ORC, która zostanie zainstalowana w zakładzie, w którym moc zainstalowana wyrażona przez energię chemiczną w paliwie jest większa niż 20 MW, podlega europejskiemu systemowi handlu emisjami CO2 (EU ETS). W związku z tym uwarunkowania systemu handlu uprawnieniami do emisji stanowią istotne ograniczenia w zadaniu optymalizacji a obrót uprawnieniami emisyjnymi stanowi jednen z
instrumentów wsparcia finansowego dla projektu inwestycyjnego. Ponadto, podstawowym mechanizmem wsparcia funkcjonowania bloku biomasowego są świadectwa pochodzenia energii elektrycznej ze źródeł odnawialnych (zielone certyfikaty). Kolejnym ważnym aspektem jest ustalenie udziału energii elektrycznej, która będzie spożytkowana na potrzeby własne zakładu oraz jaka ilość energii będzie sprzedawana na rynku bilansującym. W ten sposób właśnie zagospodarowywana jest energia wytwarzana w elektrociepłowni w Krośnie. Ważną kwestią jest również określenie potencjału lokalnego rynku biomasy i oszacowanie ceny paliwa.
Celem optymalizacji jest wybór kluczowych zmiennych decyzyjnych, by uzyskać najkorzystniejszy wskaźnik (funkcję celu), biorąc pod uwagę występujące ograniczenia. Ograniczeniami są m.in.: przebieg obciążenia cieplnego oraz parametry w sieci ciepłowniczej, minimum technologiczne poszczególnych źródeł ciepła, czy właściwości paliw. Zadaną funkcją celą jest zatem maksymalizacja wskaźnika ekonomicznego – wartości bieżącej netto (NPV):
Gdzie:
x – niezależna zmienna decyzyjna,
DCFt – zmiana przepływów pieniężnych po modernizacji,
r – stopa dyskonta,
JCHP – nakłady inwestycyjne kogeneracji.
Aby wyznaczyć różnicowy przepływ pieniężny na potrzeby analizy dyskontowej w kolejnych latach eksploatacji należy porównać stan po modernizacji ze stanem wyjściowym ciepłowni węglowej, przy niezmienionej ilości ciepła wyprowadzanego do sieci. Prowadzi to do zależności:
Gdzie:
Eel – wytworzona energia elektryczna czynna [MWh];
j – wskaźnik zużycia energii elektrycznej potrzeb własnych;
ccert – cena sprzedaży certyfikatów [zł/ MWh];
cel – średnia cena sprzedaży energii elektrycznej [zł/MWh];
DGCO2 – lokalne zmniejszenie emisji CO2 [Mg];
cUEA – cena uprawnienia emisyjnego [zł/ Mg];
DKE – zmiana kosztów emisji zanieczyszczeń [zł];
DKW – zmiana kosztów operacyjnych w układzie węglowym (wartość ujemna) [zł];
KB – koszty operacyjne kogeneracji biomasowej [zł];
DPd – zmiana podatku dochodowego [zł].
Przyjęto, że zmienną decyzyjną przy wyznaczeniu optymalnej wartości NPV jest moc nominalna generatora ORC x1 = P˙el. Zakres dopuszczalnych wartości określono tu na 200 – 5000 kW. Na podstawie dostępnych danych na rynku producentów modułów ORC opalanych biomasą podobnej konstrukcji, przy zachowaniu parametrów nominalnych wody w skraplaczu 60/80°C, opracowano zależność sprawności wytwarzania energii elektrycznej w module ORC w funkcji mocy generatora w warunkach nominalnych [6] (rys. 4a):
Całkowite nakłady inwestycyjne na biomasowe systemy kogeneracyjne z modułem ORC zostały określone na podstawie raportu IFC z 2017 roku [7]. Dane zostały zaproksymowane następującą zależnością funkcyjną (rys. 4b):
Na podstawie historycznych danych pomiarowych pozyskanych z obiektowego systemu SCADA, opracowano charakterystyki pracy kluczowych urządzeń w elektrociepłowni Krosno [5] [8] [9]. Wyznaczone funkcje odwzorowują pracę urządzeń w rzeczywistych warunkach eksploatacji. Znormalizowaną charakterystykę kotłów węglowych WR [9] przedstawia zależność:
gdzie:
hkw – chwilowa sprawność kotła węglowego;
hkw,nom – nominalna sprawność kotła węglowego (hkw,nom = 0,835);
Q˙ kw – chwilowa moc cieplna kotła węglowego [kW];
Q˙ kw,nom – nominalna moc cieplna kotła węglowego [kW].
Znormalizowana charakterystyka kotła biomasowego przedstawia się następująco [8]:
gdzie:
hkw – chwilowa sprawność kotła biomasowego;
hkw,nom – nominalna sprawność kotła biomasowego (hkw,nom = 0,854);
Q˙ ORC – chwilowa moc cieplna skraplacza ORC [kW];
Q˙ ORC,nom – nominalna moc cieplna skraplacza ORC [kW].
Charakterystyki obciążeniowe modułu ORC w Krośnie, które wyznaczono na podstawie pomiarów [8], zostały znormalizowane i zastosowane w formie bezwymiarowej. W związku z tym obliczeniową moc cieplną skraplacza ORC wyznacza się z zależności:
gdzie:
Q˙ T14,nom – nominalna moc cieplna skraplacza ORC w Elektrociepłowni Krosno, kW
Pozwala to wyznaczyć kolejne parametry modułu ORC takie jak:
–– wartość mocy elektrycznej na zaciskach generatora ORC [8]:
gdzie:
P˙T14,nom – nominalna moc elektryczna generatora ORC w Elektrociepłowni Krosno,
kW tT66 – temperatura oleju termalnego [°C];
tw2 – temperatura wylotowa ze skraplacza ORC, [°C], – sprawności elektryczna obiegu ORC brutto [8]:
gdzie:
hT14,nom – nominalna moc elektryczna generatora ORC w Elektrociepłowni Krosno [kW]; Temperatura oleju termalnego tT66 wyznaczana jest na podstawie wzoru:
Kotły węglowe oraz kogeneracja biomasowa mogą pracować tylko w zakresie dopuszczalnych obciążeń:
Chwilowe zużycie paliwa wyznaczono z następujących zależności:
–– strumień masowy węgla [kg/s]:
gdzie:
Wdw – wartość opałowa węgla, kJ/kg
–– strumień masowy biomasy [kg/s]:
gdzie:
s – wskaźnik skojarzenia definiowany jako stosunek chwilowej mocy elektrycznej do cieplnej.
Wyniki optymalizacji
Na podstawie rocznych danych z sieci ciepłowniczej (godzina po godzinie), została zasymulowana praca kogeneracji wraz z kotłami węglowymi. Chwilowy rozkład obciążeń wszystkich źródeł ciepła dobierany był tak, aby sprawność zintegrowanego układ była w każdym kroku obliczeniowym jak najwyższa. Symulacje przeprowadzane były dla całkowitego okresu eksploatacji elektrociepłowni ORC, wynoszącej 15 lat. W obliczeniach przyjęto cztery scenariusze dotyczące mechanizmów wsparcia kogeneracji biomasowej:
- Udział dotacji na projekt inwestycyjny w wysokości 30% całkowitych nakładów oraz zyski ze sprzedaży zielonych certyfikatów,
- Udział dotacji na projekt inwestycyjny w wysokości 30% całkowitych nakładów bez dodatkowych zysków ze sprzedaży zielonych certyfikatów,
- Brak dotacji na projekt inwestycyjny lecz ujęto zyski ze sprzedaży zielonych certyfikatów,
- Całkowity brak mechanizmów wsparcia.
Założono wartości cenowe, obowiązujące na początku 2019 roku. Cenę węgla przyjęto na poziomie 14,50 zł/GJ, cenę biomasy – 18,22 zł/GJ, cenę uprawnień do handlu emisjami – 103,2 zł/Mg, cenę zielonych certyfikatów – 130 zł/MWh. Cena sprzedaży energii elektrycznej założona na podstawie danych z rynku bilansującego z 2018 roku [10] (rys. 5a). Stopę dyskonta przyjęto r = 5%. W obliczeniach ekonomicznych uwzględniono również ścieżkę prognozowanych cen energii elektrycznej oraz cen uprawnień do emisji CO2 w kolejnych latach (rys 5b). Przyjęta dynamikę wzrostu cen opracowano na podstawie literatury [11][12][13].
Obliczenia zostały wykonane dla zakresu mocy generatora ORC od 200 do 5000 kW. Dobór zbyt mocy elektrycznej elektrociepłowni ORC skutkuje spadkiem rocznej dyspozycyjności nawet do 50%. Powodem tego są ograniczenia technologiczne kogeneracji i występujące sytuacje, gdy minimalne obciążenie w sieci ciepłowniczej przewyższa wymagane minimum techniczne kotła biomasowego. Na rys. 6 i 7 przedstawiono wynikowy sumaryczny czas pracy kogeneracji dla przykładowych mocy generatora oraz dyspozycyjność źródeł ciepła dla pełnego zakresu dobranej mocy generatora.
Uruchomienie bloku elektrociepłowni może wywołać zmianę sposobu eksploatacji kotłowni węglowej. Na rys. 8 przedstawiono wyniki średniorocznych sprawności kotłowni węglowej przy pracy kogeneracji. Widać wyraźnie, że nadbudowa elektrociepłowni skutkuje spadkiem efektowności kotłów węglowych względem wartością przed dokonaną modernizacją. Przy doborze modułu ORC w zakresach 1000-1400 kW, spadek sprawności kotłowni jest nieznaczny.
Po realizacji zwiększeniu ulega rzeczywista sumaryczna ilość wyemitowanego dwutlenku węgla. Warto zwrócić uwagę, że zgodnie z obowiązującymi regulacjami dobór większej mocy kogeneracji biomasowej zmniejsza udział emisji CO2 z węgla, natomiast spalana biomasa traktowana jest jako paliwo zeromisyjne. Wykres ilości emisji CO2 z podziałem na blok węglowy i biomasowy dla danej mocy kogeneracji przedstawia rys. 9.
Wyniki optymalizacji mocy kogeneracji ORC pokazują, że zakres optymalnych wartości NPV jest dosyć szeroki (rys. 10). Inaczej mówiąc, można zarekomendować kilka wartości mocy ORC, dla których funkcja celu jest zbliżona do optymalnej. Optymalny zakres mocy generatora ORC przy maksymalnym wsparciu inwestycji oscyluje w przedziale 1200 – 2000 kW. Co istotne, w przypadku układu w Krośnie dobrana w roku 2011 moc układu znajduje się obecnie w obszarze rozwiązań
optymalnych. Rekomendowane rozwiązania charakteryzują się długą dyspozycyjnością. Wielkość jednostki kogeneracyjnej powinna również dążyć do pokrycia podstawowego obciążenia cieplnego sieci ciepłowniczej w okresie letnim. Wyniki wykazały również, że projekty kogeneracji opalanych biomasą są nadal silnie uzależnione od wsparcia finansowego. Warianty bez wsparcia finansowego lub z częściowym finansowaniem są nieuzasadnione ekonomicznie lub znajdują się na granicy opłacalności. Jednym z mechanizmów wsparcia, który ma największy wpływ na opłacalność to system handlu emisjami. Do analizy założono również prognozowaną ścieżkę wartości cen energii elektrycznej i cen CO2 w systemie handlu emisjami. Na rys 9 przedstawiono uzyskane wyniki uwzględniając ceny stałe oraz prognozowane. Wyniki pokazują, że przyjęcie ścieżki cen energii elektrycznej i CO2 istotnie poprawiło wskaźniki ekonomiczne (rys. 11).
Podsumowanie
W pracy przedstawiono i omówiono kluczowe zagadnienia związane z pracą bloków elektrociepłowniczych opalanych biomasą, w których zastosowano technologię ORC. Chociaż technologia ta stanowi istotną alternatywę w zakresie wykorzystania lokalnie dostępnych zasobów biomasy, wciąż występują istotne bariery dla jej szerokiej implementacji rynkowej. Główną z nich są niskie wartości wskaźników opłacalności projektów inwestycyjnych przy obecnym poziomie cen paliw i energii. Pomimo znacznego wsparcia rozwoju tej technologii w latach minionych nie osiągnęła ona poziomu technologii konkurencyjnej, a projekty zakładające jej wykorzystanie wciąż wymagają znacznego dofinansowania. W uwarunkowaniach obecnych systemów wsparcia finansowego projektów w obszarze energetyki odnawialnej i kogeneracji istotne znacznie dla powodzenia projektu ma postać zaprojektowanego systemu technologicznego oraz sposób prowadzenia ruchu instalacji. W szczególności w studium wykonalności projektu powinna zostać uwzględniona zmienność parametrów pracy urządzeń w systemie technologicznym oraz strategia zarządzania produkcją ciepła i energii elektrycznej.
W pracy omówiono zagadnienie doboru wartości mocy kogeneracji biomasowej ORC, która jest zabudowana w istniejącej kotłowni węglowej, pracującej na potrzeby małej miejskiej sieci ciepłowniczej. Funkcją celu był wskaźnik ekonomiczny NPV. Wyniki analizy dowodzą, że przy odpowiednim doborze poszczególnych komponentów układu technologicznego oraz ich optymalnej integracji z istniejącym systemem ciepłowniczym możliwe jest uzyskanie korzystnych wskaźników opłacalności projektu. Jednakże należy mieć na uwadze, że projekt taki będzie się charakteryzować stosunkowo długim okresem zwrotu oraz podwyższonym poziomem ryzyka. Podsumowując, w
najbliższej przyszłości nie należy się spodziewać kontynuacji trendu ekspansji rozproszonych elektrociepłowni opalanych biomasą w technologii ORC, jaki obserwowano w Europie latach 2004 do 2014. Ma to głównie związek z obniżeniem intensywności wsparcia oraz utratą statusu technologii innowacyjnej. Z drugiej jednak strony, w przypadku dalszego wzrostu cen energii elektrycznej i cen uprawnień do emisji CO2, układy kogeneracyjne w technologii ORC mogą ponownie stać się istotną alternatywą technologiczną w zakresie wykorzystania lokalnie dostępnych zasobów biomasy w rozproszonych obiektach małej mocy.
L I T E R AT U R A
[1] Musiał A., Kalina J.: Wykorzystanie przemysłowej energii odpadowej do produkcji energii elektrycznej w układach ORC. INSTAL 5/2019.
[2] Tartière T.: ORC Market: A World Overview. Web project available at http://orc-worldmap. org/analysis.html
[3] Tartière T., Astolfi M.: A World Overview of the Organic Rankine Cycle Market. Proceedings of the IV International Seminar on ORC Power Systems, ORC2017, 13-15 September 2017, Milano, Italy. Energy Procedia 129 (2017), pp. 2–9.
[4] Ćwięka J.: Ciepłownictwo w poszukiwaniu najlepszych rozwiązań – Krosno postawiło na biomasę. INSTAL 5/2016.
[5] Kalina, J., Świerzewski, M., Strzałka, R., 2019, Operational experiences of municipal heating plants with biomass-fired ORC cogeneration units. Energy Conversion and Management, Vol. 181, p.544-561.
[6] Turboden srl. ORC standard units. Datasheet. www.turboden.com.
[7] International Finance Corporation (IFC) in partnership with Austrian Federal of Ministry of Finance, 2017, Converting Biomass to Energy – a Guide for Developers and Investors. Washington, D.C. June 2017. ifc.org.
[8] Kalina, J., Świerzewski M., 2019, Identification of ORC unit operation in biomass-fired cogeneration system. Renewable Energy. In Press, Accepted Manuscript, Available online 25 April 2019, doi.org/10.1016/j.renene. 2019.04.080.
[9] Świerzewski M., Kalina, J., Implementation of heat storage and network water cooler for improvement of energy and economic performance of municipal heating plant with biomass fired cogeneration module. Journal of Power Technologies 2019 vol. 99 iss. 2, s. 131-141, bibliogr. 12 poz.
[10] https://rynek-bilansujacy.cire.pl/st,17,107, me,0,0,0,0,0,ceny-energii-na-rb.html
[11] Schnell C., Roszkowski M., Instytut Jagieloński: ceny energii. Trwały wzrost cen energii zmieni dopiero modifikacja polskiego miksu energetycznego. Czy stać nas na marnowanie czasu? Warszawa 2018.
[12] Long-Term Carbon Price Forecast Report. Submitted to: Ontario Energy Board Submitted by: ICF Consulting Canada, Inc. May 31, 2017 Updated on July 19, 2017.
[13] Luckow P., Stanton E.A., Fields S., Ong W., Biewald B., Jackson S., Fisher J.: Spring 2016 National Carbon Dioxide Price Forecast. Updated March 16, 2016. Raport agencji Synapse Energy Economics Inc. www.synapse- energy.com.