劉金玉

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

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

文章編號：1004-7026（2019）17-0116-02? ? ? ? ?中國圖書分類號：F019.3 ? ? ? ? 文獻標志碼：A

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

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

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

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

1? 研究區(qū)概況

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

2? 研究方法

2.1? 多元線性回歸分析

運用多元線性回歸分析，從總?cè)丝跀?shù)、地區(qū)生產(chǎn)總值和城鎮(zhèn)化水平3個方面，進行社會經(jīng)濟發(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)濟發(fā)展與環(huán)境污染程度呈倒U型關(guān)系。對山東省環(huán)境污染物排放量與經(jīng)濟發(fā)展程度間的關(guān)系進行EKC曲線模擬，選擇人均GDP為經(jīng)濟增長指標，二氧化硫、煙粉塵、工業(yè)廢水、工業(yè)固體廢物的排放量為環(huán)境污染指標，分別使用對數(shù)函數(shù)、二次函數(shù)、三次函數(shù)進行曲線回歸模擬，根據(jù)擬合優(yōu)度和參數(shù)檢驗的顯著性結(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對二氧化硫排放量的影響是顯著的。地區(qū)生產(chǎn)總值與二氧化硫排放量呈負相關(guān)，說明隨著經(jīng)濟的發(fā)展，沒有以犧牲環(huán)境換取經(jīng)濟增長，反而加強環(huán)境治理，使生態(tài)環(huán)境質(zhì)量持續(xù)好轉(zhuǎn)。②煙粉塵排放量與城市化率呈負相關(guān)，說明隨著山東省城市化進程的加快，煙粉塵排放量減少，環(huán)境質(zhì)量變好。③工業(yè)廢水排放量與社會經(jīng)濟指標的各個變量都具有顯著得相關(guān)性。④工業(yè)固體廢物排放量與地區(qū)生產(chǎn)總值、城市化率呈正相關(guān)，說明經(jīng)濟發(fā)展導致固體廢物的增加，未來要發(fā)展循環(huán)經(jīng)濟，減少“三廢”排放。

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

3.2.1? 二氧化硫與人均GDP的EKC檢驗

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

3.2.2? 煙粉塵與人均GDP的EKC檢驗

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

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

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

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

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

4? 結(jié)論

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

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

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

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

參考文獻：

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

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

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

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

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

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

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

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