河北省霧霾污染的空間分布及其影響因素的實(shí)證分析

回 瑩，毛 培，戴宏偉

（中央財(cái)經(jīng)大學(xué) 經(jīng)濟(jì)學(xué)院，北京 100081）

近幾年來，中國霧霾天氣頻發(fā)，京津冀、長(zhǎng)三角等多地成為中國主要的霧霾污染重災(zāi)區(qū)，邢臺(tái)、邯鄲、唐山等市更是長(zhǎng)期出現(xiàn)在城市空氣質(zhì)量較差的前十城市排行榜上。霧霾天氣頻發(fā)，一部分是由于氣象因素，更深層次的原因是產(chǎn)業(yè)結(jié)構(gòu)不合理、城鎮(zhèn)化水平粗放增長(zhǎng)、汽車尾氣排放等因素密切相關(guān)。本文以中國霧霾污染最嚴(yán)重的河北省為例，探討霧霾污染的空間溢出效應(yīng)及其與產(chǎn)業(yè)結(jié)構(gòu)、城鎮(zhèn)化水平之間的相互關(guān)系。

2017年4月1日，河北雄安新區(qū)設(shè)立，未來河北雄安新區(qū)要建成一座“綠色生態(tài)宜居新城區(qū)、創(chuàng)新驅(qū)動(dòng)發(fā)展引領(lǐng)區(qū)、協(xié)調(diào)發(fā)展示范區(qū)、開放發(fā)展先行區(qū)”，而雄安新區(qū)位于河北省霧霾污染重地，未來河北省治霾形勢(shì)十分嚴(yán)峻。目前河北省已經(jīng)出臺(tái)一系列的政策措施來緩解霧霾污染，京津冀地區(qū)也采取了一系列的治霾聯(lián)防聯(lián)控措施，《大氣污染防治行動(dòng)計(jì)劃》中明確提出，2017年河北省PM2.5濃度值下降25%。但由于霧霾成因的復(fù)雜性，河北省霧霾治理形勢(shì)依然嚴(yán)峻。因此，研究河北省霧霾污染的成因，河北省霧霾污染的影響因素具有重要的現(xiàn)實(shí)意義。

一、文獻(xiàn)綜述

隨著經(jīng)濟(jì)的發(fā)展，環(huán)境問題越來越突出。越來越多的學(xué)者開始研究環(huán)境與經(jīng)濟(jì)發(fā)展及其他方面的關(guān)系。Panayoutou（1993）提出了著名的“環(huán)境庫茨涅茨曲線”（EKC）[1]，隨后Selden et al.（1995）、Brajer et al.（2011）等人從經(jīng)濟(jì)的不同角度研究了經(jīng)濟(jì)發(fā)展與環(huán)境之間的關(guān)系，認(rèn)為二者或存在“倒U型曲線”或“J型曲線”[2-3]。我國一些學(xué)者也從不同角度論證了經(jīng)濟(jì)發(fā)展與環(huán)境之間的關(guān)系：包群 等（2006）利用中國省級(jí)數(shù)據(jù)驗(yàn)證了“環(huán)境庫茨涅茨曲線”在中國的正確性[4]；張少華（2009）基于行業(yè)面板數(shù)據(jù)研究發(fā)現(xiàn)經(jīng)濟(jì)全球化顯著降低了我國的環(huán)境污染[5]；張悅 等（2014）則使用中國經(jīng)濟(jì)增長(zhǎng)和環(huán)境污染等相關(guān)指標(biāo)對(duì)“環(huán)境庫茨涅茨曲線”進(jìn)行驗(yàn)證后，發(fā)現(xiàn)我國經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染之間呈現(xiàn)出倒N型關(guān)系[6]。

伴隨著霧霾天氣的頻發(fā)，越來越多的學(xué)者開始研究霧霾污染與經(jīng)濟(jì)發(fā)展之間的關(guān)系：郭俊華 等（2014），冷艷麗 等（2015），李鵬（2015）等人論證了中國產(chǎn)業(yè)結(jié)構(gòu)與霧霾污染之間的關(guān)系，均得出中國霧霾污染與產(chǎn)業(yè)結(jié)構(gòu)具有很強(qiáng)的相關(guān)關(guān)系[7-9]；馬麗梅 等（2014）基于中國31個(gè)省份的數(shù)據(jù)研究認(rèn)為霧霾污染與產(chǎn)業(yè)結(jié)構(gòu)、能源結(jié)構(gòu)密切相關(guān)，人均GDP的增長(zhǎng)導(dǎo)致霧霾污染的加劇[10]；胡飛（2011）實(shí)證分析了中國東部和中部地區(qū)的產(chǎn)業(yè)結(jié)構(gòu)升級(jí)和霧霾污染之間的關(guān)系[11]。王家庭 等（2010）、劉伯龍 等（2015）[12-13]則探討了中國城市化與環(huán)境污染之間的反“U型”和正相關(guān)關(guān)系。

