Log In

OBNIŻENIE PARAMETRÓW TEMPERATUROWYCH MIEJSKIEJ SIECI CIEPŁOWNICZEJ


W pracy przedstawiono efekty energetyczne, ekologiczne i ekonomiczne usprawnienia miejskiej sieci ciepłowniczej poprzez wprowadzenie nowych tabel regulacyjnych dostosowujących parametry pracy systemu ciepłowniczego do aktualnych potrzeb użytkowników, tj. dostosowania systemu do pracy w warunkach dynamicznie przebiegających procesów termomodernizacyjnych w budownictwie oraz wprowadzania indywidualnych systemów rozliczeń za zużyte ciepło. Na podstawie uzyskanych wyników stwierdzono m.in., że obniżenie temperatury zasilania z 135 do 120oC oraz powrotu z 70 do 60oC pozwala na obniżenie mocy zamówionej i rocznego zużycia nieodnawialnej energii pierwotnej pokrywających straty ciepła na przesyle nośnika ciepła o 12 proc. Zakładając, że prace związane z obniżaniem parametrów temperaturowych sieci ciepłowniczych nie wymuszają wykonania prac skutkujących dodatkową emisją to wraz z efektem energetycznym uzyskany zostaje również efekt obniżenia emisji zanieczyszczeń do atmosfery. Łączna kwota szacowanych oszczędności z tytułu obniżenia strat ciepła sieci ciepłowniczych, poprzez obniżenie temperatur wody sieciowej, może wynieść w tym przypadku 3146 PLN/MW dla standardowego sezonu grzewczego.

1. Wprowadzenie
Udział strat ciepła w sieciach ciepłowniczych w Polsce wynosi od 5 proc. do 15 proc. w okresie sezonu zimowego oraz od 20 proc. do 30 proc. w czasie sezonu letniego [1]. Jednym z parametrów determinujących straty ciepła na przesyle jest temperatura czynnika grzewczego. Wraz z jej wzrostem zwiększa się udział strat, co z kolei pociąga za sobą większe koszty wytworzenia i dostarczenia zamówionej mocy cieplnej przez odbiorców. Tak jak przedstawiono w pracy [2], można wyróżnić dwa zasadnicze kierunki prac zmierzających do obniżenia strat ciepła w sieci ciepłowniczej. Pierwszy kierunek to wzrost oporu przenikania ciepła przez ściankę przewodów ciepłowniczych. Proces ten realizowany jest w trakcie wymian wyeksploatowanych sieci ciepłowniczych wykonanych w technologii kanałowej na sieć preizolowaną, przy czym są to procesy wymuszone koniecznością bieżącej eksploatacji sieci i często obarczone wysokimi kosztami inwestycyjnymi.

Z drugiej strony straty ciepła w sieci ciepłowniczej są wprost proporcjonalne do różnicy temperatury wody sieciowej i powietrza zewnętrznego. Na rysunku 1 przedstawiono przykładowy rozkład temperatur wody sieciowej w funkcji temperatury powietrza zewnętrznego dla systemu ciepłowniczego pracującego w warunkach braku dostosowania temperatur wody sieciowej do bieżącego zapotrzebowania na ciepło. Jak wynika z danych przedstawionych na rysunku 1 ilość zdarzeń, w których temperatury wody sieciowej są wyższe od wartości założonych w tabelach regulacyjnych jest znacząca. Każdy tego typu stan pracy sieci przyczynia się do wzrostu nieuzasadnionych strat ciepła. Z drugiej strony wyniki pomiarów przedstawione na rysunku wskazują, że eksploatacja prezentowanej sieci ciepłowniczej przy wymuszonej pracy z temperaturami wody zasilającej nie przekraczających 120oC, a więc niższymi w porównaniu do temperatur zgodnych z tabelą regulacyjną, nie spowodowała żadnych problemów z dostawą ciepła do odbiorców końcowych i kolizji z wymaganymi minimalnymi spadkami temperatur w węzłach cieplnych. Kolejnym problemem dotyczącym pracy sieci ciepłowniczych jest występowanie znacznej rozbieżności pomiędzy projektowym a rzeczywistym zapotrzebowaniem na moc cieplną u odbiorców końcowych.

