Hőbányászat települési szilárd hulladéklerakókból X. Fazola-napok, „Hulladék mint nyersanyag” szakmai-tudományos konferencia, Miskolc, MAB székház, 2016. szeptember 16. Magyar Tamás, PhD hallgató Dr. habil. Faitli József, egyetemi docens Miskolci Egyetem, Nyersanyagelőkészítési és Környezeti Eljárástechnikai Intézet
Szakirodalmi áttekintés – A hulladékgazdálkodás jelenlegi helyzete 1,3 Mrd tonna TSZH/év (2012) → 2,2 Mrd tonna TSZH/év (előrejelzés 2025-re). 1,2 kg/fő/nap TSZH (2012) → 1,42 kg/fő/nap (előrejelzés 2025-re). A begyűjtött TSZH közel 95 %-a kerül deponálásra (KURNIAWAN ET AL., 2006). 400
Ártalmatlanított hulladék mennyisége, [millió tonna/év]
Jelmagyarázat:
350
1. Lerakással történő ártalmatlanítás (védelemmel ellátott hulladéklerakó) 2. Újrafeldolgozással hasznosított hulladék 3. Energiahasznosítással történő égetéssel hasznosított hulladék 4. Lerakással történő ártalmatlanítás (műszaki védelem nélküli hulladéklerakó) 5. Komposztálás 6. Egyéb
300
250
200
150
100
50
0
1.
2.
3.
4.
5.
6.
1. ábra: A TSZH ártalmatlanítási lehetőségeinek aránya világviszonylatban (HOORNWEG és BHADA-TATA, 2012 alapján).
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
2/24
Szakirodalmi áttekintés – A TSZH-okban lejátszódó bomlási folyamatok A lebomlási folyamatok eredményeként keletkező úgynevezett „melléktermékek”: Csurgalékvíz Biogáz Hő
2. ábra: Hő- és gázképződés egy TSZH lerakóban (COCCIA ET AL., 2013).
3. ábra: Hőképződés az idő függvényében (HANSON ET AL., 2013).
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
3/24
Szakirodalmi áttekintés – A kinyerhető energia nagysága (YOUNG, 1992) YOUNG, 1992-es tanulmánya szerint az energiaképződésre vonatkozó képlet az alábbi formában írható fel:
Tkörnyezet
chulladék Thulladék
n
E Δt i c w M(t)i i 1
Ahol: E: Hőképződés [MJ/m3] Δti: hőmérsékletkülönbség [K] cw: hulladék hőkapacitása [MJ/m3K] M(t)i: a felszabaduló energia azon hányada, amely a lerakó fűtésére fordítódik
4. ábra: A hőképződés természettudományos megközelítése (YOUNG, 1992).
Megállapítás: A természettudományos megközelítés szerint a lebomlás következtében energia keletkezik, amely az anyag felmelegedését, - így a belső energia növekedését - eredményezi.
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
4/24
Szakirodalmi áttekintés – A kinyerhető energia nagysága (YESILLER ET AL., 2015) YESILLER ET AL., 2015-es tanulmánya szerint az energiaképződésre vonatkozó képletét az 1 éven keresztül mért hőmérséklet-eloszlás alapján becsülte:
Depóniagáz: ∆T(átlag) depógáz x c(átlag)=5,2 MJ/m3 Talaj: ∆T(átlag) talaj x c(átlag)=48 MJ/m3
Ahol: ∆Tátlag – időbeli és a teljes lerakóra térbeli átlag
6. ábra: A kinyerhető energia becslésének elméleti háttere (YESILLER ET AL., 2015).
5. ábra: A vizsgált TSZH lerakó hőmérséklet-eloszlása (2D) (YESILLER ET AL., 2015).
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
5/24
Vizsgálati helyszín – gyáli regionális hulladéklerakó Az .A.S.A. Magyarország Kft. telephelye Gyálon található. A lerakóban megközelítőleg 200 000 tonna kevert TSZH kerül deponálásra évente. A lerakó 5 ütemre osztható fel, jelenleg a beérkező TSZH-ot az 5. ütemben deponálják.
Ütem azonosító
Befogadó kapacitás [m3]
Deponálás kezdete
1.
377 596
1999
2.
426 322
2003
3.
593 059
2006
4.
400 000
2009
5.
423 900
2012
1. táblázat: Az egyes ütemekhez tartozó kapacitások és a deponálás időtartama.
7. ábra: A gyáli regionális hulladéklerakó látképe.
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
6/24
Kísérleti eszközök fejlesztése – Hőmérséklet-monitoring rendszer kiépítése A monitoring rendszer kiépítésének első lépéseként 10 fúrási pontot jelöltünk ki a depónia első 4 ütemében, amelyekben a hőmérséklet 16 méteres mélységig hosszútávon vizsgálható (1 mes osztásközzel). A hőmérsékletérzékelő szenzorok elhelyezése a depótestet alkotó agresszív közegben nagy kihívást jelentett. Számos tényező mérlegelése után úgy döntöttünk, hogy KPE csőből védőburkolatot kell a depótestbe fúrt lyukba helyezni, amelybe a hőmérő szondák leengedhetők.
Ütem azonosító
Hőmérséklet-monitoring kutak száma
Hulladék kora
1.
2
12-16
2.
3
9-12
3.
2
6-9
4.
3
0-6
2. táblázat: A kiépített hőmérséklet-monitoring kutak száma és hulladék kora az egyes ütemekre vonatkozóan.
8. ábra: Furatok készítése a depóniában.
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
7/24
Kísérleti eszközök fejlesztése – TSZH hőtani paramétereinek vizsgálata A berendezés segítségével mért és számított paraméterek: Fizikai: minta tömege (m), térfogata (V), nedvességtartalom (nm), halmazsűrűség (ρB), száraz tömeg (msg) és térfogat (Vsg), folyadék- (εl), szilárd- (εs), és gáz (εg) fázisok térfogatarányai. Hőtani: hővezetési tényező (λ), fajhő (c), hődiffuzivitás (κ). Hőmérsékletmérő szenzorok
Hőáramsűrűség mérő szenzorok
9. ábra: A TSZH hőtani paramétereinek mérésére kifejlesztett mérőberendezés.
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
8/24
Eredmények – Hőmérséklet-eloszlás (idősebb hulladék esetében) 60
50
Jelmagyarázat -15 m (T = 55,5 °C) -14 m (T = 54,6 °C) -13 m (T = 55,1 °C) -12 m (T = 53,4 °C)
30
-11 m (T = 52,2 °C) -10 m (T = 51,2 °C) -9 m (T = 49,3 °C) -8 m (T = 47,6 °C)
20
-7 m (T = 45,3 °C) -6 m (T = 43,9 °C)
Azonosító: II. ütem Hulladék kora: 9-12 év Lebomlási fázis: anaerob Hőmérséklet trend: állandósult állapot Illesztés: lineáris
Külső hőmérséklet (T = 14,5 °C) Illesztett lináris egyenes (-15 m) Illesztett lináris egyenes (-11 m) Illesztett lináris egyenes (-6 m)
10
0
Illesztés egyenlete: Y = -A × x + B
01.11.2014
01.10.2014
01.09.2014
01.08.2014
01.07.2014
01.06.2014
01.05.2014
01.04.2014
01.03.2014
01.02.2014
01.01.2014
01.12.2013
01.11.2013
01.10.2013
01.09.2013
01.08.2013
-10
19.06.2013 01.07.2013
Hőmérséklet, T [°C]
40
Dátum, [nn.hh.éé.] 10. ábra: Hőmérséklet monitoring eredmények a II/2-es kút esetében.
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
9/24
Eredmények – Hőmérséklet-eloszlás (friss hulladék esetében) 70
60
50
-14 m (T = 55,7 °C)
40
-13 m (T = 53,7 °C) -12 m (T = 52,7 °C) -11 m (T = 53,1 °C) -10 m (T = 52,1 °C)
30
-9 m (T = 50,9 °C) -8 m (T = 50,5 °C)
Azonosító: IV. ütem Hulladék kora: 0-6 év Lebomlási fázis: aerob és anaerob is (mélység függő) Hőmérséklet trend: felmelegedési szakasz Illesztés: exponenciális
20
10
0
-7 m (T = 49,8 °C) -6 m (T = 49,4 °C) Külső hőmérséklet (T = 13,8 °C) Illesztett exponenciális görbe (-15 m) Illesztett exponenciális görbe (-11 m) Illesztett exponenciális görbe (-6 m)
Illesztés egyenlete: Y = T × (1 - A × EXP (-B × x)) 17.07.2014
01.07.2014
01.06.2014
01.05.2014
01.04.2014
01.03.2014
01.02.2014
01.01.2014
01.12.2013
01.11.2013
01.10.2013
01.09.2013
01.08.2013
01.07.2013
-10
19.06.2013
Hőmérséklet, T [°C]
Jelmagyarázat -15 m (T = 55,9 °C)
Dátum, [nn.hh.éé.] 11. ábra: Hőmérséklet monitoring eredmények a IV/3-as kút esetében.
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
10/24
Eredmények – TSZH hőtani paramétereinek vizsgálata (λƩ 3 fázis-ra) A térfogathányadok (koncentrációk) összege: ε s ε l ε g 1 Soros elrendezés: elméleti minimum a λ-ra vonatkozóan.
λ
λs λl λg εs λ l λ g ε l λs λ g ε g λs λl
Párhuzamos elrendezés: elméleti maximum a λ-ra vonatkozóan.
λ λs ε s λ l ε l λ g ε g
12. ábra: A soros- és párhuzamos hővezetés elméleti háttere 3 fázis esetében.
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
11/24
Eredmények – TSZH hőtani paramétereinek becslése (3 fázisra) Fázis Folyadék (víz) Gáz (levegő) Szilárd (TSZH)
Sűrűség [kg/dm3] ρl = 1 ρg = 0,0012 ρs = 1,297
Hővezetési tényező [W/mK] λl = 0,6 λg = 0,025 λs = 3,99
Fajhő [J/gK] cl = 4,181 cg = 1,012 cs = 1,8
3. táblázat: Az egyes fázisok eredő hőtani jellemzőinek becslése.
Szabványos hulladék kategória 1. Biológiailag lebomló 2. Papír 3. Karton 4. Kompozit 5. Textil 6. Higiéniai 7. Műanyag 8. Éghető egyéb 9. Üveg 10. Fém 11. Nem éghető egyéb 12. Veszélyes Finom (< 20 mm) Szilárd fázis eredője:
Mért tömeg hányad [%] 21,6 12,7 4,7 2,1 3,6 4,4 19,9 2,9 3,6 3,6 4,4 0,7 15,7 100 %
Becsült sűrűség [kg/dm3] 0,7 1,1 1,1 1,1 0,25 1 1,15 1,2 2,7 4 2,5 3,5 1,5 1,297
Becsült hővezetési tényező [W/mK] 0,15 0,05 0,21 0,2 0,06 0,155 0,25 0,2 0,87 91,92 1,59 56,22 0,16 3,99
Becsült fajhő [J/gK] 2,721 1,34 1,34 2,19 1,29 1,48 1,67 2,19 0,837 0,483 0,858 1,1 2,112 1,8
4. táblázat: A szilárd fázis hőtani jellemzőinek becslése a szabványos (MSZ 21420/28 és 29) hulladékösszetételi vizsgálat eredményei alapján.
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
12/24
Eredmények – TSZH hőtani paramétereinek vizsgálata (λƩ) Szaturáció határértéke:
Jelmagyarázat P=440 W, B=350-450 kg/dm3
2.5
P=1760 W, 350-450 kg/dm3 2.25
P=440 W, 450-550 kg/dm3
Eredő hővezetési tényező, [ W/mK ]
P=1760 W, 450-550 kg/dm3 P=440 W, 550-750 kg/dm3
2
P=1760 W, 550-750 kg/dm3 400 kg/dm3
1.75
500 kg/dm3 650 kg/dm3
1.5
Eredő hővezetési tényező: párhuzamos elrendezés (elméleti maximum)
1.25
1
0.75
0.5
Eredő hővezetési tényező: soros elrendezés (elméleti minimum)
0.25
0 0
0.1
0.2
0.3
0.4
0.5
Folyadék fázis térfogataránya, l [ - ] 13. ábra: A mért és a számított hővezetési tényezők.
0.6
εl
ρB ρl
Amennyiben εl értéke meghaladja a szaturáció hatáértékét, úgy abban az esetben csak 2 fázisú rendszerről (szilárd-folyadék) beszélünk. Hővezetési Szakirodalmi forrás tényező, λ [W/mK] HOUI et al. (1997) 0,10 YOSHIDA et al. (1997) 0,53 ZANETTI et al. (1997) 0,0445 LEFEBVRE et al. (2000) 0,10 HANSON et al. (2006) 0,30 HANSON et al. (2008) 0,6-1,5 ROWE et al. (2010) 0,35-0,96 BONANY et al. (2013) 0,67 FAITLI et al. (2015) 0,24-1,15
Fajlagos hőkapacitás, c [J/kgK] 1900-3000 3300 2200 1900-3000 719 1200-2200 1940-2360 1400 900-2100
5. táblázat: A TSZH hőtani paraméterei a szakirodalom alapján.
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
13/24
A megépített technológia – Hőkinyerő- és hőhasznosító körök
14. ábra: A hőkinyerő- és hőhasznosító technológia elvi vázlata (H - hideg, előremenő csőág; M - meleg, visszatérő csőág).
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
14/24
A megépített technológia – Hőkinyerő körök „Hurkos” „Piskóta”
15. ábra: A hőkinyerő körök (vízszintes: balra; függőleges: jobbra).
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
15/24
A megépített technológia – Hőhasznosító körök Csurgalékvíz-gyűjtő medence Üvegház
16. ábra: A hőhasznosító körök (csurgalékvíz-gyűjtő medence: balra; üvegház: jobbra).
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
16/24
Hőcserélési kísérletek – Elvégzett vizsgálatok Mérés jelölése 2014/I.
Hőkinyerő és hőhasznosító Függőleges a csurgalékvíz medencére
Hőkinyerés stratégiája
Időtartam
intenzív hőkinyerés, hosszú regeneráció alacsony intenzitással stacionér állapot elérése
2014.08.21.2014.09.15. 2014.09.17.2014.10.13. 2014.11.01.2014.11.25. 2014.11.26.2014.12.11. 2014.12.12.2014.12.20. 2015.06.01.2015.06.08. 2015.06.09.2015.06.19. 2015.06.19.2015.07.07. 2015.07.08.2015.07.16. 2015.07.17.2015.07.31.
2014/II.
Függőleges az üvegházra
2014/III.
„Hurkos” a csurgalékvíz medencére
szakaszos működés, változó intenzitás
2014/IV.
„Piskóta” az üvegházra
szakaszos működés, változó intenzitás
2014/V.
Függőleges az üvegházra
intenzív hőkinyerés
2015/I. 2015/II. 2015/III.
„Piskóta” a csurgalékvíz medencére Függőleges a csurgalékvíz medencére „Piskóta” a csurgalékvíz medencére
2015/IV.
Függőleges az üvegházra
2015/V.
„Piskóta” a csurgalékvíz medencére
kis intenzitás, folyamatos üzem kis intenzitás, napi váltásban a kinyerés és regeneráció mérsékelt intenzitás, hosszú aktív és passzív szakasz mérséklet intenzitás, rövidebb aktív és passzív szakasz változó intenzitás, hosszú aktív és passzív szakasz
3. táblázat: Az elvégzett félüzemi méretű hőkinyerési- és hőhasznosítási vizsgálatok főbb jellemzői.
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
17/24
Hőcserélési kísérletek – Hőmérő szenzorok elhelyezkedése a kutakban
Fa apríték a hőszigetelés biztosítása érdekében (λfa apríték≈ 0,14 W/mK)
Beton a jó hővezetés és a mechanikai stabilitás érdekében (λbeton≈ 1,09 W/mK)
17. ábra: A hőmérő szenzorok elhelyezkedése függőleges hőcserélő körben („A és „B” kutak esetében).
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
18/24
Hőcserélési kísérletek – Hőkinyerés az „A” jelű függőleges kútból (2014/I.) 60
Hőkinyerési szakasz: 9085 perc (≈ 6 nap) Átlagos hőáram: 1152 W Kinyert hőenergia: 0,63 GJ A lebomlásból származó fajlagos hőteljesítmény: 0,53 W/m3
55
Hőmérséklet, T [°C]
50 45
40 35
30 25
20 15
Hőkinyerési szakasz
10
Regenerációs szakasz
8000 Jelmagyarázat A1 A2 A3 A4 A6 A7 A8 A9 A10 Hőáram
7000
5000
4000
Regenreációs szakasz: 18180 perc (≈ 13 nap) A lebomlás effektív fajlagos hőteljesítménye: 0,18 W/m3
3000
2000 1000
09.09.14.
08.09.14.
07.09.14.
06.09.14.
05.09.14.
04.09.14.
03.09.14.
02.09.14.
01.09.14.
31.08.14.
30.08.14.
29.08.14.
28.08.14.
27.08.14.
26.08.14.
25.08.14.
24.08.14.
23.08.14.
22.08.14.
0
21.08.14.
Hőáram, q [W]
6000
Dátum, [nn.hh.éé.]
18. ábra: Hőkinyerés az „A” jelű függőleges kút esetében.
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
19/24
Hőcserélési kísérletek – „Csőhéj” modell a hőfejlődés figyelembe vételével A hővezetés differenciál egyenlete (belső energia mérlegegyenlete): dQ dt q λ A dτ dr A hővezetés differenciál egyenletének megoldása hengerszimmetrikus, időben állandó hőáram esetében:
q
2πλh t n t 1 rn ln r1
Azonban a vizsgált hengeren belül folyamatos hőfejlődés van:
q x p rn2 rx2 π h 17. ábra: A „csőhéj” modell.
A fentiek alapján a hőáram- és hőmérséklet-eloszlás: 2 πλ r Δr qx Δt x h q x ln x rx Δr r x ln Δt x 2 πλh rx
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
20/24
Konklúziók
Lebomlás fajlagos hőteljesítménye:
0,53 W/m3 Lebomlás effektív fajlagos hőteljesítménye (=kinyerhető):
0,18 W/m3 1 éven keresztül 0,18 W/m3 fajlagos hőteljesítménnyel kinyerhető hőenergia mennyisége:
5,67 MJ Henger paraméterei: Henger magassága: 16 m Iterálással meghatározott sugár: 6 m Henger térfogata: 1809 m3
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
21/24
Irodalomjegyzék
BONANY, J. E., VAN GEEL, P. J., GUNAY, B. H. és ISGOR, B. O. (2013): Simulating waste temperatures in an operating landfill in Québec, Canada. Waste Management & Research, 31(7), pp. 692-699. COCCIA, C.J.R., GUPTA, R., MORRIS, J. és MCCARTNEY, J.S. (2013): Municipal solid waste landfills as geothermal heat sources. Elsevier, Renewable and Sustainable Energy Reviews 19, pp. 463-474. HANSON, J. L., LIU, W.-L. és YESILLER, N. (2008): Analytical and Numerical Methodology for Modeling Temperatures in Landfills. Selected Sessions of GeoCongress 08: Geotechnics of Waste Management and Remediation, New Orleans, Louisiana, pp. 1-8. HANSON, J. L., YESILLER, N., HOWARD, K. A., LIU, W. -L., és COOPER, S. P. (2006): Effects of Placement Conditions on Decomposition of Municipal Solid Wastes in Cold Regions. Proceedings of the 13th International Conference on Cold Regions Engineering. HANSON, J.L., YESILLER, N., ONNEN, M.T., LIU, W.-L., OETTLE, N.K. és MARINOS, J.A. (2013): Development of numerical model for predicting heat generation and temperatures in MSW landfills. Waste Management 33 (10), pp. 1993-2000. HOORNWEG, D. és BHADA-TATA, P. (2012): What a waste - A global review of solid waste management, Research report of World Bank, March 2012, No. 15, p. 116. HOUI, D., PAUL, E. és COUTURIER, C. (1997): Heat and mass transfer in landfills and biogas recovery. Proc., 6th Int. Waste Management and Landfill Symp., CHRISTENSEN, T.H. et al., eds., Vol. I, CISA, Italy, pp. 101-108. KURNIAWAN, T.A., LO, W.-H. és CHAN, G.Y.S. (2006): Degradation of recalcitrant compounds from stabilized landfill leachate using a combination of ozone-GAC adsorption treatment. J. Hazard. Mater. B137, pp. 433-455. LEFEBVRE, X., LANINI, S., HOUI, D. (2000): The role of aerobic activity on refuse temperature rise. I: Landfill experimental study. Waste Manage. Res., 18(5), pp. 444–452. ROWE, R.K., HOOR, A. és POLLARD, A. (2010): Numerical Examination of a Method for Reducing the Temperature of Municipal Solid Waste Landfill Liners. Journal of Environmental Engineering, 136(8), pp. 794-804. YESILLER, N., HANSON, J.L., KOPP, K.B. és YEE, E.H. (2015): Assessing approaches for extraction of heat from msw landfill. In: R., COSSU, P., HE, P., KJELDSEN, Y., MATSUFUJI, D., REINHART, R., STEGMANN (szerk.), 15th International Waste Management and Landfill Symposium: G13. Workshop: Heat utilization from landfills. S. Margherita di Pula: CISA Publisher. YOSHIDA, H., TANAKA, N. és HOZUMI, H. (1997): Theoretical study on heat transport phenomena in a sanitary landfill. Proc., 6th Int. Waste Management and Landfill Symp., CHRISTENSEN, T.H. et al., eds., Vol. I, CISA, Italy, pp. 109-120. YOUNG, A. (1992): Application of computer modelling to landfill processes. DoE Rep. No. CWM 039A/92, Dept. of Environment, London. ZANETTI, M. C., MANNA, L. és GENON, G. (1997): Biogas production Evaluation by Means of Thermal Balances. Proc., 6th Int. Waste Management and Landfill Symp., CHRISTENSEN, T.H. et al., eds., Vol. II, CISA, Italy, pp. 523-531.
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
22/24
Köszönetnyilvánítás
A Depóniahő-hasznosítási technológia kidolgozása (KMR_12-1-2012-0128) megnevezésű projekt a Magyar Kormány támogatásával, a Nemzeti Fejlesztési Ügynökség kezelésében, a Kutatási és Technológiai Innovációs Alap finanszírozásával valósul meg.
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
23/24
Köszönöm a figyelmet!
X. Fazola-napok - „Hulladék mint nyersanyag” szakmai-tudományos konferencia
24/24