上海市奉賢區(qū)O3污染特征與影響分析

梁星

(上海市奉賢區(qū)環(huán)境監(jiān)測(cè)站，上海 201400)

臭氧(O3)是地球大氣中一種重要的痕量氣體[1]，主要存在于大氣平流層中，O3可以吸收對(duì)人體有害的短波紫外線從而保護(hù)人類健康。但近地面對(duì)流層中高濃度的O3會(huì)損害人類健康，影響植物生長(zhǎng)，是城市大氣環(huán)境重要的污染物之一。全球O3污染呈現(xiàn)不斷惡化的趨勢(shì)，這與氣候變化、人為污染的加重以及全球范圍的污染物區(qū)域性傳播等有關(guān)。近年來(lái)我國(guó)一些特大城市 (如北京 、上海 、廣州等 )O3超標(biāo)嚴(yán)重，而且趨勢(shì)還在加劇[2]。2019年，全國(guó)337個(gè)城市的O3平均質(zhì)量濃度[ρ(O3)]為148 μg/m3，我國(guó)《環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)》(GB 3095—2012)O3日最大8小時(shí)平均質(zhì)量濃度[ρ(O3-8 h)]二級(jí)標(biāo)準(zhǔn)限值為160 μg/m3。全國(guó)以O(shè)3為首要污染物的超標(biāo)天數(shù)占總超標(biāo)天數(shù)的41.8%，僅次于以細(xì)顆粒物(PM2.5)為首要污染物的超標(biāo)天數(shù)(占比45%)[3]。近年來(lái)，O3污染已引起廣泛關(guān)注，全國(guó)各個(gè)地區(qū)從不同角度對(duì)O3污染特征進(jìn)行了研究[4-7]。已有研究表明，環(huán)境空氣的首要污染物由PM2.5逐漸轉(zhuǎn)變?yōu)镺3，大氣復(fù)合污染特征更趨明顯，以O(shè)3為主的大氣污染問(wèn)題突出。O3污染具有明顯的季節(jié)變化特征，春、夏污染最為嚴(yán)重，且與氣象因素關(guān)系密切，高溫、低濕度、低風(fēng)速的氣象條件更易造成O3污染。氮氧化物(NOX)和揮發(fā)性有機(jī)物(VOCs)等大氣污染物的光化學(xué)反應(yīng)為O3生成提供了重要的前體物。對(duì)流層 O3的來(lái)源主要有2種：平流層的垂直方向傳遞和對(duì)流層光化學(xué)過(guò)程產(chǎn)生[8-9]。對(duì)流層中O3產(chǎn)生的主要途徑為：二氧化氮(NO2)的光解生成O3；VOCs被氧化為活性自由基，尤其是過(guò)氧烷基自由基(RO2·)，可引起一氧化氮( NO)向NO2轉(zhuǎn)化，進(jìn)一步提供了生成O3的 NO2源[10]。

2019年上海市ρ(O3-8 h)為98 μg /m3，對(duì)O3日最大8 h平均濃度第90百分位數(shù)[ρ(O3-8h-90)]進(jìn)行評(píng)價(jià)，結(jié)果為達(dá)標(biāo)。奉賢區(qū)2019年ρ(O3-8 h)為104 μg /m3，對(duì)ρ(O3-8h-90)進(jìn)行評(píng)價(jià)，結(jié)果為未達(dá)標(biāo)。由此可見(jiàn)奉賢區(qū)O3污染程度比上海市平均水平更加嚴(yán)重。

現(xiàn)以上海市奉賢區(qū)近6 a來(lái)逐日O3自動(dòng)監(jiān)測(cè)數(shù)據(jù)為基礎(chǔ)，研究O3質(zhì)量濃度的時(shí)間、空間變化特征及其與氣象影響因素之間的關(guān)系，探索奉賢區(qū)O3的變化規(guī)律及影響因素，為奉賢區(qū)未來(lái)O3污染防治提供依據(jù)。

