一种基于小波变换的微能源网典型日选取与工况集构建方法与流程

1.本发明属于综合能源系统与智慧能源领域，具体涉及一种基于小波变换的微能源网典型日选取与工况集构建方法。


背景技术：

2.微能源网有冷热电多种负荷，需统筹协调多能流以达到源网荷储的系统优化。能源微网能降低间歇性电源的不利影响，提高供电可靠性和电能质量，引起了学界的广泛关注。由于微能源网通常需要接入大量的可再生能源，其规划设计需要充分考虑可再生能源的波动性，这导致了微能源网规划设计问题与其运行优化问题的耦合。
3.微能源网的优化规划一般会通过选定典型日来降低规划的复杂度，而典型日特征数据则是无监督地选取典型日的基础。使用合理的典型日选取方法不仅可以减少规划问题的计算复杂度，还可以保留尽可能多的有效信息提高规划的准确性。典型日选取的常见方法有主观选择法，抽样法，线性规划法，聚类法等。其中主观选择法主要根据季度、日期、天气情况、用能习惯等因素选定典型日，主观性较强，当数据较为复杂时难以应用。抽样法依据随机抽样的思路进行典型日的选取，容易出现数据聚集的现象，难以避免数据信息严重丢失的问题。线性规划法通过建立数学模型对典型日选取结果进行最优化求解，可以自主定制优化的目标和约束条件，但核心的优化思路与聚类法差别不大，且该方法在不同数据集上的性能难以保证。聚类法通过对数据进行分类并计算类质心，可以得到具有一定代表性的典型日集合，而且可以保留数据的时序信息。本发明提供一种基于小波变换的微能源网典型日选取与工况集构建方法，可为冷、热、电等多能流微网以及风光储综合能源系统提供一种高效、实用的典型日选取与工况集构建方法，从而指导系统分析规划与统筹调度优化。


技术实现要素：

