一种主动配电网日前调度方法、装置、电子设备和存储介质与流程

1.本发明属于主动配电网经济调度技术领域，特别提出一种主动配电网日前调度方法、装置、电子设备和存储介质。


背景技术：

2.主动配电网中蕴含海量的灵活性负荷，对提升电力系统运行效率和经济性具有巨大潜力。然而，灵活性负荷数量大、规模小、种类多，对主动配电网运营商的集中管理和优化调度带来了挑战。为充分挖掘配电系统多种海量灵活性负荷的调度潜力，解决海量灵活性负荷带来的“维度灾”问题，需要制定一种考虑海量灵活性负荷的合理的主动配电网日前调度方法。
3.目前，学术界已经针对主动配电网的优化调度问题已经展开了大量研究，《计及储能和空调负荷的主动配电网多目标优化调度》和《协调储能与柔性负荷的主动配电网多目标优化调度》中分别对空调负荷和柔性负荷进行建模，旨在提升主动配电网的可调性和灵活性。《基于模型预测控制的主动配电网多时间尺度动态优化调度》基于模型预测控制和多步动态滚动优化理论提出了主动配电网的多时间尺度动态优化调度方法，实现了分布式电源、储能及灵活性负荷的协调管理。现有研究和技术主要存在以下不足之处：
4.(1)上述方案中仅对少量灵活性负荷进行逐一建模，然而，在实际系统中，灵活性负荷具有数量大、容量小和参数各异的特点，若对灵活性负荷进行逐一建模，主动配电网运营者将面临“维度灾”难题。因此，主动配电网调度方法需要和有效的集群聚合管理方法相结合，提升系统运营管理效率。
5.(2)上述方案主要关注配电系统内部的优化调度问题。


技术实现要素：

