1. 引言
改革开放之后,为了实现经济的快速发展,我国的能源消耗量也迅速加大 [1] ,近年来由于能源过度使用以及低效率使用所带来的环境污染问题、能源紧张问题已经逐渐显现了出来。现在所推行的新能源汽车、绿色照明工程、建筑节能工程等都是贯彻国家节能减排政策、优化能源结构、既要金山银山又要绿水青山的思想的一系列有效的措施。GDP每万元碳排放量等指标也成为国家关注发展的质量的一个重要指标,但是仅仅通过这样一个比例关系并不能完全的反应二者之间存在关系,代表性较弱。因此深入探索能源消耗与经济发展的关系,可以为国家节能减排政策的制定提供科学的建议,同时也可以对上述政策的实施效果做一个侧面的评价,从而为我国经济可持续发展、经济结构的有效调整提供内在动力。
2. 数据来源于符号说明
本文从山东省统计局公布的经济发展与能源数据中选取了以下年度指标进行研究 [2] (见表1)。
对本文所出现的符号做统一的规定(见表2)。
3. 研究方法
3.1. 主成分分析
主成分分析法的本质是降维,将原来多个复杂的变量进行最佳综合,简化为少数几个相互独立的综合指标,这几个综合指标既可以包含原来多个复杂变量的绝大部分信息,又去除了原变量信息间的重叠,具有较强的客观性。
主成分是原始变量的一些线性组合,从而将大量原始变量的信息提取出来。主成分的方差,即Var(F)反应了该主成分所包含的信息量,显然应该选取方差最大的主成分做成分作为第一主成分 [3] 。通常情况下如果数据相关性较强选择第一主成分就可以对数据有一个较好的代表性,当然也可以选择第二、三、四等主成分加强对原始数据的解释。但是为了有效的反映原始数据,要保证主成分两两正交,即即F1已有的信息就不需要再出现在F2中,在确定第一主成分后找到满足
的第二个线性组合F2,就是第二主成分。以此类推可以得到以下公式:
代表前p个主成分,
是原始变量矩阵标准化后的数据,P < n。
(
)是标准化后数据的协方差阵的最大特征值对应的特征向量。
3.2. 协整性检验
3.2.1. 单整
对于非平稳的时间序列数据,通过运用差分的方式进行平稳化处理,一阶单整指的是经过一次差分处理后非平稳时间序列便趋于平稳,记为I(1)。一般而言,如果时间序列平稳,记为I(0),如果d阶差分后平稳,称为d阶单整,记为I(d)。
3.2.2. 协整
如果时间序列
,都是d阶单整,即I(d),存在一个向量
使得
,这里
则称序列
是(d, b)阶协整,记为
,
为协整向量。
协整在探索经济社会发展规律方面是具有显著意义的,经济社会发展是一个具有长期规律性波动的过程。但是其发展情况及其复杂,通常受到外界各种因素的影响,很有可能存在一个或多个与它保持长期稳定均衡关系的指标。这种协整关系的存在便可以通过其它变量的变化来影响另一变量的变化 [4] 。
3.2.3. 协整性检验(E-G两步法)
两个或多个时间序列同阶单整时才可能协整,否则直接判断为非协整。若同阶单整,用以下两种方法检验是否协整。
E-G两步法比较适合于两个变量之间协整关系的检验。xt和yt均为d阶单整的时间序列数据。
第一步,做普通最小二乘回归。
得到模型的残差
第二步,检验残差是否平稳,如果平稳,则认为时间序列xt和yt是协整的。
3.3. Bootstrap方法
Bootstrap称为自助法,是一种有放回的重复抽样。从容量为n的样本中进行有放回的抽样,并不改变样本容量 [5] [6] 。
Bootstrap优势:一是可以减小估计量的偏差,二是对很多估计能够获得精准的置信区间,三是在某些情况下极限分布难以获得,bootstrap可以用做统计推断,四是存在严重的离群点或样本容量不够大时,自助法非常实用。
4. 实证分析
4.1. 主成分分析
本文采用R语言进行山东省能源消耗与经济发展关系的实证分析,首先利用主成分分析方法,分别对反映经济发展的指标和反映能源消耗的指标进行降维,确定两个综合指标得分。
首先提取能源消耗的主成分,为了消除数据量纲的差异,首先对对数据进行了标准化处理。根据指标间的相关系数矩阵(见图1)可以发现能源消耗变量间存在很强的相关性,适用于进行主成分降维。
由主成分分析结果(见图2)可以看出第一主成分的方差贡献率到达了84.9%,第二主成分的方差贡献率为12.32%,前两个主成分的累计方差贡献率高达97.23%,损失信息非常少,所以选择前两个主成分便可以代表所有变量的信息,通过碎石图可以更加直观的说明主成分的选择(见图3)。
同理经济发展变量中我们提取了第一主成分的方差贡献率即可达到97.72%,即选取第一主成分就可以代表97.72%的信息
计算1995年到2015年各年份在第一主成分上的得分情况。下面展示了前5年的结果(见图4、图5):
![](//html.hanspub.org/file/1-3050046x23_hanspub.png)
Figure 1. Correlation coefficient matrix for energy consumption
图1. 能源消耗相关系数矩阵
![](//html.hanspub.org/file/1-3050046x24_hanspub.png)
Figure 2. Principal component analysis of energy consumption
图2. 能源消耗主成分分析
![](//html.hanspub.org/file/1-3050046x26_hanspub.png)
Figure 4. Principal component scores of energy consumption indicators
图4. 能源消耗主成分得分
![](//html.hanspub.org/file/1-3050046x27_hanspub.png)
Figure 5. Principal component scores of economic indicators
图5. 经济指标主成分得分
最后根据综合评价函数
计算主成分的综合得分如下(见表3):
由综合得分时序图(见图6)可以看出,经济发展的主成分得分是逐年提高的,这表明山东省的经济发展一直保持逐年增长的趋势,而能源消耗的主成分得分却出现了一些小的波动,特别是在2011年有了一个明显的下降。一方面,受利比亚战乱、埃及局势动荡以及欧佩克组织限产等因素的影响,世界石油价格高涨。另一方面,全省节能降耗政策取得新成效,如351户企业完成落后产能淘汰。因此,燃料油、汽油、柴油和液化石油气的消耗量较2010年有较大下降,能源消耗的主成分得分在2011年显著下降。
![](Images/Table_Tmp.jpg)
Table 3. Scores of comprehensive energy and economic indicators
表3. 能源综合指标与经济综合指标得分
以1997年、2004年、2011年为三个分割点,可以看出山东省经济发展受能源消耗的影响的发展趋势基本是相同的,但是其影响力却逐渐发生着变化。
4.2. 协整分析
由上图可以很清晰的得出结论,两个综合得分时间序列是不平稳的。分别作一阶差分后的时序图(见图7)如下:
两项指标单位根检验的P值分别为0.0252和0.0232均小于0.05 (见图8),说明经过一阶差分序列已经平稳。进而可以利用两变量的Engle-Granger检验进行协整性检验,首先建立回归方程并用最小二乘法估计参数,得到非均衡误差。
经过对误差项的单位根检验,可以确定误差项平稳,说明能源消耗与经济发展存在长期均衡关系,其回归方程为
。
![](//html.hanspub.org/file/1-3050046x31_hanspub.png)
Figure 7. Composite Score After Difference
图7. 差分后综合指标得分
但是,这里存在一个问题,模型的DW检验未能通过(p-value = 3.288e^−7),说明模型中有严重的自相关问题
所以考虑把因变量的滞后一期和滞后二期的结果引入回归模型中重新建模(见图9)。
所以经过修正后的模型变为Nengyuanscore = 0.3*jingjiscore,即经济得分每提高提高一分,能源得分提高0.3分。
4.3. Bootstrap方法与最小二乘法的比较
模型诊断图(见图10)可以看出,误差项的正态性假定并不太满足,样本量也不大,根据误差正态性假定得出的回归系数的置信区间不具有说服力。
利用bootstrap方法,重复抽样2000次,对回归方程的系数进行估计。步骤如下:
1) 有样本容量为n的样本
,从中有放回的重复抽样n次,得到了第一个样本容量为n的新样本,称为自助样本1,对于原始数据对
,在自助样本中可能出现多次,可能出现零次。
2) 用自助样本1进行最小二乘法回归分析,估计其回归参数,不妨记为t1
3) 重复上述步骤m次,得到
。取置信水平为 95%,那么把
从小到大排序,置信区间为这列数的(上0.975分位数,上0.025分位数)。
通过比较可以发现,采用Bootstrap所得到的回归方程,其参数估计的置信区间要更具说服力(见表4)。因此,采用bootstrap法建立能源消耗与经济发展的回归模型可以更好的反映二者的关系,从而有助于根据能源消耗量测算其所带来的经济增长量。经济得分每提高一分,能源消耗得分提高0.03分到0.5分。
5. 结论建议
由能源消耗与经济发展的相互拉动关系可以看出,能源消耗对于经济的贡献与在我国不同的发展时
![](//html.hanspub.org/file/1-3050046x37_hanspub.png)
Figure 9. Regression results with the introduction of hysteresis
图9. 引入滞后效应的回归结果
![](Images/Table_Tmp.jpg)
Table 4. Comparison of Least Squares with Bootstrap
表4. 最小二乘与Bootstrap比较
期具有明显的差异,在2005年到2010年,能源消耗的贡献率明显低于2010年以后,这说明我国近年来的能源结构调整和经济发展战略都取得了显著的成效,必须将这种改革深入贯彻落实下去。同时必须采取积极有效的措施解决现有经济-能源-环境关系,为此可以采取以下措施:
1) 在大力倡导使用天然气等清洁能源的同时着力解决其使用过程中存在的问题,例如天然气汽车稳定性、加气便利性的问题都亟待解决,只有让消费者感受到清洁能源要更加的便捷经济,才能够从根本上刺激消费转型,从而在不影响经济发展的基础上减少能源消耗带来的环境问题。
2) 继续推进产业结构升级,寻找新的经济增长的拉动因素,降低能源消耗特别是不可再生能源的消耗对于经济发展的贡献率,进而达到提高单位能源消耗的经济转化率。
3) 建立能源消耗对环境影响的评价指标体系,直接将能源消耗与环境问题挂钩,不仅可以提高相关部门、企业的重视程度,也可以提供一个更直观的评价标准。
致谢
感谢在本文完成过程中家人和老师对我的帮助,同时也感谢所有参考或引用的文献作者。