姜綿峰, 葉春明,*, 盛真真, 盛小星, 花宇

1. 上海理工大學(xué)管理學(xué)院, 上海 200090 2. 廣東醫(yī)學(xué)院研究生學(xué)院, 廣東湛江 524000

上海市霧霾健康經(jīng)濟損失風險評估

姜綿峰1, 葉春明1,*, 盛真真1, 盛小星1, 花宇2

1. 上海理工大學(xué)管理學(xué)院, 上海 200090 2. 廣東醫(yī)學(xué)院研究生學(xué)院, 廣東湛江 524000

為量化評估霧霾顆粒物PM2.5的區(qū)域人群健康損害, 基于2006—2015年上海市地方統(tǒng)計數(shù)據(jù), 使用暴露反應(yīng)關(guān)系函數(shù)計算上海市霧霾健康經(jīng)濟損失及占同年GDP比例, 采用Bootstrap信息擴散綜合模型分析了健康損失在不同風險水平下的超越概率、期望水平和重現(xiàn)周期, 并利用多項式擬合預(yù)測了未來5年的健康損失值以及占GDP比例。結(jié)果顯示上海市PM2.5國家標準健康經(jīng)濟損失最高為2013年的120.89億元, 占GDP比例最大為2015年的4.60‰; 目前上海市PM2.5的健康經(jīng)濟損失超過80億元的概率為56.62%, 近10年平均健康損失約94.64億元, 平均占GDP比例約5.42‰, 5年內(nèi)會出現(xiàn)144億元左右的健康損失, 預(yù)計2020年損失值達到261.85億元。上海市PM2.5經(jīng)濟健康損失出現(xiàn)快速波動上升趨勢, 需要有效措施控制霧霾污染同時設(shè)置預(yù)警機制應(yīng)對潛在的突發(fā)事件。

上海; 霧霾; 經(jīng)濟損失; 占GDP比例; 超越概率

1 前言

2000年后經(jīng)過將近5輪環(huán)保三年行動計劃, 上海市環(huán)境空氣質(zhì)量優(yōu)良率總體呈上升趨勢, 然而新型大氣污染問題凸顯, 其中主要由細顆粒物(PM2.5)造成的霧霾近幾年成為上海市最嚴重的環(huán)境問題[1]。細顆粒物(PM2.5)的物理性質(zhì)特殊且化學(xué)成分復(fù)雜, 不僅會引起可見度下降還能引發(fā)各種急慢性健康問題,已成為治理城市霧霾污染的重點和關(guān)鍵[2–3]。

針對霧霾污染的健康危害, 國內(nèi)外學(xué)者展開了不同層面的經(jīng)濟損失評估。穆泉等[4]學(xué)者計算了2013年1月全國霧霾直接經(jīng)濟損失以及健康經(jīng)濟損失的比例。QUAH E等[5]估算了霧霾對新加坡人口死亡和發(fā)病效應(yīng)的經(jīng)濟損失和占到當年GDP比例。MATUS K等[6]學(xué)者研究了空氣質(zhì)量、人口、收入水平與霧霾污染邊際損失的關(guān)系以及污染損失占GDP比例的變化趨勢。地區(qū)城市方面, 陳仁杰等[7]估算出2006年P(guān)M10對我國113個主要城市的健康經(jīng)濟總損失, 劉曉云等[8]學(xué)者分別就國家和國際兩個標準計算了PM10對珠三角9個城市的急性健康損失及占GDP比例。黃德生等[9]學(xué)者研究了京津冀地區(qū)PM2.5的健康和經(jīng)濟改善效益。潘小川等[10]學(xué)者估算了的北京上海廣州等中心城市2010年P(guān)M2.5引起的早死人數(shù)和經(jīng)濟損失。謝元博等[11]學(xué)者研究了2013年1月 PM2.5造成的北京市健康損失, 陳依等[12]分析了10年間南京市PM10的健康損失變化特征, 侯青等[13]用 meta分析定量評估了蘭州市 2002—2009年間PM10健康經(jīng)濟損失, KAN等[14]分析2001年上海市PM10的經(jīng)濟損失及占 GDP比例。上述研究利用歷史資料計算了霧霾導(dǎo)致的直接健康經(jīng)濟損失, 污染因子以 PM10為主, 健康終端以人口早逝為主, 部分研究使用舊的國家空氣標準和國外參數(shù), 對健康經(jīng)濟損失的概率分布規(guī)律沒有深入分析, 缺少對不同損失規(guī)模的預(yù)先評估。

