Analiza Activitatii Unei Statii de Betoane Asupra Calitatii Aerului

Lucrare de licență

ANALIZA ACTIVITATII UNEI STATII DE BETOANE ASUPRA CALITATII AERULUI

CUPRINS

1. INTRODUCERE

2. TRANSPORTUL ȘI DISPERSIA POLUANȚILOR ÎN ATMOSFERĂ

2.1 Forme de transport ale poluanților

2.2 Factori care influențează dispersia poluanților în atmosferă

3. LEGISLATIA PRIVIND TRANSPORTURILE POLUANTILOR ATMOSFERICI

4. ASPECTE DE MEDIU

5. BAZELE TEORETICE A DISPERSIEI POLUANTILOR IN ATMOSFERA

5.1. Ecuațiile de difuzie

5.2 Dispersia în condiții de extremă stabilitate

5.3 Înălțimea efectivă de emisi

5.4 Controlul emisiilor în aer

5.5. Fenomenul de sinergism al agenților poluanți ai atmosferei

5.6. Considerații generale privind modelarea matematică a transportului și dispersiei poluanților în atmosferă

5.7 Descrierea fenomenologică a poluării atmosferei din surse punctiforme

5.8. Elemente de aerostatica și aerodinamică care intervin la studiul și evoluția penei

5.8.1 Elemente de aerostatică

5.8.2Elemente de aerodinamică

5.9. Distribuția pe verticală a vitezei curenților de aer în stratul limită

4.10. Distribuția pe verticală a temperaturii aerului

5.11. Criterii de stabilitate a atmosferei

5.12. Influența reliefului asupra structurii curenților atmosferici și asupra formei penei

5.13. Modele matematice pentru calculul transportului si dispersiei poluantilor in atmosfera

6 PROGRAMUL DE CALCUL UTILIZAT PENTRU MODELAREA DISPERSIEI POLUANTILOR IN ATMOSFERA

6.1. Principalele variabilele folosite în program

6.2 Instrucțiuni de utilizare

6.3 Concluzii referitoare la programul de calcul

7. Studiu de caz – Statie de beton mobila amplasata pe soseua de centura Târgu Mures

8. CONCLUZII

9. BIBLIOGRAFIE

INTRODUCERE

Tema proiectului este analiza activitatii unei statii de betoane asupra calitatii aerului . În Romania statiile de betoane mobile se impun cu rapiditate pe o piata cu cerinte speciale de productie beton, necesar de reducere a costurilor de productie a betonului si nu în putine cazuri mobilitatea statiilor de beton la santierele de constructii.

Statia de beton mobilă permite să fie instalata pe santierul clientului într-o singură zi, ceea ce permite să aibă costuri reduse de relocare cât si o productivitate crescută de beton, atât prin faptul că are un timp extrem de redus în care nu productia este 0 cât si faptul că are un sistem hidraulic independent si nu este necesara o macara pentru montarea ei.

Statia de beton mobilă are nevoie pentru mutarea pe un alt santier de productie doar de un singur cap de tractor.Statia de betoane mobila oferă controlul total clientului asupra productiei de beton cât si al raportărilor de productie prin intermediul sistemului de automatizare, al casetei de comandă si al sistemului informatic cu care este echipata statia încă de la constructia în fabrica noastră.

Statia de beton mobilă prezentată este situată pe soseaua de centură a orasului Târgu-Mures în partea de E-NE la distanta de 1 km de la începerea soselei în partea de S.

Programul în care am lucrat este structurat sub formă modulară astfel încât permite perfecționarea ulterioară și în acest sens sunt în vedere adăugarea unor noi module de calcul pentru calculul concentrațiilor de poluant provenite din surse liniare, surse intermitente, surse aflate în mișcare.

Cu ajutorul programului am obtinut atât concentratile maxime ale poluantului studiat cât si evolutia spatiala norului de poluant în diferite ipoteze de calcul.

2. TRANSPORTUL ȘI DISPERSIA POLUANȚILOR ÎN ATMOSFERĂ

2.1 Forme de transport ale poluanților

Transportul poluanților în mediu la distanțe apreciabile față de sursa de poluare, se datorează unei serii de fenomene: advecție-dispersie, fenomene întârzietoare care determină mărimea ariei poluate.

Convecția (advecția) reprezintă transportul substanțelor de către curentul de aer, apă de suprafață sau apă subterană, la viteza medie a acestuia.

Dispersia este rezultatul acțiunii simultane a unui fenomen pur mecanic și a unui fenomen fizico-chimic.

Dispersie = Advecție (Transport) + Diluție (Difuzie)

Dispersia moleculară – difuzia se produce ca urmare a diferențelor de concentrații ale poluanților, evident de la concentrații mai mari către concentrații mai mici.

2.2 Factori care influențează dispersia poluanților în atmosferă

Dispersia poluanților în atmosferă, stabilirea gradientului de poluare al acestuia și în final determinarea concentrațiilor de poluanți sunt puternic influențate de condițiile climaterice și meteorologice.

Temperatura

Variația temperaturii aerului în funcție de presiune și de înălțime, este un factor important care intervine în deplasarea maselor de aer și implicit în răspândirea nocivităților în atmosferă.

Temperatura aerului variază cu înălțimea. Această variație are o valoare cuprinsă între 0.4 și 0.9 C la fiecare sută de metri.

Variația temperaturii aerului cu înălțimea se numește gradient vertical de temperatură și se exprimă prin raportul T/Z în care T este temperatura aerului în C iar Z înălțimea în metri.

Întrucât, în general, temperatura aerului scade cu înălțimea, media gradientului vertical de temperatură este pozitivă. În funcție de mărimea gradientului vertical de temperatură se pot deosebi mai multe tipuri de stratificare a atmosferei.

Instabilitatea

Aceasta este caracterizată prin scăderea temperaturii cu înălțimea la o valoare mai mare decât media gradientului vertical. Această condiție se realizează cel mai frecvent în zilele senine de vară când încălzirea rapidă a solului de către soare va influența temperatura straturilor de aer dând naștere unui amestec vertical pronunțat de curenți ascendenți puternici.Masa de aer în ascensiune se răcește și își continuă deplasarea fiind mai puțin densă decât atmosfera înconjurătoare. Această situație de instabilitate față de mediul înconjurător va favoriza dispersia substanțelor nocive.

Stabilitatea

Aceasta este caracterizată prin faptul că scăderea temperaturii cu înălțimea este mai puțin accentuată decât gradientul vertical mediu. În această situație masa de aer în ascensiune se va răci cu mai puțin decât această medie.

Atmosfera înconjurătoare, considerată în repaus, dispune de o energie termică mai mare, capabilă să frâneze mișcarea ascendentă a aerului.

Tendința de a rezista deplasărilor verticale este caracteristică stabilității, starea atmosferică în care difuzia este aproape absentă, ceea ce defavorizează răspândirea poluanților.

Starea indiferentă

Această stare a atmosferei, numită și neutră, prezintă o scădere a temperaturii cu înălțimea după o relație adiabatică, adică în repartiția verticală a temperaturii aerului din atmosfera înconjurătoare și din masa de aer asupra deplasării nu există nici o diferență.

Volumul de aer în ascensiune va avea aceeași temperatură cu straturile de aer care îl înconjoară și va fi echilibru la orice nivel, neavând tendința să revină la poziția inițială.

Stările atmosferice deosebite sunt izotermiile și inversiunile termice.

Izotermia este caracterizată de faptul că temperatura aerului este staționară cu înălțimea.

În cazul inversiunilor termice temperatura este crescătoare cu înălțimea, ceea ce implică un gradient vertical de temperatură negativă.

Aceste stări ale atmosferei pot surveni datorită răcirii aerului de lângă sol ca urmare a pierderii de căldură de către Pământ prin intermediul radiațiilor cu lungimi mari de undă în nopțile senine. [1]

Figura 2.1. Inversiune terimica

O situație foarte stabilă ar avea loc atunci când aerul rece este aproape de sol, și există un strat de aer cald deasupra lui. Aceasta este o abatere de la situația normală, și se numește inversiune de temperatură, sau inversiune termică. Acest aer rece este mai dens decât aerul cald. Mai presus de stratul de inversiune, aerul va fi răcoros cu altitudinea din nou în creștere.

Umiditatea

Umiditatea aerului se opune difuziei poluanților și respectiv micșorării concentrațiilor lor, împiedicând particulele să se deplaseze. Umiditatea crescută duce la formarea de ceață care produce concentrarea impurităților. În zonele poluate, ceața se produce frecvent.

Vântul

Vântul definit ca mișcare orizontală a aerului este cel mai important factor care conduce la răspândirea poluanților în atmosferă. El nu este un element stabil, viteza lui putând varia în limite largi între rafală și acalmie.

Direcția vântului reprezintă caracteristica principală a mișcării aerului în relația ei cu poluarea, deoarece pe direcția de mișcare se produce poluarea cea mai intensă, în timp ce în alte zone, aerul poate rămâne curat.

