賀天鋒 黃婧 李國星等

摘要：目的 評(píng)估寧波市大氣可吸入顆粒物（PM¹⁰）短期暴露造成人群死亡的歸因風(fēng)險(xiǎn)。方法 收集寧波市2011—2015年人群死亡、空氣污染和氣象資料，利用時(shí)間序列研究方法和分布滯后非線性模型，分析PM₁₀短期暴露造成人群死亡的超額危險(xiǎn)度、歸因死亡人數(shù)和人群歸因分值。結(jié)果 PM₁₀短期濃度升高與人群死亡數(shù)增加相關(guān)，效應(yīng)存在明顯的滯后性和持續(xù)性，日均濃度每升高10μg/m³，總死亡累積4 d（Lag0～3）的超額危險(xiǎn)度為0.35%（95%CI：0.16%～0.53%）。寧波市2011—2015年歸因于PM₁₀短期暴露的總死亡人數(shù)為3972例（95%CI：1877～6032），人群歸因分值為2.15%（95%CI：1.02%～3.27%）。呼吸系統(tǒng)疾病對(duì)PM¹⁰。更為敏感，其歸因死亡人數(shù)和人群歸因分值分別為1 323例（95%C/：619～1 990）和4.68%（95%CI：2.19%～7.03%）；女性和高年齡組人群的歸因風(fēng)險(xiǎn)分別高于男性和低年齡組。如果將PM¹⁰日均濃度控制在150μg/m³和100μg/m³限值以內(nèi)，人群歸因風(fēng)險(xiǎn)將分別下降25.45%和53.07%。結(jié)論 PM₁₀濃度變化對(duì)健康有較大的影響，應(yīng)立即采用有效的控制措施。

關(guān)鍵詞：空氣污染；顆粒物；歸因風(fēng)險(xiǎn)；分布滯后非線性模型 中圖分類號(hào)：R122 文獻(xiàn)標(biāo)志碼：A

隨著我國經(jīng)濟(jì)的快速發(fā)展，空氣污染問題日益嚴(yán)重，對(duì)人群健康造成了極大威脅。大氣顆粒物被認(rèn)為是我國最主要的空氣污染因子，國內(nèi)外已有大量研究證實(shí)顆粒物短期暴露與人群死亡率升高相關(guān)。但目前已有的研究多數(shù)局限于采用率的指標(biāo)評(píng)估顆粒物和死亡的關(guān)系，如相對(duì)危險(xiǎn)度或比值比等，這些指標(biāo)反映了顆粒物每單位濃度升高時(shí)的效應(yīng)，卻并未考慮到研究地區(qū)的污染強(qiáng)度和暴露人群數(shù)量。因此，無法準(zhǔn)確反映顆粒物造成的整體疾病負(fù)擔(dān)。盡管已有部分學(xué)者開展了大氣顆粒物的疾病負(fù)擔(dān)評(píng)估研究，卻并未考慮顆粒物的滯后效應(yīng)，或者缺乏對(duì)不同疾病類別和人群類別的分析。因此，對(duì)這方面還需要開展進(jìn)一步的研究。寧波市是我國長三角地區(qū)南翼經(jīng)濟(jì)中心，本研究擬利用寧波市數(shù)據(jù)，定量評(píng)估大氣可吸入顆粒物（PM₁₀）短期暴露造成的人群死亡歸因風(fēng)險(xiǎn)。

1資料與方法

1.1資料來源

死因資料來源于寧波市疾病預(yù)防控制中心，包括2011—2015年全市死因監(jiān)測個(gè)案數(shù)據(jù)。按照性別和年齡（<65歲組，≥65歲組）進(jìn)行人群分層，并依據(jù)國際疾病分類標(biāo)準(zhǔn)編碼（ICD-10）挑選出心血管疾病死亡（100-199）、呼吸系統(tǒng)疾病死亡（J00～J99）和其他疾病死亡進(jìn)行分層分析。同期大氣可吸入顆粒物（PM₁₀）的日均濃度由寧波市環(huán)境監(jiān)測中心提供，通過全市各個(gè)監(jiān)測站點(diǎn)的數(shù)據(jù)求得區(qū)域算術(shù)平均值，作為人群污染物日均暴露量。氣象資料來源于寧波市氣象局，變量包括日均氣溫、日均相對(duì)濕度和日均氣壓，同樣采用區(qū)域算術(shù)平均值作為人群氣象因子日均暴露量。