我國空氣質(zhì)量標準滯后于國際相關(guān)標準, 霧霾顆粒物 PM2.5的監(jiān)測數(shù)據(jù)缺乏共同導(dǎo)致霧霾健康損失這一方面缺乏系統(tǒng)的風險評估。作為全國經(jīng)濟發(fā)展改革中心城市, 上海市人口對土地、能源和交通需求巨大導(dǎo)致資源生態(tài)壓力指數(shù)居高不下[15–16], 霧霾對居民的健康影響急需預(yù)先的風險評估。依據(jù)最新空氣質(zhì)量標準和計算方法, 定量分析近10年上海市 PM2.5的健康經(jīng)濟影響, 并以此為基礎(chǔ)計算不同損失規(guī)模的概率分布, 就出現(xiàn)周期和未來發(fā)展趨勢做出詳細分析, 旨在為上海市霧霾污染治理和能源產(chǎn)業(yè)改革提供決策參考。

2 研究方法

在洪澇、氣象和地震等災(zāi)害風險評估主流方法中, 概率統(tǒng)計方法和信息擴散模型方法都需要參考歷史數(shù)據(jù)進行分析[17]。目前對于霧霾健康損失缺乏可信的官方資料, 尤其是各省市的健康損失資料較少。通過流行病學(xué)的暴露反應(yīng)函數(shù)計算上海市歷年的健康經(jīng)濟損失, 在既得歷年損失值的基礎(chǔ)上, 根據(jù)信息擴散模型深入研究健康損失的概率分布。

2.1 暴露反應(yīng)函數(shù)

當前大氣污染健康損失計算中, 健康損失首先通過流行病學(xué)的暴露反應(yīng)關(guān)系模型來得到健康效應(yīng),然后采用量化方法得到具體損失價值, 常用方法有醫(yī)療費用法、修改的人力資本法、支付意愿調(diào)查法、統(tǒng)計生命價值和影響途徑分析法[12,18–20]。暴露反應(yīng)函數(shù)確定了人群在不同污染物濃度下某種健康問題的實際發(fā)病率E(公式1), 流行病學(xué)依據(jù)既定的E值和暴露反應(yīng)函數(shù)逆向求取污染因子的基期發(fā)病率E0,通過基期E0與實際發(fā)病率E的差值得到污染因子導(dǎo)致的某種健康問題的發(fā)病率, 再結(jié)合暴露總?cè)丝赑ed以及各類健康效應(yīng)的經(jīng)濟價值 HCmu可得到暴露人群的總體健康損失ECal(公式2)。歷年的健康經(jīng)濟損失將作為進一步統(tǒng)計分析風險狀況的基礎(chǔ)。

其中,β為不同健康問題頻率與污染因子濃度的相關(guān)系數(shù), C為污染因子實際濃度, C0為污染因子流行病學(xué)濃度閾值, E為污染因子在濃度C下人群健康效應(yīng), E0為污染因子在濃度C0下人群健康效應(yīng)。

其中, Ped為暴露人口數(shù), ECal為污染因子造成的居民健康損失, HCmu表示暴露人群的各類健康問題的經(jīng)濟價值。

2.2 基于Bootstrap的信息擴散模型

在災(zāi)害風險評估方法中, 信息擴散模型具有所需數(shù)據(jù)量小、意義明確等優(yōu)點, 常用于水災(zāi)和旱災(zāi)損失的風險評估[21–23]。歐陽蔚等研究表明傳統(tǒng)信息擴散估計缺少精確度, 而基于Bootstrap重抽樣技術(shù)的復(fù)合信息擴散得到的區(qū)間估計結(jié)果更精確可信[24]。其基本原理如下:

(1) 構(gòu)造重抽樣本。設(shè)原始樣本X=(x1,x2,…,xm),采用隨機抽樣法從樣本X中抽取N組容量為m的樣本X'(t)=(x1',x2',…,xm'), t=1,2,…,N, N是整數(shù), 則X'(t)稱為X的Bootstrap樣本。

(2)重抽樣信息擴散。設(shè)樣本信息擴散域為U=(u1,u2,…,un), n表示擴散控制點數(shù)量, 一般小于樣本容量 m。將樣本X'(t)中各樣本點帶入公式 3中,可將所帶信息擴散到U中。

