河北經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染關(guān)系的檢驗(yàn)

杜 穎，肖榮閣，朱春華

（1.中國(guó)地質(zhì)大學(xué) （北京），北京100083；2.河北地質(zhì)職工大學(xué)，河北 石家莊 050081）

河北經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染關(guān)系的檢驗(yàn)

杜 穎1，2，肖榮閣1，朱春華1

（1.中國(guó)地質(zhì)大學(xué) （北京），北京100083；2.河北地質(zhì)職工大學(xué)，河北 石家莊 050081）

本文根據(jù)1991～2012年河北省五類(lèi)環(huán)境污染指標(biāo)和人均GDP時(shí)間序列數(shù)據(jù)，采用環(huán)境庫(kù)茲涅茨曲線、脈沖響應(yīng)分析和方差分解方法實(shí)證檢驗(yàn)了經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染之間的動(dòng)態(tài)關(guān)系。研究發(fā)現(xiàn)：隨著經(jīng)濟(jì)快速增長(zhǎng)，河北環(huán)境質(zhì)量總體趨于惡化，經(jīng)濟(jì)增長(zhǎng)加速了環(huán)境污染物的排放且持續(xù)時(shí)間較長(zhǎng)，但污染物排放量的增多并沒(méi)有大幅推動(dòng)經(jīng)濟(jì)增長(zhǎng)。經(jīng)濟(jì)增長(zhǎng)與環(huán)境治理之間的博弈是現(xiàn)階段河北發(fā)展面臨的最大難題，必須綜合利用排污權(quán)交易、稅收補(bǔ)貼等手段，加快技術(shù)革新、優(yōu)化產(chǎn)業(yè)結(jié)構(gòu)，實(shí)現(xiàn)經(jīng)濟(jì)環(huán)境協(xié)調(diào)發(fā)展。

經(jīng)濟(jì)增長(zhǎng)；環(huán)境污染；環(huán)境庫(kù)茲涅茨曲線；脈沖效應(yīng)分析；方差分解；河北

隨著經(jīng)濟(jì)增長(zhǎng)與環(huán)境保護(hù)的矛盾日益突出，國(guó)外學(xué)者大量利用計(jì)量經(jīng)濟(jì)學(xué)和統(tǒng)計(jì)學(xué)方法對(duì)經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染之間的關(guān)系展開(kāi)研究，較有代表性的是美國(guó) Grossman和 Krueger（1991）［1］提出的環(huán)境庫(kù)茲涅茨曲線（EKC），該理論將倒U型假說(shuō)引入經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染關(guān)系的研究中，認(rèn)為一個(gè)國(guó)家的污染程度會(huì)隨國(guó)民收入的增加而上升，當(dāng)國(guó)民收入達(dá)到某一轉(zhuǎn)折點(diǎn)后，環(huán)境污染水平會(huì)隨之下降。此后，各國(guó)學(xué)者利用截面數(shù)據(jù)、面板數(shù)據(jù)對(duì)EKC曲線進(jìn)行理論論證和實(shí)證檢驗(yàn)，如 Ekins（1997）［2］、Levinson等（2000）［3］以及 Dinda（2004）［4］等。針對(duì)我國(guó)經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染之間關(guān)系的研究起步較晚，包群和彭水軍（2006）［5］根據(jù)1996～2000年中國(guó)30個(gè)省市面板數(shù)據(jù)，實(shí)證檢驗(yàn)了大多省份的經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染之間存在EKC倒U曲線關(guān)系。原志華等（2008）［6］對(duì)山西進(jìn)行實(shí)證檢驗(yàn)，結(jié)果表明，除煙塵排放量外，環(huán)境指標(biāo)與經(jīng)濟(jì)增長(zhǎng)的關(guān)系呈波浪式的N型或倒N型曲線關(guān)系，山西省環(huán)境污染隨著人均GDP的增長(zhǎng)而不斷加劇。李鵬等（2009）［7］則發(fā)現(xiàn)山西經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染之間存在EKC倒U曲線關(guān)系。李琳等（2009）［8］運(yùn)用VAR模型考察了重慶經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染的相互影響機(jī)制，實(shí)證檢驗(yàn)了工業(yè)廢水、工業(yè)粉塵、工業(yè)SO2、工業(yè)固體廢棄物排放量的增加對(duì)重慶經(jīng)濟(jì)增長(zhǎng)存在著反作用力。陳桂月等（2012）［9］基于Granger因果檢驗(yàn)研究了內(nèi)蒙古經(jīng)濟(jì)發(fā)展與環(huán)境污染關(guān)系，發(fā)現(xiàn)內(nèi)蒙古各環(huán)境污染因子排放量的增多并沒(méi)有大幅促進(jìn)經(jīng)濟(jì)增長(zhǎng)，經(jīng)濟(jì)增長(zhǎng)卻帶動(dòng)了工業(yè)污染物的排放。王敏等（2015）［10］根據(jù)2003～2010年中國(guó)112個(gè)城市大氣污染濃度數(shù)據(jù)，認(rèn)為我國(guó)經(jīng)濟(jì)增長(zhǎng)和環(huán)境污染之間的呈現(xiàn)出“U形”曲線關(guān)系，但高增長(zhǎng)并不一定會(huì)導(dǎo)致高污染。上述研究各有側(cè)重，但對(duì)EKC曲線成立條件以及河北經(jīng)濟(jì)環(huán)境的相互作用機(jī)理，有待進(jìn)一步研究。

