山東省污染物排放與經(jīng)濟(jì)發(fā)展關(guān)系研究

劉金玉

摘 要：運(yùn)用多元線性回歸和環(huán)境庫茲涅茨曲線模型模擬1981—2015年山東省環(huán)境質(zhì)量與經(jīng)濟(jì)發(fā)展之間的關(guān)系。研究結(jié)果表明，山東省SO2排放量與人均GDP的關(guān)系不顯著，SO2排放量隨經(jīng)濟(jì)發(fā)展先急劇上升后波浪式下降;煙粉塵排放量與人均GDP符合U型模型，工業(yè)廢水、固體廢物排放量與人均GDP存在倒U型關(guān)系。山東省應(yīng)優(yōu)化產(chǎn)業(yè)結(jié)構(gòu)，提高公眾環(huán)境保護(hù)意識(shí)，加強(qiáng)環(huán)境治理，加大投資力度，避免環(huán)境污染程度隨經(jīng)濟(jì)發(fā)展出現(xiàn)反復(fù)惡化的情況。

關(guān)鍵詞：山東省;污染物排放;經(jīng)濟(jì)發(fā)展;多元線性回歸分析;環(huán)境庫茲涅茨曲線

文章編號(hào)：1004-7026（2019）17-0116-02? ? ? ? ?中國圖書分類號(hào)：F019.3 ? ? ? ? 文獻(xiàn)標(biāo)志碼：A

隨著我國經(jīng)濟(jì)高速發(fā)展，環(huán)境污染問題日益加重，對(duì)生態(tài)系統(tǒng)和人類社會(huì)造成直接或間接的破壞，是制約經(jīng)濟(jì)發(fā)展的重要因素。

目前，我國經(jīng)濟(jì)正處在由高速發(fā)展向高質(zhì)量發(fā)展的重要轉(zhuǎn)折時(shí)期，肩負(fù)著發(fā)展經(jīng)濟(jì)和保護(hù)環(huán)境的雙重任務(wù)。因此，研究環(huán)境污染與經(jīng)濟(jì)發(fā)展的關(guān)系具有重大意義[1-3]。

不同國家和地區(qū)的經(jīng)濟(jì)發(fā)展水平和生態(tài)環(huán)境狀況不同，所選取的研究指標(biāo)不同，EKC曲線的形式也不同，有倒U型、U型、倒U+U型、N型等[4-8]。

為了研究山東省環(huán)境質(zhì)量與經(jīng)濟(jì)發(fā)展水平之間的關(guān)系，以社會(huì)經(jīng)濟(jì)指標(biāo)和環(huán)境污染指標(biāo)進(jìn)行多元線性回歸分析和環(huán)境庫茲涅茨曲線模擬（EKC），定量研究山東省1981—2015年環(huán)境污染與社會(huì)經(jīng)濟(jì)發(fā)展之間的關(guān)系，為山東省環(huán)境治理和經(jīng)濟(jì)調(diào)控提供理論支撐。

1? 研究區(qū)概況

山東省位于中國東部沿海，是環(huán)渤海經(jīng)濟(jì)圈的重要組成部分，是中國經(jīng)濟(jì)第3大省，2015年山東省地區(qū)生產(chǎn)總值占全國GDP的9.2%。山東省經(jīng)濟(jì)發(fā)展水平不斷提升，2000年以前呈現(xiàn)平穩(wěn)緩慢增長(zhǎng)趨勢(shì)，2000年以后地區(qū)生產(chǎn)總值呈現(xiàn)高速增長(zhǎng)趨勢(shì)，1981—2015年山東省GDP年均增長(zhǎng)率為16.54%。近年來，山東省產(chǎn)業(yè)結(jié)構(gòu)表現(xiàn)為第一產(chǎn)業(yè)的比重持續(xù)下降，第二和第三產(chǎn)業(yè)的比重緩慢增長(zhǎng)并基本持平。第一、二、三產(chǎn)業(yè)結(jié)構(gòu)由1981年的38.2%、44.8%和17.0%，變?yōu)?015年的7.9%、46.8%和45.3%。山東省經(jīng)濟(jì)水平的迅速發(fā)展，使工業(yè)“三廢”的排放持續(xù)增加，加劇了環(huán)境污染，造成了環(huán)境質(zhì)量惡化。

2? 研究方法

2.1? 多元線性回歸分析

運(yùn)用多元線性回歸分析，從總?cè)丝跀?shù)、地區(qū)生產(chǎn)總值和城鎮(zhèn)化水平3個(gè)方面，進(jìn)行社會(huì)經(jīng)濟(jì)發(fā)展速度與環(huán)境污染物排放量關(guān)系的定量研究。多元回歸分析的模型如下。

Pol_i=β_1i+β_2iPopu_i+β_3iGDP_i+β_4iUrb_i+ε_i（1）