公式3中, i=1,2,…,m; j=1,2,…,n; h表示擴散系數(shù),

式中, m表示樣本個數(shù), a=max{X'(t )}, b=min{X'(t )}

(3)將N組bootstrap樣本帶入步驟2可得N組超越概率 Pt'。將uj處對應(yīng)的超越概率值Pt'(u≥uj)共 N組, 按大小重新排列, 則該點置信區(qū)間為{Pt1'(u≥uj),Pt2'(u≥uj)}, t1=round(N·(α/2)), t2=round(N·(1-α/2)), 表示顯著性水平。

3 健康經(jīng)濟損失風險評估

3.1 數(shù)據(jù)來源

表1整理了2006—2015年上海市PM2.5年平均濃度、常住人口、GDP和居民消費水平。其中, 2006—2012 PM2.5年平均濃度, 由每年的上海環(huán)境狀況公告中PM10濃度, 按世界標準0.6進行轉(zhuǎn)換得到[10],2013—2015 PM2.5年濃度數(shù)據(jù)來源于上海市環(huán)境保護局, 包含了2013年1月1日至2015年7月1日的PM2.5日均濃度值。PM2.5濃度閾值同時選用最新國家空氣質(zhì)量二級標準[10], 即 35 μg·m-3和 WHO(世界衛(wèi)生組織)空氣質(zhì)量準則值[25], 即 10 μg·m-3, 避免單一標準誤差過大。暴露人口即上海市常住人口,居民消費水平和GDP來自2007—2014年上海市統(tǒng)計年鑒, 并根據(jù)平均增長率推算了2014、2015年常住人口。

從表1看出, 2006年至2015年上海市PM2.5平均濃度最高出現(xiàn)在2015年, 達到55.7 μg·m-3, 高出國家二級標準59.14%, 最低點為2012年, 達到42.6 μg·m-3,高出了國家二級標準 21.71%, 說明上海市 PM2.5污染狀況很嚴重。這10年間上海市經(jīng)濟社會快速發(fā)展,2013年上海常住人口比5年前增加了12.82%, 2013年上海GDP比5年前增加53.53%, 居民消費水平比5年前增加46.59%。

3.2 健康損失及超越概率

選擇霧霾污染的主要因素 PM2.5作為健康損失值核算的污染因子。鑒于我國南方北方間的地理氣候、能源結(jié)構(gòu)和收入水平差別很大[26–27], 不同省市的人口對霧霾污染的敏感反應(yīng)程度不同且不同收入水平下相同健康問題的費用差異較大。為了準確計算上海市行政區(qū)域內(nèi)的霧霾健康損失, 綜合考慮研究對象和參數(shù)的完整統(tǒng)一性, 文章選用上海交通大學(xué)趙文昌博士關(guān)于上海市 PM2.5的健康效應(yīng)終點以及對應(yīng)的暴露關(guān)系系數(shù)、健康終點發(fā)病率基數(shù)以及相應(yīng)的經(jīng)濟損失值[28], 其中健康效應(yīng)終點包括: 呼吸系統(tǒng)發(fā)病率、呼吸系統(tǒng)住院率、哮喘患病率、心血管疾病死亡率及活動受限日, 以上效應(yīng)均屬于慢性效應(yīng)。霧霾對人體健康的危害可分為急性效應(yīng)和慢性效應(yīng), 流行病學(xué)研究表明霧霾顆粒物慢性效應(yīng)遠大于其急性效應(yīng)[8,12], 即上述慢性效應(yīng)能全面反映霧霾健康經(jīng)濟損失實際值。利用公式(1)和(2)可計算得到2006年到2015年上海市PM2.5引起的健康經(jīng)濟損失, 其中單位經(jīng)濟損失都以 2009年為基期, 按CPI變化做相應(yīng)調(diào)整, 兩種標準的計算結(jié)果如表2所示。

表1 2006—2015年上海市經(jīng)濟環(huán)境情況Tab. 1 Economic and environmental condition in Shanghai from 2006 to 2015

