inż. Mariusz KOPKA
Właściwy dobór dolnego źródła ciepła daje gwarancję wysokiej efektywności działania pompy ciepła. Należy pamiętać o najważniejszej roli instalacji pompy ciepła jaką jest zapewnienie wystarczającej ilości energii cieplnej do ogrzania budynku. Na opłacalność całej inwestycji w dłuższym przedziale czasu ma też duży wpływ zdolność dolnego źródła do regeneracji, jego koherentność, sposób wykonania instalacji czy konieczność uzyskania pozwoleń w przypadku pionowego GWC. Jeśli chodzi o możliwość powrotu instalacji do temperatury początkowej w następnym sezonie grzewczym, to w przypadku kolektorów pionowych może dojść do spadku temperatury.
ZAŁOŻENIA PROJEKTOWE
Obliczenia projektowe zostały wykonane dla wolnostojącego budynku jednorodzinnego znajdującego się w Malborku, którego charakterystykę cieplną przedstawiono w tabeli 1. Występujące w nim straty energii cieplnej ograniczono przez wymianę stolarki okiennej na okna z szybą podwójną oraz zastosowanie warstwy izolacji cieplnej.
Parametr | Oznaczenie | Jednostka | Typ systemu ogrzewania w budynku | ||
konwekcyjne | ścienne | podłogowe | |||
Strefa klimatyczna | [-] | 2 | |||
Projektowa temperatura zewnętrzna |
Θe | [oC] | -18 | ||
Średnia roczna temperatura zewnętrzna |
Θm,e | [oC] | 7,9 | ||
Powierzchnia ogrzewana budynku |
Ah | [m2] | 150,0 | 147,5 | 150,0 |
Kubatura ogrzewana budynku |
Vh | [m3] | 356,5 | 349,6 | 340,7 |
Projektowe obciążenie cieplne budynku |
ΦHL | [W] | 6963 | 6961 | 7236 |
Temperatura zasilania obiegu grzewczego |
THV | [oC] | 55 | 45 | 35 |
Przedstawiony przykład obliczeniowy dotyczy urządzenia charakteryzującego się temperaturą zasilania obiegu grzewczego THV= 55oC, która odpowiada konwekcyjnemu systemowi ogrzewania budynku.
OKREŚLENIE CHARAKTERYSTYKI SPRĘŻARKOWEJ POMPY CIEPŁA
W celu wyznaczenia parametrów sprężarkowej pompy ciepła oraz instalacji dolnego źródła ciepła posłużono się modelem obliczeniowym oraz wytycznymi technologii wykonania gruntowych wymienników ciepła opartymi na zaleceniach PORT PC oraz firmy Viessmann. Obliczenia teoretycznej mocy grzewczej sprężarkowej pompy ciepła odpowiadającej projektowemu obciążeniu cieplnemu budynku wykonano w programie Audytor OZC.
Wartość rzeczywistej mocy grzewczej została powiększona w porównaniu do teoretycznej ze względu na przerwy w pracy urządzenia wymuszone koniecznością regeneracji gruntu. Przyjęto, że czas przerwy w pracy sprężarkowej pompy ciepła będzie wynosił Tp= 10 h.
gdzie:
Tp – czas przerwy w pracy sprężarkowej pompy ciepła, który należy pomniejszyć o dwie godziny ze względu na bezwładność budynku [h].
Na podstawie rzeczywistej mocy grzewczej dobrano sprężarkową pompę ciepła, której charakterystykę techniczną przedstawiono w tabeli 2.
Korzystając ze współczynnika wydajności grzewczej charakteryzującego dobrane urządzenie wyznaczono moc napędową sprężarki:
gdzie:
COP – współczynnik wydajności grzewczej sprężarkowej pompy ciepła [-].
Kolejnym krokiem było wyznaczenie wydajności chłodniczej pompy ciepła:
Znając parametry sprężarkowej pompy ciepła zestawione w tabeli 2, pracującej w trybie grzania określono strumień ciepła doprowadzony do gruntu przez rewersyjną pompę ciepła:
Tabela 2. Charakterystyka techniczna sprężarkowej pompy ciepła Vitocal 200_G BWC 201.A13 firmy Viessmann [1]
Parametr | Oznaczenie | Jednostka | Wartość | ||
Czas pracy sprężarkowej pompy ciepła |
Tsp | [h] | 2900 | ||
Czas przerwy w pracy sprężarkowej pompy ciepła |
Tp | [h] | 10 | ||
Temperatura zasilania obiegu grzewczego |
THV | oC | 55 | 45 | 35 |
Teoretyczna moc grzewcza | [kW] | 6,96 | 6,96 | 7,24 | |
Rzeczywista moc grzewcza | 11,93 | 11,93 | 12,41 | ||
Znamionowa moc grzewcza | 11,98 | 12,46 | 12,95 | ||
Współczynnik wydajności grzewczej COP |
COP | [-] | 3 | 4 | 5 |
Moc sprężarki | Psp | [kW] | 3,99 | 3,12 | 2,59 |
Wydajność chłodnicza | 7,99 | 9,35 | 10,36 | ||
Moc grzewcza w trybie chłodzenia |
15,97 | 15,58 | 15,54 | ||
Spadek temperatury czynnika niskozamarzającego w parowniku |
ΔTp | [K] | 5 | ||
Opór przepływu czynnika niskozamarzającego w parowniku |
ΔPpc | [kPa] | 4,5 |
WYZNACZENIE PARA METRÓW GEO METRYCZNYCH PIONOWEGO GWC
Pierwszy krokiem w wyznaczeniu parametrów geometrycznych pionowego GWC podanych w tabeli 6, było określenie średniego współczynnika przewodzenia ciepła, który wyznaczono na podstawie struktury geologicznej gruntu opisanej w tabeli 3.
gdzie:
ui – udział i- tej warstwy gruntu w całej długości pionowego GWC [-], λi – współczynnik przewodzenia ciepła dla i- tej warstwy gruntu [W/mK].
Tabela 3. Współczynnik przewodzenia ciepła dla wybranej struktury geologicznej gruntu [2]
Rodzaj gruntu | Udział warstwy gruntu w całej długości pionowego GWC |
Współczynnik przewodzenia ciepła |
u | λ | |
[-] | [W/mK] | |
Piasek lekki (5%wody) | 0,25 | 1,35 |
Piasek lekki (15% wody) | 0,15 | 1,55 |
Glina ciężka (5% wody) | 0,40 | 1,2 |
Piasek ciężki (5% wody) | 0,20 | 2,2 |
Znając średni współczynnik przewodzenia ciepła określono jednostkową wydajność cieplną pionowego GWC na podstawie rysunku 1:
Rys. 1. Zależność między jednostkową wydajnością cieplną pionowego GWC qv a współczynnikiem przewodzenia ciepła gruntu λ [3] |
Następnie obliczono długość sondy pionowej:
Przy wyznaczaniu wymaganej długości odwiertu uwzględniono współczynnik korekcyjny zależny od czasu pracy sprężarkowej pompy ciepła.
Zgodnie z zaleceniami o maksymalnej długości odwiertu określono ilość sond w instalacji pionowego GWC.
Znając ilość sond pionowego GWC, obliczono długość pojedynczego wymiennika:
Przystępując do wyznaczenia średnicy wewnętrznej sondy pionowej przyjęto następujące założenia:
- czynnik niskozamarzający stanowi wodny roztwór glikolu propylenowego, którego charakterystykę przedstawiono w tabeli 4;
- przewody rozprowadzające oraz rury doprowadzające wykonane będą z polietylenu sieciowanego (PE-80);
- odległość między sąsiadującymi sondami powinna wynosić minimum 6 m.
Tabela 4. Właściwości cieplno-fizyczne roztworu wodnego glikolu propylenowego [3]
Parametr | Oznaczenie | Jednostka | Wartość |
Gęstość | ρcn | [kg/m3] | 1057 |
Lepkość kinematyczna | υcn | [mm2/s] | 8,30 |
Ciepło właściwe | ccn | [kJ/kgK] | 3,76 |
Zawartość koncentratu w roztworze wodnym |
zk | [%] | 36 |
W celu wyznaczenia średnicy sondy pionowej, określono wydajność chłodniczą przypadającą na jedną sekcję GWC:
oraz strumień objętościowy czynnika niskozamarzającego uwzględniając spadek temperatury nośnika ciepła w parowniku ΔTp= 5 K:
gdzie:
ρcn – gęstość czynnika niskozamarzającego [kg/m3], ccn – ciepło właściwe czynnika niskozamarzającego [kJ/ kgK],
Przy założeniu prędkości przepływu czynnika niskozamarzającego ϑ=0,4 m/s otrzymano średnicę wewnętrzną sondy pionowej, która wyniesie:
Na podstawie otrzymanej wartości oraz tabeli 5 określono średnicę nominalną sondy pionowej:
W przypadku wyznaczania średnicy przewodu doprowadzającego należało przyjąć wydajność chłodniczą przypadającą na cały obieg dolnego źródła pompy ciepła.
Przy zwiększonym strumieniu objętościowym oraz zachowanej prędkości przepływu roztworu wodnego glikolu w przewodzie doprowadzającym w odniesieniu do sondy pionowej nastąpiło zwiększenie średnicy rurociągu zasilającego:
Na podstawie otrzymanej wartości średnicy wewnętrznej oraz tabeli 5 określono średnicę nominalną przewodów doprowadzających, która wynosi:
Posługując się tabelą 5 oraz wyznaczoną średnicą nominalną sondy pionowej i przewodu doprowadzającego określono ilość płynu potrzebną do napełnienia pionowego GWC.
gdzie:
Lw – wymagana długość odwiertu [m], Xs – pojemność rur sondy pionowej [l/m], Lpd – długość przewodów doprowadzających [m], Xpd – pojemność przewodów doprowadzających [l/m].
Tabela 5. Pojemność przewodu rurowego w zależności od średnicy wewnętrznej rury [4]
Średnica zewnętrzna Dz / grubość ścianki gs[mm] |
Średnica nominalna Dn [mm] | Pojemność rury X [l/m] |
20/ 2,0 | 15 | 0,201 |
25/ 2,3 | 20 | 0,327 |
32/ 3,0 | 25 | 0,531 |
40/ 2,3 | 32 | 0,984 |
40/ 3,7 | 32 | 0,835 |
50/ 2,9 | 40 | 1,595 |
50/ 4,6 | 40 | 1,308 |
63/ 5,8 | 50 | 2,070 |
63/ 3,6 | 50 | 2,445 |
Do obliczenia ilości suchej masy wiążącej wykorzystano właściwości materiału wypełniającego odwiert. Znając średnicę nominalną określono średnicę zewnętrzną sondy pionowej posługując się tabelą 5.
gdzie:
gs – grubość ścianki sondy pionowej [m].
W kolejnym kroku wyznaczono objętość materiału wypełniającego dla wszystkich odwiertów GWC:
gdzie:
Do – średnica odwiertu sondy pionowej [m].
Następnie obliczono ilość suchej masy wiążącej przy założonej gęstości materiału ρsmw=810 kg/m3:
Tabela 6. Charakterystyka pionowego GWC dla nierewersyjnej pompy ciepła [opr.wł.]
Parametr | Oznaczenie | Jednostka | Nierewersyjna pompa ciepła |
Rewersyjna pompa ciepła |
||||
Współczynnik wydajności grzewczej COP |
Współczynnik wydajności grzewczej COP |
|||||||
3 | 4 | 5 | 3 | 4 | 5 | |||
Moc sprężarki | Psp | [kW] | 3,99 | 3,12 | 2,59 | 3,99 | 3,12 | 2,59 |
Wydajność chłodnicza/ Moc grzewcza w trybie chłodzenia |
[kW] | 7,99 | 9,35 | 10,36 | 15,97 | 15,58 | 15,54 | |
Łączna długość sondy pionowej |
Lo | [m] | 257,63 | 301,45 | 334,19 | 644,09 | 628,02 | 626,61 |
Wymagana długość odwiertu |
Lw | [m] | 373,57 | 437,10 | 484,58 | 933,92 | 910,64 | 908,59 |
Ilość sond pionowego GWC | Ns | [-] | 4 | 5 | 5 | 10 | 10 | 10 |
Długość pojedynczej sondy pionowej |
Lw,1 | [m] | 93,39 | 87,42 | 96,92 | 93,39 | 91,06 | 90,86 |
UWAGA: W przypadku wyznaczenia parametrów GWC dla rewersyjnej pompy ciepła, uwzględniając zmianę charakteru pracy urządzenia poza okresem grzewczym, obliczenia wykonano w oparciu o strumień ciepła odprowadzany do gruntu.
Po określeniu wymiarów gruntowego wymiennika ciepła dla dobranej sprężarkowej pompy ciepła pracującej w trybie grzewczym oraz w trybie chłodzenia aktywnego, wyznaczono współczynnik zwymiarowania instalacji przedstawiony w tabeli 7 oraz rysunku 2 [5].
Tabela 7. Parametry pionowego GWC dla pompy ciepła pracującej w trybie grzania oraz w trybie chłodzenia aktywnego z uwzględnieniem wskaźnika WZI [opr.wł.]
Parametr | Oznaczenie | Jednostka | Współczynnik wydajności grzewczej COP | ||
3 | 4 | 5 | |||
Wymagana długość odwiertu dla nierewersyjnej pompy ciepła |
Lw | [m] | 373,57 | 437,10 | 484,58 |
Wymagana długość odwiertu dla rewersyjnej pompy ciepła |
Lw,r | [m] | 933,92 | 910,64 | 908,59 |
Współczynnik zwymiarowania instalacji |
WZI | [-] | 2,5 | 2,1 | 1,9 |
Wzrost współczynnika wydajności grzewczej COP powoduje spadek współczynnika zwymiarowania instalacji. Wynika to z mniejszego strumienia ciepła doprowadzonego do gruntu podczas pracy SPC w trybie chłodzenia aktywnego przez ogrzewanie niskotemperaturowe w porównaniu do systemu charakteryzującego się wyższą temperaturą obiegu grzewczego.
Rys. 2. Zależność współczynnika zwymiarowania instalacji w funkcji współczynnika COP [opr.wł.] |
OKREŚLENIE PARAMETRÓW POMPY OBIEGU PIERWOTNEGO
Przed wyznaczeniem mocy pompy obiegu pierwotnego sprężarkowej pompy ciepła przedstawionej w tabeli 8, należało określić straty ciśnienia wywołane oporami hydraulicznymi w przewodach doprowadzających oraz wymiennikach gruntowych. Pierwszy krokiem prowadzącym do określenia liniowych strat ciśnienia w rurach rozprowadzających było wyznaczenie liczby Reynoldsa uwzględniając strumień objętościowy płynu przypadający na jedną sondę.
gdzie:
υcn – lepkość kinematyczna czynnika niskozamarzającego [m2/ s].
Znając charakter przepływu czynnika niskozamarzającego, wyznaczono współczynnik oporów liniowych:
Następnie obliczono prędkość przepływu nośnika ciepła dla założonej średnicy nominalnej rurociągu:
Kolejnym krokiem było wyznaczenie jednostkowych, liniowych strat ciśnienia w sondzie pionowej.
Liczba Reynoldsa określająca charakter przepływu w przewodzie doprowadzającym została obliczona po uwzględnieniu zwiększonego w porównaniu do sondy pionowej strumienia objętościowego płynu oraz średnicy nominalnej rurociągu:
Większa burzliwość przepływu objawiła się mniejszym współczynnikiem oporów linowych.
W celu obliczenia jednostkowych linowych strat ciśnienia w przewodzie doprowadzającym, wyznaczono prędkość przepływu czynnika niskozamarzającego.
Po obliczeniu strat ciśnienia przypadających na jeden metr przewodu wymiennika, wyznaczono liniowe straty ciśnienia GWC.
Kolejnym krokiem było określenie miejscowych strat ciśnienia przy uwzględnieniu współczynnika oporów miejscowych dla rury polietylenowej wynoszącego ζ=0,42.
gdzie:
Σζ – suma współczynników oporów miejscowych występujących w działce [-].
Suma oporów miejscowych została wyznaczona na podstawie liczby zagięć rury wymiennika odpowiadającej ilości sond pionowych:
Tabela 8. Dane techniczne pompy obiegowej dla nierewersyjnej i rewersyjnej pompy ciepła
Przyjmując straty ciśnienia w parowniku sprężarkowej pompy ciepła zgodnie z jej charakterystyką, określono całkowite opory hydrauliczne [1]:
Uwzględniając lepkość kinematyczną czynnika niskozamarzającego we wcześniejszych obliczeniach przy wyznaczaniu rzeczywistej wysokości podnoszenia pompy ciepła oraz strumienia objętościowego przepływu czynnika niskozamarzającego przez przewody doprowadzające, pominięto współczynnik korekcyjny:
gdzie:
g – przyspieszenie ziemskie [m/s2].
Na podstawie obliczonych parametrów dobrano pompę dla obiegu pierwotnego.
PODSUMOWANIE
Określenie obciążenia cieplnego budynku przy wykorzystaniu programu Audytor OZC pozwoliło na wyznaczenie teoretycznej mocy sprężarkowej pompy ciepła. Na podstawie rzeczywistej charakterystyki SPC uwzględniającej czas przerwy w jej pracy dobrano urządzenie. Korzystając ze współczynnika wydajności grzewczej odpowiadającego temperaturze zasilania obiegu grzewczego dla danego typu systemu ogrzewania określono moc sprężarki, wydajność chłodniczą oraz strumień ciepła doprowadzonego do gruntu przez rewersyjną pompę ciepła w trybie chłodzenia aktywnego. Wyznaczenie tych parametrów pozwoliło przejść do obliczenia wymiarów GWC. Poza wydajnością chłodniczą dla nierewersyjnej SPC oraz strumieniem ciepła doprowadzonego do gruntu w przypadku rewersyjnej pompy ciepła, wpływ na parametry geometryczne sondy pionowej mają: czas pracy urządzenia oraz jednostkowa wydajność cieplna gruntu. Zauważono, że dla nierewersyjnej pompy ciepła spadek temperatury zasilania obiegu grzewczego skutkował zwiększeniem wymiarów GWC. Odwrotne zjawisko zaobserwowano w przypadku urządzenia pracującego w trybie chłodzenia aktywnego znajdując potwierdzenie w malejącym wskaźniku zwymiarowania instalacji WZI przy zwiększającym się współczynniku wydajności grzewczej. Przed określeniem charakterystyki pompy obiegowej obiegu pierwotnego SPC, określono straty ciśnienia wywołane oporami hydraulicznymi w przewodach doprowadzających oraz wymiennikach gruntowych. Ze względu na największe wymiary instalacji DŹC rewersyjnej pompy ciepła o współczynniku COP=3, zastosowano pompę obiegową o największym poborze mocy.
cdn...
LITERATURA
1. Wytyczne projektowe pomp ciepła Vitocal Viessmann, 05/ 2012.
2. Raport o oddziaływaniu przedsięwzięcia na środowisko – Zmiana koncecesji dla przedsięwzięcia polegającego na poszukiwaniu i rozpoznawaniu złóż gazu ziemnego w obrębie obszaru koncesji Malbork, URS Polska Sp. z o.o, 02/2013.
3. Wytyczne projektowania, wykonania i odbioru instalacji z pompami ciepła, Wyd. PORT PC, Kraków, 2013.
4. Wytyczne projektowe – Systemy pomp ciepła, Viessmann Sp. z o.o., Wrocław, 05/2005.
5. Koszorek P.: Analiza możliwości regeneracji gruntu za pomocą rewersyjnej pompy ciepła. „Technika Chłodnicza i Klimatyzacyjna”, 10/ 2015.