TÜRKiYE’DE BİR ÇİMENTO FABRİKASI İÇİN HAVA KALİTESİ MODELLEME ÇALIŞMASI

Barış R.CANPOLAT*, Aysel T.ATIMTAY*, Ismet MUNLAFALIOĞLU**, Ersan KALAFATOĞLU***

*ODTÜ, Çevre Mühendisliği Bölümü, 06531 Ankara

**T.Ç.M.B. Eskişehir Yolu, 9.Km, 06530 Ankara

***TÜBİTAK-Marmara Araştırma Merkezi, P.K.21, 41470 Gebze, Kocaeli

Özet: Bu çalışmada, Türkiye Çimento Müstahsilleri Birliği (TÇMB) tarafından bir çimento fabrikası için yapılan emisyon ölçüm sonuçları kullanılarak ortam havası kalitesi için bir modelleme  çalışması  yapılmıştır.  Modelleme  çalışması  sonunda  bulunan  yer  seviyesi konsantrasyonları ile imisyon ölçüm sonuçları karşılaştırılmıştır. Matematiksel model olarak Industrial Source Complex-ISC Version 3.0 kullanılmıştır. Çimento fabrikası için bir yıllık hava kalitesi  ölçüm sonuçları TÇMB tarafından sağlanmıştır. Model  sonuçları ile ölçüm sonuçları arasındaki farkın nedenleri açıklanmaya çalışılmıştır. Bundan sonra yapılacak hava kalitesi ölçümleri yapılmadan önce nelere dikkat edilmesi konusunda öneriler geliştirilmiştir. Anahtar kelimeler: Emisyonlar, modelleme, çimento fabrikası toz dağılımı, ISC modeli

1.   GİRİŞ

Çimento endüstrisi, tüm dünyada hava kirliliğine sebep olan en önemli endüstrilerden biri olarak  düşünülmüştür.  Yüksek  enerji  tüketimi  ve  partikül  madde  emisyon  potansiyeli açısından hava kirliliği kaynakları arasında öncelikli listeye girmesine neden olmuştur. Bugün Türkiye çapında 40 çimento fabrikası ile 16 öğütme-paketleme fabrikası mevcuttur. Çimento fabrikaları 1998 yılı sonu itibarı ile 29 148 000 ton/yıl klinker ve 37 488 000 ton/yıl çimento üretme kapasitesine sahiptir . Fabrikaların üretim kapasiteleri ortalama 600 000-650 000 ton/yıl’dır.  Çimento  sektöründe  sadece  üretim  açısından  değil,  çevresel  kirliliği  önleme

açısından da en son teknoloji takip edilmektedir. Çimento fabrikaları HKKY’de verilen toz emisyon değerlerinden çok daha düşük emisyon değerlerini sağlamaktadır. Türk çimento

sektörü çimento üretiminde son yıllarda Avrupa’da birinci ve dünyada 8. sırada gelmektedir. Kapasite açısından da İtalya ve Almanyadan sonra Avrupada 3. sıradadır.

Çalışmanın amacı, seçilen bir çimento fabrikası için yapılmış emisyon ve imisyon ölçümlerinden yararlanarak model çalışması yoluyla yer konsantrasyonlarını hesaplamak ve bulunan değerleri ölçüm sonuçları ile karşılaştırmaktır. Hava kirlenmesine katkı değerlerinin hesaplanmasında pahalı ve zaman alıcı ölçüm sürecinin model sonuçları ile bulunmasının ne kadar olası olduğu da bu çalışmanın amaçları arasındadır.

2.   KULLANILAN MATEMATİKSEL MODEL

Kısa dönemli Industrial Source Complex (ISCST) modeli, tipik bir endüstriyel kaynaktan yayılan emisyonların geniş bir alanda dağılımını hesaplayabilmektedir. Modelin temeli Gauss Dağılımına dayanmaktadır (1). Emisyon kaynakları 4 ayrı grupta toplanmaktadır. Bunlar; nokta kaynaklar, hacimsel kaynaklar, alansal kaynaklar ve açık alan kaynaklarıdır. Hacimsel ve alansal kaynak opsiyonları çizgisel kaynakların simulasyonu için kullanılabilmektedir.