由表2可知, 2006—2015年期間, 國家標準下,上海市PM2.5最大健康損失在2013年, 為120.89億元, 最小損失值在2012年, 為23.62億元; 2015年占GDP比例最大, 為4.60‰; 2012年占GDP比例最小,為 1.57‰; 2006—2015年, WHO 標準下, 上海市PM2.5健康損失值最大為 218.23億元, 最小值為86.76億元; 占 GDP比例最大為 9.30‰, 最小為5.77‰。健康經(jīng)濟損失及占GDP比例上都表現(xiàn)出波動上升規(guī)律, 這與 PM2.5年平均濃度增加密不可分,另一方面, 上海市人口、居民消費水平和國內(nèi)生產(chǎn)總值的快速增長起到了放大作用。

采用復(fù)合信息擴散模型, 對上述經(jīng)濟損失計算超越概率及可信區(qū)間。其中, 經(jīng)濟損失擴散域設(shè)定為[20, 240], 風險控制點數(shù)量 m=23, 原始樣本為兩i種標準下的 20個經(jīng)濟損失值, 抽樣次數(shù) N=10000,顯著水平α取0.25和0.1。由表1可知, 兩種標準得到的經(jīng)濟損失差額達到 1倍以上, 這種極端的結(jié)果不能完全反映健康損失的實際分布, 為此將同年的兩種結(jié)果平均值也納入原始數(shù)據(jù)中, 用來計算上海市PM2.5健康損失的信息擴散超越概率以及Bootstrap抽樣的區(qū)間估計, 結(jié)果見圖1, 2和表3, 4。

表2 2006—2015年上海市PM2.5健康經(jīng)濟損失Tab. 2 Economic loss of PM2.5health damage in Shangha from 2006 to 2015

圖1 上海市PM2.5健康經(jīng)濟損失風險曲線Fig. 1 Risk curve of economic loss of PM2.5health damage in Shanghai

圖2 上海市PM2.5健康經(jīng)濟損失占GDP比例風險曲線Fig. 2 Risk curve of GDP proportion for economic loss of PM2.5health damage in Shanghai

表3 上海市PM2.5健康經(jīng)濟損失超越概率Tab. 3 Exceeding probability of economic loss of PM2.5health damage in Shanghai

表4 上海市PM2.5健康損失占GDP比例超越概率Tab. 4 Exceeding probability of GDP proportion for economic loss of PM2.5health damage in Shanghai

以信息擴散模型計算的超越概率為基礎(chǔ), 利用bootstrap重抽樣技術(shù)來確定超越概率的置信區(qū)間能提高分析的精確度, 避免單一概率值的誤差過大。由步驟(2)計算得到信息擴散超越概率, 通過步驟(3)可以得到不同置信水平下的置信區(qū)間, 結(jié)果如圖 1所示。圖 1中, “區(qū)間估計上下限”指各控制點處超越概率均值的置信區(qū)間[μ- σ, μ+ σ], “超越概率上下限1、2”表示兩種置信水平時超越概率的區(qū)間估計值。一般情況下置信水平越高, 結(jié)果越可信。但研究結(jié)果表明, 在區(qū)間估計時置信水平過高會導(dǎo)致區(qū)間估計過于寬泛。為此, 借鑒歐陽蔚等學(xué)者對置信水平的研究結(jié)論, 設(shè)定置信水平分別為 0.75、0.9。以經(jīng)濟損失100億元為例做進一步說明, 0.9置信水平下的區(qū)間估計值為[29.11%, 54.57%], 0.7置信水平下的區(qū)間估計值為[32.54%, 50.51%], 而信息擴散模型計算的超越概率為41.58%時, 對應(yīng)的均值區(qū)間估計為[33.96%, 49.34%]。由此可知, 置信水平為0.9時的估計區(qū)間包含了置信水平為0.75時的估計區(qū)間, 置信概率 0.75的區(qū)間估計比置信概率為 0.9的區(qū)間估計更接近信息擴散超越概率及其正態(tài)區(qū)間估計值, 可見置信概率0.75的區(qū)間估計更精確。

圖1顯示了上海市10年來PM2.5健康經(jīng)濟損失的逐漸增加, 對應(yīng)的超越概率逐次降低, 超越概率的置信區(qū)間先逐漸變大, 然后逐漸變小, 起點和終點都收斂于同一處。表3顯示了不同風險水平下的健康經(jīng)濟損失、對應(yīng)的超越概率以及對應(yīng)的置信區(qū)間。經(jīng)濟損失超過60億元時, 超越概率為73.52%,置信區(qū)間為[66.32%, 80.36%]; 當經(jīng)濟損失超過100億元時, 超越概率為 41.58%, 置信區(qū)間為[32.49%,50.28%]; 而經(jīng)濟損失超過160億元時, 超越概率為15.46%, 置信區(qū)間為[8.46%, 22.36%]; 最后經(jīng)濟損失超過200億元時, 超越概率為6.47%, 置信區(qū)間為[2.28%, 10.59%]。