京津冀地區(qū)作為中國霧霾污染的重災(zāi)區(qū)，也成為學(xué)者們的研究對(duì)象之一：潘慧峰 等（2015）對(duì)京津冀地區(qū)7個(gè)城市的PM2.5進(jìn)行估計(jì)后，發(fā)現(xiàn)京津冀地區(qū)的霧霾污染存在溢出效應(yīng)[14－15]；馮博 等（2015）則認(rèn)為京津冀地區(qū)考慮霧霾效應(yīng)之后能源效率有所下降[16]；王穎 等（2016）運(yùn)用跨界治理理論破解京津冀地區(qū)治理困境[17]；李瑞 等（2014）分析了河北省霧霾污染嚴(yán)重的成因并提出相關(guān)政策建議[18]。

空間計(jì)量學(xué)的發(fā)展為區(qū)域經(jīng)濟(jì)增長(zhǎng)問題又提供了一種新的研究方法，冷艷麗 等（2015）、馬麗梅 等（2014）、邵帥 等（2016）學(xué)者均使用中國省級(jí)數(shù)據(jù)對(duì)中國霧霾污染進(jìn)行了空間分析[8，10，19]；胡宗義 等（2015）使用面板數(shù)據(jù)建立空間模型分析了能源消耗、污染與區(qū)域經(jīng)濟(jì)增長(zhǎng)之間的動(dòng)態(tài)關(guān)系[20]。

從以上文獻(xiàn)可以看出，大量學(xué)者開始重視環(huán)境問題，研究環(huán)境污染與經(jīng)濟(jì)增長(zhǎng)之間關(guān)系的文獻(xiàn)較多；近年來，霧霾污染與產(chǎn)業(yè)結(jié)構(gòu)、霧霾污染與城鎮(zhèn)化等不同影響因素之間關(guān)系的文獻(xiàn)增多，但鮮有使用空間計(jì)量方法研究霧霾污染及其影響因素的文獻(xiàn)，研究河北省霧霾污染的空間特征及其影響因素的實(shí)證分析更加少有。本文對(duì)河北省2000—2012年的PM2.5濃度值進(jìn)行了空間相關(guān)性分析，并進(jìn)一步構(gòu)建空間杜賓模型進(jìn)行實(shí)證分析得出相關(guān)結(jié)論。

二、霧霾污染的空間相關(guān)性分析

（一）空間權(quán)重矩陣的選擇

1.空間地理權(quán)重（W1）。該權(quán)重基于車相鄰（rook contiguity）方式設(shè)定河北省11個(gè)地級(jí)市的二進(jìn)制鄰接空間權(quán)重矩陣，以反映河北省各地級(jí)市之間的空間關(guān)系，具體設(shè)定為：

2.地理－經(jīng)濟(jì)空間權(quán)重矩陣（W2）。為了分析更具穩(wěn)健性，在空間權(quán)重矩陣的基礎(chǔ)上增加經(jīng)濟(jì)因素，此種空間權(quán)重矩陣的設(shè)定方法是將地理因素和經(jīng)濟(jì)因素結(jié)合在一起（林光平 等，2005）[21]，具體設(shè)定如下：

式（2）中，W1為空間地理權(quán)重為 i地區(qū)的人均GDP的平均值為j地區(qū)的人均GDP的平均值。

（二）莫蘭指數(shù)分析

1.PM2.5數(shù)據(jù)來源。河北省各地級(jí)市對(duì)PM2.5數(shù)據(jù)的監(jiān)測(cè)始于2013年，為了得到更早的PM2.5的相關(guān)數(shù)據(jù)，本文借鑒馬麗梅 等（2014）[10]、邵帥 等（2016）[19]的文章，采用哥倫比亞大學(xué)社會(huì)經(jīng)濟(jì)數(shù)據(jù)和應(yīng)用中心（CIESIN）公布的關(guān)于PM2.5濃度值的全球年均值的遙感地圖①，并利用Arcgis等相關(guān)軟件提取出河北省11個(gè)地級(jí)市的PM2.5濃度值數(shù)據(jù)。

2.河北省PM2.5的Moran’s I。本文使用全局Moran’s I和局部Moran’s I來分析河北省霧霾污染的空間相關(guān)性，并利用GeoDa軟件得出結(jié)果。全局Moran’s I的計(jì)算公式為：

式（3）中，I表示全局 Moran’s I，S2表示樣本方差，Wij表示空間權(quán)重矩陣，Xi和Xj分別表示河北省城市i和城市j的霧霾污染程度，n表示河北省地級(jí)市數(shù)量，I的取值范圍為［－1，1］。當(dāng) I＞0 時(shí)，表示存在空間正相關(guān)，當(dāng)I＜0時(shí)，表示存在空間負(fù)相關(guān)，當(dāng)I接近0，表示不存在空間自相關(guān)。

利用GeoDa軟件，可以得出2000—2012年河北省 PM2.5濃度值的 Moran’s I，并繪制出 Moran’s I根據(jù)時(shí)間變化的趨勢(shì)圖（見圖1）。由圖1可以看出，河北省PM2.5濃度值的Moran’s I均大于0小于1，顯著為正，具有明顯的正自相關(guān)。這說明河北省的霧霾污染在空間上呈現(xiàn)出正相關(guān)性，即PM2.5濃度高的與濃度高的地級(jí)市相鄰，PM2.5濃度值低的與濃度值低的地級(jí)市相鄰，這也表明河北省各地級(jí)市的霧霾污染具有明顯的空間相互作用。