Na rysunku 2 przedstawiono wyniki badań przeprowadzonych na 50-ciu budynkach edukacyjnych na terenie Częstochowy. Jak wskazują badania [3, 4] oraz dane przedstawione na rysunku 2 rozbieżności te mogą wynosić nawet do 60 proc. Dlatego też, modernizacja sieci ciepłowniczych obejmować powinna opracowanie metod określania czynników kształtujących zapotrzebowanie na ciepło przez odbiorców komunalnych w nowych warunkach ich eksploatacji oraz zagadnienia dotyczące optymalizacji obliczeniowych temperatur wody sieciowej.

Rys. 1. Temperatura wody sieciowej w zależności od temperatury zewnętrznej
dla przykładowego systemu ciepłowniczego
Rys. 2. Rozbieżność pomiędzy teoretycznym zapotrzebowaniem na ciepło a rzeczywistym zużyciem ciepła dla 50 budynków edukacyjnych w Częstochowie

Wyniki tych prac pozwolą na uzyskanie najkorzystniejszego cieplnego i hydraulicznego stanu pracy centrali, sieci cieplnej oraz podłączonych do niej węzłów cieplnych, a co za tym idzie dostawę wymaganej przez użytkownika ilości ciepła przy zwiększonej efektywności systemu ciepłowniczego.

2. Zmienność temperatur powietrza zewnętrznego
Zapotrzebowanie na moc cieplną silnie uzależnione jest m.in. od warunków atmosferycznych występujących w danym sezonie grzewczym. Stąd też, w celu określenia właściwych parametrów wody sieciowej wyznacza się wykres regulacyjny, czyli zależność regulowanego parametru (temperatury wody sieciowej, strumienia masy wody sieciowej) od temperatury zewnętrznej. Ocenę sezonów grzewczych dokonać można na podstawie porównania trzech zasadniczych wielkości: długości trwania sezonu, średniej temperatury otoczenia tego okresu, względnie liczby stopniodni ogrzewania. W zależności od przyjętej wielkości sezon może być: długi lub krótki, o wysokiej lub niskiej temperaturze, mniej lub bardziej energochłonny. Stosując więc do oceny minionych sezonów metodę stopniodni w tabeli 1 zestawiono kolejne sezony grzewcze od 1963/1964 do 2009/2010 dla Częstochowy [5] oraz niezbędne dane w postaci: długości sezonu, sezonowej temperatury otoczenia oraz liczby stopniodni sezonu.

Tabela 1. Zestawienie sezonów grzewczych [5]

Wielkości te uzupełniono względną liczbą stopniodni, wyznaczoną jako iloraz aktualnej liczby stopniodni do średniej liczby stopniodni obliczonej dla wszystkich sezonów. Takie ujęcie liczby stopniodni umożliwiło wyznaczenie względnych wartości ekstremalnych. Odczytana wartość maksymalna względnej liczby stopniodni wynosi 1,16, zaś minimalna osiągnęła poziom 0,78. Na podstawie zestawionych danych liczbowych stwierdzić można, iż kryterium energochłonności jest wskaźnikiem uniwersalnym pozwalającym jednoznacznie ocenić analizowany sezon grzewczy. Uzasadnieniem tego stwierdzenia są dane zawarte w tabeli 2, gdzie przedstawiono sezony o małej (2006/2007) i dużej (1978/1979) energochłonności. Sezon 2006/2007 spełniał również kryterium wysokiej temperatury otoczenia, zaś sezon 1978/1979 kryterium długości, sezon długi. Nadmienić należy, że w analizowanych sezonach występował również sezon krótki (1991/1992) o liczbie dni 194, lecz niższej średniej temperaturze powietrza zewnętrznego.

Tabela 2. Porównanie sezonów grzewczych w latach 1963/1964 – 2009/2010 [5]

Można więc przyjąć, zachowując pewien margines dokładności, i pomijając sezony najcieplejsze, że liczba stopniodni, w porównaniu do wartości średniej, zmienia się dla analizowanych sezonów (sezony od 1963/1964 do 2009/2010) o ± 15 proc. Na rysunku 3 zamieszczono przebieg zmian liczby dni ogrzewania, sezonowej temperatury otoczenia oraz, wyrażonej w tysiącach, liczby stopniodni w kolejnych sezonach grzewczych od 1963/1964 do 2009/2010. Dane te ilustrują wzajemną zależność tych trzech parametrów w analizowanym przedziale czasu, jednak trudnym jest określenie systematycznej zależności pomiędzy nimi. Jak przedstawiono na rysunku 3, obserwowany jest trend wzrostu wartości sezonowej temperatury otoczenia wraz z kolejnymi sezonami grzewczymi. Przebieg zmian względnej wartości stopniodni w kolejnych sezonach grzewczych zamieszczono na rysunku 4, który zawiera dwa zestawy krzywych.