長(zhǎng)期以來(lái)，河北經(jīng)濟(jì)快速增長(zhǎng)的同時(shí)，環(huán)境狀況日益惡化，工業(yè)“三廢”特別是工業(yè)廢水和廢氣排放量明顯上升，帶來(lái)了嚴(yán)重的水污染和大氣霧霾，造成河北經(jīng)濟(jì)可持續(xù)發(fā)展的巨大桎梏。本文結(jié)合河北經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染之間的環(huán)境庫(kù)茲涅茨曲線的實(shí)證檢驗(yàn)，深入分析了河北經(jīng)濟(jì)環(huán)境的相互作用機(jī)制，通過(guò)脈沖響應(yīng)分析和方差分解得出工業(yè)三廢與經(jīng)濟(jì)增長(zhǎng)之間的動(dòng)態(tài)影響軌跡，為河北經(jīng)濟(jì)環(huán)境協(xié)調(diào)發(fā)展提供參考。

1 研究方法

1.1 數(shù)據(jù)來(lái)源

通常使用工業(yè)廢水排放量、工業(yè)廢氣排放量和工業(yè)固體廢物產(chǎn)生量三個(gè)指標(biāo)，即工業(yè)“三廢”指標(biāo)來(lái)反映某國(guó)（地區(qū)）的環(huán)境污染程度。本文延續(xù)這一慣例，并結(jié)合河北大氣污染嚴(yán)重的區(qū)域性特點(diǎn)，加入工業(yè)SO2排放量和工業(yè)煙粉塵排放量，共五類(lèi)污染指標(biāo)（工業(yè)廢水、工業(yè)廢氣、工業(yè)固體廢棄物、工業(yè)SO2、工業(yè)煙粉塵）反映環(huán)境污染狀況，選取實(shí)際人均GDP（1978年基期）反映經(jīng)濟(jì)增長(zhǎng)水平，數(shù)據(jù)時(shí)序長(zhǎng)度為1991～2012年，符號(hào)表示及數(shù)據(jù)來(lái)源詳見(jiàn)表1。鑒于時(shí)間序列數(shù)據(jù)取對(duì)數(shù)能較好消除序列的異方差，因此后續(xù)檢驗(yàn)均采用各變量的對(duì)數(shù)值。文中出現(xiàn)的dlny、dlnfs、dlnfq、dlnfw、dlnSO2、dlnyc分別表示一階差分后的lny、lnfs、lnfq、lnfw、lnSO2、lnyc；plny、llny分別表示lny的二次方和三次方。

表1 變量選取與數(shù)據(jù)來(lái)源

1.2 研究方法

1.2.1 環(huán)境庫(kù)茲涅茨曲線

環(huán)境庫(kù)茲涅茨曲線（EKC）表明，一國(guó)環(huán)境質(zhì)量最初隨著經(jīng)濟(jì)增長(zhǎng)而不斷惡化，越過(guò)拐點(diǎn)后，環(huán)境質(zhì)量才會(huì)隨之改善。本文采用經(jīng)典的GK模型［1］，引用三次多項(xiàng)式對(duì)數(shù)函數(shù)對(duì)河北省環(huán)境庫(kù)茲涅茨曲線進(jìn)行實(shí)證檢驗(yàn)，為使數(shù)據(jù)平穩(wěn)，函數(shù)關(guān)系中被解釋變量（工業(yè)廢水、工業(yè)廢氣、工業(yè)固體廢棄物、工業(yè)SO2、工業(yè)煙粉塵）和解釋變量（實(shí)際人均GDP）均以對(duì)數(shù)形式表示，見(jiàn)式（1）［11］。