圖1 2000—2012年河北省PM2.5濃度值的Moran's I

局部Moran’s I的計(jì)算公式為：

式（4）中各指標(biāo)含義與式（1）相同。局部Moran’s I用來考察各地區(qū)被周圍地區(qū)包圍的情況。圖2－a和圖2－b以2006年和2012年為例，描繪了河北省PM2.5濃度值的局部Lisa圖。從圖2－a和圖2－b可以看出，河北省各地級(jí)市的PM2.5高濃度值的被高濃度值的包圍，低濃度值的被低濃度值的包圍，即呈現(xiàn)出高濃度值與高濃度值地區(qū)集聚、低濃度值與低濃度值集聚的空間分布特征。其中，邢臺(tái)市、邯鄲市為霧霾污染高-高型集聚區(qū)，承德市、秦皇島市為霧霾污染的低-低型集聚區(qū)。

圖2 2006年、2012年河北省PM2.5濃度值局部集聚（Lisa）

三、實(shí)證分析

（一）變量選取及數(shù)據(jù)來源

1.產(chǎn)業(yè)結(jié)構(gòu)。霧霾污染與產(chǎn)業(yè)結(jié)構(gòu)密切相關(guān)，尤其是河北省第二產(chǎn)業(yè)比重偏高，并多以煤炭開采和洗選業(yè)、黑色金屬礦采選業(yè)、黑色金屬冶煉和壓延加工業(yè)、金屬制品業(yè)等“高污染、高能耗、高成本”產(chǎn)業(yè)為主，這些“三高”產(chǎn)業(yè)產(chǎn)生大量的煙粉塵、含氮有機(jī)顆粒物等，這也是造成河北省霧霾污染嚴(yán)重的主要原因之一。本文選取河北省各地級(jí)市第二產(chǎn)業(yè)的區(qū)位商值代表河北省產(chǎn)業(yè)結(jié)構(gòu)，區(qū)位商值能更好的反映出第二產(chǎn)業(yè)在河北省所占的主導(dǎo)地位。相關(guān)數(shù)據(jù)主要來源于《中國區(qū)域經(jīng)濟(jì)統(tǒng)計(jì)年鑒》《中國城市統(tǒng)計(jì)年鑒》《中國統(tǒng)計(jì)年鑒》。

2.城鎮(zhèn)化水平。城鎮(zhèn)化水平也會(huì)對(duì)霧霾污染產(chǎn)生一定影響，隨著城市人口的增加，城市生產(chǎn)、生活產(chǎn)生的“三廢”等污染物增加，勢(shì)必會(huì)加重城市的霧霾污染。河北省城鎮(zhèn)化水平近幾年呈不斷上升趨勢(shì)，但呈粗放型增長(zhǎng)，城鎮(zhèn)化水平質(zhì)量不高。由于河北省各地級(jí)市數(shù)據(jù)所限，本文采取各地級(jí)市非農(nóng)人口占總?cè)丝诒戎貋肀硎境擎?zhèn)化水平。河北省各地級(jí)市人口數(shù)據(jù)主要來源于《全國分縣市人口統(tǒng)計(jì)資料》。

3．控制變量。除了產(chǎn)業(yè)結(jié)構(gòu)、城鎮(zhèn)化水平之外，還有大量因素會(huì)對(duì)霧霾污染產(chǎn)生影響。（1）實(shí)際人均GDP。Panayoutou（1993）[1]認(rèn)為環(huán)境質(zhì)量與人均收入之間存在著“倒U型”關(guān)系，這就是著名的“環(huán)境庫茨涅茨曲線”；Krugman（1990）也論證了人均收入與環(huán)境質(zhì)量的關(guān)系[22]；國內(nèi)學(xué)者張悅 等（2014）[6]也對(duì)經(jīng)濟(jì)增長(zhǎng)與環(huán)境之間的關(guān)系進(jìn)行了論證。本文借鑒相關(guān)文獻(xiàn)，選取實(shí)際人均GDP來表示河北省各地級(jí)市的經(jīng)濟(jì)增長(zhǎng)，該指標(biāo)使用各地級(jí)市的國民生產(chǎn)總值指數(shù)和GDP平減指數(shù)計(jì)算得來。數(shù)據(jù)主要來源于《中國城市統(tǒng)計(jì)年鑒》和《中國統(tǒng)計(jì)年鑒》。（2）外商直接投資額?！拔廴颈茈y所假說”（Copeland et al.，1994）[23]認(rèn)為環(huán)境管制嚴(yán)厲的發(fā)達(dá)國家會(huì)把產(chǎn)業(yè)轉(zhuǎn)移到環(huán)境標(biāo)準(zhǔn)低的發(fā)展中國家，從而把污染轉(zhuǎn)移到這些國家進(jìn)行投資和生產(chǎn)，使其成為發(fā)達(dá)國家的“污染避難所”；List et al.（2000）[24]、沙文兵等（2006）[25]等都論證了該假說的成立。冷艷麗 等（2015）[26]認(rèn)為外商直接投資與霧霾污染呈現(xiàn)出正相關(guān)關(guān)系。本文選取河北省11個(gè)地級(jí)市的外商直接投資額數(shù)據(jù)，并使用固定資產(chǎn)價(jià)格指數(shù)進(jìn)行平減，得到實(shí)際的外商直接投資額。該指標(biāo)數(shù)據(jù)主要來源于《中國城市統(tǒng)計(jì)年鑒》《中國區(qū)域經(jīng)濟(jì)統(tǒng)計(jì)年鑒》《新河北60年》等。（3）其他控制變量。為了增加統(tǒng)計(jì)結(jié)果的穩(wěn)健性，本文還選取了對(duì)外貿(mào)易依存度、人口密度、公路里程數(shù)等指標(biāo)。童玉芬 等（2014）[27]認(rèn)為人口規(guī)模與霧霾污染之間存在著雙向關(guān)系；秦蒙 等（2016）[28]認(rèn)為人口規(guī)模的增長(zhǎng)會(huì)提高PM2.5濃度，本文選取河北省各地級(jí)市的人口密度來表示人口規(guī)模。此外，公路交通機(jī)動(dòng)車尾氣等也是霧霾污染嚴(yán)重的主要原因之一，本文選取各地級(jí)市的公路里程數(shù)表示交通運(yùn)輸強(qiáng)度。以上指標(biāo)數(shù)據(jù)主要來源于《中國城市統(tǒng)計(jì)年鑒》《中國區(qū)域經(jīng)濟(jì)統(tǒng)計(jì)年鑒》《新河北60年》等。