1.2研究方法

1.2.1基本假設(shè) 根據(jù)世界衛(wèi)生組織（WHO）發(fā)布的《空氣質(zhì)量準(zhǔn)則》，PM₁₀年均濃度準(zhǔn)則值為20μg/m³，低于該限值時(shí)認(rèn)為PM₁₀對(duì)健康幾乎沒有影響。同時(shí)根據(jù)已有文獻(xiàn)中通常將PM₁₀與死亡的關(guān)系設(shè)為線性，因此，本研究假設(shè)PM₁₀與死亡的暴露-反應(yīng)關(guān)系為以20μg/m³為閾值的線性函數(shù)。

1.2.2 PM₁的短期效應(yīng) 時(shí)間序列分析是當(dāng)前研究空氣污染短期影響的主要方法，其原理是在控制時(shí)間趨勢等混雜因素的基礎(chǔ)上，比較暴露因素的短期變化和健康結(jié)局的變化是否存在關(guān)聯(lián)?；灸Ｐ蜆?gòu)建如下：

利用赤池信息準(zhǔn)則（Akaike information criteri-on，AIC）確定交叉基函數(shù)中的自由度，選擇相對(duì)危險(xiǎn)度（RR）和超額危險(xiǎn)度（ER）作為PM₁₀短期效應(yīng)評(píng)價(jià)指標(biāo)，本研究中ER表示為PM₁₀每升高10μg/m³造成死亡增加的百分比，計(jì)算方法如下：ER=（RR-1）×100%=[exp（β）-1]×100%

1.2.3 PM₁₀的歸因風(fēng)險(xiǎn) 歸因風(fēng)險(xiǎn)的基礎(chǔ)指標(biāo)是歸因人數(shù)（attribute number，AN）和人群歸因分值（attribute fraction，AF）。本研究中AF表示由于PM₁₀短期暴露造成的超額死亡人數(shù)占總死亡的比例，也可理解為如果人群不再暴露于PM₁₀時(shí)，相應(yīng)死亡減少的數(shù)量占總死亡的比例。兩者的計(jì)算方法分別如下：

根據(jù)WHO《空氣質(zhì)量準(zhǔn)則》中不同過渡時(shí)期的目標(biāo)值，本研究分別選擇150、100、75、50μg/m³作為PMI。日均濃度控制限值，并分析了當(dāng)超出該限值時(shí)所對(duì)應(yīng)的PM₁₀短期暴露對(duì)總死亡影響的歸因風(fēng)險(xiǎn)大小。

1.2.4敏感性分析 由于WHO和環(huán)境空氣質(zhì)量標(biāo)準(zhǔn)（GB 3095—2012）均將50μg/m³設(shè)為PM¹⁰日均濃度準(zhǔn)則值，同時(shí)部分研究也將0μg/m³作為顆粒物危害評(píng)估基線。因此，本研究也采用了50μg/m³和0μg/m³作為PM₁₀閾值進(jìn)行敏感性分析。為驗(yàn)證模型的穩(wěn)定性，本研究改變了模型中PM¹⁰最長滯后時(shí)間（5～10 d）、PM₁₀滯后時(shí)間的自由度（3～5）、時(shí)間變量的自由度（6～8）、溫度的自由度（3～6）、濕度和氣壓的自由度（3～6）進(jìn)行敏感性分析，結(jié)果均使用ER表示。

1.3統(tǒng)計(jì)學(xué)分析

研究中采用R3.1.1軟件進(jìn)行統(tǒng)計(jì)分析，分布滯后非線性模型采用“dlnm”程序包。檢驗(yàn)水準(zhǔn)α=0.05。

2結(jié)果

2.1一般情況

