Modelowanie rozprzestrzeniania zanieczyszczeń przy zakładzie

Modelowanie rozprzestrzeniania zanieczyszczeń przy zakładzie to symulacja dyspersji substancji z emitorów, która wyznacza stężenia w siatce obliczeniowej i zasięg oddziaływania. Uwzględnia emisję, parametry komina, meteorologię i teren. Dla NO2, SO2, PM10 i opadu pyłu pozwala porównać maksymalne stężenia z normami. W analizowanym przypadku wartości pozostają z dużym marginesem poniżej dopuszczalnych poziomów.

Czym jest modelowanie rozprzestrzeniania zanieczyszczeń przy zakładzie i dlaczego jest niezbędne

Modelowanie rozprzestrzeniania zanieczyszczeń przy zakładzie to matematyczno-fizyczny opis tego, jak substancje uwolnione z emitorów trafiają do powietrza i jak ulegają dyspersji. Zamiast zgadywać, czy dym z komina dotrze do najbliższej zabudowy, model symuluje stężenia w siatce obliczeniowej i wyznacza zasięg oddziaływania zakładu.

Więcej informacji znajdziesz tutaj: skontaktuj się z nami.

Punktem wyjścia jest emisja, czyli ilość substancji wprowadzanej do powietrza w jednostce czasu, oraz parametry źródła. Dla komina w zakładzie Edison Next Poland Sp. z o.o. uwzględnia się m.in. wysokość 150 m i średnicę 2,6 m. Wysokość pozornego punktu emisji opisuje wzór H = h + Δh, gdzie h to rzeczywista wysokość komina, a Δh — wyniesienie smugi.

Symulacja wymaga też danych meteorologicznych — kierunku i prędkości wiatru, stabilności atmosfery, temperatury — a także ukształtowania terenu i charakterystyki substancji. Obliczenia mogą obejmować NO2, SO2, PM10 i opad pyłu, z uśrednianiem 1-godzinnym i rocznym. Wykorzystuje się wówczas programy takie jak Ek100W (Atmoterm), oparty na modelu Pasquilla, zgodnie z metodyką referencyjną z rozporządzenia z 26 stycznia 2010 r. (Dz.U. 2010 nr 16 poz. 87).

Dane emisyjne: parametry komina, rodzaj i ilość zanieczyszczeń.
Dane meteorologiczne: wiatr, stabilność, temperatura, wilgotność.
Dane terenowe: siatka obliczeniowa, zabudowa, odległości od emitora.

Modelowanie pozwala ocenić, czy maksymalne stężenia NO2, SO2 i PM10 mieszczą się w normach. W analizowanym przypadku nawet najwyższe wartości pozostają daleko poniżej dopuszczalnych poziomów.

Rodzaje zanieczyszczeń emitowanych przez zakłady przemysłowe i ich wpływ na otoczenie

Zanieczyszczenia z zakładów przemysłowych dzieli się zwykle na gazowe, pyłowe oraz związane z cząstkami stałymi — metale ciężkie i WWA. Każda grupa inaczej rozprzestrzenia się w powietrzu i oddziałuje na ludzi, rośliny, glebę oraz wodę.

Tlenki azotu (NOx, w tym NO2)— powstają głównie podczas spalania paliw w wysokiej temperaturze. Drażnią układ oddechowy, współtworzą smog fotochemiczny i ozon troposferyczny, przyczyniają się także do eutrofizacji i zakwaszenia gleb.
Dwutlenek siarki (SO2)— pochodzi ze spalania paliw zawierających siarkę. Podrażnia drogi oddechowe, nasila choroby płuc, a w środowisku sprzyja kwaśnym deszczom.
Pyły zawieszone PM10 i PM2.5— mieszaniny cząstek stałych i kropelek. PM10 zatrzymuje się w górnych drogach oddechowych, natomiast PM2.5 przenika głębiej, do pęcherzyków płucnych; oba sprzyjają chorobom serca i płuc. Opad pyłu osadza się na roślinach i glebie, zmieniając ich właściwości.
Metale ciężkie(np. ołów, kadm, rtęć, arsen) — są toksyczne i kumulują się w glebie, osadach oraz organizmach, skąd mogą przenikać do łańcucha pokarmowego.
WWA— wielopierścieniowe węglowodory aromatyczne, w tym benzo(a)piren, powstają przy niepełnym spalaniu. Mają działanie rakotwórcze i mutagenne.

