陶 燕,羊德容 ,蘭 嵐 ,王洪新 ,王式功 (.蘭州大學資源環(huán)境學院,甘肅 蘭州 730000；.蘭州大學大氣科學學院,甘肅 蘭州 730000)

近年來,空氣污染對人類健康的影響逐漸成為國內(nèi)外研究的熱點之一.WHO對全球3211個城市的研究表明,2000年全球室外空氣可吸入顆粒物污染引起的早死人數(shù)約為79.9萬人,其中亞太地區(qū)為48.7萬人,空氣污染的疾病負擔主要發(fā)生在發(fā)展中國家,亞洲的空氣污染疾病負擔占全球的 2/3[1].眾多研究均表明[2-5]低劣的空氣會刺激呼吸系統(tǒng)使其出現(xiàn)病變.Pikhart等[6]調(diào)查了3680名兒童,多水平分析結(jié)果顯示,NO2每升高10μg/m3,兒童喘息發(fā)作的危險性增加 1.16倍;SO2每升高 10μg/m3,兒童喘息發(fā)作的危險性增加 1.08倍.周燕榮等[7]分析了重慶大氣污染與上呼吸道疾病與病癥的關(guān)系,結(jié)果發(fā)現(xiàn)在高、中、低污染區(qū),上呼吸道疾病與病癥的發(fā)病率分別為22.0%、12.9%和9.1%.

蘭州市是一個以石油化工、冶金和機械工業(yè)為主的工業(yè)城市,空氣污染較為嚴重,再加上地理條件的制約,特有的盆地地形、特殊的氣象條件等因素使蘭州空氣污染物不易擴散.因此本研究利用蘭州市2007年1月1日至2009年12月31日空氣污染物和呼吸系統(tǒng)疾病日入院人數(shù)數(shù)據(jù),運用時間序列方法,通過廣義相加模型分析蘭州市空氣污染對居民呼吸系統(tǒng)疾病的影響,具有明確的地域代表性,同時也為疾病的預防控制和有關(guān)環(huán)境的決策提供科學依據(jù).

1 資料與方法

1.1 疾病資料

呼吸系統(tǒng)疾病日入院資料取自蘭州市 3所大型綜合醫(yī)院2007年1月1日至2009年12月31日的住院登記表.呼吸系統(tǒng)疾病由醫(yī)院專業(yè)人員根據(jù)國際疾病分類標準第 10版(ICD-10:J00-J99)進行分類[8],并剔除了因人為原因造成的疾病.3所醫(yī)院呼吸系統(tǒng)疾病入院人數(shù)占蘭州市區(qū)呼吸系統(tǒng)疾病總?cè)朐喝藬?shù)的30%～40%.

1.2 空氣污染物資料

收集了2007年1月1日至2009年12月31日蘭州市3種主要污染物PM10、SO2、NO2日均濃度,資料來自于蘭州市環(huán)境監(jiān)測站.采用線性內(nèi)插法對缺失數(shù)據(jù)進行填補得到的完整的日均濃度序列.

1.3 氣象資料

收集了2007年1月1日至2009年12月31日蘭州市逐日地面氣象資料,包括平均氣溫、最高氣溫、最低氣溫、氣壓、相對濕度等氣象因子,資料來自于甘肅省氣象局.

1.4 方法

采用半?yún)?shù)廣義相加模型(GAM)[9]分析空氣污染物與呼吸系統(tǒng)疾病日入院人數(shù)的影響,并分別按性別和年齡的不同分層建立模型.相對于總?cè)丝趤碚f,呼吸系統(tǒng)疾病患者入院屬于小概率事件,其實際分布近似 Poisson分布[10],故本研究擬合 Poisson回歸的 GAM 模型.在排除星期幾(day of the week)效應,用樣條平滑函數(shù)排除長期趨勢、氣象等混雜因素影響的基礎(chǔ)上,將當天(lag0)、1d前至 7d前污染物 PM10、SO2、NO2濃度逐一引入模型.具體模型如下:

Log[E(Yk)]=α+ DOW+βXk+s(time,df)+ s(Zk,df)式中:Yk為第k日呼吸系統(tǒng)疾病入院人數(shù);E(Yk)為第k日呼吸系統(tǒng)疾病入院人數(shù)的期望值;α為殘差;β為回歸系數(shù);Xk為第k日污染物濃度;s為樣條平滑函數(shù);df為自由度;DOW為星期啞元虛擬變量;time為日歷時間;Zk為第k日氣象要素.

確定單污染物的最優(yōu)模型后,引入 PM10、SO2和 NO2,分別建立單污染物模型和多污染物模型,分析多污染物協(xié)同作用下空氣污染物對人群健康造成的損害如何變化,由此確定研究期間主要影響呼吸系統(tǒng)入院人數(shù)的危險因子.并且采用 Akaike信息標準(AIC)作為評判標準[11],對回歸模型進行因子選擇與優(yōu)度檢驗.

根據(jù) GAM 模型估算出空氣污染物的回歸系數(shù)β,計算相對危險度(RR)RR=EXP[IQR(inter-quartile range)·β],在此基礎(chǔ)上計算出 RR 的95%的可信區(qū)間(95%CI) EXP[(β±1.96SE)·IQR],污染物濃度每增加 10μg/m3疾病日入院人數(shù)變化的百分比為 10/IQR×(RR-1)×100%.統(tǒng)計分析由R軟件和SPSS軟件實現(xiàn).

2 結(jié)果

2.1 描述性分析

2007年1月1日至2009年12月31日蘭州市大氣中主要污染物PM10、NO2、SO2的日均濃度年均值、主要氣象因素和呼吸系統(tǒng)疾病日入院人數(shù)的描述性統(tǒng)計結(jié)果見表 1.可以看出,蘭州市大氣中主要污染物PM10、NO2、SO2的日均濃度年均值均高于國家空氣質(zhì)量二級標準[12].3年期間呼吸系統(tǒng)疾病入院人數(shù)合計為 19009人,男女性別比將近 2:1;老年(年齡≥65)人群占總呼吸系統(tǒng)疾病入院人數(shù) 25.7%,年齡≤15歲的人群占38.3%;年齡介于15～65歲之間的人群占36.0%.

2.2 空氣污染物及氣象因素間的相關(guān)性

Pearson 相關(guān)性分析結(jié)果見表2,3種污染物之間均存在明顯的正相關(guān),且有統(tǒng)計學意義,說明污染物之間具有一定的協(xié)同作用,其中 NO2與SO2的相關(guān)性較強,相關(guān)系數(shù)為0.430.氣溫、能見度與 3種污染物之間均呈現(xiàn)負相關(guān),且具有統(tǒng)計學意義.

2.3 空氣污染物與日入院人數(shù)的關(guān)系

2.3.1 單污染物模型 圖1顯示了在單污染模型下空氣污染物PM10、SO2和 NO2日平均濃度對居民呼吸系統(tǒng)疾病日入院人數(shù)的 RR.可以看出,PM10對不同人群的影響具有一定的滯后效應,且滯后時間略有差異,全人群和女性人群在滯后2d有統(tǒng)計學意義,其他人群均在滯后 5d具有統(tǒng)計學意義;PM10每升高1個IQR,全人群、男性、女性、年齡≥65歲的老年人群、15～65歲以及≤15歲人群的RR依次為1.015、1.022、1.026、1.039、1.033和1.048.SO2對不同人群的均在滯后3d RR達到最大值,且除女性人群和15～65歲人群外,均具有統(tǒng)計學意義;SO2每升高1個IQR,全人群、男性、女性、年齡≥65歲的老年人群、15～65歲人群以及≤15歲人群的RR依次為1.049、1.061、1.052、1.085、1.072和1.101.NO2對全人群、女性及≤15歲人群的影響在當天具有統(tǒng)計學意義,且每升高 1個 IQR,其 RR分別為 1.040、1.050和1.078;對男性、年齡≥65歲的老年人群及15～65歲人群的影響在滯后3d具有統(tǒng)計學意義,且每升高1個IQR,其RR分別為1.065、1.074和1.066.

表1 蘭州市2007～2009年空氣污染物濃度、氣象因素和呼吸系統(tǒng)疾病入院人數(shù)統(tǒng)計Table 1 Statistics of air pollutants, meteorological factors and hospital admissions for respiratory diseases in Lanzhou,from 2007 to 2009