寧波市2011—2015年共監(jiān)測到死亡人數(shù)184 542例，其中心血管疾病共死亡52415例，呼吸系統(tǒng)疾病共死亡28 302例，其余人群死亡和氣象資料的一般情況見表1。研究期間PM₁₀日均濃度為（80.48±49.69）μg/m³，其中超過WHO《空氣質(zhì)量準(zhǔn)則》年均濃度準(zhǔn)則值（20μg/m³）共1814 d，占總天數(shù)的99.34%。PM₁₀日均濃度的時(shí)間分布趨勢呈現(xiàn)明顯的季節(jié)性，冬季污染嚴(yán)重，夏季污染較輕。見圖1。

2.2 PM₁₀的短期效應(yīng)（圖2）

PM₁₀對(duì)人群死亡的影響具有明顯的滯后性。在滯后1～3 d時(shí)PM₁₀和總死亡的關(guān)聯(lián)有統(tǒng)計(jì)學(xué)意義（P<0.05），滯后1 d時(shí)效應(yīng)最強(qiáng)；PM₁₀每升高10μg/m³，累積4 d（Lag0～3）可造成總死亡增加0.35%（95%CI：0.16%～0.53%）。根據(jù)研究假設(shè)，PM₁₀與總死亡累積4 d的暴露一反應(yīng)關(guān)系如圖3所示，閾值（20μg/m³）以下PM₁₀效應(yīng)為零，超過閾值后死亡的相對(duì)危險(xiǎn)度呈線性上升。

2.3 PM₁₀短期暴露的歸因風(fēng)險(xiǎn)

當(dāng)以20μg/m³作為PM₁₀閾值時(shí)，2011—2015年寧波市歸因于PM₁₀短期暴露的總死亡人數(shù)為3972例（95%CI：1 877～6 032），人群歸因分值為2.15%（95%CI：1.02～3.27%）。不同人群組別中，女性和I>65歲年齡組的歸因風(fēng)險(xiǎn)分別高于男性和<65歲年齡組，兩者的人群歸因分值分別為2.91%和2.84%。呼吸系統(tǒng)疾病的歸因風(fēng)險(xiǎn)高于心血管疾病，歸因死亡人數(shù)和人群歸因分值分別為1323例和4.68%。其他疾病歸因于PM₁₀短期暴露的死亡人數(shù)共為1505例，人群歸因分值為1.45%。見表2。

2011—2015年歸因于PM₁₀短期暴露的總死亡人數(shù)分別為795、671、1145、672、689例，其時(shí)間分布趨勢與PM₁₀日均濃度分布一致。冬季歸因死亡人數(shù)較多，共為1 663例，占總歸因人數(shù)的41.86%；夏季人數(shù)較少，共457例，占總歸因人數(shù)的11.50%。見圖4。

寧波市2011—2015年P(guān)M₁₀日均濃度超過150μg/m³共148 d，占總天數(shù)的8.11%，其所對(duì)應(yīng)的歸因人數(shù)為1011例，占總歸因人數(shù)的25.45%；超過100μg/m³共444 d，占總天數(shù)的24.32%，所對(duì)應(yīng)的歸因死亡人數(shù)為2108例，占總歸因人數(shù)的53.07%；超過50μg/m³共1313 d，占總天數(shù)的71.91%，所對(duì)應(yīng)的歸因死亡人數(shù)為3671例，占總歸因人數(shù)的92.42%。見表3。

2.4敏感性分析

當(dāng)改變PM₁₀閾值、時(shí)間變量的自由度、氣象變量的自由度等參數(shù)時(shí)，其效應(yīng)均未發(fā)生較大改變，結(jié)果顯示模型穩(wěn)定性較好。見表4。

3討論

空氣污染造成的人群健康風(fēng)險(xiǎn)已成為當(dāng)前政府及公眾高度關(guān)注的問題，及時(shí)開展本地區(qū)空氣污染對(duì)人群健康的影響研究，為決策者制定相應(yīng)政策提供理論依據(jù)，具有重要的公共衛(wèi)生學(xué)意義。本研究在以往分析顆粒物每單位濃度的效應(yīng)基礎(chǔ)上，結(jié)合寧波市本地空氣污染濃度和人群暴露數(shù)量，定量評(píng)估了寧波市近幾年歸因于PM₁₀短期暴露的人群死亡數(shù)量和人群歸因分值。結(jié)果發(fā)現(xiàn)，PM₁₀每升高10μg/m³，累積4 d（Lag0～3）可造成總死亡增加0.35%（95%CI：0.16%～0.53%），此效應(yīng)大小與中國其他城市的研究結(jié)果接近。寧波市2011—2015年歸因于PM₁₀短期暴露的死亡人數(shù)共3972例，占總死亡人數(shù)的2.15%，表明空氣污染已經(jīng)給人群健康造成了較大損失。