Przy zakładzie modeluje się najczęściej substancje o ustalonych normach i typowe dla danego procesu. W analizach zakładowych powtarza się zestaw: NO2, SO2, PM10 oraz opad pyłu. Dla Edison Next Poland Sp. z o.o. właśnie te zanieczyszczenia objęto obliczeniami, a najwyższe wartości mieściły się z dużym marginesem poniżej dopuszczalnych poziomów. PM2.5, metale ciężkie i WWA również bywają modelowane — zwłaszcza gdy zakład spala paliwa stałe lub prowadzi procesy wysokotemperaturowe.

Modele dyspersji i metody obliczeniowe stosowane w ocenie rozprzestrzeniania

W praktyce zakładowej wybór modelu wynika z rodzaju emisji, wymaganej dokładności i charakteru terenu. Dla typowych emitorów punktowych, takich jak kominy, najczęściej stosuje się modele gaussowskie. Opisują one smugę jako rozkład normalny wokół osi wiatru; stężenie zależy wówczas od emisji, odległości, wysokości pozornego punktu emisji oraz warunków meteorologicznych. Model Gaussa jest szybki i dobrze sprawdza się w ocenach zgodności z normami, gdy teren nie jest silnie zabudowany ani pofałdowany.

Model Lagrange’a śledzi trajektorie pakietów cząstek w zmiennym polu wiatru. Sprawdza się przy zmiennej meteorologii, dłuższych okresach uśredniania i złożonych źródłach, ale wymaga większej liczby danych wejściowych. Z kolei CFD, czyli obliczeniowa mechanika płynów, rozwiązuje równania przepływu w siatce 3D. Znajduje zastosowanie w terenie zurbanizowanym, przy zabudowie, przeszkodach i emisjach niskich, gdzie modele uproszczone zawodzą. CFD jest jednak kosztowny obliczeniowo i wymaga kalibracji.

Parametry wejściowe obejmują geometrię emitora, prędkość i temperaturę spalin, strumień masy zanieczyszczeń, współczynnik dyfuzji, klasę stabilności atmosfery, profil wiatru i temperatury, a także rzeźbę terenu i przeszkody.

Emisja punktowa, teren otwarty— model Gaussa.
Emisja zmienna, duży obszar— model Lagrange’a.
Zabudowa, wąwozy, emisja niska— CFD.

Dobór modelu to kompromis między dokładnością, czasem obliczeń i dostępnością danych. W ocenie przy zakładzie najpierw stosuje się model gaussowski, a gdy wyniki są bliskie granicom lub teren jest skomplikowany, weryfikuje się je modelem Lagrange’a albo CFD.

Rola warunków meteorologicznych i ukształtowania terenu w symulacji rozprzestrzeniania

Warunki meteorologiczne decydują o zasięgu i kształcie smugi. Największe znaczenie mają prędkość i kierunek wiatru: prędkość steruje wynoszeniem, rozcieńczaniem i czasem kontaktu z podłożem, a kierunek wyznacza sektor, w którym mogą pojawić się podwyższone stężenia. Dla komina o wysokości 150 m i średnicy 2,6 m nawet niewielka zmiana kierunku przesuwa maksimum obliczeń, natomiast spadek prędkości wiatru sprzyja kumulacji zanieczyszczeń. Temperatura i wilgotność wpływają na pionowe mieszanie oraz wyniesienie gazów odlotowych.

