一种基于CVMD-GRU-DenseNet混合模型的短期电力负荷预测方法

一种基于cvmd-gru-densenet混合模型的短期电力负荷预测方法
技术领域
1.本发明涉及电力负荷预测技术领域，具体涉及一种基于cvmd-gru-densenet混合模型的短期电力负荷预测方法。

背景技术：

2.负荷预测是电力系统运行决策的重要前提，准确的短期负荷预测不仅能为电力公司制定发电计划提供可靠依据，也有助于制定适当的备用和维护计划，降低运维成本。
3.短期电力负荷预测方法分为统计学方法和机器学习方法。统计学方法难以准确预测高复杂度负荷时间序列的未来变化趋势。而早期的机器学习方法如支持向量机(support vector machine,svm)虽对非线性负荷序列具有拟合能力强的优点，但难以挖掘负荷的时序性规律。近年来，深度学习方法因其强大的数据特征挖掘能力成为负荷预测的研究热点。其中长短期记忆网络(long short term memory,lstm)通过增加门控结构，解决了传统循环神经网络学习长时间序列的梯度消失或爆炸问题，被广泛应用于长时序预测。但受制于门控参数多以及时序信息的顺序传递机制，lstm的收敛速度较慢。门控循环单元神经网络(gated recurrent unit,gru)将lstm的遗忘门和输入门合并为重置门，通过简化结构，提高了收敛速度。卷积神经网络(convolutional neural networks,cnn)因其强大的局部特征提取能力，在时间序列预测中具有较为明显的优势，但随着网络深度的增加，存在训练结果不稳定、梯度消失等问题。残差网络(deep residual network,resnet)通过构造恒等映射设置旁路进行残差学习减少了梯度消失和模型退化问题。在resnet的基础上，密集连接卷积网络(densely connected convolutional networks,densenet)通过在各层间建立连接关系，减少了网络参数和计算成本，加快了特征信息传递速度。
4.通过信号分解技术准确分析和挖掘负荷变化规律是提高短期负荷预测精度的另一个重要手段，主要方法有小波分解(wavelet decomposition,wd)、经验模态分解(empirical mode decomposition,emd)和变分模态分解(variational mode decomposition,vmd)等。wd的分解效果受小波基函数与分解层的影响，适用性不强。emd可自适应地分解负荷序列但易发生模态混叠现象。vmd采用非递归和变分模态分解来处理原始信号，具有更强的鲁棒性，有效减少了模态混叠发生。
5.申请公布号cn114626581a公开了一种基于vmd分解机制的tcn-svm组合用电负荷预测方法，该方法未考虑如何确定vmd分解次数，分解次数设置过小，负荷功率时间序列并未分解完全，各子序列呈现规律不够明显，影响后续预测的精度；分解次数设置过大，会导致分解过度，相邻的子序列分量重复的问题。
6.申请公布号cn113902206a公开了一种基于vmd-bigru的短期负荷预测方法，该方法对经vmd分解后的子序列采用相同的预测模型，未针对各子序列的差异选择合适的预测模型。

技术实现要素：