如圖2所示, 置信概率0.75時的區(qū)間估計比置信概率為00.9時更精確。同時健康損失占GDP比例逐漸降低, 對應(yīng)的超越概率也逐漸降低, 對應(yīng)的超越概率估計區(qū)間先逐漸變大然后變小, 起點和終點都收斂于同一點。健康損失 GDP比例的風險水平、具體數(shù)值、超越概率以及區(qū)間估計, 如表4所示。健康損失占GDP比例超過3‰時的超越概率為88.98%, 置信區(qū)間為[84.28, 93.43]。健康損失占GDP比例超過 5‰的概率為 60.41%, 置信區(qū)間為[51.23, 69.23]; 健康損失占GDP比例超過8‰的概率為 15.41%, 置信區(qū)間為[9.03, 21.78]; 健康損失占 GDP比例超過 10‰的概率為 1.36%, 置信區(qū)間為[9.03, 21.78]。

3.3 數(shù)學(xué)期望和重現(xiàn)周期

依據(jù)災(zāi)害損失的風險期望計算方法[29], 可得2006-2015年上海市霧霾顆粒物PM2.5的健康損失值數(shù)學(xué)期望為94.64億元, 按照2013年上海常住人口計算, 人均健康損失約392元。而2006—22015年上海PM2.5的健康損失占GDP比例的數(shù)學(xué)期望為5.42‰。

為進一步分析上海市 PM2.5健康經(jīng)濟損失的風險狀況, 使用超越概率的倒數(shù)表示不同風險水平的重現(xiàn)周期[21,29], 即不同規(guī)模的損失多少年一遇。由上述超越概率求倒數(shù), 可得重現(xiàn)周期的不同健康損失規(guī)模分布, 結(jié)果如表5所示。

由表5和圖3可以看出, 隨著重現(xiàn)周期增加, 損失規(guī)模不斷上升, 但增加速度不斷減小, 逐漸趨向穩(wěn)定。具體地講, 當遇到5年一遇的霧霾污染時, 健康經(jīng)濟損失規(guī)模達到144億元左右, 按照2013年上海常住人口計算, 人均健康損失約 596元; 當遇到10年一遇的霧霾污染時, 健康經(jīng)濟損失規(guī)模達到183億元左右, 按 2013年上海常住人口計算, 人均健康損失約758元; 當遇到15年一遇的霧霾污染時,健康經(jīng)濟損失規(guī)模達到199億元規(guī)模, 按照2013年上海常住人口計算, 人均健康損失達824元。

3.4 健康損失預(yù)測

健康經(jīng)濟損失的變化總體呈現(xiàn)出快速增加的趨勢。為了準確預(yù)測未來的損失情況, 采用數(shù)據(jù)擬合對未來幾年的情況進行預(yù)測, 歷史數(shù)據(jù)使用表 2中的損失值。考慮到數(shù)據(jù)量大小, 分別使用線性、二次多項式、三次多項式擬合健康損失序列[30]。在擬合前, 需要對數(shù)據(jù)預(yù)處理。常用的0—1標準化處理能解決量綱不統(tǒng)一的問題, 但嚴重壓縮了數(shù)據(jù)中的有用信息, 為此將2006年設(shè)定為橫坐標基期, 依次處理年份。擬合參數(shù)如表6、7所示。

表5 PM2.5健康損失重現(xiàn)周期Tab. 5 Recurrence period of economic loss of PM2.5health damage in Shanghai

表6 國家標準霧霾損失擬合效果Tab. 6 Polynomial fitting of economic loss of haze health in Shanghai under national standard

表7 WHO標準霧霾損失擬合效果Tab. 7 Polynomial fitting of economic loss of PM2.5health damage under WHO standard