1 研究方法

1.1 監(jiān)測(cè)時(shí)間

2014年1月1日—2019年12月31日。

1.2 監(jiān)測(cè)點(diǎn)位

上海市奉賢區(qū)位于上海西南部，南臨杭州灣，處于長(zhǎng)江三角洲沖積平原，臨江瀕海，屬于亞熱帶季風(fēng)氣候，常年主導(dǎo)風(fēng)為東南風(fēng)，四季分明，冬、夏長(zhǎng)，春、秋短。奉賢區(qū)共有4 個(gè)環(huán)境空氣自動(dòng)監(jiān)測(cè)站點(diǎn)，分別為南橋站、海灣站、新城站和奉浦站。4個(gè)站點(diǎn)分別分布在奉賢區(qū)的老城區(qū)、背景區(qū)、新城區(qū)和工業(yè)區(qū)，能夠真實(shí)反映奉賢區(qū)的空氣質(zhì)量狀況(圖1)。

1.3 監(jiān)測(cè)儀器

O3數(shù)據(jù)由Thermo 49i型O3分析儀[量程為(0～0.05)×10-6～200 ×10-6，(0～0.1)～400 mg/m3，美國(guó)賽默飛世爾科技公司]采集，NO2數(shù)據(jù)由Thermo 42i型氮氧化物分析儀[量程為(0～0.05)×10-6～20×10-6，(0～0.1)～30 mg/m3，

美國(guó)賽默飛世爾科技公司]采集，CO數(shù)據(jù)由Thermo 48i紅外吸收CO分析儀[量程為(0～1)×10-6～10 000×10-6，美國(guó)賽默飛世爾科技公司]采集。設(shè)備運(yùn)行參照《環(huán)境空氣氣態(tài)污染物(SO2、NO2、O3、CO)連續(xù)自動(dòng)監(jiān)測(cè)系統(tǒng)運(yùn)行和質(zhì)控技術(shù)規(guī)范》(HJ 818—2018)進(jìn)行。

1.4 數(shù)據(jù)來(lái)源

ρ(O3-8 h)數(shù)據(jù)來(lái)自上海市奉賢區(qū)4個(gè)環(huán)境空氣自動(dòng)監(jiān)測(cè)站點(diǎn)，所有數(shù)據(jù)均經(jīng)過(guò)嚴(yán)格質(zhì)量控制。季節(jié)劃分為春季(3—5月)，夏季(6—8月)，秋季(9—11月)，冬季(12—次年2月)。每周對(duì)O3分析儀進(jìn)行1次零點(diǎn)/跨度漂移檢查，每月對(duì)NO2、CO分析儀進(jìn)行1次零點(diǎn)/跨度漂移檢查，每季度進(jìn)行1次儀器多點(diǎn)線性校準(zhǔn)，每年對(duì)各站點(diǎn)O3校準(zhǔn)儀開(kāi)展1次O3標(biāo)準(zhǔn)傳遞。

地面氣象觀測(cè)數(shù)據(jù)包括氣溫、相對(duì)濕度、風(fēng)速和風(fēng)向，數(shù)據(jù)來(lái)自國(guó)家氣象局奉賢區(qū)氣象監(jiān)測(cè)點(diǎn)位(CN1010210000,東經(jīng)121°45′84.7″，北緯30°91′23.5″)。

2 結(jié)果與討論

2.1 O3的時(shí)間變化特征

2.1.1 O3的年際變化特征