目前認(rèn)為敏感人群數(shù)小時(shí)或數(shù)天的高濃度暴露即可引起心肌梗死或心力衰竭，與以往文獻(xiàn)一致。本研究發(fā)現(xiàn)，呼吸系統(tǒng)疾病和心血管疾病對(duì)顆粒物污染均較其他疾病更為敏感，兩者由于PM₁₀短期暴露造成的人群歸因分值均高于其他疾病。其中呼吸系統(tǒng)疾病受到的影響更強(qiáng)，大約每20個(gè)呼吸系統(tǒng)疾病死亡病例中就有1個(gè)是因PM₁₀短期暴露造成的。已有較多的證據(jù)表明，老年人是空氣污染的敏感人群，本研究結(jié)論與以往一致。同時(shí)，研究發(fā)現(xiàn)女性受到PM₁₀的影響比男性更強(qiáng)，兩者的歸因分值分別為2.91%和1.56%，此機(jī)制尚不清楚，可能與兩者生活習(xí)慣、肺功能、內(nèi)分泌改變等差異有關(guān)。

污染物濃度差異同樣會(huì)導(dǎo)致整體疾病負(fù)擔(dān)不同，寧波市2012年和2013年的PM₁₀年均濃度分別為71.46μg/m³和107.43μg/m³，造成的歸因死亡人數(shù)分別為671和1145例，相差近1倍。據(jù)此類推，由于中國各地顆粒物濃度和暴露人數(shù)不同，造成的疾病負(fù)擔(dān)會(huì)存在較大差異，如李湉湉等的研究發(fā)現(xiàn)，2013年1月北京、上海、廣州、西安PM_2.5日均濃度分別為240、83、94、210μg/m³，歸因死亡人數(shù)分別為725、296、310、85例。與其他國家相比，中國具有眾多的暴露人口和較高的顆粒物濃度，因此整體歸因疾病負(fù)擔(dān)會(huì)更加嚴(yán)重。本研究同時(shí)發(fā)現(xiàn)，如果采取及時(shí)有效措施降低空氣污染濃度，將會(huì)獲得較大的人群健康收益。如將PMI。日均濃度均控制在150μg/m³以下時(shí)，即能夠降低25.45%的歸因風(fēng)險(xiǎn)；而當(dāng)濃度控制在100μg/m³以下時(shí)，歸因風(fēng)險(xiǎn)將下降一半以上。

本研究方法建立在一系列的假設(shè)基礎(chǔ)上，其結(jié)果仍存在一些局限性。首先研究中只評(píng)估了PM₁₀短期暴露的歸因風(fēng)險(xiǎn)，而國外隊(duì)列研究顯示顆粒物的長期效應(yīng)可能更強(qiáng)。因此，本研究發(fā)現(xiàn)的歸因風(fēng)險(xiǎn)僅僅是PM₁₀疾病負(fù)擔(dān)中的一小部分。此外，以往文獻(xiàn)發(fā)現(xiàn)高溫可能會(huì)增強(qiáng)PM₁₀的效應(yīng)，如果考慮PM₁₀與氣溫的交互作用，夏季的歸因風(fēng)險(xiǎn)可能比研究中估計(jì)的更強(qiáng)。但與以往研究相比，本研究不再局限于分析PM₁₀單位濃度的效應(yīng)，而是定量評(píng)估了PM₁₀短期暴露的整體疾病，并比較了不同疾病和人群的差異，研究結(jié)果能夠?yàn)闆Q策者制定相關(guān)環(huán)境防控政策提供理論依據(jù)，具有一定的創(chuàng)新性和實(shí)用性。