（一）指標(biāo)說(shuō)明

第一，碳排放指標(biāo)。能源部門通常是溫室氣體排放中的最重要部門，在發(fā)達(dá)國(guó)家，其貢獻(xiàn)一般占CO2排放量的90%以上和溫室氣體總排放量的75%。根據(jù)《中華人民共和國(guó)氣候變化初始國(guó)家信息通報(bào)》中的數(shù)據(jù)，能源活動(dòng)是中國(guó)最主要的CO2排放源，1994年中國(guó)能源活動(dòng)的CO2排放量在全國(guó)CO2排放總量中占90.95%。因此，本文主要研究中國(guó)與能源活動(dòng)有關(guān)的CO2排放，包括煤炭、石油、天然氣等三種化石燃料消耗。

中國(guó)碳排放量采用以下公式進(jìn)行估算：

（1）

其中，C為碳排放量；E為中國(guó)一次能源的消費(fèi)總量；Fi為i類能源的碳排放系數(shù)；Si為i類能源在總能源中所占的比重。其中一次能源消費(fèi)總量E和比重Si的計(jì)算均采用發(fā)電煤耗計(jì)算法，數(shù)據(jù)取自于《中國(guó)能源統(tǒng)計(jì)年鑒1989》和《中國(guó)能源統(tǒng)計(jì)年鑒2009》。通過(guò)查閱有關(guān)文獻(xiàn)，收集有關(guān)能源消耗的碳排放系數(shù)并進(jìn)行比較計(jì)算，最終取平均值確定為各能源消耗碳排放系數(shù)Fi。

（二）數(shù)據(jù)來(lái)源及處理

本文以1978年至2008年我國(guó)重工業(yè)產(chǎn)值和碳排放量為分析的基礎(chǔ)數(shù)據(jù)，重工業(yè)產(chǎn)值和碳排放量分別用ZGY和TPF表示。由于數(shù)據(jù)的自然對(duì)數(shù)變換不改變?cè)瓉?lái)的協(xié)整關(guān)系，并能使其趨勢(shì)線性化，消除時(shí)間序列中存在的異方差現(xiàn)象，對(duì)兩個(gè)變量取自然對(duì)數(shù)，分別表示為L(zhǎng)ZGY和LTPF。圖1是1978年到2008年間我國(guó)重工業(yè)產(chǎn)值和碳排放量變化趨勢(shì)。從圖中可以看出，考察期內(nèi)，我國(guó)重工業(yè)產(chǎn)值和碳排放量的走勢(shì)基本一致，尤其是2000年以后，兩者的變化趨勢(shì)幾乎完全一致，都經(jīng)歷了一個(gè)上升的過(guò)程。這說(shuō)明，在我國(guó)重工業(yè)發(fā)展和碳排放量之間存在著長(zhǎng)期對(duì)應(yīng)關(guān)系。本文通過(guò)運(yùn)用時(shí)間序列動(dòng)態(tài)均衡的分析方法，在單位根檢驗(yàn)、協(xié)整檢驗(yàn)的基礎(chǔ)上建立了誤差修正模型，以此說(shuō)明我國(guó)產(chǎn)業(yè)結(jié)構(gòu)重型化與碳排放的長(zhǎng)期均衡關(guān)系和短期動(dòng)態(tài)關(guān)系。

（一）單位根檢驗(yàn)

在實(shí)際經(jīng)濟(jì)運(yùn)行中，經(jīng)濟(jì)變量的時(shí)間序列數(shù)據(jù)很少是平穩(wěn)的，若采用傳統(tǒng)普通最小二乘法進(jìn)行回歸分析，很容易導(dǎo)致“謬誤回歸”。從圖1我們也可以看出，1978-2008年間，我國(guó)重工業(yè)產(chǎn)值增長(zhǎng)和碳排放總量均具有明顯的上升趨勢(shì)，顯示出不平穩(wěn)性。因此，為了保證回歸的真實(shí)性，本文使用ADF法對(duì)變量LZGYt和LTPFt進(jìn)行平穩(wěn)性檢驗(yàn)，檢驗(yàn)形式由序列圖確定，滯后階數(shù)根據(jù)SC最小準(zhǔn)則確定，具體檢驗(yàn)結(jié)果如表2所示。

表2的結(jié)果顯示，對(duì)數(shù)序列LZGY、LTPF的ADF值均大于1%的顯著性水平的臨界值，不能拒絕存在單位根的假設(shè)，即序列都是非平穩(wěn)的；但其一階差分在1%的置信性水平下拒絕了原假設(shè)，說(shuō)明是平穩(wěn)的。因此可以判定LZGY、LTPF都是一階單整的，它們都是I（1）序列，滿足變量協(xié)整的條件，即LZGY和LTPF很可能存在協(xié)整關(guān)系。

（二）協(xié)整關(guān)系檢驗(yàn)

雖然一些經(jīng)濟(jì)變量的本身是非平穩(wěn)序列，但是它們的線性組合卻有可能是平穩(wěn)序列，這種平穩(wěn)的線性組合被稱為協(xié)整方程且可被解釋為變量之間的長(zhǎng)期穩(wěn)定的均衡關(guān)系。目前最常用的協(xié)整分析方法主要有Engle-Granger（EG）的兩步法以及Johansen和Juselius（JJ）的極大似然法。前者適合于雙變量的協(xié)整檢驗(yàn)，本文采用此方法進(jìn)行檢驗(yàn)。

