趙俊粟 金圣博（對外經(jīng)濟(jì)貿(mào)易大學(xué)保險(xiǎn)學(xué)院，北京 100029）

1 引言

1.1 選定研究對象

作為空氣污染的主要污染物，大氣中的二氧化硫污染，不僅會帶來酸雨等環(huán)境問題，如果人體過量攝入，還可能引發(fā)過敏反應(yīng)，出現(xiàn)呼吸困難、嘔吐等癥狀。因此，做好二氧化硫污染的防治非常重要。本文針對我國大氣中二氧化硫污染的主要影響因素進(jìn)行了實(shí)證分析，從技術(shù)、政策等方面，提出了相應(yīng)的解決對策，希望能夠?yàn)槎趸蛭廴镜姆乐翁峁┮恍﹨⒖肌?/p>

1.2 確定被解釋變量與解釋變量

本文引入第二工業(yè)產(chǎn)值，道路面積，城市建筑用地面積，機(jī)動車擁有量四個變量.其中第二工業(yè)產(chǎn)值作為目標(biāo)變量，剩下三個作為控制變量，分析二氧化硫排放與這些結(jié)構(gòu)性因素的關(guān)聯(lián)。從各地的數(shù)據(jù)情況來看，相關(guān)性如下：

一是第二工業(yè)產(chǎn)值對二氧化硫排放的影響是顯而易見的，作為工業(yè)生產(chǎn)的主要廢料之一，工業(yè)產(chǎn)值高必然伴隨著二氧化硫排放增加。二是機(jī)動車擁有量和道路面積顯然對二氧化硫排放有影響，而這兩個變量本身有一定的相關(guān)關(guān)系，所以本文會引入一個交互項(xiàng)。三是城市建筑用地面積，城市建筑用地面積大，在一定程度上也會加大工廠的用地面積，二氧化硫排放也會有所增加。

2 數(shù)據(jù)的收集和整理

本文數(shù)據(jù)來源于萬得數(shù)據(jù)庫，以及中國統(tǒng)計(jì)年鑒1994——2017年的地級市面板數(shù)據(jù)。

其中，共有287個地級市數(shù)據(jù)，年份跨度23年，本文做面板回歸時(shí)選取了其中十三年的數(shù)據(jù)，本文引入第二工業(yè)產(chǎn)值，道路面積，城市建筑用地面積，機(jī)動車擁有量四個變量.其中第二工業(yè)產(chǎn)值作為解釋變量，剩下三個作為控制變量。

3 截面數(shù)據(jù)回歸結(jié)果分析和檢驗(yàn)

采用的回歸方程是y=β0+β1lnGDP+β2CAR+β3CONST+β4ROAD+β5(ROAD*CAR)

根據(jù)截面數(shù)據(jù)的回歸情況，發(fā)現(xiàn)解釋變量是99%顯著的，說明第二工業(yè)產(chǎn)值是會顯著影響二氧化硫排放的。同時(shí)根據(jù)回歸系數(shù)，對回歸方程的解釋為：當(dāng)?shù)诙I(yè)產(chǎn)值增加1%元，二氧化硫排放平均增加0.16847千噸。

根據(jù)截面數(shù)據(jù)的回歸情況，發(fā)現(xiàn)回歸結(jié)果的均方誤差大，但由于截面數(shù)據(jù)本身的特點(diǎn)，將該數(shù)據(jù)納入可以接受的范圍內(nèi)。

（1）根據(jù)散點(diǎn)圖的情況，決定采用非線性模型。

原因如下：第一，采用對數(shù)的形式可以很好的擬合散點(diǎn)圖的形狀。

第二，由于數(shù)據(jù)的本身的特性，可以得知，一個城市開始發(fā)展的時(shí)候必然是先發(fā)展工業(yè)，二氧化硫的排放主要來源于化石燃料的燃燒，而發(fā)展工業(yè)的首要步驟就是獲取能源，所以當(dāng)?shù)诙a(chǎn)業(yè)產(chǎn)值由0開始增長的時(shí)候二氧化硫一定會迅速增加，當(dāng)其他工業(yè)企業(yè)形成完整體系的時(shí)候，二氧化硫增加量會逐漸變緩慢，非常符合對數(shù)函數(shù)的圖像。

（2）根據(jù)回歸方程，幾個解釋變量與控制變量之間不存在明顯的可以相互表示的情況，根據(jù)回歸結(jié)果得到了進(jìn)一步證實(shí)。

4 面板回歸分析

（1）采用了加入虛擬變量的方式對個體固定效應(yīng)進(jìn)行了回歸，回歸方程為：

回歸結(jié)果發(fā)現(xiàn)R2顯著接近于1，而且第二工業(yè)產(chǎn)值也在1%的顯著性水平下顯著，說明在研究第二產(chǎn)業(yè)產(chǎn)值對二氧化硫排放量的影響過程中個體差異比較明顯，而時(shí)間差異不太明顯。根據(jù)常理也可以推斷，每個城市都有不一樣的城市規(guī)劃，城市功能劃分不盡相同，但對于城市化進(jìn)程來說所有城市基本都遵循工業(yè)化帶動第三產(chǎn)業(yè)發(fā)展的模式，所以在探究這個問題的過程中，個體差異要顯著大于時(shí)間差異。