（二）空間面板數(shù)據(jù)模型

本文選取2000—2012年河北省11個(gè)地級(jí)市空間面板數(shù)據(jù)來分析河北省霧霾污染的空間溢出效應(yīng)。使用（Robust）LM檢驗(yàn)對(duì)比結(jié)果如表1，發(fā)現(xiàn)空間地理權(quán)重矩陣和地理-經(jīng)濟(jì)空間權(quán)重矩陣下空間滯后模型均優(yōu)于空間誤差模型。

表1 LM檢驗(yàn)結(jié)果

本文選取更具廣泛意義的空間杜賓模型，模型設(shè)定如下：

（5）式中，PMit代表河北省i市在t時(shí)期的PM2.5的濃度值，str代表河北省產(chǎn)業(yè)結(jié)構(gòu)，urb代表城鎮(zhèn)化水平，Xit代表各控制變量，即實(shí)際人均GDP、外商直接投資額、對(duì)外貿(mào)易依存度、人口密度和公路里程數(shù)，W代表空間權(quán)重，α為常數(shù)項(xiàng)，β、γ、φ代表各參數(shù)，ρ為空間自回歸系數(shù)，θ、φ、η代表各變量空間滯后系數(shù)，u為服從正態(tài)分布的擾動(dòng)項(xiàng)。本文使用Matlab 2012b軟件進(jìn)行實(shí)證估計(jì)。

（三）實(shí)證結(jié)果與分析

依據(jù)式（3）按照空間計(jì)量方法進(jìn)行實(shí)證估計(jì)。根據(jù)Hausman檢驗(yàn)，兩種權(quán)重矩陣下的p值分別為0.004 2，0.000 0，均小于原假設(shè)條件0.025，拒絕如果p小于0.025接受隨機(jī)模型的原假設(shè)，因此選取固定效應(yīng)模型。

由表2可以看出，在空間地理權(quán)重下，混合模型和固定效應(yīng)模型中河北省的霧霾污染具有明顯的空間溢出效應(yīng)，修正后的R2在0.9以上，說明數(shù)據(jù)擬合的比較好。在W1的空間權(quán)重矩陣下，兩種模型中的產(chǎn)業(yè)結(jié)構(gòu)和城鎮(zhèn)化水平均在1%的置信水平上顯著；W2權(quán)重矩陣下，兩種模型的產(chǎn)業(yè)結(jié)構(gòu)在1%的置信水平上顯著，混合模型的城鎮(zhèn)化水平在5%的置信水平上顯著，固定效應(yīng)模型的城鎮(zhèn)化水平在1%的置信水平上顯著；兩種空間權(quán)重矩陣下變量系數(shù)均為正，說明產(chǎn)業(yè)結(jié)構(gòu)和城鎮(zhèn)化水平對(duì)河北省各市霧霾污染具有顯著正影響。隨著河北省第二產(chǎn)業(yè)比重的上升，主導(dǎo)地位的提升，城鎮(zhèn)化速度的上升，都會(huì)不利于河北省霧霾污染的緩解：混合模型下第二產(chǎn)業(yè)區(qū)位商值每提高1%，霧霾污染上升1%左右，固定效應(yīng)模型下第二產(chǎn)業(yè)區(qū)位商每提高1%，霧霾污染上升0.6%左右；城鎮(zhèn)化水平每提高1%，霧霾污染上升1.3%左右。實(shí)際人均GDP在1%的置信水平上顯著，系數(shù)為負(fù)，說明隨著經(jīng)濟(jì)的增長(zhǎng)，更加注重綠色增長(zhǎng)，非常有利于河北省霧霾污染的減輕。在兩種權(quán)重矩陣下，外商直接投資額分別在5%和1%的置信水平上顯著，系數(shù)均為負(fù)，說明河北省外商直接投資額的增加并沒有增加霧霾污染，這與“污染避難所”假說并不一致。對(duì)外貿(mào)易依存度對(duì)河北省霧霾污染具有正效應(yīng)，但影響并不顯著。人口密度和公路里程在兩種空間權(quán)重矩陣下在1%的水平上顯著，系數(shù)為正，說明人口的集聚和交通運(yùn)輸狀況對(duì)河北省霧霾污染具有顯著正效應(yīng)，不利于霧霾污染的緩解，人口密度和公路里程每提高1%，霧霾污染分別上升0.7%和0.1%左右。