ln（waste）t＝c＋β1lnyt＋β2lny21＋β3lny31＋ut（1）

式中，waste代表五類(lèi)污染排放指標(biāo)，y代表實(shí)際人均GDP，下標(biāo)t表示時(shí)間指標(biāo)，c為常數(shù)，u為隨機(jī)擾動(dòng)項(xiàng)。參數(shù)β1、β2、β3分別表示人均GDP的一次、二次和三次項(xiàng)系數(shù)，其取值范圍不同，則模型反映的EKC曲線關(guān)系也不同，具體包括七種關(guān)系（見(jiàn)表2）。如果擬合結(jié)果顯示三次項(xiàng)系數(shù)不顯著，則選用二次多項(xiàng)式對(duì)數(shù)函數(shù)重新回歸。

1.2.2 脈沖響應(yīng)分析和方差分解

脈沖響應(yīng)函數(shù)（impulse response functions，IRF）側(cè)重分析當(dāng)系統(tǒng)受到一個(gè)內(nèi)生變量的一個(gè)標(biāo)準(zhǔn)差沖擊時(shí)對(duì)系統(tǒng)內(nèi)各變量的動(dòng)態(tài)影響路徑，用于反映河北經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染之間的長(zhǎng)期相互作用關(guān)系。系統(tǒng)內(nèi)某變量的一個(gè)沖擊，不僅影響變量本身，還會(huì)傳遞給其他內(nèi)生變量，脈沖響應(yīng)分析可以較好的刻畫(huà)這些影響的動(dòng)態(tài)軌跡，展現(xiàn)一個(gè)變量的變動(dòng)是如何影響到系統(tǒng)內(nèi)其他變量的。

表2 EKC曲線判定標(biāo)準(zhǔn)

方差分解（variance decomposition）用于進(jìn)一步分析河北經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染之間影響的貢獻(xiàn)度。方差分解通過(guò)分析每個(gè)環(huán)境污染變量沖擊對(duì)人均GDP變化的貢獻(xiàn)度，進(jìn)而評(píng)價(jià)各污染變量沖擊對(duì)經(jīng)濟(jì)增長(zhǎng)的重要性，反之亦然。脈沖響應(yīng)和方差分解的具體方法見(jiàn)高鐵梅（2006）［12］。

2 協(xié)整檢驗(yàn)

為避免計(jì)量經(jīng)濟(jì)學(xué)中的偽回歸問(wèn)題，需要各變量具備相同的單整階數(shù)，同時(shí)，變量之間存在協(xié)整關(guān)系，因此，要對(duì)各環(huán)境污染指標(biāo)和經(jīng)濟(jì)指標(biāo)進(jìn)行平穩(wěn)性檢驗(yàn)和協(xié)整關(guān)系檢驗(yàn)。

當(dāng)各變量具有相同的單整階數(shù)時(shí)，則其為平穩(wěn)序列。本文采用單位根檢驗(yàn)法中的ADF檢驗(yàn)，它以“變量是非平穩(wěn)序列”為原假設(shè)、“變量是平穩(wěn)序列”為備選假設(shè)的左側(cè)單邊檢驗(yàn)。借助于eviews7.0對(duì)河北各環(huán)境變量和人均GDP進(jìn)行ADF檢驗(yàn)，結(jié)果如表3所示，除lny，各環(huán)境變量的ADF檢驗(yàn)值的絕對(duì)值均小于5%臨界值，變量序列不平穩(wěn)。一階差分后，各變量在5%顯著水平下均拒絕接受存在單位根的假設(shè)，都是平穩(wěn)序列。因此，各變量是一階單整序列，可以繼續(xù)進(jìn)行協(xié)整檢驗(yàn)。

在通過(guò)平穩(wěn)性檢驗(yàn)的前提下，進(jìn)行變量間的協(xié)整關(guān)系檢驗(yàn)。本文采用Johansen檢驗(yàn)，較好適用于多變量檢驗(yàn)。檢驗(yàn)結(jié)果見(jiàn)表4，特征根跡檢驗(yàn)和最大特征值檢驗(yàn)均表明，在5%顯著性水平下，各環(huán)境污染指標(biāo)和經(jīng)濟(jì)指標(biāo)變量之間存在協(xié)整關(guān)系，不存在偽回歸問(wèn)題。

3 實(shí)證分析