Kısa dönemli ISC modeli, dumanın yükselmesi, taşınması, difüzyonu ve çökelmesi ile ilgili koşulların tanımı için saatlik meterolojik verilere ihtiyaç duymaktadır. Model her kaynak (source)  ve  alıcı  (receptor)  kombinasyonu  için  çökelme  veya  konsantrasyon  değerlerini, meterolojik girdilerin verildiği her saat için hesaplamaktadır. Ayrıca, kullanıcı tarafından belirlenen zaman aralıklarındaki ortalama konsantrasyonu da hesaplayabilmektedir. Model, nokta  kaynakların  yakın  civarındaki  binaların  partikül  maddeler  üzerindeki  aerodinamik etkilerini belirleyebilmek, partikül maddelerin birim alandaki kuru ve yaş çökelme hızlarını ve ayrıca toplam çökelme hızını hesaplamak için değişik algoritmalara sahiptir. Model, birden fazla emisyon kaynağını da (nokta, alan ve yer seviyesindeki açık kaynak olarak) kullanabilir. Modelde ele alınan kaynaklardan çıkan emisyonların hızı sabit kabul edilebileceği gibi ay, mevsim, ya da başka bir zaman periyodu için değişken olarak alınabilir. Emisyon kaynakları, tek olarak, ya da grup olarak belirlenebilir. Kısa dönemli model (ISCST) için kullanıcı, bazı emisyon kaynakları veya kaynakların tamamı için saatlik emisyon hızı da tanımlayabilir (1).

Uzun dönemli ISC Modeli (ISCLT) için uzun süre içinde toplanmış meteorolojik veriler kullanılmaktadır. Uzun dönem için kullanılan ISC modelinin kısa dönem için kullanılan

modelden tek farkı, kullanılan meteorolojik verilerdedir. ISCST saatlik meteorolojik veriler, ISCLT yıllık meteorolojik veriler kullanmaktadır.

Alıcıların Belirlenmesindeki Opsiyonlar

ISC Modelinin, alıcıların yerlerinin belirlenmesinde önemli bir esnekliği vardır. Kullanıcı, model çalışması için çok sayıda alıcı noktasından oluşan bir grid sistemi belirleyebildiği gibi, grid sisteminde alıcıların yerlerinin belirlenmesinde kartezyen ve polar koordinatları aynı sistem içerisinde kullanabilir (1).

Meteorolojik Opsiyonlar

Kısa dönemli model (ISCST3), ön işlemcisi tarafından oluşturulan formatlanmamış sıralı dosyaları kullanabilmektedir. Bu durum, çökelme algoritmalarının kullanılmadığı durumlar için geçerli değildir. Kısa dönemli model, kuru ve yaş çökelme algoritmalarını içermektedir. Kuru çökelme algoritması, ISCST3 ön işlemcisi tarafından sağlanabilecek Monin-Obukhov uzunluğu ve yüzey sürtünme hızı gibi ilave meteorolojik girdilere ihtiyaç duymaktadır. Yaş çökelme algoritması ise ayrıca yağış verilerine gereksinim duymaktadır (1).

3.   MODEL ÇALIŞMASI

Örnek olarak seçilen bir çimento fabrikasının yakın civarındaki havada asılı partikül madde (PM10) ve çöken tozların yer seviyesi konsantrasyonları, kısa dönemli ISCST3 modeli kullanılarak hesaplanmıştır. Bu çalışmadaki çimento fabrikasının seçilme nedenleri:

a)  Fabrikanın yakın civarında 1 yıl süre ile TÇMB tarafından hava kalitesi ve çöken toz ölçümlerinin yapılması