Rys. 3. Rozkład liczby dni ogrzewania, liczby stopniodni oraz średniej temperatury otoczenia dla sezonów grzewczych od 1963/1964 do 2009/2010 [5]

Pierwszy z nich reprezentuje względną liczbę stopniodni występującą w poszczególnych sezonach grzewczych, natomiast zestaw drugi przedstawia przebieg tego parametru, uporządkowany malejąco, uzupełniony udziałem procentowym sezonów grzewczych dla założonych przedziałów względnej liczby stopniodni. Przedstawione wyniki analiz wskazują na obniżanie się wartości względnej liczby stopniodni wraz z kolejnymi sezonami grzewczymi. Z analizy zebranych danych wynika, że w rozpatrywanych sezonach grzewczych wystąpiło osiem sezonów o średniej dobowej temperaturze otoczenia niższej od te = -18oC, jednak występowanie minimalnych dobowych temperatur, na poziomie niższym niż te = – 20oC, dotyczy tylko pięciu z nich. Wyznaczone średnie pięciodobowe temperatury otoczenia dla tych sezonów grzewczych wykazały, że w żadnym przypadku nie przekroczona została temperatura obliczeniowa III strefy klimatycznej (te = – 20oC). Przeprowadzona analiza wykazała dodatkowo, że największą różnicę pomiędzy temperaturami maksymalnymi i minimalnymi zaobserwowano dla miesiąca stycznia na poziomie Δt = 13,6oC oraz lutego Δt = 13,5oC, najniższą różnicą na poziomie Δt = 7,1oC charakteryzował się miesiąc wrzesień. Po odrzuceniu dwóch sezonów z minimalną oraz dwóch sezonów z maksymalną wartością względnej liczby stopniodni, przedział wartości tej wielkości mieści się pomiędzy 0,9 a 1,1, tj. odchylenie od wartości średniej wynosi ± 10 proc. Dokonane oceny wskazały także wyraźną tendencję wzrostu wartości sezonowej temperatury otoczenia oraz tendencje obniżania się wartości względnej liczby stopniodni wraz z kolejnymi sezonami grzewczymi.

Rys. 4. Względna liczba stopniodni:
a) w bieżących sezonach grzewczych, b) uporządkowana malejąco [5]

3. Efekty aktualizacji tabel regulacyjnych
Optymalizację temperatur wody sieciowej prowadzi się poprzez obliczenia zdążające do utworzenia nowych tabel regulacyjnych dla źródeł, sieci i użytkowników systemów ciepłowniczych, umożliwiających ich dostosowanie do aktualnych potrzeb odbiorców przy jednoczesnym zagwarantowaniu niezawodności dostaw ciepła. Właściwie ustalona temperatura wody sieciowej, w zależności od parametrów powietrza zewnętrznego, jest zatem bardzo istotnym parametrem decydującym o możliwościach zwiększania efektywności energetycznej systemów ciepłowniczych. Umożliwia ona bowiem, przy jednoczesnym zapewnieniu komfortu cieplnego odbiorców, minimalizację strat ciepła, oszczędność surowców energetycznych, jak i również poprawę wskaźników ekologicznych. Tabele regulacyjne wody sieciowej dla sieci ciepłowniczej przy regulacji centralno-jakościowej systemu ciepłowniczego opracowywane były na podstawie wytycznych ministerialnych [6]. Zgodnie z tymi wytycznymi temperatury wody sieciowej w źródle ciepła były zależne od obliczeniowej wartości współczynnika obciążenia cieplnego i udziału ciepłej wody użytkowej Ucw. Tabele posiadały dziewięć kolumn i uwzględniały temperaturę powietrza zewnętrznego, nasłonecznienie i prędkość wiatru. Ta podstawowa metoda ustalania parametrów wody sieciowej wykorzystywana przy regulacji systemu ciepłowniczego jest powszechnie stosowana dla sieci ciepłowniczych obecnie pracujących, gdzie jak wcześniej stwierdzono występują znaczące rozbieżności pomiędzy projektowym a rzeczywistym zapotrzebowaniem na ciepło. Należy dodać, że z uwagi na zastosowanie w szerokim zakresie automatyzacji węzłów cieplnych, wykres centralnej regulacji ma w rzeczywistości charakter ilościowo – jakościowy, gdyż zmiany zapotrzebowania na moc cieplną w poszczególnych węzłach ciepłowniczych kompensowane są poprzez zmianę przepływu wody sieciowej.