2014—2019年奉賢區(qū)4個(gè)站點(diǎn)的O3質(zhì)量濃度統(tǒng)計(jì)見(jiàn)表1。由2014—2019年平均氣溫可知，2017年平均氣溫為6 a內(nèi)最高(18.0 ℃)，光化學(xué)反應(yīng)增強(qiáng)，上海地區(qū)副熱帶高壓環(huán)流指數(shù)異常偏強(qiáng)、面積異常偏大、西伸脊點(diǎn)異常偏西[11]，導(dǎo)致該年出現(xiàn)了O3質(zhì)量濃度最大值(352 μg/m3)，與郭欣曈等[12]研究結(jié)果一致。其他年份O3質(zhì)量濃度年均值趨勢(shì)大體平穩(wěn)，且略微下降，ρ(O3)最大值呈現(xiàn)逐步上升的趨勢(shì)。

2.1.2 O3的月際變化特征

2014—2019年4個(gè)站點(diǎn)O3質(zhì)量濃度月均值變化見(jiàn)圖2(a)—(f)。由圖2可見(jiàn)，1—4月ρ(O3)呈現(xiàn)逐步上升趨勢(shì)；5—8月ρ(O3)達(dá)到峰值，每年6月為上海的梅雨季，降水量增多，期間會(huì)出現(xiàn)ρ(O3)下降的現(xiàn)象；9月ρ(O3)逐步降低。研究結(jié)果與易睿等[13]關(guān)于長(zhǎng)三角地區(qū)城市O3污染特征的研究結(jié)果一致。其中2014，2017和2018年出現(xiàn)明顯的雙峰型變化趨勢(shì)；每年4—10月均會(huì)出現(xiàn)ρ(O3)的最高點(diǎn)和次高點(diǎn)。研究結(jié)果與相鄰的金山區(qū)的ρ(O3)峰值變化類似[14]。4個(gè)站點(diǎn)ρ(O3)變化總趨勢(shì)相似，海灣站點(diǎn)作為背景點(diǎn)，其值最低，奉浦站點(diǎn)由于企業(yè)排放點(diǎn)比較密集，其值最高。

2.1.3 O3的季節(jié)變化特征

2014—2019年奉賢區(qū)ρ(O3)季節(jié)均值和超標(biāo)天數(shù)變化見(jiàn)圖3。由圖3可見(jiàn)，奉賢區(qū)近地面ρ(O3)的季節(jié)變化規(guī)律為：夏季>春季>秋季>冬季，研究結(jié)果與長(zhǎng)三角地區(qū)的淮安市、南京市以及嘉興市(嘉善縣善西)類似[15-17]，夏季和春季的ρ(O3)明顯大于秋季和冬季，這主要與氣象條件有關(guān)。夏季溫度高，太陽(yáng)輻射強(qiáng)，光化學(xué)反應(yīng)強(qiáng)烈，導(dǎo)致ρ(O3)升高。而春季ρ(O3)出現(xiàn)小高峰的原因尚存在爭(zhēng)議，有研究表明，對(duì)流層層頂折疊會(huì)造成平流層O3向?qū)α鲗虞斔?，較強(qiáng)的光化學(xué)反應(yīng)導(dǎo)致ρ(O3)升高[18-19]。也有研究認(rèn)為，可能是由于冬季的NOX和VOCs等前體物積累，隨著溫度上升而導(dǎo)致光化學(xué)反應(yīng)加快造成ρ(O3)升高[18]。秋季和冬季溫差大、日照時(shí)間短，導(dǎo)致ρ(O3)明顯下降[20]。O3季節(jié)超標(biāo)趨勢(shì)與季節(jié)濃度變化趨勢(shì)基本一致，春季和夏季超標(biāo)天數(shù)占比較大，其中春季超標(biāo)天數(shù)占總超標(biāo)天數(shù)的29.6%，夏季超標(biāo)天數(shù)占總超標(biāo)天數(shù)的49.7%，秋季超標(biāo)天數(shù)占總超標(biāo)天數(shù)的20.1%，冬季超標(biāo)天數(shù)僅占總超標(biāo)天數(shù)的0.6%。

圖2 2014—2019年奉賢區(qū)ρ(O3)月均值變化

圖3 2014—2019年奉賢區(qū)ρ(O3)季節(jié)均值和超標(biāo)天數(shù)變化