擬合參數(shù)中, R2(確定系數(shù))和Adjusted-R2越接近 1, SSE(和方差)和 RMSE(均方根)越接近0, 擬合效果越好。由表6、7中可以看出二次多項式回歸的擬合效果最理想。圖3、4顯示了二次多項式擬合及預(yù)測的具體情況。二次多項式擬合結(jié)果如圖3、4所示。設(shè)定預(yù)測區(qū)間為[2016,2020], 根據(jù)圖3可知國家標準下, 5年后的健康經(jīng)濟損失(億元)分別為140.02,166.17, 195.19, 227.08, 261.85; WHO標準下, 未來5年的健康經(jīng)濟損失(億元)分別為 249.22, 283.49,320.82, 361.22, 404.68; 2020年的損失額度將是2015年的3.5倍。結(jié)果表明如果沒有過多干預(yù), 霧霾產(chǎn)生的健康損失將會出現(xiàn)迅速增加。目前需要對于霧霾的治理進度亟需加快, 對于未來可能出現(xiàn)的大規(guī)模公共衛(wèi)生問題, 也需要相關(guān)部門做好預(yù)防準備工作。

圖3 國家標準下上海霧霾健康損失預(yù)測Fig. 3 Forecast of economic loss of PM2.5health damage in Shanghai under national standard

圖4 WHO標準下上海霧霾健康損失預(yù)測Fig. 4 Forecast of economic loss of PM2.5health damage in Shanghai under WHO standard

4 結(jié)論與討論

4.1 結(jié)論

文章通過暴露反應(yīng)關(guān)系計算了國家標準和 WHO標準下上海市2006年至2015年P(guān)M2.5的健康經(jīng)濟損失及占同年GDP的比例。兩種標準條件下健康損失值都表現(xiàn)出波動增長趨勢, 健康損失占 GDP比例也表現(xiàn)出波動上漲趨勢, 但健康損失增長速度明顯高于占 GDP比例的增加速度。健康損失值的波動幅度越來越大, 而健康損失占 GDP比例變化幅度逐漸減小。

在健康經(jīng)濟損失數(shù)據(jù)基礎(chǔ)上, 聯(lián)合Bootstrap重抽樣技術(shù)和信息擴散模型對經(jīng)濟健康損失仿真實驗,得到霧霾健康經(jīng)濟損失值、超越概率和置信區(qū)間。并且根據(jù)超越概率計算得到不同損失額度的數(shù)學(xué)期望和重現(xiàn)周期。其中, 上海市PM2.5的健康經(jīng)濟損失超過80億元的概率為56.62%, 置信區(qū)間為[47.96%,65.34%]; 超過100億元的發(fā)生概率為41.58%, 置信區(qū)間為[32.94%, 50.51%]; 2006年至 2015年期間,PM2.5健康經(jīng)濟損失平均值為94.64億元, 占GDP比例平均水平為5.42‰; 上海PM2.5健康經(jīng)濟損失5年內(nèi)會出現(xiàn)144億元的規(guī)模, 10年內(nèi)會出現(xiàn)183億元的水平。

最后根據(jù)健康損失值2006年至2015年的數(shù)據(jù),利用多項式擬合得到了健康損失值隨年份變化的函數(shù)關(guān)系, 然后依據(jù)擬合函數(shù)預(yù)測了未來 5年上海市PM2.5的健康經(jīng)濟損失值。二次多項式擬合函數(shù)的預(yù)測結(jié)果顯示: 未來 5年上海的霧霾污染對人體健康危害會進一步擴大, 2020年國家標準下上海市PM2.5健康經(jīng)濟損失將達到261.85億元。

4.2 討論

在計算健康經(jīng)濟損失時, 污染物時間空間分布差異、人口年齡結(jié)構(gòu)和暴露人口空間分布差異都會使經(jīng)濟損失值結(jié)果不同。計算模型中污染物的暴露反應(yīng)系數(shù)、閾值有多種標準, 不同標準下的計算結(jié)果也不同。考慮到上述問題, 運用信息擴散模型計算了不同風險水平的超越概率值并通過重抽樣技術(shù)確定了置信區(qū)間, 同時完成情景分析和模擬預(yù)測,實現(xiàn)了霧霾健康損失的風險評估。霧霾健康損失有一定的時滯性, 而且對經(jīng)濟社會產(chǎn)生許多間接損失,這些都需要做進一步的研究。

參考文獻

[1] 周馮琦, 湯慶合, 任文偉. 上海資源環(huán)境發(fā)展報告(2014)[M]. 北京: 社會科學(xué)文獻出版社, 2014.