3.1 河北環(huán)境庫(kù)茲涅茨曲線的回歸分析

利用河北1991～2012年各污染數(shù)據(jù)和人均GDP數(shù)據(jù)，借助于eviews7.0對(duì)式（1）進(jìn)行環(huán)境庫(kù)茲涅茨曲線回歸分析，參數(shù)估計(jì)結(jié)果見(jiàn)表5，擬合曲線見(jiàn)圖1。

表3 ADF單位根檢驗(yàn)結(jié)果

表4 Johansen協(xié)整檢驗(yàn)結(jié)果

在5%顯著性水平下參數(shù)β1＞0、β2＜0、β3＞0，可知lny與lnfs呈N型曲線關(guān)系，圖1的擬合曲線也驗(yàn)證了兩者的N型關(guān)系。從擬合曲線看，工業(yè)廢水排放量先降后升，且當(dāng)前處于右半部分上升階段，說(shuō)明河北省工業(yè)廢水排放量正隨著人均GDP的增加而增加，環(huán)境污染加劇。

lny與lnfq呈倒N型曲線關(guān)系，工業(yè)廢氣排放量先降后升，且當(dāng)前處于右半部分上升階段，預(yù)計(jì)未來(lái)河北工業(yè)廢氣排放量將隨著人均GDP的增加而增加，但趨勢(shì)減緩，大氣污染壓力有所減輕。

lny與lnfw呈U型曲線關(guān)系，工業(yè)廢棄物產(chǎn)生量先降后升，且當(dāng)前處于右半部分上升階段，說(shuō)明河北工業(yè)廢棄物排放量將隨著人均GDP的增加而增加，環(huán)境質(zhì)量趨于惡化。

lny與lnSO2呈倒U型曲線關(guān)系，工業(yè)SO2排放量先升后降，且當(dāng)前已進(jìn)入右半部分下降階段，說(shuō)明河北工業(yè)SO2排放量將隨著人均GDP的增加而減緩，環(huán)境壓力有所減輕。

lny與lnyc呈倒U型曲線關(guān)系，工業(yè)煙粉塵排放量先升后降，且當(dāng)前已進(jìn)入右半部分下降階段，說(shuō)明河北工業(yè)煙粉塵排放量將隨著人均GDP的增加而減緩，環(huán)境壓力有所減輕。

綜上可知，隨著河北省人均GDP的不斷增長(zhǎng)，環(huán)境質(zhì)量總體趨于惡化，河北經(jīng)濟(jì)增長(zhǎng)的環(huán)境污染代價(jià)較大。

表5 河北省環(huán)境庫(kù)茲涅茨曲線參數(shù)估計(jì)結(jié)果

圖1 河北省各環(huán)境污染指標(biāo)與人均GDP的EKC擬合曲線

3.2 脈沖響應(yīng)分析結(jié)果

分析河北人均GDP對(duì)各環(huán)境污染指標(biāo)排放的沖擊響應(yīng)，得到的脈沖響應(yīng)函數(shù)圖（圖2），橫軸表示沖擊作用的滯后期間，本文設(shè)定為10期，縱軸表示各環(huán)境污染指標(biāo)的響應(yīng)程度。

圖2（a）表明，當(dāng)本期給人均GDP一個(gè)正沖擊后，工業(yè)廢水排放量在前4期為正，第4～6期為負(fù)，隨后上下波動(dòng)于第7期達(dá)到最高點(diǎn)0.066520，累計(jì)響應(yīng)值0.17975。這表明人均GDP受外部某一沖擊后，經(jīng)傳遞后給工業(yè)廢水排放量帶來(lái)同向沖擊，且這一沖擊具有較強(qiáng)的正效應(yīng)。

圖2（b）表明，當(dāng)本期給人均GDP一個(gè)正沖擊后，工業(yè)廢氣排放量在前4期小幅波動(dòng)，第4期達(dá)到最高點(diǎn)0.064449，隨后緩慢下降但始終為正影響，累計(jì)響應(yīng)值0.451022。這表明人均GDP受外部某一沖擊后，經(jīng)傳遞后給工業(yè)廢氣排放量帶來(lái)同向沖擊，且這一沖擊對(duì)工業(yè)廢氣排放量具有顯著的促進(jìn)作用，并且這一顯著促進(jìn)作用具有較長(zhǎng)的持續(xù)效應(yīng)。

