一种新能源出力引起的反调峰风险的电网峰谷时段划分方法与流程

本发明涉及一种电网峰谷时段划分方法。

背景技术：

国内能源需求的急剧增长，使化石能源的消耗急剧增加，但我国化石能源的产量相对不足。而且大量使用化石能源导致了非常严重的环境问题。所以大力发展可再生的新能源成为应对能源不足问题的主要途径。但是持续增加的新能源负荷峰谷差，间歇性和大波动性，使电网出现反调峰现象增多，其可预测性降低。这使新能源消纳难，并网难，并引起大量的弃风弃光现象。为了能够大量消纳新能源，必须抑制反调峰现象的发生。反调峰现象多发生在峰荷区和谷荷区，因此对峰荷区和谷荷区进行合理的划分对抑制反调峰的发生有重要意义。

目前对峰荷区和谷荷区的划分都是为了适应计算分时电价的需要，利用需求侧实现“削峰填谷”的目的。但这种划分方式并没有考虑到新能源出力特性，没有从源头上解决反调峰现象的发生。

技术实现要素：

针对现有技术以上存在的问题，本发明提出一种基于新能源出力所引起反调峰风险的电网峰谷时段划分方法。本发明以新能源出力所引起的反调峰风险为标准对峰谷时段进行合理的划分，通过对新能源出力特性的分析以抑制反调峰现象的发生，减少弃风弃光现象的发生，从电源端实现“削峰填谷”的目的。

本发明的技术方案为：

步骤1：利用模糊半梯形隶属度函数法确定负荷曲线上各点的峰、谷隶属度。

1.1利用新疆电网典型负荷日的负荷数据，把一天24小时按15分钟为一个时间段，共划分为96个时间段。时段集合为t＝{t1t2,...,t96}，对应的负荷集合为绘制负荷曲线图。其中，集合t中的元素表示每个时间段的起点时刻；集合q中的元素为对应时刻的负荷值。

1.2一天的负荷曲线中必有一个最低谷点a，和一个最高峰点b，采用以下两个原则，并根据最高点和最低点进行峰谷时段划分：

①最低谷点a在谷时段的概率为100％，在峰时段的概率为0；

②最高峰点b在谷时段的概率为0，在峰时段的概率为100％。

1.3用半梯形隶属度函数确定负荷曲线上其余各点处于峰时段和谷时段的可能性，由于该方法只需确定负荷曲线上其余各点相对于最高峰点与最低谷点的距离大小，与具体数值无关，因此可以只采用隶属度函数的线性部分。采用偏大型半梯形隶属度函数计算时点ti的峰隶属度采用偏小型隶属度函数计算时点ti的谷隶属度

式中：为ti时刻的负荷值，b为负荷值中的最大值，a为负荷值中的最小值。

时点的隶属度属性通常指的是该时点至下一时点之间的峰、谷隶属度属性问题，即“以点盖段”。

步骤2：利用半梯形隶属度函数得出负荷曲线上各点的峰、谷隶属度，通过对新能源出力特性的分析得到反调峰风险指标。

该反调峰指标由负荷与新能源和传统能源合成出力之差以及负荷波动率与合成出力波动率之差确定。设典型负荷日当天的新能源和传统能源合成出力集合为：则

式中，为i时点的负荷值，为i时点的新能源和传统能源合成出力，表示i时点的反调峰风险指标；表示负荷与新能源和传统能源合成出力的差值；表示i时点的前一个时点的负荷与新能源和传统能源合成出力的差值；为负荷波动率与合成出力波动率的差值，i＝1,2,…,95，上式中的值采用的值表示。

步骤3：对时点进行模糊聚类，确定各时点所在的时段。

3.1以各时刻负荷数据为分类对象，各时刻的峰、谷隶属度和反调峰风险为统计指标：

则可得到日负荷曲线上各时点组成的时刻点集合的特征指标矩阵x：

该特征指标矩阵中，分别为i时点的峰、谷隶属度指标，i＝1,2,…,96；表示i时点的反调峰风险指标,i＝1,2,…,96。

式中，xi为i时点的统计指标向量；x表示日负荷曲线上各时点组成的时刻点集合的特征指标矩阵x；和分别表示i时点的峰、谷隶属度；表示i时刻的负荷与新能源和传统能源合成出力之差；表示i时刻负荷波动率与合成出力波动率之差。

3.2一般不同的数据有不同的量纲，为了使不同量纲的量也能进行比较，需要对数据作适当的变换。采用平移—标准差变化的方法对矩阵x进行处理，具体方法如下：

式中，xi表示i时点的统计指标向量，为统计向量xi的均值，s为统计向量xi的均方根值；x′i表示数据变换后i时点的统计指标向量。

3.3依照传统聚类方法中的绝对值减数法确定相似系数，建立模糊相似矩阵来表示x′i和x'j之间的关系。

式中，i＝1,2,...,96，j＝1,2,...,96，rij表示模糊相似矩阵中的元素；x′ik和x'jk分别表示经数据变换后i时点与j时点的第k个指标，其中i时点与j时点表示两个不相同的时点；c为约束rij的约束变量；k表示指标的序号。

为了保证rij∈[0,1]，在聚类过程中将c作为约束变量得到模糊相似矩阵r。

3.4采用模糊传递闭包法进行聚类，对于模糊相似矩阵r，依次求其二次方，即r²,r⁴,r⁶,...,r²ⁱ,...，直到第一次出现r²ⁿ＝rⁿ*rⁿ＝rⁿ时，rⁿ便为一个传递闭包t(r)。在传统闭包t(r)＝tij中，0≤tij≤1,i,j＝1,2,…,96。从高到低依次给定传递闭包的阈值λ，求取t(r)的λ—截矩阵rλ：

