Politechnika Wrocławska Rok akademicki: 2013/2014
Instytut Inżynierii Ochrony Środowiska
Zakład Ochrony Atmosfery
Ochrona powietrza 2
Projekt
Temat nr 20
Prowadząca ćwiczenia: Mgr inż. Maria Skrętowicz Wykonała:
Magdalena Grabowska
Spis treści
1. Cel zadania projektowego
2. Zakres projektu.
2.1 Obliczenia wstępne
2.2 Obliczenia rozprzestrzeniania zanieczyszczeń
3. Charakterystyka obiektu.
4. Charakterystyka otoczenia
5. Tło zanieczyszczeń
6. Obliczenia wstępne
6.1.Wydajność cieplna źródła kotła
6.2. Wydajność cieplna kotłowni
6.3. Zużycie paliwa dla kotłowni w poszczególnych okresach
6.4 Unos zanieczyszczeń (godzinowy)
6.5 Emisja substancji gazowych
6.6 Prędkość swobodnego opadania ziaren pyłu
6.7 Określenie składu spalin suchych (przy braku analizy elementarnej paliwa)
6.8 Określenie gęstości spalin suchych w warunkach normalnych
6.9 Określenie ilości powstających spalin przy spalaniu węgla kamiennego lub brunatnego (przy braku analizy elementarnej paliwa)
6.10 Określenie stężeń zanieczyszczeń w spalinach w przeliczeniu na spaliny suche przy zawartości 6% O2
6.11 Prędkości wylotowe spalin
7. Parametry emisyjne przyjęte do obliczeń.
8. Sprawdzenie warunków skróconego zakresu obliczeń.
9. Obliczenia rozkładów przestrzennych stężeń maksymalnych i średniorocznych
10. Wymagane redukcje emisji
11. Obliczenia stężeń maksymalnych na wysokości zabudowy mieszkalnej
12. Obliczenia opadu pyłu
13. Wymagane redukcje emisji
14. Wyznaczenie emisji granicznych według kryterium imisyjnego i emisyjnego oraz proponowane emisje dopuszczalne
15. Podsumowanie i proponowane rozwiązania
16. Izolinie stężeń maksymalnych i średniorocznych przed i po redukcji emisji
Zestawienie tabel
Tabela 1. Dane emitora
Tabela 2. Wartości odniesienia dla tła zanieczyszczeń
Tabela 3. Zestawienie parametrów frakcji emitowanego pyłu
Tabela 4. Emisja średnia dla poszczególnych substancji
Tabela 5. Emisja wariantowa dla poszczególnych substancji
Tabela 6. Pył całkowity
Tabela 7. Wyniki obliczeń
Tabela 8. Podokresy emisji dla poszczególnych substancji
Tabela 9. Przekroczenia stężeń dopuszczalnych
Tabela 10. Analiza najwyższych wartości stężeń
Tabela 11. Emisja w wariantach po redukcji
Tabela 12. Emisja średnia po redukcji
Tabela 13. Najwyższe wartości stężeń maksymalnych I średniorocznych po redukcji emisji
Tabela 14. Stężenia maksymalne na wysokości budynku
Tabela 15.Parametry frakcji pyłu
Tabela 16. Przekroczenia
Tabela 17. Ładunek pyłu całkowitego po redukcji
Tabela 18. Określenie emisji granicznych oraz proponowanych emisji dopuszczalnych
Zestawienie rysunków
Rys. 1. Schemat sieci obliczeniowej.
Spis oznaczeń
Ar – zawartość popiołu w paliwie; %
Ari - liczba Archimedesa dla ziaren danej frakcji pyłu
Bki – zużycie paliwa dla kotłowni w poszczególnych wariantach; kg/h
C - stała Sutherlanda, dla powietrza C=112
CO2max - maksymalna zawartość CO2 w spalinach; dla węgla kamiennego wynosi 18,7%
D1 – wartość odniesienia substancji w powietrzu lub dopuszczalny jej poziom, uśrednienie dla jednej godziny, μg/m3
Da – wartość odniesienia substancji w powietrzu lub dopuszczalny jej poziom, uśrednione dla 1 roku, μg/m3
de - wewnętrzna średnica emitora; m
dpi - średnica przedziałowa ziaren pyłu; m
Ei – emisja zanieczyszczenia w danym wariancie; kg/h
Ered – emisja zredukowana, kg/h
Edop – emisja dopuszczalna, kg/h
Fe - pole przekroju wylotowego emitora; m2
g - przyspieszenie ziemskie; m/s2
he – wysokość geometryczna emitora
Hmax – najwyższa efektywna wysokość emitora w zespole z obliczonych dla wszystkich sytuacji meteorologicznych, m
ki - obciążenie kotła; %
M – masa cząsteczkowa substancji
n – liczba kotłów pracujących w danym wariancie z obciążeniem ki
P – zawartość części palnych w pyle; %
R – tło substancji gazowej, μg/m3
Rp – tło opadu pyłu, g/(m2rok)
Sa – wartość dyspozycyjna, μg/m3
QKi - wydajność cieplna kotłowni; kW
QN – nominalna wydajność kotła, MW
Repszi - szacunkowa wartość liczby Reynoldsa dla ziaren danej frakcji pyłu
Repi - liczba Reynoldsa dla ziaren danej frakcji pyłu
s - zawartość siarki w paliwie; %
ti – czas trwania wariantu; h
Tot - temperatura otoczenia; K
Ts – temperatura spalin; K
upi - prędkość opadania ziaren danej frakcji pyłu; m/s
Wd - wartość opałowa paliwa; kJ/kg
v – prędkość gazów odlotowych na wylocie emitora, m/s
x- zawartość wilgoci w spalinach suchych; g/kgss
Oznaczenia greckie:
η - sprawność cieplna kotła
ηSO2 – skuteczność odsiarczania; %
ηp – całkowita skuteczność odpylania; %
ηpi – frakcyjna skuteczność odpylania; %
ηpow - dynamiczny współczynnik lepkości powietrza w warunkach rzeczywistych; Pa*s
η0pow - dynamiczny współczynnik lepkości powietrza w warunkach normalnych =17,04*10-6 Pa*s
λ - współczynnik nadmiaru powietrza
µUi – udział poszczególnej frakcji w unosie pyłu ogółem
ρpow - gęstość powietrza w warunkach rzeczywistych; kg/m3
ρ0pow - gęstość powietrza w warunkach normalnych (powietrze suche)=1,293 kg/m3n
ρp - gęstość pyłu; kg/m3
ρSSn - gęstość spalin suchych w warunkach normalnych; kg/m3n
Cel zadania projektowego
Celem zadania projektowego jest ocena stanu zanieczyszczenia powietrza przez kotły oraz wyznaczenie wymaganego stopnia i sposobu redukcji zanieczyszczeń dla zadanej kotłowni.
Zakres zadania projektowego
Projekt został wykonany według poniższego schematu:
Obliczenia wstępne
Obliczenia emisji godzinowych SO2, NO2, CO, pyłu ogółem, pyłu zawieszonego PM10 i pozostałych frakcji pyłu dla wszystkich wariantów pracy i współpracy kotłów
Obliczenia prędkości opadania ziaren pyłu dla poszczególnych frakcji pyłu
Określenie podstawowego składu i gęstości spalin suchych
Obliczenia ilości spalin wyprowadzanych przez emitor dla poszczególnych wariantów pracy i współpracy kotłów:
- strumienie spalin wilgotnych w warunkach normalnych i rzeczywistych
- strumienie spalin suchych w warunkach normalnych
- strumienie spalin suchych w warunkach umownych przyjętych przy określaniu standardów emisyjnych
Określenie stężeń zanieczyszczeń w spalinach dla warunków umownych przyjętych przy określaniu standardów emisyjnych
Określenie prędkości wylotowych spalin z emitora dla wszystkich wariantów pracy i współpracy kotłów
Obliczenia rozprzestrzeniania zanieczyszczeń
Sprawdzenie warunku skróconego zakresu obliczeń dla substancji gazowych oraz sprawdzenie warunków kryterium opadu pyłu –w przypadku niespełnienia warunku skróconego zakresu dla jakiejkolwiek substancji gazowej przechodzimy do pełnego zakresu obliczeń
Projekt siatki obliczeniowej
Obliczenia stężeń wybranych zanieczyszczeń na poziomie zabudowy mieszkalnej wyższej niż parterowej zlokalizowanej w zasięgu 10 wysokości emitora
Obliczenia opadu pyłu
Określenie emisji granicznych ze względu na kryterium imisyjne
Określenie emisji granicznych ze względu na kryterium emisyjne
Propozycja emisji dopuszczalnych przedstawiona tabelarycznie
Określenie niezbędnego stopnia i sposobu ograniczania emisji – można uwzględnić możliwość zmiany parametrów paliwa i parametrów emitora
Wnioski i zalecenia
Charakterystyka obiektu
Zadanym obiektem jest kotłownia, w której umieszczone są 2 kotły opalane węglem kamiennym o wydajności cieplnej kotła wynoszącej 10 MW z rusztem mechanicznym o sprawności cieplnej wynoszącej η= 82% . W kotłach spalany jest węgiel kamienny o wartości opałowej
Wd = 21 800 kJ/kg. Zawartość popiołu w paliwie wynosi: Ar = 18%, zawartość siarki w paliwie wynosi: s = 1,9%. Spaliny o temp = 425 K, odprowadzane z kotła zawierają wilgoć o wartości wynoszącej
X=50 g/kgs Gęstość pyłu odprowadzanego do atmosfery wynosi 1850 kg/m3. Skład frakcyjny pyłu wynosi:
0-10μm –15 %
10-20μm – 20%
20-40μm – 16%
40-60μm – 19%
60-100μm – 12%
>100μm – 18%
Kotłownię wyposażono w instalację do odpylania o skuteczności odpylania ηuj wynoszącej dla poszczególnych frakcji:
0-10μm – 3%
10-20μm – 6%
20-40μm – 25%
40-60μm – 52%
> 60 μm – 97%
Praca kotłów odbywa się w 3 wariantach. W sezonie letnim w jednym wariancie pracuje jeden kocioł ze średnim obciążeniem w czasie 3790 godzin. W sezonie zimowym praca kotłów odbywa się w dwóch wariantach: dwa kotły z obciążeniem średnim w czasie 3280 godzin lub 2 kotły z obciążeniem maksymalnym w czasie 1300 godzin.
Zanieczyszczenia są emitowane za pomocą emitora o średnicy wylotowej 1 m oraz wysokości 30m. Dane dotyczące emitora przedstawiono w tabeli 1.
Tabela 1. Dane emitora
------------------------------------------------------------------------
ATMOTERM Opole EK100W
------------------------------------------------------------------------
MAGDALENA GRABOWSKA
------------------------------------------------------------------------
DANE EMITORÓW
Obiekt: CWICZENIE PROJEKTOWE Z ROZPRZESTRZENIANIA ZANIECZYSZCZEN NR 20
Identyfikator obiektu GRAM
Wybrane emitory: od: 1 do: 1
------------------------------------------------------------------------
| Emitor | Współrzędne | Wysokość| Wymiar |
lp. | Nr | x [m], y [m] | h [m] |d[m], a[m]| Typ
------------------------------------------------------------------------
1| 1 | KOTLOWNIA |
| | 0,0 0,0 | 30,0 | 1,00 | OTWARTY
------------------------------------------------------------------------
Charakterystyka otoczenia
Kotłownia zlokalizowana jest na obszarze miejskim o liczbie mieszkańców do 10 tysięcy. W odległości 100 m w kierunku zachodnim od emitora zlokalizowany jest XI kondygnacyjny budynek mieszkalny. Różę wiatrów przyjęto dla miejscowości Częstochowa. Dla danego terenu wg. Tabeli 4 z Rozporządzenia Dz.U z 2010r. Nr 16 poz. 87 przyjęto współczynnik aerodynamicznej szorstkości podłoża wynoszący: z0 = 1
Tło zanieczyszczeń
Tło zanieczyszczeń czyli aktualny stan jakości powietrza określony jako stężenie uśrednione dla roku. Dla danego tematu wynosi ono R = 60%Da, zaś dla pyłu R = 10% Dp. Dopuszczalny poziom substancji w powietrzu (dla 1 roku) Da przyjmowany jest zgodnie z Rozporządzeniem Ministra Środowiska z dnia 26 stycznia 2010 r. w sprawie wartości odniesienia dla niektórych substancji w powietrzu.
Da ≥ Sa + R
Dp ≥ Op +Rp
Tło zanieczyszczeń powietrza – zawartość gazów lub pyłów w zanieczyszczonym powietrzu atmosferycznym, określana jako wartość średnia odniesiona do roku. Jest wyznaczana i porównywana przede wszystkim dla tych zanieczyszczeń, dla których ustalono dopuszczalne wartości stężeń.
Przyjęte oznaczenie dla tła zanieczyszczeń: R
Do celów projektowych przyjęto tło na poziomie 60 % wartości odniesienia uśrednionej dla roku. Przyjęte wartości odniesienia i tła przestawia tabela 2.
Tabela 2. Wartości odniesienia dla tła zanieczyszczeń
----------------------------------------------------------------------------------
ATMOTERM Opole EK100W
----------------------------------------------------------------------------------
MAGDALENA GRABOWSKA
----------------------------------------------------------------------------------
WARTOŚCI ODNIESIENIA
Obiekt: CWICZENIE PROJEKTOWE Z ROZPRZESTRZENIANIA ZANIECZYSZCZEN NR 20
Identyfikator obiektu: GRAM
---------------------------------------------------------------------------------
Substancja | Numer CAS | D1[ug/m3]| Da[ug/m3]|R[ug/m3] | Dp* Rp*
---------------------------------------------------------------------------------
70 ditl. azotu | 10102-44-0 | 200,0000| 40,0000| 24,0000 | - -
72 ditl. siarki | 7446-09-5 | 350,0000| 20,0000| 12,0000 | - -
137 pył zaw. PM10 | | 280,0000| 40,0000| 24,0000 | 200,00 20,000
150 tlenek węgla | 630-08-0 |30000,0000| - | - | - -
* - [g/m2*rok]
Obliczenia wstępne
Wydajność cieplna źródła (kotła):
Qkmax= QN
$$Q_{ksr} = Q_{N} \bullet \frac{k}{100},\ MW$$
QN – nominalna wydajność cieplna kotła, MW
k – obciążenie kotła, %
Qkmax= 10 MW
$$Q_{ksr} = 10 \bullet \frac{75}{100} = 7,5\ MW$$
Wydajność cieplna kotłowni:
Okres letni:
QhI= Qksr
QhI= 7,5 MW
Okres zimowy:
QhII= 2 * Qksr
QhII= 2 * 7,5
QhII= 15 MW
Okres zimowy maksymalny:
QhIII= 2 * Qkmax
QhIII= 2 * 10
QhIII= 20 MW
Zużycie paliwa dla kotłowni w poszczególnych podokresach:
$$B_{\text{hi}} = \frac{Q_{\text{hi}}}{W_{d} \cdot \eta}\ \ ,\ kg/h$$
Qhi – wydajność cieplna kotła,
Wd – wydajność opałowa paliwa, kJ/kg
η – sprawność cieplna kotła, %
Okres letni:
$$B_{\text{hI}} = \frac{7,5 \bullet 10^{3} \bullet 3600}{21800 \cdot 0,82} = 1510,4\ \ ,\ kg/h$$
Okres zimowy:
$$B_{\text{hII}} = \frac{15 \bullet 10^{3} \bullet 3600\ }{21800 \cdot 0,82} = 3020,8\ \ ,\ kg/h$$
Okres zimowy maksymalny:
$$B_{\text{hIII}} = \frac{20 \bullet 10^{3} \bullet 3600\ }{21800 \cdot 0,82} = 4027,7\ \ ,\ kg/h$$
Unos zanieczyszczeń (godzinowy)
Unos jest to wartość zanieczyszczeń wytworzonych w źródle.
Ilość energii SN wprowadzonej w paliwie
SN= $\frac{Q_{N}}{\eta}$
SN= $\frac{10}{0,82}$= 12,2 MW
• wskaźnik unosu SO2 wynosi 17,0 kg/Mg%,
• wskaźnik unosu NO2 wynosi 4,0 kg/Mg%,
• wskaźnik unosu CO wynosi 5,0 kg/Mg%,
• wskaźnik unosu pyłu całkowitego wynosi 3,0 kg/Mg%.
Zawartość części palnych w pyle wynosi 15 %
Unos SO2
$$U_{SO2} = \frac{W_{\text{SO}_{2}} \bullet s \bullet B_{\text{hi}}}{1000},\ \ kg/h$$
$$U_{SO2} = \frac{17 \bullet s \bullet B_{\text{hi}}}{1000},\ \ kg/h$$
Bhi – zużycie paliwa, kg/h
s – zawartość siarki w paliwie, %
Okres letni:
$$U_{SO2I} = \frac{17 \bullet 1,9 \bullet 1510,4}{1000} = 48,8\ \ kg/h$$
Okres zimowy:
$$U_{SO2II} = \frac{17 \bullet 1,9 \bullet 3020,8}{1000} = 97,6\ \ kg/h$$
Okres zimowy maksymalny:
$$U_{SO2III} = \frac{17 \bullet 1,9 \bullet 4027,7}{1000} = 130,1\ \ kg/h$$
Unos NO2
UNO2 = Bhi • WNO2, kg/h
Bhi – zużycie paliwa, kg/h
WNO2 – wskaźnik unosu NO2, kg/kg
Okres letni:
UNO2I = 1510, 4 • 0, 004 = 6, 0 kg/h
Okres zimowy:
UNO2II = 3020, 8 • 0, 004 = 12, 1 kg/h
Okres zimowy maksymalny:
UNO2III = 4027, 7 • 0, 004 = 16, 1 kg/h
Unos CO
UCO = Bhi • WCO, kg/h
Bhi – zużycie paliwa, kg/h
WCO – wskaźnik unosu CO, kg/kg
Okres letni:
UCOI = 1510, 4 • 0, 005 = 7, 6 kg/h
Okres zimowy:
UCOII = 3020, 8 • 0, 005 = 15, 1 kg/h
Okres zimowy maksymalny:
UCOIII = 4027, 7 • 0, 005 = 20, 1 kg/h
Unos pyłu całkowitego ze spalania węgla kamiennego
$$U_{p} = B_{\text{hi}} \bullet W_{p} \bullet Ar \bullet \frac{100}{100 - P}\ \ kg/h$$
Bhi – zużycie paliwa, kg/h
Wp – wskaźnik unosu pyłu, kg/Mg%
Ar – zawartość popiołu w paliwie, %
P – zawartość części palnych w pyle, % (P=15)
Okres letni:
$$U_{\text{pI}} = \frac{1510,4 \bullet 3 \bullet 18 \bullet 10^{- 3}}{1 - 15/100} = 96,0\ \ kg/h$$
Okres zimowy:
$$U_{\text{pII}} = \frac{3020,8 \bullet 3 \bullet 18 \bullet 10^{- 3}}{1 - 15/100} = 191,9\ \ kg/h$$
Okres zimowy maksymalny:
$$U_{\text{pIII}} = \frac{4027,7 \bullet 3 \bullet 18 \bullet 10^{- 3}}{1 - 15/100} = 255,9\ \ kg/h$$
Unos poszczególnych frakcji pyłu
Upji = μuj • Upi
μuj – udział j-tej frakcji w pyle całkowitym
Obliczenia wartości unosu poszczególnych zanieczyszczeń dla podanych wariantów pracy
kotłów.
Okres letni:
Frakcja 0-10μm
$$UpI = \frac{15}{100} \bullet 96,0 = 14,4\ \ kg/h$$
Frakcja 10-20μm
$$UpI = \frac{20}{100} \bullet 96,0 = 19,2\ \ kg/h$$
Frakcja 20-40μm
$$UpI = \frac{16}{100} \bullet 96,0 = 15,4\ \ kg/h$$
Frakcja 40-60μm
$$UpI = \frac{19}{100} \bullet 96,0 = 18,2\ \ kg/h$$
Frakcja 60-100μm
$$UpI = \frac{12}{100} \bullet 96,0 = 11,5\ \ kg/h$$
Frakcja >100μm
$$UpI = \frac{18}{100} \bullet 96,0 = 17,3\ \ kg/h$$
Okres zimowy:
Frakcja 0-10μm
$$UpII = \frac{15}{100} \bullet 191,9 = 28,8\ \ kg/h$$
Frakcja 10-20μm
$$UpII = \frac{20}{100} \bullet 191,9 = 38,4\ \ kg/h$$
Frakcja 20-40μm
$$UpII = \frac{16}{100} \bullet 191,9 = 30,7\ \ kg/h$$
Frakcja 40-60μm
$$UpII = \frac{19}{100} \bullet 191,9 = 36,5\ \ kg/h$$
Frakcja 60-100μm
$$UpII = \frac{12}{100} \bullet 191,9 = 23,0\ \ kg/h$$
Frakcja >100μm
$$UpII = \frac{18}{100} \bullet 191,9 = 34,5\ \ kg/h$$
Okres zimowy maksymalny:
Frakcja 0-10μm
$$UpIII = \frac{15}{100} \bullet 255,9 = 38,4\ \ kg/h$$
Frakcja 10-20μm
$$UpIII = \frac{20}{100} \bullet 255,9 = 51,2\ \ kg/h$$
Frakcja 20-40μm
$$UpIII = \frac{16}{100} \bullet 255,9 = 40,9\ \ kg/h$$
Frakcja 40-60μm
$$UpIII = \frac{19}{100} \bullet 255,9 = 48,6\ \ kg/h$$
Frakcja 60-100μm
$$UpIII = \frac{12}{100} \bullet 255,9 = 30,7\ \ kg/h$$
Frakcja >100μm
$$UpIII = \frac{18}{100} \bullet 255,9 = 46,1\ \ kg/h$$
Emisja substancji gazowych:
Obliczenia wartości emisji poszczególnych zanieczyszczeń dla podanych wariantów pracy
kotłów.
Emisje substancji gazowych:
Emisja SO2
$$E_{SO2} = U_{SO2} \bullet \left( 1 - \frac{\eta_{SO2}}{100} \right)\ \ \ kg/h$$
ηSO2 -skuteczność odsiarczania, %
Okres letni:
$$E_{SO2I} = 48,8 \bullet \left( 1 - \frac{30}{100} \right) = 34,2\ \ \ kg/h$$
Okres zimowy:
$$E_{SO2II} = 97,6 \bullet \left( 1 - \frac{30}{100} \right) = 68,3\ \ \ kg/h$$
Okres zimowy maksymalny:
$$E_{SO2III} = 130,1 \bullet \left( 1 - \frac{30}{100} \right) = 91,1\ \ \ kg/h$$
Emisja NO2
ENO2 = UNO2
Okres letni:
ENO2I = 6, 0 kg/h
Okres zimowy:
ENO2II = 12, 1 kg/h
Okres zimowy maksymalny:
ENO2III = 16, 1 kg/h
Emisja CO
ECO = UCO
Okres letni:
ECOI = 7, 6 kg/h
Okres zimowy:
ECOII = 15, 1 kg/h
Okres zimowy maksymalny:
$$\text{\ E}_{\text{COIII}} = 20,1\ \frac{\text{kg}}{h}$$
Obliczenie emisji pyłu całkowitego
$$E_{p} = U_{p} \bullet \left( 1 - \frac{\eta_{\text{pc}}}{100} \right),\ \ kg/h$$
Skuteczność całkowita odpylania
ηpc - całkowita skuteczność odpylania
ηuj- frakcyjna skuteczność odpylania
ηpc = (0,15 ∙ 0,03) + ( 0,2 ∙ 0,06) +( 0,16∙ 0,25) +(0,19 ∙ 0,52) +(0,12 ∙ 0,97) +(0,18 ∙ 0,97) = 0,4463
Obliczenie emisji pyłu całkowitego
Okres letni:
$$E_{\text{pI}} = 96,0 \bullet \left( 1 - \frac{44,63}{100} \right) = 53,13\ kg/h$$
Okres zimowy:
$$E_{\text{pII}} = 191,9 \bullet \left( 1 - \frac{44,63}{100} \right) = 106,26\ kg/h$$
Okres zimowy maksymalny:
$$E_{\text{pIII}} = 255,9 \bullet \left( 1 - \frac{44,63}{100} \right) = 141,68\ kg/h$$
Obliczenie emisji frakcyjnej pyłu
$E_{\text{pi}} = U_{\text{pi}} \bullet \left( 1 - \frac{\eta_{i}}{100} \right),\ \ kg/h$ ,
ηi – skuteczność odpylania
Okres letni:
Frakcja pyłu 0 – 10 μm :
$$E_{\text{pI}} = 14,4 \bullet \left( 1 - \frac{3}{100} \right) = 13,9615\ kg/h$$
Frakcja pyłu 10 – 20 μm :
$$E_{\text{pI}} = 19,2 \bullet \left( 1 - \frac{6}{100} \right) = 18,0396\ kg/h$$
Frakcja pyłu 20 – 40 μm :
$$E_{\text{pI}} = 15,4 \bullet \left( 1 - \frac{25}{100} \right) = 11,5146\ kg/h$$
Frakcja pyłu 40 – 60 μm :
$$E_{\text{pI}} = 18,2 \bullet \left( 1 - \frac{52}{100} \right) = 8,7511\ kg/h$$
Frakcja pyłu 60 – 100 μm :
$$E_{\text{pI}} = 11,5 \bullet \left( 1 - \frac{97}{100} \right) = 0,3454\ kg/h$$
Frakcja pyłu >100 μm :
$$E_{\text{pI}} = 17,3 \bullet \left( 1 - \frac{97}{100} \right) = 0,5182\ kg/h$$
Okres zimowy:
Frakcja pyłu 0 – 10 μm :
$$E_{\text{pII}} = 28,8 \bullet \left( 1 - \frac{3}{100} \right) = 27,9229\ kg/h$$
Frakcja pyłu 10 – 20 μm :
$$E_{\text{pII}} = 38,4 \bullet \left( 1 - \frac{6}{100} \right) = 36,0791\ kg/h$$
Frakcja pyłu 20 – 40 μm :
$$E_{\text{pII}} = 30,7 \bullet \left( 1 - \frac{25}{100} \right) = 23,0292\ kg/h$$
Frakcja pyłu 40 – 60 μm :
$$E_{\text{pII}} = 36,5 \bullet \left( 1 - \frac{52}{100} \right) = 17,5022\ kg/h$$
Frakcja pyłu 60 – 100 μm :
$$E_{\text{pI}} = 23,0 \bullet \left( 1 - \frac{97}{100} \right) = 0,6909\ kg/h$$
Frakcja pyłu >100 μm :
$$E_{\text{pI}} = 34,5 \bullet \left( 1 - \frac{97}{100} \right) = 1,0363\ kg/h$$
Okres zimowy maksymalny:
Frakcja pyłu 0 – 10 μm :
$$E_{\text{pIII}} = 38,4 \bullet \left( 1 - \frac{3}{100} \right) = 37,2306\ kg/h$$
Frakcja pyłu 10 – 20 μm :
$$E_{\text{pIII}} = 51,2 \bullet \left( 1 - \frac{6}{100} \right) = 48,1055\ kg/h$$
Frakcja pyłu 20 – 40 μm :
$$E_{\text{pIII}} = 40,9 \bullet \left( 1 - \frac{25}{100} \right) = 30,7056\ kg/h$$
Frakcja pyłu 40 – 60 μm :
$$E_{\text{pIII}} = 48,6 \bullet \left( 1 - \frac{52}{100} \right) = 23,3363\ kg/h$$
Frakcja pyłu 60 – 100 μm :
$$E_{\text{pI}} = 30,7 \bullet \left( 1 - \frac{97}{100} \right) = 0,9212\ kg/h$$
Frakcja pyłu >100 μm :
$$E_{\text{pI}} = 46,1 \bullet \left( 1 - \frac{97}{100} \right) = 1,3818\ kg/h$$
Frakcje powyżej 60 μm są usuwane w 97%.
Udział frakcji w pyle emitowanym
μEj= $\frac{E_{\text{pj}}}{E_{p}}$
Dla wszystkich okresów udział frakcji w emitowanym pyle jest taki sam. Uzyskane wyniki przedstawiono dla okresu letniego.
Dla d=5μm
$$\mu_{E} = \ \frac{13,9615}{53,13\ } = 0,2628$$
Dla d=15μm
$$\mu_{E} = \ \frac{18,0396\ }{53,13\ } = 0,3395$$
Dla d=30μm
$$\mu_{E} = \ \frac{11,5146}{53,13\ } = 0,2167$$
Dla d=50μm
$$\mu_{E} = \ \frac{8,7511\ }{53,13\ } = 0,1647$$
Dla d=80μm
$$\mu_{E} = \ \frac{0,3454\ \ }{53,13\ } = 0,0065$$
Dla d=150μm
$$\mu_{E} = \ \frac{0,5182\ }{53,13\ } = 0,0098$$
6.6. Prędkość swobodnego opadania ziaren pyłu
Dynamiczny współczynnik lepkości powietrza w warunkach rzeczywistych
gdzie:
-dynamiczny współczynnik lepkości powietrza dla warunków rzeczywistych;
C – stała Sutherlanda, dla powietrza C= 112
Tot - temperatura otoczenia; 280,7 K
$\eta_{\text{pow}} = 17,04\ \bullet 10^{- 6}\ \bullet \ \frac{273\ + 112\ }{280,7\ + \ 112}\ \bullet \ \left( \frac{280,7}{273} \right)^{1,5} = 1,74176\ \bullet 10^{- 5}\ Pa \bullet s$
Gęstość powietrza dla warunków rzeczywistych
ρpow= ρ0pow * $\frac{273}{T_{\text{ot}}}$
ρ0 pow - gęstość powietrza dla warunków normalnych dla powietrza suchego= 1,293 kg/m3
$$\rho_{\text{pow}} = 1,293 \bullet \ \frac{273}{280,7} = 1,2575\ \frac{\text{kg}}{m^{3}}$$
Do obliczeń przyjęto gęstość powietrza dla warunków rzeczywistych równą ρpow = 1,2571 kg/m3
Liczba Archimedesa
gdzie:
dp- średnica ziarna pyłu; m
g- przyspieszenie ziemskie; m/s2
gęstość pyłu; 1850 kg/m3
Wartości średnicy ziarna przyjęto jako wartości średnie przedziałów rozpatrywanych.
Zakres przedziału [μm] | Średnica ziarna pyłu [μm] | Średnica ziarna pyłu [m] |
---|---|---|
0 – 10 | 5 | 0,000005 |
10 – 20 | 15 | 0,000015 |
20 – 40 | 30 | 0,00003 |
40 – 60 | 50 | 0,00005 |
60 – 100 | 80 | 0,00008 |
>100 | 150 | 0,00015 |
Dla d=5μm
$$Ar = \frac{{0,000005}^{3} \bullet \ 1,2575\ \bullet 9,81 \bullet \ \left( 1850 - \ 1,2575 \right)}{\left( 1,74176\ \bullet 10^{- 5} \right)^{2}} = 0,0094\ $$
Dla d=15μm
$$Ar = \frac{{0,000015}^{3} \bullet \ 1,2575\ \bullet 9,81 \bullet \ \left( 1850 - \ 1,2575 \right)}{\left( 1,74176\ \bullet 10^{- 5} \right)^{2}} = 0,2537\ $$
Dla d=30μm
$$Ar = \frac{{0,00003}^{3} \bullet \ 1,2575\ \bullet 9,81 \bullet \ \left( 1850 - \ 1,2575 \right)}{\left( 1,74176\ \bullet 10^{- 5} \right)^{2}} = 2,0298\ $$
Dla d=50μm
$$Ar = \frac{{0,00005}^{3} \bullet \ 1,2575\ \bullet 9,81 \bullet \ \left( 1850 - \ 1,2575 \right)}{\left( 1,74176\ \bullet 10^{- 5} \right)^{2}} = 9,3971\ $$
Dla d=80μm
$$Ar = \frac{{0,00008}^{3} \bullet \ 1,2575\ \bullet 9,81 \bullet \ \left( 1850 - \ 1,2575 \right)}{\left( 1,74176\ \bullet 10^{- 5} \right)^{2}} = 38,4907\ $$
Dla d=150μm
$$Ar = \frac{{0,00015}^{3} \bullet \ 1,2575\ \bullet 9,81 \bullet \ \left( 1850 - \ 1,2575 \right)}{\left( 1,74176\ \bullet 10^{- 5} \right)^{2}} = 253,7230\ $$
Szacunkowa liczba Reynoldsa
ruch laminarny Repsz ≤ 2 => Rep= $\frac{\text{Ar}}{18}$
ruch przejściowy 2 < Repsz ≤ 500 => Rep= 0,152 * Ar0, 715
ruch turbulencyjny Repsz > 500 => Rep= 0,74 * Ar0, 5
Dla d=5μm
$$\text{Re}_{p_{\text{sz}}} = \ \frac{0,0094\ }{18 + 0,61 \bullet \ \sqrt{0,0094}} = 0,0005\ \ \ \ \ \Rightarrow \text{Re}_{p} = \ \frac{0,0094}{18} = \ 0,0005$$
Dla d=15μm
$$\text{Re}_{p_{\text{sz}}} = \ \frac{0,2537\ }{18 + 0,61 \bullet \ \sqrt{0,2537}} = 0,014\ \ \ \ \ \Rightarrow \text{Re}_{p} = \ \frac{0,2537}{18} = 0,0141$$
Dla d=30μm
$$\text{Re}_{p_{\text{sz}}} = \ \frac{2,0298\ }{18 + 0,61 \bullet \ \sqrt{2,0298}} = 0,1076\ \ \ \ \ \Rightarrow \text{Re}_{p} = \ \frac{2,0298}{18} = 0,1128$$
Dla d=50μm
$$\text{Re}_{p_{\text{sz}}} = \ \frac{9,3971\ }{18 + 0,61 \bullet \ \sqrt{9,3971}} = 0,4729\ \ \ \ \Rightarrow \text{Re}_{p} = \ \frac{9,3971}{18} = 0,5221\ $$
Dla d=80μm
$$\text{Re}_{p_{\text{sz}}} = \ \frac{38,4907\ }{18 + 0,61 \bullet \ \sqrt{38,4907}} = 1,7669\ \ \ \ \ \Rightarrow \text{Re}_{p} = \ \frac{38,4907}{18} = 2,1384$$
Dla d=150μm
$$\text{Re}_{p_{\text{sz}}} = \ \frac{253,7230\ }{18 + 0,61 \bullet \ \sqrt{253,7230}} = 9,1542\ \ \ \ \ \ \Rightarrow \text{Re}_{p} = 0,152 \bullet \ {275,4293}^{0,715} = 7,9609$$
Prędkość opadania ziarna pyłu
$$u_{p} = \frac{\text{Re}_{p}*\ \eta_{\text{pow}}}{d_{p} \bullet \rho_{\text{pow}}}$$
Dla d=5μm
$$u_{p} = \frac{0,0005 \bullet 1,74176\ \bullet 10^{- 5}}{0,000005 \bullet \ 1,2575} = 0,0014\ \frac{m}{s}$$
Dla d=15μm
$$u_{p} = \frac{0,0141 \bullet 1,74176 \bullet 10^{- 5}}{0,000015 \bullet \ 1,2575} = 0,0130\ \frac{m}{s}$$
Dla d=30μm
$$u_{p} = \frac{0,1128 \bullet 1,74176 \bullet 10^{- 5}}{0,00003 \bullet \ 1,2575} = 0,0521\ \frac{m}{s}$$
Dla d=50μm
$$u_{p} = \frac{0,5221 \bullet 1,74176 \bullet 10^{- 5}}{0,00005 \bullet \ 1,2575} = 0,1446\ \frac{m}{s}$$
Dla d=80μm
$$u_{p} = \frac{2,1384 \bullet 1,74176 \bullet 10^{- 5}}{0,00008 \bullet \ 1,2575} = 0,3702\ \frac{m}{s}$$
Dla d=150μm
$$u_{p} = \frac{7,9609 \bullet 1,74176 \bullet 10^{- 5}}{0,00015 \bullet \ 1,2575} = 0,7351\ \frac{m}{s}$$
Tabela 3. Zestawienie parametrów frakcji emitowanego pyłu
frakcja | Parametry frakcji pyłu oraz unosy i emisje frakcyjne | Prędkości opadania ziaren pyłu |
---|---|---|
∖nμuj |
∖nηuj |
|
μm | - | - |
0-10 | 0,15 | 0,03 |
10-20 | 0,2 | 0,06 |
20-40 | 0,16 | 0,25 |
40-60 | 0,19 | 0,52 |
60-100 | 0,12 | 0,97 |
>100 | 0,18 | 0,97 |
- | Σ=1 | - |
Określenie składu spalin suchych (przy braku analizy elementarnej paliwa)
Korzystamy z formuły na obliczanie współczynnika nadmiaru powietrza:
Λ=$\frac{21}{21 - O_{2}}$
O2= $\frac{21*(\lambda - 1)}{\lambda}$
O2= = $\frac{21*(1,5 - 1)}{1,5}$= 7%
Zawartość CO2
maksymalna zawartość CO2 w spalinach
CO2= 18,7% * (1 - $\frac{7}{21}$)= 12,5%
Zawartość N2
N2= 100 – 7 – 12,5 = 80,5 %
6.8. Określenie gęstości spalin suchych w warunkach normalnych
ρssn= $\frac{7}{100}$* $\frac{32}{22,4}$ + $\frac{80,5}{100}$ * $\frac{28}{22,4}$ + $\frac{18,7}{100}$ * $\frac{44}{22,4}$ = 1,35 kg/ m3
6.9. Określenie ilości powstających spalin przy spalaniu węgla kamiennego lub brunatnego (przy braku analizy elementarnej paliwa)
Minimalne jednostkowe zapotrzebowanie powietrza (dla λ=1)
vpmin= $\frac{1,012*\ W_{d}}{4,19*\ 10^{3}}$ + 0,5
vpmin= $\frac{1,012*\ 21800}{4,19*\ 10^{3}}$ + 0,5 =5,7653 m3/kg
Minimalna jednostkowa ilość powstających spalin wilgotnych
VSWmin= $\frac{0,89*\ W_{d}}{4,19*\ 10^{3}}$ + 1,65
VSWmin= $\frac{0,89*\ 21800}{4,19*\ 10^{3}}$ + 1,65 = 6,2805 m3/kg
Rzeczywista jednostkowa ilość spalin wilgotnych w warunkach normalnych
VSW=VSWmin + (λ−1) • vpmin
$$V_{\text{SW}}\mathbf{=}6,2805 + \left( 1,5 - 1 \right) \bullet 5,7653 = 9,1632\ \frac{m^{3}}{\text{kg}}$$
Λ- współczynnik nadmiaru powietrza
Strumień spalin wilgotnych w warunkach normalnych
$\dot{V}$SWn = VSw ·Bh, m3n /h
Okres letni:
$\dot{V}$SWI = 9,1632 · 1510,4 = 13840,1403 m3/h
Okres zimowy:
$\dot{V}$SWII = 9,1632 · 3020,8 = 27680,2806 m3/h
Okres zimowy maksymalny:
$\dot{V}$SWIII = 9,1632 · 4027,7 = 36907,0409 m3/h
Strumień spalin wilgotnych w warunkach rzeczywistych
$\dot{V}$SWrz= $\dot{V}$SWn· , m3n /h
Okres letni:
$\dot{V}$SWrz= 27680,2806·$\frac{425}{273}$ = 43092,0120 m3n /h
Okres zimowy:
$\dot{V}$SWrz= 13840,1403·$\frac{425}{273}$ = 21546,0060 m3n /h
Okres zimowy maksymalny:
$\dot{V}$SWrz= 36907,0409·$\frac{425}{273}$ = 57456,0160 m3n /h
Strumień spalin suchych w warunkach normalnych
VSSn= , m3n /h
x – zawartość wilgoci w spalinach; kg/kg suchych spalin
ρSSn – gęstość spalin suchych w warunkach normalnych, kg/m3n
Okres letni:
VSSn= $\frac{13840,1403}{1 + \ \frac{0,05*\ 1,35*22,4}{18}}$ = 12766,5230 m3/h
Okres zimowy:
VSSn= $\frac{27680,2806}{1 + \ \frac{0,05*\ 1,35*22,4}{18}}$ = 25533,0460 m3/h
Okres zimowy maksymalny:
VSSn= $\frac{36907,0409}{1 + \ \frac{0,05*\ 1,35*22,4}{18}}$ = 34044,0614 m3/h
Strumień spalin suchych w warunkach normalnych przy zawartości 6% O2 w spalinach.
, m3n /h
Okres letni:
$$\left( {\dot{V}}_{\text{SS}_{n}} \right)_{6\%} = 12766,5230 \bullet \frac{21 - 7}{21 - 6} = 11915,4215\ \frac{m^{3}}{h}$$
Okres zimowy:
$$\left( {\dot{V}}_{\text{SS}_{n}} \right)_{6\%} = 25533,0460 \bullet \frac{21 - 7}{21 - 6} = 23830,8430\ \frac{m^{3}}{h}$$
Okres zimowy maksymalny:
$$\left( {\dot{V}}_{\text{SS}_{n}} \right)_{6\%} = 34044,0614 \bullet \frac{21 - 7}{21 - 6} = 31774,4573\ \frac{m^{3}}{h}$$
6.10. Określenie stężeń zanieczyszczeń w spalinach w przeliczeniu na spaliny suche przy zawartości 6% O2
Dla wszystkich okresów stężenia zanieczyszczeń w spalinach w przeliczeniu na spaliny suche przy zawartości 6% O2 są takie same. Uzyskane wyniki przedstawiono dla okresu letniego.
Obliczenie stężenia SO2:
(SSO2 )SS6% = , mg/m3n
$$\left( S_{\text{SO}_{2}} \right)_{\text{SS}_{6\%}} = \frac{34,2\ \bullet 10^{6}}{11915,4215} = 2866,0553\ \frac{\text{mg}}{m^{3}}\ $$
Obliczenie stężenia NO2:
(SNO2 )SS6% = , mg/m3n
$$\left( S_{\text{NO}_{2}} \right)_{\text{SS}_{6\%}}\frac{6,0 \bullet 10^{6}}{11915,4215} = 507,0421\ \frac{\text{mg}}{m^{3}}$$
Obliczenie stężenia CO:
$$\left( S_{\text{CO}} \right)_{\text{SS}_{6\%}} = \frac{E_{\text{CO}}}{\left( {\dot{V}}_{\text{SS}_{n}} \right)_{6\%}},\ \frac{\text{mg}}{m^{3}}$$
$$\left( S_{\text{CO}} \right)_{\text{SS}_{6\%}} = \frac{7,6 \bullet 10^{6}}{11915,4215} = 633,8026\ \frac{\text{mg}}{m^{3}}$$
Obliczenie stężenia pyłu całkowitego:
(Sp)SS6% = , mg/m3n
Obliczenie stężenia frakcji pyłu PM10:
$$\left( S_{PM10} \right)_{\text{SS}_{6\%}} = \frac{13,9615 \bullet 10^{6}}{11915,4215} = 1171,7146\ \frac{\text{mg}}{m^{3}}$$
6.11. Prędkości wylotowe spalin
V= , m/s
- pole przekroju wylotowego emitora; m2
de2 - wewnętrzna średnica przekroju wylotowego emitora;
Okres letni:
$$v = \frac{4 \bullet \frac{21546,0060}{3600}}{\pi \bullet 1^{2}} = 7,6203\ \frac{m}{s}$$
Okres zimowy:
$$v = \frac{4 \bullet \frac{43092,0120}{3600}}{\pi \bullet {0,7}^{2}} = 15,2407\ \frac{m}{s}$$
Okres zimowy maksymalny:
$$v = \frac{4 \bullet \frac{57456,0160}{3600}}{\pi \bullet {0,7}^{2}} = 20,3209\ \frac{m}{s}$$
7. Parametry emisyjne przyjęte do obliczeń
W projekcie przyjęto 3 warianty pracy kotłowni wynikające z różnej pracy kotłów. Emisje w wariantach oraz roczny ładunek zanieczyszczeń na emitorze, przyjęte do obliczeń rozprzestrzeniania się zanieczyszczeń, zobrazowane są w poniższych tabelach.
Poniższa tabela zawiera dane o emisji w poszczególnych wariantach.
Tabela 4. Emisja wariantowa dla poszczególnych substancji
-------------------------------------------------------------------------------------
ATMOTERM Opole EK100W
-------------------------------------------------------------------------------------
MAGDALENA GRABOWSKA
-------------------------------------------------------------------------------------
EMISJA W WARIANTACH
Obiekt: CWICZENIE PROJEKTOWE Z ROZPRZESTRZENIANIA ZANIECZYSZCZEN NR 20
Identyfikator obiektu: GRAM
Wybrane emitory: od: 1 do: 1
-------------------------------------------------------------------------------------
Emitor |War.| Czas trwania [h] | Substancja |
Nr | Nr | Zima Lato Rok | kod nazwa CAS | Emisja [kg/h]
-------------------------------------------------------------------------------------
1 | 1 | 0,0 0,0 3790,0 Vwyl[m/s]= 7,6, Tsp[K]= 425,0
| | | 70 ditl. azotu , 10102-44-0| 6,0000000
| | | 72 ditl. siarki , 7446-09-5 | 34,2000000
| | |137 pył zaw. PM10, | 13,9615000
| | |150 tlenek węgla, 630-08-0 | 7,6000000
1 | 2 | 0,0 0,0 3280,0 Vwyl[m/s]=15,2, Tsp[K]= 425,0
| | | 70 ditl. azotu , 10102-44-0| 12,1000000
| | | 72 ditl. siarki , 7446-09-5 | 68,3000000
| | |137 pył zaw. PM10, | 27,9229000
| | |150 tlenek węgla, 630-08-0 | 27,9229000
1 | 3 | 0,0 0,0 1300,0 Vwyl[m/s]=20,3, Tsp[K]= 425,0
| | | 70 ditl. azotu , 10102-44-0| 16,1000000
| | | 72 ditl. siarki , 7446-09-5 | 91,1000000
| | |137 pył zaw. PM10, | 37,2306000
| | |150 tlenek węgla, 630-08-0 | 20,1000000
Tabela 5. Emisja średnia dla poszczególnych substancji
---------------------------------------------------------------------------------------
ATMOTERM Opole EK100W
---------------------------------------------------------------------------------------
MAGDALENA GRABOWSKA
---------------------------------------------------------------------------------------
EMISJA ŚREDNIA
Obiekt: CWICZENIE PROJEKTOWE Z ROZPRZESTRZENIANIA ZANIECZYSZCZEN NR 20
Identyfikator obiektu: GRAM
W - obliczona z wariantów, P - podana przez użytkownika
--------------------------------------------------------------------------------------
Emitor |Substancja, numer CAS | EmSr[kg/h] | EmSr[Mg/rok] |Tsp[K]|Vwyl[m/s]
--------------------------------------------------------------------------------------
1W| | | | 425,0| 12,58
| 72 ditl. siarki , 7446-09-5 | 56,40048 | 472,0720 | |
| 70 ditl. azotu , 10102-44-0| 9,959140 | 83,35800 | |
|150 tlenek węgla , 630-08-0 | 17,50551 | 146,5211 | |
|137 pył zaw. PM10, | 23,04671 | 192,9010 | |
Tabela 6. Pył całkowity
-------------------------------------------------------------------------------------
ATMOTERM Opole EK100W
-------------------------------------------------------------------------------------
MAGDALENA GRABOWSKA
-------------------------------------------------------------------------------------
ŁADUNEK SUBSTANCJI NA POSZCZEGÓLNYCH EMITORACH
( rok )
Obiekt: CWICZENIE PROJEKTOWE Z ROZPRZESTRZENIANIA ZANIECZYSZCZEN NR 20
Identyfikator obiektu: GRAM
------------------------------------------------------------------------------------
Nr | Substancja | Ładunek [Mg]
------------------------------------------
emitora | Kod i nazwa CAS |Gaz, pył zawieszony | Pył całkowity
------------------------------------------------------------------------------------
1| 70 ditl. azotu 10102-44-0| 83,358000 |
| 72 ditl. siarki 7446-09-5 | 472,072000 |
| 137 pył zaw. PM10 | 192,90100 | 734,079500
| 150 tlenek węgla 630-08-0 | 146,521100 |
8. Sprawdzenie warunków skróconego zakresu obliczeń
Jeżeli z obliczeń wstępnych wynika, ze spełnione są następujące warunki:
dla pojedynczego emitora lub zespołu emitorów, z których został utworzony emitor zastępczy:
oraz:
kryterium opadu pyłu ,
to na tym kończą się wymagane dla tego zakresu obliczenia.
Jeżeli nie jest spełniony warunek określony w punkcie a) i b) zakresu skróconego, to na całym obszarze to należy wykonać obliczenia opadu substancji pyłowych w sieci obliczeniowej rozkład maksymalnych stężeń substancji w powietrzu uśrednionych dla 1 godziny, z uwzględnieniem statystyki warunków meteorologicznych w celu sprawdzenia czy został spełniony warunek:
Smm ≤ D1
W przypadku naszego projektu obliczanego w programie komputerowym, zostały przez niego wyznaczone na określonym obszarze terenu stężenia maksymalne (w oparciu o emisje maksymalne) w osi wiatru dla poszczególnych sytuacji meteorologicznych. W kolejnym kroku sprawdza on, czy w żadnej sytuacji nie zostały przekroczone wartości stężeń odpowiadające 0, 1Da co oznacza, że konieczne jest spełnienie nierówności Smm ≤ 0, 1 • Da. Jeżeli tak nie jest, to wówczas warunek zakresu skróconego obliczeń nie jest spełniony.
Wyniki obliczeń z zakresu skróconego przedstawia tabela nr 7.
Tabela 7. Wyniki obliczeń
-------------------------------------------------------------------------
ATMOTERM Opole EK100W
-------------------------------------------------------------------------
MAGDALENA GRABOWSKA
-------------------------------------------------------------------------
ZAKRES OBLICZEŃ
Obliczenia dla wariantów emisji
Obiekt: CWICZENIE PROJEKTOWE Z ROZPRZESTRZENIANIA ZANIECZYSZCZEN NR 20
Identyfikator obiektu: GRAM
Wysokość anemometru: 14,0 Wektor szorstkości: 1,00000
Obszar: Obszar zwykły sezon: ROK
------------------------------------------------------------------
Substancja | Nr CAS | Smm[ug/m3] | 0,1*D1 | Zakres
------------------------------------------------------------------
70 ditl. azotu 10102-44-0 122,84077 20,00000 pełny
72 ditl. siarki 7446-09-5 693,39050 35,00000 pełny
137 zaw. PM10 141,73845 28,00000 pełny
150 tlenek węgla 630-08-0 283,47690 3000,00000 skrócony
Zakres skrócony oznacza, że substancja nie powoduje przekroczeń
10% dopuszczalnego poziomu w powietrzu lub 10% wartości odniesienia
------------------------------------------------------------------
Ponieważ nie są dotrzymane warunki zakresu skróconego należy wykonać pełne obliczenia w odniesieniu do trzech substancji: SO2, NO2 i frakcji pyłu PM10.
9. Obliczenia rozkładów przestrzennych stężeń maksymalnych i średniorocznych:
Do przeprowadzenia obliczeń rozkładów stężeń substancji należy zaprojektować siatkę obliczeniową. Projektując siatkę obliczeniową należało założyć zakres; założono, że emitor będzie znajdował się w środku układu (0,0). Parametry siatki:
- promień siatki:
R= 50 ∙ he = 50 ∙ 30 = 1500 m
gdzie: he- wysokość emitora
- długość kroku :
k= 3 ∙ he = 3 ∙ 30 = 90 m
- ilość kroków obliczeniowych:
n= $\frac{R}{k}$ = $\frac{1500}{90}$ = 17
Siatka obliczeniowa będzie miała wymiary 1500 m x 1500 m (w punkcie o współrzędnych (0,0) będzie zlokalizowany emitor). Prawy górny róg siatki będzie miał współrzędne (1500, 1500) a lewy dolny róg (-1500, -1500).
Do obliczeń rozkładów stężeń średniorocznych substancji wymagane jest określenie podokresów emisji (wariantów). W analizowanym przypadku są one równoważne wariantom pracy emitora.
Tabela 8. Podokresy emisji dla poszczególnych substancji
-----------------------------------------------------------------
ATMOTERM Opole EK100W
-----------------------------------------------------------------
MAGDALENA GRABOWSKA
-----------------------------------------------------------------
PODOKRESY NA SUBSTANCJĘ
Obiekt: CWICZENIE PROJEKTOWE Z ROZPRZESTRZENIANIA ZANIECZYSZCZEN NR 20
Identyfikator obiektu GRAM
Substancja: 70 ditl. azotu , 10102-44-0
----------------------------------------------
Podokres | Czas przyjęty | Emitor |Wariant
Nr |do obliczeń [h] | Nr | Nr
----------------------------------------------
R01 | 3790,0 | 1| 1
R02 | 3280,0 | 1| 2
R03 | 1300,0 | 1| 3
Substancja: 72 ditl. siarki , 7446-09-5
----------------------------------------------
Podokres | Czas przyjęty | Emitor |Wariant
Nr |do obliczeń [h] | Nr | Nr
----------------------------------------------
R01 | 3790,0 | 1| 1
R02 | 3280,0 | 1| 2
R03 | 1300,0 | 1| 3
Substancja: 137 pył zaw. PM10,
----------------------------------------------
Podokres | Czas przyjęty | Emitor |Wariant
Nr |do obliczeń [h] | Nr | Nr
----------------------------------------------
R01 | 3790,0 | 1| 1
R02 | 3280,0 | 1| 2
R03 | 1300,0 | 1| 3
Substancja: 150 tlenek węgla, 630-08-0
----------------------------------------------
Podokres | Czas przyjęty | Emitor |Wariant
Nr |do obliczeń [h] | Nr | Nr
----------------------------------------------
R01 | 3790,0 | 1| 1
R02 | 3280,0 | 1| 2
R03 | 1300,0 | 1| 3
Wyniki obliczeń przedstawiono w Tabeli 9:
Tabela 9. Przekroczenia stężeń dopuszczalnych
--------------------------------------------------------------------------------------------------
ATMOTERM Opole EK100W
--------------------------------------------------------------------------------------------------
MAGDALENA GRABOWSKA
--------------------------------------------------------------------------------------------------
ANALIZA STĘŻEŃ UŚREDNIONYCH DLA 1 GODZINY
Punkty z przekroczeniami dopuszczalnych norm stężeń
dla wybranych substancji
Obiekt: CWICZENIE PROJEKTOWE Z ROZPRZESTRZENIANIA ZANIECZYSZCZEN NR 20
Identyfikator obiektu: GRAM Zbiór wyników: T01GRAM.DBF
--------------------------------------------------------------------------------------------------
Współrzędne St. maksymalne Percentyl Stężenie średnioroczne Częstość
Z[m] X[m] Y[m] [µg/m3] [µg/m3] [µg/m3] przekroczeń
--------------------------------------------------------------------------------------------------
Współczynnik szorstkości z0 = 1,00000
----------------------------------------------------------------------------------------------
70 ditl. azotu (gaz) D1=200,000 Obszar zwykły
CAS 10102-44-0 percentyl 99,800 0,2%
Nie ma przekroczeń
----------------------------------------------------------------------------------------------
72 ditl. siarki (gaz) D1=350,000 Obszar zwykły
CAS 7446-09-5 percentyl 99,726 0,274%
0,0 210,0 -600,0 500,74268 (3) 301,18451 8,08693 0,17
0,0 -240,0 -510,0 531,17548 (3) 310,72934 8,62912 0,19
0,0 -150,0 -510,0 540,28265 (3) 324,29880 8,68089 0,22
0,0 -60,0 -510,0 543,53412 (3) 331,38336 8,38368 0,22
0,0 30,0 -510,0 544,07794 (3) 333,04443 8,60515 0,22
0,0 120,0 -510,0 542,44812 (3) 339,98669 9,53958 0,24
0,0 210,0 -510,0 534,90674 (3) 345,40686 10,00285 0,26
0,0 300,0 -510,0 520,13715 (3) 323,66956 9,53869 0,22
0,0 390,0 -510,0 497,74725 (3) 298,18765 8,27717 0,18
0,0 -330,0 -420,0 539,74268 (3) 333,25449 9,15590 0,23
0,0 -240,0 -420,0 546,25854 (3) 369,04056 10,54732 0,41
0,0 -150,0 -420,0 560,44049 (2) 392,65149 11,20432 0,50
0,0 -60,0 -420,0 568,91040 (2) 402,91248 10,83165 0,45
0,0 30,0 -420,0 570,61969 (2) 405,74277 11,05792 0,45
0,0 120,0 -420,0 564,37732 (2) 401,30551 12,41409 0,56
0,0 210,0 -420,0 549,34302 (2) 387,96780 12,67825 0,51
0,0 300,0 -420,0 543,53418 (3) 364,20728 11,20748 0,32
0,0 390,0 -420,0 527,47021 (3) 327,97308 9,23568 0,23
0,0 -420,0 -330,0 539,74268 (3) 315,04623 8,47986 0,20
0,0 -330,0 -330,0 550,99359 (2) 373,05261 10,50905 0,38
0,0 -240,0 -330,0 573,47998 (2) 425,26791 12,84070 0,65
0,0 -150,0 -330,0 620,34503 (3) 467,44232 14,57558 0,82
0,0 -60,0 -330,0 644,37183 (3) 469,99548 14,38644 0,81
0,0 30,0 -330,0 646,95447 (3) 467,18399 14,87070 0,78
0,0 120,0 -330,0 631,61237 (3) 478,62741 16,80316 0,95
0,0 210,0 -330,0 590,09052 (3) 460,32022 15,92348 0,81
0,0 300,0 -330,0 560,44049 (2) 397,70996 12,75867 0,56
0,0 390,0 -330,0 544,07794 (3) 338,22772 9,73634 0,24
0,0 -420,0 -240,0 546,25854 (3) 352,38327 9,44046 0,33
0,0 -330,0 -240,0 573,47998 (2) 417,07092 11,72371 0,52
0,0 -240,0 -240,0 641,15802 (3) 481,18716 15,16281 0,84
0,0 -150,0 -240,0 657,81415 (2) 524,92548 18,84096 1,38
0,0 -60,0 -240,0 681,24524 (2) 523,63525 19,35855 1,51
0,0 30,0 -240,0 681,24524 (2) 513,76959 19,98730 1,44
0,0 120,0 -240,0 671,77429 (2) 532,57269 22,52427 1,69
0,0 210,0 -240,0 654,11023 (3) 500,05560 18,52564 1,13
0,0 300,0 -240,0 597,81177 (3) 439,77197 13,76657 0,60
0,0 390,0 -240,0 555,41913 (2) 379,52573 10,92299 0,42
0,0 480,0 -240,0 538,66425 (3) 319,34000 9,39917 0,22
0,0 570,0 -240,0 508,81903 (3) 270,40543 8,40919 0,14
0,0 -510,0 -150,0 540,28265 (3) 326,98761 9,15988 0,22
0,0 -420,0 -150,0 560,44049 (2) 390,43295 10,92497 0,43
0,0 -330,0 -150,0 620,34503 (3) 460,47397 13,53952 0,68
0,0 -240,0 -150,0 657,81415 (2) 501,05673 17,06355 1,10
0,0 -150,0 -150,0 682,25330 (2) 544,36279 21,84374 1,89
0,0 -60,0 -150,0 666,40790 (3) 576,96112 22,86483 2,09
0,0 30,0 -150,0 666,40796 (3) 554,83112 22,81903 1,69
0,0 120,0 -150,0 693,25720 (2) 564,90851 26,12756 2,39
0,0 210,0 -150,0 678,52576 (2) 509,64749 20,33927 1,27
0,0 300,0 -150,0 644,37183 (3) 469,99548 17,09565 0,82
0,0 390,0 -150,0 571,76215 (2) 408,59293 14,35832 0,57
0,0 480,0 -150,0 545,16718 (3) 370,82233 12,59547 0,33
0,0 570,0 -150,0 521,17847 (3) 349,00717 11,45974 0,26
0,0 660,0 -150,0 481,10831 (3) 308,24854 10,03921 0,19
0,0 750,0 -150,0 437,06958 (3) 278,79272 8,47189 0,13
0,0 -600,0 -60,0 515,99261 (3) 298,65106 8,57014 0,17
0,0 -510,0 -60,0 543,53412 (3) 355,64120 10,35808 0,29
0,0 -420,0 -60,0 568,91040 (2) 414,84561 12,76546 0,59
0,0 -330,0 -60,0 644,37183 (3) 485,76230 16,18938 1,01
0,0 -240,0 -60,0 681,24524 (2) 534,70728 20,42136 1,67
0,0 -150,0 -60,0 666,40790 (3) 567,75818 22,06012 1,92
0,0 -60,0 -60,0 364,93338 (3) 298,03055 8,31886 0,09
0,0 120,0 -60,0 626,34534 (3) 516,10583 23,59980 1,36
0,0 210,0 -60,0 674,78973 (2) 534,34528 30,12868 2,23
0,0 300,0 -60,0 658,04675 (3) 519,42340 26,06394 1,57
0,0 390,0 -60,0 585,38867 (3) 481,98184 21,77233 1,13
0,0 480,0 -60,0 546,25854 (3) 442,34607 17,08383 0,72
0,0 570,0 -60,0 527,47021 (3) 373,86935 14,20831 0,36
0,0 660,0 -60,0 487,89120 (3) 342,31216 11,70341 0,23
0,0 750,0 -60,0 443,23160 (3) 297,29443 10,08366 0,14
0,0 840,0 -60,0 402,29721 (2) 258,97281 8,61858 0,06
0,0 -600,0 30,0 516,50885 (3) 300,14810 8,67336 0,17
0,0 -510,0 30,0 544,07794 (3) 357,78149 10,42379 0,29
0,0 -420,0 30,0 570,61969 (2) 422,38718 12,97780 0,62
0,0 -330,0 30,0 646,95447 (3) 497,33127 16,48208 1,06
0,0 -240,0 30,0 681,24524 (2) 534,27271 20,93391 1,82
0,0 -150,0 30,0 666,40796 (3) 576,89703 21,85216 1,87
0,0 120,0 30,0 599,53766 (1) 503,21411 33,24818 1,88
0,0 210,0 30,0 682,25330 (2) 569,24548 40,00212 3,29
0,0 300,0 30,0 659,36407 (3) 521,78540 31,95307 2,05
0,0 390,0 30,0 590,09052 (3) 485,86087 24,49850 1,18
0,0 480,0 30,0 546,25854 (3) 442,34604 19,13600 0,74
0,0 570,0 30,0 528,52625 (3) 375,74341 15,33107 0,36
0,0 660,0 30,0 488,86798 (3) 343,34067 12,53483 0,23
0,0 750,0 30,0 444,11899 (3) 297,88962 10,50277 0,14
0,0 840,0 30,0 402,69965 (2) 259,23190 8,95413 0,06
0,0 -510,0 120,0 542,44812 (3) 335,78418 9,38581 0,23
0,0 -420,0 120,0 564,37732 (2) 398,90344 11,35783 0,50
0,0 -330,0 120,0 631,61237 (3) 470,71667 13,90466 0,75
0,0 -240,0 120,0 671,77429 (2) 507,10565 17,22079 1,13
0,0 -150,0 120,0 693,25720 (2) 532,06995 21,25386 1,82
0,0 -60,0 120,0 626,34534 (3) 545,74884 25,49908 2,04
0,0 30,0 120,0 599,53766 (1) 527,26099 34,14540 2,34
0,0 120,0 120,0 664,07898 (2) 590,31946 42,84605 3,43
0,0 210,0 120,0 681,24524 (2) 534,27271 36,08693 2,59
0,0 300,0 120,0 651,49908 (3) 508,63116 28,70195 1,58
0,0 390,0 120,0 573,47998 (2) 476,14529 22,60263 1,03
0,0 480,0 120,0 546,25854 (3) 425,34750 18,11711 0,56
0,0 570,0 120,0 524,31488 (3) 369,98724 14,72967 0,35
0,0 660,0 120,0 484,00363 (3) 336,59152 12,22609 0,21
0,0 750,0 120,0 439,69992 (3) 293,74823 10,27965 0,14
0,0 840,0 120,0 401,89508 (2) 257,10089 8,82629 0,06
0,0 -510,0 210,0 534,90674 (3) 309,28287 8,16190 0,17
0,0 -420,0 210,0 549,34302 (2) 367,20578 9,69939 0,33
0,0 -330,0 210,0 590,09052 (3) 432,79163 12,01500 0,54
0,0 -240,0 210,0 654,11023 (3) 499,32452 16,39849 1,01
0,0 -150,0 210,0 678,52576 (2) 542,55511 22,78632 1,82
0,0 -60,0 210,0 674,78973 (2) 562,71594 31,44062 2,66
0,0 30,0 210,0 682,25330 (2) 574,65790 39,43855 3,28
0,0 120,0 210,0 681,24524 (2) 547,45264 41,53162 3,42
0,0 210,0 210,0 659,36407 (3) 524,92554 33,62520 2,05
0,0 300,0 210,0 616,63422 (3) 479,01074 25,26703 1,09
0,0 390,0 210,0 562,68677 (2) 442,26376 20,01850 0,80
0,0 480,0 210,0 542,44812 (3) 395,80893 16,27495 0,39
0,0 570,0 210,0 513,93274 (3) 344,99014 13,54638 0,27
0,0 660,0 210,0 472,99860 (3) 309,85129 11,41414 0,20
0,0 750,0 210,0 430,13211 (3) 277,02075 9,76350 0,13
0,0 840,0 210,0 400,69122 (2) 257,34711 8,45765 0,04
0,0 -420,0 300,0 543,53418 (3) 328,66245 8,96324 0,21
0,0 -330,0 300,0 560,44049 (2) 392,65149 11,63438 0,48
0,0 -240,0 300,0 597,81177 (3) 479,49002 16,01375 0,88
0,0 -150,0 300,0 644,37183 (3) 499,15781 20,86909 1,29
0,0 -60,0 300,0 658,04675 (3) 520,98590 26,70823 1,80
0,0 30,0 300,0 659,36407 (3) 526,22192 31,45337 2,12
0,0 120,0 300,0 651,49908 (3) 528,42798 32,77839 2,02
0,0 210,0 300,0 616,63422 (3) 503,66254 29,16404 1,40
0,0 300,0 300,0 568,91040 (2) 471,40756 23,28644 0,98
0,0 390,0 300,0 546,25854 (3) 424,65448 18,14531 0,51
0,0 480,0 300,0 529,58435 (3) 370,89035 14,70143 0,32
0,0 570,0 300,0 497,24979 (3) 329,01141 12,31429 0,25
0,0 660,0 300,0 457,18698 (3) 297,70184 10,53772 0,19
0,0 750,0 300,0 416,16943 (3) 272,79565 9,12432 0,06
0,0 -420,0 390,0 527,47021 (3) 303,36063 8,74183 0,19
0,0 -330,0 390,0 544,07794 (3) 366,76562 11,68194 0,33
0,0 -240,0 390,0 555,41913 (2) 414,85223 14,42370 0,66
0,0 -150,0 390,0 571,76215 (2) 463,08304 17,64088 0,86
0,0 -60,0 390,0 585,38867 (3) 496,66827 22,01685 1,22
0,0 30,0 390,0 590,09052 (3) 501,15848 24,29627 1,27
0,0 120,0 390,0 573,47998 (2) 478,23032 25,33582 1,25
0,0 210,0 390,0 562,68677 (2) 464,38925 23,84887 1,02
0,0 300,0 390,0 546,25854 (3) 439,97723 20,49575 0,62
0,0 390,0 390,0 534,90674 (3) 389,51657 16,92008 0,38
0,0 480,0 390,0 508,81903 (3) 343,06036 13,68424 0,26
0,0 570,0 390,0 474,41977 (3) 310,78223 11,31965 0,19
0,0 660,0 390,0 436,19635 (3) 281,48874 9,71412 0,12
0,0 750,0 390,0 401,89508 (2) 257,68115 8,48163 0,06
0,0 -420,0 480,0 500,24225 (3) 300,28229 8,92587 0,18
0,0 -330,0 480,0 524,31488 (3) 342,13013 11,02308 0,26
0,0 -240,0 480,0 538,66425 (3) 365,07123 12,94343 0,33
0,0 -150,0 480,0 545,16718 (3) 393,43130 15,20385 0,46
0,0 -60,0 480,0 546,25854 (3) 444,39911 17,76369 0,80
0,0 30,0 480,0 546,25854 (3) 445,73431 19,43140 0,82
0,0 120,0 480,0 546,25854 (3) 438,65930 19,91140 0,64
0,0 210,0 480,0 542,44812 (3) 418,33221 19,40400 0,47
0,0 300,0 480,0 529,58435 (3) 389,73557 17,71172 0,38
0,0 390,0 480,0 508,81903 (3) 362,29855 15,22654 0,28
0,0 480,0 480,0 480,14713 (3) 318,68115 12,95985 0,21
0,0 570,0 480,0 446,79169 (3) 289,77142 10,63789 0,14
0,0 660,0 480,0 412,02847 (3) 269,27231 9,07998 0,06
0,0 -420,0 570,0 465,49084 (3) 284,83395 8,82958 0,16
0,0 -330,0 570,0 489,84668 (3) 307,13150 10,03634 0,17
0,0 -240,0 570,0 508,81903 (3) 317,64410 11,55757 0,21
0,0 -150,0 570,0 521,17847 (3) 369,04980 13,39561 0,32
0,0 -60,0 570,0 527,47021 (3) 393,35913 14,96317 0,42
0,0 30,0 570,0 528,52625 (3) 394,54099 15,97810 0,42
0,0 120,0 570,0 524,31488 (3) 387,50275 16,24544 0,39
0,0 210,0 570,0 513,93274 (3) 365,01248 15,95884 0,32
0,0 300,0 570,0 497,24979 (3) 352,38455 15,10868 0,28
0,0 390,0 570,0 474,41977 (3) 327,57697 13,72195 0,22
0,0 480,0 570,0 446,79169 (3) 300,58273 12,00740 0,16
0,0 570,0 570,0 417,00259 (3) 277,53934 10,34110 0,06
0,0 660,0 570,0 400,69122 (2) 257,60458 8,56524 0,04
0,0 -420,0 660,0 428,84366 (3) 275,32788 8,24493 0,06
0,0 -330,0 660,0 450,38034 (3) 270,27115 9,17246 0,12
0,0 -240,0 660,0 468,29221 (3) 290,69540 10,26580 0,16
0,0 -150,0 660,0 481,10831 (3) 334,52884 11,72465 0,23
0,0 -60,0 660,0 487,89120 (3) 342,31216 12,68066 0,27
0,0 30,0 660,0 488,86798 (3) 343,34067 13,32902 0,27
0,0 120,0 660,0 484,00363 (3) 337,89090 13,58968 0,25
0,0 210,0 660,0 472,99860 (3) 326,59570 13,45183 0,23
0,0 300,0 660,0 457,18698 (3) 310,04672 12,92196 0,21
0,0 390,0 660,0 436,19635 (3) 290,82541 12,09989 0,14
0,0 480,0 660,0 412,02847 (3) 278,09497 10,99128 0,06
0,0 570,0 660,0 400,69122 (2) 281,45221 9,79954 0,04
0,0 660,0 660,0 395,12061 (2) 259,49130 8,51515 0,02
0,0 -330,0 750,0 410,79425 (3) 267,66150 8,31965 0,05
0,0 -240,0 750,0 425,85226 (3) 276,28302 9,23656 0,07
0,0 -150,0 750,0 437,06958 (3) 291,40762 10,19972 0,16
0,0 -60,0 750,0 443,23160 (3) 321,45905 11,09439 0,16
0,0 30,0 750,0 444,11899 (3) 321,78067 11,46658 0,16
0,0 120,0 750,0 439,69992 (3) 315,29724 11,69149 0,16
0,0 210,0 750,0 430,13211 (3) 285,35190 11,54199 0,15
0,0 300,0 750,0 416,16943 (3) 277,81702 11,19337 0,06
0,0 390,0 750,0 401,89508 (2) 280,32867 10,66918 0,06
0,0 480,0 750,0 399,09167 (2) 275,81302 9,94925 0,04
0,0 570,0 750,0 393,93704 (2) 257,86227 9,02237 0,02
0,0 660,0 750,0 386,90967 (2) 256,57620 8,01183 0,01
0,0 -240,0 840,0 400,29077 (2) 281,45221 8,39957 0,04
0,0 -150,0 840,0 401,89511 (2) 281,30112 9,17224 0,07
0,0 -60,0 840,0 402,29721 (2) 295,51602 9,72685 0,07
0,0 30,0 840,0 402,69965 (2) 296,10764 10,02136 0,07
0,0 120,0 840,0 401,89508 (2) 292,86829 10,14692 0,07
0,0 210,0 840,0 400,69122 (2) 281,58255 10,06475 0,04
0,0 300,0 840,0 399,09167 (2) 274,98682 9,80381 0,04
0,0 390,0 840,0 396,30777 (2) 262,36145 9,42476 0,04
0,0 480,0 840,0 391,18909 (2) 257,60455 8,95264 0,01
0,0 570,0 840,0 384,59509 (2) 255,68924 8,34899 0,01
0,0 -150,0 930,0 393,93704 (2) 273,71640 8,18103 0,02
0,0 -60,0 930,0 395,51599 (2) 277,29794 8,63277 0,02
0,0 30,0 930,0 395,51599 (2) 277,57538 8,87495 0,05
0,0 120,0 930,0 394,72571 (2) 275,08838 8,98420 0,02
0,0 210,0 930,0 393,14993 (2) 269,91107 8,89825 0,02
0,0 300,0 930,0 389,62750 (2) 257,34711 8,68619 0,01
0,0 390,0 930,0 385,75061 (2) 256,20108 8,39828 0,01
0,0 480,0 930,0 380,00754 (2) 249,81560 8,07678 0,01
----------------------------------------------------------------------------------------------
137 pył zaw. PM10(pył) D1=280,000 Obszar zwykły
CAS percentyl 99,800 0,2%
Nie ma przekroczeń
----------------------------------------------------------------------------------------------
150 tlenek węgla (gaz) D1=30000,0 Obszar zwykły
CAS 630-08-0 percentyl 99,800 0,2%
Nie ma przekroczeń
W żadnym punkcie stężenie nie przekracza
10% wartości odniesienia
Na podstawie wykonanych obliczeń można odczytać wartość maksymalnego stężenia substancji w punkcie. Porównując tą wartość z maksymalną dopuszczalną wartością można sprawdzić czy występują przekroczenia.
Dla PM10 oraz dwutlenku azotu nie odnotowano przekroczeń maksymalnych.
Maksymalna wartość stężenia dwutlenku siatki wynosi 590,31946 µg/m3 i przekracza dopuszczalną wartość 350 µg/m3 stanowiąc 168,66 % dopuszczalnej wartości. Przekroczenie dopuszczalnej wartości maksymalnej powoduje konieczność przeprowadzenia redukcji emisji. Wartość, do jakiej powinno się obniżyć stężenie powinna odpowiadać wartości emisji dopuszczalnej wyliczonej zgodnie z najnowszymi aktami prawnymi: Rozporządzeniem Ministra Środowiska w sprawie wartości odniesienia niektórych substancji w powietrzu oraz Rozporządzeniem Ministra Środowiska w sprawie standardów emisyjnych z instalacji.
10. Wymagane redukcje emisji :
Ze względu na występujące przekroczenia, zostały przeprowadzone obliczenia stopni redukcji emisji na podstawie najwyższych wartości stężeń.
Tabela 10. Analiza najwyższych wartości stężeń.
--------------------------------------------------------------------------------------------------
ATMOTERM Opole EK100W
--------------------------------------------------------------------------------------------------
MAGDALENA GRABOWSKA
--------------------------------------------------------------------------------------------------
ANALIZA STĘŻEŃ UŚREDNIONYCH DLA 1 GODZINY
Punkty z maksymalnymi wartościami
Obiekt: CWICZENIE PROJEKTOWE Z ROZPRZESTRZENIANIA ZANIECZYSZCZEN NR 20
Identyfikator obiektu: GRAM Zbiór wyników: T01GRAM.DBF
* - wartosc maksymalna
--------------------------------------------------------------------------------------------------
Współrzędne St. maksymalne Percentyl Stężenie średnioroczne Częstość
Z[m] X[m] Y[m] [µg/m3] [µg/m3] [µg/m3] przekroczeń
--------------------------------------------------------------------------------------------------
Współczynnik szorstkości z0 = 1,00000
----------------------------------------------------------------------------------------------
70 ditl. azotu (gaz) D1=200,000 Obszar zwykły
CAS 10102-44-0 percentyl 99,800 0,2%
0,0 -150,0 120,0 122,81718* 99,65326 3,75299 0,00
0,0 30,0 210,0 120,86772 107,24877* 6,96402 0,00
0,0 120,0 120,0 117,64796 106,54928 7,56571* 0,00
0,0 1470,0 1470,0 47,22558 30,91552 0,50214 0,00*
Wymagane obliczenia rozkładu stężeń uśrednionych dla roku, ponieważ
maksymalne stężenie 1-godz. przekracza 10% wartości odniesienia
i 10% dopuszczalnego poziomu substancji w powietrzu
----------------------------------------------------------------------------------------------
72 ditl. siarki (gaz) D1=350,000 Obszar zwykły
CAS 7446-09-5 percentyl 99,726 0,274%
0,0 -150,0 120,0 693,25720* 532,06995 21,25386 1,82
0,0 120,0 120,0 664,07898 590,31946* 42,84605* 3,43*
Wymagane obliczenia rozkładu stężeń uśrednionych dla roku, ponieważ
maksymalne stężenie 1-godz. przekracza 10% wartości odniesienia
i 10% dopuszczalnego poziomu substancji w powietrzu
----------------------------------------------------------------------------------------------
137 pył zaw. PM10(pył) D1=280,000 Obszar zwykły
CAS percentyl 99,800 0,2%
0,0 -150,0 120,0 141,71121* 114,98379 4,34244 0,00
0,0 30,0 210,0 139,46185 124,00422* 8,05781 0,00
0,0 120,0 120,0 135,74678 122,94113 8,75401* 0,00
0,0 1470,0 1470,0 54,60363 35,67166 0,58101 0,00*
Wymagane obliczenia rozkładu stężeń uśrednionych dla roku, ponieważ
maksymalne stężenie 1-godz. przekracza 10% wartości odniesienia.
----------------------------------------------------------------------------------------------
150 tlenek węgla (gaz) D1=30000,0 Obszar zwykły
CAS 630-08-0 percentyl 99,800 0,2%
0,0 -150,0 120,0 283,42242* 210,02635 6,59675 0,00
0,0 30,0 210,0 278,92371 234,59233* 12,24089 0,00
0,0 120,0 120,0 271,49356 230,49591 13,29850* 0,00
0,0 1470,0 1470,0 108,81500 52,63937 0,88263 0,00*
W żadnym punkcie stężenie nie przekracza
10% wartości odniesienia
Dla NO2 zanotowano, że stężenie maksymalne wynosi 122,81718 a maksymalna częstość przekroczeń wynosi 0,00. Dla SO2 i PM10 odczytujemy wartości w taki sam sposób, zaś dla CO nie zanotowano aby stężenie przekraczało 10% wartości odniesienia.
Stopień redukcji emisji :
st. Red. ze względu na stężenie max= $\frac{max.wartosc\ percentyla - D_{1}}{max.\ wartosc\ percentyla}$
st.red. ze względu na stęż. Średnioroczne= $\frac{max.\ wartosc\ stezenia\ sredniorocznego - W}{max.\ wartosc\ stezenia\ sredniorocznego}$
gdzie: W=Da∙ (1-tło)
Dla NO2:
st.red= $\frac{107,24877 - 200}{107,24877}$= -0,865
st.red= $\frac{7,56571 - 16}{7,56571}$= -1,115 , W= 40∙ (1-0,6) = 16
Dla SO2:
st.red= $\frac{590,31946 - 350}{590,31946}$= 0,407
st.red = $\frac{42,84605 - 8}{42,84605}$= 0,813 , W= 20∙ (1-0,6) = 8
PM10:
st.red= $\frac{124,00422 - 280}{124,00422}$= -1,258
st.red= $\frac{8,75401 - 16}{8,75401}$= -0,828, W= 40∙ (1-0,6) = 16
W=16 > W=8
Tabela 11. Emisja w wariantach po redukcji
-------------------------------------------------------------------------------------
ATMOTERM Opole EK100W
-------------------------------------------------------------------------------------
MAGDALENA GRABOWSKA
-------------------------------------------------------------------------------------
EMISJA W WARIANTACH
Obiekt: CWICZENIE PROJEKTOWE Z ROZPRZESTRZENIANIA ZANIECZYSZCZEN NR 20
Identyfikator obiektu: GRAM
Wybrane emitory: od: 1 do: 1
-------------------------------------------------------------------------------------
Emitor |War.| Czas trwania [h] | Substancja |
Nr | Nr | Zima Lato Rok | kod nazwa CAS | Emisja [kg/h]
-------------------------------------------------------------------------------------
1 | 1 | 0,0 0,0 3790,0 Vwyl[m/s]= 7,6, Tsp[K]= 425,0
| | | 70 ditl. azotu , 10102-44-0| 6,0000000
| | | 72 ditl. siarki , 7446-09-5 | 6,2000000
| | |137 pył zaw. PM10, | 13,9615000
| | |150 tlenek węgla, 630-08-0 | 7,6000000
1 | 2 | 0,0 0,0 3280,0 Vwyl[m/s]=15,2, Tsp[K]= 425,0
| | | 70 ditl. azotu , 10102-44-0| 12,1000000
| | | 72 ditl. siarki , 7446-09-5 | 12,4000000
| | |137 pył zaw. PM10, | 27,9229000
| | |150 tlenek węgla, 630-08-0 | 27,9229000
1 | 3 | 0,0 0,0 1300,0 Vwyl[m/s]=20,3, Tsp[K]= 425,0
| | | 70 ditl. azotu , 10102-44-0| 16,1000000
| | | 72 ditl. siarki , 7446-09-5 | 16,8000000
| | |137 pył zaw. PM10, | 37,2306000
| | |150 tlenek węgla, 630-08-0 | 20,1000000
Eredi= Ei∙ (1-st.red)
EredI= 34,2∙ (1-0,813)= 6,386
EredII= 68,3∙ (1-0,813)= 12,753
EredIII= 91,1∙ (1-0,813)= 17,010
Tabela12. Emisja średnia po redukcji
---------------------------------------------------------------------------------------
ATMOTERM Opole EK100W
---------------------------------------------------------------------------------------
MAGDALENA GRABOWSKA
---------------------------------------------------------------------------------------
EMISJA ŚREDNIA
Obiekt: CWICZENIE PROJEKTOWE Z ROZPRZESTRZENIANIA ZANIECZYSZCZEN NR 20
Identyfikator obiektu: GRAM
W - obliczona z wariantów, P - podana przez użytkownika
--------------------------------------------------------------------------------------
Emitor |Substancja, numer CAS | EmSr[kg/h] | EmSr[Mg/rok] |Tsp[K]|Vwyl[m/s]
--------------------------------------------------------------------------------------
1W| | | | 425,0| 12,58
| 72 ditl. siarki , 7446-09-5 | 10,27598 | 86,00999 | |
| 70 ditl. azotu , 10102-44-0| 9,959140 | 83,35800 | |
|150 tlenek węgla , 630-08-0 | 17,50551 | 146,5211 | |
|137 pył zaw. PM10, | 23,04671 | 192,9010 | |
Tabela13. Najwyższe wartości stężeń maksymalnych I średniorocznych po redukcji emisji
--------------------------------------------------------------------------------------------------
ATMOTERM Opole EK100W
--------------------------------------------------------------------------------------------------
MAGDALENA GRABOWSKA
--------------------------------------------------------------------------------------------------
ANALIZA STĘŻEŃ UŚREDNIONYCH DLA 1 GODZINY
Punkty z maksymalnymi wartościami
Obiekt: CWICZENIE PROJEKTOWE Z ROZPRZESTRZENIANIA ZANIECZYSZCZEN NR 20
Identyfikator obiektu: GRAM Zbiór wyników: T02GRAM.DBF
* - wartosc maksymalna
--------------------------------------------------------------------------------------------------
Współrzędne St. maksymalne Percentyl Stężenie średnioroczne Częstość
Z[m] X[m] Y[m] [µg/m3] [µg/m3] [µg/m3] przekroczeń
--------------------------------------------------------------------------------------------------
Współczynnik szorstkości z0 = 1,00000
----------------------------------------------------------------------------------------------
70 ditl. azotu (gaz) D1=200,000 Obszar zwykły
CAS 10102-44-0 percentyl 99,800 0,2%
0,0 -150,0 120,0 122,81718* 99,65326 3,75299 0,00
0,0 30,0 210,0 120,86772 107,24877* 6,96402 0,00
0,0 120,0 120,0 117,64796 106,54928 7,56571* 0,00
0,0 1470,0 1470,0 47,22558 30,91552 0,50214 0,00*
Wymagane obliczenia rozkładu stężeń uśrednionych dla roku, ponieważ
maksymalne stężenie 1-godz. przekracza 10% wartości odniesienia
i 10% dopuszczalnego poziomu substancji w powietrzu
----------------------------------------------------------------------------------------------
72 ditl. siarki (gaz) D1=350,000 Obszar zwykły
CAS 7446-09-5 percentyl 99,726 0,274%
0,0 -150,0 120,0 125,86221* 97,07301 3,87239 0,00
0,0 120,0 120,0 121,54967 108,86242* 7,80641* 0,00
0,0 1470,0 1470,0 49,27887 26,59475 0,51812 0,00*
Wymagane obliczenia rozkładu stężeń uśrednionych dla roku, ponieważ
maksymalne stężenie 1-godz. przekracza 10% wartości odniesienia
i 10% dopuszczalnego poziomu substancji w powietrzu
----------------------------------------------------------------------------------------------
137 pył zaw. PM10(pył) D1=280,000 Obszar zwykły
CAS percentyl 99,800 0,2%
0,0 -150,0 120,0 141,71121* 114,98379 4,34244 0,00
0,0 30,0 210,0 139,46185 124,00422* 8,05781 0,00
0,0 120,0 120,0 135,74678 122,94113 8,75401* 0,00
0,0 1470,0 1470,0 54,60363 35,67166 0,58101 0,00*
Wymagane obliczenia rozkładu stężeń uśrednionych dla roku, ponieważ
maksymalne stężenie 1-godz. przekracza 10% wartości odniesienia.
----------------------------------------------------------------------------------------------
150 tlenek węgla (gaz) D1=30000,0 Obszar zwykły
CAS 630-08-0 percentyl 99,800 0,2%
0,0 -150,0 120,0 283,42242* 210,02635 6,59675 0,00
0,0 30,0 210,0 278,92371 234,59233* 12,24089 0,00
0,0 120,0 120,0 271,49356 230,49591 13,29850* 0,00
0,0 1470,0 1470,0 108,81500 52,63937 0,88263 0,00*
W żadnym punkcie stężenie nie przekracza
10% wartości odniesienia
11. Obliczenie stężeń maksymalnych na wysokości zabudowy mieszkalnej:
Obliczenia stężeń uśrednionych dla godziny na poziomie zabudowy mieszkalnej wykonujemy wtedy, kiedy w odległości mniejszej niż 10∙he znajdują się wyższe niż parterowe budynki mieszkalne, biurowe lub użyteczności publicznej. Należy wtedy sprawdzić czy nie są one narażone na przekroczenia wartości odniesienia substancji w powietrzu. W tym celu należy obliczyć maksymalne stężenia substancji w powietrzu dla odpowiednich wysokości. Rozróżniamy dwa przypadki:
gdy geometryczna wysokość emitora (najniższego) jest większa od wysokości najwyżej położonej kondygnacji budynku Z w zadanej odległości od emitora -> obliczenia stężeń wykonuje się na wysokości Z
gdy geometryczna wysokość emitora (najniższego) jest mniejsza od wysokości najwyżej położonej kondygnacji budynku w zadanej odległości od emitora -> obliczenia wykonuje się co 1m zaczynając od geometrycznej wysokości najniższego emitora do wysokości:
gdy Hmax≥ Z (Hmax – najwyższa efektywna wysokość emitora w zespole) – obliczenia wykonujemy od wysokości Z
gdy Hmax< Z – obliczenia wykonujemy do wysokości Hmax
Wysokość budynku: Z= (11∙ 2,8)+1= 32m
Współrzędne położenia budynku: (-100;0)
W programie EKAW100 wykonano odpowiednie obliczenia co 1m.
Tabela 14. Stężenia maksymalne na wysokości budynku
---------------------------------------------------------------------------
ATMOTERM Opole EK100W
---------------------------------------------------------------------------
MAGDALENA GRABOWSKA
---------------------------------------------------------------------------
ANALIZA STĘŻEŃ UŚREDNIONYCH DLA 1 GODZINY
Punkty z przekroczeniami dopuszczalnych norm stężeń
dla wybranych substancji
Obiekt: CWICZENIE PROJEKTOWE Z ROZPRZESTRZENIANIA ZANIECZYSZCZEN NR 20
Identyfikator obiektu: GRAM Zbiór wyników: T03GRAM.DBF
---------------------------------------------------------------------------
Współrzędne St. maksymalne Percentyl
Z[m] X[m] Y[m] [µg/m3] [µg/m3]
---------------------------------------------------------------------------
Współczynnik szorstkości z0 = 1,00000
-----------------------------------------------------------------------
70 ditl. azotu (gaz) D1=200,000 Obszar zwykły
CAS 10102-44-0 percentyl 99,800
30,0 -100,0 0,0 574,83069 (2) 422,01788
31,0 -100,0 0,0 648,11981 (2) 489,22073
32,0 -100,0 0,0 719,87354 (2) 574,67950
-----------------------------------------------------------------------
72 ditl. siarki (gaz) D1=350,000 Obszar zwykły
CAS 7446-09-5 percentyl 99,726
30,0 -100,0 0,0 598,58398 (3) 424,36850
31,0 -100,0 0,0 674,90155 (3) 480,15457
32,0 -100,0 0,0 749,62030 (3) 541,09888
-----------------------------------------------------------------------
137 pył zaw. PM10(pył) D1=280,000 Obszar zwykły
CAS percentyl 99,800
30,0 -100,0 0,0 1326,52625 (3) 982,00043
31,0 -100,0 0,0 1495,65405 (3) 1128,96375
32,0 -100,0 0,0 1661,23889 (3) 1337,22632
-----------------------------------------------------------------------
150 tlenek węgla (gaz) D1=30000,0 Obszar zwykły
CAS 630-08-0 percentyl 99,800
Nie ma przekroczeń
W żadnym punkcie stężenie nie przekracza
10% wartości odniesienia
Po wykonaniu obliczeń w programie nie odnotowano przekroczeń na wysokości zabudowy dla tlenku węgla CO. Natomiast w przypadku SO2, NO2 i CO odnotowano przekroczenia.
12. Obliczenia opadu pyłu:
W celu sprawdzenia warunku (III) należy wykonać obliczenia w programie EKAW100 opadu pyłu w sieci obliczeniowej. Obliczenia wykonane zostały w tej samej siatce obliczeniowej, co obliczenia emisji substancji gazowych.
Parametry frakcji pyłu przyjęte do obliczeń przedstawiono w Tabeli 14.
Tabela 15. Parametry frakcji pyłu.
---------------------------------------------------------------------------
ATMOTERM Opole EK100W
---------------------------------------------------------------------------
MAGDALENA GRABOWSKA
---------------------------------------------------------------------------
PARAMETRY FRAKCJI PYŁU
Obiekt: CWICZENIE PROJEKTOWE Z ROZPRZESTRZENIANIA ZANIECZYSZCZEN NR 20
Identyfikator obiektu: GRAM
Substancja: pył ogółem
--------------------------------------------------------------------
Nr | Nr | | Udział | Emisja |
emitora |wariantu| Frakcje[µm] | proc.[%] | [Mg/rok] | Wf[m/s]
--------------------------------------------------------------------
1 | 1 | 0-10 | 26,28 | 52,91812 | 0,001
| | 10-20 | 33,95 | 68,36264 | 0,013
| | 20-40 | 21,67 | 43,63530 | 0,052
| | 40-60 | 16,47 | 33,16444 | 0,145
| | 60-100 | 0,65 | 1,30886 | 0,370
| | 100-200 | 0,98 | 1,97335 | 0,735
| 2 | 0-10 | 26,28 | 91,59442 | 0,001
| | 10-20 | 33,95 | 118,32690 | 0,013
| | 20-40 | 21,67 | 75,52706 | 0,052
| | 40-60 | 16,47 | 57,40335 | 0,145
| | 60-100 | 0,65 | 2,26546 | 0,370
| | 100-200 | 0,98 | 3,41562 | 0,735
| 3 | 0-10 | 26,28 | 48,40355 | 0,001
| | 10-20 | 33,95 | 62,53047 | 0,013
| | 20-40 | 21,67 | 39,91267 | 0,052
| | 40-60 | 16,47 | 30,33510 | 0,145
| | 60-100 | 0,65 | 1,19720 | 0,370
| | 100-200 | 0,98 | 1,80500 | 0,735
Wyniki obliczeń przedstawiono w Tabeli 15:
Tabela 16. Przekroczenia
---------------------------------------------------------------------------
ATMOTERM Opole EK100W
---------------------------------------------------------------------------
MAGDALENA GRABOWSKA
---------------------------------------------------------------------------
ANALIZA OPADU PYŁU
Punkty z wartościami opadu przekraczającymi progi
Obiekt: CWICZENIE PROJEKTOWE Z ROZPRZESTRZENIANIA ZANIECZYSZCZEN NR 20
Identyfikator obiektu: GRAM Zbiór wyników: D01GRAM.DBF
---------------------------------------------------------------------------
Współrzędne Pył ogółem
X[m] Y[m] [g/m2*rok]
---------------------------------------------------------------------------
Współczynnik szorstkości z0 = 1,00000
--------------------------------------------------------------------------
Obszar zwykły Dp-Rp=180,000
-150,0 -240,0 218,81821
-60,0 -240,0 222,17070
30,0 -240,0 231,58420
120,0 -240,0 299,93674
210,0 -240,0 211,77307
-150,0 -150,0 371,68253
-60,0 -150,0 547,85822
30,0 -150,0 456,30417
120,0 -150,0 532,70709
210,0 -150,0 228,19080
-240,0 -60,0 294,09277
-150,0 -60,0 429,93646
-60,0 -60,0 557,06714
30,0 -60,0 359,34723
120,0 -60,0 563,05316
210,0 -60,0 308,27448
300,0 -60,0 362,60114
390,0 -60,0 210,06877
-240,0 30,0 306,35226
-150,0 30,0 602,63623
-60,0 30,0 201,52541
120,0 30,0 1244,10327
210,0 30,0 714,01727
300,0 30,0 373,61465
390,0 30,0 214,19354
-150,0 120,0 301,57285
-60,0 120,0 691,57233
30,0 120,0 986,95117
120,0 120,0 1030,96875
210,0 120,0 493,66586
300,0 120,0 282,07962
-150,0 210,0 292,27829
-60,0 210,0 386,16031
30,0 210,0 561,92401
120,0 210,0 605,02264
210,0 210,0 411,79184
300,0 210,0 216,87282
-150,0 300,0 181,41264
-60,0 300,0 296,45752
30,0 300,0 304,78082
120,0 300,0 347,92023
210,0 300,0 268,74994
300,0 300,0 196,58858
-60,0 390,0 180,47437
30,0 390,0 183,64505
120,0 390,0 213,73524
Z tabel wynika, że kryterium opadu pyłu nie zostało spełnione.
13. Wymagane redukcje emisji:
Ponieważ odnotowano wiele przekroczeń. W związku z tym należy zredukować emisję pyłu do wartości dopuszczalnej (w każdym wariancie i w każdej frakcji pyłu redukcję wykonano o ten sam udział procentowy).
St.red= $\frac{max.\ wartosc\ opadu\ pylu - (D_{p - \ }R_{p)\ }}{max.\ wartosc\ opadu\ pylu}$, gdzie Dp − Rp = 180
St.red= $\frac{1244,10327 - 180}{1244,10327}$= 0,855
Wyniki po redukcji przedstawiono w Tabeli16:
Tabela 17. Ładunek pyłu całkowitego po redukcji
-----------------------------------------------------------------------------------
ATMOTERM Opole EK100W
-----------------------------------------------------------------------------------
MAGDALENA GRABOWSKA
-----------------------------------------------------------------------------------
PARAMETRY FRAKCJI PYŁU
Obiekt: CWICZENIE PROJEKTOWE Z ROZPRZESTRZENIANIA ZANIECZYSZCZEN NR 20
Identyfikator obiektu: GRAM Przedział emitorów: od: 1 do: 1
-----------------------------------------------------------------------------------
Nr | Nr | | Frakcje | Udział | Emisja |
emitora |wariantu|Substancja | [um] | proc.[%] | [Mg/rok] |Wf[m/s]
-----------------------------------------------------------------------------------
1| 1 | pył ogółem | 0-10 | 26,28 | 7,65535| 0,001
| | | 10-20 | 33,95 | 9,88962| 0,013
| | | 20-40 | 21,67 | 6,31246| 0,052
| | | 40-60 | 16,47 | 4,79770| 0,145
| | | 60-100| 0,65 | 0,18934460 | 0,370
| | | 100-200| 0,98 | 0,28547340 | 0,735
| 2 | pył ogółem | 0-10 | 26,28 | 13,25128| 0,001
| | | 10-20 | 33,95 | 17,11876| 0,013
| | | 20-40 | 21,67 | 10,92676| 0,052
| | | 40-60 | 16,47 | 8,30474| 0,145
| | | 60-100| 0,65 | 0,32775240 | 0,370
| | | 100-200| 0,98 | 0,49414970 | 0,735
| 3 | pył ogółem | 0-10 | 26,28 | 7,00294| 0,001
| | | 10-20 | 33,95 | 9,04679| 0,013
| | | 20-40 | 21,67 | 5,77449| 0,052
| | | 40-60 | 16,47 | 4,38883| 0,145
| | | 60-100| 0,65 | 0,17320810 | 0,370
| | | 100-200| 0,98 | 0,26114450 | 0,735
---------------------------------------------------------------------------
ATMOTERM Opole EK100W
---------------------------------------------------------------------------
MAGDALENA GRABOWSKA
---------------------------------------------------------------------------
ANALIZA OPADU PYŁU
Punkty z wartościami opadu przekraczającymi progi
Obiekt: CWICZENIE PROJEKTOWE Z ROZPRZESTRZENIANIA ZANIECZYSZCZEN NR 20
Identyfikator obiektu: GRAM Zbiór wyników: D02GRAM.DBF
---------------------------------------------------------------------------
Współrzędne Pył ogółem
X[m] Y[m] [g/m2*rok]
---------------------------------------------------------------------------
Współczynnik szorstkości z0 = 1,00000
--------------------------------------------------------------------------
Obszar zwykły Dp-Rp=180,000
Nie ma przekroczeń
14. Wyznaczenie emisji granicznych według kryteriów imisyjnego i emisyjnego oraz proponowane emisje dopuszczalne:
Tabela 18. Określenie emisji granicznych oraz proponowanych emisji dopuszczalnych
15. Podsumowanie i proponowane rozwiązania:
Podczas wykonywania ćwiczenia projektowego oceniono stan zanieczyszczenia powietrza. Kierowano się kryteriami: emisyjnym (badając stężenia zanieczyszczeń w gazach odlotowych tuż przy wylocie z emitora) i imisyjnym (badając stężenia zanieczyszczeń na poziomie terenu i zabudowy w określonych odległościach od emitora).
Po wykonaniu obliczeń wstępnych, wprowadzono je do programu, który obliczył jak emisja będzie rozprzestrzeniać się na danym obszarze. Następnie określono stopień redukcji emisji i co za tym idzie redukcję zredukowaną jaka odpowiada standardom. Stężenia substancji gazowych SO2, NO2 i substancji stałe tj. pył ogółem zostały zmniejszone. Wyliczona emisja zredukowana spełnia kryteria emisyjne i imisyjne.
W celu osiągnięcia poziomu emisji zredukowanej należy zastosować niżej wymienione metody:
usuwanie SO2: zastosowanie instalacji o lepszej skuteczności odsiarczania
usuwanie NO2: wybrano metodę pierwotną ograniczenia powstawania zanieczyszczeń poprzez optymalizację pracy kotła połączoną ze sterowaniem procesami spalania poprzez regulację temperatury w kotle
usuwanie pyłu całkowitego: ze względu na dużą zawartość popiołu paliwie, zalecana jest zmiana węgla na taki, którego zawartość popiołu jest mniejsza; możliwe jest również zastosowanie filtrów instalacji odpylającej o większej skuteczności
PM10: Elektrofiltry - odpylacze elektrostatyczne spalin. W takim odpylaczu wykorzystano zjawisko ujemnego ładownia się cząsteczek pyłu w gazie zjonizowanym silnym polem elektrycznym. Pod wpływem sił pola elektrycznego zjonizowane cząsteczki dążą do dodatniej elektrody, gdzie neutralizują się i odpadają do leja zbiorczego. Elektrofiltr składa się z: odpylacza elektrostatycznego, zespołu zasilającego. Skuteczność odpylania współczesnych elektrofiltrów osiąga 99,9%. W razie potrzeby można zastosować: odpylacze tkaninowe lub cyklonowe, jednakże trzeba pamiętać, że nie mają one tylu zalet.
16. Izolinie stężeń maksymalnych i średniorocznych przed i po redukcji emisji:
Izolinie stężeń uśrednionych dla godziny i średniorocznych na poziomie terenu oraz opadu pyłu.
Przed:
- CO_M
- NO2_M
- NO2_R
- PM10_M
- PM10_R
- PYL_CALK
- SO2_M
- SO2_R
Po:
- CO_M
- NO2_M
- NO2_R
- PM10_M
- PM10_R
- PYL_CALK
- SO2_M
- SO2_R