Inwersja termiczna zwiększa stężenia, ponieważ odwraca normalny profil temperatury: powietrze przy ziemi jest chłodniejsze, a wyżej — cieplejsze. Powstaje warstwa blokująca, która ogranicza pionową dyfuzję. Smuga nie może się wznieść i rozprzestrzenia się poziomo, często tuż przy ziemi. Przy słabym wietrze rozcieńczanie jest minimalne, więc rosną zarówno stężenia 1-godzinne, jak i średnie roczne — szczególnie w obniżeniach terenu.

Ukształtowanie terenu modyfikuje przepływ lokalny. Chropowatość zwiększa turbulencję mechaniczną, spowalnia wiatr i wydłuża zaleganie zanieczyszczeń. Wzgórza, doliny, wąwozy i zabudowa tworzą kanały oraz strefy cienia aerodynamicznego, dlatego model musi uwzględniać rzeźbę, a nie tylko dominującą różę wiatrów. Ten sam emitor w terenie otwartym i w dolinie da inne izolinie; pominięcie inwersji lub topografii może zaniżyć wyniki w lokalnych ogniskach, nawet jeśli maksymalne stężenia pozostają daleko poniżej dopuszczalnych poziomów.

Interaktywny kalkulator stężeń zanieczyszczeń i przykładowe obliczenia krok po kroku

Kalkulator stężeń służy do wstępnego oszacowania stężenia zanieczyszczeń i zasięgu wpływu zakładu. Nie zastępuje pełnej oceny, ale pozwala szybko sprawdzić, które parametry emisji i meteorologiczne najmocniej zmieniają wynik. Poniżej zestawiono dane wejściowe oraz przykładowe obliczenia.

Parametry do wprowadzenia:

geometria emitora: wysokość 150 m i średnica 2,6 m,
prędkość i temperatura spalin oraz strumień masy dla NO2, SO2, PM10 i opadu pyłu,
czasy uśredniania: 1-godzinny i roczny,
meteorologia: prędkość i kierunek wiatru, klasa stabilności atmosfery według Pasquilla, temperatura powietrza,
teren: chropowatość, przeszkody, zabudowa.

Przykładowe obliczenia krok po kroku:

Przyjmij rzeczywistą wysokość komina h = 150 m. Na podstawie prędkości wiatru i różnicy temperatur spalin oszacuj wyniesienie Δh. Oblicz wysokość pozornego punktu emisji H = h + Δh.
Wprowadź średnicę 2,6 m i prędkość wylotową, aby wyznaczyć strumień objętości gazów.
Wybierz substancje: NO2, SO2, PM10 oraz opad pyłu. Ustaw czasy uśredniania: 1-godzinny i roczny.
Zdefiniuj warunki meteorologiczne, w tym klasę stabilności Pasquilla i profil wiatru.
Uruchom kalkulator oparty na modelu Pasquilla — np. Ek100W (Atmoterm), zgodny z rozporządzeniem Ministra Środowiska z 26 stycznia 2010 r. (Dz.U. 2010 nr 16 poz. 87).
Odczytaj izolinie stężeń i zasięg wpływu. Porównaj maksymalne wartości z obowiązującymi normami.

W analizowanym przypadku nawet najwyższe obliczone stężenia NO2, SO2 i PM10 w siatce obliczeniowej mieszczą się z dużym marginesem poniżej dopuszczalnych poziomów. Kalkulator pozwala jednak testować scenariusze: zmianę kierunku wiatru, inwersję termiczną, spadek prędkości wiatru czy modyfikację emisji. Dzięki temu można szybko porównać warianty jeszcze przed pełną symulacją i oceną zgodności.

Wskazówka: kalkulacja_pl — dla wiarygodnego wyniku nie pomijaj Δh ani danych o terenie. Zasięg wpływu rośnie, gdy wiatr słabnie, a smuga zalega w obniżeniach.