4.针对现有技术中存在的问题，本发明提供一种基于小波变换的微能源网典型日选取与工况集构建方法。
5.本发明采用以下技术方案实现：
6.一种基于小波变换的微能源网典型日选取与工况集构建方法，包括以下步骤：
7.步骤s1，根据综合能源系统结构特征、冷热电多能流主体差异、主要环境影响因素，选取微能源网关键数据构成典型日的特征向量；
8.步骤s2，通过haar小波变换对自然日特征数据进行分解，提取时频特征并构建自然日初始特征矩阵；
9.步骤s3，考虑微能源网设备个体的容量配置信息，对自然日初始特征矩阵进行加权，得到自然日动态加权特征矩阵；
10.步骤s4，将自然日特征数据样本采用步骤s2、s3方法处理得到自然日动态加权特征矩阵集合，并基于聚类方法选取典型日，典型日的工况构成工况集。
11.在步骤s1中：
12.根据结构特征及冷、热、电多能流生产主体的差异，考虑微能源网的主要环境影响因素，选取特征数据构成典型日的特征向量。
13.其特征数据选取原则为：对微网系统考虑冷、热、电多能流，供给侧考虑出力与能源生产量，输配侧考虑环境与能量耗散，用户侧考虑负荷需求与实际负荷。其特征数据采样时间要求以自然日为单位，尽可能涵盖冷、热、电等多种微网所有能流。通常情况下每种数据的要求采样间隔时间为1小时。若选取n个特征数据，以原始序列的数值或标值直接构造特征向量，每个特征向量的维度为24
×
n维。
14.通常，对微网系统采用电负荷、冷负荷和热负荷3种数据，构成每个自然日的特征向量，其维度为72。
15.在步骤s2中：
16.s21，采用haar小波变换对自然日原始数据进行特征提取和降维重构。对每类数据的原始序列zr进行一次上采样，获得2n(n∈n*)个采样点的序列zu，然后对zu进行n层的小波分解后可得到各个子空间上的尺度分量和细节分量，如式(1)所示。
17.zu→
(c0,d0,d1,...,di,...,d
n-1
)
ꢀꢀꢀ
(1)
18.式中，c0为尺度分量，为细节分量，其中是第i层分解中第2i个母小波的系数。
19.s22,计算尺度分量和细节分量对原始序列能量的绝对贡献值，得到处理后的每类数据的自然日特征向量。
20.以尺度分量和细节分量对原始序列的贡献值组合来表示自然日的特征，其中尺度分量可以描述序列的整体信息，而细节分量则能有效表征细节信息。通过计算尺度分量和细节分量对原始序列能量的绝对贡献值，可以得到每类数据的自然日特征向量，如式(2)-(3)所示。
21.f＝(c
02
,f0,f1,...,fi,...,f
n-1
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
22.fi＝||di||2ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
23.式中，f为冷、热、电负荷及微能源网特征的自然日特征向量。通过特征提取及处理后，原始序列维度至少可以从接近2
n-1
1维降低至n 1维。
24.s23，归一化处理得到自然日初始特征矩阵。为了保证每类数据在聚类过程中都不受自身量纲的影响，对每类数据的自然日特征向量进行归一化处理，得到自然日初始特征矩阵，以采用微能源网中电负荷、冷负荷和热负荷为特征向量为例，自然日初始特征矩阵如式(4)所示。
[0025][0026]
式中，fr为自然日初始特征矩阵；f
ie
，f
ih
，f
ic
分别为第i个自然日的电负荷、热负荷和冷负荷的归一化特征列向量。
[0027]
在步骤s3中：
[0028]
对式(4)中归一化后的自然日初始特征矩阵进行不同维度的动态加权，以体现特征成分在典型日选取过程中的重要程度。其原则是：
[0029]
当某一能流在系统中起主导作用时，典型日选取过程中应更多地考虑该能流的特征，如对以电供暖、电制冷为主体的综合能源系统，电在冷、热、电多能流中起主要，应赋予其更大的特征数据选取权重。
[0030]
当某一设备的额定功率越大时，典型日选取过程中应更多地考虑该设备对应的特征，因而给对应特征赋予的权值应该更大；如式(5)所示。
[0031][0032]
式中，f为自然日动态加权特征矩阵；和为电、热、冷负荷的最大值或接近实际运行的出力额定值。
[0033]
在步骤s4中：
[0034]
s41，根据自然日特征数据样本，获取自然日动态加权特征矩阵集合。
[0035]
根据步骤s1确定的典型日特征向量采集特征数据，获取待处理的微能源网数据样本。进行数据清洗与预处理后，根据步骤s2、s3所述，将数据样本以自然日为单位进行计算，最终得到自然日动态加权特征矩阵的集合。
[0036]
s42，确定聚类簇数并进行聚类
[0037]
自然日动态加权特征矩阵的集合涵盖了数据样本中的信息，代表每个自然日的特征，需进一步处理。
[0038]
聚类在能源系统典型日选取中具有广泛应用，其核心思想是把n个自然日划分为k个簇，使簇内尽量紧凑，而簇间尽量分散，并以簇的质心作为典型日代表所有簇内的自然日，一般采用均方差作为目标函数，如式(5)所示。
[0039][0040]
式中，k为簇的数量；p为自然日对象空间中的一点；mk为簇μk的质心。其中簇的数量k是需要指定的参数。
[0041]
k聚类以最小化样本与质点平方误差作为目标函数，将每个簇的质点与簇内样本点的平方距离误差和称为畸变程度，那么，对于一个簇，它的畸变程度越低，代表簇内成员越紧密，畸变程度越高，代表簇内结构越松散。畸变程度会随着类别的增加而降低，但对于有一定区分度的数据，在达到某个临界点时畸变程度会得到极大改善，之后缓慢下降，这个临界点就可以考虑为聚类性能较好的点，可采用肘部准则决定k值。
[0042]
s43，根据总量偏差、分布偏差指标，比较典型日选取效果的好坏。
[0043]
为了比较典型日选取效果的好坏，定义总量偏差和分布偏差两个指标进行评估，如式(6)和式(7)所示。
[0044][0045][0046]
式中，δtot为典型日集合的总量偏差；d为典型日集合；ωd和sd为典型日d代表的天数和该典型日的全天聚类数据(电负荷、太阳辐射强度或风速)总量；sa为全年聚类数据总量；δdis为典型日集合的分布偏差；d0为所有自然日的集合；为典型日d在t时刻的聚类数据值；为自然日d在t时刻的聚类数据值。
[0047]
步骤s44，根据步骤s41与s42的计算结果，聚类得出典型日合集。
[0048]
求解完全后，总量偏差越小，则典型日选取效果越好；分布偏差越小，则典型日选取效果越好。由此，可聚类得出典型日合集，可用典型日的工况(即典型日的特征数据所代表的系统状态)指导后续规划优化。
[0049]
本发明的有益效果是：
[0050]
本发明通过haar小波变换对自然日特征数据进行分解，结合聚类方法选取典型日，本发明可为冷、热、电等多能流微网以及风光储综合能源系统提供一种高效、实用的典型日选取与工况集构建方法，从而指导系统分析规划与统筹调度优化。
附图说明
[0051]
图1为本发明的一种基于小波变换的微能源网典型日选取与工况集构建方法；
[0052]
图2为本发明的一种具体实施案例中的社区全年负荷及气象数据曲线图；
[0053]
图3为本发明的一种具体实施案例中的聚类目标函数值随聚类簇数变化曲线图；
具体实施方式
[0054]
下面结合附图和具体实施例对本发明做进一步详细描述，其内容是对本发明的解释而不是限定：
[0055]
本发明公开了一种基于小波变换的微能源网典型日选取与工况集构建方法，包括以下步骤：
[0056]
步骤s1，根据风光储能源生产主体的差异，考虑微能源网的主要环境影响因素，选取数据类型构成典型日的特征向量；步骤s2，通过haar小波变换对自然日特征数据进行分解，提取时频特征并构建自然日初始特征矩阵；步骤s3，考虑微能源网设备个体的容量配置信息，对自然日初始特征矩阵进行加权，得到自然日动态加权特征矩阵；步骤s4，将自然日特征数据样本作为输入获取自然日动态加权特征矩阵集合，并基于聚类方法选取典型日。
[0057]
在步骤s1中：
[0058]
考虑以风-光-储为能源生产主体的微能源网的主要环境影响因素，电负荷为系统的核心参数，太阳辐射强度数据为光电关键参数，风速数据为风电关键参数，可构成选取典
型日的特征向量。
[0059]
要求典型日的选取以自然日为单位，由上采用电负荷、太阳辐射强度和风速3种数据构成每个自然日的特征向量。
[0060]
通常情况下每种数据的采样间隔时间为1小时，若以原始序列的数值或标值直接构造特征向量，每个特征向量的维度为72维。
[0061]
采用美国东北部某社区的实测电负荷数据，以及当地的风速和太阳辐射强度数据，如图2所示。该社区的全年平均电负荷为135.39kw，峰值负荷为379.91kw；全年平均太阳辐射强度为195.90w/m-2，最大太阳辐射强度为1146.22w/m-2；全年平均风速为6.40m/s，最大风速为24.94m/s。
[0062]
在步骤s2中：
[0063]
s21，采用haar小波变换对自然日原始数据进行特征提取和降维重构。对每类数据的原始序列zr进行一次上采样，获得2n(n∈n*)个采样点的序列zu，然后对zu进行n层的小波分解后可得到各个子空间上的尺度分量和细节分量，如下式所示。
[0064]zu
→
(c0,d0,d1,...,di,...,d
n-1
)
[0065]
式中，c0为尺度分量，为细节分量，其中是第i层分解中第2i个母小波的系数。
[0066]
s22,计算尺度分量和细节分量对原始序列能量的绝对贡献值，得到处理后的每类数据的自然日特征向量。
[0067]
以尺度分量和细节分量对原始序列的贡献值组合来表示自然日的特征，其中尺度分量可以描述序列的整体信息，而细节分量则能有效表征细节信息。通过计算尺度分量和细节分量对原始序列能量的绝对贡献值，可以得到每类数据的自然日特征向量，如下式所示。
[0068]
f＝(c
02
,f0,f1,...,fi,...,f
n-1
)
[0069]fi
＝||di||2[0070]
式中，f为电负荷、太阳辐射强度或风速的自然日特征向量。通过特征提取及处理后，原始序列维度至少可以从接近2
n-1
1维降低至n 1维。
[0071]
s23，归一化处理得到自然日初始特征矩阵。为了保证每类数据在聚类过程中都不受自身量纲的影响，对每类数据的自然日特征向量进行归一化处理，得到自然日初始特征矩阵，如下式所示。
[0072][0073][0074][0075][0076]
式中，fr为自然日初始特征矩阵；分别为第i个自然日的电负荷、太阳辐射强度和风速的归一化特征列向量。
[0077]
在步骤s3中：
[0078]
对步骤s23式中归一化后的自然日初始特征矩阵进行不同维度的动态加权，以体现特征成分在典型日选取过程中的重要程度。
[0079]
当某一机组的额定功率越大时，典型日选取过程中应更多地考虑该机组对应的特征，因而给对应特征赋予的权值应该更大，如下式所示。
[0080][0081]
式中，f为自然日动态加权特征矩阵；为电负荷的最大值；pe
pv
和pe
wt
分别为上层模型决定的光伏和风电机组的额定功率，这两个参数的数值接近实际运行过程中光伏和风电机组的出力最大值。以容量配置信息如下表的具体项目为例：
[0082][0083]
在步骤s4中：
[0084]
s41，根据自然日特征数据样本，获取自然日动态加权特征矩阵集合。
[0085]
根据步骤s1确定的典型日特征向量采集特征数据，获取待处理的微能源网数据样本。进行数据清洗与预处理后，根据步骤s2、s3所述，将数据样本以自然日为单位进行计算，最终得到自然日动态加权特征矩阵的集合。
[0086]
s42，采用肘部准则决定聚类簇数。
[0087]
自然日动态加权特征矩阵的集合涵盖了数据样本中的信息，代表每个自然日的特征，需进一步处理。聚类在能源系统典型日选取中具有广泛应用，其核心思想是把n个自然日划分为k个簇，使簇内尽量紧凑，而簇间尽量分散，并以簇的质心作为典型日代表所有簇内的自然日，一般采用均方差作为目标函数，如下式所示。
[0088][0089]
式中，k为簇的数量；p为自然日对象空间中的一点；mk为簇μk的质心。其中簇的数量
k是需要指定的参数。可采用肘部准则决定k值。
[0090]
对聚类算法的目标函数值随着聚类簇数的变化关系进行研究，如图3所示，展示了聚类簇数在2到100的区间内对应的聚类目标函数值。为了平衡计算时间与计算的准确性，根据肘部法则选取聚类簇数为14。
[0091]
s43，根据总量偏差、分布偏差指标，比较典型日选取效果的好坏，聚类得出典型日合集。为了比较典型日选取效果的好坏，定义总量偏差和分布偏差两个指标进行评估，如下式所示。
[0092][0093][0094]
式中，δtot为典型日集合的总量偏差；d为典型日集合；ωd和sd为典型日d代表的天数和该典型日的全天聚类数据(电负荷、太阳辐射强度或风速)总量；sa为全年聚类数据总量；δdis为典型日集合的分布偏差；d0为所有自然日的集合；为典型日d在t时刻的聚类数据值；为自然日d在t时刻的聚类数据值。
[0095]
步骤s44，求解完全后，总量偏差越小，则典型日选取效果越好；分布偏差越小，则典型日选取效果越好。由此，可聚类得出典型日评估结果如下表所示。
[0096][0097]
需要说明的是，以上所述仅为本发明实施方式的一部分，根据本发明所描述的系统所做的等效变化，均包括在本发明的保护范围内。本发明所属技术领域的技术人员可以对所描述的具体实例做类似的方式替代，只要不偏离本发明的结构或者超越本权利要求书所定义的范围，均属于本发明的保护范围。