Ochrona powietrza projekt

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

  1. 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.

  1. Zakres zadania projektowego

Projekt został wykonany według poniższego schematu:

  1. Obliczenia wstępne

- 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


  1. 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:

Kotłownię wyposażono w instalację do odpylania o skuteczności odpylania ηuj wynoszącej dla poszczególnych frakcji:

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

------------------------------------------------------------------------

  1. 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

  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]

  1. Obliczenia wstępne

    1. 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$$

  1. 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

  1. 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$$

  1. 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 %


$$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$$


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


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


$$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$$


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$$

  1. Emisja substancji gazowych:

Obliczenia wartości emisji poszczególnych zanieczyszczeń dla podanych wariantów pracy

kotłów.

Emisje substancji gazowych:


$$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$$


ENO2 = UNO2

Okres letni:


ENO2I = 6, 0  kg/h

Okres zimowy:


ENO2II = 12, 1  kg/h

Okres zimowy maksymalny:


ENO2III = 16, 1  kg/h


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}$$


$$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

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$$

$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%.

μ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

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$

ρ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

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\ $$

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$$


$$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 -
  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%

maksymalna zawartość CO2 w spalinach

CO2= 18,7% * (1 - $\frac{7}{21}$)= 12,5%

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)

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

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


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

$\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

$\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

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

, 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.

(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}}\ $$

(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}}$$


$$\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}}$$

(Sp)SS6% = , mg/m3n


$$\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:

  1. dla pojedynczego emitora lub zespołu emitorów, z których został utworzony emitor zastępczy:

oraz:

  1. 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 :

  1. st. Red. ze względu na stężenie max= $\frac{max.wartosc\ percentyla - D_{1}}{max.\ wartosc\ percentyla}$

  2. 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:

  1. st.red= $\frac{107,24877 - 200}{107,24877}$= -0,865

  2. st.red= $\frac{7,56571 - 16}{7,56571}$= -1,115 , W= 40∙ (1-0,6) = 16

Dla SO2:

  1. st.red= $\frac{590,31946 - 350}{590,31946}$= 0,407

  2. st.red = $\frac{42,84605 - 8}{42,84605}$= 0,813 , W= 20∙ (1-0,6) = 8

PM10:

  1. st.red= $\frac{124,00422 - 280}{124,00422}$= -1,258

  2. 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:

  1. 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

  2. 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:

  1. gdy Hmax≥ Z (Hmax – najwyższa efektywna wysokość emitora w zespole) – obliczenia wykonujemy od wysokości Z

  2. 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:

  1. usuwanie SO2: zastosowanie instalacji o lepszej skuteczności odsiarczania

  2. usuwanie NO2: wybrano metodę pierwotną ograniczenia powstawania zanieczyszczeń poprzez optymalizację pracy kotła połączoną ze sterowaniem procesami spalania poprzez regulację temperatury w kotle

  3. 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

  4. 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.

  1. Przed:

- CO_M

- NO2­_M

- NO2_R

- PM10_M

- PM10_R

- PYL_CALK

- SO2_M

- SO2_R

  1. Po:

- CO_M

- NO2_M

- NO2_R

- PM10_M

- PM10_R

- PYL_CALK

- SO2_M

- SO2_R


Wyszukiwarka

Podobne podstrony:
Ochrona powietrza projekt
op projekt ochrona powietrza
projekt ochrona powietrza
Ochrona Powietrza 2[P] MagdalenaG TEMAT
OCHRONA POWIETRZA, inżynieria ochrony środowiska kalisz, a pwsz kalisz ioś, IV ochrona powietrza
instrumenty ochrony powietrza oraz metody ich wykorzystania
ochrona powietrza
ochrona powietrza
Ochrona powietrza (zaliczenie ćwiczeń), UMK, Ochrona środowiska
Oczyszczanie Gazow Odlotowych, Ochrona Środowiska studia, 4 rok (2009-2010), Semestr VII (Rok 4), Oc
PYTANIA!!!(2), Ochrona Środowiska studia, 4 rok (2009-2010), Semestr VII (Rok 4), Ochrona Powietrza
Miareczkowanie konduktometryczne kwasu solnego, Studia - IŚ - materiały, Semestr 05, Ochrona powietr
ochrona powietrza
PDiU VI semestr ochrona srodowiska w projektowaniu drog
Ochrona powietrza 2, studia mgr rok 2, semestr II, Prawo Ochrony środowiska
Ochrona powietrza
ochrona powietrza, Technik BHP, Semestr II

więcej podobnych podstron