式中：i=1表示二氧化硫，i=2表示煙粉塵，i=3表示工業(yè)廢水，i=4表示工業(yè)固體廢物;Pol為環(huán)境污染物排放量;Popu為人口數(shù)，GDP為地區(qū)生產(chǎn)總值，Urb代表城市化率。

2.2? EKC曲線模擬

EKC曲線描述了經(jīng)濟(jì)發(fā)展與環(huán)境污染程度呈倒U型關(guān)系。對(duì)山東省環(huán)境污染物排放量與經(jīng)濟(jì)發(fā)展程度間的關(guān)系進(jìn)行EKC曲線模擬，選擇人均GDP為經(jīng)濟(jì)增長(zhǎng)指標(biāo)，二氧化硫、煙粉塵、工業(yè)廢水、工業(yè)固體廢物的排放量為環(huán)境污染指標(biāo)，分別使用對(duì)數(shù)函數(shù)、二次函數(shù)、三次函數(shù)進(jìn)行曲線回歸模擬，根據(jù)擬合優(yōu)度和參數(shù)檢驗(yàn)的顯著性結(jié)果選擇最佳模型。

3? 結(jié)果與分析

3.1? 多元回歸結(jié)果

二氧化硫排放量的回歸方程為Pol₁=-263.300+1.199Popu₁-1.202GDP₁

煙粉塵排放量的回歸方程為Pol₂=132.496-0.535Urb₂

工業(yè)廢水排放量的回歸方程為Pol₃=375230.296-0.830Popu₃+0.385GDP₃+1.354Urb₃

工業(yè)固體廢物排放量的回歸方程為Pol₄=1393.864+0.828GDP₄+0.177Urb₄

模型輸出結(jié)果表明：①總?cè)丝诤虶DP對(duì)二氧化硫排放量的影響是顯著的。地區(qū)生產(chǎn)總值與二氧化硫排放量呈負(fù)相關(guān)，說明隨著經(jīng)濟(jì)的發(fā)展，沒有以犧牲環(huán)境換取經(jīng)濟(jì)增長(zhǎng)，反而加強(qiáng)環(huán)境治理，使生態(tài)環(huán)境質(zhì)量持續(xù)好轉(zhuǎn)。②煙粉塵排放量與城市化率呈負(fù)相關(guān)，說明隨著山東省城市化進(jìn)程的加快，煙粉塵排放量減少，環(huán)境質(zhì)量變好。③工業(yè)廢水排放量與社會(huì)經(jīng)濟(jì)指標(biāo)的各個(gè)變量都具有顯著得相關(guān)性。④工業(yè)固體廢物排放量與地區(qū)生產(chǎn)總值、城市化率呈正相關(guān)，說明經(jīng)濟(jì)發(fā)展導(dǎo)致固體廢物的增加，未來要發(fā)展循環(huán)經(jīng)濟(jì)，減少“三廢”排放。

3.2? EKC曲線模擬結(jié)果

3.2.1? 二氧化硫與人均GDP的EKC檢驗(yàn)

二氧化硫與人均GDP間沒有擬合較好的模型，伴隨經(jīng)濟(jì)的發(fā)展，二氧化硫排放量呈現(xiàn)先急劇上升后波浪式下降的趨勢(shì)，具有非常明顯的波動(dòng)特征，這說明政府對(duì)二氧化硫排放的治理存在“短視”現(xiàn)象。

3.2.2? 煙粉塵與人均GDP的EKC檢驗(yàn)

煙粉塵排放量與人均GDP關(guān)系可取二次擬合曲線Y=7.389×10^-8X²-0.005X+122.698，說明隨著經(jīng)濟(jì)的不斷發(fā)展，煙粉塵排放量與人均GDP存在U型關(guān)系，當(dāng)前煙粉塵排放處在U型曲線的右半邊呈現(xiàn)上升趨勢(shì)，結(jié)合山東統(tǒng)計(jì)年鑒統(tǒng)計(jì)口徑的變化和數(shù)據(jù)本身分析，實(shí)際情況應(yīng)與擬合曲線有所不同。目前，隨著環(huán)境政策的細(xì)化和人們環(huán)保意識(shí)的增強(qiáng)，煙粉塵的排放量基本能維持當(dāng)前水平或進(jìn)一步下降。

3.2.3? 工業(yè)廢水與人均GDP的EKC檢驗(yàn)

工業(yè)廢水排放與人均GDP關(guān)系為二次函數(shù)，函數(shù)方程為Y=-3.260×10^-5X²+3.668X+8.471×10⁴。工業(yè)廢水排放量與經(jīng)濟(jì)水平的關(guān)系為倒U型曲線，2012年左右拐點(diǎn)在人均GDP50 000元。2012年以前，工業(yè)廢水排放量隨人均GDP的增長(zhǎng)而增長(zhǎng)，到達(dá)50 000元臨界點(diǎn)后，再隨人均GDP的增長(zhǎng)逐漸減少。由于近幾年第三產(chǎn)業(yè)的快速發(fā)展及科學(xué)技術(shù)的進(jìn)步，工業(yè)廢水排放量得到控制，并在未來隨經(jīng)濟(jì)發(fā)展有可能會(huì)呈下降趨勢(shì)，應(yīng)積極采取措施進(jìn)行治理廢水，加強(qiáng)管理工業(yè)廢水排放源頭，增加廢水處理設(shè)備的資金投入等。