圖2（c）表明，當(dāng)本期給人均GDP一個(gè)正沖擊后，工業(yè)廢棄物產(chǎn)生量在前3期大幅攀升，第5期達(dá)到最高點(diǎn)0.055353，隨后保持穩(wěn)定，累計(jì)響應(yīng)值0.463803。這表明人均GDP受外部某一沖擊后，經(jīng)傳遞后給工業(yè)廢棄物產(chǎn)生量帶來(lái)同向沖擊，且這一沖擊對(duì)工業(yè)廢棄物產(chǎn)生量具有顯著的促進(jìn)作用，并且這一顯著促進(jìn)作用具有較長(zhǎng)的持續(xù)效應(yīng)。

圖2（d）表明，當(dāng)本期給人均GDP一個(gè)正沖擊后，工業(yè)SO2排放量第一期為0，前3期上升，第4期達(dá)到最高點(diǎn)0.006933，隨后保持穩(wěn)定，累計(jì)響應(yīng)值0.051913。這表明人均GDP受外部某一沖擊后，經(jīng)傳遞后給工業(yè)SO2排放量帶來(lái)同向沖擊，且這一沖擊對(duì)工業(yè)SO2排放量具有較弱的促進(jìn)作用，但持續(xù)效應(yīng)較長(zhǎng)。

圖2（e）表明，當(dāng)本期給人均GDP一個(gè)正沖擊后，工業(yè)煙粉塵排放量隨之下降，第1期為0，第3期達(dá)到最低點(diǎn)－0.06616，之后大幅回彈至接近0點(diǎn)，累計(jì)響應(yīng)值－0.27844。這表明人均GDP受外部某一沖擊后，經(jīng)傳遞后給工業(yè)煙粉塵排放量帶來(lái)反向沖擊，且這一沖擊對(duì)工業(yè)煙粉塵排放量具有較強(qiáng)的抑制作用，但這一抑制作用的持續(xù)時(shí)間較短。

反過(guò)來(lái)分析各環(huán)境污染指標(biāo)的變動(dòng)對(duì)人均GDP的影響，圖3橫軸表示沖擊作用的滯后期間，縱軸表示人均GDP的響應(yīng)程度。

圖3（a）表明，當(dāng)本期給工業(yè)廢水排放量一個(gè)正沖擊后，人均GDP持續(xù)上升，在第4～7期小幅上下波動(dòng)之后于第9期達(dá)到最高點(diǎn)0.039531，累計(jì)響應(yīng)值0.226038。這表明工業(yè)廢水排放量受外部某一沖擊后，經(jīng)傳遞后給人均GDP帶來(lái)同向沖擊，且這一沖擊具有持續(xù)顯著的促進(jìn)作用。

由圖3（b）表明，當(dāng)本期給工業(yè)廢氣排放量一個(gè)正沖擊后，人均GDP大幅上升，第7期達(dá)到最高點(diǎn)0.033107，隨后小幅下滑，累計(jì)響應(yīng)值0.278383。這表明工業(yè)廢氣排放量受外部某一沖擊后，經(jīng)傳遞后給人均GDP帶來(lái)同向沖擊，且這一沖擊具有極強(qiáng)的促進(jìn)作用和較長(zhǎng)的持續(xù)效應(yīng)。

圖3（c）表明，當(dāng)本期給工業(yè)廢棄物產(chǎn)生量一個(gè)正沖擊后，人均GDP持續(xù)上升，累計(jì)響應(yīng)值0.173458。這表明工業(yè)廢棄物產(chǎn)生量受外部某一沖擊后，經(jīng)傳遞后給人均GDP帶來(lái)同向沖擊，且這一沖擊具有非常顯著的促進(jìn)作用和較長(zhǎng)的持續(xù)效應(yīng)。

圖3（d）表明，當(dāng)本期給工業(yè)SO2排放量一個(gè)正沖擊后，人均GDP持續(xù)下滑，第2期降為負(fù)值，累計(jì)響應(yīng)值－0.09274。這表明工業(yè)SO2排放量受外部某一沖擊后，經(jīng)傳遞后給人均GDP帶來(lái)反向沖擊，且這一沖擊具有較弱的抑制作用和較長(zhǎng)的持續(xù)效應(yīng)。

圖3（e）表明，當(dāng)本期給工業(yè)煙粉塵排放量一個(gè)正沖擊后，人均GDP持續(xù)下滑，第4期降為最低點(diǎn)－0.007，隨后略有上升，累計(jì)響應(yīng)值－0.05439。這表明工業(yè)煙粉塵產(chǎn)生量受外部某一沖擊后，經(jīng)傳遞后給人均GDP帶來(lái)反向沖擊，且這一沖擊具有較弱的抑制作用和較長(zhǎng)的持續(xù)效應(yīng)。