Prin urmare vântul reprezintă procesul prin care are loc deplasarea poluanților și creează fenomenul prin care aceștia se difuzează în bazinul aerian. Difuzia este direct proporțională cu viteza vântului.

Vântul uniform și de viteză mică menține concentrațiile ridicate de poluanți în stratul de aer în care a ajuns.

Cu cât viteza vântului este mai mare, cu atât înălțimea la care ajung poluanții este mai mică. Vitezele mari ale vântului împing spre sol pana de gaze sau poluanți eliminați de surse, din cauza rezistenței de frecare mai mare pe care o întimpină aerul din apropierea solului și a vitezei mai mari de la înălțime.

Turbulența

Mișcarea curenților de aer este în general o mișcare turbulentă deci cu o structură haotică, întâmplătoare, volumele mici de aer aflate în mișcare nepăstrându-și individualitatea.

Din punct de vedere al modului cum se produce turbulența, aceasta poate fi mecanică și termică.

Turbulența mecanică are ca sursă principală frecarea aerului cu suprafața pământului și a altor obstacole. Ea se propagă în sus până la înălțimi mari și crește cu rugozitatea pământului și cu intensitatea vântului.

Turbulența termică sau convecția ia naștere ca urmare a diferenței de temperatură dintre pământ și curentul de aer, deci ca urmare a repartizării neuniforme a temperaturii.

Ambele feluri de turbulență se găsesc concomitent în atmosferă, în diverse proporții, funcție de condițiile meteo.Criteriile după care se apreciază turbulența aerului sunt următoarele:

Criteriul dinamic (Reynolds)

Criteriul termic (Richardson)

Turbulența este un factor care favorizează amestecurile și de aceea difuzia impurităților în masa de aer turbulent se face mai repede.

3. LEGISLATIA PRIVIND TRANSPORTURILE POLUANTILOR ATMOSFERICI

In Romania legislatia privind calitatea aerului si nivelul maxim de poluare este reglementat de urmatoarea lege :

Legea nr. 104/15.06.2011 privind calitatea aerului înconjurător (publicată în Monitorul Oficial nr. 452/28.06.2011) si are ca scop protejarea sanatatii omului si a mediului prin reglementari destinate mentinerii calitatii aerului sub limita maxima admisa de poluare [2]

Tabelul nr. 3.1 Dioxidul de sulf – SO2 conform [3]

Tabelul nr. 3.2 Oxizi de azot -NOx – conform [4]

Tabelul nr 3.3 Plumb- Pb- conform [5]

Tabelul nr 3.4 Nivelul metalelor grele ( Pb, As, Cd si Ni )- conform [6]

Tabelul 3.5 Nivelul pulberilor in suspensie ( PM10 / PM 2,5 ) – conform [7]

Acte legislative in vederea producerii betonului :

1. Ordinul ministrului dezvoltării lucrărilor publice și locuinței nr.577/2008 privind aprobarea Normativului pentru producerea betonului și executarea lucrărilor din beton, beton armat și beton precomprimat – Partea 1: Producerea betonului, indicativ NE 012/1-2007 publicat in Monitorul Oficial al României, Partea I, nr. 374 din 16 /05/2008 [8]

2. Ordinul ministrului dezvoltării regionale și turismului nr. 2514/2010 privind aprobarea Normativului pentru producerea betonului și executarea luct se face mai repede.

3. LEGISLATIA PRIVIND TRANSPORTURILE POLUANTILOR ATMOSFERICI

In Romania legislatia privind calitatea aerului si nivelul maxim de poluare este reglementat de urmatoarea lege :

Legea nr. 104/15.06.2011 privind calitatea aerului înconjurător (publicată în Monitorul Oficial nr. 452/28.06.2011) si are ca scop protejarea sanatatii omului si a mediului prin reglementari destinate mentinerii calitatii aerului sub limita maxima admisa de poluare [2]

Tabelul nr. 3.1 Dioxidul de sulf – SO2 conform [3]

Tabelul nr. 3.2 Oxizi de azot -NOx – conform [4]

Tabelul nr 3.3 Plumb- Pb- conform [5]

Tabelul nr 3.4 Nivelul metalelor grele ( Pb, As, Cd si Ni )- conform [6]

Tabelul 3.5 Nivelul pulberilor in suspensie ( PM10 / PM 2,5 ) – conform [7]

Acte legislative in vederea producerii betonului :

1. Ordinul ministrului dezvoltării lucrărilor publice și locuinței nr.577/2008 privind aprobarea Normativului pentru producerea betonului și executarea lucrărilor din beton, beton armat și beton precomprimat – Partea 1: Producerea betonului, indicativ NE 012/1-2007 publicat in Monitorul Oficial al României, Partea I, nr. 374 din 16 /05/2008 [8]

2. Ordinul ministrului dezvoltării regionale și turismului nr. 2514/2010 privind aprobarea Normativului pentru producerea betonului și executarea lucrărilor din beton, beton armat și beton precomprimat-Partea 2: Executarea lucrărilor din beton, indicativ NE 012/2-2010 publicat in Monitorul Oficial al României, Partea I, nr. 853 din 20 decembrie 2010 [9]

3. Ordinul ministrului transporturilor construcțiilor și turismului/ administrației și internelor nr.1.822/394/2004 pentru aprobarea Regulamentului privind clasificarea si încadrarea produselor pentru construcții pe baza performanțelor de comportare la foc, cu modificările și completările ulterioare publicat in Monitorul Oficial al României, Partea I nr. 90 din 27/01/2005 [10]

4. ASPECTE DE MEDIU

Urmatoarele acte normative reglementeaza realizarea studiului de impact asupra mediului :

Ordonanța de Urgență nr. 195/30 decembrie 2005, aprobată de legea 265/2006, privind protecția mediului; [11]

Ordinul 184/1997, emis de Ministrul Apelor, Pădurilor și Protecției Mediului, pentru aprobarea procedurii de realizare a bilanțurilor de mediu; [12]

Ordinul 860/2002, emis de Ministrul Apelor și Protecției Mediului, pentru aprobarea procedurii de evaluare a impactului asupra mediului; [13]

Ordinul 863/2002, emis de Ministrul Apelor și Protecției Mediului, pentru aprobarea ghidurilor metodologice de evaluare a impactului asupra mediului. [14]

Statia de beton mobila asupra caruia este facut studiul are o capacitate maxima de producere a betonului de 80mc3/h. Este de productie italiana construita in anul 2008 si se afla pe o platforma de 250 mp2

Alimentarea energetica se face din reteaua municipiului Tg. Mures si in caz de defectiuni sau caderi de curent se foloseste un generator de curent

Structura statiei de beton :

1 – 2 silozuri de ciment – capacitate de 140 tone

2 – Schip transportor – transfera agregate uscate la malaxor

3 – Platforma cu cantar pentru schip – capacitatate 2 tone

4 – 4 containare tip palnie – se folosesc pentru depozitarea agregatelor

5 – Supapa de siguranta – amplasata inaintea intrarii cimentului sub presiune in filtru

6 – Filtru cu cartuse – retine pulberile fine emanate

7 – Structura de ancorare pentru silozuri –fixeaza silozurile pentru siguranta

8 – 2 snecuri melcate – capacitate de 80 tone pe ora

9 – Cantar pentru ciment – capacitate de o tona

10- Cantar pentru apa – capacitate de 800 litrii

11 – Instalatie de pompare a apei – putere de 8 cp

12 – Cantar dozator de apa – bazat pe impulsuri de 1 litru

13 – Malaxor cu ax dublu orizontal – 2 mc

14- Palnie de descarcare a malaxorului

15- Sistem de curatare manuala a malaxorului – cu apa sub presiune la temperatura ridicata

16- Instalatie pneumatica cu compresor – se foloseste la functionarea pneumatica a utilajelor

17- Panou electric

18- Panou de control

Figura 4.1 Schema procesului tehnologic de producere a betonului

Tabelul 4.1 EVALUAREA IMPACTULUI ASPECTELOR DE MEDIU

Dupa evaluarea impactului aspectelor de mediu consideram ca toate valorile peste 60 de puncte sunt considerate aspecte semnificative de mediu

Tabelul 4.2 LISTA ASPECTELOR SEMNIFICATIVE DE MEDIU

Deseurile provenite din activitatea de productie a betonului :

Din procesul de producere a betonului nu rezulta deseuri iar apa care apare din diferite scurgeri sau spalari a malaxorului sau a diferitelor utilaje se capteaza printr-un sistem de canale intr-un bazin colector de decantare. Apa se refoloseste la producerea betonului de a doua zi. Deseurile menajere nu sunt o problema pentru ca zona de producere a betonului este dotata cu cosuri de cunoi iar acestea sunt colectate saptamanal de o firma specializata.

Zgomotul provenit de la statia mobila de beton:

Zgomotul este un factor important al statiei de beton, el provine de la motorul electric al statiei, de la utilajele pneumatice, malaxor, pompe. De asemenea un factor important este si vibratia. Avand in vedere ca zona de amplasare a statiei este foarte putin populata si nu exista case locuibile in jurul ei, rezulta ca nici zgomotul si nici vibratia nu cauzeaza o poluare fonica a locuitorilor. Cu toate acestea constructorii statiei de betoane au inconjurat malaxorul cu pereti de tip senvis ( spuma poliuretanică si tabla) cu toate ca fisa tehnica a utilajului arata o poluare fonica provenita de la malaxor mai mica de 80 db.

5. BAZELE TEORETICE A DISPERSIEI POLUANTILOR IN ATMOSFERA

În prezentările obișnuite ale dinamicii mediilor continue compoziția fluidelor este considerată omogenă în întreg domeniul spațio-temporal în care mișcarea este studiată. Pentru un amestec de lichide sau gaze, a căror compoziție este variabilă în timp și spațiu, ecuațiile clasice ale mecanicii fluidelor trebuie modificate.

Se tratează un amestec de două componente fluide. Compoziția fluidelor va fi descrisă cu ajutorul concentrației ci, care se definește ca raportul masei unui component la masa unui volum de control (evident egală cu suma masei tuturor constituenților amestecului):

; (i = 1,2)

Odata cu variația timpului se va modifica și concentrația componentelor amestecului de fluide, pe care în cele ce urmează îl vom numi simplu fluid, atunci când nu există pericol de confuzie. Variația concentrației ia două forme:

variație datorată mișcării macroscopice a fluidului, la care participă fiecare particulă elementară de compoziție invariabilă. Prin această mișcare are loc un amestec numit amestec mecanic pur. Cu toate că, pentru fiecare particulă elementară de amestec fluid în mișcare, compoziția nu variază, ea este variabilă de la punct la punct în spațiul fizic bine definit și fixat al mișcării. În același timp au loc procese de transfer ale căldurii și de frecare internă în amestecul de fluid, ceea ce determină o variație a concentrației, variație considerată ca proces termodinamic ireversibil, care are loc fără disipare de energie.

pe de altă parte compoziția din interiorul unei particule de fluid în mișcare se modifică prin procese moleculare de transport al masei, având tendința de omogenizare a concentrației. Tendința de omogenizare a concentrației în acest tip de variație a concentrației în exteriorul fiecărui volum elementar de fluid este cunoscută sub numele de difuzie. Difuzia este un proces ireversibil. Impreună cu transportul de căldură și vâscozitate, este una din cauzele disipării energiei într-un amestec de gaze.

Dacă se noteaza cu densitatea totală a amestecului de fluid, atunci ecuația de continuitate pentru masa totală de fluid iși păstrează forma cunoscută:

Această ecuatie exprimă faptul că masa totală a amestecului de fluide, într-un volum arbitrar ales, suferă o variație numai prin fluxurile de masă care intră sau părăsesc acest volum prin suprafața sa. Trebuie însă, de la bun început, să evidențiem faptul că noțiunea de viteză trebuie redefinită pentru amestecul de fluide.

Scriind ecuația de continuitate de mai sus, am considerat viteza în sensul definiției, Ecuațiile Navier – Stokes nu se modifică pentru amestec. Aceasta înseamnă că derivata totală dc/dt este identic nulă, adică:

Aceasta ecuație, tinând cont de relația de continuitate, se mai poate scrie:

adică, sub forma unei ecuații de continuitate pentru unul din constituenții amestecului (ci este masa constituentului pe unitatea de volum).

Deoarece difuzia contribuie suplimentar la fluxul constituentului considerat, ea va genera în fluxul total o componentă distinctă de flux, care actionează ca transport de materie, chiar dacă miscarea generală încetează. Fie densitatea acestui flux de difuzie, adică o cantitate din constituentul considerat, care este transportată prin unitatea de suprafață în unitatea de timp.

Variația aceastei cantități de substanță într-un volum dat, este dată de relația diferențială:

Folosind ecuația de continuitate, ultima relație de mai sus se poate rescrie:

5.1. Ecuațiile de difuzie

Concentrațiile, notate cu , de aerosoli sau gaze (aerosolii fiind particule cu diametrul mai mic de 20 de microni) provenite de la o sursă continuă cu o înălțime efectivă de emisie H, calculate pe direcțiile x, y și z sunt date de ecuația: unde: (x, y, z; H) = concentrația

H = înălțimea centrului penei, considerată în momentul când aceasta ajunge să aibă un volum substanțial, și este suma dintre înălțimea fizică a coșului de emisie, notată cu h și H care reprezintă înălțarea penei.

Se va ține cont de următoarele ipoteze:

pana va dezvolta o dispersie cu o distribuție de tip Gauss atât în plan orizontal cât și în cel vertical. Deviațiile distribuțiilor de concentrații, pe orizontală și verticală (Ϭysi Ϭz ) sunt considerate a fi standard;

u = viteza medie a vântului care afectează pana;

emisia de poluant este uniformă, debitul ei fiind Q;

se presupune că solul reflectă în totalitate pana, deci nu intervin depuneri sau reacții cu solul;

y, z= coeficenți de difuzie după axele y respectiv z.

Ecuația anterioară este valabilă în cazul în care difuzia pe direcția deplasării penei poate fi neglijată, deci dacă nu există dispersie pe direcția x.

Ipoteza de mai sus este valabilă în cazurile în care:

emisia se face continuu;

durata emisiei este mai mare sau egală cu timpul de deplasare (x/u) de la sursă până la zona de interes.

Figura 5.1.1.Distributia gaussiană bidimensională a penelor de poluant generate de surse aflate în originea sistemului de coordonate

Concentratia la sol:

Pentru concentrațiile calculate la nivelul solului (z=0), ecuația se simplifică

Pentru explicarea metodei reflexiei penelor de poluant pe sol, care au ca rezultat suprapunerea distributiei concentratiei unei surse la cota H, cu distributia concentratiei unei surse la cota H, cu distributia concentratiei generate de o sursa fictiva aflata la cota -H (sub nivelul solului).

Acolo unde concentrația se calculează de-a lungul axei penei (y=0) intervin noi schimbări (simplificări):

Pentru o sursă la nivelul solului, fără a exista o înălțare a penei (H=0):

Efectele stabilității. Clase de stabilitate.Valorile lui y și z variază cu:

gradul de turbulență al atmosferei;

înălțimea în raport cu suprafața;

“rugozitatea” (neregularitățile) suprafeței;

eșantionul de timp pe care se va estima concentrația;

viteza vântului;

distanța față de sursă.

Valorile lui y și z sunt estimate conform stabilității atmosferei, care la rândul ei este determinată în funcție de viteza vântului (măsurată la o înălțime de aproximativ 10 metri), de intensitatea radiației solare (în timpul zilei) și de nebulozitate (pe timpul nopții).

Clasele de stabilitate, șase la număr, sunt date în tabelul următor:

Tabel 5.1.Clase de stabilitate

Clasa F este cea mai stabilă iar A cea mai instabilă.

Perioada de noapte este considerată ca fiind intervalul de timp care să conțină noaptea propiu-zisă plus o oră înainte de asfințitul și o oră după răsăritul soarelui.

De reținut că, pentru condiții de vreme închisă, indiferent de viteza vântului, se poate folosi clasa neutră D.

Radiația solară “puternică” corespunde cazului cu cer senin și altitudinea solară mai mare de 60. Insolație slabă înseamnă cer senin și altitudine solară între 15 și 35.

Nebulozitatea va diminua radiația solară receptată și va influența, împreună cu altitudinea solară, încadrarea în cele trei categorii (P=Puternică, M=Moderată, S=Slabă).

Radiația solară, puternică în zilele cu cer senin, va fi moderată în perioadele cu nori sparți, aflați la o înălțime medie și care asigură o nebulozitate de 5/8 până la 7/8, și chiar slabă în cazul norilor sparți și joși.

Aceste metode vor furniza informații reprezentative asupra stabilității deasupra unui câmp deschis sau asupra zonelor rurale, dar sunt irelevante pentru zonele urbane. Diferența se datorează efectului de insulă de căldură și faptului că zonele urbane prezintă o “rugozitate” mult mai mare a suprafeței. Diferența se accentuează în nopțile calde și senine, când, spre deosebire de zona rurală unde există stabilitate, în zonele urbane se manifestă condiții de ușoară instabilitate sau de aproape neutralitate până la o înălțime de câteva ori mai mare decât media înălțimii construcțiilor, cu un strat de stabilitate deasupra ei.

5.2 Dispersia în condiții de extremă stabilitate

În situațiile cu vânt slab și cer senin noaptea, având surse la nivelul solului neinfluențate de neregularități topografice, apar de obicei schimbări frecvente ale direcției vântului, ceea ce servește la împrăștierea orizontală a volumului de dispersat.

Pentru surse de înălțime, în aceste condiții de extremă stabilitate, nu se vor semnala concentrații semnificative la nivelul solului până la schimbarea condițiilor de stabilitate.