7.本发明的目的是为了解决现有技术存在的利用vmd对负荷进行分解难以确定最佳分解次数以及根据分解后的各子序列特点及差异难以选择合适的预测模型的技术问题。
8.为了解决上述技术问题，本发明采用的技术方案为：
9.一种基于cvmd-gru-densenet混合模型的短期电力负荷预测方法，它包括以下步骤：
10.步骤1：采集目标地区一定时间段为单位的负荷功率时间序列数据集；
11.步骤2：将负荷功率时间序列分解为k个本征模态函数子序列分量和1个残差分量；
12.步骤3：筛选与预测分量关系密切的输入特征，依据各分量序列周期及复杂度的差异选用不同的预测模型；
13.步骤4：将各子序列分量的预测结果叠加，重构得到未来日时序负荷预测结果。
14.在步骤2中，采用依据相关熵确定最佳分解数的变分模态分解技术，将原始负荷时间序列分解为k个本征模态函数imf子序列分量和1个残差分量；具体包括以下步骤：
15.步骤2.1：设置分解次数k；
16.步骤2.2：输入步骤1中采集的负荷功率时间序列数据集，构造约束变分问题：
[0017][0018]
其中，{uk}是imf集合，{ωk}是中心频率集合，δ(t)是狄拉克分布，表示卷积，uk(t)代表第k个imf子序列分量，t为采样时刻点。
[0019]
步骤2.3：引入二次罚因子α和拉格朗日乘法算子λ，将式(1)转化为无约束变分问题，表示为：
[0020][0021]
其中，α是二次罚因子，λ是拉格朗日乘法算子，《,》表示内积；
[0022]
步骤2.4：采用交替方向乘子法求解式(2)，uk、ωk和λ迭代更新过程为：
[0023][0024][0025]
[0026]
其中，n为迭代次数，f(ω)、分别代表时域下的f(t)、经过傅里叶变换到频域的结果，τ为噪声容忍度；
[0027]
步骤2.5：根据公式(6)、(7)计算imf分量间的相关熵最大值，通过检测imf分量两两间的相关熵是否近似为1，从而判断对应分解数k条件下是否出现模态混叠现象，以此为依据确定最佳分解数k
opt
。
[0028][0029][0030]
步骤2.6：若v》1-ε1则输出最佳分解数k
opt
＝k；否则k＝k+1，转入步骤2.2。
[0031]
在步骤3中，针对不同分量序列的特点，利用最大信息系数筛选与预测分量关系密切的输入特征，并使用近似熵依据各分量序列周期及复杂度的差异选用不同的预测模型；其中，低频子序列分量采用gru预测模型；高频子序列分量则采用densenet预测模型；
[0032]
在步骤3中，具体包括以下子步骤：
[0033]
步骤3.1：针对不同分量序列的特点，利用最大信息系数筛选与预测分量关系密切的输入特征；
[0034]
步骤3.2：依据各分量序列周期及复杂度将子序列分为高频子序列和低频子序列；
[0035]
步骤3.3：低频子序列分量使用gru预测模型预测；高频子序列分量则使用densenet预测模型预测。
[0036]
在步骤3.1中，具体包括以下步骤：
[0037]
步骤3.1.1：设负荷功率的影响因素集合为e＝{eq}(q＝1,2,
…
,q)，q为全部影响因素个数。eq与uk之间的mic值计算如下：
[0038][0039][0040]
式中，mi(,)为互信息，p(eq,uk)是eq和uk的联合概率密度分布函数，p()表示边缘概率密度，r*s为网格尺度，b通常取样本总数的0.6次方；
[0041]
步骤3.1.2：设定一个阈值mic
opt
，当f
mic
(eq,uk)》mic
opt
时，则影响因素eq为uk的相关外部影响因素；将子序列分量uk与筛选得到的影响因素集合ek共同作为预测模型的输入变量xk＝[uk,ek](k＝1,2,
…
,k
opt
+1)。
[0042]
在步骤3.2中，具体包括以下步骤：
[0043]
步骤3.2.1：选取合适的嵌入维数m将该时间序列重构为m维矢量：
[0044][0045]
其中，sk(t)为重构后的m维矢量；
[0046]
步骤3.2.2：计算sk(i)和sk(j)的距离dm[sk(i),sk(j)]:
[0047][0048]
其中，i＝1,2,
…
,t-m+1，j＝1,2,
…
,t-m+1。
[0049]
步骤3.2.3：给定相似容限r，定义其中num代表数目，求的自然对数并取平均值，记作φm(r)：
[0050][0051]
步骤3.2.4：维数增加为m+1，重复1)-3)得到φ
m+1
(r)，近似熵表示为：
[0052]apen
(m,r)＝φm(r)-φ
m+1
(r)
ꢀꢀꢀꢀꢀꢀꢀ
(13)
[0053]
通常m＝2，r为0.15～0.2s，s为时间序列的标准差。
[0054]
通过上述方法计算出各子序列uk的周期lk和近似熵a
pen
后，根据设定的阈值t
l
和a，当子序列uk同时满足周期lk＞t
l
且近似熵a
pen
＜a时则为低频子序列，否则为高频子序列。假设按此标准归类为低频子序列的共有k1个，其集合表示为归类为高频子序列的共有k
opt
+1-k1个，其集合表示为其中k
l
和kh分别表示高频和低频子序列下标集合。
[0055]
在步骤3.3中，具体包括以下步骤：
[0056]
步骤3.3.1：低频分量预测模型为：其中每个时刻信息为k∈k
l
，对应预测输出为y
k,t
，f
l
为低频分量预测模型。gru神经网络对于信息的处理公式如下：
[0057]rt
＝σ(wr·
[h
t-1
,x
t
])
ꢀꢀꢀꢀ
(14)
[0058]zt
＝σ(wz·
[h
t-1
,x
t
])
ꢀꢀꢀꢀ
(15)
[0059][0060][0061]
其中，z
t
、r
t
代表更新门和重置门输出，x
t
为当前时刻输入，h
t-1
为上一时刻隐含层信息，wr、wz、代表权重矩阵，σ是sigmoid激活函数,[]为两个向量间的连接；
⊙
为hadamard乘积。
[0062]
步骤3.3.2：输出低频子序列预测结果。
[0063]
步骤3.3.3：高频分量预测模型为：输入数据由uk的ti个历史时刻数据及未来to时刻的相关影响因素ek共同组成：k∈kh；
[0064]
步骤3.3.4：首先密集连接模块从输入数据xk中提取局部特征，第m个密集连接块中第i层输出该层前所有层的信息，其公式如下:
[0065][0066]
其中，表示第k个子序列第m个密集连接块中第i层卷积核权重向量，表示卷积运算符，表示偏移项，f表示非线性函数，包括采用的归一化处理函数bn、relu激活函
数。
[0067]
步骤3.3.5：将m个密集连接模块级联组成densenet网络，如图3所示。在输入的每个密集连接模块前添加卷积层将密集连接块中卷积核通道数统一变换为输入层的通道数。通过m个密集连接块后，得到输出如下：
[0068][0069]
步骤3.3.6：为了避免过拟合，在m-1个的密集连接模块后添加了一个dropout层。最后通过线性全连接层得到预测值：
[0070][0071]
其中，wk表示线性全连接层权重向量，bk表示偏移向量。
[0072]
一种将负荷时间序列进行分解的方法，采用依据相关熵确定最佳分解数的变分模态分解技术，将原始负荷时间序列分解为k个本征模态函数imf子序列分量和1个残差分量；具体包括以下步骤：
[0073]
步骤2.1：设置分解次数k＝1；
[0074]
步骤2.2：输入步骤1中采集的负荷功率时间序列数据集，构造约束变分问题：
[0075][0076]
其中，{uk}是imf集合，{ωk}是中心频率集合，δ(t)是狄拉克分布，表示卷积，uk(t)代表第k个imf子序列分量，t为采样时刻点。
[0077]
步骤2.3：引入二次罚因子α和拉格朗日乘法算子λ，将式(1)转化为无约束变分问题，表示为：
[0078][0079]
其中，α是二次罚因子，λ是拉格朗日乘法算子，《,》表示内积；
[0080]
步骤2.4：采用交替方向乘子法求解式(2)，uk、ωk和λ迭代更新过程为：
[0081][0082][0083]
[0084]
其中，n为迭代次数，f(ω)、分别代表时域下的f(t)、经过傅里叶变换到频域的结果，τ为噪声容忍度；
[0085]
步骤2.5：根据公式(6)、(7)计算imf分量间的相关熵最大值，通过检测imf分量两两间的相关熵是否近似为1，从而判断对应分解数k条件下是否出现模态混叠现象，以此为依据确定最佳分解数k
opt
。
[0086][0087][0088]
步骤2.6：若v》1-ε1则输出最佳分解数k
opt
＝k；否则k＝k+1，转入步骤2.2。
[0089]
一种将子序列划分为高频子序列、低频子序列的方法，它包括以下步骤：
[0090]
步骤1)针对不同分量序列的特点，利用最大信息系数筛选与预测分量关系密切的输入特征；
[0091]
步骤2)依据各分量序列周期及复杂度将子序列分为高频子序列和低频子序列。
[0092]
在步骤1)中，具体包括以下步骤：
[0093]
步骤(1)设负荷功率的影响因素集合为e＝{eq}(q＝1,2,
…
,q)，q为全部影响因素个数。eq与uk之间的mic值计算如下：
[0094][0095][0096]
式中，mi(,)为互信息，p(eq,uk)是eq和uk的联合概率密度分布函数，p()表示边缘概率密度，r*s为网格尺度，b通常取样本总数的0.6次方；
[0097]
步骤(2)设定一个阈值mic
opt
，当f
mic
(eq,uk)》mic
opt
时，则影响因素eq为uk的相关外部影响因素；将子序列分量uk与筛选得到的影响因素集合ek共同作为预测模型的输入变量xk＝[uk,ek](k＝1,2,
…
,k
opt
+1)。
[0098]
在步骤2)中，具体包括以下步骤：
[0099]
步骤(1)选取合适的嵌入维数m将该时间序列重构为m维矢量：
[0100][0101]
其中，sk(t)为重构后的m维矢量；
[0102]
步骤(2)计算sk(i)和sk(j)的距离dm[sk(i),sk(j)]:
[0103][0104]
其中，i＝1,2,
…
,t-m+1，j＝1,2,
…
,t-m+1。
[0105]
步骤(3)给定相似容限r，定义其中num代表数目，求的自然对数并取平均值，记作φm(r)：
[0106][0107]
步骤(4)维数增加为m+1，重复1)-3)得到φ
m+1
(r)，近似熵表示为：
[0108]apen
(m,r)＝φm(r)-φ
m+1
(r)
ꢀꢀꢀꢀ
(13)
[0109]
通常m＝2，r为0.15～0.2s，s为时间序列的标准差。
[0110]
通过上述方法计算出各子序列uk的周期lk和近似熵a
pen
后，根据设定的阈值t
l
和a，当子序列uk同时满足周期lk＞t
l
且近似熵a
pen
＜a时则为低频子序列，否则为高频子序列。假设按此标准归类为低频子序列的共有k1个，其集合表示为归类为高频子序列的共有k
opt
+1-k1个，其集合表示为其中k
l
和kh分别表示高频和低频子序列下标集合。
[0111]
与现有技术相比，本发明具有如下技术效果：
[0112]
1)本发明所提出的一种cvmd-gru-densenet混合模型短期电力负荷预测方法；在分解阶段，提出改进变分模态分解技术依据相关熵确定最佳子序列数量，提高分解质量；在预测阶段，针对各子序列特点分别选择输入特征和预测模型，其中平滑且规律性强的低频子序列分量采用门控循环单元神经网络预测模型，强随机性的高频子序列分量采用密集卷积神经网络预测模型；最后将各分量在未来预测时段的预测结果叠加，重构为预测时段的负荷预测曲线。以湖北某城市四季的实际负荷数据为算例，验证结果表明，所提方法能有效提高短期负荷预测的精度，并且具有较强的泛化能力；
[0113]
2)本发明第一次将相关熵与vmd相结合，为确定vmd最佳分解次数提供了新的方法；
[0114]
2)分解后的负荷数据更有利于预测，可以显著提高预测模型的运行效率和预测精度；
[0115]
3)针对分解后各分量序列周期及复杂度的差异选用不同的预测模型，为负荷预测提供了新思路；
[0116]
4)实际算例结果表明，该方法对不同季节短期电力负荷预测的有效性，可满足实际电力系统短期负荷预测的需求。
附图说明
[0117]
下面结合附图和实施例对本发明作进一步说明：
[0118]
图1为cvmd-gru-densenet的整体结构框架图；
[0119]
图2为gru神经网络结构图；
[0120]
图3为densenet神经网络结构图；
[0121]
图4为实施例中a城市2016年3月-2017年2月日均负荷；
[0122]
图5为各时期不同k时相关熵最大值；
[0123]
图6为时期2的cvmd分解结果；
[0124]
图7为mic相关性分析统计；
[0125]
图8为时期1的cvmd-gru-densenet与单一预测方法对比图；
[0126]
图9为时期2的cvmd-gru-densenet与单一预测方法对比图；
[0127]
图10为时期3的cvmd-gru-densenet与单一预测方法对比图；
[0128]
图11为时期4的cvmd-gru-densenet与单一预测方法对比图；
[0129]
图12为时期1的cvmd-gru-densenet与组合预测方法对比图；
[0130]
图13为时期2的cvmd-gru-densenet与组合预测方法对比图；
[0131]
图14为时期3的cvmd-gru-densenet与组合预测方法对比图；
[0132]
图15为时期4的cvmd-gru-densenet与组合预测方法对比图。
具体实施方式
[0133]
如图1所示，一种基于cvmd-gru-densenet混合模型的短期电力负荷预测方法，包括以下步骤：
[0134]
步骤1：采集某地区以15min为单位的负荷功率时间序列数据集；
[0135]
步骤2：采用依据相关熵确定最佳分解数的变分模态分解技术(correntropy variational mode decomposition,cvmd)，将原始负荷时间序列分解为k个本征模态函数(intrinsic mode function,imf)子序列分量和1个残差分量(residual,res)。
[0136]
步骤2.1：设置分解次数k＝1；
[0137]
步骤2.2：输入步骤1中采集的负荷功率时间序列数据集，构造约束变分问题：
[0138][0139]
其中，{uk}是imf集合，{ωk}是中心频率集合，δ(t)是狄拉克分布，表示卷积，uk(t)代表第k个imf子序列分量，t为采样时刻点。
[0140]
步骤2.3：引入二次罚因子α和拉格朗日乘法算子λ，将式(1)转化为无约束变分问题，表示为：
[0141][0142]
其中，α是二次罚因子，λ是拉格朗日乘法算子，《,》表示内积；
[0143]
步骤2.4：采用交替方向乘子法求解式(2)，uk、ωk和λ迭代更新过程为：
[0144][0145][0146]
[0147]
其中，n为迭代次数，f(ω)、分别代表时域下的f(t)、经过傅里叶变换到频域的结果，τ为噪声容忍度。
[0148]
步骤2.5：根据公式(6)、(7)计算imf分量间的相关熵最大值，通过检测imf分量两两间的相关熵是否近似为1，从而判断对应分解数k条件下是否出现模态混叠现象，以此为依据确定最佳分解数k
opt
。
[0149][0150][0151]
步骤2.6：若v》1-ε1则输出最佳分解数k
opt
＝k；否则k＝k+1，转入步骤2.2。
[0152]
步骤3：针对不同分量序列的特点，利用最大信息系数(maximal information coefficient,mic)筛选与预测分量关系密切的输入特征，并使用近似熵(approximate entropy,apen)依据各分量序列周期及复杂度的差异选用不同的预测模型。其中，低频子序列分量采用gru预测模型，它可快速预测子序列的整体时序变化趋势；高频子序列分量则采用densenet预测模型，它更善于学习子序列复杂的非线性局部行为；
[0153]
步骤3.1：针对不同分量序列的特点，利用最大信息系数(maximal information coefficient,mic)筛选与预测分量关系密切的输入特征。
[0154]
步骤3.1.1：设负荷功率的影响因素集合为e＝{eq}(q＝1,2,
…
,q)，q为全部影响因素个数。eq与uk之间的mic值计算如下：
[0155][0156][0157]
式中，mi(,)为互信息，p(eq,uk)是eq和uk的联合概率密度分布函数，p()表示边缘概率密度，r*s为网格尺度，b是变量，网格规模需要小于b，通常取样本总数的0.6次方；。
[0158]
步骤3.1.2：设定一个阈值mic
opt
，当f
mic
(eq,uk)》mic
opt
时，则影响因素eq为uk的相关外部影响因素。将子序列分量uk与筛选得到的影响因素集合ek共同作为预测模型的输入变量xk＝[uk,ek](k＝1,2,
…
,k
opt
+1)。
[0159]
步骤3.2：依据各分量序列周期及复杂度将子序列分为高频子序列和低频子序列。
[0160]
步骤3.2.1：选取合适的嵌入维数m将该时间序列重构为m维矢量：
[0161][0162]
其中，sk(t)为重构后的m维矢量；
[0163]
步骤3.2.2：计算sk(i)和sk(j)的距离dm[sk(i),sk(j)]:
[0164][0165]
其中，i＝1,2,
…
,t-m+1，j＝1,2,
…
,t-m+1。
[0166]
步骤3.2.3：给定相似容限r，定义其中num代表数目，求的自然对数并取平均值，记作φm(r)：
[0167][0168]
步骤3.2.4：维数增加为m+1，重复1)-3)得到φ
m+1
(r)，近似熵表示为：
[0169]apen
(m,r)＝φm(r)-φ
m+1
(r)
ꢀꢀꢀꢀꢀ
(13)
[0170]
通常m＝2，r为0.15～0.2s，s为时间序列的标准差
[20]
。
[0171]
通过上述方法计算出各子序列uk的周期lk和近似熵a
pen
后，根据设定的阈值t
l
和a，当子序列uk同时满足周期lk＞t
l
且近似熵a
pen
＜a时则为低频子序列，否则为高频子序列。假设按此标准归类为低频子序列的共有k1个，其集合表示为归类为高频子序列的共有k
opt
+1-k1个，其集合表示为其中k
l
和kh分别表示高频和低频子序列下标集合。
[0172]
步骤3.3：低频子序列分量采用gru预测模型预测；高频子序列分量则采用densenet预测模型预测。
[0173]
步骤3.3.1：低频分量预测模型为：其中每个时刻信息为k∈k
l
，对应预测输出为y
k,t
，f
l
为低频分量预测模型。gru神经网络对于信息的处理公式如下：
[0174]rt
＝σ(wr·
[h
t-1
,x
t
])
ꢀꢀꢀꢀ
(14)
[0175]zt
＝σ(wz·
[h
t-1
,x
t
])
ꢀꢀꢀꢀ
(15)
[0176][0177][0178]
其中，z
t
、r
t
代表更新门和重置门输出，x
t
为当前时刻输入，h
t-1
为上一时刻隐含层信息，wr、wz、代表权重矩阵，σ是sigmoid激活函数,[]为两个向量间的连接；
⊙
为hadamard乘积。
[0179]
步骤3.3.2：输出低频子序列预测结果。
[0180]
步骤3.3.3：高频分量预测模型为：输入数据由uk的ti个历史时刻数据及未来to时刻的相关影响因素ek共同组成：k∈kh；
[0181]
步骤3.3.4：首先密集连接模块从输入数据xk中提取局部特征，第m个密集连接块中第i层输出该层前所有层的信息，其公式如下:
[0182][0183]
其中，表示第k个子序列第m个密集连接块中第i层卷积核权重向量，表示卷积运算符，表示偏移项，f表示非线性函数，包括采用的归一化处理函数bn、relu激活函数。
[0184]
步骤3.3.5：将m个密集连接模块级联组成densenet网络，如图3所示。在输入的每个密集连接模块前添加卷积层将密集连接块中卷积核通道数统一变换为输入层的通道数。
通过m个密集连接块后，得到输出如下：
[0185][0186]
步骤3.3.6：为了避免过拟合，在m-1个的密集连接模块后添加了一个dropout层。最后通过线性全连接层得到预测值：
[0187][0188]
其中，wk表示线性全连接层权重向量，bk表示偏移向量。
[0189]
步骤4：最后将各子序列分量的预测结果叠加，重构得到未来日时序负荷预测结果。
[0190]
步骤4.1：将步骤3中得到的预测数据反归一化；
[0191]
步骤4.2：将步骤4.1的结果叠加，重构得到未来日时序负荷预测结果。
[0192]
在本实施例中，某城市2016年3月-2017年2月15min采样频率的实际负荷数据作为算例，验证所提方法的有效性。由于图4所示的日平均负荷曲线显示出显著的负荷季节性特点，为了更好地评价预测性能，将负荷样本划分为春夏秋冬四个季节时期，测试集分别为四个时期的最后1周。
[0193]
利用cvmd对四个时期负荷数据进行分解，不同k值下的相关熵最大值如图5所示，ε1取阈值0.02时，四个时期对应的最佳分解次数k
opt
分别为6、5、5、5。以时期2的分解结果为例，分解后的各分量如图6所示。将气温(日均气温、日最高气温、日最低气温)、降雨量以及日类型(节假日、星期类型)与各负荷子序列分量进行mic相关性分析，柱状图见图7。图8为四个时期测试集实际值和预测值的散点图，横坐标为预测值，纵坐标为实际负荷。为验证本模型的精确性选取svm、cnn、lstm、gru和densenet五种常见的单一模型，以及emd-gru-densenet、vmd-gru-densenet两种组合模型和本文模型进行对比研究，图8-图5给出了四个时期测试集中周三的日预测曲线，其中图8-图11为四个时期本文模型与单一预测模型对比，图12-图15为四个时期本文模型与组合预测模型对比。相较于常见单一的预测模型，本文所提组合预测模型通过对负荷数据合理分解，并针对分解后各子序列的特点选择合适的输入变量及预测模型，可使mape预测误差降低0.63％～1.81％。相较于基于emd和传统vmd分解技术的组合预测模型，本文通过相关熵确定vmd最佳分解次数k
opt
可提高分解质量，使mape预测误差降低0.31％～1.33％，并且对不同季节的负荷预测具有较好的泛化性。夏季负荷波动剧烈，预测难度高于其它季节，mape误差比其它季节高约0.68％～0.99％。

