王德征，江國虹，顧清，張輝，徐忠良、宋桂德，張穎、沈成鳳

采用時(shí)間序列泊松回歸分析天津市大氣污染物對心腦血管疾病死亡的急性影響*

王德征，江國虹**，顧清，張輝，徐忠良、宋桂德，張穎、沈成鳳

目的：了解天津市大氣污染物對居民心腦血管疾病死亡的急性影響，探索早期預(yù)警方法，為心腦血管疾病的預(yù)防控制提供依據(jù)。

方法：采用天津市疾病預(yù)防控制中心收集的居民全死因監(jiān)測數(shù)據(jù)， 氣象資料和大氣污染物資料來源于天津市氣象局和天津市環(huán)境監(jiān)測中心。采用時(shí)間序列的泊松回歸GAM進(jìn)行天津市每日大氣污染物與居民心腦血管疾病死亡危險(xiǎn)度分析， 同時(shí)控制氣象因素、長期趨勢、星期幾效應(yīng)以及人口數(shù)等混雜因素的影響，進(jìn)行單污染物和多污染物分析。

結(jié)果：單因素GAM分析顯示大氣中SO2、NO2、PM10效應(yīng)在當(dāng)日達(dá)到最大，風(fēng)險(xiǎn)RR值分別增加1.13% [95%CI（0.76～1.51） %]、 0.78% [95%CI（ 0.41～1.15）%]和 0.61% [95%CI（0.51～0.71）%]；SO2、NO2、PM10日均濃度每升高10 μg/m3，0～5日后心腦血管疾病死亡的風(fēng)險(xiǎn)RR值分別增加0.70% [95%CI（0.47～0.94）]%， 0.51% [95%CI（0.27～0.74）%]和 0.16% [95%CI（0.06～0.27）%]；多因素分析顯示，SO2、NO2、PM10日均濃度每增加10 μg/m3，心腦血管疾病死亡的風(fēng)險(xiǎn)RR值分別增加0.77% [95%CI（0.58～0.97）%]，0.41% [95%CI（0.05～0.78）%]和 0.38% [95%CI（0.12～0.64）%]。

結(jié)論：天津市大氣污染物能增加人群心腦血管疾病死亡風(fēng)險(xiǎn)，可能對心腦血管疾病急性發(fā)作起到誘導(dǎo)作用，有必要開展大氣污染物治理，建立預(yù)警系統(tǒng)，合理分配衛(wèi)生資源，降低心腦血管疾病死亡危害。

心腦血管疾?。凰劳雎?；廣義可加模型; 大氣污染物; 急性影響

Methods: The Mortality for patients with cardio cerebral vascular disease from 2001-01 to 2009-12 was from Tianjin Centers for Disease Control and Prevention, the meteorological and air pollution data were from Tianjin meteorological bureau and Tianjin environmental monitoring station respectively. The Single and multiple generalized additive model (GAM) extended poisson regression analysis was performed to calculate the relationship between air pollution and cardio-cerebral vascular disease mortality by controlling the time trends, weather, the day of week and air population.

Results: Results: Single GAM analysis indicated that when the effect of SO2, NO2and PM10in the air reached the maximum in the day, the RR values for the mortality in patients with cardio-cerebral vascular disease increased 1.13%

[95%CI (0.76-1.51)%], 0.78% [95%CI ( 0.41-1.15)%] and 0.61% [95%CI ( 0.51-0.71)%] respectively; when the average concentration of SO2, NO2and PM10 increasing 10μg/m3per day, after 0-5 days, the RR values for the mortality elevated 0.70% [95%CI (0.47-0.94)%], 0.51% [95%CI (0.27-0.74)%] and 0.16% [95%CI (0.06-0.27)%] respectively. Multiple GAM analysis presented that when SO2, NO2and PM10increasing 10 μg/m3per day, the RR values for the mortality elevated 0.77% [95%CI (0.58-0.97)%], 0.41% [95%CI (0.05-0.78) %] and 0.38% [95%CI (0.12-0.64%)%] respectively. Conclusion: The air pollution could increase the mortality risk in patients with cardio-cerebral vascular disease, it is necessary to establish the prevention system in order to decrease the mortality risk in those patients.