Existența unui strat stabil deasupra unuia instabil va avea ca rezultat restricționarea difuziei verticale. Evaluarea dispersiei va fi modificată în acest caz, introducând în calcule un termen nou, acesta fiind înălțimea de bază a stratului stabil, notat cu L.

5.3 Înălțimea efectivă de emisie

În majoritatea problemelor, se dovedește a fi necesară estimarea înălțimii efective a coșului, notat cu H, înălțimea la care pana poate fi considerată ca întinsă.Rareori această înălțime corespunde celei fizice a coșului, h.

Dacă pana este antrenată în turbulența stârnită de coș sau de clădirile învecinate, ea va fi proiectată rapid spre sol datorită fenomenului de trenă aerodinamică.

Dacă pana este emisă într-o zonă lipsită de turbulențe, atunci înălțarea ei va fi influențată de o serie de factori de emisie și meteorologici.

1) Factorii de emisie reprezentați de datele de proiect ale coșului:

a) diametrul interior dc al coșului la cota maximă construită (hc);

b) viteza Vc de evacuare a gazelor din coș;

c) emisia de căldură pe unitatea de timp Qh(cal/sec), care deseori este înlocuită cu temperatura Tc de evacuare a gazelor, folosită în calculul termenului de flotabilitate a penei.

2) Factorii meteorologici (condițiile din atmosferă).

a) stabilitatea termică a atmosferei;

b) temperatura aerului Ta;

c) viteza și direcția vântului;

d) gradientul vertical al vitezei ;

e) intensitatea turbulenței (caracteristica pentru procesul de amestec poluant – aer);

f) gradul și etajele de nebulozitate.

3) Natura fizico – chimică a efluentilor emisi în atmosferă:

a) mărimea particulelor (particulele cu diametrul > 20 sunt rapid depuse la sol);

b) densitatea gazului;

c) vâscozitate, capacitate calorică și călduri specifice ale efluen-ților;

d) potențialul reactivității chimice și fotochimice a efluenților.

Marea majoritate a ecuațiilor formulate pentru calculul înălțimii efective de emisie sunt semiempirice.

Ecuația lui Holland a fost desprinsă dintr-un set de date experimentale provenite prin studierea unor coșuri mai mari în diametru (între 1,7 și 4,3 m) și cu temperaturi mai mari ale efluentului (între 82 și 204 C). Ecuația poate subestima uneori înălțimea efectivă de emisie dar prezintă de obicei un ușor grad de siguranță. Ea este:

H = înălțarea penei deasupra coșului [m];

Vs= viteza gazului [m/s];

d= diametrul interior al coșului [m];

p= presiunea atmosferică [mb];

Ts= temperatura gazului [K];

Ta= temperatura aerului [K];

2,6810-3 este o constantă care se măsoară în [mb-1 m-1].

Calitatea aerului în așezările umane se determină prin măsurarea concentrațiilor medii orare, zilnice sau lunare ale diferiților poluanți și compararea acestora cu valorile limită sau după caz cu valorile maxime admisibile prevăzute în actele normative în vigoare (Ordinul MAPM 592/2002) pentru poluanții relevant, respectiv STAS 12574-87 pentru restul poluanților reglementați.

Concentrația emisiilor de poluanți în aerul ambient poate varia în funcție de condițiile meteorologice, mai mult sau mai puțin favorabile unei bune dispersii.

Supravegherea calității aerului la nivel local se face prin măsurarea:

Poluanților comuni (SO2,NO2,NH3,pulberi în suspensie) în toate rețelele locale din zonele industriale și/ sau urbane

Poluanților specifici (HCl,fenoli,aldehide,Cl2,metale grele:Pb,Cd), în zone industriale în funcție de activitate [14]

5.4 Controlul emisiilor în aer

Controlul în aer al emisiilor care provin din procese industriale, au început să fie realizate prin intermediul unor reglementări specific importante.

Reglementările mai vechi au solicitat în general controlul emisiilor de noxe generate de un numar limitat de procese de producție din sectorul industriei. Autoritățile din domeniu au avut datoria de a verifica aplicarea “celor mai bune practici pentru împiedicarea scăpărilor de noxe sau de gaze agresive și pentru transformarea acestor emisii de gaze în emisii inofensive si nedaunatoare”.[17]

Evoluțiile legislative recente s-au axat în special pe acele activități industriale care generează volume mari de poluanți sau cu toxicitate ridicată. Au fost stabilite valori limită pentru emisiile generate de:

Marile centrale cu ardere

Incineratoarele de deșeuri periculoase

Incineratoare de deșeuri orășenești

Incineratoare de deșeuri

În situații adecvate, astfel de limite nu se aplică autorizațiile IPPC.

5.5. Fenomenul de sinergism al agenților poluanți ai atmosferei

Agenți poluanți ai atmosferei reacționează prin mecanisme complexe: reacții fotochimice, de polimerizare, de oxidoreducere, de hidroliză sau cataliză, agravându-se astfel efectul factorilor de mediu.

Efecte:

Compuși ai azotului:

NO + O2 => NO3 ; NO3 + O2 => NO2 + O3

NO2 => NO + O ( în atmosferă superioară)

Compuși ai sulfului:

SO2 + H2O => H2SO3

H2SO3 + H2O => H2SO4 +H2

SO2 +2NO2+ 2H2O => H2SO3 + HNO2

5.6. Considerații generale privind modelarea matematică a transportului și dispersiei poluanților în atmosferă

Cu toate progresele făcute, modelarea matematică a proceselor de transport și dispersie a poluanților în atmosferă prezintă dificultăți mari, cel puțin din două puncte de vedere:

Cunoașterea limitată a fizicii fenomenului, a proceselor chimice, biologice care se produc în timpul transportului și dispersiei poluanților;

Modelele matematice sunt în curs de perfecționare întrucât fenomenul pe care-l descriu este în curs de cunoaștere și întrucât aparatul matematic folosit încă prezintă schematizări atât de complexe care au loc în evoluția poluanților.

Modelele matematice sunt axate pe transportul și dispersia poluanților:

Pe distanțe mici și medii la microscară și pe distanțe mari, la macroscara, care însoțește fenomenele meteorologice în troposferă și care țin seama de câmpul baric, de accelerația Coriolis, de câmpul termic, de circulația norilor.[16]

5.7 Descrierea fenomenologică a poluării atmosferei din surse punctiforme

Exemple de surse punctiforme evidente și des întâlnite sunt coșurile de evacuare a gazelor și a pulberilor de la diverse industrii.

Figura 5.7.1 Schemă de evoluție a unei pene de la o sursă punctiformă

În figura de mai sus se prezintă schema generală a unei pene de gaze de pulberi evacuate dintr-un coș al unei uzine sau centrale termice. Pe această schemă sunt ilustrate convențional și diverse fenomene care se pot produce în timpul evoluției penei: sedimentare uscată de particule solide, depuneri umede, reacții chimice și fotochimice, reflexia de sol, reflexia de plafonul de nori sau de straturile superioare ale atmosferei, formare de vârtejuri ș.a.

În cazul unei atmosfere stabile în care viteza u a vântului este mai mică de 2 – 3 ori decât viteza w a gazelor la iesirea din coș ( w> 2u ), pana parcurge patru faze:

Faza 1, la ieșirea din coș, caracterizată printr-o poziție aproape verticală, cu viteza w de ieșire a gazelor mai mare decât viteza vântului u ( w> u ) și se consideră pana în secțiunea în care w ≅ u

Faza 2, în care pana se curbează pâna în secțiunea unde viteza w este apropiată ca mărime de mărimea pulsațiilor vitezei vântului, u’ ( w≅u’ );

Faza 3, în care densitatea gazului emis este apropiată de densitatea aerului, cantitatea de mișcare proprie a gazelor este practic nulă; procesul de transport prin convecție este dominant pe direcția vântului; procesul de difuzie turbulent este important în secțiunea transversală, pe direcția vântului fiind în general neglijabil în raport cu convecția.

Figura 5.7.2.Evoluția penei de gaze de emisie pentru h>0.

Faza 4, cand pana atinge suprafața pământului și se reflectă de aceasta , parțial sau total, ca și cum ar fi doua surse de poluare, una reala cu H= h + H>0 și alta imaginară cu H <0.

În cazul în care viteza gazelor la iesirea din coș, w, este mai mică decât viteza vântului, sau 1,5 ori viteza vântului, ( w< 1,5 u) supraînălțarea penei H poate fi negativă. Procesele fizico-chimice, chimice, biochimice care au loc în interiorul penei sunt greu de identificat, unele pot fi sinergice.

Figura 5.7.3.Evoluția penei de gaze emise când h<0.