3.2.4? 工業(yè)固體廢物與人均GDP的EKC檢驗(yàn)

工業(yè)固體廢物排放與人均GDP的關(guān)系為二次方程Y=-1.819×10^-6X²+0.394X+2.466×10³。工業(yè)固體廢物排放與人均GDP近似呈倒U型關(guān)系，當(dāng)前固體廢物排放呈上升的趨勢(shì)，但上升趨勢(shì)有所放緩，未來應(yīng)提高固體廢物的綜合利用率，發(fā)展循環(huán)經(jīng)濟(jì)，使拐點(diǎn)提前出現(xiàn)。

4? 結(jié)論

（1）二氧化硫排放量隨人均GDP增長(zhǎng)的變化趨勢(shì)為先急劇上升后波浪式下降;煙粉塵排放量隨人均GDP增長(zhǎng)先減小后增加，符合U型曲線;工業(yè)廢水排放量隨經(jīng)濟(jì)的發(fā)展表現(xiàn)為先上升后下降的變化趨勢(shì)，符合倒U型曲線;工業(yè)固體廢物排放量隨經(jīng)濟(jì)發(fā)展表現(xiàn)為持續(xù)增長(zhǎng)趨勢(shì)，未來即將越過拐點(diǎn)改變趨勢(shì)，近似于倒U型曲線。

（2）環(huán)境庫茲涅茨曲線只是一種現(xiàn)象，人口、經(jīng)濟(jì)政策、產(chǎn)業(yè)結(jié)構(gòu)、科學(xué)技術(shù)、環(huán)保投入等都會(huì)影響EKC曲線的變化。統(tǒng)計(jì)口徑的變化和數(shù)據(jù)的準(zhǔn)確性會(huì)影響污染物排放量變化趨勢(shì)的預(yù)測(cè)。

曲線顯示山東省的某些污染物排放量已經(jīng)越過拐點(diǎn)，但如果未來不繼續(xù)加強(qiáng)環(huán)境保護(hù)仍然有反彈的可能性。出現(xiàn)拐點(diǎn)說明環(huán)境污染物排放量與經(jīng)濟(jì)發(fā)展的關(guān)系發(fā)生變化，但并不意味著環(huán)境質(zhì)量變好或變壞。

（3）建設(shè)可持續(xù)發(fā)展的經(jīng)濟(jì)強(qiáng)省的必然要求，就是要協(xié)調(diào)好環(huán)境質(zhì)量和經(jīng)濟(jì)發(fā)展的關(guān)系。在高質(zhì)量發(fā)展經(jīng)濟(jì)的同時(shí)，進(jìn)一步調(diào)整產(chǎn)業(yè)布局、調(diào)整產(chǎn)業(yè)結(jié)構(gòu)、完善環(huán)保法規(guī)、加大環(huán)境治理投資和加強(qiáng)環(huán)境保護(hù)意識(shí)等，從而減少污染物的排放。

參考文獻(xiàn)：

[1]彭水軍，包群.經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染—環(huán)境庫茲涅茨曲線假說的中國檢驗(yàn)[J].財(cái)經(jīng)問題研究，2006（8）：3-17.

[2]周澤輝，趙娜.中國環(huán)境污染與經(jīng)濟(jì)發(fā)展水平的關(guān)系研究[J].中國城市經(jīng)濟(jì)，2010（6）：88-90.

[3]楊珂玲，葛翔宇，董利民.洱海流域社會(huì)經(jīng)濟(jì)發(fā)展速度與污染物排放量關(guān)系研究[J].生態(tài)經(jīng)濟(jì)，2011（12）：37-40.

[4]虞依娜，陳麗麗.中國環(huán)境庫茲涅茨曲線研究進(jìn)展[J].生態(tài)環(huán)境學(xué)報(bào)，2012，21（12）：2018-2023.

[5]張曉.中國環(huán)境政策的總體評(píng)價(jià)[J].中國社會(huì)科學(xué)，1999（3）：88-99.

[6]吳玉萍，董縮成，宋健峰.北京市經(jīng)濟(jì)增長(zhǎng)與環(huán)境污染水平計(jì)量模型研究[J].地理研究，2002，21（3）：239-246.

[7]段曉峰，許學(xué)工.山東省污染物排放與經(jīng)濟(jì)發(fā)展水平的關(guān)系[J].地理科學(xué)進(jìn)展，2010，29（3）：342-346.

[8]邢秀鳳，劉穎宇.山東省經(jīng)濟(jì)發(fā)展與環(huán)境保護(hù)關(guān)系的計(jì)量分析[J].中國人口.資源與環(huán)境，2006，（1）：58-61.