Sposoby zapobiegania zanieczyszczeniu powietrza i redukcji emisji z zakładu

Zapobieganie zanieczyszczeniu powietrza powinno być wbudowane w proces, a nie dodawane na końcu. Najpierw eliminuje się lub ogranicza powstawanie zanieczyszczeń, następnie hermetyzuje instalację, a dopiero potem oczyszcza gazy odlotowe. W instalacjach emitujących NO2, SO2, PM10 i opad pyłu kluczowe znaczenie mają technologie redukcji dopasowane do rodzaju substancji.

Odsiarczanie — usuwa SO2 metodami mokrymi, półsuchymi lub suchymi; reagenty wiążą siarkę w produkty możliwe do zagospodarowania.
Odazotowanie — ogranicza NOx/NO2 przez niskoemisyjne palniki, recyrkulację spalin oraz selektywną redukcję katalityczną (SCR) lub niekatalityczną (SNCR).
Filtry workowe — skutecznie wychwytują pył, w tym PM10 i drobne frakcje; wymagają doboru tkaniny, właściwego zakresu temperatur i regularnej wymiany worków.
Odpylanie wstępne i scrubbery — elektrofiltry, cyklony i płuczki ograniczają unos pyłów oraz emisję gazów kwaśnych i odorowych.
Najlepsze dostępne techniki BAT — porządkują wybór technologii według procesu, skali, paliwa i kosztów; obejmują także monitorowanie, konserwację i zarządzanie awaryjne.

Nie istnieje jedna uniwersalna technologia. W zakładach przemysłowych najlepiej sprawdza się układ łączony: niskoemisyjne spalanie, odazotowanie, odpylanie i odsiarczanie. Taki ciąg redukuje stężenia NO2, SO2, PM10 oraz opadu pyłu, a po modernizacji efekt warto potwierdzić pomiarami i ponowną symulacją rozprzestrzeniania.

Co najbardziej zanieczyszcza środowisko – porównanie źródeł emisji i udział zakładu

Zanieczyszczenie środowiska w skali lokalnej i globalnej tworzy suma różnych źródeł, nie pojedynczy emitor. W miastach dominuje transport drogowy (tlenki azotu, pyły ze ścierania), a za nim przemysł i energetyka. Emisja niska — kotły domowe, piece węglowe — odpowiada za znaczną część PM10 i PM2.5 w sezonie grzewczym, szczególnie w zabudowie rozproszonej, gdzie smuga z niskich kominów nie rozprasza się efektywnie.

Jak porównać udział zakładu z innymi źródłami? Punktem wyjścia jest bilans emisji: ilość NO2, SO2, PM10 i opadu pyłu wprowadzana do powietrza przez zakład w skali roku w zestawieniu z emisją transportu, energetyki i emisji niskiej na tym samym obszarze. Drugi krok to modelowanie rozprzestrzeniania — dla Edison Next Poland Sp. z o.o. maksymalne obliczone stężenia NO2, SO2 i PM10 w siatce obliczeniowej pozostają znacznie poniżej dopuszczalnych poziomów. Emitor o wysokości 150 m i średnicy 2,6 m wprowadza zanieczyszczenia wysoko, więc ich udział w stężeniach przyziemnych jest relatywnie mniejszy niż w przypadku niskich źródeł.

W praktyce porównanie sprowadza się do trzech pytań:

Jaka jest roczna masa emisji z zakładu na tle transportu, energetyki i emisji niskiej w promieniu kilku kilometrów?
Jak wysokość emitora wpływa na stężenia przyziemne w porównaniu z niskimi źródłami?
Czy modelowanie pokazuje przekroczenia norm, czy — jak dla Edison Next Poland — wartości pozostają daleko poniżej dopuszczalnych poziomów?

Przykładowe dane wejściowe i wyniki pomiarów do demonstracji kalkulatora