b)  Fabrikanın konumu itibariyle fabrika civarındaki diğer toz emisyon kaynaklarının azlığı (fabrikanın çevresinde yoğun bir trafiğin ve yerleşim yerinin olmaması, yine yakın çevrede başka bir sanayi tesisinin bulunmaması gibi)

c)  Fabrikanın imisyon ölçümlerinin yapıldığı zaman aralığında emisyon verilerinin mevcut olmasıdır.

Seçilen çimento fabrikasında emisyon ve hava kalitesi (imisyon) ölçümleri, TÇMB Kalite ve Çevre Kontrol Müdürlüğü tarafından 1 yıl süre ile gerçekleştirilmiştir. Bu çalışma sırasında, çimento fabrikasından kaynaklanan emisyonun sabit olduğu varsayılmıştır.

Model çalışmaları için, Hava Kalitesinin Korunması Yönetmeliğine (HKKY) göre (2) çimento fabrikasının merkezinde yer aldığı 4 km x 4 km’lik bir alan seçilmiştir. TÇMB’nin gerçekleştirdiği hava kalitesi ölçümleri çimento fabrikasına eşit uzaklıkta 4 ayrı noktada yapılmıştır. Ölçüm noktaları Şekil 1’de gösterilmiştir.

2 km

Şekil 1. Ölçüm Noktalari

Modelde, 4 km x 4 km’lik alan, 200 m x 200 m’lik kare alanlara bölünmüştür. Bu alanların kesiştiği noktaların yükseltileri, 1/25000’lik topografik haritadan okunmuştur. Modelleme çalışmasında kullanılan alanın 3 boyutlu görüntüsü Şekil 2’de verilmiştir. Bu alan üzerinde

21 x 21 olmak üzere toplam 441 alıcı nokta belirlenmiştir.

Kısa dönemli model; rüzgar hızı, rüzgar yönü, karışım yüksekliği (sabah ve öğleden sonra), sıcaklık ve basınç gibi parametreleri içeren saatlik meteorolojik verilere ihtiyaç duymaktadır. Bu veriler Devlet Meteoroloji İşleri Genel Müdürlüğü’nden temin edilmiştir. Çalışma yapılan alandaki  meteorolojik  veriler  Tablo  1’de  verilmiştir.  Daha  önce  42  yıllık  meteorolojik verilerle hazırlanan rüzgar gülü, baskın rüzgar yönünün NE ve NNE olduğunu göstermektedir.

ISCST3 programı için yeterli veri girildikten sonra, yer seviyesinde PM10 ve çöken toz konsantrasyonları bulunmuştur. TÇMB, havada asılı partikül madde ölçümünü (PM10), HKKY’ye göre her ay ayda 5 gün olmak üzere yılda 60 gün gerçekleştirmiş; çöken toz ölçümünü ise, sürekli olarak her gün yapmıştır.

4.    BULGULAR VE YORUM

PM10 Sonuçları

1 no’lu ölçüm yerinde model ile elde edilen sonuçlar, Tablo 1’den görüleceği üzere, son üç sonuç hariç, ölçüm sonuçlarından yüksek çıkmıştır. Nisan 1999’da tamamlanan by-pass bacası revizyonu sonrasında, modelde girdi olarak verilen emisyon ölçümlerinin düşük olması program sonuçlarının da düşük çıkmasına sebep olmuştur. Özellikle, baskın rüzgar yönünün N ve NNE olması, çimento fabrikasının güney-batı tarafında kalan 1 no’lu ölçüm noktasında elde edilen model sonuçlarının yüksek çıkmasına neden olmuştur. Ağustos, Kasım, Aralık, Ocak  ve  Şubat  aylarında  rüzgar  hızı  oldukça  yüksektir.  Böylece,  havada  asılı  partikül maddeler, çimento fabrikasından uzak noktalara taşınmıştır. HKKY’de havada asılı partikül maddeler için 150     g/m3 sınır değer belirlenmiştir. 1 no’lu noktada yapılan hava kalitesi ölçümleri,  sınır  değerinin  %   8-43’ne  karşılık   gelen  bir   kirlilik  göstermiştir.  Model sonuçlarında ise, çimento fabrikasının hava kirliliğine % 0-84 arasında bir katkı getirdiği görülmektedir. Kalafatoğlu ve diğerleri (3), 7 ayrı çimento fabrikasında yaptıkları çalışmada, hava kirlenmesine katkı değerini % 4-43 arasında bulmuşlardır. 1 yıllık model sonuçları irdelendiğinde maksimum hava kirlenmesine katkı değerinin % 475 olduğu tespit edilmiştir.

