有效利用低碳技术是应对气候变化的重要途径。以2006~2015年中国30个省级面板数据为研究范围,采用空间杜宾模型,实证研究低碳创新及其组合对碳排放的作用机制。研究显示,低碳创新与碳排放在空间分布上存在一定偏差;更具突破性的低碳创新作用高于增量创新,但两类创新的组合作用最显著;城市化、产业结构以及FDI等引致碳排放增加,而环境规制对碳排放的抑制作用仍限于本地。这就要求既要协同不同类型创新活动,又要协同不同区域低碳治理。 The effective use of low-carbon technologies is an important way to combat climate change. The paper adopts spatial Durbin model to empirically study the mechanism of low-carbon innovation and its portfolio on carbon emissions with 30 provincial panel data of China from 2006 to 2015. The research shows that the spatial distribution of low-carbon innovation and carbon emissions has a certain deviation. And the breakthrough low-carbon innovation plays a bigger role than the incremental innovation, but the combination of the two types of innovation plays the most signifi-cant role; Urbanization, industry structure and FDI lead to increased carbon emissions, while the impact of environmental regulations on carbon emissions is still limited to local areas. This requires the collaboration of different types of innovation activities and the collaboration of different regional low-carbon governance.
王淼,王为东,卢娜,张财经
江苏大学财经学院,江苏 镇江
收稿日期:2019年2月3日;录用日期:2019年2月19日;发布日期:2019年2月26日
有效利用低碳技术是应对气候变化的重要途径。以2006~2015年中国30个省级面板数据为研究范围,采用空间杜宾模型,实证研究低碳创新及其组合对碳排放的作用机制。研究显示,低碳创新与碳排放在空间分布上存在一定偏差;更具突破性的低碳创新作用高于增量创新,但两类创新的组合作用最显著;城市化、产业结构以及FDI等引致碳排放增加,而环境规制对碳排放的抑制作用仍限于本地。这就要求既要协同不同类型创新活动,又要协同不同区域低碳治理。
关键词 :碳排放,突破式,增量式,低碳创新组合,空间杜宾
Copyright © 2019 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
气候变化正在逐渐危及到人类赖以生存的地球环境,是涉及地球各层级及人类社会各方面的综合性环境问题 [
专利是衡量创新能力的重要指标 [
关于气候变化影响因素的识别及其作用研究一直以来深受学术界的关注,不同的学者从不同的方面,采用不同的方法对影响碳排放的因素进行考察,目前已识别出经济增长、产业结构、城市化率、外商投资以及技术创新等因素 [
但以上研究只是了采用了一般的普通面板回归模型,忽略了地区之间的关联性、区域间的信息和资源的交换性、以及二氧化碳在地区间的流动性,即地区间不存在任何的交换和溢出作用。然而空间计量专家Anselin经过探讨发现空间因素对环境问题有重要影响 [
综上,本文的创新点在于:一是依据已有的较为成熟的份额估计的测度方法将低碳创新活动进行分类,在评测出合理的具有突破式作用和增量式作用的低碳创新技术的同时,考虑其组合作用对碳排放的影响;二是从空间的角度对低碳创新活动对碳排放的作用进行分析,考虑了要素之间的空间溢出性,更完整、准确的揭示影响气候变化的机制。
根据文献综述,碳排放的研究主要围绕着经济发展水平、结构性因素和技术因素对缓解气候变化的促进作用。本文在检验低碳创新活动对碳排放的促进作用时,综合考虑经济增长、产业结构、城市化率、环境规制以及外商直接投资等因素对碳排放的影响。因此,设定的初始计量模型如式(1)所示。
L n c a r b o n i t = α 0 + α 1 Y 02 i t + δ L n X + ε i t (1)
式中,i代表各省市;t代表时间;Lncarbon为被解释变量,代表碳排放量;Y02代表低碳创新专利(包含突破式创新、增量式创新及其组合);X为一组控制变量,包含了可能对碳排放产生影响的若干变量; δ 代表控制变量的弹性系数; ε i t 代表残差项。
由于碳排放和低碳创新活动具有空间溢出型,普通的计量模型可能存在偏差,本文将采用空间计量模型从空间视觉验证低碳创新及其组合对碳排放的作用。本文具体采用更具一般性的空间杜宾模型(SDM),它是空间滞后模型(SLM)和空间误差模型(SEM)的一般形式,表达式如式(2)所示。
y i t = c + ρ ∑ i = 1 n W i t y i t + α x i t + ∑ i = 1 n W i t x i t γ + μ i + δ t + ε i t (2)
式中, y i t 为因变量, x i t 为自变量,α为自变量系数,c为常数项, ρ 为因变量空间自回归系数, γ 为自变量空间滞后系数, μ i 、 δ t 分别表示空间效应和时间效应, ε i t 为残差项,为空间权重矩阵,表示各个空间要素之间的关联性和相互影响程度,本文选取空间邻接矩阵,即相邻地区为1,不相邻地区为0。表示本地区自变量对因变量的影响;表示本地区自变量对邻近地区因变量的影响,即空间溢出效应。
本文预期低碳创新及其组合对气候变化具有正作用,即随着突破式创新专利的不断增加,碳排放量不断减少,同时低碳创新的组合的影响机制显著与单项低碳创新活动。碳排放的核算参考IPCC提供的核算方法,公式如式(3)所示。
c a r b o n i = 12 44 [ ∑ n E n , i × α n × β n ] E n , i (3)
式中,i指不同省市; c a r b o n i 值i省市碳排放,单位为万t; E n , i 指i省市第n种能源终端消费量,焦炉煤气、其他煤气、天然气的单位为m3;其他能源为万t; α n 指第n种能源折标准煤系数; β n 指第n种能源CO2排放系数,单位万t/万t ce。
低碳创新专利(IP)以2013年美国和英国联合颁布的CPC(合作专利分类法)中Y02分类申请专利数来衡量 [
考虑到碳排放还会受经济结构、人口结构、外商需求等其他变量的影响,参照已有成果 [
本文以2006~2013年为研究区间,30个省份为研究对象(因西藏及港澳台数据缺失,不予以考虑)。各省份的Y02专利数据来自incopat专利数据库;能源消费、折标煤系数以及碳排放系数来源于《中国能源统计年鉴》;各省份的人口数据来源于国家统计局;地区产值以及外商直接投资来源于《中国统计年鉴》。为了统一货币单位,采用人民币兑美元年均价汇率将美元为单位的外商直接投资换成人民币。同时为了消除价格波动的影响,以2006年为基期,对经济变量做不变价处理,为了消除异方差,对部分变量进行取对数处理。
空间相关性是指邻近空间分布对象属性值之间的统计相关性,在进行空间计量模型之前需要对核心变量进行空间自相关检验,通常采用Moran’s I指数检验区域变量是否存在自相关。全局空间自相关反映了变量在全部样本区域的空间依赖程度,检验结果如表1所示。可以看出,在邻接权重矩阵下,碳排放与低碳创新专利的全局Moran’s指数均为正,绝大多数通过5%水平检验,说明这三项指标均存在空间正相关性。
碳排放强度 | 突破式低碳创新专利 | 增量式低碳创新专利 | |||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
年份 | Moran’s I | Z值 | p值 | Moran’s I | Z值 | p值 | Moran’s I | Z值 | p值 | ||
2006 | 0.273 | 1.934 | 0.027 | 0.097 | 0.901 | 0.184 | 0.254 | 1.803 | 0.036 | ||
2007 | 0.309 | 2.116 | 0.017 | 0.227 | 1.794 | 0.036 | 0.300 | 2.072 | 0.019 | ||
2008 | 0.292 | 2.012 | 0.022 | 0.009 | 0.294 | 0.385 | 0.265 | 1.879 | 0.030 | ||
2009 | 0.309 | 2.118 | 0.017 | 0.011 | 0.307 | 0.379 | 0.288 | 2.025 | 0.021 | ||
2010 | 0.338 | 2.312 | 0.010 | 0.200 | 1.548 | 0.061 | 0.323 | 2.248 | 0.012 | ||
2011 | 0.268 | 1.903 | 0.029 | 0.250 | 1.858 | 0.032 | 0.328 | 2.281 | 0.011 | ||
2012 | 0.304 | 2.115 | 0.017 | 0.124 | 1.018 | 0.154 | 0.379 | 2.658 | 0.004 | ||
2013 | 0.299 | 2.087 | 0.018 | 0.249 | 1.876 | 0.030 | 0.366 | 2.568 | 0.005 | ||
2014 | 0.243 | 1.749 | 0.040 | 0.227 | 1.721 | 0.043 | 0.264 | 1.873 | 0.031 | ||
2015 | 0.245 | 1.756 | 0.040 | 0.239 | 1.165 | 0.049 | 0.307 | 2.166 | 0.015 | ||
表1. 2006~2013年中国各省市碳排放及低碳创新全局Moran’I表
局部空间自相关采用Moran’s I散点图用来验证各省市碳排放及低碳创新专利与其邻近地区之间的关系。局部自相关共有四种模式,对应分布在四个象限。分别为第一象限的高值区周围是高值区(HH型),表明两者空间差异小;第二象限的低值区周围是高值区(LH型),表明两者空间差异大;第三象限的低值区周围是低值区(LL型),表明两者空间差异小;第四象限的高值区周围是低值区(HL型),表明两者空间差异大。(由于篇幅有限)本文选取碳排放的2006年和2012年的散点图进行展示,结果如图1所示。可以看出碳排放局部空间集聚特征较明显,总体以HH集聚和LL集聚类型为主,并随着时间的变化,各省份的类型呈现动态变化特征。
图1. 碳排放Moran’s I散点图
通过前文的空间相关性检验可知,碳排放和低碳创新活动的Moran’s I均为正的空间自相关性,且绝大多数通过了5%显著性检验,因此在构造模型时,有必要纳入空间因素。首先根据普通面板数据模型(OLS)回归残差的空间自相关检验显示,LMlag、Robust-LM lag和LMerror、Robust-LM error的p统计值均在1%显著性水平下显著,说明构建空间计量模型比较合理;接着对SDM模型进行LR检验,结果显示均通过了1%显著性水平检验,说明SDM模型不可退化为SLM或SEM模型;最后Hausman检验结果,以及空间杜宾模型的空间自相关系数(rho)显示应该选用随机效应形式(具体模型结果见表2)。
SDM模型 | SDM模型分解 | |||||||||
---|---|---|---|---|---|---|---|---|---|---|
固定效应 | 随机效应 | 直接效应 | 间接效应 | 总效应 | ||||||
disr | −0.00265** (−0.0011) | −0.00252** (−0.00126) | −0.00251* (−0.00134) | −0.0005 (−0.00306) | −0.0031 (−0.00368) | |||||
incr | −3.64e-05** (1.43e-05) | −3.42e-05** (1.59e-05) | −3.64e-05** (1.63e-05) | −0.0000 (3.24e-05) | −0.0001 (3.92e-05) | |||||
lnpgdp | −0.0212 (−0.0166) | −0.0290* (−0.0173) | −0.0397** (−0.0179) | −0.188*** (−0.0503) | −0.228*** (−0.0613) | |||||
ur | 3.565*** (−0.526) | 3.272*** −0.5630 | 3.308*** (−0.518) | 0.5960 (−0.606) | 3.904*** (−0.435) | |||||
is | 1.150*** (−0.25) | (1.220***) −0.2530 | 1.283*** (−0.246) | 0.954* (−0.526) | 2.237*** (−0.601) | |||||
lner | 2.2730 (−2.284) | (−1.199) −2.6000 | −1.0720 (−2.689) | 0.0716 (−7.8) | −1.0000 (−9.032) | |||||
lnfdi | 0.0572*** (−0.0216) | −0.0374 (−0.0249) | 0.0597** (−0.0267) | 0.334*** (−0.118) | 0.393*** (−0.129) | |||||
C | 10.73*** (−1.127) | |||||||||
rho | −0.1330 (−0.0852) | 0.351*** (−0.0631) | |||
---|---|---|---|---|---|
Log-likelihood | 229.1803 | 77.6818 | |||
R-sq | 0.2030 | 0.5692 | |||
Hausman | −16.6600 | ||||
LM-lag | 37.556*** | ||||
Robust LM-lag | 36.422*** | ||||
LM-error | 24.25*** | ||||
Robust LM-error | 12.90*** | ||||
LR-lag | 28.64*** | ||||
LR-error | 28.14*** |
表2. 空间Dubin模型(SDM)回归结果及模型分解(1)
注:*、**、***分别表示通过10%、5%、1%的显著性检验;括号内的值为标准误。
空间杜宾模型将变量的滞后因子纳入研究范围,因此解释变量的系数不能直接拿来反映对因变量的影响,但从模型的回归结果中可以看出突破时低碳创新与增量式低碳创新对碳排放估计系数显著为负,表明低碳创新活动对碳排放具有明显的抑制作用,与预期相符。
从直接效应来看:突破式低碳创新(disr)对本地区碳排放的影响系数为−0.00251,通过了10%的显著性检验,增量式低碳创新(incr)对本地区的碳排放的影响系数为−0.0000364,通过了5%的显著性检验,尽管影响力度不是很理想,可能的原因是专利从开始到被发掘、大量被引,发挥作用需要一定的时间 [
从间接效应来看:突破式低碳创新与增量式低碳创新的空间溢出效应为负,说明一个地区的低碳创新活动对其他地区的碳排放也具有抑制作用,但尚不显著,表明目前一个地区的低碳创新活动还主要局限于本地,区域的联动性与技术发展水平不高,共同推进低碳创新的机制尚未形成。
从控制变量的直接效用和间接效应来看:1) 经济增长(lnpgdp)。人均GDP对无论本地区还是临近地区的碳排放都具有显著的反向作用,即经济越发达,碳排放反而越少。说明随着人们对生态问题的不断关注,经济越是发达的地区越是加快调整产业结构、改变发展模型。2) 城市化率(ur)。城市化率对本地区的碳排放具有显著的正向作用,即城市化水平的提高会促进碳排放。出现这种现象说明随着城镇人口不断增加,生活方式不断发生改变,对城镇的资源结构和城市公共设施的需求的也不断提高,进而增加对能源的需求,从而造成了碳排放的增加。而本地区的城市化率对其他地区的碳排放的作用效果不显著,即空间溢出效应不显著。3) 产业结构(is)。产业结构的直接效应更大,且通过了1%水平的显著性检验,间接效应为较微弱。由于第二产业是我国经济高速发展时期的主导产业,且第二产业是碳排放的主要来源之一,因此产业结构对碳排放的增加具有显著的正向影响,产业结构对本地区的碳排放具有显著的正向作用,但每个地区经济实力不通、产业结构也各有差异,主要受本地区因素影响制约,对相邻地区的影响就较微弱。4) 环境规制(lner)。环境规制对碳排放的直接效应为负,间接效应为正,尽管尚不显著,表明环境规制对碳排放的抑制作用主要限于本地区,形成对本地碳排放的压力和减排动力,但可能由于地区之间大多各自为政,未形成政策的扩散与流动,所以对其他地区的碳排放不明显。5) 外商直接投资(lnfdi)。外商直接投资对本地区的直接效应与其他地区的空间溢出效应均为正向作用,且均通过了1%的显著性水平检验。这表明增加人均外商投资,虽然可能会给地区带来经济的增长,同时低质量的外商投资也会将碳排放转移到我国,造成本地区碳排放的增加。同时由于市场的临近性,一个地区的人均外商投资不仅影响本地区市场的调整,也会引起跨地区整合的供应链活动。
以上研究发现突破式低碳创新与增量式低碳创新对碳排放都具有显著的抑制作用,但以上研究只将两者看作独立的个体,尚未考虑两者同作为低碳创新活动的子项目,组合起来是否会对碳排放产生更大的效应?基于此,将突破式低碳创新与增量式低碳创新进行交叉相乘作为组合项,带入原模型,考量低碳创新组合是否会对碳排放产生更大的作用机制。模型的选择过程同上文,检验结果显示选择随机效应下的空间杜宾模型。模型的结果见表3。
SDM模型 | SDM模型分解 | ||||
---|---|---|---|---|---|
固定效应 | 随机效应 | 直接效应 | 间接效应 | 总效应 | |
disr | −0.0010 (−0.00124) | −0.0004 (−0.0014) | −0.000313 (−0.00145) | 0.00104 (−0.00314) | 0.000726 (−0.00359) |
incr | 0.0000 (1.80e-05) | 0.0000 (2.01e-05) | 7.15E-06 (1.95e-05) | 3.15E-05 (3.92e-05) | 3.86E-05 (4.46e-05) |
disr*incr | −9.28e-07*** (2.81e-07) | −1.13e-06*** (3.21e-07) | −1.22e-06*** (3.25e-07) | −1.75e-06** (7.99e-07) | −2.96e-06*** (9.51e-07) |
lnpgdp | −0.0194 (−0.0165) | −0.0284* (−0.017) | −0.0395** (−0.0174) | −0.177*** (−0.0455) | −0.216*** (−0.0548) |
ur | 3.495*** (−0.528) | 3.172*** (−0.56) | 3.163*** (−0.517) | −0.23 (−0.628) | 2.933*** (−0.495) |
is | 1.154*** (−0.246) | 1.160*** (−0.249) | 1.218*** (−0.244) | 0.675 (−0.504) | 1.893*** (−0.56) |
lner | 2.6860 (−2.259) | −0.7460 (−2.559) | −0.702 (−2.768) | 0.641 (−7.371) | −0.0613 (−8.529) |
lnfdi | 0.0527** (−0.0214) | 0.0309 (−0.0245) | 0.0496** (−0.0244) | 0.313*** (−0.12) | 0.363*** (−0.13) |
C | 11.38*** (−1.151) | ||||
rho | −0.1370 −0.0852 | 0.329*** (−0.0639) | |||
Log-likelihood | 234.7048 | 85.0275 | |||
R-sq | 0.1934 | 0.5893 |
Hausman | −24.3700 | |||
---|---|---|---|---|
LM-lag | 28.718*** | |||
Robust LM-lag | 27.259*** | |||
LM-error | 15.037*** | |||
Robust LM-error | 13.578*** | |||
LR-lag | 35.19*** | |||
LR-error | 34.28*** |
表3. 空间Dubin模型(SDM)回归结果及模型分解(2)
注:*、**、***分别表示通过10%、5%、1%的显著性检验;括号内的值为标准误。
从模型的回归结果可以明显看出低碳创新组合对碳排放的估计系数为负,且通过了1%的显著性检验,表明低碳创新组合对碳排放具有抑制作用。从模型的分解结果来看,无论是直接效应还是间接效应以及总效应,低碳创新组合对碳排放都具有反向作用,且分别在1%、5%、1%的水平上显著,表示低碳创新组合不仅能够有效的抑制本地区的碳排放,对改善其他地区的碳排放问题具有辅助作用,即空间溢出效应明显。
将两次模型的回归结果进行对比分析,在未加入创新组合前,突破式低碳创新、增量式低碳创新对碳排放具有显著的抑制作用,在加入创新组合变量后,突破式低碳创新对本地区的碳排放虽仍为负向作用,但作用效果不显著,且对其他地区的作用效果由反向作用变为正向作用,尚不显著;增量式低碳创新的作用机制在加入创新组合后,影响效果几乎未零,说明相较于低碳创新对碳排放的单独作用机制,创新组合对碳排放的影响效果更大,即要实现碳减排问题不仅要看创新某一单方面要素的影响,要尝试不同要素之间的组合,最大化的发挥要素之间的协同作用。其他控制变量的变化不是很显著。
本文在鉴别和筛选突破式与增量式低碳专利的基础上,以2006~2015年30个省份为研究对象,采用空间杜宾模型实证检验低碳创新专利及其组合对碳排放的直接效应和空间溢出效应。研究结论如下:1) 样本期间,中国的碳排放与低碳创新整体空间分布具有较强的集聚性,局部空间集聚性具有一定的差异,并且随着时间的变化,局部空间集聚性呈现动态演变特征;2) 突破式低碳创新对碳排放的影响强于增量式低碳创新,低碳创新组合的作用机制强于突破式与增量式低碳创新,且创新组合对其他创新的作用机制具有一定的排它性。但由于知识创新活动作用可能的滞后性,低碳创新活动及其组合的作用效果不是很突出。3) 城市化、产业结构升级以及外商直接投资等都会导致中国碳排放的增长,而经济增长并没有显著引致碳排放增加,环境规制则对本地碳排放起到了抑制作用,尽管对其他地区作用不显著。
以上研究结论得出的政策建议如下:1) 由于中国具有突破式的低碳创新专利作用效果仍较微弱,科技政策的重心应该进一步放在具有重要影响的突破式创新专利上,继续加入研发投入,突破常规的研发。2) 在发挥低碳创新活动的作用时,应注重不同的创新活动之间的关联性与互补性,通过不同类型创新的组合最大化创新活动的作用。3) 政府充分发挥宏观调控作用、整合优化资源配置,各地区积极构建信息共享机制,共同促进地方经济的持续发展和节能减排目标的实现。
教育部人文社科基金(16YJC630125):面向持续生态创新的我国制造业集群动态能力形成机制研究(201607-201912);
江苏省社科基金项目(17GLB020):生态创新与构建江苏制造业集群竞争新优势研究(201711-202012);
江苏省高校自然科学研究面上资助项目(17KJB170004):中国建筑业碳排放效率的空间格局演变及影响因素研究(201711-202012);
江苏省大创重点项目:中国如何进行颠覆式低碳创新?基于绿色专利数据挖掘的研究(201810299018Z);
大创项目:我国突破式创新对工业绿色增长影响的空间溢出效应研究(201810299240W)。
王 淼,王为东,卢 娜,张财经. 低碳创新及其组合对碳排放作用的空间外部性研究 Study of Spatial Externalities of Low-Carbon Innovation and Its Portfolio on Carbon Emissions[J]. 低碳经济, 2019, 08(01): 19-29. https://doi.org/10.12677/JLCE.2019.81003