3.3 方差分解分析結(jié)果

通過(guò)方差分解，可知各環(huán)境污染變量對(duì)河北經(jīng)濟(jì)增長(zhǎng)的貢獻(xiàn)率，如圖4所示，不考慮人均GDP自身的貢獻(xiàn)率，工業(yè)廢氣排放量和工業(yè)廢棄物產(chǎn)生量對(duì)人均GDP的貢獻(xiàn)度最大，且具有長(zhǎng)期貢獻(xiàn)，第9期分別達(dá)到最高點(diǎn)22.7%和22.6%；其次是工業(yè)煙粉塵，第1期為0，隨后大幅上升到第5期的15.5%，第10期達(dá)到最大值16.8%，貢獻(xiàn)逐步增大且持續(xù)時(shí)間長(zhǎng)；工業(yè)廢水排放量的貢獻(xiàn)度從第1期的0上升到第3期的最高點(diǎn)9.3%，隨后圍繞6%上下波動(dòng)；工業(yè)SO2的貢獻(xiàn)度最低，10期內(nèi)基本穩(wěn)定在2%左右。

反過(guò)來(lái)分析河北經(jīng)濟(jì)增長(zhǎng)對(duì)各環(huán)境污染指標(biāo)的貢獻(xiàn)度。圖5表明，不考慮環(huán)境污染指標(biāo)自身的貢獻(xiàn)率，人均GDP對(duì)工業(yè)廢氣排放量的貢獻(xiàn)率最大，且長(zhǎng)期貢獻(xiàn)較強(qiáng)，從第1期60%大幅上升到第2期70%，隨后基本穩(wěn)定；其次是工業(yè)廢水排放量，人均GDP對(duì)其貢獻(xiàn)率在前4期穩(wěn)定上升，第5期達(dá)到最高點(diǎn)40%左右。隨后在33%～40%之間上下波動(dòng)；人均GDP對(duì)工業(yè)廢棄物產(chǎn)生量的貢獻(xiàn)率逐步上升，第10期達(dá)到最高點(diǎn)30.8%；對(duì)工業(yè)SO2和工業(yè)煙粉塵排放量的貢獻(xiàn)率相對(duì)較小，最大值分別為10%和4.8%。

圖3 河北人均GDP對(duì)各污染指標(biāo)的沖擊響應(yīng)函數(shù)

圖4 河北各環(huán)境污染指標(biāo)對(duì)人均GDP的貢獻(xiàn)率

圖5 河北人均GDP對(duì)各環(huán)境污染指標(biāo)的貢獻(xiàn)率

4 結(jié) 論

本文采用EKC、脈沖響應(yīng)分析和方差分解方法，利用河北省1991～2012年時(shí)間序列數(shù)據(jù)，實(shí)證檢驗(yàn)了河北環(huán)境污染與經(jīng)濟(jì)增長(zhǎng)之間的EKC曲線關(guān)系、相互影響軌跡和貢獻(xiàn)度，結(jié)論總結(jié)如下。

1）環(huán)境庫(kù)茲涅茨曲線回歸檢驗(yàn)表明：1991～2012年間，河北工業(yè)廢水排放量與人均GDP呈N型曲線關(guān)系，且正處于上升階段，環(huán)境質(zhì)量趨于惡化。工業(yè)廢氣排放量與人均GDP呈倒N型曲線關(guān)系，且正處于上升階段，環(huán)境質(zhì)量惡化。工業(yè)廢棄物產(chǎn)生量與人均GDP呈U型曲線關(guān)系，且正處于上升階段，環(huán)境質(zhì)量趨于惡化。工業(yè)SO2排放量與人均GDP之間呈倒U型曲線關(guān)系，且正處于下降階段，環(huán)境壓力有所減輕。工業(yè)煙粉塵與人均GDP之間呈倒U型曲線關(guān)系，且正處于下降階段，環(huán)境壓力有所減輕。綜合來(lái)看，隨著經(jīng)濟(jì)的快速增長(zhǎng)，河北省環(huán)境質(zhì)量總體趨于惡化，環(huán)境污染代價(jià)較大。