Tablo 1. Ölçüm Noktalarında Hava Kalitesi Ölçüm ve Model Sonuçları (PM10 için) (4)

Ay	1 no’lu ölçüm   noktası	2 no’lu ölçüm   noktası	3 no’lu ölçüm   noktası	4 no’lu ölçüm   noktası
Ö (  g/m3)	M (  g/m3)	Ö (  g/m3)	M (  g/m3)	Ö (  g/m3)	M (  g/m3)	Ö (  g/m3)	M (  g/m3)
7/98	31.4	126.1	28.6	72.2	35.0	195.3	37.4	0.1
8/98	62.4	115.2	62.6	208.9	54.0	7.9	55.6	0.1
9/98	39.2	87.4	34.0	27.1	41.4	16.9	36.0	0.7
10/98	30.6	59.3	31.6	17.0	31.0	39.0	37.6	27.1
11/98	47.6	123.0	49.4	50.4	64.8	0.6	123.4	0.2
12/98	55.2	120.7	52.0	1.3	63.6	0.0	98.8	0.0
1/99	64.8	19.7	69.2	7.6	72.0	14.8	131.2	50.8
2/99	35.4	75.4	34.6	3.1	39.4	5.0	84.0	13.1
3/99	37.4	10.7	38.2	4.8	37.0	82.6	65.2	28.7
4/99	12.6	0.0	10.2	0.4	8.0	0.2	16.4	0.2

5/99	23.8	0.3	22.2	0.3	23.8	0.0	23.2	0.0
6/99	48.2	0.3	43.0	0.2	48.2	0.3	43.4	0.3

Ö:Ölçüm         M: Model

2 no’lu ölçüm noktasında elde edilen model sonuçları, 1 no’lu ölçüm noktasında elde edilen sonuçlara göre daha düşük seviyededir. Bu noktada yapılan hava kalitesi ölçümleri ise bütün yıl boyunca pek değişiklik göstermemiştir. En yüksek model sonucunun alındığı Ağustos

1998 ayında, ikinci en yüksek ölçüm sonucu elde edilmiştir. Ölçüm sonuçlarına göre hava kirlenmesine katkı değeri, % 7-46 arasında, model sonuçlarına göre ise % 0-138 arasında bulunmuştur.

3 no’lu ölçüm noktasında elde edilen model sonuçları, diğer ölçüm noktalarında elde edilen sonuçlara göre daha düşüktür. Ölçüm sonuçlarına göre hava kirlenmesi katkı değeri % 5-48 arasında, model sonuçlarına göre ise % 0-130 arasında bulunmuştur. Bu noktadaki model sonuçlarının düşük çıkmasının sebebi baskın rüzgar yönünün N ve NNE olmasıdır.

En yüksek ölçüm değerlerinin bulunduğu yer olan 4 no’lu ölçüm noktası, model sonuçlarına göre en düşük PM10 konsantrasyonunun görüldüğü noktadır. Bu durum, ölçüm noktasının

çok yakınında mevcut olan yerleşim yerinin etkisiyle açıklanabilir. Ölçüm sonuçlarına göre %

11-87 arasında bulunan hava kirlenmesi katkı değeri, model sonuçlarına göre % 0-34 arasında bulunmuştur. Şekil 4’te model çıktıları 3 boyutlu olarak gösterilmiştir.