Problemele importante care se pun constau în cunoașterea cantitativă a acestor procese, cunoașterea concentrației poluantului la suprafața solului, în zona de contact cu solul sau cu obiectele de la suprafața pământului, deci la imisie, în funcție de concentrația poluantului la sursă, deci la emisie, în funcție de viteza de ieșire a gazelor, de diametrul și înălțimea coșului, în funcție de viteza vântului, de diferența de temperatură și densitate între poluant și aerul atmosferic, în funcție de gradientul termic al aerului și de alți factori. De asemenea, prezintă importanță cunoașterea depunerilor de poluant la suprafața solului, pe învelișul vegetal, pe suprafața apelor, cunoașterea efectului (impactului) pe care îl au poluanții transportați asupra factorilor de mediu,asupra omului și asupra ecosistemelor, cunoașterea distanței până la care concentrația poluanților este dăunătoare.

Sursa de poluare din figura de mai sus poate fi considerată ca sursă punctiformă virtual, situată undeva în punctul S, foarte aproape de verticala coșului, dar la o înăltime față de secțiunea de iesire din coș, numită supraînălțare, respectiv la înălțimea H= h + față de suprafața solului, unde h este înălțimea coșului.

Fenomenul de transport și dispersie a poluanților se complică, dacă în calea penei apar forme de relief care modifică forma acesteia.

5.8. Elemente de aerostatica și aerodinamică care intervin la studiul și evoluția penei

5.8.1 Elemente de aerostatică

În condiții de echilibru al atmosferei există relații de stare între presiune, temperatură și masa specifică a aerului. Se știe că, în hidrostatica aplicată la gaze compresibile, există relația de echilibru:

unde n este coeficientul de politropie, care la procese adiabatice devine:

unde cp este caldura specifică la presiune constantă, iar cv căldura specifică la volum constant. Pentru aer k = 1,4.

Pe de altă parte se poate apela la ecuatia de stare Clapeyron

sau

cu R0= 29,27 m/grad

p – presiunea, z – cota fața de un plan de referință, de exemplu, planul la suprafața pământului, – masa specifică, T – temperatura absolută ( în grade Kelvin) , R constanta universală:

R = 288,41 m2/s2grad – pentru aer cu umiditate medie;

R = 287,14m2/s2grad – pentru aer uscat

Din relațiile de mai sus, ținând seama de relația , unde g este accelerația gravitațională, obținem:

sau

sau încă

Pentru n= 1,2÷1,4 rezultă:

unde reprezintă gradientul de echilibru adiabatic (politropic). De cele mai multe ori se adopta

Dacă distribuția temperaturii pe verticală corespunde cazului în care solul este mai rece decât aerul, iar temperatura crește pe verticală, este evident ca stratificația atmosferei este stabilă: aerul mai rece este mai greu, aerul mai cald este mai ușor și orice perturbație se amortizează.

Dacă, dimpotrivă, temperatura aerului este mai mică decat temperatura solului, gradientul de temperatură , apar curenți de aer cald ascendenți de la sol spre exterior și curenți descendeți de aer rece care vin sa ia locul eliberat de aerul cald.

Figura 3.8.1.Atmosfera instabila

5.8.2Elemente de aerodinamică

Atmosfera este construită din mai multe straturi. Pentru ilustrare în figura de mai jos (fig. 5.) este prezentată variația temperaturii și variația presiunii în funcție de altitudine.

Fenomenul de poluare, ca și fenomenele meteorologice, se produc în principal în troposferă, care, convențional, constituie un inveliș de circa 10km grosime ( 8 – 9 km la ecuator si 13 – 15 km la poli).

Figura 5.8.2.1.Variația temperaturii și variația presiunii în funcție de altitudine.

La contactul cu suprafața pământului, curenții de aer formează stratul limită de suprafață de circa 100 – 200 m și stratul limită de tranziție, deasupra stratului limită de suprafață, cu o grosime de circa 1000m. Este evident ca frecările și schimbul limită nu se mai resimte frecarea de suprafața pământului, vântul se racordează la vântul geografic care are directive paralela cu liniile izobare.

5.9. Distribuția pe verticală a vitezei curenților de aer în stratul limită

Se disting două cazuri:

atmosferă stabilă

atmosferă instabilă

În cazul atmosferei stabile, când stratificația orizontală a aerului se păstrează, distribuția pe verticală a vitezei vântului este de tip logaritmic, având expresia:

Dacă suprafața pământului este netedă și

Dacă suprafața pământului este rugoasă, unde k = 0,4 este constanta lui Karman, u0= este viteza dinamică, – efortul unitar de forfecare, – masa specifică, – coeficient de vâscozitate cinematic a aerului, z0 – rugozitatea, C – constantă. Mărimile u*, z0, C se determină prin măsurători. De exemplu, dacă la înăltimea z1 se măsoară viteza u1, și la înaltimea z2 se măsoara viteza u2, din ecuație obtinem:

de unde rezultă u*și z0.

Este important de reținut faptul că gradientul vitezei pe verticală, atât din relatiile de mai sus are expresia:

În cazul în care gradientul de temperatură diferă de gradientul adiabatic, pentru descrierea distribuției vitezei pe verticală, se folosește lungimea caracteristică L introdusă de Monin și Obukhov.

Unde T este temperatura absolută măsurată în grade Kelvin, cp- căldura specifică la presiune constantă, q – fluxul de căldură pe direcția verticală. Dacă atmosfera este stabilă

fluxul de caldură este negativ (q<0) și L>0. Dacă atmosfera este instabilă , fluxul de căldură este pozitiv (q>0) și L<0.Pentru q=0, L=∞ ceea ce se întamplă la

În cazurile în care, , legea de distribuție pe verticală a vitezei vântului se deduce din relația

unde este o funcție care ține seama de gradientul termic.

Pentru determinarea funcției au fost efectuate măsurători și prelucrări, care în conditii de <0 (L>0)

Figura 5.9.1.Graficul functiei

4.10. Distribuția pe verticală a temperaturii aerului

Se constată că distribuția temperaturii aerului pe verticală este descrisă de relații de tip ca și distribuția pe verticală a vitezei curenților orizontali:

unde

Este temperatura de frecare, iar

Prin analogie cu km ( coeficientul de transfer de cantitate de mișcare) se numește coeficient de transfer de căldură.

La atmosfera neutră se consideră

km= si

La atmosfera stabilă < km, la atmosferă instabilă > km.

Relațiile de mai sus, ca și relațiile de distribuție pe verticală a vitezei curenților orizontali, necesită măsurători pentru a fi utilizate. În lipsă de măsuratori se consideră distributia log – liniară a vitezei si km =, deci , cu valorile approximate 4<<7.

5.11. Criterii de stabilitate a atmosferei

Se poate face o analiză detaliată a stabilității atmosferice. În condiții de calm atmosferic (u0 ) un criteriu de stabilitate termică se poate obține din compararea gradientului termic cu gradientul adiabatic. Avem:

atmosferă instabilă

atmosferă stabilă

Relațiile de mai sus evidențiază apariția și dezvoltarea unor curenți verticali de aer cald (ascendenți) și de aer rece (descendenți), cu formarea unor vârtejuri, deci a unei turbulențe atmosferice de natura termică.

În condiții de vânt, frecarea aerului de suprafața pământului conduce la o distribuție neuniformă a vitezei curentului, cu formarea de vârtejuri de natură mecanică.

Cu cât diferența de viteză a curentului de aer între două puncte pe aceeași verticală e mai mare, cu atât amestecul paturilor de aer este mai puternic și temperatura mai apropiată de cea uniform distribuită. Dar diferența de viteză pe verticală a curenților orizontali este dată de gradientul de viteză. Asadar, turbulența atmosferică la scară mică (microscopic) are două origini:

termică, exprimată prin gradientul de temperatura

mecanică, exprimată prin gradientul de viteză al curentilor

Criteriul de stabilitate a atmosferei în condiții dinamice este dat de numărul Richardson:

Care combină cei doi gradienți de temperatură și de viteză ( T- temperatura absolută, g – accelerație gravitatională).

Când vântul este de intensitate mare gradientul de viteză este de asemenea foarte mare și mărimea de la numitorul raportului are același efect ca apropierea gradientului temperaturii de gradientul adiabatic. Efectul vântului puternic, este din punct de vedere fizic, dublu: pe de o parte, creșterea gradientului vitezei din cauza frecării aerului de suprafața pământului, pe de altă parte intensificarea schimbului de cantitate de mișcare între stratele de aer,ceea ce conduce la uniformizarea temperaturii.

Criteriul Richardson se scrie:

Ri> 0 – atmosfera stabilă;

Ri= 0 – atmosfera neutră;

Ri< 0 – atmosfera instabilă; Întrucat numărul Richardson necesită măsurarea gradientului vitezei, deci a distribuției vitezei pe verticală, se numește numărul Richardson global sub forma:

unde viteza u este masurată doar la cota z.

Se poate completa analiza fenomenologică a dispersiei poluanților în atmosferă luând în considerare și fenomenul de inversiune termică, respectiv de schimbare a semnului gradientului termic pe verticală.

Figura 5.11.1.Condițiile de instabilitate foarte accentuate (Ri<0).