2）脈沖響應(yīng)分析結(jié)果表明：河北經(jīng)濟(jì)增長(zhǎng)對(duì)工業(yè)廢水、廢氣排放和工業(yè)廢棄物產(chǎn)生具有持續(xù)較強(qiáng)的正面影響；對(duì)工業(yè)SO2排放具有持續(xù)較弱的正面影響；對(duì)工業(yè)煙粉塵排放具有較強(qiáng)的負(fù)面影響。工業(yè)廢水排放、工業(yè)廢氣排放、工業(yè)廢棄物產(chǎn)生對(duì)經(jīng)濟(jì)增長(zhǎng)產(chǎn)生持續(xù)較強(qiáng)的正面影響，工業(yè)SO2排放量對(duì)經(jīng)濟(jì)增長(zhǎng)具有持續(xù)較弱的負(fù)面影響，工業(yè)煙粉塵排放量對(duì)經(jīng)濟(jì)增長(zhǎng)具有持續(xù)較弱的負(fù)面影響。

3）方差分解分析表明：河北工業(yè)廢氣和廢棄物產(chǎn)生量對(duì)經(jīng)濟(jì)增長(zhǎng)的貢獻(xiàn)度最大，達(dá)22%左右，且具有長(zhǎng)期貢獻(xiàn)；其次是工業(yè)煙粉塵，貢獻(xiàn)度在15%左右，影響逐步增大且持續(xù)時(shí)間長(zhǎng)；工業(yè)廢水和SO2排放量對(duì)經(jīng)濟(jì)增長(zhǎng)的長(zhǎng)期貢獻(xiàn)度分別在6%和2%左右。反過(guò)來(lái)看經(jīng)濟(jì)增長(zhǎng)對(duì)各環(huán)境污染指標(biāo)的貢獻(xiàn)度，河北經(jīng)濟(jì)增長(zhǎng)對(duì)廢氣排放量的貢獻(xiàn)率，且長(zhǎng)期貢獻(xiàn)最大，高達(dá)70%；對(duì)工業(yè)廢水和工業(yè)廢棄物的貢獻(xiàn)率分別為40%和31%；對(duì)工業(yè)SO2和工業(yè)煙粉塵排放量的貢獻(xiàn)率則相對(duì)較小。簡(jiǎn)而言之，河北各環(huán)境污染指標(biāo)排放量（產(chǎn)生量）的增大并沒(méi)有大幅推動(dòng)經(jīng)濟(jì)增長(zhǎng)，但經(jīng)濟(jì)增長(zhǎng)對(duì)環(huán)境污染的貢獻(xiàn)度較高、解釋力較強(qiáng)。

4）綜合脈沖效應(yīng)分析和方差分解結(jié)果可知，河北經(jīng)濟(jì)的快速增長(zhǎng)加劇了環(huán)境污染，這也在一定程度上解釋了河北面臨的大氣污染治理困局。由于工業(yè)廢氣排放與經(jīng)濟(jì)增長(zhǎng)相互作用的持續(xù)性較強(qiáng)，意味著短期內(nèi)難以解決河北霧霾問(wèn)題，大氣治理將是一項(xiàng)長(zhǎng)期任務(wù)，但鑒于工業(yè)廢氣與經(jīng)濟(jì)增長(zhǎng)之間的正效應(yīng)較強(qiáng)，河北若減少工業(yè)廢氣排放，將會(huì)明顯抑制經(jīng)濟(jì)增長(zhǎng)。

河北工業(yè)發(fā)展在推動(dòng)經(jīng)濟(jì)增長(zhǎng)的同時(shí)，導(dǎo)致環(huán)境污染不斷加劇，工業(yè)“三廢”指標(biāo)均呈上升趨勢(shì)，目前呈現(xiàn)以犧牲環(huán)境為代價(jià)的經(jīng)濟(jì)增長(zhǎng)局面，經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染之間的博弈將是河北未來(lái)經(jīng)濟(jì)環(huán)境協(xié)調(diào)發(fā)展面臨的最大難題。在河北環(huán)境尚未遭到完全破壞的當(dāng)下，必須采取有力措施，加快技術(shù)革新、優(yōu)化產(chǎn)業(yè)結(jié)構(gòu)、實(shí)施節(jié)能減排戰(zhàn)略，將環(huán)境保護(hù)與經(jīng)濟(jì)可持續(xù)發(fā)展有機(jī)結(jié)合起來(lái)，既要將環(huán)境污染納入經(jīng)濟(jì)發(fā)展質(zhì)量評(píng)測(cè)之中，盡力減輕經(jīng)濟(jì)增長(zhǎng)對(duì)環(huán)境的破壞程度，又要通過(guò)資源定價(jià)、排污權(quán)交易、稅收補(bǔ)貼等手段多管齊下，盡快達(dá)到環(huán)境庫(kù)茲涅茨曲線的轉(zhuǎn)折點(diǎn)，早日實(shí)現(xiàn)河北省經(jīng)濟(jì)與環(huán)境協(xié)調(diào)發(fā)展。