（2）然后采用了stata 自帶的面板數(shù)據(jù)回歸，得到的結(jié)果和個體固定效應(yīng)回歸差別不大，可以看到，第二產(chǎn)業(yè)產(chǎn)值仍然在1%的顯著性水平下顯著，并且第二工業(yè)產(chǎn)值每增加1%，二氧化硫的排放平均會增加5.563千噸。

5 模型選擇與評價(jià)

5.1 模型評價(jià)

5.1.1 遺漏變量偏差

研究目標(biāo)是為了探尋城市第二工業(yè)產(chǎn)值對該城市二氧化硫排放的關(guān)系，但是進(jìn)行回歸研究時(shí)，總會存在很多遺漏變量。

為了解決遺漏變量偏差的問題，本文搜集到了道路面積，城市建筑用地面積，機(jī)動車擁有量的數(shù)據(jù)，將它們作為控制變量加入回歸模型，盡可能地充分控制遺漏變量的影響。

5.1.2 回歸函數(shù)形式的誤設(shè)

觀察散點(diǎn)圖發(fā)現(xiàn)，被解釋變量二氧化硫與目標(biāo)變量第二工業(yè)產(chǎn)值的關(guān)系，很可能不是線性的。觀察回歸結(jié)果，發(fā)現(xiàn)采用對數(shù)形式的回歸結(jié)果的修正R2更大，回歸效果更好。

5.1.3 雙向因果關(guān)系

對于本文研究的問題，雙向因果效應(yīng)不是特別明顯。

工業(yè)產(chǎn)值增加會導(dǎo)致二氧化硫排放量增加，工業(yè)排放二氧化硫占總二氧化硫污染的百分之九十五以上，而我國主要的能源即化石燃料的燃燒，同時(shí)我國的高能源消耗型工業(yè)生產(chǎn)，工業(yè)產(chǎn)值增加一定會導(dǎo)致能源需求的增加，二氧化硫排放必然會提升。

但二氧化硫排放量上升并不一定導(dǎo)致工業(yè)產(chǎn)值增加，有可能是該地區(qū)政策原因沒有限制企業(yè)的污染排放，企業(yè)沒有合適的廢料回收機(jī)制，某些無良企業(yè)使用劣質(zhì)原材料導(dǎo)致等等，即解釋變量與被解釋變量之間具有單項(xiàng)因果關(guān)系。

5.2 最優(yōu)模型與結(jié)論

最終的回歸結(jié)果的模型形式

該模型表明：

截距項(xiàng)：199.7296，

5.56：當(dāng)?shù)诙I(yè)產(chǎn)值增加1%，二氧化硫排放平均增加0.0556千噸

-0.026：汽車擁有量增加一千輛，二氧化硫排放平均減少26噸

-0.003：道路面積增加一百萬平方米，二氧化硫排放平均減少3噸，

-0.7004：建筑面積增加十萬平方米，二氧化硫排放平均減少700.4噸。

（后三個變量，其回歸系數(shù)在一定程度上違反常理，但它們作為控制變量，OLS 估計(jì)得到的系數(shù)是不準(zhǔn)確的，因此經(jīng)濟(jì)意義可能不夠有效）

研究結(jié)論：

第二工業(yè)產(chǎn)值對二氧化硫的排放量有顯著的影響，如果想要控制二氧化硫的排放量，從根本上來說還是應(yīng)該控制工業(yè)生產(chǎn)。

可以從三個方面降低工業(yè)產(chǎn)值，進(jìn)而減少二氧化硫的排放量，打到污染治理的目的。

其一，特大城市進(jìn)行產(chǎn)業(yè)升級和轉(zhuǎn)型，北京已經(jīng)將大多數(shù)工業(yè)企業(yè)遷至郊區(qū)或者周邊省市，這是減少北京工業(yè)產(chǎn)值的一大關(guān)鍵舉措。但不意味著僅僅遷移工廠，在遷移的同時(shí)，北京這樣的大城市應(yīng)該注重第三產(chǎn)業(yè)的發(fā)展，補(bǔ)充第二產(chǎn)業(yè)造成的gdp缺口。

其二，工業(yè)產(chǎn)業(yè)為支柱的省市著重推行能源高效利用，開啟清潔能源試點(diǎn)，盡量從源頭減少二氧化硫的產(chǎn)生，這樣不會減少太多工業(yè)產(chǎn)值，同時(shí)能達(dá)到減少污染的目的。

其三，重污染生產(chǎn)外包，類似稀土開采加工這樣的高污染行業(yè)完全可以轉(zhuǎn)包給其他國家，將污染的負(fù)外部性轉(zhuǎn)嫁給其他國家，即便付出一些成本，也是值得的。