(Chinese circulation Journal, 2014,29:453.)

在全球范圍內(nèi)， 60%的由心腦血管疾病導(dǎo)致的負(fù)擔(dān)都發(fā)生在發(fā)展中國家[1]。我國城市和農(nóng)村冠心病在全死因順位中居第三或第四，腦卒中在全死因順位中居第二或第三[2]，是導(dǎo)致我國居民死亡的主要原因之一。心腦血管疾病是一種多因素疾病，是遺傳因素和環(huán)境因素共同作用的結(jié)果[3]，特別是隨著工業(yè)化、城市化和人口快速增長，人類活動導(dǎo)致的環(huán)境變化特別是可吸入顆粒物（PM10）、二氧化硫（SO2）、二氧化氮（ NO2）、一氧化碳（CO） 、臭氧（O3）等的增加對心腦血管疾病發(fā)病、死亡的影響已引起人們廣泛關(guān)注。天津市是典型的北方工業(yè)城市, 探討大氣污染物對天津市居民心腦血管疾病死亡的急性影響，探索早期預(yù)警方法，為心腦血管疾病的預(yù)防控制提供依據(jù)。

1 材料和方法

材料：采用2001-01-01至2009-12-31天津市疾病預(yù)防控制中心所收集的天津市居民全死因監(jiān)測數(shù)據(jù)。心腦血管疾病診斷分類按照“疾病和有關(guān)健康問題的國際統(tǒng)計(jì)分類 （ICD） ”：ICD-9: 390-459；ICD-10: I00-I99。以實(shí)際死亡日期為分析時(shí)間點(diǎn)。采用天津市氣象局和天津市環(huán)境監(jiān)測中心收集的天津市同時(shí)段氣象指標(biāo)（包括每日平均溫度、平均相對濕度、平均大氣壓、平均風(fēng)速等，以及大氣污染物資料包括8個(gè)國控測點(diǎn)的可吸入顆粒物PM10、SO2及NO2的每日平均值，每日24小時(shí)24次監(jiān)測值的算術(shù)平均值）。對大氣污染物資料、氣象資料中的缺失值，采用前后各兩天數(shù)值的中位數(shù)為缺失值當(dāng)天的數(shù)值，缺失值填補(bǔ)比例分別為0.06%和0.91%。

統(tǒng)計(jì)方法：統(tǒng)計(jì)軟件為SPSS11.5軟件和SAS8.1。對于總?cè)丝趤碚f，每日心腦血管疾病死亡為小概率事件，作為一種時(shí)間序列資料，其實(shí)際分布近似泊松分布。應(yīng)采用時(shí)間序列的泊松回歸廣義可加模型（GAM）。具體模型為:

[Yi ：觀察日；i 當(dāng)天的死亡人數(shù)， E（ Yi ） ：觀察日i日死亡數(shù)的預(yù)期值；α：截距；X：對應(yīng)變量產(chǎn)生線性影響指示變量；β：通過回歸模型估計(jì)的指示變量系數(shù)： f ：非參數(shù)自然三次光滑樣條函數(shù)；Z ：對應(yīng)變量發(fā)生非線性影響的變量，∑：加和]。

通過擬合單因素GAM并進(jìn)行偏差性檢驗(yàn)，同時(shí)考慮氣象因素、人口結(jié)構(gòu)、星期幾效應(yīng)以及長期趨勢等混雜因素，分析大氣單污染物對心腦血管疾病死亡的影響，并用相對危險(xiǎn)度（RR）評價(jià)危害的大小?？紤]到污染物對人體危害可能的滯后效應(yīng)，選取遲滯日（滯后天數(shù)）1～5天，分析大氣污染物對心腦血管疾病死亡影響的變化情況。為檢驗(yàn)?zāi)Ｐ偷姆€(wěn)定性,采用廣義交叉確認(rèn)（ GCV） 法進(jìn)行不同統(tǒng)計(jì)學(xué)模型的擬合優(yōu)度檢驗(yàn)[4]。