［1］Grossman M，Krueger B.Environmental impacts of a North American free trade agreement ［D］.National Bureau of Economic Research Working Paper 3914，London：Cambridge MA，1991.

［2］Ekins P.The Kuznets curve for the environment and economic growth：examining the evidence［J］.Environment and planning，1997，29：805－829.

［3］Levinson A，Hilton F.Factoring the environmental Kuznets curve：evidence from automotive lead emission［J］.Journal of environmental economics and management，2000，35：126－141.

［4］Dinda S.Environmental Kuznets curve hypothesis：a survey［J］.Ecological economics，2004，49：431－455.

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

［6］原志華，延軍平，劉宇峰.山西省經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染水平計(jì)量模型研究［J］.農(nóng)業(yè)系統(tǒng)科學(xué)與綜合研究，2008，24（2）：158－162.

［7］李鵬，蘭宜生.經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染關(guān)系的研究——基于山西省數(shù)據(jù)的實(shí)證分析［J］.對(duì)外經(jīng)濟(jì)貿(mào)易大學(xué)學(xué)報(bào)，2009，6：54－62.

［8］李琳，朱金山，高潤(rùn)霞.基于JLU模型的重慶市經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染的關(guān)系研究［J］.西南大學(xué)學(xué)報(bào)：自然科學(xué)版，2009，31（11）：92－96.

［9］陳桂月，李海濤，梁濤.內(nèi)蒙古經(jīng)濟(jì)發(fā)展與環(huán)境污染之間關(guān)系的檢驗(yàn)［J］.自然資源學(xué)報(bào)，2012，27（11）：1845－1859.

［10］王敏，黃瀅.中國(guó)的環(huán)境污染與經(jīng)濟(jì)增長(zhǎng)［J］.經(jīng)濟(jì)學(xué)，2015，2：557－578.

［11］周璇，孫慧.中國(guó)工業(yè)廢水排放量與經(jīng)濟(jì)增長(zhǎng)關(guān)系的區(qū)域分異研究［J］.干旱區(qū)資源與環(huán)境，2013，27（12）：15－19.

［12］高鐵梅.計(jì)量經(jīng)濟(jì)分析方法與建模［M］.北京：清華大學(xué)出版社，2006：277－314.

An empirical analysis on the relationship between economic growth and environmental pollution in Hebei

DU Ying1，XIAO Rong－ge1，ZHU Chun－h(huán)ua1
（1.China University of Geosciences（Beijing），Beijing 100083，China；2.Hebei University of Geological Workers，Shijiazhuang 050081，China）

Based on the time series data of five pollution indices and per capita GDP in Hebei from 1991 to 2012，this paper investigated the interactions between economic growth and environmental pollution through EKC，impulse response functions and variance decomposition.Results indicate：The environment becomes deteriorated with the rapid economic growth of Hebei.The economic growth could push up pollution level for a long time while the growth of pollutant emissions cannot significantly promote the economy.At this stage，the most challenging problem of Hebei economy is the competition between economic growth and environmental management.To solve this problem and realize the coordinated development of economy and environment，we should speed up technical innovation and optimize economic structure using a combination of emission trading and tax subsidies.

economic growth；environmental pollution；environmental Kuznets curve；impulse response analysis；variance decomposition；Hebei province

肖榮閣（1949－），男，漢族，吉林舒蘭人，中國(guó)地質(zhì)大學(xué)（北京）教授，博士生導(dǎo)師，主要從事礦業(yè)勘查方面的教學(xué)與研究工作。E－mail：rgxiao＠cugb.edu.cn。

X196

A

1004－4051（2015）09－0062－07

2014－11－28

河北省教育廳社科研究2014年度基金項(xiàng)目資助（編號(hào)：SZ141209）

杜穎（1986－），女，漢族，河北藁城人，中國(guó)地質(zhì)大學(xué)（北京）資源產(chǎn)業(yè)經(jīng)濟(jì)專(zhuān)業(yè)博士研究生，主要從事資源經(jīng)濟(jì)的研究工作。E－mail：148121960＠qq.com。