Dlatego też, dla potrzeb poprawy efektywności energetycznej systemu ciepłowniczego coraz częściej w źródłach tych systemów stosuje się regulacje parametrów sieci w oparciu o prognozowane zapotrzebowanie na ciepło. W tym celu wykorzystuje się rozkład parametrów pracy sieci uzyskany z pomiarów przeprowadzonych we wcześniejszych sezonach grzewczych w funkcji temperatury radiacyjno-efektywnej (TRE)1.. W przypadku regulacji jakościowej dostosowanie mocy cieplej do zapotrzebowania przy zmiennych warunkach atmosferycznych, przy stałym przepływie wody sieciowej w ciągu całego sezonu grzewczego, osiągane jest poprzez odpowiednio dostosowane temperatury na zasilaniu (wyjściu ze źródła ciepła). Natomiast temperatura wody sieciowej na powrocie jest wynikiem ilości odebranego ciepła przez użytkowników.

Obecnie brak jest jasnych metod określania temperatur wody sieciowej. Nowe koncepcje sterowania siecią realizowane są w oparciu o indywidualne modele matematyczne opisujące pracę sieci w rzeczywistych warunkach pogodowych na danym terenie, tzn. powstałe modele uwzględniają specyfikę zachowania się odbiorców końcowych dla danej sieci ciepłowniczej, dając tym samym podwaliny pod inteligentne sieci ciepłownicze. Obecnie dąży się do pełnej automatyzacji węzłów, według planu lub idei pracy sieci. W dostępnej literaturze naukowo-technicznej nie proponuje się żadnych jasnych i precyzyjnych metod tworzenia nowych tabel regulacyjnych wody sieciowej w systemach ciepłowniczych pracujących w obecnych realiach rynku ciepła. Dlatego też, w bieżącej pracy zaproponowano zestaw zależności, które można stosować w rozwiązywaniu problematyki aktualizacji tabel regulacyjnych. Rozkład temperatur wody sieciowej opisany zależnościami 1-6 opracowano w oparciu o pracę [7] oraz analizę rozkładu rzeczywistych temperatur wody sieciowej w funkcji TRE dla kliku krajowych systemów ciepłowniczych. Podstawowe dane wymagane do przeprowadzenia obliczeń rozkładu temperatur wody sieciowej przedstawiono na rysunku 5.

gdzie:

gdzie:

tws.z – temperatura wody w rurociągu zasilającym w źródle ciepła, oC,
tws.p – temperatura wody w rurociągu powrotnym w źródle ciepła, oC,
twi.z.obl – obliczeniowa temperatura wody instalacyjnej na zasilaniu w węźle cieplnym, oC,
twi.p.obl – obliczeniowa temperatura wody instalacyjnej na powrocie w węźle cieplnym, oC,
tws.z.obl – obliczeniowa temperatura wody sieciowej na zasilaniu w węźle cieplnym, oC,
tws.p.obl – obliczeniowa temperatura wody sieciowej na powrocie w węźle cieplnym, oC,
tw – temperatura wewnętrzna w pomieszczeniach ogrzewanych, oC,

u – współczynnik redukcji temperatury w węźle cieplnym, -,
αo – względny przepływ wody sieciowej, -,
δtobl – obliczeniowa różnica temperatur zasilania i powrotu instalacji centralnego ogrzewania w węźle cieplnym, K,
∆tobl – obliczeniowy spadek temperatury w instalacji centralnego ogrzewania, oC,
∆tws,z – obniżenie temperatury wody sieciowej dostarczanej ze źródła do danego przyłącza wskutek strat ciepła podczas przesyłania, K,
∆tws,p – obniżenie temperatury wody sieciowej dostarczanej od danego przyłącza do źródła ciepła wskutek strat ciepła podczas przesyłania, K, ∆tw,p – różnica temperatur wody sieciowej i instalacyjnej na powrocie wymagana warunkami techniczno-eksploatacyjnymi wymiennika ciepła, K,
φo – względne zapotrzebowanie na ciepło, -,