技术特征：

1.一种基于cvmd-gru-densenet混合模型的短期电力负荷预测方法，其特征在于，它包括以下步骤：步骤1：采集目标地区一定时间段为单位的负荷功率时间序列数据集；步骤2：将负荷功率时间序列分解为k个本征模态函数分量和1个残差分量；k为：分解后本征模态函数分量的个数，为大于0的整数步骤3：筛选与预测分量关系密切的输入特征，依据各分量序列周期及复杂度的差异选用不同的预测模型；步骤4：将各子序列分量的预测结果叠加，重构得到未来日时序负荷预测结果。2.根据权利要求1所述的方法，其特征在于，在步骤2中，采用依据相关熵确定最佳分解数的变分模态分解技术，将原始负荷时间序列分解为k个本征模态函数imf子序列分量和1个残差分量；具体包括以下步骤：步骤2.1：设置分解次数k；步骤2.2：输入步骤1中采集的负荷功率时间序列数据集，构造约束变分问题：其中，{u
k
}是imf集合，{ω
k
}是中心频率集合，δ(t)是狄拉克分布，表示卷积，u
k
(t)代表第k个imf子序列分量，t为采样时刻点；步骤2.3：引入二次罚因子α和拉格朗日乘法算子λ，将式(1)转化为无约束变分问题，表示为：其中，α是二次罚因子，λ是拉格朗日乘法算子，<,>表示内积；步骤2.4：采用交替方向乘子法求解式(2)，u
k
、ω
k
和λ迭代更新过程为：和λ迭代更新过程为：和λ迭代更新过程为：其中，n为迭代次数，f(ω)、分别代表时域下的f(t)、经过傅里叶变换到频域的结果，τ为噪声容忍度；
步骤2.5：根据公式(6)、(7)计算imf分量间的相关熵最大值，通过检测imf分量两两间的相关熵是否近似为1，从而判断对应分解数k条件下是否出现模态混叠现象，以此为依据确定最佳分解数k
opt
；；步骤2.6：若v>1-ε1则输出最佳分解数k
opt
＝k；否则k＝k+1，转入步骤2.2。3.根据权利要求1所述的方法，其特征在于，在步骤3中，针对不同分量序列的特点，利用最大信息系数筛选与预测分量关系密切的输入特征，并使用近似熵依据各分量序列周期及复杂度的差异选用不同的预测模型；其中，低频子序列分量采用gru预测模型；高频子序列分量则采用densenet预测模型；在步骤3中，具体包括以下子步骤：步骤3.1：针对不同分量序列的特点，利用最大信息系数筛选与预测分量关系密切的输入特征；步骤3.2：依据各分量序列周期及复杂度将子序列分为高频子序列和低频子序列；步骤3.3：低频子序列分量使用gru预测模型预测；高频子序列分量则使用densenet预测模型预测。4.根据权利要求3所述的方法，其特征在于，在步骤3.1中，具体包括以下步骤：步骤3.1.1：设负荷功率的影响因素集合为e＝{e
q
}(q＝1,2,
…
,q)，q为全部影响因素个数，e
q
与u
k
之间的mic值计算如下：之间的mic值计算如下：式中，mi(,)为互信息，p(e
q
,u
k
)是e
q
和u
k
的联合概率密度分布函数，p()表示边缘概率密度，r*s为网格尺度，b是变量，网格规模需要小于b；步骤3.1.2：设定一个阈值mic
opt
，当f
mic
(e
q
,u
k
)>mic
opt
时，则影响因素e
q
为u
k
的相关外部影响因素；将子序列分量u
k
与筛选得到的影响因素集合e
k
共同作为预测模型的输入变量x
k
＝[u
k
,e
k
](k＝1,2,
…
,k
opt
+1)。5.根据权利要求3所述的方法，其特征在于，在步骤3.2中，具体包括以下步骤：步骤3.2.1：选取合适的嵌入维数m将该时间序列重构为m维矢量：其中，s
k
(t)为重构后的m维矢量；步骤3.2.2：计算s
k
(i)和s
k
(j)的距离d
m
[s
k
(i),s
k
(j)]:
其中，i＝1,2,
…
,t-m+1，j＝1,2,
…
,t-m+1；步骤3.2.3：给定相似容限r，定义其中num代表数目，求的自然对数并取平均值，记作φ
m
(r)：步骤3.2.4：维数增加为m+1，重复步骤3.2.1-步骤3.2.3得到φ
m+1
(r)，近似熵表示为：a
pen
(m,r)＝φ
m
(r)-φ
m+1
(r)
ꢀꢀꢀꢀꢀꢀꢀ
(13)通常m＝2，r为0.15～0.2s，s为时间序列的标准差；通过上述方法计算出各子序列u
k
的周期l
k
和近似熵a
pen
后，根据设定的阈值t
l
和a，当子序列u
k
同时满足周期l
k
＞t
l
且近似熵a
pen
＜a时则为低频子序列，否则为高频子序列；假设按此标准归类为低频子序列的共有k1个，其集合表示为归类为高频子序列的共有k
opt
+1-k1个，其集合表示为其中k
l
和k
h
分别表示高频和低频子序列下标集合。6.根据权利要求3所述的方法，其特征在于，在步骤3.3中，具体包括以下步骤：步骤3.3.1：低频分量预测模型为：其中每个时刻信息为对应预测输出为y
k,t
，f
l
为低频分量预测模型；gru神经网络对于信息的处理公式如下：r
t
＝σ(w
r
·
[h
t-1
,x
t
])
ꢀꢀꢀꢀꢀꢀꢀꢀ
(14)z
t
＝σ(w
z
·
[h
t-1
,x
t
])
ꢀꢀꢀꢀꢀꢀꢀ
(15)(15)其中，z
t
、r
t
代表更新门和重置门输出，x
t
为当前时刻输入，h
t-1
为上一时刻隐含层信息，w
r
、w
z
、代表权重矩阵，σ是sigmoid激活函数,[]为两个向量间的连接；
⊙
为hadamard乘积；步骤3.3.2：输出低频子序列预测结果；步骤3.3.3：高频分量预测模型为：输入数据由u
k
的t
i
个历史时刻数据及未来t
o
时刻的相关影响因素e
k
共同组成：步骤3.3.4：首先密集连接模块从输入数据x
k
中提取局部特征，第m个密集连接块中第i层输出该层前所有层的信息，其公式如下:其中，表示第k个子序列第m个密集连接块中第i层卷积核权重向量，表示卷积运算符，表示偏移项，f表示非线性函数，包括采用的归一化处理函数bn、relu激活函数；步骤3.3.5：将m个密集连接模块级联组成densenet网络；在输入的每个密集连接模块前添加卷积层将密集连接块中卷积核通道数统一变换为输入层的通道数；通过m个密集连
接块后，得到输出如下：步骤3.3.6：为了避免过拟合，在m-1个的密集连接模块后添加了一个dropout层；最后通过线性全连接层得到预测值：其中，w
k
表示线性全连接层权重向量，b
k
表示偏移向量。7.一种将负荷时间序列进行分解的方法，其特征在于，采用依据相关熵确定最佳分解数的变分模态分解技术，将原始负荷时间序列分解为k个本征模态函数imf子序列分量和1个残差分量；具体包括以下步骤：步骤2.1：设置分解次数k＝1；步骤2.2：输入步骤1中采集的负荷功率时间序列数据集，构造约束变分问题：其中，{u
k
}是imf集合，{ω
k
}是中心频率集合，δ(t)是狄拉克分布，表示卷积，u
k
(t)代表第k个imf子序列分量，t为采样时刻点；步骤2.3：引入二次罚因子α和拉格朗日乘法算子λ，将式(1)转化为无约束变分问题，表示为：其中，α是二次罚因子，λ是拉格朗日乘法算子，<,>表示内积；步骤2.4：采用交替方向乘子法求解式(2)，u
k
、ω
k
和λ迭代更新过程为：和λ迭代更新过程为：和λ迭代更新过程为：其中，n为迭代次数，f(ω)、分别代表时域下的f(t)、经过傅里叶变换到频域的结果，τ为噪声容忍度；步骤2.5：根据公式(6)、(7)计算imf分量间的相关熵最大值，通过检测imf分量两两间的相关熵是否近似为1，从而判断对应分解数k条件下是否出现模态混叠现象，以此为依据
确定最佳分解数k
opt
；；步骤2.6：若v>1-ε1则输出最佳分解数k
opt
＝k；否则k＝k+1，转入步骤2.2。8.一种将子序列划分为高频子序列、低频子序列的方法，其特征在于，它包括以下步骤：步骤1)针对不同分量序列的特点，利用最大信息系数筛选与预测分量关系密切的输入特征；步骤2)依据各分量序列周期及复杂度将子序列分为高频子序列和低频子序列。9.根据权利要求8所述的方法，其特征在于，在步骤1)中，具体包括以下步骤：步骤(1)设负荷功率的影响因素集合为e＝{e
q
}(q＝1,2,
…
,q)，q为全部影响因素个数；e
q
与u
k
之间的mic值计算如下：之间的mic值计算如下：式中，mi(,)为互信息，p(e
q
,u
k
)是e
q
和u
k
的联合概率密度分布函数，p()表示边缘概率密度，r*s为网格尺度，b通常取样本总数的0.6次方；步骤(2)设定一个阈值mic
opt
，当f
mic
(e
q
,u
k
)>mic
opt
时，则影响因素e
q
为u
k
的相关外部影响因素；将子序列分量u
k
与筛选得到的影响因素集合e
k
共同作为预测模型的输入变量x
k
＝[u
k
,e
k
](k＝1,2,
…
,k
opt
+1)。10.根据权利要求8所述的方法，其特征在于，在步骤2)中，具体包括以下步骤：步骤(1)选取合适的嵌入维数m将该时间序列重构为m维矢量：其中，s
k
(t)为重构后的m维矢量；步骤(2)计算s
k
(i)和s
k
(j)的距离d
m
[s
k
(i),s
k
(j)]:其中，i＝1,2,
…
,t-m+1，j＝1,2,
…
,t-m+1；步骤(3)给定相似容限r，定义其中num代表数目，求的自然对数并取平均值，记作φ
m
(r)：步骤(4)维数增加为m+1，重复1)-3)得到φ
m+1
(r)，近似熵表示为：a
pen
(m,r)＝φ
m
(r)-φ
m+1
(r)
ꢀꢀꢀꢀꢀꢀꢀꢀ
(13)
通常m＝2，r为0.15～0.2s，s为时间序列的标准差；通过上述方法计算出各子序列u
k
的周期l
k
和近似熵a
pen
后，根据设定的阈值t
l
和a，当子序列u
k
同时满足周期l
k
＞t
l
且近似熵a
pen
<a时则为低频子序列，否则为高频子序列，假设按此标准归类为低频子序列的共有k1个，其集合表示为归类为高频子序列的共有k
opt
+1-k1个，其集合表示为其中k
l
和k
h
分别表示高频和低频子序列下标集合。

技术总结

一种基于CVMD-GRU-DenseNet混合模型的短期电力负荷预测方法，它包括以下步骤：步骤1：采集目标地区一定时间段为单位的负荷功率时间序列数据集；步骤2：将负荷功率时间序列分解为K个本征模态函数分量和1个残差分量；步骤3：筛选与预测分量关系密切的输入特征，依据各分量序列周期及复杂度的差异选用不同的预测模型；步骤4：将各子序列分量的预测结果叠加，重构得到未来日时序负荷预测结果。本发明的目的是为了解决现有技术存在的利用变分模态分解VMD对负荷进行分解难以确定最佳分解次数的技术问题。术问题。术问题。