表2 蘭州市2007～2009年空氣污染物濃度與氣象因素的Pearson相關(guān)系數(shù)Table 2 Pearson correlation coefficients between atmospheric pollution concentrations and meteorological factors in Lanzhou from 2007 to 2009

2.3.2 多污染物模型 表3顯示了在雙污染模型中,分別引入SO2或NO2后,PM10對全人群、男性、老年人群和年齡≤15歲人群呼吸系統(tǒng)疾病入院人數(shù)的影響不大,且仍具有統(tǒng)計學意義,對女性和 15～65歲人群的影響的 RR明顯下降,但不具統(tǒng)計學意義.引入PM10后,SO2對全人群、老年人群和年齡≤15歲人群呼吸系統(tǒng)疾病日入院人數(shù)的影響不大,對男性和15～65歲人群的影響的RR明顯增加,且仍均具有統(tǒng)計學意義,對女性人群影響的RR明顯下降,但不具統(tǒng)計學意義;引入NO2后,SO2對全人群、女性、15～65歲和和年齡≤15歲人群影響的 RR明顯增加,除女性人群外均具有統(tǒng)計學意義,對男性和老年人群的影響的 RR明顯下降,但均不具有統(tǒng)計學意義.引入 PM10后,NO2對全人群、男性、老年人群和年齡≤15歲人群呼吸系統(tǒng)疾病日入院人數(shù)影響的影響RR明顯下降,但僅男性具有統(tǒng)計學意義,對女性和15～65歲人群的影響的RR明顯增加,且仍均具有統(tǒng)計學意義;引入SO2后,NO2對男性、老年人群和年齡≤15歲人群呼吸系統(tǒng)疾病日入院人數(shù)影響的影響RR明顯下降,但均不具有統(tǒng)計學意義,對女性和15～65歲人群的影響的RR明顯增加,且仍均具有統(tǒng)計學意義,對全人群影響變化不大.在多污染模型中,同時引入PM10、SO2和NO2,對全人群的影響仍具有統(tǒng)計學意義,對老年人群的影響不具有統(tǒng)計學意義.

圖1 空氣污染影響呼吸系統(tǒng)疾病入院人數(shù)的相對危險度(單污染模型)Fig.1 Relative risks of the association between atmospheric pollution and hospital admissions for respiratory diseases ( single-pollutant model)

表3 空氣污染影響呼吸系統(tǒng)疾病入院人數(shù)的相對危險度(多污染模型)Table 3 Relative risks of the association between atmospheric pollution and hospital admissions for respiratory diseases( multi-pollutant model)

3 討論

自20世紀90年代以來,基于廣義相加模型的時間序列分析已被國際上廣泛用于空氣污染急性暴露對各種健康效應終點的研究.時間序列方法對同一研究人群反復觀察暴露條件改變后的健康效應,因此與時間變化相關(guān)的一些變量,如年齡改變、吸煙、社會經(jīng)濟因素等,就不再成為研究空氣污染健康效應的混雜因素[13],這是時間序列的一個關(guān)鍵優(yōu)點.

國內(nèi)外研究證明,從輕微的呼吸系統(tǒng)癥狀的產(chǎn)生到心肺疾病的門診人數(shù)、入院人數(shù)和死亡率的增加都與空氣污染有著密切關(guān)系[14-18].本研究收集了2007年1月1日至2009年12月31日蘭州市大氣污染物(PM10、SO2、NO2)日均濃度和呼吸系統(tǒng)疾病日入院人數(shù)資料,探討了蘭州市大氣污染與呼吸系統(tǒng)疾病日入院人數(shù)變化的暴露-反應關(guān)系.由于研究時間較國外同類研究相對較短,而且,不同大氣污染物之間存在著明顯的共線性,具體區(qū)分某個污染物的健康效應相對困難,多污染物模型由于增加模型擬合結(jié)果的標準差[19],因而在統(tǒng)計學上的意義較低.本研究亦與國內(nèi)外報道[20-21]的空氣污染與呼吸系統(tǒng)疾病日入院人數(shù)關(guān)系的研究一致.經(jīng) GAM模型分析得知,當空氣污染物PM10、SO2和NO2濃度每增加1個IQR,全人群呼吸系統(tǒng)疾病日入院人數(shù)分別增加1.5%、4.9%和4.0%.PM10對人群的影響低于其他地區(qū)和國家的研究結(jié)果,原因可能為:一方面,雖然蘭州市的 PM10濃度較其他地區(qū)高,人群暴露水平高,但人群對其敏感程度反而較低,可能與當?shù)厝巳簩諝馕廴疽赘行缘牟顒e、人口的年齡結(jié)構(gòu)差異;另一方面,可能由于暴露-反應關(guān)系曲線在高暴露時存在近似“平臺期”,使斜率較低而導致[22].發(fā)達國家及地區(qū)的 PM10主要來自機動車尾氣排放,含有大量的二次氣溶膠,毒性高,對人體影響大,而蘭州地區(qū)顆粒物源解析表明,城市揚塵和土壤風沙貢獻率較高,其主要成分為無機礦物質(zhì),毒性弱,對人體健康影響小[22].這說明由于顆粒物特征、來源不同,顆粒物毒性存在很大差異,從而對人群健康差別很大.