Rys. 5. Podstawowe dane wymagane do obliczeń rozkładu temperatur czynnika grzewczego dla sieci ciepłowniczej

Poniżej przedstawiono możliwe do uzyskania efekty energetyczne, ekologiczne i ekonomiczne obniżenia temperatur wody sieciowej na zasilaniu i powrocie dla wybranego systemu ciepłowniczego o mocy cieplnej 340 MW oraz długości łącznej magistral 28 km. Na rysunku 6 przedstawiono uzyskany efekt energetyczny i ekonomiczny. Na podstawie badań stwierdzono, że:
• obniżenie temperatury zasilania z 135 do 120oC oraz powrotu z 70 do 60oC pozwala na obniżenie mocy zamówionej i rocznego zużycia nieodnawialnej energii pierwotnej pokrywających straty ciepła na przesyle nośnika ciepła o 12 proc.,
• łączna kwota szacowanych oszczędności z tytułu obniżenia strat ciepła sieci ciepłowniczych, poprzez obniżenie temperatur wody sieciowej, może wynieść w tym przypadku 3146 PLN/MW dla standardowego sezonu grzewczego.

Na rysunku 7 przedstawiono jednostkowy efekt ekologiczny2 obniżenia strat przesyłu sieci ciepłowniczej pracującej przy temperaturach 135/70oC oraz wyjściowych stratach ciepła na przesyle wynoszących 8 proc. i 15 proc. W obliczeniach jednostkowego efektu ekologicznego przyjęto emisje jednostkowe, które przedstawiono w tabeli 3. Szacowany jednostkowy efekt ekologiczny przedstawiony jako emisja równoważnawynikający z obniżenia temperatury zasilania z 135 do 130oC oraz temperatur powrotu z 70do 65oC wyniósł 5,7 kg/MW mocy zainstalowanej oraz przy obniżeniu temperatury zasilania z 135 do 120oC oraz powrotu z 70 do 60oC wyniósł 14,7 kg/MW mocy zainstalowanej dla standardowego sezonu grzewczego oraz stratach wyjściowych sieci na poziomie 8 proc.



Rys. 6. Efekt energetyczny i ekonomiczny obniżenia strat przesyłu sieci ciepłowniczej:
a) Nominalne parametry sieci: 135/70
oC, straty na przesyle 8 proc;
b) Nominalne parametry sieci: 135/70
oC, straty na przesyle 15 proc.
Tabela 3. Emisje jednostkowe przyjęte do obliczeń efektu ekologicznego

W przypadku strat sieci wynoszących 15 proc. emisja równoważna wyniosła odpowiednio 10,7 kg/MW oraz 27,5 kg/MW. Zmiana rzeczywistych warunków pracy sieci ciepłowniczej po procesach termomodernizacji budynków przyłączonych do niej prowadzi także do możliwości obniżenia emisji ditlenku węgla. Zakładając, że prace związane z obniżaniem parametrów temperaturowych sieci ciepłowniczych nie wymuszają w większości przypadków wykonania prac skutkujących dodatkową emisją CO2, to uzyskany efekt redukcji emisji ditlenku węgla jest też wskaźnikiem poprawy efektywności energetycznej sieci ciepłowniczej.

W analizowanym przypadku (rysunek 7) uzyskano obniżenie emisji CO2 odpowiednio od 3,4 do 8,7 kg/MW dla strat ciepła w sieci 8 proc. oraz od 6,4 do 16,4 kg/MW dla strat ciepła w sieci 15 proc. w standardowym sezonie grzewczym.



Rys. 7. Efekt ekologiczny obniżenia strat przesyłu sieci ciepłowniczej;
a) Nominalne parametry sieci: 135/70
oC, straty na przesyle 8 proc.;
b) Nominalne parametry sieci: 135/70
oC, straty na przesyle 15 proc.