2.2 O3與氣象影響因素的關(guān)系

2.2.1 O3與日均氣溫的關(guān)系

氣溫(T)是影響O3光化學(xué)反應(yīng)的重要因素，有研究表明，溫度越高，太陽(yáng)輻射越強(qiáng)，光化學(xué)反應(yīng)也越強(qiáng)烈，使得ρ(O3)升高[21]。2014—2019年奉賢區(qū)ρ(O3)日均值和日均氣溫的線性關(guān)系見(jiàn)表2。由表2可見(jiàn)，2014—2019年奉賢區(qū)近地面ρ(O3)日均值(y)與日均氣溫(x)均呈顯著正相關(guān)關(guān)系，相關(guān)系數(shù)(r)為0.664～0.868。

表2 2014—2019年奉賢區(qū)ρ(O3)日均值和日均氣溫的線性關(guān)系

2014—2019年奉賢區(qū)不同日均氣溫下ρ(O3)日均值和超標(biāo)率見(jiàn)圖4。由圖4可見(jiàn)，O3的超標(biāo)率隨著溫度的升高呈現(xiàn)上升趨勢(shì)，日均氣溫>10 ℃時(shí)，O3超標(biāo)現(xiàn)象逐漸顯現(xiàn)，超標(biāo)率和日均氣溫呈顯著正相關(guān)關(guān)系。2014—2019年日均氣溫為30～35 ℃的天數(shù)為130 d，其中O3超標(biāo)天數(shù)總計(jì)為63 d，超標(biāo)率達(dá)48.5%。同時(shí)，隨著日均氣溫的上升，ρ(O3)上升幅度也越來(lái)越大。

圖4 2014—2019年奉賢區(qū)不同日均氣溫下ρ(O3)日均值和超標(biāo)率

2.2.2 O3與相對(duì)濕度的關(guān)系

2014—2019年奉賢區(qū)不同相對(duì)濕度(RH)下ρ(O3)日均值和超標(biāo)率見(jiàn)圖5。由圖5可見(jiàn)，相對(duì)濕度>30%時(shí)，O3開(kāi)始出現(xiàn)超標(biāo)現(xiàn)象；相對(duì)濕度為40%～50%時(shí)，超標(biāo)率出現(xiàn)最大值；之后隨著相對(duì)濕度的增加，超標(biāo)率緩慢下降；相對(duì)濕度≥80%時(shí)，超標(biāo)率急劇下降。造成相對(duì)濕度越高近地面ρ(O3)越低的原因是：一方面，相對(duì)濕度較高時(shí)，大氣中的水汽影響太陽(yáng)紫外線輻射強(qiáng)度，導(dǎo)致大氣光化學(xué)反應(yīng)減弱[13]；另一方面，相對(duì)濕度較高時(shí)，空氣中的水汽所含的HO2·、OH·等自由基可快速消耗O3，從而使ρ(O3)降低[22]。

圖5 2014—2019年奉賢區(qū)不同相對(duì)濕度下ρ(O3)日均值和超標(biāo)率

2.2.3 O3與日均風(fēng)速的關(guān)系