Figura 5.11.2.Condiții de stabilitate aproape neute (Ri=0).

Figura 5.11.3.Condiții de stabilitate accentuate (Ri>0).

Figura 5.11.4.Inversiune termica tip A.

Figura 5.11.5.Inversiune termica tip B.

5.12. Influența reliefului asupra structurii curenților atmosferici și asupra formei penei

Pana de poluant se racordează la structura curenților de aer care se formează în condiții de relief accidentat sau în zona construcțiilor.

Figura 4.12.1.Ocolirea unui munte(deal)

Figura 4.12.2.Denivelare pozitivă

5.13. Modele matematice pentru calculul transportului Și dispersiei poluantilor in atmosfera

Modele matematice pentru studiul dispersiei poluanților pot fi clasificate în trei tipuri:

modele matematice de simulare numerică atât a câmpului vitezelor și presiuni aerului, cât și a evoluției jetului de fluid poluant emis în atmosferă, în condiții concrete de spațiu și timp;

modele aproximative, de tip Gaussian, pentru care câmpul de viteză al vântului este dat de fiecare moment, în regim nepermanent sau permanent de evoluție a jetului de gaze poluante;

modele aproximative, tot de tip Gaussian, dar în valori prelucrate statistic ale câmpului de viteză, ca mărime, directie și frecvență.Se mai numesc și modele statistice.

Aceste din urmă modele matematice conțin ipoteza că media în timp a rezultatelor calculelor privind concentrația poluanților, în orice punct din spațiu, obținută prin simulare numerică a fenomenului dispersiei, este egală cu rezultatul calculului de dispersie a poluanților în condițiile în care datele de intrare sunt mediate în timp, statistic. Funcția de transformare a datelor de intrare în rezultat final este liniară:

Unde x0 este media aritmetică a valorilor xi.

Modele de calcul Pasquill – Gifford – Turner

Este de tip Gaussian cu parametrii estimați pe baza unor elemente mai mult apreciative decât analitice, mai mult calitative decât cantitative. Este vorba în special de stabilirea unor clase de stabilitate a atmosferei și a unor clase de variație a vitezei vântului cu ajutorul carora se calculează parametrii și

Modelul propus de Pasqiull în anul 1959 are la baza considerarea a șase clase de stabilitate a atmosferei si anume:

clasa A – foarte instabilă

clasa B – instabilă

clasa C – ușor instabilă

clasa D – neutră

clasa E – stabilă

clasa F – foarte stabilă

Clasele de stabilitate au fost caracterizate de doi parametri: viteza vântului u ,la înăltimea de 10 m deasupra terenului și radiația solară în timpul zilei sau nebulozitatea, respectiv gradul de acoperire cu nori a cerului în timpul nopții.

Pentru cer senin radiația solară incidentă este puternică dacă înălțimea soarelui este 600; este moderată dacă este cuprinsă între 600 și 350 și este slabă dacă această înălțime este cuprinsă între 350 și 150.

Noaptea înseamnă o ora înainte de apusul soarelui până la o oră după răsăritul soarelui. Hino a cuantificat radiația solară incidentă puternică, moderată și slabă în felul următor:

Puternică RD > 600w/m2 sau RD >123 cal/m2s

Moderată 300 w/m2 <RD <123cal/m2s

Slabă RD < 300w/m2 sau RD < 61 cal/m2s

Clase de stabilitate a atmosferei după Pasquill

Tabel 5.13.1Clase de stabilitate dupa Pasquill

Gifford a propus următorul tabel:

Tabel 5.13.2. Clase de stabilitate a atmosferei

6 PROGRAMUL DE CALCUL UTILIZAT PENTRU MODELAREA DISPERSIEI POLUANTILOR IN ATMOSFERA

6.1. Principalele variabilele folosite în program

Parametri de emisie pentru fiecare sursă:

Q – 14,130,40 kg/h – debitul de poluant (vector);

h – 10 m -înălțimea coșului (vector);

t – 120 ° C temperatura efluentului (vector);

VS – 1.36, 12.82, 3.96 m/s -viteza efluentului (vector);

d – 0.6 m-diametrul coșului (vector);

Parametrii meteorologici

an, luna, zi – data de calcul;

ora – ora de calcul;

Radiație – tipul radiației solare (puternică, moderată, slabă);

TOA – 20°C – temperatura mediului ambiant;

.u – 0.1, 2.4 m/s – viteza vântului;

.NEBUL – 5*1/8 – nebulozitatea;

p – 1 barr -presiunea atmosferică;

Parametrii rețelei și setări grafic

nSurse – numărul de surse (max. 10);

EditR – matrice ce conține poziția fiecărei surse în rețea;

Definirea curbelor de izoconcentrație

FIzoConc.ListBox1 – listă ce conține titlul și valoarea fiecărei curbe;

culori – vector ce conține culoarea fiecărei curbe;

ElLeg – vector ce conține fiecare element al legendei;

Variabile de calcul

Hef – vector ce conține înălțimea efectivă a fiecărei surse

cl – păstrează clasa de stabilitate necesară calculului;

IzoConc – vector de serii de tip TPointSeries necesar graficului;

pps – numărul de pixeli pe segment (necesar la afișarea la scară a graficului);

dx, dy – numărul de pixeli pe axa Ox, respectiv Oy;

6.2 Instrucțiuni de utilizare

În figura de mai jos este prezentată fereastra de dialog principală cu ajutorul căreia utilizatorul va putea executa următoarele operații:

Figura 6.2.1.Fereastra principală a programului

Salvarea (înregistrarea pe disc) a setărilor graficului

În cazul în care se dorește păstrarea setărilor graficului se selectează opțiunea Fișiere|Salvare setări grafic din meniu și se alege numele și calea fișierului ce va păstra setările respective.

Deschiderea unui fișier cu parametri deja existent

Această operație se execută în orice moment prin selectarea opțiunii din meniul Fișiere|Deschidere fișier parametri de intrare din meniu sau prin apăsarea celui de-al doilea buton din bara de butoane (o carte deschisă). Datele existente în fișierul ales vor fi citite și afișate în fereastra de setare a parametrilor de intrare.

Execuția acestei operații va conduce la pierderea tuturor parametrilor de intrare ce nu au fost salvați în prealabil precum și la necesitatea reluării calculelor de concentrații

Tipărirea la imprimantă a graficelor

Se va putea face numai dacă pe ecran există un grafic, acesta fiind trimis la imprimantă sub format A4 prin selectarea opțiunii Fișiere|Tipărire grafic din meniu sau prin apăsarea butonului corespunzător (imprimanta).

Introducerea datelor necesare setării graficului

Datele vor fi introduse în fereastra de dialog (figura 16) afișată la apăsarea butonului ce simbolizează un grafic sau la selectarea opțiunii Setări|Grafic din meniu în orice moment.

Figura 6.2.2.Forma grafic

Se va introduce distanța și numărul punctelor de calcul (după direcția x) pentru care se va calcula dispersia poluanților. La apăsarea butonului de afișare a graficului calculatorul va obține distanța după direcția y în funcție de dimensiunile ecranului (astfel încât graficul să fie la scară) și le va afișa în partea dreaptă a ferestrei în grupul Date rezultate prin calcul.

Prin apăsarea butonul Renunțare se revine în meniul principal fără a se păstra pe ecran graficul și setările făcute. Continuitatea programului se realizează prin apăsarea butonului Ok, acesta activând butonul de calcul.

Introducerea datelor necesare rulării programului

Se realizează prin selectarea opțiunii Setări|Parametrii de intrare din meniul principal sau prin apăsarea butonului corespunzător din bara de butoane (al cincilea). La apariția pe ecran a ferestrei din figura 17.a și în cazul în care graficul nu a fost afișat pe ecran, va fi activat numai butonul Ok.

Butoanele Renunțare și Salvare date se vor activa numai la apariția unei modificări a datelor din fereastră, iar butoanele Figurare surse și Figurare rețea vor fi activate numai dacă pe ecran este afișat graficul.

În fereastra de introducere a datelor, mărimile utilizate în procedurile de calcul și a căror valori trebuiesc date de utilizator au fost grupate pe pagini și categorii pentru a evita eventualele confuzii datorate asemănării notațiilor și numărului mare al acestora. Astfel fereastra de dialog Parametri de intrare conține 3 pagini: Parametri rețelei, Parametri meteorologici și Parametri de emisie ce se derulează prin apăsarea tab-urilor corespunzătoare.

În pagina Parametri rețelei, inițial afișată (vezi figura 4.2.3), se va seta numărul de surse prin folosirea controlului UpDown1 din dreapta, iar poziția surselor se va modifica prin scrierea în controalele Edit a distanțelor dorite. Trecerea de la o sursă la alta se va face la fel prin folosirea controlului UpDown2 din dreapta grupului (acest control este afișat numai dacă numărul de surse este mai mare decât 1).

Figura 6.2.3. Parametrii rețelei