6.本发明的目的是为克服已有技术的不足之处，提出一种主动配电网日前调度方法、装置、电子设备和存储介质。本发明可以实现主动配电网中海量灵活性负荷协调调度，具有计算效率高、兼容性广、实用性强的优点。
7.本公开第一方面实施例提出了一种主动配电网日前调度方法，包括：对主动配电网各节点下的灵活性负荷进行集群负荷聚类分群；根据分群结果，建立主动配电网日前调度优化模型并求解，得到日前调度的预优化结果，根据预优化结果获得各集群的参考功率；对各集群的参考功率进行分解优化，得到集群中各灵活性负荷的实际功率需求；利用所述实际功率需求，完成主动配电网日前调度终优化，得到最终的日前调度优化方案。
8.在本公开的一个实施例中，所述集群负荷聚类分群具体方法如下：
9.获取主动配电网各节点下灵活性负荷的需求参数，根据该需求参数基于njw谱聚类算法对各节点下的灵活性负荷进行集群负荷聚类分群，得到各节点的灵活性负荷集群划分结果，并根据该结果计算每个集群的需求参数。
10.在本公开的一个实施例中，所述计算每个集群的需求参数，具体包括：
[0011]1‑
1)获取各灵活性负荷的需求参数；
[0012]
所述灵活性负荷的需求参数为x
l,j
＝(t
start,l,j
,t
end,l,j
,p
l,j,max
,p
l,j,min
,e
l,j,max
‑
de
,e
l,j,min
‑
de
)，该需求参数满足如下需求模型：
[0013]
p
l,j,min
≤p
l,j,t
≤p
l,j,max
,t
start,l,j
≤t≤t
end,l,j
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0014]
p
l,j,t
＝0,t＜t
start,l,j or t＞t
end,l,j
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
[0015]
e
l,j,t
＝e
l,j,t
‑1 δt
·
p
l,j,t
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
[0016][0017]
其中，p
l,j,min
和p
l,j,max
分别为灵活性负荷j的功率需求下限和上限；e
l,j,min
‑
de
和e
l,j,max
‑
de
分别为灵活性负荷j的能量需求下限和上限；t
start,l,j
和t
end,l,j
分别为灵活性负荷j的需求起始时间和终止时间；δt为调度时间间隔；p
l,j,t
为灵活性负荷j在t时刻的有功功率需求；e
l,j,t
为灵活性负荷j在t时刻的能量状态；
[0018]1‑
2)基于njw谱聚类算法的灵活性负荷集群划分；包括：
[0019]1‑2‑
1)获取各电网节点期望获得的集群数目，其中将节点i处期望获得的集群的数目记为n
agg,i
；
[0020]1‑2‑
2)将各节点下每个灵活性负荷的各需求参数分别进行极值归一化；
[0021]1‑2‑
3)利用极值归一化后的需求参数，构建节点下的各灵活性负荷之间的差异矩阵；
[0022]1‑2‑
4)将差异矩阵转换成邻接矩阵；
[0023]1‑2‑
5)将邻接矩阵归一化，计算归一化后邻接矩阵k前n
agg,i
个最大特征值及其对应的特征向量，利用该特征向量构建特征向量空间矩阵；
[0024]1‑2‑
6)对所述特征向量空间矩阵u进行聚类，得到对应节点下灵活性负荷聚类分群结果；
[0025]1‑
3)根据所述所述灵活性负荷聚类分群结果，计算集群的需求参数包括：
[0026]
节点i处集群k的起始需求时间和终止需求时间和节点i处集群k的功率需求上限和功率需求下限和节点i处集群k的能量需求上限和能量需求下限和各参数表达式分别为：
[0027][0028][0029][0030]
其中，n
l,k,i
为节点i处集群k内部灵活性负荷的数目。
[0031]
在本公开的一个实施例中，所述各集群的参考功率计算方法如下：
[0032]2‑
1)建立主动配电网日前调度优化模型的目标函数，表达式如下：
[0033][0034]
其中，n
s
和ρ
s
分别为场景数目和各场景权重；n
dg
和n
ess
分别为主动配电网中分布式发电机组和储能设备的数量；t为主动配电网运营总时段数；为t时刻场景s下日前能量市场预测电价；为t时刻场景s下旋转备用市场预测电价；为t时刻场景s下主动配电网向能量市场的投标功率；为t时刻场景s下主动配电网向旋转备用市场的投标容量；为分布式发电机组i在t时刻场景s下的有功出力；为分布式发电机组i在t时刻场景s下的备用容量；为分布式发电机组i在t时刻的启停状态；和分别为储能设备i在t时刻场景s下的有功功率输出和有功功率输入；为储能设备i在t时刻场景s下的备用容量；
[0035]
f
idg
(
·
)和f
iess
(
·
)分别为分布式发电机组i和储能设备i的运行成本，表达式如下：
[0036][0037][0038]
其中，为t时刻旋转备用的调用概率，参考取值范围为(0,1)；a
i
，b
i
和c
i
分别为分布式发电机组i的运行成本系数；ap
i
为储能设备i的运行成本系数；
[0039]2‑
2)确定主动配电网日前调度优化模型的约束条件，具体如下：
[0040]
潮流方程模型约束；
[0041][0042][0043][0044][0045][0046][0047]
其中，n
b
为主动配电网中的节点总数；r
ij
和x
ij
分别为节点i和节点j之间线路的电阻和电抗；θ
i,s,t
，v
i,s,t
，p
i,s,t
，q
i,s,t
分别为节点i处t时刻场景s下的相
位、电压幅值、注入有功功率和注入无功功率；和分别为风电场i和光伏电站i在t时刻场景s下的有功出力；和分别为风电场i和光伏电站i在t时刻场景s下的无功出力；为节点i处集群k在t时刻场景s下的参考功率需求；和分别为节点i处t时刻的有功负荷和无功负荷；p
ij,s,t
和q
ij,s,t
分别为节点i和节点j之间线路在t时刻场景s下的有功潮流和无功潮流；
[0048][0049]
其中，s
ij,n
为节点i和节点j之间线路的额定容量；v
i,max
和v
i,min
分别为节点i电压幅值上限和下限；
[0050]
分布式发电机组模型约束；
[0051][0052][0053][0054][0055]
其中，和分别为分布式发电机组i的出力上限和下限；为分布式发电机组i的最大爬坡速率；为分布式发电机组i的额定容量；
[0056]
可再生能源出力约束：
[0057][0058][0059]
其中，和分别为风电场i和光伏电站i在t时刻场景s下的最大预测出力；和分别为风电场i和光伏电站i的额定容量；
[0060]
储能模型约束：
[0061][0062][0063][0064][0065]
其中，分别为储能设备i的最大输入功率和最大输出功率，分别为储能设备i的最小输入功率和最小输出功率；和分别为储能设备i的储能上限和储能下限；η
in
和η
out
分别为功率输入能量转换系数和功率输出能量转换系数；λ为储能的耗散系数；δt为调度时间间隔；为储能设备i的初始能量；
[0066]
系统旋转备用平衡约束：
[0067][0068]
其中，|δr
l
|为主动配电网内部负荷，和分别为t时刻风电场i和光伏电站i的备用需求；n
wp
和n
pv
分别为主动配电网中风电场和光伏电站的数量。
[0069]
灵活性负荷集群模型约束：
[0070][0071][0072][0073]
其中，为节点i处集群k在t时刻场景s下的参考能量状态；
[0074]2‑
3)式(11)
‑
(34)构成的优化调度模型中的决策变量包括(34)构成的优化调度模型中的决策变量包括对该模型求解得到日前调度的预优化结果，并利用各场景权重根据计算各决策变量的总的优化值；
[0075]
根据优化模型的预优化结果，获得各灵活性负荷集群在各场景下每个时刻的参考功率
[0076]
在本公开的一个实施例中，所述集群中各灵活性负荷的实际功率需求计算方法如下：
[0077]3‑
1)各集群根据对应的参考功率分别建立各自的灵活性负荷集群日前分解优化模型，表达式如下：
[0078][0079][0080][0081]
其中，集群中每个灵活性负荷在各时刻的功率满足如式(1)
‑
(4)所示的需求模型；
[0082]
f
1dis
(p
l,j,t
)表示灵活性负荷j所在集群的优化目标1，f
2dis
(p
l,j,t
)表示灵活性负荷j所在集群的优化目标2；m为一正数；
[0083]3‑
2)求解步骤3
‑
1)建立的日前分解优化模型，获得集群中各灵活性负荷的p
l,j,t
的最优解，从而计算得到该灵活性负荷集群的实际功率需求表达式如下：
[0084][0085]
其中，为节点i处集群k在t时刻的实际功率需求。
[0086]
在本公开的一个实施例中，所述主动配电网日前调度终优化，包括：
[0087]
将主动配电网日前调度优化模型中的决策变量替换为此时模型的决策变量更新为再次求解该日前调度优化模型，得到各决策变量的最优解，进而利用得到各决策变量总的最终优化值，从而获取主动配电网的最终日前调度方案。
[0088]
为达上述目的，本公开第二方面实施例提出了一种主动配电网日前调度装置，所述装置位于主动配电网日前调度系统中。所述装置包括：灵活性负荷集群聚类模块，用于对主动配电网各节点下灵活性负荷进行集群负荷聚类分群；主动配电网日前调度预优化模块，用于根据分群结果，建立主动配电网日前调度优化模型并求解，得到日前调度的预优化结果，根据预优化结果获得各集群的参考功率；参考功率分解优化模块，用于对各集群的参考功率进行分解优化，得到集群中各灵活性负荷的实际功率需求；主动配电网日前调度终优化模块，用于利用所述实际功率需求，完成主动配电网日前调度终优化，得到最终的日前调度优化方案。
[0089]
本公开第三方面实施例提出了一种电子设备，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述程序，以用于：对主动配电网各节点下的灵活性负荷进行集群负荷聚类分群；根据分群结果，建立主动配电网日前调度优化模型并求解，得到日前调度的预优化结果，根据预优化结果获得各集群的参考功率；对各集群的参考功率进行分解优化，得到集群中各灵活性负荷的实际功率需求；利用所述实际功率需求，完成主动配电网日前调度终优化，得到最终的日前调度优化方案。
[0090]
本公开第四方面实施例提出了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行，以用于：对主动配电网各节点下的灵活性负荷进行集群负荷聚类分群；根据分群结果，建立主动配电网日前调度优化模型并求解，得到日前调度的预优化结果，根据预优化结果获得各集群的参考功率；对各集群的参考功率进行分解优化，得到集群中各灵活性负荷的实际功率需求；利用所述实际功率需求，完成主动配电网日前调度终优化，得到最终的日前调度优化方案。
[0091]
本发明的特点及有益效果：
[0092]
(1)考虑灵活性负荷聚合和解聚合过程，本发明可以实现主动配电网中海量灵活性负荷协调调度，具有高效的计算效率，能够充分挖掘电力系统负荷侧的灵活性，有效地解决了海量小规模的灵活性负荷给系统优化调度带来的“维度灾”难题；(2)本发明所提的主动配电网优化调度模型兼顾了配电网与上级电力市场的交互，可以实现内部调度计划和外部交互方案的协同优化，在电力市场背景下的主动配电网优化调度领域具有很强的普适性和实用价值。
附图说明
[0093]
图1为本实施例所提供的一种主动配电网日前调度方法的流程示意图；
[0094]
图2为本实施例所提供的一种主动配电网日前调度装置的结构示意图；
[0095]
图3为本实施例所提供的一种主动配电网日前调度方法的硬件结构示意图。
具体实施方式
[0096]
本发明提出一种主动配电网日前调度方法，流程如图1所示，包括以下步骤：
[0097]
1)灵活性负荷集群聚合；
[0098]
主动配电网代理应统计主动配电网每条总线上意愿参与主动配电网聚合运营的灵活性负荷，根据各灵活性负荷的需求参数对各节点下意愿参与主动配电网聚合运营的灵活性负荷进行集群负荷聚类分群，得到各节点的灵活性负荷集群划分结果，并根据该结果计算每个集群的需求参数。
[0099]
2)建立基于灵活性负荷集群的主动配电网日前调度优化模型并求解，得到日前调度的预优化结果，根据预优化结果计算得到各集群的参考功率；
[0100]
根据灵活性负荷集群聚合结果以及可再生能源和各类电价的预测信息，主动配电网代理需要进行预优化，并将预优化调度结果下发至各集群。在预优化步骤结束之后，每个集群都会收到参考功率。
[0101]
3)对各集群建立灵活性负荷集群日前分解优化模型并求解；
[0102]
各灵活性负荷集群将收到的参考功率分解给其内部每个灵活性负荷，通过优化，得到各集群中各灵活性负荷的实际功率需求并进行上报。
[0103]
4)利用步骤3)得到的各集群中各灵活性负荷的实际功率需求，完成基于灵活性负荷集群的主动配电网日前调度终优化；
[0104]
根据所有灵活性负荷的实际功率需求，再次求解基于灵活性负荷集群的主动配电网日前调度优化模型，得到最终的日前调度优化方案，包括：各可再生能源电站、分布式发电机组及储能设备的调度计划和向电力市场投标计划。
[0105]
所述步骤1)具体方法实施例如下：
[0106]
通过灵活性负荷集群，海量灵活性负荷可以给主动配电网代理提供巨大的灵活性，而且不会增加过大的计算负担。本发明提出了一种基于ng
‑
jordan
‑
weiss(njw)谱聚类算法的灵活性负荷集群聚合方法，将具有相似特征的灵活性负荷聚类划分为同一组，形成一个集群统一管理，从而参与主动配电网调度和运营。其中灵活性负荷的需求参数可以基于历史统计数据或概率分布获取。具体步骤如下：
[0107]1‑
1)获取各灵活性负荷的需求参数；
[0108]
记每个灵活性负荷的需求参数为x
l,j
＝(t
start,l,j
,t
end,l,j
,p
l,j,max
,p
l,j,min
,e
l,j,max
‑
de
,e
l,j,min
‑
de
)用于反映其需求信息，灵活性负荷的需求模型可表示如下：
[0109]
p
l,j,min
≤p
l,j,t
≤p
l,j,max
,t
start,l,j
≤t≤t
end,l,j
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0110]
p
l,j,t
＝0,t＜t
start,l,j or t＞t
end,l,j
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
[0111]
e
l,j,t
＝e
l,j,t
‑1 δt
·
p
l,j,t
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
[0112][0113]
其中，p
l,j,min
和p
l,j,max
分别为灵活性负荷j的功率需求下限和上限；e
l,j,min
‑
de
和e
l,j,max
‑
de
分别为灵活性负荷j的能量需求下限和上限。t
start,l,j
和t
end,l,j
分别为灵活性负荷j的需求起始时间和终止时间；δt为调度时间间隔；p
l,j,t
为灵活性负荷j在t时刻的有功功率需求；e
l,j,t
为灵活性负荷j在t时刻的能量状态。
[0114]1‑
2)基于njw谱聚类算法的灵活性负荷集群划分；
[0115]
采用njw谱聚类算法，根据灵活性负荷的特点将其划分到不同的集群。具体步骤如下：
[0116]1‑2‑
1)输入灵活性负荷需求参数和各电网节点期望获得的集群数目，其中将节点i处期望获得的集群的数目记为n
agg,i
。
[0117]1‑2‑
2)将各节点下每个灵活性负荷的各需求参数分贝进行极值归一化；
[0118]
其中，对于节点i下的灵活性负荷j，归一化方法如下：
[0119][0120]
其中，n
l,i
为节点i处灵活性负荷的数目；x
l,j
＝(t
start,l,j
,t
end,l,j
,p
l,j,max
,p
l,j,min
,e
l,j,max
‑
de
,e
l,j,min
‑
de
)为灵活性负荷j归一化前的原始需求参数；x'
l,j
＝(t'
start,l,j
,t'
end,l,j
,p'
l,j,max
,p'
l,j,min
,e'
l,j,max
‑
de
,e'
l,j,min
‑
de
)为灵活性负荷j经过归一化之后的需求参数。
[0121]1‑2‑
3)采用欧式距离来衡量每个节点下不同灵活性负荷之间的差异；
[0122]
对于任一节点i，将该节点下的各灵活性负荷之间的差异矩阵定义为d，其中该矩阵中各元素表达式如下：
[0123][0124]
其中，d(l,j)为矩阵d中第l行第j列的元素，其代表对应节点下灵活性负荷l和j之间的差异；α1，α2，α3，α4，α5和α6为权重系数，其取值范围为(0,1)，在实际实施中，应根据对应参数的标准差确定，其参考值设置为：α1＝α2＝0.2，α3＝α4＝0.5，α5＝α6＝0.8。
[0125]1‑2‑
4)利用高斯核函数将差异矩阵d转换成邻接矩阵k，其中矩阵k中各元素表达式如下：
[0126][0127]
其中，为高斯核函数参数，参考取值范围为[1,2]。
[0128]1‑2‑
5)将邻接矩阵k归一化，计算归一化后矩阵k前n
agg,i
个最大特征值及其对应的特征向量，将得到的特征向量分别记为利用该特征向量构建特征向量空间矩阵
[0129]1‑2‑
6)利用k
‑
means算法对步骤1
‑2‑
5)得到的特征向量空间矩阵u进行聚类，节点i下灵活性负荷聚类分群结果。
[0130]1‑
3)根据步骤1
‑2‑
6)得到的灵活性负荷聚类分群结果，对于每个灵活性负荷集群，根据该集群内部所有灵活性负荷的需求参数，可以获取该集群需求参数如下：
[0131][0132]
[0133][0134]
其中，n
l,k,i
为节点i处集群k内部灵活性负荷的数目；和分别为节点i处集群k的起始需求时间和终止需求时间；和分别为节点i处集群k的功率需求上限和功率需求下限；和分别为节点i处集群k的能量需求上限和能量需求下限。
[0135]
所述步骤2)具体步骤如下：
[0136]2‑
1)建立基于灵活性负荷集群的主动配电网日前调度优化模型的目标函数；
[0137]
本发明适用于作为市场电价接收者的主动配电网运营商。采用基于典型场景集的随机优化方法来对主动配电网运营优化中可再生能源及电价不确定性进行建模，主动配电网日前调度优化模型以主动配电网代理的总运营成本最低为目标，模型中考虑的运营成本应包括以下四项：主动配电网从上层能量市场的电能购买成本，从旋转备用市场的备用购买成本，主动配电网内部分布式发电机组运行成本以及储能运行成本，具体的目标函数如下：
[0138][0139]
其中，n
s
和ρ
s
分别为场景数目和各场景权重；n
dg
和n
ess
分别为主动配电网中中分布式发电机组和储能设备的数量；t为主动配电网运营总时段数；为t时刻场景s下日前能量市场预测电价；为t时刻场景s下旋转备用市场预测电价；为t时刻场景s下主动配电网向能量市场的投标功率；为t时刻场景s下主动配电网向旋转备用市场的投标容量；为分布式发电机组i在t时刻场景s下的有功出力；为分布式发电机组i在t时刻场景s下的备用容量；为分布式发电机组i在t时刻的启停状态；和分别为储能设备i在t时刻场景s下的有功功率输出和有功功率输入；为储能设备i在t时刻场景s下的备用容量。
[0140]
其中，f
idg
(
·
)和f
iess
(
·
)分别为分布式发电机组i和储能设备i的运行成本，表示如下：
[0141][0142][0143]
其中，γ
trm
为t时刻旋转备用的调用概率，参考取值范围为(0,1)；a
i
，b
i
和c
i
分别为分布式发电机组i的运行成本系数；ap
i
为储能设备i的运行成本系数。
[0144]2‑
2)确定基于灵活性负荷集群的主动配电网日前调度优化模型的约束条件，具体如下：
[0145]2‑2‑
1)潮流方程模型约束；
[0146]
本发明采用一种线性化的潮流模型来描述主动配电网内部的有功和无功潮流平
衡如下：
[0147][0148][0149][0150][0151][0152][0153]
其中，n
b
为主动配电网中的节点总数；r
ij
和x
ij
分别为节点i和节点j之间线路的电阻和电抗；θ
i,s,t
，v
i,s,t
，p
i,s,t
，q
i,s,t
分别为节点i处t时刻场景s下的相位、电压幅值、注入有功功率和注入无功功率；和分别为风电场i和光伏电站i在t时刻场景s下的有功出力；和分别为风电场i和光伏电站i在t时刻场景s下的无功出力；为节点i处集群k在t时刻场景s下的参考功率需求；和分别为节点i处t时刻的有功负荷和无功负荷；p
ij,s,t
和q
ij,s,t
分别为节点i和节点j之间线路在t时刻场景s下的有功潮流和无功潮流。
[0154]
此外，各支路潮流和节点电压应限制在一定的范围内：
[0155][0156]
其中，s
ij,n
为节点i和节点j之间线路的额定容量；v
i,max
和v
i,min
分别为节点i电压幅值上限和下限。
[0157]2‑2‑
2)分布式发电机组模型约束；
[0158]
分布式发电机组的模型约束包括机组出力、爬坡出力、启停机时长等：
[0159][0160][0161][0162][0163]
其中，和分别为分布式发电机组i的出力上限和下限；为分布式发电机组i的最大爬坡速率；为分布式发电机组i的额定容量。
[0164]2‑2‑
3)可再生能源出力约束：
[0165]
[0166][0167]
其中，和分别为风电场i和光伏电站i在t时刻场景s下的最大预测出力；和分别为风电场i和光伏电站i的额定容量。
[0168]2‑2‑
4)储能模型约束：
[0169][0170][0171][0172][0173]
其中，分别为储能设备i的最大输入功率和最大输出功率，分别为储能设备i的最小输入功率和最小输出功率；和分别为储能设备i的储能上限和储能下限；η
in
和η
out
分别为功率输入能量转换系数和功率输出能量转换系数；λ为储能的耗散系数；δt为调度时间间隔；为储能设备i的初始能量；。
[0174]2‑2‑
5)系统旋转备用平衡约束：
[0175][0176]
其中，|δr
l
|为主动配电网内部负荷，和分别为t时刻风电场i和光伏电站i的备用需求；n
wp
和n
pv
分别为主动配电网中风电场和光伏电站的数量。
[0177]2‑2‑
6)灵活性负荷集群模型约束：
[0178][0179][0180][0181]
其中，为节点i处集群k在t时刻场景s下的参考能量状态。
[0182]2‑
3)上述式(11)
‑
(34)构成的优化调度模型中的决策变量包括(34)构成的优化调度模型中的决策变量包括对该模型求解得到日前调度的预优化结果，其中，最终的调度决策结果应根据各场景权重计算各决策变量的总的优化值
[0183]
然后根据优化模型的预优化结果，可以获得各灵活性负荷集群在各场景下每个时刻的参考功率并下发至各对应集群。
[0184]
所述步骤3)具体方法如下：
[0185]
各集群根据主动配电网预优化阶段获取的参考功率各灵活性负荷集群需制定其内部灵活性负荷的运行计划。对于某个特定的集群k，有以下两个优化目标，目标1：集
群整体输出功率应尽量追踪预优化阶段获取的参考功率f
1dis
(p
l,j,t
)；目标2：减小每个灵活性负荷的启停次数，促进灵活性负荷稳定运行f
2dis
(p
l,j,t
)；此外，通过设置分解优化的目标函数中系数m的取值可以调整不同优化目标之间的权重。一般来说，分解模型应优先保证集群的功率追踪精度，因此m通常取一个较大数值，参考取值范围为[1000,10000]。
[0186]
对每个集群，分别建立各自的灵活性负荷集群日前分解优化模型，表达式如下：
[0187][0188][0189][0190]
此外每个灵活性负荷在各时刻的功率需求应满足其需求模型，如式(1)
‑
(4)所示；
[0191]
该优化模型为典型的二次规划模型，可直接利用cplex进行求解。
[0192]
求解该优化模型，可以获得集群中各灵活性负荷的调度计划(p
l,j,t
)，以及该灵活性负荷集群的实际功率需求具体计算公式如下：
[0193][0194]
其中，为节点i处集群k在t时刻的实际功率需求。
[0195]
值得注意的是，每个集群的参考功率和实际功率之间可能会存在一定的偏差，因此主动配电网代理需要一个终优化过程来计算最终的日前调度结果。
[0196]
所述步骤4)具体方法如下：
[0197]
终优化与预优化的模型完全相同，不过，在终优化阶段，所有灵活性负荷集群的功率需求已经确定，因此，需将步骤2)的模型式(11)
‑
(34)中的决策变量替换为步骤3)计算得到的此时模型的决策变量更新为此时模型的决策变量更新为再次求解步骤2)的模型，，得到各决策变量的最优解，，进而利用得到各决策变量总的最终优化值，从而获取主动配电网内部其他可控设备(可再生能源电站，分布式发电机组和储能设备)的最终日前调度方案，以及主动配电网向能量市场及辅助服务市场的日前投标方案。
[0198]
本实施例还提出一种主动配电网日前调度装置，结构如图2所示，
[0199]
图2为本实施例所提供的一种主动配电网日前调度装置的结构示意图。
[0200]
其中，主动配电网日前调度装置位于主动配电网日前调度系统中，其中，主动配电网日前调度系统位于电子设备。
[0201]
如图2所示，该主动配电网日前调度装置的实现装置可以包括灵活性负荷集群聚类模块、主动配电网日前调度预优化模块、参考功率分解优化模块和主动配电网日前调度终优化模块；其中：
[0202]
灵活性负荷集群聚类模块，用于对主动配电网各节点下灵活性负荷进行集群负荷
聚类分群；
[0203]
主动配电网日前调度预优化模块，用于根据分群结果，建立基于灵活性负荷集群的主动配电网日前调度优化模型并求解，得到日前调度的预优化结果，根据预优化结果获得各集群的参考功率；
[0204]
参考功率分解优化模块，用于对各集群的参考功率进行分解优化，得到集群中各灵活性负荷的实际功率需求；
[0205]
主动配电网日前调度终优化模块，用于利用所述实际功率需求，完成基于灵活性负荷集群的主动配电网日前调度终优化，得到最终的日前调度优化方案。
[0206]
为了实现上述实施例，本发明实施例还提出一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行，用于执行上述实施例的主动配电网日前调度方法。
[0207]
下面参考图3，其示出了适于用来实现本公开实施例的电子设备900的结构示意图。其中，需要说明的是，该电子设备900中包括主动配电网日前调度系统，其中，本公开实施例中的电子设备可以包括但不限于诸如移动电话、笔记本电脑、数字广播接收器、pda(个人数字助理)、pad(平板电脑)、pmp(便携式多媒体播放器)、车载终端(例如车载导航终端)等等的移动终端以及诸如数字tv、台式计算机、服务器等等的固定终端。图3示出的电子设备仅仅是一个示例，不应对本公开实施例的功能和使用范围带来任何限制。
[0208]
如图3所示，电子设备900可以包括处理装置(例如中央处理器、图形处理器等)901，其可以根据存储在只读存储器(rom)902中的程序或者从存储装置908加载到随机访问存储器(ram)903中的程序而执行各种适当的动作和处理。在ram 903中，还存储有电子设备900操作所需的各种程序和数据。处理装置901、rom 902以及ram 903通过总线904彼此相连。输入/输出(i/o)接口905也连接至总线904。
[0209]
通常，以下装置可以连接至i/o接口905：包括例如触摸屏、触摸板、键盘、鼠标、摄像头、麦克风等的输入装置906；包括例如液晶显示器(lcd)、扬声器、振动器等的输出装置907；包括例如磁带、硬盘等的存储装置908；以及通信装置909。通信装置909可以允许电子设备900与其他设备进行无线或有线通信以交换数据。虽然图3示出了具有各种装置的电子设备900，但是应理解的是，并不要求实施或具备所有示出的装置。可以替代地实施或具备更多或更少的装置。
[0210]
特别地，根据本公开的实施例，上文参考流程图描述的过程可以被实现为计算机软件程序。例如，本实施例包括一种计算机程序产品，其包括承载在计算机可读介质上的计算机程序，该计算机程序包含用于执行流程图中所示方法的程序代码。在这样的实施例中，该计算机程序可以通过通信装置909从网络上被下载和安装，或者从存储装置908被安装，或者从rom 902被安装。在该计算机程序被处理装置901执行时，执行本公开实施例的方法中限定的上述功能。
[0211]
需要说明的是，本公开上述的计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质或者是上述两者的任意组合。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件，或者任意以上的组合。计算机可读存储介质的更具体的例子可以包括但不限于：具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机访问存储器(ram)、只读存储器(rom)、可擦式可编程只读存储器(eprom或闪存)、光纤、便携式紧凑磁盘只读存储器(cd
‑
rom)、光存储器件、磁存储器件、
或者上述的任意合适的组合。在本公开中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。而在本公开中，计算机可读信号介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读信号介质还可以是计算机可读存储介质以外的任何计算机可读介质，该计算机可读信号介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。计算机可读介质上包含的程序代码可以用任何适当的介质传输，包括但不限于：电线、光缆、rf(射频)等等，或者上述的任意合适的组合。
[0212]
上述计算机可读介质可以是上述电子设备中所包含的；也可以是单独存在，而未装配入该电子设备中。
[0213]
上述计算机可读介质承载有一个或者多个程序，当上述一个或者多个程序被该电子设备执行时，使得该电子设备：对主动配电网各节点下的灵活性负荷进行集群负荷聚类分群；根据分群结果，建立基于灵活性负荷集群的主动配电网日前调度优化模型并求解，得到日前调度的预优化结果，根据预优化结果获得各集群的参考功率；对各集群的参考功率进行分解优化，得到集群中各灵活性负荷的实际功率需求；利用所述实际功率需求，完成基于灵活性负荷集群的主动配电网日前调度终优化，得到最终的日前调度优化方案。
[0214]
可以以一种或多种程序设计语言或其组合来编写用于执行本公开的操作的计算机程序代码，上述程序设计语言包括面向对象的程序设计语言—诸如java、smalltalk、c 、python，还包括常规的过程式程序设计语言—诸如“c”语言或类似的程序设计语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络——包括局域网(lan)或广域网(wan)—连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
[0215]
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本技术的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
[0216]
此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本技术的描述中，“多个”的含义是至少两个，例如两个，三个等，除非另有明确具体的限定。
[0217]
流程图中或在此以其他方式描述的任何过程或方法描述可以被理解为，表示包括一个或更多个用于实现特定逻辑功能或过程的步骤的可执行指令的代码的模块、片段或部分，并且本技术的优选实施方式的范围包括另外的实现，其中可以不按所示出或讨论的顺
序，包括根据所涉及的功能按基本同时的方式或按相反的顺序，来执行功能，这应被本技术的实施例所属技术领域的技术人员所理解。
[0218]
在流程图中表示或在此以其他方式描述的逻辑和/或步骤，例如，可以被认为是用于实现逻辑功能的可执行指令的定序列表，可以具体实现在任何计算机可读介质中，以供指令执行系统、装置或设备(如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统)使用，或结合这些指令执行系统、装置或设备而使用。就本说明书而言，"计算机可读介质"可以是任何可以包含、存储、通信、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。计算机可读介质的更具体的示例(非穷尽性列表)包括以下：具有一个或多个布线的电连接部(电子装置)，便携式计算机盘盒(磁装置)，随机存取存储器(ram)，只读存储器(rom)，可擦除可编辑只读存储器(eprom或闪速存储器)，光纤装置，以及便携式光盘只读存储器(cdrom)。另外，计算机可读介质甚至可以是可在其上打印程序的纸或其他合适的介质，因为可以例如通过对纸或其他介质进行光学扫描，接着进行编辑、解译或必要时以其他合适方式进行处理来以电子方式获得程序，然后将其存储在计算机存储器中。
[0219]
应当理解，本技术的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中，多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。例如，如果用硬件来实现，和在另一实施方式中一样，可用本领域公知的下列技术中的任一项或他们的组合来实现：具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路，具有合适的组合逻辑门电路的专用集成电路，可编程门阵列(pga)，现场可编程门阵列(fpga)等。
[0220]
本技术领域的普通技术人员可以理解实现上述实施例方法携带的全部或部分步骤，可以通过程序来指令相关的硬件完成，所开发的程序可以存储于一种计算机可读存储介质中，该程序在执行时，包括方法实施例的步骤之一或其组合。
[0221]
此外，在本技术各个实施例中的各功能单元可以集成在一个处理模块中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个模块中。上述集成的模块既可以采用硬件的形式实现，也可以采用软件功能模块的形式实现。集成的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。
[0222]
上述提到的存储介质可以是只读存储器，磁盘或光盘等。尽管上面已经示出和描述了本技术的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本技术的限制，本领域的普通技术人员在本技术的范围内可以对上述实施例进行变化、修改、替换和变型。