Çöken Tozlar

Çöken tozlar için yapılan ölçüm sonuçları ile model çıktıları Tablo 2 de verilmiştir. Elde edilen sonuçlar arasında büyük farklar mevcuttur. Bu durumun 3 sebebi olabilir:

Modelde  yağmur  verileri  kullanılmamıştır.  Özellikle,  yağmurun  tozun  çökelmesinde büyük bir etkisi mevcuttur.

4 no’lu ölçüm noktasında bulunan değerlerde, yerleşim yerinin etkisi bulunmaktadır. Özellikle çöken toz için, çimento fabrikası dışındaki etkenlerin etkisi çok olmaktadır.

Modelde yağmur verilerinin kullanılamamasının nedeni, bu bölge için saatlik yağmur ölçümü sonuçlarının mevcut olmamasıdır. Kalafatoğlu ve diğerleri (3) çimento fabrikalarının, çöken tozlar için hava kirlenmesine katkı değerini  % 1-33 olarak bulmuşlardır. HKKY’de çöken tozlar için belirlenen sınır değeri 350 mg/m2/gün’dür.  Bu çalışmada, hava kirlenmesine katkı

değeri ölçüm sonuçlarına göre % 55-1750, model sonuçlarına göre ise %0-43 arasında bulunmuştur. Fabrikadaki revizyondan sonra, çöken toz miktarında önemli bir azalma tesbit edilmiştir.

Tablo 2. Ölçüm ve Model Sonuçlarına göre Çöken Tozlar (mg/m2/gün) (4)

Ay	1 no’lu ölçüm   noktası	2 no’lu ölçüm   noktası	3 no’lu ölçüm   noktası	4 no’lu ölçüm   noktası
Ö	M	Ö	M	Ö	M	Ö	M
9/98	3559	54	1694	49	3053	32	4469	6
10/98	3962	34	840	51	1700	38	3174	14
11/98	2450	44	666	24	494	11	2395	11
12/98	1628	50	483	5	1071	4	6124	59
1/99	3571	28	1509	16	2018	7	5038	22
2/99	4851	35	2461	11	912	22	3980	150
3/99	1500	52	637	43	1868	26	2032	22
4/99	209	0	522	0	381	0	191	0
5/99	592	0	362	0	358	0	458	0
6/99	213	0	498	0	845	0	141	0

Ö:Ölçüm       M: Model

Sonuç

ISCST3  ile  yapılan  model  çalışmaları,  çimento  fabrikasının  konumunun  hava  kirliliği açısından önemini ortaya koymuştur. Türkiye’de bugüne kadar yapılan çalışmalar, hava kirlenmesine katkının özellikle fabrika merkezli 1000 m yarıçaplı bir alan üzerinde olduğunu göstermiştir. Yapılan model çalışmalarında elde edilen sonuçlar ile ölçüm sonuçları arasında önemli bir ilişki tespit edilememiştir. Bu duruma, özellikle ölçüm noktalarının fabrikaya yakınlığı ve yakın civardaki kirlilik kaynaklarının faaliyetleri sebep olmuştur.

5.   KAYNAKLAR

1.   Users’ Guide for the Industrial Source Complex (ISC3) Dispersion Models, Vol.1-2, EPA, 1995.

2.   Hava Kalitesinin Korunması Yönetmeliği, 2.11.1986, Resmi Gazete No. 19269, Ankara.

3.   Kalafatoğlu, E., Örs, N., Gözmen, T., Sain, S., Munlafalıoğlu, İ., 1997, “Air Pollution

Contribution of Some Cement Plants in Turkey”, Environmental Research Forum, Cilt.

7-8, s.191-196.

4.   Canpolat, B.R., 1999, “Calculation of Emission Factors in Cement Plants in Turkey and

A Modelling Study for a Cement Plant”, Y. Lisans Tezi, ODTÜ, Ankara.