Do odtworzenia obliczeń przygotuj anonimizowany zestaw: dane wejściowe, próbki emisji, profile wiatrów oraz wyniki symulacji i pomiary stężeń. Dla zakładu Edison Next Poland Sp. z o.o. punktem odniesienia są emitor o wysokości 150 m i średnicy 2,6 m, substancje NO2, SO2, PM10 i opad pyłu oraz czasy uśredniania: 1-godzinny i roczny.

Dane wejściowe: geometria emitora, prędkość i temperatura gazów odlotowych, strumień masowy zanieczyszczeń, współrzędne, rzędna terenu i chropowatość podłoża.
Próbki emisji: wyniki pomiarów albo bilansu dla NO2, SO2, PM10 i opadu pyłu; podaj jednostkę, czas próbkowania i niepewność.
Profile wiatrów: kierunek, prędkość i klasa stabilności dla godzin symulacji. Szukaj ich w IMGW-PIB, reanalizie ERA5, modelach WRF/CALMET, a najlepiej weryfikuj własnymi pomiarami (dane_wlasne) z anemometru lub sodaru.
Wyniki symulacji: siatka obliczeniowa, izolinie i maksima stężeń. W przykładzie dla Edison Next Poland maksymalne wartości NO2, SO2 i PM10 nie zbliżają się do poziomów dopuszczalnych.
Pomiary stężeń: jeśli są dostępne, zestaw je z symulacją w tych samych punktach i okresach; pokaże to, czy kalkulator nie zaniża ani nie zawyża wyników.

Obliczenia wykonaj w programie Ek100W (Atmoterm), opartym na modelu Pasquilla, zgodnie z rozporządzeniem Ministra Środowiska z 26 stycznia 2010 r. (Dz.U. 2010 nr 16 poz. 87). Przy wyniesieniu smugi stosuj H = h + Δh, a wynik zapisz jako plik CSV z kolumnami: x, y, stężenie, czas uśredniania.

Lista pytań technicznych i niezaadresowanych luk w analizie do uzupełnienia

Poniższe pytania porządkują braki, które trzeba domknąć, zanim wyniki modelowania trafią do publikacji lub FAQ.

Parametry emisji: Czy dla każdej substancji podano strumień masowy, stężenie na wylocie, prędkość i temperaturę gazów, wilgotność oraz zawartość tlenu? Czy dane pochodzą z pomiarów czy bilansu? Jak wygląda zmienność w rozruchu, postoju i awarii? Czy uwzględniono niepewność pomiaru i reprezentatywność próbek?
Model dyspersji: Dlaczego wybrano właśnie ten model i program do wizualizacji izolinii, a nie model Lagrange’a, CFD lub hybrydowy? Jak uzasadniono współczynniki dyspersji, pole wiatru i wyniesienie smugi? Czy model obsługuje chemię NO/NO2, depozycję, obecność budynków i przeszkód terenowych? Czy wyniki porównano z alternatywnym modelem?
Scenariusze meteorologiczne: Które lata, sezony i pory dnia są reprezentatywne? Czy uwzględniono klasy stabilności, cisze, inwersje, wiatry przeciwne do receptorów i epizody smogowe? Czy użyto danych lokalnych, reanalizy, WRF/CALMET, a może pomiarów własnych? Jak wybrano rok referencyjny?
Monitoring i walidacja: Gdzie i jak często mierzyć stężenia? Czy planowany jest monitoring ciągły czy kampanie? Jak zestawić wyniki z modelem w tych samych punktach i okresach? Jakie wskaźniki walidacji zastosować: korelację, FAC2, NMSE, błąd systematyczny? Jak odejmować tło regionalne?
Niepewność modelowania: Czy wykonano analizę wrażliwości na parametry emisji, meteorologię, siatkę i rzędną terenu? Jak propagować niepewność do izolinii i wartości maksymalnych? Czy wniosek o dużym marginesie do norm jest odporny na założenia pesymistyczne? Czy przedziały ufności będą raportowane?