首先用LZGY對(duì)LTPF進(jìn)行回歸，即建立如下計(jì)量方程：

（2）

采用OLS法對(duì)式（2）進(jìn)行估計(jì)，得：

（3）

R2=0.977586；Adj-R2=0.976813；F=1264.835***；D-W=0.188229

式中，括號(hào)內(nèi)數(shù)字為t統(tǒng)計(jì)量值，***表示在0.01水平上顯著，et為回歸殘差。

式（3）的D-W值指示誤差項(xiàng)存在正自相關(guān)，因此，參數(shù)估計(jì)值顯著性檢驗(yàn)失效，從而預(yù)測(cè)值盡管是無(wú)偏的，但卻不是有效的。為了消除誤差項(xiàng)的序列相關(guān)性，我們建立如下（p，q）階回歸分布滯后模型：

(4）

其中，p、q分別為L(zhǎng)TPF和LZGY的最大滯后階數(shù)。依據(jù)英國(guó)計(jì)量經(jīng)濟(jì)學(xué)家亨德瑞等人倡導(dǎo)的“從一般到簡(jiǎn)單”的Hendary/LSE建模方法，考慮到采用的是年度數(shù)據(jù)以及自由度的約束，我們選擇最大滯后階數(shù)p、q均為2，建立一般模型，然后利用Wald檢驗(yàn)和t檢驗(yàn)，依次剔出一個(gè)或多個(gè)系數(shù)不顯著的變量，進(jìn)行模型簡(jiǎn)化，并結(jié)合模型選擇統(tǒng)計(jì)量（AIC和SC）選擇一個(gè)最優(yōu)的模型作為最終模型，估計(jì)結(jié)果如下：

（5）

R2=0.996901；Adj-R2=0.996384；F=1929.861***；D-W=1.441374

式中，括號(hào)內(nèi)數(shù)值為t統(tǒng)計(jì)量值，**、***分別表示在0.01、0.05水平上顯著。

接下來(lái)需要對(duì)（5）式的殘差項(xiàng)進(jìn)行序列相關(guān)性檢驗(yàn)。由于方程的解釋變量存在被解釋變量的滯后項(xiàng)，此時(shí)D-W值就不能作為判斷回歸方程的殘差是否存在序列相關(guān)的標(biāo)準(zhǔn)，如果殘差序列存在序列相關(guān)，那么，顯著性水平、擬合優(yōu)度，F(xiàn)統(tǒng)計(jì)量將不再可信。這種情況下可以利用LM統(tǒng)計(jì)量進(jìn)行檢驗(yàn)。采用LM統(tǒng)計(jì)量進(jìn)行檢驗(yàn)（p=2），得到結(jié)果如下：

LM統(tǒng)計(jì)量顯示，F(xiàn)統(tǒng)計(jì)量和T×R2統(tǒng)計(jì)量均小于在10%的顯著性水平下的臨界值，因此無(wú)法拒絕原假設(shè)“直到p階后不存在序列相關(guān)”，即可以認(rèn)為回歸方程的殘差序列不存在序列相關(guān)性，因此用方程（5）來(lái)分析產(chǎn)業(yè)結(jié)構(gòu)重型化對(duì)碳排放的影響。

由式（5）計(jì)算得誤差項(xiàng)：

（6）

對(duì)該誤差項(xiàng)進(jìn)行平穩(wěn)性檢驗(yàn)，如果et是平穩(wěn)的，則說(shuō)明LTPF和LZGY之間具有協(xié)整關(guān)系。仍采用ADF單位根檢驗(yàn)方法，選擇既無(wú)趨勢(shì)項(xiàng)也無(wú)常數(shù)項(xiàng)，檢驗(yàn)時(shí)采用SC最小準(zhǔn)則自動(dòng)選擇滯后階數(shù)為5。最終檢驗(yàn)結(jié)果見表4。

殘差序列et的ADF檢驗(yàn)結(jié)果表明，ADF檢驗(yàn)統(tǒng)計(jì)量小于1%顯著水平下的臨界值，說(shuō)明殘差et為平穩(wěn)序列。由此可以得出結(jié)論，我國(guó)產(chǎn)業(yè)結(jié)構(gòu)重型化和碳排放之間存在協(xié)整關(guān)系。而為了得到二者之間的長(zhǎng)期均衡關(guān)系，只需令：LTPFt=LTPFt-1=LTPFt-2=LTPF*，LZGYt=LZGYt-1=LZGY*，并代入（5）式，整理后得到：

（6）

（6）式即為產(chǎn)業(yè)結(jié)構(gòu)重型化與碳排放之間的長(zhǎng)期均衡關(guān)系，它表明，重工業(yè)產(chǎn)值每提高1%，碳排放量就會(huì)提高0.278141%。這表明產(chǎn)業(yè)結(jié)構(gòu)重型化加大了我國(guó)的碳排放量。