-
权威、可靠、全面的数据是定性与定量分析的基础。由于目前文献中并无可靠的广东各类能源与经济指标数据集,首先需使用数据挖掘方法收集整理相关指标数据为后续开展量化分析提供条件。
广东省官方统计数据主要由广东省统计局发布,包括各年底发布至其前一年为止各种统计数据的广东省统计年鉴[16]和与在广东省统计信息网[17]发布,包含部分部门月度、季度数据的各种统计报告。综合考虑数据可用性与相关文献,选取表1中指标进行量化分析,并编写程序从对应原始数据发布页中收集整理。其中,广东工业分行业增加值[18]包括了精确到41个工业门类的月度数据。
数据名称 时间范围 周期 单位 广东省地区生产总值 2004-2023 年/季 亿元 广东省工业分行业增加值 2002-2023 年/月 亿元 广东省总固定资产投资 2006-2023 年/月 亿元 广东省能源行业投资 2006-2023 年/月 亿元 广东省总就业 2014-2022 年/季 万人 Table 1. Data used in quantitative analysis
根据国民经济行业分类(GB/T 4574—2017),能源行业由包括煤炭开采和洗选业、石油和天然气开采业的能源开采业,由包括石油、煤炭及其他燃料加工业、电力、热力生产和供应业、燃气生产和供应业的能源加工业组成。分析中的能源行业增加值、能源行业投资由上述分行业数据加总计算。
2000年以来,广东省能源行业总增加值持续增长(图1),从21世纪初的689亿元增长到2021年的3 593亿元,累计增长4.2倍,年均增速9.6%。分阶段看,2000年至2013年能源行业增加值增速快速攀升后震荡回落,总体增速保持在较高水平,年均增速12.6%;2013年后增速总体呈缓慢下降趋势,2021、2022年由于煤炭、天然气等大宗能源商品价格大幅上涨,增速明显回升。
从能源行业增加值的构成(见图2)看,电力热力生产和供应业一直是广东能源行业的支柱,占比最高达70%,近年来占比有所下降,但仍在50%左右。煤炭开采与洗选业自2005年广东退出煤炭生产后占比极低。随着2000年以来广东一批石化项目陆续投产、达产和天然气消费的快速增长,石油煤炭及其他燃料加工业、燃气生产和供应业的占比总体呈上升趋势,至2023年分别提升至约25%、10%。
能源产业的发展与国民经济整体结构有密切的关系。Simon Kuznets基于产业结构的划分标准[19]与Chenery的一般标准工业化模型中以人均GDP的判断标准[20]均表明,广东省在2012年前后进入工业化后期。从能源行业增加值占当年工业增加值、当年地区生产总值的比重变化趋势来看(见图3),2012年前能源行业占比呈震荡上升趋势,同期广东的电力供应较为紧张。2013年后随着经济进入新常态,能源需求增速降低,能源基础设施建设进度放缓,能源行业增加值占比也逐渐下降。“十四五”后期,“双碳”目标推动能源结构转型速度加快,海上风电、光伏与清洁煤电建设加快,加上同期高耗能产业投资大幅上升[17],拉动了能源行业增加值占比再次上涨。
-
本节利用广东省2014-2023年能源行业发展与经济增长指标数据,通过限制VAR模型与Feder生产函数对工业化后期阶段广东能源行业的经济贡献进一步展开定量分析。首先对文中采用的理论、方法及研究步骤进行简单介绍。
-
由于原始数据均为时间序列,在进行分析前需要使各数据变量的统计频率、时间范围和统计口径全部一致。在广东公布的官方数据中,整体GDP数据按季度公布,工业部门的增加值、电力消费、投资等数据却精确到月度,各指标变量的统计周期间存在较大差异。为保持一致,将月度数据均按对应季节加总后形成季度数据。由于广东省没有发布精确到能源或电力行业的价格指数,所有经济变量均保留现价计算方式。由于2018年后不再公布各部门固定资产投资绝对值,2018-2023年的固投数据均按公布的增速推算。
进行回归模型的拟合之前,一般需要通过平稳性检验来判断变量时间序列的平稳性,以避免出现“伪回归”,即由于变量随着时间变化的不平稳成分在拟合结果中占据权重过大,出现虽然拟合优度和显著性指标足够好,但实际无法真实反映变量间均衡关系的现象,造成结论可靠性较低。首先使用较为通用的增广式迪基-富勒(Augmented Dickey-Fuller,ADF)单位根检验测试各变量的平稳性,根据结果对变量统一进行差分、去除线性分量、取对数等平稳化处理。
协整检验是一种通过分析时间序列间协方差和相关系数来描述时间序列间长期相关性和稳定性的关系的统计工具。在分析中协整检验具有两方面的作用:(1)通过验证变量间长期稳定均衡关系的数量,为Granger因果性检验结果提供支撑,从而明确限制VAR模型中的约束项;(2)作为生产函数动态模型中建立误差修正项的参考[21]。本研究中由于同一面板内变量数量较多,选择使用更适合用于存在多对长期关系的情况下的Johansen迹检验。
-
VAR模型是一个多变量随机过程,用于概括多个随时间变化的量之间的线性相关关系[22]。在VAR模型中,每个变量的状态变化由模型中所有变量最大n期的滞后值和固定截距项决定。模型的阶数通常根据不同可能滞后值的多次试算得出的赤池信息准则(Akaike Information Criteria,AIC)选择。阶数为n的VAR模型的一般形式如式(1)所示,其中X是所有变量的状态向量,C和Ai分别是时不变的截距与各阶迟滞的系数矩阵,$ {\mathrm{e}}_{t} $是回归误差。确定VAR模型的参数需要计算所有系数矩阵Ai与截距矩阵C,可以等同于一个超定的线性系统参数识别问题,并使用回归方法求解。
$$ {\boldsymbol{X}}_{{t}}=\boldsymbol{C}+\displaystyle \sum _{{i}=1}^{{n}}{\boldsymbol{A}}_{{i}}{\boldsymbol{X}}_{{t}-{i}}+{{{\mathrm{e}}}}_{{t}} $$ (1) 限制向量自回归模型(Restricted-VAR)是VAR模型的一种特殊情况[23]。通过在参数回归过程中加入约束,可以间接规定系数矩阵C和Ai中指定元素的符号与大小,从而提高回归模型的现实意义。一般使用Granger因果性检验结果判断各变量间关系的显著性,并对VAR模型中对应不显著关系的系数施加零约束,使Granger检验判定不存在因果关系的变量间系数为0。
完成了VAR模型的参数计算后,可以使用冲激响应函数分析(Impulse Response Function,IRF)研究VAR模型的动态特性。在经典系统与控制理论中,一个线性系统的冲激响应函数是该系统的输入变量为单位冲激函数时,整个系统输出的经时间衰减的自然响应。由于VAR模型也属于线性系统,可以使用类似理论通过冲激响应对VAR模型进行研究[24]。通过对代表能源行业发展的指标施加冲激函数并观察经济增长指标的动态表现,可以定量分析能源行业与经济增长间的相互作用关系。
预测误差方差分解(Forecast Error Variance Decomposition,FEVD)是一种常用于分析时间序列模型中各变量对未来波动的贡献程度的技术。通过计算并比较VAR模型预测结果中的方差中来源于各输入的贡献,可以得出模型中各自变量对因变量影响力程度的相对大小[25]。
-
Feder模型是一种同等技术水平下生产要素的投入和产出间的数量关系进行定量分析的方法。该方法首次由世界银行的经济学家Gershon Feder于1982年提出[26],用于测算并对比各发展中国家10年间出口产业的经济贡献。Feder两部门模型理论将出口产业对各国整体经济的贡献分为直接贡献和间接贡献,将一国或地区的所有国民经济部门根据部门产品是否对外出口划分为出口与非出口两大部门,将国民经济领域的生产活动全部划分在这两个经济部门之一中进行。
基于Feder两部门模型构建生产函数,将整体经济产出分为能源和非能源部门,能源行业的产出作为非能源产业的生产要素建模,分析能源部门增长对整体经济增长的间接贡献。根据Feder两部门分析思想,建立如下生产函数:
$$ {Y}={f}({{K}}_{\mathrm{N}},{{L}}_{\mathrm{N}},{E})\mathrm{ }+\mathrm{ }{g}({{K}}_{\mathrm{E}},{{L}}_{\mathrm{E}}) $$ (2) 其中$ {Y} $代表整体经济产出,由代表全部非能源部门(下标N)的生产函数$ f $和代表能源部门(下标E)的生产函数$ g $组成。能源部门的产出由能源部门投资$ K_{{\mathrm{E}}} $和劳动力$ L_{{\mathrm{E}}} $决定,而能源部门的产出$ E $作为非能源部门的生产要素,与非能源部门的投资和劳动力决定整体经济产出。由于各指标数据单位不一、数量级相差较大,将原始数据归一化并统一化为百分比增速,根据式(2)中生产函数并参考Feder理论[21]建立如下线性回归模型:
$$ \begin{split} &\frac{\Delta {Y}}{{Y}}=C+\alpha \frac{\Delta {K}}{{K}}+\beta \frac{\Delta {L}}{{L}}+\gamma \frac{\Delta {E}}{{Y}}+\delta \frac{\Delta {E}}{{E}}+\\&\varphi {E}{C}{M}+\varepsilon \end{split} $$ (3) 其中,整体经济指标$Y$由地区生产总值代表;总资本投入$K$由广东省固定资产完成情况月报中总投资额代表;劳动力指标$L$由总就业人口代表;能源部门产出$E$由各能源部门增加值总和代表。在变量间存在协整关系时,使用原始数据回归试算得出残差的一阶滞后序列建立误差修正模型(Error Correction Model,ECM)加入模型中以增强短期动态回归精度。模型系数$ \alpha $、$ \beta $分别代表资本、劳动力的边际生产弹性,$ \gamma $可用于计算能源与非能源部门边际生产率的差异,$ \mathbf{\delta } $代表了能源部门对整体经济产出的溢出效应。
-
首先对未平稳化的原始数据进行Johansen协整检验验证各指标间的长期均衡关系。判定各指标间存在协整关系,可以支撑VAR模型建立前进行Granger因果性检验的结论。协整关系检验结果也可用于在生产函数回归中建立误差修正项,利用变量间的长期均衡关系根据预测误差动态修正每一期的回归与预测结果,从而提高模型的优度。
表2总结了对进行协整检验的结果。在5%显著性水平下,最多可以拒绝协整秩为0~3的零假设,表明数据横截面中的5个变量之间存在4个长期均衡关系。
虚无假设(H0) 迹统计量 显著性(5%) 结论 没有协整关系 136.901 69.063 拒绝 不大于1个协整 75.467 40.175 拒绝 不大于2个协整 37.280 24.276 拒绝 不大于3个协整 13.531 12.321 拒绝 不大于4个协整 2.9878 4.129 无法拒绝 Table 2. Johansen cointegration test results for raw data panel used in quantitative analysis
对原始数据及其一阶差分进行平稳性检验(表3),发现除能源投资外的4个变量均在5%显著性水平下与一阶差分后平稳。由于此4个变量达到平稳所需要的差分阶数相同,均为一阶单整,支持了协整性存在的结论。在后续建模分析中,由于VAR模型与生产函数模型的性质不同,在VAR模型中,对全部数据同时取对数处理以平稳化;在Feder模型拟合中,由于式(3)中的生产要素不包括能源行业投资,而其余变量均已平稳,故采用一阶差分的原始数据结合ECM进行生产函数拟合。
数据名称 前处理 ADF 显著性(5%) P值 结论 地区生产总值 原始数据 0.066 −2.964 0.964 非平稳 一阶差分 −4.294 0.005 平稳 总固定资产投资 原始数据 −1.882 −2.954 0.340 非平稳 一阶差分 −3.078 0.028 平稳 总就业 原始数据 −1.232 −2.949 0.660 非平稳 一阶差分 −5.492 0.000 平稳 能源行业增加值 原始数据 0.943 −2.964 0.966 非平稳 一阶差分 −3.117 0.025 平稳 能源行业投资 原始数据 −0.404 −2.968 0.909 非平稳 一阶差分 −1.638 0.463 非平稳 Table 3. Unit root test results for data used in quantitative analysis
-
为确定VAR模型中非零参数的数量与位置,首先对平稳化后的数据进行Granger因果关系检验。每次检验的虚无假设都是:原因变量的过去值对结果变量的未来值没有显著影响。表4中整理了各指标变量间的Granger检验结果,检验结果以矩阵的形式呈现,其中每个元素的行名和列名分别对应Granger因果关系的原因与结果,每个元素的值为拒绝虚无假设的p值。
结果变量 原因变量 生产总值 总投资 能源投资 能源增加值 总就业 生产总值 1 0.005 4 0.091 7 0 0 总投资 0 1 0.011 9 0 0 能源投资 0.012 8 0.057 3 1 0 0 能源增加值 0.001 1 0.002 4 0 1 0.141 4 总就业 0.016 9 0.031 5 0.010 3 0.174 0 1 注:p-值小于0.01%时记为0。 Table 4. Granger causality matrix for data used in restricted VAR model
检验结果表明数据横截面中存在多对Granger因果关系。除能源行业投资外,地区生产总值与其他变量的检验均在高显著性水平下通过,说明了此发展阶段下能源行业发展与整体经济发展间的互相促进关系。能源增加值与除就业以外的指标间都存在双向的Granger关系,说明能源行业与国民经济间关联较为明显。此数据横截面反映了广东经济结构持续升级和改革的过程中,不断提高的工业生产力创造了日益增长的能源需求,并刺激了能源基础设施的发展,这反过来又促成了对可靠和低成本能源供应具有强烈依赖性的高科技产业的建立。从长远来看,由于广东省政府制定了以实体经济为导向的发展战略,目前区域发展规划中包括了多个超超临界燃煤电厂和能源密集型计算中心与石化、化工项目,省内能源需求在“十五五”可能将持续增长。为了在满足能源需求的同时实现减排目标,未来的能源政策应重点推进清洁能源一体化和节能技术。
表5中Granger检验的结论被用来建立限制VAR模型。结果显示,能源投资到生产总值、总投资到能源投资和能源增加值与总就业之间的因果关系在5%的水平上不显著,因此在VAR模型中剔除相应项,通过在回归过程中加入零约束,强制其对应的系数为零。完成回归后,使用单位根检验模型的数值稳定性,即系数矩阵Ai的特征值全部位于复平面的单位圆内。
期数 生产总值 总投资 能源投资 能源增加值 总就业 1 0.553 0.190 0.042 0.155 0.058 2 0.435 0.238 0.036 0.200 0.089 3 0.389 0.255 0.032 0.219 0.104 4 0.370 0.262 0.029 0.228 0.111 5 0.359 0.265 0.028 0.232 0.115 6 0.354 0.267 0.027 0.234 0.117 7 0.351 0.268 0.026 0.235 0.118 8 0.349 0.268 0.026 0.236 0.118 Table 5. FEVD result for total regional production in VAR model
接下来使用冲激响应分析参数确定后VAR模型中关键变量间的动态表现。图4中展示了对各个自变量施加单位冲激后生产总值的短期累积动态响应。在每幅图中,蓝色实线代表各自变量施加冲击后10期内生产总值的累计变化,黑色实线是10期间的平均值。
从图4中可以看出,能源增加值的单位冲激对生产总值产生较为明显的正向影响,冲激发生后10期内收敛至约0.255个单位。表明如下一季度能源行业增长1个单位,会在之后数月内对整体经济产生相当于当季生产总值0.255个单位的累积影响。投资对生产总值的拉动效应则较不明显。从经济增长对能源行业的影响看,生产总值对能源投资、能源增加值的累计影响均小于0,在一定程度上侧面印证了广东能源消费与经济发展的脱钩趋势。在数据横截面的时间范围上,广东正处于能源消费与GDP增长脱钩的过程中。同期,越来越多的经济增长来自于第三产业的扩张,而第三产业的能源强度较低,并不会引起能源需求的剧烈增长。
最后对VAR模型进行FEVD分析研究各因素对经济增长指标变化的相对影响力。表6中展示了对最多8期生产总值对于各变量FEVD分析的结果。从表中看出,生产总值的演变动态在短期主要来源于自回归项,受其他变量影响相对较小。在各影响因素中,总投资与能源增加值变化对生产总值变化的解释率较明显,分别达到25%与22%左右。总就业对于生产总值的影响在初期较小,但4期后即上升至10%以上。
回归模型指标 数值 R2 0.869 F-统计量 13.23 P值 0.000 385 AIC −65.86 变量 系数 显著性(P值) 固定资产投资 0.437 9 0.017 总就业 −0.935 7 0.087 能源增长/生产总值 −2.079 7 0.651 溢出效应系数 0.387 1 0.05 误差修正项 −0.878 0 0.005 Table 6. Regression parameters of Feder production function
-
将历史数据代入(3)式中建立的Feder生产函数回归模型,使用生产总值、总投资、总就业、能源行业增加值的增长率以及能源增加值增长与生产总值的比值等变量进行直接回归并计算残差。所用变量均为一阶单整且存在协整关系,需要将残差的一阶滞后作为额外误差修正变量进行重新拟合并使用包含误差修正项的结果作为分析参考。
使用最小二乘法回归计算各项参数并列入表6中。观察各模型指标可以发现,反映拟合优度的R2指标为0.869,说明几个自变量指标的变化可以解释86.9%生产总值的变化;P值为0.000385,说明模型的置信度在0.0385%水平,优于一般经验的5%标准。分析模型中各项参数,其中代表总就业的系数显著性较低,对整体经济影响较为不明显,支持了表5中对VAR模型的方差分解分析结论。总投资、能源行业溢出效应指标与误差修正项系数显著性超过5%水平,能源行业溢出效应的系数$ \mathbf{\delta } $为0.3871,说明在其他条件不变的情况下,能源产业增加值增长1个单位,可以通过溢出效应促进非能源部门产出增长0.387个单位。与文献中类似结论[27]对比发现,此结果稍高,推测可能由于广东省外电、一次能源输入占比较高,而这部分能源的增加值未纳入广东省统计数据,导致能源行业的对应系数偏高。最后,误差修正项参数(即调整系数)小于0并在0.5%水平显著,说明短期波动偏离均衡时,以87.8%的速率将因变量拉回均衡状态,从回归角度印证了各变量间长期均衡关系的存在。
图5使用回归后的模型代入各自变量数据计算生产总值(红色虚线)并与真值(蓝色圆点)进行比较。观察模型与真值对比结果,可以发现此次回归精确度较高,较为有效的预测了因变量的变化趋势,对几个波动较大的数据点预测较为精准,侧面说明使用能源增加值对生产总值建模估算有较明显的作用。
Economic Contribution Measurement of Energy Sector Development During Late Industrialization Period in Guangdong Province
doi: 10.16516/j.ceec.2024-058
- Received Date: 2024-03-04
- Accepted Date: 2024-04-16
- Rev Recd Date: 2024-03-26
-
Key words:
- energy planning /
- energy policy /
- energy economic /
- data science /
- statistical analysis
Abstract:
Citation: | WU Dongqi, GAN Junwen, GUO Jingtao, YU Xinmei, LIAO Yi. Economic Contribution Measurement of Energy Sector Development During Late Industrialization Period in Guangdong Province[J]. SOUTHERN ENERGY CONSTRUCTION. doi: 10.16516/j.ceec.2024-058 |