從年齡與性別分層來看,年齡≤15歲的人群對空氣污染物最為敏感,PM10、SO2和 NO2每增加1個IQR,年齡≤15歲人群的呼吸系統(tǒng)疾病日入院人數(shù)分別增加4.8%、10.1%和7.8%;其次為年齡≥65歲的老年人群,PM10、SO2和NO2每增加1個 IQR,呼吸系統(tǒng)疾病日入院人數(shù)分別增加3.9%、8.5%和7.4%.原因在于:年齡≤15歲的人群生理結(jié)構(gòu)發(fā)育不夠完全,免疫系統(tǒng)不夠完善;年齡≥65歲的老年人群隨著年齡的增長,各種機能開始減退,體質(zhì)較差,機體免疫力下降,抵御不良環(huán)境的能力較差,大氣污染物中攜帶的細菌、病毒及有毒有害物質(zhì)容易進入體內(nèi),導致疾病的發(fā)生.王黎華等[23]觀察了北京市大氣污染與兒童最大呼氣流量(PEF)的關(guān)系,發(fā)現(xiàn)空氣中 NO2和 SO2的濃度與兒童 PEF的變化顯著相關(guān).與此類似,馬洪寶等[24]比較了上海大氣 SO2污染程度輕而顆粒物污染程度不同地區(qū)兒童的肺功能,結(jié)果發(fā)現(xiàn)大氣顆粒物污染可引起用力呼吸氣流(FEF)25%～75%和FEF 75%～85%的減少,提示長期暴露在顆粒物污染下,可損害兒童的肺功能.同時除PM10外,SO2和NO2對男性人群的影響較女性人群影響的 RR值大,可能是男性從事戶外作業(yè)的比例較女性高,且防護意識比女性淡薄,不良的生活習慣(如抽煙、酗酒)使男性更易受外界污染物的影響而誘發(fā)呼吸系統(tǒng)疾病.

4 結(jié)論

4.1 蘭州市空氣污染物對年齡≤15歲人群的呼吸系統(tǒng)疾病入院人數(shù)的影響最大,其次為老年人群;對男性的影響較女性明顯.

4.2 PM10滯后 2或 5d,SO2滯后3d,NO2在當天或滯后3d對呼吸系統(tǒng)疾病入院人數(shù)的RR值達到最大.

[1]Cohen A J, Ross Anderson H, Ostro B, et al. The global burden of disease due to outdoor air pollution [J]. Journal of Toxicology and Environmental Health Part A, 2005,68(13-14):1301-1307.

[2]Oberd?rster Gunter. Pulmonary effects of inhaled ultrafine particles [J]. International Archives of Occupational Environmental Health, 2001,74(1):1-8.

[3]朱 悅,楊小南,李宇斌.沈陽市某地區(qū)大氣污染與兒童呼吸系統(tǒng)疾病門診人次的關(guān)系 [J]. 環(huán)境與健康雜志, 2008,25(7):575-578.

[4]Zeghnoun A, Czernichow P, Beaudeau P, et al. Short-term effect of air pollution on mortality in the cities of Rouen and Le Havre,France,1990-1995 [J]. Archives of Environmental Health, 2001,56(4):327-335.