Pagina Parametri meteorologici (figura 4.2.4) va conține date referitoare la data de calcul și condițiile atmosferice existente la acea dată. Sunt cerute radiația solară și nebulozitatea atmosferei pentru ca în funcție de ora de calcul (noaptea este între orele 20 și 06) programul să stabilească clasa de stabilitate.

Figura 6.2.4.Parametrii meteorologici

În figura 4.2.5este prezentată ultima pagină a ferestrei de dialog Parametrii de intrare în care de asemenea vor fi specificați parametrii de emisie ai fiecărei surse. Trecerea de la o sursă la alta se face cu ajutorul controlului ScrollBar de pe marginea superioară a grupului (este vizibil numai dacă numărul surselor este mai mare decât 1).

Figura 6.2.5Parametrii de emisie

La apăsarea butonului Salvare date, toți parametri specificați în fereastra Parametrii de intrare vor fi salvați pe disc sub forma unui fișier (*.fpi), al cărui nume va fi dat în fereastra dialog ulterioară. Datele vor putea fi reîncărcate în memorie în orice moment la cerere (vezi deschiderea unui fișier cu parametri deja existent).

În cazul în care utilizatorul se răsgândește asupra modificării parametrilor, el poate apăsa butonul Renunțare, în acest caz toate datele fiind restaurate din memorie iar fereastra se va închide. Restaurarea datelor din memorie nu este posibilă și la deschiderea unui fișier de date. Aici se vor reface numai datele modificate ulterior deschiderii fișierului.

La apăsarea butonului Ok toate datele existente în fereastra dialog Parametri de intrare vor fi trecute în variabilele de lucru din memorie și programul poate rula în continuare cu efectuarea calculelor (sau setarea graficului după caz).

Efectuarea calculelor și trasarea curbelor de izoconcentrație

Aceste operații se vor putea face numai dacă pe ecran este afișat graficul respectiv pentru trasare numai dacă a fost efectuat calculul. Deci efectuarea calculelor va determina activarea butonului de trasare a izoconcentrațiilor, a butonului de citire a graficului și dezactivarea butonului de calcul (opțiunii Execuție program|Calculul concentrațiilor din meniul principal).

Efectuarea calculelor va trebui refăcută de fiecare dată când apare vreo modificare în setările graficului sau în parametrii de intrare. La apariția unei modificări eventualele curbe de pe grafic vor dispare, se vor dezactiva butoanele de legendă, de citire a graficului respectiv de trasare a izoconcentrațiilor și se va activa butonul de calcul (aceleași operații se vor efectua și în meniuri).

Trasarea curbelor de izoconcentrație se execută la selectarea opțiunii Execuție program|Trasare izoconcentrații din meniul principal sau apăsarea butonului corespunzător din bara de butoane. Fereastra de dialog ce va apare este prezentată în figura 4.2.6.

Figura 6.2.6. Fereastra de creare a listei de izohipse

Izo 0.01,a fost adăugată în listă prin apăsarea butonului Adaugă izoconcentrație și scrierea valorii curbei în controlul Edit al ferestrei apărute.Procedeul se poate aplica în continuare pentru toate curbele dorite.

La selectarea unui element al listei, în colțul din dreapta jos va apare un grup de controale ce conțin caracteristicile curbei respective (valoarea și culoarea). Valoarea nu se poate modifica decât prin ștergerea curbei prin apăsarea pe butonul Ștergere izoconcentrație (în momentul în care aceasta a fost în prealabil selectată din listă). Culoarea se poate modifica prin apăsarea pe butonul colorat (implicit roșu), moment în care pe ecran va apare o paletă de culori din care utilizatorul își va alege culoarea preferată. Se poate de asemenea ca utilizatorul să își creeze propria culoare urmând indicațiile din fereastră.

Pentru rapiditatea ștergerii curbelor a fost prevăzut butonul Ștergere listă izoconcentrații care va distruge toate curbele create și va goli lista.

La apăsarea butonului Ok programul va începe baleierea graficului, va afișa curbele cerute și va activa opțiunea Ajutor|Legendă – Izoconcentrații din meniul principal. În continuare utilizatorul poate cere trimiterea la imprimantă a graficului cu curbele de izoconcentrație trasate sau poate continua cu modificarea graficului pe o altă distanță sau modificarea parametrilor de intrare.

Vizualizarea unei ferestre ce conține poziția și concentrația fiecărui punct din grafic indicat

Dacă pe ecran este afișat graficul și dacă calculele au fost efectuate utilizatorul are opținea de a citi datele din grafic prin intermediul mouse-ului și a unei ferestre ce va apare la selectarea opțiunii Ajutor|Citirea datelor din grafic sau la apăsarea butonului ce simbolizează o lupă. Datele citite din grafic vor fi afișate în această fereastră și se vor modifica prin simpla mișcare a mouse-ului pe grafic.

Fereastra se închide la cerere prin apăsarea butonului Ok.

Afișarea legendei graficului

Această opțiune este disponibilă în meniu (Ajutor|Legenda – Izoconcentrații) sau bara de butoane (ultimul buton) în momentul în care curbele au fost deja trasate pe grafic. Elementele legendei sunt curbele definite în fereastra de creare a listei de izoconcentrații.

Legenda se închide la cerere prin apăsarea butonului Ok.

Obținerea unor informații despre program

Se poate face în orice moment prin selectarea opțiunii Ajutor|Despre program

6.3 Concluzii referitoare la programul de calcul

Programul realizat prezintă următoarele avantaje:

introducerea ușoară a parametrilor de intrare;

oferă posibilitatea modificării și salvării acestora;

realizează în mod automat scările de reprezentare determinând distanța pe Oy în funcție de distanța cerută pe Ox;

permite citirea concentrațiilor de poluant în orice punct indicând atât valoarea concentrației cât și poziția respectivului punct în spațiu;

permite trasarea curbelor de izoconcentrație provenite de la una sau mai multe surse (numărul maxim de surse este 10);

curbele trasate se pot identifica prin afișarea legendei;

realizează calculul concentrațiilor având în vedere 3 variante: normală, cu fumigație, cu blocarea penei;

Programul este structurat sub formă modulară astfel încât permite perfecționarea ulterioară și în acest sens sunt în vedere adăugarea unor noi module de calcul pentru calculul concentrațiilor de poluant provenite din surse liniare, surse intermitente, surse aflate în mișcare.

7. Studiu de caz – STAtIE DE BETON MOBILĂ AMPLASATĂ PE soseua de centura Târgu Mures

Pe baza modelului prezentat anterior, s-a realizat o evaluare a emisiilor de gaze pentru statia de betoane.

Orașul Târgu Mureș este situat în centrul Transilvaniei fiind traversat de râul Mureș și înconjurat de dealuri împădurite.

Drumul de ocolire, la est de Târgu Mureș, realizează legătura dintre DN13 (E60), km 164+850, și DN15, km 83+900 intrare în localitatea Ernei, la nord de Târgu Mureș, aproximativ la 4 km de oraș.

Târgu Mureș este traversat de drumuri naționale și județene în direcția nord- sud, DN13, DN15, DJ152 și în direcția est – nord, DN15, DJ135.

Statia de beton mobilă prezentată este situată pe soseaua de centură a orasului Târgu-Mures în partea de E-NE la distanta de 1 km de la începerea soselei în partea de S.

Figura 7.1.Plan de incadrare in zona a soselei de centură Târgu Mures

Din punct de vedere climatic teritoriul județului Târgu-Mures se încadrează în sectorul de climă continental-moderată, zona studiată aparținând ținutului climatic al podișului Transilvaniei.

Numărul mediu anual al zilelor de îngheț totalizează 130-140 zile în zona premontană.

Aerul atmosferic in judetul Mures este influentat într-o cantitate moderată de emisiile din activitatile economico – sociale. Sursele antropice de emisie în atmosfera cu potential semnificativ sunt amplasate în Târgu Mures, Iernut si Târnaveni în timp ce în zone, precum Reghin, Ludus , Sovata, sursele de emisie antropice nu produc poluare puternica.[18]

A.) O scurtă caracterizare a surselor fixe se prezintă astfel:

surse industriale:

industria chimică, industria de prelucrare a lemnului, producerea

materialelor de constructie;

industria energetica si termoficare care utilizeaza gazele naturale în scopuri

energetice si deseuri din lemn în centralele termice recent modernizate;

stocarea si distributia carburantilor;

utilizarea solventilor;

prepararea hranei si consumul casnic;

În general, combustibilul folosit este gazul natural (peste 97 %). Sursele de emisie degajate in atmosfera din lucrari agricole sunt slabe dar nu de neglijat. Deseurile depozitate reprezintă surse moderate de emisie în atmosfera, fiind amplasate la distante mari fata de localitati.

B). Surse mobile:

traficul rutier care se desfasoară in mare parte pe DN13 (E60) și DN15 (Târgu Mures – Toplita) si care tranziteaza localitatile urbane si rurale din judetului Mures;