[2] 郭新彪, 魏紅英. 大氣 PM2.5對健康影響的研究進展[J].科學(xué)通報, 2013, 58(13) ): 1171–1177.

[3] 闞海東, 陳仁杰. PM2.5對人體危害有多大[J]. 中國經(jīng)濟報告, 2015(4): 114–116.

[4] 穆泉, 張世秋. 2013年1月中國大面積霧霾事件直接社會經(jīng)濟損失評估[J]. 中國環(huán)境科學(xué), 2013, 33(11):2087–2094.

[5] QUAH E, BOON T L. The economic cost of particulate air pollution on health in Singapore[J]. Journal of Asian Economics, 2002, 14(1): 73–90.

[6] MATUS K, NAM K, SELIN N E., et al. Health damages from air pollution in China[J]. Global Environmental Change, 2011, 22(1): 55–66.

[7] 陳仁杰, 陳秉衡, 闞海東. 我國113個城市大氣顆粒物污染的健康經(jīng)濟學(xué)評價[J]. 中國環(huán)境科學(xué), 2010, 30(3):410–415.

[8] 劉曉云, 謝鵬, 劉兆榮, 等. 珠江三角洲可吸入顆粒物污染急性健康效應(yīng)的經(jīng)濟損失評價[J]. 北京大學(xué)學(xué)報(自然科學(xué)版), 2010, 46(5): 829–834.

[9] 黃德生, 張世秋. 京津冀地區(qū)控制PM2.5污染的健康效益評估[J]. 中國環(huán)境科學(xué), 2013, 33(1): 166–174.

[10] 潘小川, 李國星, 高婷. 危險的呼吸––PM2.5的健康危害和經(jīng)濟損失評估研究[M]. 北京: 中國環(huán)境科學(xué)出版社,2012.

[11] 謝元博, 陳娟, 李巍. 霧霾重污染期間北京居民對高濃度 PM2.5持續(xù)暴露的健康風險及其損害價值評估[J]. 環(huán)境科學(xué), 2014, 35(1): 1–8.

[12] 陳依, 柏益堯, 錢新. 南京市 PM10人群健康經(jīng)濟損失評估[J]. 安徽農(nóng)業(yè)科學(xué), 2015, 43(19): 248–250+303.

[13] 侯青, 安興琴, 王自發(fā), 等．2002–2009年蘭州 PM10人體健康經(jīng)濟損失評估[J]. 中國環(huán)境科學(xué), 2011, 31(8):1398–1402.

[14] KAN Haidong, CHEN Bingheng. Particulate air pollution in urban areas of Shanghai, China: health-based economic assessment[J]. Science of the Total Environment, 2004, 322:71–79.

[15] 姜綿峰, 葉春明. 上海城市生態(tài)足跡動態(tài)研究––基于ARIMA模型[J]. 華東經(jīng)濟管理, 2015, 29(1): 18–24.

[16] 姜綿峰, 葉春明. 上海市建設(shè)用地時空演變的驅(qū)動力研究[J]. 資源開發(fā)與市場, 2015, 31(2): 160–165.

[17] 毛熙彥, 蒙吉軍, 康玉芳. 信息擴散模型在自然災(zāi)害綜合風險評估中的應(yīng)用與擴展[J]. 北京大學(xué)學(xué)報(自然科學(xué)版), 2012, 48(3): 513–518.

[18] 過孝民, 於方, 趙越, 環(huán)境污染成本評估理論與方法[M].北京: 中國環(huán)境科學(xué)出版社, 2009.

[19] 曾先峰, 王天瓊, 李印. 基于損害的西安市大氣污染經(jīng)濟損失研究[J]. 干旱區(qū)資源與環(huán)境, 2015, 29(1):105–110.

[20] WEI Jiuchang, GUO Xiumei, MARINOVA D, et al.Industrial SO2pollution and agricultural losses in China:evidence from heavy air polluters[J]. Journal of Cleaner Production, 2014, 64: 404–413.

[21] 歐朝敏, 尹輝, 張磊. 洞庭湖區(qū)不同情景下農(nóng)業(yè)水旱災(zāi)害風險損失評估[J]. 農(nóng)業(yè)現(xiàn)代化研究, 2011, 32(6):691–694.

[22] 張竟竟. 河南省農(nóng)業(yè)水旱災(zāi)害風險評估與時空分布特征[J]. 農(nóng)業(yè)工程學(xué)報, 2012, 28(18): 98–106.