對單因素分析有意義、相對危險(xiǎn)度較大且模型擬合優(yōu)度較好的污染物、滯后天數(shù)等，并考慮可能的混雜因素，建立共同擬合多因素GAM，并進(jìn)行偏差性檢驗(yàn)和擬合優(yōu)度檢驗(yàn)。

以2000年世界標(biāo)準(zhǔn)人口計(jì)算年齡、性別標(biāo)化死亡率，采用趨勢卡方檢驗(yàn)對心腦血管疾病的年死亡率進(jìn)行趨勢分析。

采用Spearman等級相關(guān)分析大氣污染物、氣象因素指標(biāo)與心腦血管疾病死亡的相關(guān)性。

2 結(jié)果

2001～2009年天津市每日心腦血管疾病死亡、氣象因素和大氣污染物濃度分布：心腦血管疾病死

亡、氣象因素和大氣污染物日間分布存在明顯的周期性（圓分布檢驗(yàn)，P＜0.01），且均不符合正態(tài)分布（P＜0.05）。2001～2009年天津市每日中位心腦血管疾病死亡數(shù)為77例；氣象因素中，中位氣溫、濕度、氣壓、風(fēng)速分別為16 ℃、69 %、1 015 hPa、10 Km/h；大氣污染物中的SO2、NO2和PM10中位濃度分別為4.90 10 μg/m3、 4.68 10 μg/m3和10.32 10 μg/m3。表1

2001～2009年天津市每日心腦血管疾病死亡數(shù)與氣象因素和大氣污染物因素Spearman相關(guān)分析：SO2、NO2每日濃度同心腦血管疾病每日死亡數(shù)間呈正相關(guān)（r 分別為0.479、0.250），差異有統(tǒng)計(jì)學(xué)意義（ P＜0.01）；PM10每日濃度同心腦血管疾病每日死亡數(shù)間呈負(fù)相關(guān)（r=-0.023），但差異無統(tǒng)計(jì)學(xué)意義（ P=0.182）。每日氣象因素中，每日溫度、每日濕度與心腦血管疾病每日死亡數(shù)呈負(fù)相關(guān)（r 分別為-0.572和-0.164），差異有統(tǒng)計(jì)學(xué)意義（ P＜0.01）；每日氣壓、每日風(fēng)速與心腦血管疾病每日死亡數(shù)呈正相關(guān)（r 分別為0.479和0.063），差異有統(tǒng)計(jì)學(xué)意義（P＜0.01）。同時(shí)每日大氣污染物、氣象因素內(nèi)部及相互之間呈現(xiàn)不同程度的相關(guān)關(guān)系，具體結(jié)果見表2。

表1 2001～2009年天津市每日心腦血管疾病死亡例數(shù)、氣象因素和大氣污染物濃度分布

表2 2001～2009年天津市每日心腦血管疾病死亡數(shù)與氣象因素和大氣污染物因素Speaman相關(guān)分析

大氣污染物對不同遲滯日心腦血管疾病死亡的單因素GAM分析：經(jīng)過可能的混雜因素調(diào)整后大氣污染物日均濃度對當(dāng)日及不同滯天數(shù)的心腦血管疾病死亡的泊松回歸GAM分析顯示: SO2、NO2、PM10日均濃度每升高10 μg/m3，對心腦血管疾病死亡的風(fēng)險(xiǎn)在當(dāng)日達(dá)到最大，其RR值及95%可信區(qū)間分別增加 1.13% [95%CI（0.76～1.51 ） %]、0.78% [95%CI（ 0.41～1.15）%]和 0.61% [95%CI（0.51～0.71）%]； 5日后SO2、NO2、PM10日均濃度每升高10 μg/m3，心腦血管疾病死亡的風(fēng)險(xiǎn)RR值分別增加0.70% [95%CI（0.47～0.94 ） %], 0.51% [95%CI（0.27～0.74）%]和 0.16% [95%CI（0.06～0.27）%]，其滯后作用隨著時(shí)間延長而逐漸減小。表3