traficul feroviar este reprezentat slab.[19]

Degajarile de praf, ce survin în timpul lucrarilor de construcției, sunt in mare parte din cauza lucrărilor de excavare și mutare a pământului, dar si lucrari necesere precum taluzarea sau nivelarea degaja mult praf in aerul de la nivelul solului

Degajările de praf în aer difera adesea substanțial de la o zi la alta, depinzând de nivelul activității, de specificul operațiilor și de condițiile meteorologice.[20]

Statia de betoane de ciment

Se face ipoteza ca, în perioadele de vârf, vor fi necesare cantităti de betoane de 400 m3/zi respectiv cca 1000 t/zi, ce se vor realiza în baza de productie.

Pentru evaluarea emisiilor s-a folosit metodologia US-EPA/AP42/1999.

Factorii de emisie si debitele masice ale emisiilor de particule sunt prezentate în tabelul 5.1.Emisiile cuprind în principal praf de ciment si particule fine din agregatele minerale utilizate la prepararea betonului.

Tabelul 7.1. Emisiile statiei/statiilor de betoane

*Particule materiale.

Debitele masice de PM rezultate din statia de preparare a betonului au valori importante si pot depăsi concentratiile maxime admise în aer de 0,5 mg/mc în condiții meteorologice nefavorabile si concentrări a activitatii pe arii restranse.

Activitatea în stația/stațiile de preparare a betoanelor de ciment

Poluarea specifica a acestei activități include exclusiv producerea betonului.Sunt luate în considerare emisiile de PM, inclusiv ciment, de la producerea betonului. Nu se iau în considerare emisiile de particule rezultate prin eroziunea vântului din depozitele de agregate,din circulația mijloacelor de transport și activitatea utilajelor, aceste emisii fiind apreciate global în cadrul activității utilajelor de construcție și mijloacelor de transport.

Cele mai mari emisii sunt de particule de ciment ce pot fi aprox. 1 ‰ din cantitatea mutata. Prevederea cu filtre textile in silozurile folosite pentru stocarea cimentului si verificarea atenta a instalatiei pneumatice de incarcare/descarcare a cimentului sunt obligatorii pentru micsorarea pierderilor de pulberi de ciment si încadrarea concentratiilor de PM în aer în standardul legal.

Pentru viteza efluentului s-au obtinut urmatoarele rezultate:

ρ=

W=

ρaer=1,023kg/m3

Qm=14==1,4kg/h

Qm=1,4kg/h*1000/3600=0,39 g/s

ρ====0,38

Q=A*V

V=

ρ=densitate[kg/3]

m=masa[kg]

w=volum[]

Qm=debit masic[g/s]

Qw=debit volumic[g/s]

V=viteza efluentului[m/s]

Q=debit[kg/zi]

A=aria[]

Tabelul 7.2.Calcului vitezei efluentului

Pe baza acestor valori s-a realizat modelarea transportului poluantilor în atmosferă, poluanti care sunt presupusi a fi evacuati în atmosferă prin cotul de fum a statiei de betoane. Statia de betoane este amplasată pe soseaua de centura Târgu Mures Est.

În figurile urmatoare se prezintă distributia spatială izoconcentratiilor celor 2 variante analizate, pe directia predominantă a vântului (NV-SE).

Figura 7.2.Evolutia spatiala a concentratiei de PM în varianta de radiatie solară puternică, fumigatie si viteza vântului de 0.1 m/s.

Figura 7.3.Evolutia spatiala a concentatiei de PM în varianta de radiatie solară puternică cu variantă cu fumigatie si viteza de radiatie solară puternică si viteza vântului de 0,1 m/s

Figura 7.4.Evolutia spatiala a concentratiei de PM în varianta de radiatie solară slabe cu variantă cu fumigatie si viteza vântului de 0,1

Figura 7.5.Evolutia spatiala a concentratiei de PM în varianta de radiatie solară slabă cu variană normală si viteza vântului de 0,1

Figura 7.6. Evolutia spatiala a concentratiei de PM în varianta de radiatie solară slabă cu variantă cu fumigatie si viteza vântului de 2,4

Figura 7.7.Evolutia spatiala a concentratiei de PM în variant de radiatie solară slabă cu varianta normală si viteza vântului de 2,4

Figura 7.8. Evolutia spatiala a concentratiei de PM în varianta de radiatie solară puternică cu variantă cu fumigatie si viteza vântului de 2,4

Figura 7.9. Evolutia spatiala a concentratiei de PM în varianta de radiatie solară puternică cu variantă cu normală si viteza vântului de 2,4

În concluzie, conform figurilor, în tabelul urmator se prezintă concentratiile maxime de poluanti, distantele la care se depun penele de poluant, lungimile si lătimile acestora.

Tabel 7.3. Concentratiile maxime atinse de poluanti, pe directia predominantă a vântului

8. CONCLUZII

Obiectivul prezentei lucrări de licență a fost aplicarea unui program ce poate fi adaptat pentru orice poluant rezultat în atmosfera din diferite activitati și cu ajutorul căruia să poată fi modelat comportamentul acestuia pentru a permite evaluarea nivelului de poluare al atmosferei si solului.

Cu ajutorul programului s-au obsinut concentratiile maxime de poluanti emisi de statia de beton, precum si evolutia acestora reprezentate sub forma hărtilor de distributie spatiala a izoconcentratiilor.

Este important să fie identificate zonele afectate de norul de poluant, în diferite ipoteze, pentru stabilirea masurilor ce trebuie adoptate în vederea reducerii concentratiilor aferente

9. BIBLIOGRAFIE

[1] http://mmut.mec.upt.ro

[2] http://www.calitateaer.ro/parametrii.php

[3] LEGEA nr. 104 din 15 iunie 2011, www.calitateaer.ro/parametri.php‎

[4] LEGEA nr. 104 din 15 iunie 2011, www.calitateaer.ro/parametri.php‎

[5] LEGEA nr. 104 din 15 iunie 2011, www.calitateaer.ro/parametri.php‎

[6] LEGEA nr. 104 din 15 iunie 2011, www.calitateaer.ro/parametri.php‎

[7] LEGEA nr. 104 din 15 iunie 2011, www.calitateaer.ro/parametri.php‎

[8] www.mdrl.ro/_documente/constructii/punct_info/Ordin_577_2008.pdf

[9] http://www.scribd.com/doc/197210872/NE-012-2

[10] http://www.scribd.com/doc/197210872/NE-012-2

[11] http://www.mmediu.ro/protectia_mediului/legislatie_orizontala

[12] http://lege5.ro/Gratuit/geytgobt/ordinul-nr-184-1997

[13] http://www.mmediu.ro/protectia_mediului/legislatie_orizontala

[14] Transportul și dispersia poluanților – Simion Hâncu, Gabriela Marin.

Editura Carte Universitara 2008

[15] http://www.statii-beton.ro/statie-betoane-mobila-eca-1000-statii-beton- mobile.html

[16] Diagrama evoluției diferitelor substanțe chimice în atmosferă – www.qedoc.org

[17] http://www.mediu.ro/legislatie/calitate_aer.htm

[18] http://apmms.anpm.ro/pages/categorie/53

[19] http://www.cjmures.ro:1880/urbanism/patj/parte_scrisa

. BIBLIOGRAFIE

[1] http://mmut.mec.upt.ro

[2] http://www.calitateaer.ro/parametrii.php

[3] LEGEA nr. 104 din 15 iunie 2011, www.calitateaer.ro/parametri.php‎

[4] LEGEA nr. 104 din 15 iunie 2011, www.calitateaer.ro/parametri.php‎

[5] LEGEA nr. 104 din 15 iunie 2011, www.calitateaer.ro/parametri.php‎

[6] LEGEA nr. 104 din 15 iunie 2011, www.calitateaer.ro/parametri.php‎

[7] LEGEA nr. 104 din 15 iunie 2011, www.calitateaer.ro/parametri.php‎

[8] www.mdrl.ro/_documente/constructii/punct_info/Ordin_577_2008.pdf

[9] http://www.scribd.com/doc/197210872/NE-012-2

[10] http://www.scribd.com/doc/197210872/NE-012-2

[11] http://www.mmediu.ro/protectia_mediului/legislatie_orizontala

[12] http://lege5.ro/Gratuit/geytgobt/ordinul-nr-184-1997

[13] http://www.mmediu.ro/protectia_mediului/legislatie_orizontala

[14] Transportul și dispersia poluanților – Simion Hâncu, Gabriela Marin.

Editura Carte Universitara 2008

[15] http://www.statii-beton.ro/statie-betoane-mobila-eca-1000-statii-beton- mobile.html

[16] Diagrama evoluției diferitelor substanțe chimice în atmosferă – www.qedoc.org

[17] http://www.mediu.ro/legislatie/calitate_aer.htm

[18] http://apmms.anpm.ro/pages/categorie/53

[19] http://www.cjmures.ro:1880/urbanism/patj/parte_scrisa