霧霾污染還會(huì)受到上一期霧霾的影響，加入滯后變量對(duì)空間杜賓模型進(jìn)行估計(jì)后發(fā)現(xiàn)（見表2），兩種空間權(quán)重下，河北省霧霾污染仍然明顯的空間溢出效應(yīng)，尤其是加入經(jīng)濟(jì)因素后，參數(shù)ρ在1%的水平上顯著，空間溢出效應(yīng)明顯。在混合模型中，兩種空間權(quán)重矩陣的空間溢出系數(shù)為正，說明霧霾污染具有一定的連續(xù)性，上一期較高的PM2.5濃度值會(huì)導(dǎo)致當(dāng)期較高的PM2.5濃度值；考慮時(shí)間和空間因素后，即在時(shí)空固定效應(yīng)模型中，兩種空間權(quán)重矩陣的空間溢出系數(shù)為負(fù)，且分別在10%和1%的水平上顯著，說明上一期的霧霾污染對(duì)當(dāng)期具有較強(qiáng)的“警示效應(yīng)”（邵帥 等，2016）[19]，即如果河北省前一期的霧霾污染嚴(yán)重，會(huì)對(duì)當(dāng)期產(chǎn)生警示作用，從而采取一系列措施來減少霧霾污染，因此系數(shù)估計(jì)值為負(fù)?；旌夏Ｐ椭?，滯后一期的產(chǎn)業(yè)結(jié)構(gòu)和城鎮(zhèn)化水平對(duì)河北省霧霾污染具有減緩作用，但在固定效應(yīng)模型中不具有緩解作用。實(shí)際人均GDP和對(duì)外貿(mào)易依存度對(duì)河北省各市的霧霾污染仍具有一定的促增作用，對(duì)外直接投資額、人口密度、公路里程的滯后項(xiàng)系數(shù)估計(jì)值為負(fù)，說明會(huì)減緩霧霾污染。

表2 空間杜賓面板固定效應(yīng)模型估計(jì)結(jié)果

表3 空間杜賓模型中各解釋變量的空間效應(yīng)估計(jì)

下面進(jìn)一步分析各影響因素的直接效應(yīng)、間接效應(yīng)和總效應(yīng)。由于霧霾污染空間溢出效應(yīng)的存在，一個(gè)城市的霧霾污染的嚴(yán)重程度會(huì)受到該市產(chǎn)業(yè)結(jié)構(gòu)、城鎮(zhèn)化水平等影響的同時(shí)，該市的產(chǎn)業(yè)結(jié)構(gòu)、城鎮(zhèn)化水平等狀況也會(huì)影響其相鄰城市的PM2.5的濃度值。根據(jù)表3可以看出，在W1和W2兩種空間權(quán)重矩陣下，核心變量產(chǎn)業(yè)結(jié)構(gòu)和城鎮(zhèn)化水平的直接效應(yīng)的估計(jì)系數(shù)均顯著為正，說明第二產(chǎn)業(yè)區(qū)位商值和城鎮(zhèn)化水平的提高會(huì)加重該市PM2.5濃度值；間接效應(yīng)在W1權(quán)重矩陣下系數(shù)估計(jì)為負(fù)，在W2權(quán)重矩陣下系數(shù)估計(jì)為正，說明河北省各市第二產(chǎn)業(yè)區(qū)位商值和城鎮(zhèn)化水平的提高也會(huì)使在空間上相鄰、經(jīng)濟(jì)發(fā)展水平接近的地區(qū)霧霾污染有可能變得嚴(yán)重，從而使得河北省的霧霾污染呈現(xiàn)出高濃度值與高濃度值集聚和低濃度值與低濃度值集聚的狀況。總效應(yīng)上，產(chǎn)業(yè)結(jié)構(gòu)和城鎮(zhèn)化水平對(duì)河北省霧霾污染的影響顯著為正，在對(duì)自身影響的直接效應(yīng)和對(duì)鄰近地區(qū)的間接效應(yīng)雙重影響下，兩種空間權(quán)重矩陣下第二產(chǎn)業(yè)區(qū)位商值每提高1%，霧霾污染分別上升0.553、0.637，城鎮(zhèn)化水平每提高1%，霧霾污染分別上升0.855、1.548。由此可見，加入經(jīng)濟(jì)因素后，產(chǎn)業(yè)結(jié)構(gòu)和城鎮(zhèn)化水平對(duì)河北省霧霾污染的貢獻(xiàn)分別提高15.2%和81.1%，說明目前河北省產(chǎn)業(yè)結(jié)構(gòu)和城鎮(zhèn)化水平都處于粗放型增長(zhǎng)階段，未來還有很大的改進(jìn)空間和潛力。

