最大用电负荷与经济变量的关联研究

    作者简介:朱平芳,上海社会科学院数量经济研究中心研究员(上海 200020)
    作者简介:谢婼青,上海社会科学院数量经济研究中心、上海社会科学院经济研究所助理研究员(上海 200020)
    作者简介:刘盼盼,国网上海市电力公司经济技术研究院经济师(上海 200120)

摘要: 用电负荷的估计与预测是经济社会稳定发展的重要因素。以用电负荷数据为导向,发现在高温天气当降温用电设备处于满负荷运行的状态下,用电负荷基本上由经济因素所决定,从而分离经济和气温两个对电力需求影响的因素,提出一种对不可观测最大基准用电负荷的估计方法以准确预测最大用电负荷。运用特征指数方法构建最大基准用电负荷指数,探索其与经济变量之间的关系。以上海为例构建年度和季度的最大基准用电负荷指数,发现其与经济发展高度拟合,反映近年来上海在产业结构调整方面卓有成效。本文以数据为导向,挖掘数据的充分信息,从新的视角估计最大基准用电负荷,以其为桥梁达到准确预测最大用电负荷的目标,为经济形势与转型的预判提供有效分析和数据支持。

English Abstract

    全文HTML

    一.   问题提出与文献述评
    • 十八大以来,中国经济逐渐从高速增长转向高质量发展阶段,推动经济发展质量变革、效率变革、动力变革的抓手便是生产方式、产业结构的不断调整和优化升级,电力消费是产业结构调整的客观数据,也是衡量经济高质量发展阶段性的重要指标,可以说,电力指标是经济发展的晴雨表。用电负荷,即电力负荷(Electrical Load),是用电设备在某一时刻从电力系统取用所消耗的电功率总和,依据用户使用电能的用途不同,可分为城乡居民生活负荷、商业负荷、工业负荷、农业负荷等。随着经济社会的发展,基础设施建设的完善,用电设备的数量和种类逐渐增加,用电负荷也相应提升,因此,用电负荷的数据客观上反映了经济社会的发展情况。最大用电负荷是指用电负荷的最大值,不仅受到气温的影响,还受经济社会发展的影响。通过分解出经济因素对最大用电负荷的影响,能为经济形势的分析提供客观数据的支持。

      但是,用电负荷是复杂的数据,不仅反映经济运行状况,还包括气温、人口等对电力需求的影响;此外,用电负荷是滞后的数据,电力部门需要利用历史数据来准确地预测短期、甚至是中长期用电负荷,尤其是最大用电负荷,才能保证满足城市化进程中居民生活和生产的电力需求、减少过度电能的浪费和损耗。因此,最大用电负荷的预测是影响整个电网安全稳定运行的基础,是电力系统调度的关键。尤其面对千万人口的超大城市,如何预测最大用电负荷便成为关乎人们正常生活和城市稳定和谐的重要问题。

      用电负荷的预测按时间分类,可以分为长期预测、中期预测和短期预测,一般而言,长期预测是指一年至十年,中期预测是指一月、季度或半年,短期预测是指一天或一周。用电负荷的估计主要分为传统的统计回归方法和数据挖掘的方法,传统的统计回归方法包括线性回归、时间序列等;数据挖掘的方法包括支持向量机、人工神经网络等。Amjady提出用人工神经网络和Box-Jenkins方法预测峰值负荷

      自2014年11月国务院对城市规模划分标准调整以来,常住人口超过1000万的超大城市充分体现了中国改革开放四十年来城市化进程的成果。目前,中国已经有以北上广为代表的十个超大城市。随着经济进入新常态,传统产业不断优化升级,加快发展先进制造业,超大城市的经济社会业态也逐渐发生转变,这种转变如何通过经济变量来体现和解读,本文提供了新的视角。以上海为例,2018年三产增加值22842.96亿元,同比增长8.7%,占生产总值的比重是69.9%,占比逐年增加,战略新兴产业和制造业的层级逐步提升以及新的商业模式与服务业态的不断形成,对城市化特性所彰显的用电负荷趋势产生了重要影响。

      在研究电力需求的文献中,都提到人口增长是决定电力需求的重要因素,较高的人口增长率会增加电力消费。随着中国长期以来的快速工业化和城市化进程,电力需求也呈现指数型增长。但是,就超大城市来说,2018年上海市常住人口达2423.78万人,同比增长0.02%,近三年的人口增长率都低于1%,由于超大城市的人口规模已接近其承载力的极限边缘,人口基本处于饱和状态,人口增长趋于平缓,人口增长率非常低,因此,对电力需求的增加也较为平缓,在本文中假定其并非重要因素。

      用电负荷的估计与预测是电力调度的重要参考依据,是经济社会稳定发展的重要因素;用电负荷数据的客观性为我们解读经济发展和转型分析提供了新的视角。但是,用电负荷数据的复杂性使得学术界和电力部门一直难以对其准确估计和预测,其中的主要困难来源于用电负荷受到经济因素和气温因素的双重影响,且不可观测。在高温天气当降温用电设备处于满负荷运行的状态下,用电负荷基本上由经济因素所决定。因此,本文尝试在高温条件下选择合适的温度(称之为基准温度),达到估计不可观测的最大基准用电负荷及其与经济变量的关联性,提供较为有效的预测分析经济与电力互动方法的主要目标;由于大量数据分析认为两者的差别仅由温差引起,将最大基准用电负荷及其指数作为桥梁,能达到准确估计和预测最大用电负荷的目标。

      因此,本文不同于以往针对用电负荷估计方法的讨论,以用电负荷数据为导向,挖掘数据中的充分信息,考察城市的经济和气温两个因素对电力需求的影响机制,提出一种对不可观测最大基准用电负荷(Vr)的估计方法以准确预测最大用电负荷(Vm),用特征指数(Hedonic)方法建立最大基准用电负荷指数,为经济形势与转型的预判提供有效分析方法和数据支持。

    二.   经济与气温因素对用电负荷影响的内在机制
    • 本文假定超大城市的电力需求主要由经济与气温两个因素影响。自改革开放以来,上海市电力消费随着GDP的增长而增长,2018年的电力消费达1566.66亿千瓦时,同比增长2.61%,用电负荷也随之正向变化,2018年的上海电网整点最高负荷达3094万千瓦,同比下降5.33%;而气温因素有两个方向,当气温升高(如夏季),用电负荷与气温的变化是同向的,当气温降低(如冬季),用电负荷与气温的变化是反向的。在实际社会中,因生活、生产等用电的复杂性使得主要宏观经济变量与气温对用电负荷的作用也呈现复杂的非线性特征,是经济和气温的多值函数。由于对用电负荷影响的作用机制不明晰,经济与气温各自对用电负荷的作用互相融合、相互干扰,使得对最大用电负荷的数据分析带来了实质性的困难。

      首先,我们分析经济和气温因素对电力变量的函数关系:

      其中,Y是可观测的电力变量,如负荷、电量等;E是宏观经济变量,如GDP、二产或三产的产出、工业增加值等;T是气温;F是未知的函数关系。由统计数据和实际经验易得出,YE同向变化,而T在夏季时与Y同向,T在冬季时与Y反向,从而构成与经济与气温的非线性关系。

      用电负荷可以分为非空调用电负荷与空调用电负荷,定义为“基础用电负荷”和“降温(取暖)用电负荷”,基础用电负荷与经济发展程度有关,经济越发达,工业、服务业所需电量越大,在一定阶段,基础用电负荷是趋于稳定的;降温(取暖)用电负荷与温度有关,夏季温度越高(冬季温度越低),空调需要开启来达到降温(取暖)的功效;但是随着经济的发展,基础设施建设不断完善,空调装机数量的增加,降温(取暖)用电负荷也会随之增加。我们将这种关系简化为数学表达式,记Yi是时点i的用电负荷值,Xi是时点i的基础用电负荷值,Zi是时点i的降温(取暖)用电负荷值(i=1, 2, 3, …n)。所以,某一时点的用电负荷值包含基础用电负荷和降温(取暖)用电负荷两部分:

      假定春季和秋季气温凉爽,不使用空调(取暖设备),则用电负荷不包含降温(取暖)用电负荷,故可用春季和秋季的用电负荷平均值来替代Xi,从而我们能估算该年的基础用电负荷和降温(取暖)用电负荷。

      我们以日最大用电负荷为分析对象,可知日最大用电负荷的关系如下:

      以上海为例,包括第四部分的实证分析,本文所采用的用电数据是来自上海市电力公司的每天在二十四个时点上的用电负荷数据,时间段是2000年3月1日—2017年8月15日,共6377天,其中包括横表24时点负荷(共149040个时点负荷值,非闰年每年有8760个数据)、日最大负荷、日最小负荷、日电量和一、二、三产及居民生活用电的月度数据;温度数据是来自上海市电力公司的日最高温度、日最低温度数据;经济变量的数据是来自上海市统计局2001—2018年统计年鉴的GDP数据,一、二、三产季度数据和人口数据等。其中,本文的四季划分如下:冬季是12月11日—2月,春季是3月—5月,夏季是6月—9月,秋季是10月—12月10日,故全样本共七十季,即2000年第一季度—2017年第二季度。

      图1是上海市十八年的四季日最大用电负荷图,从图中我们可以看到,每个夏季或冬季的日最大用电负荷在季初和季末与春季或秋季的日最大用电负荷基本呈平行关系,这意味着YiXi的峰值之间是近似平行的关系;随着时间的推移,日最大用电负荷同比呈增长的趋势。从以上数据结构我们总结两个实际用电的规律:(1)在一定时段内,基础用电负荷峰值序列存在一种固定的结构(图1中春季和秋季的序列);(2)在特定时间段和一定气温范围内,降温(取暖)用电负荷值变化不大。

      图  1  四季日最大用电负荷图

      图2是选取近三年的四季日最大用电负荷图,从图中我们可以看到,每季的日最大用电负荷数据具有明显的季节性特征结构,春季和秋季的日最大用电负荷基本处于平稳的态势,进一步说明将春季和秋季的用电负荷作为该年基础用电负荷估计降温(取暖)用电负荷是可行和合理的。

      图  2  近三年的四季日最大用电负荷图

      用春季和秋季的平均值作为基础用电负荷的估计,使用横表24个时点的负荷数据,我们估算了2000—2017年夏季一天24个时点上的最大降温用电负荷(如图3所示)。

      图  3  每年夏季降温用电负荷最大值图

      图3中可以直观地看到自2000—2017年以来上海市每年夏季降温用电负荷在各个时点上的最大值曲线始终呈现出“二谷二峰”的变化态势。夏季下午14点左右的峰时最大,降温用电负荷值最高,其次是晚上21点左右;早晨6点左右的谷时最大降温用电负荷值最低,其次是傍晚18时左右。十八年来,夏季各时点的上海最大降温用电负荷值都逐步上升,体现了经济的发展增加了电力需求。但是,金融危机时期的2009年、经济转型时期的2013—2015年,上海的最大降温用电负荷表现异常,其中的原因是复杂的,很大程度上与这几年的极端高温天气有关,也有经济转型因素的干扰。下文在估计最大用电负荷和最大基准用电负荷中会进一步处理两者的影响。

      分析日最大用电负荷年平均值与GDP对数值,我们对两者分别与前一期的差分做相关系数,发现春季的相关系数高达70.94%,秋季的相关系数是27.69%,春季日最大用电负荷年平均增长率与GDP增长率的相关度更高,秋季较低。从图4中可以看到,春季和秋季的基础用电负荷随着时间的推移而呈增长趋势。春季和秋季的气温适宜,经验可知,在这两个季节里,居民生活中基本不开空调,故受气温的影响非常小。

      图  4  上海历年春秋季最大用电负荷与经济变量

      因此,我们总结出以下结论:(1)在温度适宜的情况下,基础用电负荷主要受经济变量的影响;(2)降温(取暖)用电负荷不仅受到气温的影响,还受到经济变量的影响,即随着经济的发展带来的人民收入的增长,空调装机的数量增加,因此,用电负荷也随之增加;(3)在特定时间段内,当温度一定(夏季气温足够高)时,空调(降温)设备接近完全启动,用电负荷(降温)基本趋于稳定。因此,当气温足够高时,我们可选择适当的温度值定义为基准温度,其对应的用电负荷称之为最大基准用电负荷(下面会给出精确定义)。最大用电负荷围绕最大基准用电负荷的波动仅与温差有关,这就可用于解决下文的最大基准用电负荷的估计和最大用电负荷的预测问题。

    三.   估计模型与指数编制方法

      一.   最大用电负荷估计和预测模型

    • 考虑到式(1)中经济变量E取定值,当温度T足够高,最大用电负荷趋于基本稳定。因此,可以认为此时的最大基准用电负荷是稳定的。

      为此,假定经济变量取定值E=Er,温度达到较高的值Tr(本文称之为“基准温度”),存在一个客观的最大用电负荷Vr,其含义是经济发展到一定阶段,温度足够高时所需配置的用电设备最大用电功率的加总,本文称其为“最大基准用电负荷”。Vr存在但是不可观测。实际可观测的最大温度是Tm,最大用电负荷是Vm,由于基准温度足够大,在基准温度附近的负荷值趋于稳定,因此,最大用电负荷(Vm)围绕最大基准用电负荷(Vr)上下波动,这种波动幅度小,且可认为仅由温差TmTr引起的。

      基于上述分析,最大基准用电负荷与最大用电负荷的函数关系式如下表示:

      由此我们可构建最大基准用电负荷与最大用电负荷的对应关系如下:

      其中,$H\left( {{T_m} - {T_r}} \right)$是以Tm- Tr为自变量的未知函数,ε是随机误差项,将eε看成一个正向的扰动项。假定$H\left( {{T_m} - {T_r}} \right)$能在原点附近以幂级数展开,则可以得到:

      若考虑n期时间上最大用电负荷值的变化(t=1, 2, 3, …, n.),则在每期的Er, VmTm不同,但是基准温度Tr相同,故可将式(5)写成时间序列形式的n个方程:

      式(6)中已知Vm,t, Tm,tTr,求解Vr,tKi,t存在困难,为求解方便,我们简化式(6)。(1)由于函数$H\left( {{T_{m,t}} - {T_r}} \right)$的性质与Er无关,并且时间的变化仅表现在最大温度Tm,t上,即有Ki,t=Ki,对所有的i,t均成立;(2)$H\left( {{T_{m,t}} - {T_r}} \right)$事实上是一个很小的量,将其以幂级数完全展开是没有必要的,一般使用一次或二次项就能保证精度,故而取一次型(i=1)近似得到:

      K1可看作单位温差Tm,t-Tr下最大用电负荷和最大基准用电负荷之间的变化率,比较Tm,tTr,将K1分为两段,以提高精度。因此,若得到K1的值,通过式(7)便可得到最大基准用电负荷Vr,t的序列。

      如果用直接的代数迭代算法来计算式(7)中n+2 个未知参数,需要假定诸多前提条件,这会造成计算的繁琐并且精度较差。本文用统计拟合方法估计最大基准用电负荷序列,算法思想直观,且计算简洁方便。

      首先,将式(7)取对数变形如下:

      这里,K1同式(7),εt是误差项。式(8)是一个线性回归模型,待估参数是Vr,t, K11K12。进入21世纪以来,上海市的经济和城市化进程都得到了快速发展,所以,直观上可知Vr,t是在n期上随着时间t逐年增加,且呈指数型增长。

      其次,为估计式(8),我们需要选择 ln Vr,t的初始序列 {at},使其与 {lnVr,t} 的趋势保持一致,考虑社会经济实际情况,我们认为 {lnVr,t} 的趋势蕴含在 {lnVm,t} 中,为了更好地去除噪声对 {lnVr,t} 的影响,对{lnVm,t}作三期移动平均(实证部分基于数据特征采用两期移动平均),记为:

      再将$\left\{ {\ln {V_{m,t}}} \right\}$bt做回归以提高计算精度,故而得到回归系数$\hat \beta $,我们令初始序列为:

      则将模型式(8)改成如下形式:

      其中,${\rm{ln}}\;\;{V_{r,t}} = \alpha + \gamma {a_t}$K1在线性模型式(9)中是哑变量,同式(7)。因此,式(9)中的待估参数是α, γ, K11K12,共4个,只要样本量足够大,用OLS方法就能估计这些参数。此外,我们在估计过程中将时间t*上的异常值已剔除。

    • 二.   最大基准用电负荷指数

    • 最大用电负荷估计和预测模型是一种有效分析电力数据的实用算法,但是其分析的对象是每年的最大用电负荷,因此,数据使用效率大大降低,导致损失了部分有用信息。该估计方法没有考虑负荷的堆积效应,即持续几天高温后常会有高峰负荷;用一天的数据解释全年的经济变化有点牵强;用模型预测日用电负荷的误差较大等。为此,本文在最大用电负荷估计方法的基础上,用特征指数(Hedonic)技术建模,从而形成最大基准用电负荷指数,该指数可以估计月度、季度甚至频率更高的数据。

      特征指数方法通常应用在编制住房价格指数方面,价格指数用以下形式表示:

      其中,Indext是价格指数,Pt表示交易价格,${\hat P_t}$表示通过模型拟合的价格。特征价格函数FtX1, X2, …, Xk)表示价格是特征(X1, X2, …, Xk)的载体。重复交易模型(Repeat Sale Method)是将一定时期内同一住房单元交易的两次及两次以上的价格都计入指数的数据范围,保证指数编制的同质性。本文在构建最大基准用电负荷指数时采用特征指数方法和重复交易模型的思想,将时间序列的样本点作为特征处理,从而构建模型。

      最大用电负荷估计模型已经证明了最大基准用电负荷曲线是客观存在的,且在一定条件下趋于稳定。据此,对峰值负荷而言,当气温足够高时,存在一条基准用电负荷指数,也就是说当气温大于临界值时,峰值负荷也是某种形式的“平均”。

      首先,界定临界气温T0。定义样本内的用电负荷值是指在该期最大用电负荷值附近,且气温大于临界值的数据。我们以季度数据为例,在统计出每季最大用电负荷样本后,构建如下乘法模型:

      其中,t = 1, 2, 3, …, n(年)。这里,t表示年度标记,j表示样本点标记,Vtj表示进入样本内的最大用电负荷值。$\alpha _1^{T_{tj}^1} \times \cdots \times \alpha _l^{T_{tj}^l}$表示积温的影响,其中α是待估参数,$T_{tj}^l$$V_t^j$样本点所对应的时点气温值,$T_{tj}^k\left( {k = 1, \cdots ,l - 1} \right)$Vtj样本点前若干天的气温值,以考虑用电负荷的堆积效应,l是积温期长度,一般是二到四日;$\beta _1^{{h_1}} \times \cdots \times \beta _6^{{h_6}}$表示每周不同日的影响,通常工作日与周末的用电负荷数量和结构都不尽相同,因此,本文将周一至周六不同日用哑变量huu=1, …,6)表示,Vtj所对应星期的hu取1,其余取0,其中,周日全部取0;It是待估参数,即Vtj样本点所对应的最大基准用电负荷指数值,n是总年度数,st是哑变量,当Vtj是第t年的最大用电负荷,则为1,否则是0;ε是误差项。

      模型(10)式是利用基期和多个样本期所建立的特征指数(Hedonic)模型,在模型中还以星期的哑变量标识样本所属的不同日期,以年份的哑变量标识最大基准用电负荷指数值的不同形式,从而在原有的最大用电负荷估计模型的基础上添加更多有效的信息特征,构成完整的最大基准用电负荷指数模型。

      为求解模型式(10),首先对其取对数简化如下:

      其中,t = 1, 2, 3, …, n(年)。从而将式(10)转换成典型的线性模型,待估参数包括αu, βvIw,通过线性OLS估计方法,可以估计$\dfrac{{{I_w}}}{{{I_{w - 1}}}}$,从而得到最大基准用电负荷指数。

      下文是对最大用电负荷估计模型和最大基准用电负荷指数方法的实证应用。

    四.   实证结果分析
    • 实证分析使用的数据与上文使用的数据相同。本文以上海市2000—2017年电力数据、温度数据和经济变量的数据为例,用最大用电负荷估计模型估计Vr, tK11K12,预测Vm , t,用特征指数(Hedonic)方法构建最大基准负荷指数,并对其给予经济学意义的解读和分析。

    • 一.   数据描述性统计

    • 下文使用上海市电力公司的日用电负荷数据表为主要研究对象,为了对上海市的电力数据、温度数据有直观的认识,我们分析了上海市2000年3月1日—2017年8月15日电力公司共6377天的用电负荷数据和日温度数据,变量包括日电量、日最大负荷、日最小负荷、日最高温度和日最低温度,表1对全样本和各年份的数据分别进行了描述性统计。

      电力数据
      日电量(亿千瓦时)日最大负荷(万千瓦)日最小负荷(万千瓦)
      年份平均数中位数标准差平均数中位数标准差平均数中位数标准差
      全样本2.983.020.981463.901467.60487.48948.24973.10307.48
      20001.501.450.19744.24723.35103.45473.24461.1554.08
      20011.581.530.21785.58761.10115.92498.15485.6055.79
      20021.721.700.23857.91841.20127.10540.20533.0059.70
      20031.991.910.29989.06958.70152.90626.09603.4081.39
      20042.202.130.291083.481051.40148.94699.87679.3579.86
      20052.472.390.341214.751175.80174.12785.26762.1091.54
      20062.692.580.421321.521272.50216.65857.28826.00114.22
      20072.922.800.461436.361385.60242.58931.12903.30127.84
      20083.092.960.471520.581446.65245.22990.15952.40132.74
      20093.143.070.471544.851509.20246.671001.10983.90133.02
      20103.533.330.571727.291629.30298.741137.281080.70164.08
      20113.663.530.531794.481740.90283.991174.521139.00149.32
      20123.683.590.551821.171772.10291.651171.921138.35154.95
      20133.853.590.721894.931774.40376.321220.781141.00207.40
      20143.733.670.501820.761789.90269.921183.281168.70136.23
      20153.833.740.551877.141832.70295.011219.601196.00149.15
      20164.033.840.711977.621877.05378.351290.491227.35202.47
      20174.183.980.882037.561945.70453.691337.441275.00264.59
      数据来源:笔者依据上海市电力公司用电负荷数据整理所得。

      表 1  日用电负荷数据的描述性统计

      表1中可知,随着时间的推移,日电量、日最大负荷、日最小负荷三个变量都逐年增加,主要是因为经济社会环境的发展,对电力的需求逐年增加,从一定程度上反映了上海市十八年来的工业、服务业的发展和城市化进程。三个指标的波动程度也随之更加剧烈,说明经济社会的工作时间和非工作时间的特征更加明显,而特征指数算法下的最大基准用电负荷指数法很好地使用了数据背后更为丰富的信息,将不同星期作为哑变量放进模型中。此外,日最大负荷与日最小负荷的差距逐年扩大,这差距中不仅包含了经济因素还有温度因素的混合。

      本文用“用电负荷旋风图”来表示以上数据特征。我们将用电负荷数据为横轴,以季节序号为纵轴,将6377个观测点画到同一个图内(如图5)。图5中每一横条上的点是当季日最大用电负荷(深色点)和日最小用电负荷(浅色点)的排序,可以清晰地看到样本内七十季日用电负荷的分布结构像一个旋风(Cyclone)的投影,直接反映了上海市2000年以来经济社会的迅猛发展,即随着时间的推移,日最大用电负荷和日最小用电负荷的最大值和最小值不断增加,最大值与最小值的波动也在不断增加。

      图  5  用电负荷旋风图

      表2反映了温度具有季节性变化,较少受到经济社会的发展影响,与经济因素基本是相互独立的。

      温度数据
      日最高温度(摄氏度)日最低温度(摄氏度)
      年份平均数中位数标准差平均数中位数标准差
      全样本21.5922.709.3215.0615.88.74
      200023.5624.607.7617.1118.007.41
      200121.3723.009.1014.8915.508.51
      200221.6621.608.5315.0415.607.91
      200321.5522.509.8514.7915.009.07
      200421.9623.159.0715.1215.658.62
      200521.5823.2010.2614.6516.109.61
      200622.1723.409.3715.5516.208.65
      200722.2722.408.9815.6815.308.33
      200821.3822.209.6914.6515.009.15
      200921.5524.209.4914.9216.708.90
      201021.2321.009.4914.6414.109.07
      201121.0022.009.8314.4316.009.28
      201220.6322.459.7414.5415.808.97
      201321.8122.409.9915.0815.209.33
      201420.8421.908.3314.8515.808.08
      201520.9523.208.5414.8016.008.00
      201621.4822.109.2615.4616.058.72
      201722.3022.9010.0115.3716.309.14
      数据来源:笔者依据上海市电力公司温度数据整理所得。

      表 2  日温度数据的描述性统计

    • 二.   最大用电负荷估计和预测

    • 以年最大用电负荷为研究对象,假定最高基准温度Tr为38.7摄氏度,表3分别是采用两期、三期和四期的序列{bt}移动平均所估计的参数值及标准误,移动平均的期数不宜过长,否则会造成趋势太光滑,使得精度变差。根据估计得到的结果可看到四期移动平均的调整的R2最低;三期移动平均次高;而两期移动平均最高。表4是用最大用电负荷估计模型预测得到的MSE,可以发现两期移动平均的MSE最低,精度最好;三期移动平均的MSE次低;而四期移动平均的MSE最高,精度最差。因此,采用两期移动平均更为合适。

      两期移动平均的估计值三期移动平均的估计值四期移动平均的估计值
      α0.2553[0.1504]0.6860***[0.1872]1.1076***[0.2166]
      γ1.0067***[0.0210]1.0022***[0.0274]1.0072***[0.0334]
      K11−0.0015 [0.0112]0.0077[0.0130]−0.0008[0.0143]
      K120.0287**[0.0128]0.0293*[0.0150]0.0399**[0.0166]
      调整的R20.99480.99110.9865
      p-value1.151e-153.926e-133.915e-11
      注:* 表示90%水平显著,**表示95%水平显著,***表示99%水平显著。

      表 3  最大用电负荷估计和预测模型

      两期移动平均三期移动平均四期移动平均
      MSE2662.1622919.6493318.067

      表 4  最大用电负荷估计模型的预测均方误差

      图6是在表3模型基础上的预测图,最大用电负荷预测值基本与最大用电负荷值(Vm)一致,并且估计得到的最大基准用电负荷(Vr)序列也能反映最大用电负荷值的趋势变化,其中两期移动平均的最大用电负荷估计值与其观测值基本重合,预测效果最好。

      图  6  最大用电负荷和最大基准用电负荷序列预测图

      本质上,最大用电负荷估计模型利用的是在高温气候和很高用电负荷的“极端”条件下,最大用电负荷仅与温度有关,而剥离出的最大基准用电负荷主要受经济变量影响的特点。这为我们对经济变量的解读、经济形势的预判和研究提供了新的视角和数据支持。以最大基准用电负荷为桥梁预测最大用电负荷的精准度较高。在以上的模型框架下,我们估计出不可观测的最大基准用电负荷(Vr)序列,该序列在固定基准温度Tr下,主要由经济变量所影响。我们在改变基准温度的情况下,分别用38.7°C、39°C和39.4°C三个基准温度,采用两期的移动平均计算最大基准用电负荷序列(如图7所示)。

      图  7  不同基准气温下的最大基准用电负荷Vr,t的估计序列

      图7中可知,随着基准温度的提高,整条曲线向上移动,这与常识相符,即最大基准用电负荷随着温度的上升而增加。但是三条曲线的趋势和形状几乎不变,这是因为曲线的形状是由经济变量所决定的。因此,最大基准用电负荷Vr,t的估计序列不仅能用来预测最大用电负荷(Vm,t),而且作为经济分析的变量,具有实用性和可操作性。由式(1)可知,选取的经济变量不同,最大基准用电负荷Vr,t的估计也会不同,即确定经济变量的类型后,需要构造相应的Vr,t序列,这依赖于初始序列 {at} 的选择。其中,经济变量是指全社会的经济状况,包含生产、生活、办公和商业等各个方面,下面我们选择生产过程中的经济变量为GDP来考察,由于GDP本身是平滑的,故不再做移动平均的处理。

      分别对应经济变量GDP,一产、二产、三产的GDP值作为最大基准用电负荷估计初始序列的“原材料”的 {bt},进而估算Vr,t序列[${V_{r,\;t}} \!=\! {V_r}\left( {GD{P_t},\;{T_r}} \right)\!$] 。我们定义GDP单位产值电耗估计值(${\hat V_{r,t}}$/GDP)。图8是在基准温度38.7°C计算的GDP单位产值电耗估计值。从图8(1)可以看到,十八年来上海市节能减排和产业转型是很有成效的,但是从整体图形来看,这种下降趋势逐步放缓,意味着降低GDP产值电耗会愈加困难,其中,三产的GDP单位产值电耗估计值的趋势与GDP的趋势一致,并且下降趋势更加缓慢。二产在2008年、2009年单位产值电耗最低,在2015年以后出现上升的趋势,2017年的单位产值电耗估计值出现“小峰”,推断可能是因为经济危机后期和经济转型过程中,全社会国民生产总值增速有所下滑造成的。一产在2012年之前电耗增加趋势明显,而近几年有显著的下降,体现上海市在乡村振兴战略方面得到了实质性的成效。从一定角度说明,上海市的三产受经济危机的影响不明显,二产在经济转型过程中影响显著。

      图  8  最大基准用电负荷与经济变量的关联关系

    • 三.   最大基准用电负荷指数

    • 我们以上海市日最大用电负荷数据为研究对象,数据总样本有6377个,选取三天的积温期,并将夏季和冬季积温期分别计算,要求进入计算的日最大负荷不小于该年日最大负荷的95%。通过OLS模型估计积温参数αu(6个),星期参数βv(6个),以及负荷指数IwIw-1的比值(70个),由于待估参数较多,文中不报告数值结果,可以计算出季度最大基准用电负荷指数(如图9)。

      图  9  最大基准用电负荷指数及光滑线

      在最大基准用电负荷指数的估计中,待估参数非常多,正如最大用电负荷估计和预测模型结果中表现的那样,存在着不显著的情况,我们通常在分析线性关系时看调整的R2。最大基准用电负荷指数模型中得到的R2=0.9984,均方误差为0.00025,模型拟合结果较好。

      图9是估计的最大基准用电负荷指数,可以明显地看到,该指数有四季的周期性特征,因此,对指数做四期的移动平均平滑处理。

      从平滑曲线中能直观地看到上海市全社会最大用电负荷的变化趋势,其中,2009年第二季度由于经济危机的显著影响,经济不景气明显;2014年以后,全国整体的经济形势从重视GDP增长到重视经济结构的优化升级,从要素和投资驱动逐渐转向创新驱动,最大基准用电负荷指数很好地反映出上海整体的经济逐渐趋向稳健;尤其在2016年末2017年初,经济结构的调整到了攻坚的阶段,高耗能低效落后产能的逐渐淘汰必然造成了这一指数趋向平缓甚至下滑。另外,这种指数方法还能构建月度等更高频率的指数。

      我们还将该指数与经济变量GDP,一产、二产、三产GDP做回归分析,从表5中可以看到该指数对经济变量的解释都具有显著性影响,其中,与GDP,二产、三产GDP回归的调整的R2较大,反映了三产所占比重最大、二产所占比重其次的上海市产业结构特征。因此,用最大基准用电负荷指数来分析经济形势和产业结构的变化是可行的和客观的。

      GDP一产GDP二产GDP三产GDP
      估计值0.7662***[0.0774]0.3819***[0.1113]0.7057***[0.0853]0.7794***[0.0754]
      调整的 R20.58100.13350.49080.6018
      p-value6.953e-150.00106.291e-121.18e-15
      注:* 表示90%水平显著,**表示95%水平显著,***表示99%水平显著。

      表 5  最大基准用电负荷指数与经济变量的关系

      我们用“最大基准用电负荷指数法”研究夏季用电负荷,构建夏季日最大用电负荷指数、夏季日最大降温负荷指数、基础负荷指数和夏季日用电量指数(如图10)。

      图  10  最大基准用电负荷指数光滑线

      图10中能直观地看到,日用电量指数与日最大用电负荷指数基本同趋势,这两个指数是经济与温度变量的综合反映;基础负荷指数较为平稳,只受经济因素的影响,不受温度的影响;最大降温负荷指数的波动与日最大用电负荷指数的波动一致,可以得出,日最大用电负荷的波动主要原因在于降温负荷指数的波动。以上四个指数对预测电力负荷和分析经济变量起着重要的作用。

    五.   结论与启示
    • “最大用电负荷估计和预测模型”和“最大基准用电负荷指数方法”是在电力峰值极端的情况下,固化了经济因素和气温因素的作用,以达到分离因素变量的目的。本文首先揭示了当夏季温度足够高时空调降温用电负荷趋于基本稳定,此时的最大基准用电负荷也相对稳定的事实,发现了经济变量一定时,只要温度足够高,最大用电负荷趋于基本稳定,从而能够在“极端”气温和用电负荷的条件下分析基础负荷和降温(取暖)用电负荷,提出一种对不可观测最大基准用电负荷(Vr)序列的估计方法,不仅能准确地预测最大用电负荷(Vm),还能分析经济转型过程中,一产、二产和三产分别的表现。

      本文通过特征指数(Hedonic)法建立季度甚至月度最大基准用电负荷指数,分析其与经济变量之间的关系;此外构建夏季日最大用电负荷指数、夏季日最大降温负荷指数、基础负荷指数和夏季日用电量指数,从而剥离出基础负荷指数平稳的特性,说明其与经济因素的关联关系,而不受气温因素的影响,日最大用电负荷指数与最大降温负荷指数的波动一致,体现日最大用电负荷指数与日用电量指数是由经济与温度变量的综合反映,其波动主要原因在于降温负荷指数的波动,这四个指数对于我们分析和预测电力负荷和经济变量具有良好的作用。

      以上海为例,我们发现上海在过去十八年来节能减排和产业转型是很有成效的,尤其是二产的GDP单位产值电耗估计值在经济转型过程中在2015年以后出现上升趋势,而三产受经济危机的影响不明显。2014年以后,全国整体的经济形势从重视GDP到重视经济结构的优化升级,最大基准用电负荷指数正好地反映出上海整体的经济逐渐趋向稳健;尤其在2016年末2017年初,经济结构的调整到了攻坚的阶段,高耗能低效落后产能的逐渐淘汰造成最大基准用电负荷指数趋向平缓。

      十八大以来,尤其是随着中国经济发展进入新常态的阶段,以北上广为代表的超大城市面临转型升级的要求,全国经济发展的“排头兵”是超大城市的责任和使命。但是,面对人口趋于饱和的现状,经济结构的优化使得高耗能低效率的产能逐步淘汰,从而出现经济发展的阵痛期。就上海而言,正处于转变发展方式、优化经济结构、转换增长动力的攻关期,需要寻找新的经济增长点,发展低耗环保的新经济业态,推动新型工业化、信息化、农业现代化同步发展,提高全要素生产率,推动经济发展质量变革、效率变革和动力变革。

      用电负荷的研究不仅可以应用在分析经济变量的趋势预测和分析上,还可以研究消费者行为等;用电负荷的数据庞大和复杂,未来我们会对用电负荷的数据所蕴含的信息通过数据挖掘和机器学习的方法进行更深的处理和分析。一方面,进一步挖掘和研究“用电负荷旋风图”,分析日最大用电负荷和日最小用电负荷分别与经济、气温因素的关联关系;另一方面,对夏季日最大用电负荷指数、夏季日最大降温负荷指数、基础负荷指数和夏季日用电量指数四个指数进行分析,探索其与经济变量之间的关系;此外,将来可以尝试用现代计量经济学的方法构建用电负荷与经济变量的模型,将本文提到的经济变量与气温对用电负荷的复杂非线性特征的多值函数加以模型化。

    参考文献
(10)  表(5)