[5]李 寧,彭曉武,張本延.廣州市居民呼吸系統(tǒng)疾病每日死亡人數(shù)與大氣污染的時間序列分析 [J]. 華中科技大學學報, 2010,39(6):863-867.

[6]Pikhart H, Bobak M, Kriz B, et al. Outdoor air pollution of nitrogen dioxide and sulfur dioxide and prevalence of wheezing school children [J]. Epidemiology, 2000,11(2):153-160.

[7]周燕榮,曾 慶,徐 放.重慶地區(qū)大氣污染對居民呼吸系統(tǒng)疾病的影響 [J]. 重慶醫(yī)科大學學報, 1996,21(T10):118-122.

[8]衛(wèi)生部衛(wèi)生統(tǒng)計信息中心,北京協(xié)和醫(yī)院世界衛(wèi)生組織疾病分類合作中心.國際疾病分類(ICD-10)應用指導手冊 [M]. 北京:中國協(xié)和醫(yī)科大學出版社, 2001:411-435.

[9]Dominici F, McDermott A, Zeger S L, et al. On the use of generalized additive models in time-series studies of air pollution and health [J]. American Journal of Epidemiology, 2002,156(3):193-203.

[10]Schwartz J, Spix C, Touloumi G, et al. Methodological issues in studies of air pollution and daily counts of deaths or hospital admissions [J]. Journal of Epidemiology and Community Health,1996,50(S1):3-11.

[11]Akaike H. Factor Analysis and AIC [J]. Psychometric, 1987,52(3):317-332.

[12]GB3095-1996 環(huán)境空氣質(zhì)量標準 [S].

[13]闞海東,陳秉衡,賈 健.大氣污染物暴露與人群健康效應的暴露-反應關(guān)系分析 [J]. 環(huán)境與健康雜志, 2004,21(7):253-254.

[14]陳仁杰,陳秉恒,闞海東.上海市近地面臭氧污染的健康影響評價 [J]. 中國環(huán)境科學, 2010,30(5):603-608.

[15]李 寧,彭曉武,張本延,等.大氣污染與呼吸系統(tǒng)疾病日門診量的時間序列分析 [J]. 環(huán)境與健康雜志, 2009,26(12):1077-1080.

[16]常桂秋,潘小川,謝學琴,等.北京市大氣污染與城區(qū)居民死亡率關(guān)系的時間序列分析 [J]. 衛(wèi)生研究, 2003,32(6):565-568.

[17]謝 鵬,劉曉云,劉兆榮,等.珠江三角洲地區(qū)大氣污染對人群健康的影響 [J]. 中國環(huán)境科學, 2010,30(7):997-1003.

[18]王 艷,張宜升,李欣鵬.濟南城區(qū)空氣污染對呼吸道疾病門診量的影響 [J]. 中國環(huán)境科學, 2008,28(6):571-576.

[19]Peters A, D?ring A, Wichmann H E, et al. Increased plasma viscosity during an air pollution episode: a link to mortality? [J].Lancet, 1997,349(9065):1582-1587.

[20]黃文珊,胡衡生,劉 航.南寧市大氣污染與呼吸系統(tǒng)疾病的關(guān)系 [J]. 環(huán)境與健康雜志, 2004,21(6):384-386.

[21]Leitte A M, Petrescu C, Franck U, et al. Respiratory health,effects of ambient air pollution and its pollution and its medication by air humidity in Drobeta-Turnu Severin, Romania[J]. Science of the Total Environment, 2009,407(13):4004-4011.

[22]孫兆彬,李棟梁,陶 燕,等.蘭州市大氣PM10與呼吸系統(tǒng)疾病入院人數(shù)的時間序列研究[J].環(huán)境與健康雜志, 2010,27(12):1049-1052.

[23]王黎華,徐希平,周玉芬,等.空氣污染與兒童最大呼氣流速變化的關(guān)系 [J]. 環(huán)境與健康雜志, 1994,11(6):243-246.

[24]馬洪寶,洪傳潔.大氣細顆粒物污染對兒童肺功能影響的多因素分析 [J]. 中國公共衛(wèi)生學報, 1991,10(2):75-77.