實(shí)際人均GDP直接效應(yīng)、間接效應(yīng)的系數(shù)估計(jì)為負(fù)，說明隨著人均GDP的提高，經(jīng)濟(jì)的增長(zhǎng)，人們對(duì)環(huán)境質(zhì)量訴求提高，環(huán)境治理等方面投資增加，會(huì)使得本地區(qū)和相鄰地區(qū)的霧霾污染得到一定程度的緩解。從總效應(yīng)上看，僅考慮空間因素，河北省實(shí)際人均GDP的增加依然會(huì)使得PM2.5濃度值上升，加入經(jīng)濟(jì)因素后，PM2.5濃度值會(huì)有所下降，實(shí)際人均GDP每提高1%，PM2.5濃度值下降0.294個(gè)單位。

總體來看，外商直接投資并沒有加重河北省霧霾污染，對(duì)外貿(mào)易依存度對(duì)河北省霧霾污染具有正效應(yīng)，但影響均不顯著。人口密度和公路里程數(shù)對(duì)河北省霧霾污染均具有顯著的正效應(yīng)，說明河北省的人口集聚和交通運(yùn)輸狀況不僅不利于河北省PM2.5濃度值的降低，也不利于相鄰地區(qū)的PM2.5濃度值的降低。由此可見，人口密度和公路里程也是造成河北省霧霾污染較大部分地區(qū)嚴(yán)重的原因之一。

四、結(jié)論及政策建議

通過對(duì)河北省霧霾污染的空間分布及影響因素的實(shí)證研究，發(fā)現(xiàn)以下結(jié)論：

1.河北省各地級(jí)市的霧霾污染存在著明顯的空間正相關(guān)性，并且呈現(xiàn)出高濃度值與高濃度值集聚和低濃度值與低濃度值集聚的空間分布特征。

2.河北省霧霾污染具有明顯的空間溢出效應(yīng)；河北省霧霾污染既具有一定的連續(xù)性又具有一定的“警示效應(yīng)”。產(chǎn)業(yè)結(jié)構(gòu)和城鎮(zhèn)化水平對(duì)河北省霧霾污染具有正效應(yīng)，不僅不利于本地區(qū)霧霾污染的緩解，也不利于鄰近地區(qū)霧霾污染的緩解。

3.河北省霧霾污染并不符合“污染避難所”假說，說明河北省的外商直接投資沒有促進(jìn)霧霾污染的加重；對(duì)外貿(mào)易依存度、人口密度和公路里程數(shù)對(duì)霧霾污染均具有較強(qiáng)的正效應(yīng)。

根據(jù)上述研究結(jié)果，河北省的霧霾污染防治必須要從源頭上進(jìn)行控制，倫敦的“煙霧事件”后采取了一系列的治理措施，出臺(tái)了《清潔空氣法》《污染控制法》《環(huán)境法》等一系列法律，改善產(chǎn)業(yè)結(jié)構(gòu)，推廣清潔能源等，這些都為河北省的霧霾治理提供了借鑒經(jīng)驗(yàn)，但仍需進(jìn)一步細(xì)化標(biāo)準(zhǔn)，提高環(huán)境標(biāo)準(zhǔn)，配合稅收、確權(quán)等市場(chǎng)激勵(lì)措施，加大懲處力度，杜絕相關(guān)污染事件的發(fā)生。此外，河北省的霧霾污染治理必須與產(chǎn)業(yè)結(jié)構(gòu)調(diào)整、城鎮(zhèn)化進(jìn)程結(jié)合在一起，河北省的霧霾污染具有顯著的空間溢出效應(yīng)，必須制定區(qū)域性的聯(lián)防聯(lián)控措施，才能在“源頭”和“末端”科學(xué)合理的解決該地區(qū)的霧霾污染問題。目前河北省已投入大量資金，實(shí)施了一系列嚴(yán)格的霧霾污染防控防治措施，霧霾污染治理取得了很大成就，空氣質(zhì)量有所好轉(zhuǎn)，但仍需進(jìn)一步從產(chǎn)業(yè)結(jié)構(gòu)調(diào)整、城鎮(zhèn)化水平質(zhì)量提高等方面緩解霧霾污染。