奉賢區(qū)日均風(fēng)速(V)為 0.5～4 m/s，其季節(jié)規(guī)律為：春季>夏季>冬季>秋季。2014—2019年奉賢區(qū)不同日均風(fēng)速下ρ(O3)日均值和超標(biāo)率見(jiàn)圖6。由圖6可見(jiàn)，隨著日均風(fēng)速的增加，超標(biāo)率表現(xiàn)為先升后降的規(guī)律。當(dāng)風(fēng)速>0.5 m/s時(shí)，O3開(kāi)始出現(xiàn)超標(biāo)現(xiàn)象；當(dāng)風(fēng)速達(dá)到1.5～2 m/s時(shí)，超標(biāo)率達(dá)到峰值；之后隨著風(fēng)速的增大，超標(biāo)率呈現(xiàn)下降的趨勢(shì)。風(fēng)速增大可使大氣垂直擴(kuò)散和水平擴(kuò)散作用同步增強(qiáng)。垂直擴(kuò)散會(huì)使O3向地面混合，水平擴(kuò)散會(huì)稀釋O3，這2種作用同時(shí)發(fā)生。當(dāng)風(fēng)速較低時(shí)，O3的混合作用強(qiáng)于擴(kuò)散作用，風(fēng)速增加會(huì)使O3超標(biāo)率增加；當(dāng)風(fēng)速超過(guò)特定值時(shí)，擴(kuò)散作用又占主導(dǎo)地位，風(fēng)速增加反而會(huì)使O3超標(biāo)率降低[13]。

圖6 2014—2019年奉賢區(qū)不同日均風(fēng)速下ρ(O3)日均值和超標(biāo)率

2.2.4 O3與風(fēng)向的關(guān)系

2014—2019年奉賢區(qū)不同風(fēng)向下O3超標(biāo)率風(fēng)玫瑰圖見(jiàn)圖7。由圖7可見(jiàn)，風(fēng)向?yàn)檎巷L(fēng)時(shí)，O3超標(biāo)率最高，為41%；風(fēng)向?yàn)槲髂巷L(fēng)時(shí)，超標(biāo)率為30%；風(fēng)向?yàn)槲髂掀巷L(fēng)(南西南)時(shí)，超標(biāo)率為29%。這3個(gè)風(fēng)向時(shí)ρ(O3)日均值分別為129～150 μg/m3，高于其他風(fēng)向的ρ(O3)日均值為87～126 μg/m3。奉賢區(qū)地處上海市的西南部，南鄰杭州灣，杭州灣沿岸分布著涂裝、化工、紡織印染等多種類型的工業(yè)企業(yè)，同時(shí)南部分布著上海市星火開(kāi)發(fā)區(qū)和上海市化學(xué)工業(yè)區(qū)，西南方向分布有金山石化和金山二工工業(yè)區(qū)，因此奉賢區(qū)的ρ(O3)變化對(duì)風(fēng)向的變化較敏感， 工業(yè)區(qū)排放的VOCs破壞了 NO2-NO -O3的光解循環(huán)，使O3累積，導(dǎo)致輸送效應(yīng)較為明顯[23]，對(duì)下風(fēng)向的環(huán)境空氣自動(dòng)監(jiān)測(cè)站點(diǎn)的O3污染有著很大的影響(測(cè)定結(jié)果為129～150 μg/m3)，奉賢區(qū)各環(huán)境空氣自動(dòng)監(jiān)測(cè)站點(diǎn)在風(fēng)向?yàn)槟巷L(fēng)、西南風(fēng)和西南偏南風(fēng)時(shí)，O3超標(biāo)率和濃度均較高。

圖7 2014—2019年奉賢區(qū)不同風(fēng)向下O3超標(biāo)率風(fēng)玫瑰圖

2.3 前體污染物對(duì)O3的影響

2.3.1 NOx對(duì)O3的影響

NO2的光解是光化學(xué)反應(yīng)產(chǎn)生O3最根本的反應(yīng)[24]。NO容易和氧化劑發(fā)生反應(yīng)轉(zhuǎn)化成NO2，NO2發(fā)生光解形成NO-NO2的動(dòng)態(tài)平衡。奉賢區(qū)的NOX主要來(lái)自工業(yè)源和交通源，呈現(xiàn)夏季濃度低，冬季濃度高的趨勢(shì)。2014—2019年奉賢區(qū)O3、NO、NO2與NOX質(zhì)量濃度月度變化見(jiàn)圖8，相關(guān)系數(shù)(r)及P值見(jiàn)表3。

圖8 2014—2019年奉賢區(qū)ρ(O3)、ρ(NO)、ρ(NO2)和ρ(NOX)濃度月度變化