[23] 丁青云, 艾萍, 吳軍斕, 等. 基于信息擴散理論的干旱災(zāi)害風險評估[J]. 中國農(nóng)村水利水電, 2015(3): 99–102.

[24] 歐陽蔚, 于艷青, 金菊良, 等. 基于信息擴散與自助法的旱災(zāi)風險評估模型––以安徽為例[J]. 災(zāi)害學(xué), 2015, 30(1):228–234.

[25] World Health Organization. Air quality guidelines for particulate matter, ozone, nitrogen dioxide and sulfur dioxide, Global update 2005, Summary of risk assessment[R]: WHO Press, 2006.

[26] CHENG Zhen, JIANG Jingkun, FAJARDO O, et al.Characteristics and health impacts of particulate matter pollution in China (2001-2011)[J]. Atmospheric Environment, 2013, 65: 186–194.

[27] AN Xingqin, HOU Qing, LI Nan, et al. Assessment of human exposure level to PM10 in China[J]. Atmospheric Environment, 2013, 70: 376–386.

[28] 趙文昌. 空氣污染對城市居民的健康風險與經(jīng)濟損失的研究[D]. 上海: 上海交通大學(xué), 2012.

[29] 杜鵑, 汪明, 史培軍. 基于歷史事件的暴雨洪澇災(zāi)害損失概率風險評估––以湖南省為例[J]. 應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報, 2014, 22(5): 916–927.

[30] 付艷茹. 基于MATLAB曲線擬合的應(yīng)用研究[J]. 吉林師范大學(xué)學(xué)報(自然科學(xué)版), 2010(2): 55–58.

Risk assessment for economic loss of haze health damage in Shanghai

JIANG Mianfeng1, YE Chunming1,*, SHENG Zhenzhen1, SHENG Xiaoxing1, HUA Yu2
1. Business School, University of Shanghai for Science and Technology, Shanghai 200090, China 2. Graduate College, Guangdong Medical University, Zhanjiang, 524000, Guangdong, China

For quantitative evaluation of health damage towards local people of haze particles PM2.5, exposure response unction was used to calculate the economic loss of haze health damage in Shanghai and its proportion of GDP based on available data of Shanghai in 2006-2015. The paper analyzed exceeding probability, expectation and recurrence period of different risk levels by using the bootstrap information diffusion comprehensive model, and predicted the economic loss and its hare of GDP in the next five years by applying polynomial fitting. Results indicated that the highest loss was 120.89 billion yuan in 2013 under the national standard, and the highest proportion of GDP in 2015 was 4.60‰, and the probability of exceeding 8 billion yuan was 56.62%. Average economic loss reached 94.64 billion yuan in Shanghai from 2006 to 2015, and he average proportion of GDP was 5.42‰. It will happen once or more that the economic loss exceeds 144 billion yuan in next ive years. The economic loss in Shanghai is going to exceed 261.85 billion yuan in 2020. Considering that economic loss of haze health damage in Shanghai had risen rapidly, effective measures should be adopted to control pollution, and the warning mechanism is also necessary in case of potential emergence.

Shanghai; haze; economic loss; proportion of GDP; exceeding probability

10.14108/j.cnki.1008-8873.2017.03.013

F205; X821

A

1008-8873(2017)03-090-08

姜綿峰, 葉春明, 盛真真, 等. 上海市霧霾健康經(jīng)濟損失風險評估[J]. 生態(tài)科學(xué), 2017, 36(3): 90-97.

JIANG Mianfeng, YE Chunming, SHENG Zhenzhen, et al. Risk assessment of economic loss for haze health damage in Shanghai[J].Ecological Science, 2017, 36(3): 90-97.

2015-11-14;

2015-12-11

國家自然科學(xué)基金項目(71271138); 上海市哲學(xué)社會科學(xué)規(guī)劃項目(2014BGL024); 上海市一流學(xué)科資助項目(S1201YLXK); 滬江基金資助項目(A14006)

姜綿峰(1991—), 男, 河南商城人, 碩士研究生, 主要從事生態(tài)經(jīng)濟及風險管理研究, E-mail: jiangmianfeng@yeah.net

*通信作者:葉春明(1964—), 男, 安徽宣城人, 博士, 教授, 主要從事技術(shù)經(jīng)濟及系統(tǒng)工程研究, E-mail: yechm6464@163.com