1.進(jìn)一步調(diào)整產(chǎn)業(yè)結(jié)構(gòu)，區(qū)域聯(lián)防聯(lián)控需進(jìn)一步加強(qiáng)。重工業(yè)比重偏高，產(chǎn)業(yè)結(jié)構(gòu)不合理，是河北省霧霾污染嚴(yán)重的主要原因之一。因此，中央應(yīng)對(duì)河北省進(jìn)一步弱化GDP考核，強(qiáng)化綠色環(huán)保的考核，不斷降低河北第二產(chǎn)業(yè)比重，加快傳統(tǒng)產(chǎn)業(yè)的改造升級(jí)，改變能源結(jié)構(gòu)，大力推動(dòng)汽車制造業(yè)、交通運(yùn)輸設(shè)備、計(jì)算機(jī)、通信和其他電子設(shè)備制造業(yè)等新興產(chǎn)業(yè)的發(fā)展，依托京津，發(fā)展物流、會(huì)展經(jīng)濟(jì)、金融業(yè)、信息服務(wù)、旅游休閑等現(xiàn)代服務(wù)業(yè)，提高第三產(chǎn)業(yè)比重。同時(shí)，進(jìn)行工業(yè)產(chǎn)業(yè)結(jié)構(gòu)內(nèi)部調(diào)整，引進(jìn)新技術(shù)減少污染排放量，淘汰落后產(chǎn)能，節(jié)能減排，關(guān)停高污染高耗能產(chǎn)業(yè)，盡快完成脫硫脫銷除塵改造。加強(qiáng)與北京、天津的經(jīng)濟(jì)聯(lián)系強(qiáng)度，在進(jìn)行產(chǎn)業(yè)轉(zhuǎn)移和產(chǎn)業(yè)承接時(shí)，河北省各市與北京、天津形成梯度，主要承接現(xiàn)代制造業(yè)，避免霧霾污染從一個(gè)區(qū)域傳染到另一個(gè)區(qū)域，雄安新區(qū)重點(diǎn)發(fā)展高端高新產(chǎn)業(yè)和服務(wù)業(yè)，生態(tài)優(yōu)先，以雄安新區(qū)的設(shè)立帶動(dòng)河北省各市產(chǎn)業(yè)結(jié)構(gòu)的優(yōu)化升級(jí)，以區(qū)域產(chǎn)業(yè)結(jié)構(gòu)調(diào)整優(yōu)化、協(xié)同發(fā)展推動(dòng)河北省乃至京津冀霧霾的聯(lián)合治理。

2.進(jìn)一步采取措施提高城鎮(zhèn)化水平的質(zhì)量。改變以煤炭為主的生產(chǎn)、生活能源消費(fèi)結(jié)構(gòu)，大力推廣清潔能源，提高油品質(zhì)量，控制汽車尾氣排放，增加城市綠地面積。河北省應(yīng)完善城市功能，提高城鎮(zhèn)化增速，注重城鎮(zhèn)化水平的質(zhì)量，增強(qiáng)各市的吸引力，促進(jìn)城市要素的轉(zhuǎn)移，加快城市交通路網(wǎng)的建設(shè)，增強(qiáng)與京津的城市聯(lián)系度。

3.控制人口密度集中地區(qū)的人口，繼續(xù)采取汽車尾號(hào)限號(hào)等措施減少汽車尾氣的排放，同時(shí)加大公共交通設(shè)施的建設(shè)。加大環(huán)境治理投入，建立合理的區(qū)域補(bǔ)償機(jī)制，在引進(jìn)外資和發(fā)展對(duì)外貿(mào)易時(shí)，要提高環(huán)境的準(zhǔn)入門檻，借助跨國公司的環(huán)保、污染控制技術(shù)優(yōu)勢(shì)，促進(jìn)河北省各市產(chǎn)業(yè)結(jié)構(gòu)的優(yōu)化，改善霧霾污染狀況。

注釋：

①哥倫比亞大學(xué)社會(huì)經(jīng)濟(jì)數(shù)據(jù)和應(yīng)用中心（CIESIN）公布的關(guān)于PM2.5濃度值的全球年均值的遙感地圖只截止到2012年，而我國各地級(jí)市PM2.5數(shù)據(jù)的監(jiān)測(cè)始于2013年末，之前數(shù)據(jù)缺失嚴(yán)重，為了保證研究數(shù)據(jù)統(tǒng)計(jì)口徑的一致性，因此本文選取的PM2.5濃度值數(shù)據(jù)的時(shí)間范圍為2000—2012年。

參考文獻(xiàn)：

［1］PANAYOUTOU T.Empirical tests and policy analysis of environmental degradation at different stages of economic development［Z］.Ilo Working Papers，1993，4.

［2］SELDEN T M，SONG D Q.Neoclassical growth the J curve for abatement，and the inverted-U curve for pollution［J］.Journal of environmental economics and management，1995（29）：162-168.

［3］BRAJER V，MEAD R W，XIAO F.Searching for an environmental Kuznets curve in China's air pollution［J］.China economic review，2011，22（3）：383-397.

［4］包群，彭水軍.經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染：基于面板數(shù)據(jù)的聯(lián)立方程估計(jì)［J］.世界經(jīng)濟(jì)，2006（11）：48-58.

［5］張少華.經(jīng)濟(jì)全球化對(duì)我國環(huán)境污染影響的實(shí)證研究——基于行業(yè)面板數(shù)據(jù)［J］.國際貿(mào)易問題，2009（11）：68-79.

［6］張悅，趙曉丹.中國經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染關(guān)系的研究——基于環(huán)境庫茲涅茨曲線的實(shí)證分析［J］.經(jīng)濟(jì)研究導(dǎo)刊，2014（26）：13-16.