大氣污染物及相關(guān)因素對心腦血管疾病死亡的多因素GAM分析：由污染物對心腦血管疾病死亡影響的單因素分析，選取當(dāng)日（Lag0）SO2、NO2、PM10及可能的混雜因素進(jìn)行多因素泊松GAM分析，SO2、NO2、PM10日均濃度每增加10 μg/m3，心腦血管疾病死亡的風(fēng)險(xiǎn)RR值及95% IC 分別增加0.77%（0.58～0.97）%，0.41% （0.05～0.78）%和 0.38%（0.12～0.64）%。與此同時(shí)，低溫、低氣壓是當(dāng)日（Lag0）心腦血管疾病死亡的危險(xiǎn)因素（P＜0.05）; 隨著每年人口數(shù)的增加，心腦血管疾病死亡數(shù)也在增加，P=0.0000。結(jié)果見表4。

表3 大氣污染物對不同遲滯日心腦血管疾病死亡的單因素GAM分析

表4 大氣污染物及相關(guān)因素對心腦血管疾病死亡的多因素GAM分析

3 討論

本研究結(jié)果表明，天津市2001至2009年SO2、NO2和PM10平均水平，高于國家環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)[5]二級水平，其中NO2平均水平低于北京和上海，SO2平均水平高于北京和上海，PM10平均水平在上海和北京之間[6]。天津市心腦血管疾病死亡、大氣污染物和氣象因素日間分布存在明顯的相關(guān)性，低氣溫能增加心腦血管疾病死亡的風(fēng)險(xiǎn)。Scarborough 等[7]對英格蘭1999～2004年冠心病死亡數(shù)據(jù)的分析表明，人群冠心病的死亡率隨著溫度的降低和日照時(shí)間的縮短而增加。此前關(guān)于天津市心腦血管病死亡周期性的時(shí)間序列分析[8]也表明，冠心病、腦卒中死亡率在冬季，特別是1月份呈現(xiàn)高水平。天津市四季分明, 采暖期污染物的水平明顯高于非采暖期污染物水平。Choi等[9]研究認(rèn)為大氣污染物與心腦血管疾病的關(guān)系, 必須要考慮氣象因素的影響，其中對血壓產(chǎn)生影響的大氣污染物因季節(jié)不同而不同, 在夏季主要是PM10和NO2, 而在冬季則是SO2和O3。

本研究大氣污染物單因素分析與許多研究結(jié)果相似，多因素分析模型同時(shí)引入各污染物，危險(xiǎn)度分析結(jié)果穩(wěn)健[10,11]。

許多研究認(rèn)為大氣污染物對人體健康的急性影響存在滯后效應(yīng) ，Analitis等[12]對歐洲29個(gè)城市的Meta分析結(jié)果說明大氣污染物對心腦血管系統(tǒng)疾病的滯后效應(yīng)可能比對呼吸系統(tǒng)疾病明顯。本研究選取了大氣污染物當(dāng)日及滯后1～5天的心腦血管疾病數(shù)據(jù)進(jìn)行分析，結(jié)果顯示心腦血管疾病死亡的滯后作用隨著時(shí)間延長而逐漸減小。

目前，大氣污染物對人體心血管系統(tǒng)的影響和機(jī)制研究是初步的和很不深入的。有研究表明，大氣污染物對心血管系統(tǒng)健康的影響中, 氧化應(yīng)激和炎性反應(yīng)是至關(guān)重要的, 可以引起動脈粥樣斑塊、血管內(nèi)皮、血小板和心律的變化[13]，大氣污染物暴露會損傷微血管功能，進(jìn)而引起NO介導(dǎo)的交感神經(jīng)興奮[14]。提示大氣污染物可能對心血管病急性發(fā)病起到誘導(dǎo)作用。

綜上所述，本研究發(fā)現(xiàn)SO2、NO2和PM10日平均濃度短期升高是人群心腦血管疾病死亡的危險(xiǎn)因素，有必要開展大氣污染物治理，合理分配衛(wèi)生資源，降低心腦血管死亡危害。