式中，rλ为λ—截矩阵，tij(λ)为截矩阵中的元素；tij为传统闭包中的元素。

将rλ看成由行向量组成的矩阵，即yi是由时点i的指标组成的行向量，只有yi＝yj(i≠j),(i,j＝1,2,...,96)时，yi与yj才为一类，与峰隶属度为1的时点构成的一类为峰时段，与谷隶属度为1的时点构成的一类为谷时段，其余的各类归为各个平时段。随着λ大小变化便可得到动态聚类，λ越大分类数目就越多，各类中的元素较少即表示峰谷时段持续时间短；λ越小分类数目就越少，各类中的元素较多即表示峰谷时段的持续时间长，因此λ的大小可以反映峰谷时段的持续时间。由于分类时考虑了新能源出力引起反调峰风险和负荷的峰谷隶属度指标，因此对最佳阈值λ的大小可以用类内与类间中这三个指标的距离来确定。因为峰谷时段持续时间的大小充分考虑了反调峰风险，所以划分后的峰谷时段其反调峰风险也是相近的，这为抑制反调峰的研究提供了帮助。

步骤4：构建统计量f表示类间距离与类内距离的比值，由统计量f确定最佳阈值λ，取统计量f最大时的λ值所对应的分类为最佳分类，即为划分的峰谷时段。

设其中为所有时点第k个指标的中心向量组成的行向量，为所有时点的中心向量。对应λ的分类数为r，第i类样本记为：第i类的聚类中心向量为其中为第i类第k个特征的平均值，即：

建立统计量f，

其中，

式中，x为行向量组成矩阵rλ中的元素，表示第i类样本中的第k个特征，为第i类样本中第j个样本向量，表示第i类中第j个样本的第k个特征；f为构建的统计量；ni表示第i类中含有的时点数；n为时点总数，n＝96；r为分类数；其中为总体样本第k个指标的中心向量。

f统计量是服从自由度为r-1，n-r的f分布，其分子表征类与类之间的距离，按照聚类思想，类与类之间的差异越大越好，所以分子要尽量大；分母表征类内各元素之间的距离，按照聚类思想类内各元素差异越小越好，所以分母要尽量小。因此统计量f越大，表示类与类之间的差异就越明显，分类效果就越好。若f＞fα(r-1,n-r),(α＝0.05)，其中α为f分布上分位点，则根据方差显著性理论知道类与类之间的差异是显著的，说明分类合理。

附图说明

图1是本发明实施基于新能源出力特性所引起反调峰风险的峰谷时段划分方法的流程图；

图2是偏小型半梯形隶属度函数图；

图3是偏大型半梯形隶属度函数图。

具体实施方式

以下结合附图及具体实施方式进一步说明本发明。

如图1所示，本发明新能源出力引起的电网反调峰风险的峰谷时段划分方法的流程如下：

1、利用模糊半梯形隶属度函数法确定负荷曲线上各点的峰、谷隶属度指标。

峰、谷隶属度指标分别由半梯形偏大型隶属度函数和半梯形偏小型隶属度函数得出，图2和图3为其函数图像。

2、利用半梯形隶属度函数得出负荷曲线上各点的峰、谷隶属度，通过对新能源出力特性的分析。得到反调峰风险指标，该反调峰风险指标通过式(3)，(4)确定。

3、以各时段负荷数据为分类对象，把各时刻的峰、谷隶属度和反调峰风险作为统计指标，利用模糊传递闭包法进行聚类。

一般不同的数据有不同的量纲，为了使不同量纲的量也能进行比较，需要对数据作适当的变换。而且模糊聚类要求数据在[0，1]区间内，采用平移—标准差变化的方法对指标矩阵x进行处理即将数据代入到式(7)和式(8)中，得到符合聚类要求的数据。然后通过式(9)和式(10)所示的传统聚类方法确定相似系数rij，得到模糊相似矩阵r；

对于模糊相似矩阵r，依次求其二次方，即r²,r⁴,r⁶,...,r²ⁱ,...，直到第一次出现r²ⁿ＝rⁿ*rⁿ＝rⁿ时，rⁿ便为一个传递闭包t(r)。得出传递闭包后需要确定阈值λ来得到t(r)的λ—截距阵rλ，由式(11)可以得到rλ，将rλ看成由行向量组成的矩阵，即：只有yi＝yj(i≠j)，i,j＝1,2,...,96时，yi与yj才为一类。与峰隶属度为1的时点构成的一类为峰时段，与谷隶属度为1的时点构成的一类为谷时段，其余的各类归为各个平时段。随着最佳阈值λ大小变化就可以得到动态聚类。

4、不同的λ值对应不同的分类，构建统计量f来确定最佳阈值λ。由式(12)和式(13)确定的统计量f，其分子表征类与类之间的距离，按照聚类思想，类与类之间的差异越大越好，所以分子要尽量大；分母表征类内各元素之间的距离，按照聚类思想类内各元素差异越小越好，所以分母要尽量小。因此统计量f越大，表示类与类之间的差异就越明显，分类效果就越好。由于λ∈[0，1]，可以从高到低给定多个λ值，然后分别计算出统计量f的值，以统计量f最大为目标来确定最佳阈值λ，由此得到最合理的分类。

本发明在划分峰荷区和谷荷区时加入对新能源出力特性的考虑，并对新能源所引起的反调峰风险进行计算，得出更为合理的峰谷时段划分方法。