4. Podsumowanie
Poprawa efektywności energetycznej powinna być realizowana na każdym etapie systemu, przy czym należy pamiętać aby te działania były bezwzględnie ze sobą skoordynowane. Przykładowo rozpoczynając proces modernizacji źródeł wytwarzania należy brać pod uwagę również kierunki rozwoju technologii w sektorze użytkowników.
Zwiększenie efektywności energetycznej sieci ciepłowniczej można zrealizować poprzez procesy, które nie wymagają znacznych kosztów inwestycyjnych, z jednej strony, i są wymuszane bieżącą eksploatacją sieci, z drugiej strony. Wraz z dynamicznym rozwojem budownictwa energooszczędengo (procesy termomodernizacyjne istniejących budynków oraz nowe standardy energetyczne dla budynków nowych) istnieje pilna potrzeba obniżania temperatur pracy sieci ciepłowniczych, dostosowując je w ten sposób do aktualnych potrzeb użytkowników.
Należy jednak podkreślić, że charakter prac optymalizacyjnych jest zależny od indywidualnych cech systemu ciepłowniczego. Dlatego też, działań tych nie można uogólniać podając rozwiązania, które można zastosować w wielu systemach ciepłowniczych. Są to rozwiązania adekwatne do danej specyfiki systemu ciepłowniczego. Ponadto, obniżenie temperatury wody sieciowej może powodować konieczność wprowadzenia zmian w układach pompowych w ciepłowniach i przepompowniach sieciowych, czy konieczność doregulowania strumienia czynnika grzejnego w węzłach u użytkowników.

PRZYPISY:
1TRE – temperatura efektywno-radiacyjna. Wskaźnik uwzględniający wpływ natężenia całkowitego promieniowaniasłonecznego oraz temperaturę, wilgotność względną powietrza i prędkość wiatru.
2Jednostkowy efekt ekologiczny – różnica emisji zanieczyszczeń przed i po usprawnieniu odniesiona do jednostki mocy zainstalowanej.
3Emisja równoważna – emisja zastępcza, która wynika z zsumowania rzeczywistych emisji poszczególnych rodzajów zanieczyszczeń pochodzących z danego źródła pomnożonych przez ich współczynniki toksyczności.

LITERATURA:
1. JURKIEWICZ A.: Decentralizacja systemów ciepłowniczych, jako metoda na obniżanie kosztów produkcji i dostawy ciepła i możliwość zastosowania odnawialnych źródeł energii i urządzeń rozproszonej energetyki (OZE/URE). XV Forum Ciepłowników Polskich, 18-21 września 2011, Międzyzdroje, 2011
2. KOPICA J., SEKRET R.: Efektywność energetyczna miejskich sieci ciepłowniczych. Rynek ciepła 2011. Materiały i studia, Wydawnictwo KAPRINT, Lublin, 2011
3. LIS P., SEKRET R.: Analiza porównawcza sezonowego zużycia ciepła do ogrzewania budynków edukacyjnych. Budownictwo niskoenergetyczne, IX Międzynarodowe Seminarium Naukowo-Techniczne. Problemy projektowania, realizacji i eksploatacji budynków o niskim zapotrzebowaniu na energię ENERGODOM’2008, str. 305-312, Kraków, 2008
4. LIS P., SEKRET R.: Analiza porównawcza sezonowego zużycia ciepła do ogrzewania budynków edukacyjnych. Czasopismo Techniczne – Budownictwo, nr 1-B/2009 Zeszyt 5, str. 167-174, Wydawnictwo Politechniki Krakowskiej, Kraków, 2009
5. SEKRET R., WILCZYŃSKI W.: Analiza zmian temperatury powietrza zewnętrznego oraz długości sezonu grzewczego na liczbę stopniodni na przykładzie miasta Częstochowa. Rynek Energii, Nr 4 (95), str. 58-63, 2011
6. DECYZJA Nr 4 Ministra Gospodarki Materiałowej i Paliwowej z dnia 4 czerwca 1987 r. w sprawie ustalenia temperatury wody sieciowej w źródłach ciepła i systemach ciepłowniczych.
7. SZKARKOWSKI A., ŁATKOWSKI L.: Ciepłownictwo. Wydawnictwo Naukowo-Techniczne, Warszawa, 2006
prof. dr hab. inż. Robert Sekret, Wydział Inżynierii Środowiska i Biotechnologii; Katedra Ciepłownictwa, Ogrzewnictwa i Wentylacji, Politechnika Częstochowska

źródło: nowoczesnecieplownictwo.pl

Warsaw

Banner 468 x 60 px