[1]Kim AS, Johnston SC. Global variation in the relative burden of stroke and ischemic heart disease. Circulation, 2011, 124: 314-323.

[2]中華人民共和國衛(wèi)生部. 2010中國衛(wèi)生統(tǒng)計(jì)年鑒. 北京: 中國協(xié)和醫(yī)科大學(xué)出版社, 2010.

[3]高閱春, 何繼強(qiáng), 姜騰勇, 等. 冠心病患者冠狀動脈病變嚴(yán)重程度與冠心病危險(xiǎn)因素的相關(guān)分析. 中國循環(huán)雜志, 2012, 27: 178-181.

[4]馮國雙, 陳景武, 張國英. 用GAM 程序擬合光滑樣條非參數(shù)回歸.數(shù)理醫(yī)藥學(xué)雜志, 2005, 18: 403-406.

[5]中華人民共和國環(huán)境保護(hù)部. 《環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)》（GB3095-2012）. 北京：中國環(huán)境科學(xué)出版社, 2012.

[6]宋桂香, 江莉莉, 陳國海, 等. 上海市大氣氣態(tài)污染物與居民每日死亡關(guān)系的時(shí)間序列研究. 環(huán)境與健康雜志, 2006, 23: 390-393.

[7]Scarborough P, Allender S, Rayner M, et al. Contribution of climate and air pollution to variation in coronary heart disease mortality rates in England. PLoS One, 2012, 7: e32787. Epub 2012 Mar 12.

[8]王德征, 江國虹, 宋桂德, 等. 求和自回歸滑動平均模型結(jié)合圓分布法分析腦卒中死亡率動態(tài)規(guī)律. 中華流行病學(xué)雜志, 2009, 30：82-84.

[9]Choi JH, Xu QS, Park SY, et al . Seasonal variation of eff ect of air pollution on blood pressure. J Epidemiol Community Health, 2007, 61: 314-318.

[10]Yorifuji T, Kawachi I, Sakamoto T, et al. Associations of outdoor air pollution with hemorrhagic stroke mortality. J Occup Environ Med, 2011, 53: 124-126.

[11]Scaife A, Barclay J, Hillis GS, et al. Lack of effect of nitrogen dioxide exposure on heart rate variability in patients with stable coronary heart disease and impaired left ventricular systolic function. Occup Environ Med, 2012, 69: 587-91.

[12]Analitis A, Katsouyanni K, Dimakopoulou K, et al. Short-term effects of ambient particles on cardiovascular and respiratory mortality. Epidemiology, 2006, 17: 230-233.

[13]黃新金, 李菊香, 顏素娟, 等. 冠狀微血管痙攣性心絞痛患者血內(nèi)源性一氧化氮合酶抑制物的變化. 中國循環(huán)雜志, 2007, 22: 100-102.

[14]Knuckles TL, Stapleton PA, Minarchick VC, et al. Air pollution particulate matter collected from an Appalachian mountaintop mining site induces microvascular dysfunction. Microcirculation, 2013, 20: 158-169.

A Time-series Study for Acute Effect of Air Pollution on Mortality in Patients With Cardio-cerebral Vascular Disease in Tianjin City

WANG De-zheng, JIANG Guo-hong, GU Qing, ZHANG Hui, XU Zhong-liang, SONG Gui-de, ZHANG Ying, SHEN Cheng-feng. Tianjin Center for Diseases Control and Prevention, Tianjin (300011), China

Objective: To explore the acute effect of air pollution on mortality for patients with cardio cerebral vascular disease and to provide the basis for disease prevention and control.

Cardio-cerebral vascular disease; Mortality; Generalized additive model; Air pollution; Acute effect

2013-09-16）

（助理編輯：曹洪紅）

國家973計(jì)劃課題（2011CB503801）

300011 天津市，天津市疾病預(yù)防控制中心

王德征 主治醫(yī)師 碩士 主要從事慢性病流行病學(xué) Email：wangdezheng1001@163.com**為共同第一作者 通訊作者：顧清

Email：guqing315@126.com

R541

A

1000-3614（2014）06-0453-05

10.3969/j.issn.1000-3614.2014.06.015