［7］郭俊華，劉奕瑋.我國城市霧霾天氣治理的產(chǎn)業(yè)結(jié)構(gòu)調(diào)整［J］.西北大學(xué)學(xué)報(bào)（哲學(xué)社會(huì)科學(xué)版），2014（3）：85-89.

［8］冷艷麗，杜思正.產(chǎn)業(yè)結(jié)構(gòu)、城市化與霧霾污染［J］.中國科技論壇，2015（9）：49-55.

［9］李鵬.產(chǎn)業(yè)結(jié)構(gòu)調(diào)整與環(huán)境污染之間存在倒U型曲線關(guān)系嗎?［J］.經(jīng)濟(jì)問題探索，2015（12）：56-67.

［10］馬麗梅，張曉.中國霧霾污染的空間效應(yīng)及經(jīng)濟(jì)、能源結(jié)構(gòu)影響［J］.中國工業(yè)經(jīng)濟(jì)，2014（4）：19-31.

［11］胡飛.產(chǎn)業(yè)結(jié)構(gòu)升級(jí)、對(duì)外貿(mào)易與環(huán)境污染的關(guān)系研究——以我國東部和中部地區(qū)為例［J］.經(jīng)濟(jì)問題探索，2011（7）：113-118.

［12］王家庭，王璇.我國城市化與環(huán)境污染的關(guān)系研究——基于28個(gè)省市面板數(shù)據(jù)的實(shí)證分析［J］.城市問題，2010（11）：9-15.

［13］劉伯龍，袁曉玲，張占軍.城鎮(zhèn)化推進(jìn)對(duì)霧霾污染的影響——基于中國省級(jí)動(dòng)態(tài)面板數(shù)據(jù)的經(jīng)驗(yàn)分析［J］.城市發(fā)展研究，2015（9）：23-27.

［14］潘慧峰，王鑫，張書宇.重霧霾污染的溢出效應(yīng)研究——來自京津冀地區(qū)的證據(jù)［J］.科學(xué)決策，2015（2）：1-15.

［15］潘慧峰，王鑫，張書宇.霧霾污染的持續(xù)性及空間溢出效應(yīng)分析——來自京津冀地區(qū)的證據(jù)［J］.中國軟科學(xué)，2015（12）：134-143.

［16］馮博，王雪青.考慮霧霾效應(yīng)的京津冀地區(qū)能源效率實(shí)證研究［J］.干旱區(qū)資源與環(huán)境，2015（10）：1-7.

［17］王穎，楊利花.跨界治理與霧霾治理轉(zhuǎn)型研究——以京津冀區(qū)域?yàn)槔跩］.東北大學(xué)學(xué)報(bào)（社會(huì)科學(xué)版），2016（7）：388-393.

［18］李瑞，蔡軍.河北工業(yè)結(jié)構(gòu)、能源消耗與霧霾關(guān)系探討［J］.宏觀經(jīng)濟(jì)管理，2014（5）：79-80.

［19］邵帥，李欣，曹建華，等.中國霧霾污染治理的經(jīng)濟(jì)政策選擇——基于空間溢出效應(yīng)的視角［J］.經(jīng)濟(jì)研究，2016（9）：73-88.

［20］胡宗義，劉亦文.能源消耗、污染排放與區(qū)域經(jīng)濟(jì)增長(zhǎng)關(guān)系的空間計(jì)量分析［J］.數(shù)理統(tǒng)計(jì)與管理，2015（1）：1-9.

［21］林光平，龍志，吳梅.我國地區(qū)經(jīng)濟(jì)收斂的空間計(jì)量實(shí)證分析：1978—2002 年［J］.經(jīng)濟(jì)學(xué)（季刊），2005（10）：67-82.

［22］KRUGMAN P.Increasing returns and economic geography［Z］.NBER Working Papers，1990，99（3）：483-499.

［23］COPELAND B R，TAYLOR M S.Trade，growth and environment［J］.Journal of economic literature，1994，42（1）：7-17.

［24］LIST J A，CO C Y.The effects of environmental regulations on foreign direct investment［J］.Journal of environmental economics&management，2000，40（1）：1-20.

［25］沙文兵，石濤.外商直接投資的環(huán)境效應(yīng)——基于中國省級(jí)面板數(shù)據(jù)的實(shí)證分析［J］.世界經(jīng)濟(jì)研究，2006（6）：76-81，89.

［26］冷艷麗，冼國明，杜思正.外商直接投資與霧霾污染——基于中國省際面板數(shù)據(jù)的實(shí)證分析［J］.國際貿(mào)易問題，2015（12）：74-84.

［27］童玉芬，王瑩瑩.中國城市人口與霧霾：相互作用機(jī)制路徑分析［J］.北京社會(huì)科學(xué)，2014（5）：4-10.

［28］秦蒙，劉修巖，仝怡婷.蔓延的城市空間是否加重了霧霾污染——來自中國PM2.5數(shù)據(jù)的經(jīng)驗(yàn)分析［J］.財(cái)貿(mào)經(jīng)濟(jì)，2016（11）：146-160.