表3 2014—2019年O3與NO、NO2、NOX的相關(guān)系數(shù)(r)及P值

由圖8和表3可見(jiàn)，ρ(O3)月均值與ρ(NO)、ρ(NO2)和ρ(NOX)月均值均呈負(fù)相關(guān)關(guān)系，ρ(O3)月均值與ρ(NO)月均值相關(guān)性最好，這與許波和胡正華等[15-16]的研究結(jié)果一致。其原因是隨著太陽(yáng)輻射的增強(qiáng)，光化學(xué)反應(yīng)加劇，O3濃度升高，NO更易轉(zhuǎn)化為NO2。

2.3.2 CO對(duì)O3的影響

CO是O3的前體物，CO參與光化學(xué)反應(yīng)，促進(jìn)O3生成的機(jī)理可以解釋為：CO、O2和 NO反應(yīng)生成 NO2，NO2、O2和第3體M 生成 O3。因此隨著 O3濃度的增加，CO 濃度下降[25]。2014—2019年奉賢區(qū)ρ(O3)和ρ(CO)月均值變化見(jiàn)圖9，計(jì)算spearman相關(guān)系數(shù)(r)為-0.888，P值<0.001，由此可見(jiàn)O3和CO的濃度變化呈顯著負(fù)相關(guān)，與趙長(zhǎng)民和應(yīng)誠(chéng)威等[26-27]研究結(jié)果一致。

圖9 2014—2019年奉賢區(qū)ρ(O3)和ρ(CO)月均值變化

2.4 O3的空間分布

2014—2019年奉賢區(qū)4個(gè)監(jiān)測(cè)站點(diǎn)的ρ(O3)日均值及超標(biāo)天數(shù)見(jiàn)圖10。

由圖10可見(jiàn)，奉浦站點(diǎn)ρ(O3)日均值最高，且超標(biāo)天數(shù)最多。該站點(diǎn)臨近工業(yè)區(qū)，企業(yè)分布密集，排放的O3前體物NOX和VOCs較多，經(jīng)過(guò)光化學(xué)反應(yīng)后對(duì)近地面ρ(O3)有較大貢獻(xiàn)。

3 結(jié)論

(1)2014—2019年奉賢區(qū)(除2017年外)ρ(O3)年均值趨勢(shì)大體平穩(wěn)且略微下降,ρ(O3-8 h)中位值和ρ(O3)最大值呈現(xiàn)逐步上升的趨勢(shì)。ρ(O3)的最高點(diǎn)和次高點(diǎn)均出現(xiàn)在4—10月。ρ(O3)的季節(jié)變化規(guī)律為：夏季>春季>秋季>冬季，春季和夏季O3超標(biāo)天數(shù)占比較大。

(2)2014—2019年奉賢區(qū)近地面ρ(O3)日均值與日平均氣溫均呈顯著正相關(guān)關(guān)系；隨著風(fēng)速和濕度的上升，ρ(O3)日均值和超標(biāo)率呈現(xiàn)先升后降的趨勢(shì)。風(fēng)向?yàn)槟巷L(fēng)、西南風(fēng)和西南偏南風(fēng)時(shí)，O3超標(biāo)率和濃度均較高。

(3)O3與CO、NO、NO2、NOX的相關(guān)系數(shù)(r)中，O3與NO的相關(guān)性最好，r=-0.902。奉賢區(qū)4個(gè)空氣自動(dòng)站中，奉浦站點(diǎn)ρ(O3)日均值最高，且超標(biāo)天數(shù)最多，主要由于該站點(diǎn)臨近工業(yè)區(qū)，企業(yè)分布密集，排放的O3前體物NOX和VOCs較多所致。奉賢區(qū)VOCs主要來(lái)源于本地及周邊的工業(yè)企業(yè)排放，其對(duì)O3的影響程度和規(guī)律還須進(jìn)一步研究。