一种含风储联合系统的配电网运行风险动态评估方法

摘要

本发明提供一种含风储联合系统的配电网运行风险动态评估方法，实现配电网运行状态快速在线评估与量化。其具体步骤为：首先，采集配电网的网络拓扑及原始数据，设置风速、光照强度、负荷的概率分布参数；其次，运用虚拟节点解耦法实现风储联合系统风电出力与储能设备出力的解耦；再次，建立配电网运行风险关键评估指标，并运用基于半不变量法的含QV节点概率潮流解析算法计算各指标数值；最后，输出配电网运行风险评估结果。本发明能够解决含并网点电压控制的风储联合系统配电网的运行风险动态评估问题，可为配电网的运行调度，风储联合系统并网点电压整定等提供决策指导，有利于管控电网运行风险、提高系统安全运行水平。

权利要求

1.一种含风储联合系统的配电网运行风险动态评估方法，其特征在于包括以下步骤：(1)采集含风储联合系统的配电网的网络拓扑结构，在该拓扑结构中对各节点进行编号并确定节点类型；

(2)输入原始数据；

(3)采集配电网中风机、光伏电源的技术参数；所述的技术参数包括：接于节点iWT的风机的额定容量 切入风速 额定风速 和切出风速 接于节点iPV的光伏电源的额定容量 额定光照强度 最大光照强度 其中，iWT∈ΩWT，iPV∈ΩPV，ΩWT和ΩPV分别为配电网内风电、光伏电源接入节点的集合；

(4)采用Weibull分布表征节点iWT的风电的风速的概率分布，采用Beta分布表征节点iPV的光伏电源的光照强度的概率分布，采用正态分布表征节点iL的负荷的概率分布，并基于历史统计数据设置风速、光照强度以及负荷的概率分布参数；

(5)由风速与风电出力映射关系计算节点iWT的风电出力，以及由光照强度与光伏电源出力映射关系计算节点iPV的光伏电源出力，具体计算表达式为：上式中， 为节点iWT处的风电的输出功率， 为节点iWT处的风速， 为节点iPV处的光伏电源的输出功率， 为节点iPV处的光照强度；

(6)基于历史统计数据，设置配电网内各风电，光伏电源出力之间的相关系数矩阵R；

(7)运用虚拟节点解耦法对风储联合系统的风电出力与储能设备出力进行解耦；

(8)运用基于半不变量法的含QV节点概率潮流解析算法，计算含风储联合系统的配电网概率潮流；

(9)构建含风储联合系统的配电网运行风险关键评估指标；

(10)基于步骤(8)的概率潮流结果，计算步骤(9)所建立的运行风险关键评估指标值；

(11)输出含风储联合系统的配电网运行风险评估结果；

其中，步骤(7)所述的虚拟节点解耦法是一种通过引入虚拟节点，实现所述风储联合系统的风电出力与储能设备出力解耦的方法，具体步骤如下：(1)引入虚拟PQ节点用以表示所述风储联合系统中风电的并网点；

(2)引入虚拟QV节点用以表示所述风储联合系统中储能设备的并网点。

2.根据权利要求1所述的一种含风储联合系统的配电网运行风险动态评估方法，其特征在于：步骤(1)所述的确定节点类型是根据配电网中各节点的运行特性，对各节点在潮流计算中的节点类型进行设定，包括：在含风储联合系统的配电网的网络拓扑结构中，把并网点电压控制的风储联合系统所在的节点设定为QV节点，把负荷、以及未配有储能设备的光伏电源、未配有储能设备的风电接入的节点设定为PQ节点，把与上层电网连接的节点设定为平衡节点。

3.根据权利要求1所述的一种含风储联合系统的配电网运行风险动态评估方法，其特征在于：步骤(2)所述的原始数据包括：在含风储联合系统的配电网的网络拓扑结构中，相邻节点i、节点j之间线路电阻rij、线路电抗xij、线路接地电纳bi0和bj0、系统平衡节点电压幅值Ubalance、风储联合系统所在节点e的并网点控制电压Upcc,e。

4.根据权利要求1所述的一种含风储联合系统的配电网运行风险动态评估方法，其特征在于：步骤(4)所述的基于历史统计数据设置风速、光照强度以及负荷的概率分布参数，包括：通过各节点采集的风速、光照强度、负荷的历史统计信息，设置风速、光照强度以及负荷的概率分布参数。

5.根据权利要求4所述的一种含风储联合系统的配电网运行风险动态评估方法，其特征在于：所述概率分布参数包括：位于节点iWT的风速在Weibull分布函数中的形状参数 位于节点iWT的风速在Weibull分布函数中的尺度参数 位于节点iPV的光照强度在Beta分布函数中的形状参数及 节点iL的有功负荷的期望值 节点iL的有功负荷标准差 节点iL的无功负荷的期望值 节点iL的无功负荷的标准差 其中，iL∈ΩL，ΩL为配电网内负荷接入节点的集合。

6.根据权利要求1所述的一种含风储联合系统的配电网运行风险动态评估方法，其特征在于：步骤(9)所述的含风储联合系统的配电网运行风险关键评估指标是用以反映和量化含风储联合系统的配电网运行风险的关键指标集合，包括：系统节点电压越限概率，置信度电压波动范围，线路传输功率越限概率。

7.根据权利要求6所述的一种含风储联合系统的配电网运行风险动态评估方法，其特征在于：所述的系统节点电压越限概率采用配电网内所有节点电压越限概率的最大值进行表示，具体表达式为：上式中：λvoltage为系统节点电压越限概率，ΩB为配电网内全部节点集合，Ui为节点i的电压幅值，ProbUi为节点i电压越限概率，Uimax和Uimin分别为节点i电压安全运行的上下限，Pr{·}定义为事件{·}发生的概率。

8.根据权利要求6所述的一种含风储联合系统的配电网运行风险动态评估方法，其特征在于：所述的置信度电压波动范围分别采用在一给定的置信度下系统所有节点电压波动范围的平均值和最大值进行表示，具体表达式为：上式中： 和 分别为在置信度α下配电网内所有节点电压波动范围的平均值和最大值， 和 分别为在置信度α下节点i的电压的上边界和下边界，n为配电网内节点总数，ΩB为配电网内全部节点集合。

9.根据权利要求6所述的一种含风储联合系统的配电网运行风险动态评估方法，其特征在于：所述的线路传输功率越限概率采用配电网内所有支路传输功率越限概率的最大值进行表示，具体表达式为：上式中：λpower为线路传输功率越限概率，l为配电网内支路编号，ΩS为配电网内全部支路的集合，Pl为支路l的传输功率，ProbPl为支路l传输功率越限概率，Plmax和Plmin分别为支路l允许传输功率的上下限，Pr{·}定义为事件{·}发生的概率。

说明书

一种含风储联合系统的配电网运行风险动态评估方法

技术领域

[0001] 本发明涉及含风储联合系统的配电网运行风险动态评估方法，特别涉及含并网点电压控制的风储联合系统的配电网运行风险动态评估方法。

背景技术

[0002] 我国作为世界上能源消费大国，经济发展面临资源和环境等多方面挑战。随着“清洁低碳、安全高效”能源体系的提出与建立，以风电、光伏发电为代表的可再生能源将成为我国未来能源转型与发展的重点方向。电力系统作为能源生产和消费的重要载体，势必需要对可再生能源的兼容包并。

[0003] 近年来，分布式风电在配电网的渗透率快速提升。然而，风电出力固有的波动性与间歇性，大规模接入配电系势必给电网的调度运行带来巨大挑战。具体表现为：1)分布式风电的接入，使配电网潮流形态由传统单向流动向双向流动的特点，由此引起局部馈线传输功率阻塞问题不容忽视；2)风电出力的不确定性，会对配电网节点电压造成扰动，进而引起电压越限和频繁波动等问题，直接影响配电系统的电压质量。

[0004] 为解决上述问题，在风电并网点处配置储能设备构建风储联合系统成为有效解决方案。通过对风储联合系统的并网点运行电压进行合理控制已被证明在含规模化分布式可再生能源的配电网具有良好的电压质量控制效果。

[0005] 对于含高比例风电的配电网，准确评估配电网运行风险，能够为系统调度运行决策等提供辅助信息，有助于量化与管控电网运行风险。此外，运用风险指标识别电网中的高运行风险区域与薄弱环节，有助于提出针对性的风险防范措施，提高配电网运行的安全性与可靠性水平。

[0006] 现有技术中，广泛采用概率潮流对含可再生能源(风电、光伏)配电网的运行风险进行评估。然而，对于含风储联合系统的配电网运行风险的评估，现有技术与方法存在以下几个方面的不足。1)在评估对象方面，现有针对含可再生能源配电网的运行风险评估通常仅研究可再生能源出力不确定性对系统运行风险的影响，缺乏考虑储能设备的引入对可再生能源出力行为的改变以及对评估结果的影响，现有技术尚未有针对并网点电压控制的风储联合系统的配电网的运行风险评估方法。2)在评估指标方面，现有的配电网运行风险评估指标主要关注节点电压和支路功率越限的问题，而对于衡量电压质量的另一个重要指标电压波动的刻画不足。3)在评估方法方面，现有的运行风险评估所涉及的概率潮流计算方法多基于蒙特卡洛模拟法，存在计算量巨大和计算机时过长等问题，不利于在线动态风险评估。

[0007] 针对上述问题，本发明提出了一种含风储联合系统的配电网运行风险评估方法。该方法通过引入虚拟节点实现风电出力和储能设备运行功率的解耦，并建立了反映配电网运行风险的关键指标集合，最后，运用基于半不变量法的含QV节点概率潮流解析算法求解含风储联合系统的配电网概率潮流。该方法能够实现含并网点电压控制的风储联合系统配电网的运行风险动态评估，可为配电网的运行调度，风储联合系统并网点电压整定等提供决策指导，有利于识别电网薄弱环节，管控电网运行风险，提高系统运行安全性，促进配电网对高比例可再生能源的综合消纳。

发明内容

[0008] 本发明的目的在于提供一种含风储联合系统的配电网运行风险动态评估方法，实现含并网点电压控制的风储联合系统的配电网运行风险动态评估，有利于识别配电网中的高运行风险区域，提高含可再生能源配电网的运行风险的管控能力。本发明既可为含风储联合系统的配电网的风险调度提供决策依据，也能为含风储联合系统的配电网薄弱环节的技术改造提供参考。

[0009] 本发明提出一种含风储联合系统的配电网运行风险动态评估方法，包括以下步骤：

[0010] (1)采集含风储联合系统的配电网的网络拓扑结构，在该拓扑结构中对各节点进行编号并确定节点类型；

[0011] (2)输入原始数据；

[0012] (3)采集配电网中风机、光伏电源的技术参数；所述的技术参数包括：接于节点iWT的风机的额定容量 切入风速 额定风速 和切出风速 接于节点iPV的光伏电源的额定容量 额定光照强度 最大光照强度 其中，iWT∈ΩWT，iPV∈ΩPV，ΩWT和ΩPV分别为配电网内风电、光伏电源接入节点的集合；

[0013] (4)采用Weibull分布表征节点iWT的风电的风速的概率分布，采用Beta分布表征节点iPV的光伏电源的光照强度的概率分布，采用正态分布表征节点iL的负荷的概率分布，并基于历史统计数据设置风速、光照强度以及负荷的概率分布参数；

[0014] (5)由风速与风电出力映射关系计算节点iWT的风电出力，以及由光照强度与光伏电源出力映射关系计算节点iPV的光伏电源出力，具体计算表达式为：

[0015]

[0016]

[0017] 式中， 为节点iWT处的风电的输出功率，， 为节点iWT处的风速， 为节点iPV处的光伏电源的输出功率， 为节点iPV处的光照强度；

[0018] (6)基于历史统计数据，设置配电网内各风电，光伏电源出力之间的相关系数矩阵R；

[0019] (7)运用虚拟节点解耦法对风储联合系统的风电出力与储能设备出力进行解耦；

[0020] (8)运用基于半不变量法的含QV节点概率潮流解析算法，计算含风储联合系统的配电网概率潮流；

[0021] (9)构建含风储联合系统的配电网运行风险关键评估指标；

[0022] (10)基于步骤(8)的概率潮流结果，计算步骤(9)所建立的运行风险关键评估指标值；

[0023] (11)输出含风储联合系统的配电网运行风险评估结果。

[0024] 进一步的，步骤(1)所述的确定节点类型是根据配电网中各节点的运行特性，对各节点在潮流计算中的节点类型进行设定，包括：在含风储联合系统的配电网的网络拓扑结构中，把并网点电压控制的风储联合系统所在的节点设定为QV节点，把负荷、以及未配有储能设备的光伏电源、未配有储能设备的风电接入的节点设定为PQ节点，把与上层电网连接的节点设定为平衡节点。

[0025] 进一步的，步骤(2)所述的原始数据包括：在含风储联合系统的配电网的网络拓扑结构中，相邻节点i、节点j之间线路电阻rij、线路电抗xij、线路接地电纳bi0和bj0、系统平衡节点电压幅值Ubalance、风储联合系统所在节点e的并网点控制电压Upcc,e。

[0026] 进一步的，步骤(4)所述的基于历史统计数据设置风速、光照强度以及负荷的概率分布参数，包括：通过各节点采集的风速、光照强度、负荷的历史统计信息，设置风速、光照强度以及负荷的概率分布参数。

[0027] 进一步的，所述概率分布参数包括：位于节点iWT的风速在Weibull分布函数中的形状参数 位于节点iWT的风速在Weibull分布函数中的尺度参数 位于节点iPV的光照强度在Beta分布函数中的形状参数 及 节点iL的有功负荷的期望值 节点iL的有功负荷标准差 节点iL的无功负荷的期望值 节点iL的无功负荷的标准差 其中，iL∈ΩL，ΩL为配电网内负荷接入节点的集合。

[0028] 进一步的，步骤(7)所述的虚拟节点解耦法是一种通过引入虚拟节点，实现风储联合系统的风电出力与储能设备出力解耦的方法，具体步骤如下：

[0029] (1)引入虚拟PQ节点用以表示风储联合系统中风电的并网点；

[0030] (2)引入虚拟QV节点用以表示风储联合系统中储能设备的并网点。

[0031] 进一步的，步骤(9)所述的含风储联合系统的配电网运行风险关键评估指标是用以反映和量化含风储联合系统的配电网运行风险的关键指标集合，包括：系统节点电压越限概率，置信度电压波动范围，线路传输功率越限概率。

[0032] 进一步的，所述的系统节点电压越限概率采用配电网内所有节点电压越限概率的最大值进行表示，具体表达式为：

[0033]

[0034] 式中：λvoltage为系统节点电压越限概率，ΩB为配电网内全部节点集合，Ui为节点i的电压幅值，ProbUi为节点i电压越限概率，Uimax和Uimin分别为节点i电压安全运行的上下限，Pr{·}定义为事件{·}发生的概率。

[0035] 进一步的，所述的置信度电压波动范围分别采用在一给定的置信度下系统所有节点电压波动范围的平均值和最大值进行表示，具体表达式为：

[0036]

[0037] 式中： 和 分别为在置信度α下配电网内所有节点电压波动范围的平均值和最大值， 和 分别为在置信度α下节点i的电压的上边界和下边界，n为配电网内节点总数，ΩB为配电网内全部节点集合。

[0038] 进一步的，所述的线路传输功率越限概率采用配电网内所有支路传输功率越限概率的最大值进行表示，具体表达式为：

[0039]

[0040] 式中：λpower为线路传输功率越限概率，l为配电网内支路编号，ΩS为配电网内全部支路的集合，Pl为支路l的传输功率，ProbPl为支路l传输功率越限概率，Plmax和Plmin分别为支路l允许传输功率的上下限，Pr{·}定义为事件{·}发生的概率。

[0041] 与现有技术相比，本发明的有益效果在于：

[0042] (1)提出针对含并网点电压控制的风储联合系统的配电网运行风险动态评估方法，引入虚拟节点实现风储联合系统的风电出力与储能设备出力的解耦，考虑了储能设备的引入对可再生能源出力行为改变的影响，解决了现有技术中尚未有针对并网点电压控制的风储联合系统的配电网的运行风险评估的难题。

[0043] (2)所提出的配电网运行风险评估方法通过建立反映配电网运行风险的关键指标集，并基于半不变量法的含QV节点概率潮流解析算法求解计算各指标值，实现了含风储联合系统的配电网运行风险的动态评估与量化，为配电网的运行调度，风储联合系统并网点电压整定等提供决策指导，有利于管控电网运行风险、提高系统安全运行水平。

附图说明

[0044] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。

[0045] 图1是本发明提供的一种含风储联合系统的配电网运行风险动态评估方法的流程示意图。

[0046] 图2是虚拟节点解耦示意图。

[0047] 图3是本发明实施例提供的含风储联合系统的配电网拓扑结构示意图。

[0048] 图4是本发明实施例提供的引入虚拟节点解耦示意图。

[0049] 图5是本发明实施例提供的场景1节点电压期望值和90％置信区间。

[0050] 图6是本发明实施例提供的场景2节点电压期望值和90％置信区间。

[0051] 图7是本发明实施例提供的两种方法下节点15的电压PDF和CDF对比。

具体实施方式

[0052] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

[0053] 本发明实施例提供一种含风储联合系统的配电网运行风险动态评估方法，图1为本发明实施例提供的一种含风储联合系统的配电网运行风险动态评估方法的流程示意图，参照图1，该方法包括：

[0054] (1)采集含风储联合系统的配电网的网络拓扑结构，在该拓扑结构中对各节点进行编号并确定节点类型；

[0055] 其中，所述的确定节点类型是根据配电网中各节点的运行特性，对各节点在潮流计算中的节点类型进行设定，包括：在含风储联合系统的配电网的网络拓扑结构中，把并网点电压控制的风储联合系统所在的节点设定为QV节点，把负荷、以及未配有储能设备的光伏电源、未配有储能设备的风电接入的节点设定为PQ节点，把与上层电网连接的节点设定为平衡节点。

[0056] (2)输入原始数据；

[0057] 其中，所述的原始数据包括：在含风储联合系统的配电网的网络拓扑结构中，相邻节点i、节点j之间线路电阻rij、线路电抗xij、线路接地电纳bi0和bj0、系统平衡节点电压幅值Ubalance、风储联合系统所在节点e的并网点控制电压Upcc,e。

[0058] (3)采集配电网中风机、光伏电源的技术参数；所述的技术参数包括：接于节点iWT的风机的额定容量 切入风速 额定风速 和切出风速 接于节点iPV的光伏电源的额定容量 额定光照强度 最大光照强度 其中，iWT∈ΩWT，iPV∈ΩPV，ΩWT和ΩPV分别为配电网内风电、光伏电源接入节点的集合。

[0059] (4)可选的，采用Weibull分布表征节点iWT的风电的风速的概率分布，采用Beta分布表征节点iPV的光伏电源的光照强度的概率分布，采用正态分布表征节点iL的负荷的概率分布，并基于历史统计数据设置风速、光照强度以及负荷的概率分布参数。

[0060] 其中，所述的采用Weibull分布表征节点iWT的风电的风速的概率分布，具体函数表达式为：

[0061]

[0062] 式中， 为位于节点iWT的风速， 分别为位于节点iWT的风速在Weibull分布函数中的形状参数和尺度参数；

[0063] 其中，所述的采用Beta分布表征节点iPV的光伏电源的光照强度的概率分布，具体函数表达式为：

[0064]

[0065] 式中， 为位于节点iPV的光照强度， 分别为位于节点iPV的光照强度在Beta分布函数中的形状参数；

[0066] 其中，所述的采用正态分布表征节点iL的负荷的概率分布，具体函数表达式为：

[0067]

[0068] 式中， 和 分别为节点iL的有功负荷和无功负荷的随机量， 和 分别为节点iL的有功负荷和无功负荷的期望值， 和 分别为节点iL的有功负荷和无功负荷的标准差，iL∈ΩL，ΩL为配电网内负荷接入节点的集合；

[0069] 其中，所述的基于历史统计数据设置风速、光照强度以及负荷的概率分布参数，包括：通过各节点采集的风速、光照强度、负荷的历史统计信息，设置风速、光照强度以及负荷的概率分布参数；

[0070] 其中，所述的概率分布参数，包括：位于节点iWT的风速在Weibull分布函数中的形状参数 位于节点iWT的风速在Weibull分布函数中的尺度参数 位于节点iPV的光照强度在Beta分布函数中的形状参数 及 节点iL的有功负荷的期望值 节点iL的有功负荷标准差 节点iL的无功负荷的期望值 节点iL的无功负荷的标准差 其中，iL∈ΩL，ΩL为配电网内负荷接入节点的集合。

[0071] (5)由风速与风电出力映射关系计算节点iWT的风电出力，以及由光照强度与光伏电源出力映射关系计算节点iPV的光伏电源出力，具体计算表达式为：

[0072]

[0073]

[0074] 式中， 为节点iWT处的风电的输出功率， 为节点iWT处的风速， 为节点iPV处的光伏电源的输出功率， 为节点iPV处的光照强度。

[0075] (6)基于历史统计数据，设置配电网内各风电，光伏电源出力之间的相关系数矩阵R。

[0076] 其中，所述的相关系数矩阵R可基于步骤(4)所获取的关于各风电、光伏电源接入节点的风速、光照强度的历史统计信息进行设置。

[0077] (7)运用虚拟节点解耦法对风储联合系统的风电出力与储能设备出力进行解耦。

[0078] 参照图2，所述的虚拟节点解耦法是一种通过引入虚拟节点，实现风储联合系统的风电出力与储能设备出力解耦的方法，具体步骤如下：

[0079] 1)引入虚拟PQ节点用以表示风储联合系统中风电的并网点；

[0080] 2)引入虚拟QV节点用以表示风储联合系统中储能设备的并网点。

[0081] (8)运用基于半不变量法的含QV节点概率潮流解析算法，计算含风储联合系统的配电网概率潮流。

[0082] 例如将相关参数代入评估方法中，在MATLAB平台上编写基于半不变量法的含QV节点概率潮流解析算法，进而计算含风储联合系统的配电网概率潮流。

[0083] (9)构建含风储联合系统的配电网运行风险关键评估指标。

[0084] 其中，所述的含风储联合系统的配电网运行风险关键评估指标是用以反映和量化含风储联合系统的配电网运行风险的关键指标集合，包括系统节点电压越限概率、置信度电压波动范围，线路传输功率越限概率；

[0085] 其中，所述的系统节点电压越限概率采用配电网内所有节点电压越限概率的最大值进行表示，具体表达式为：

[0086]

[0087] 式中：λvoltage为系统节点电压越限概率，ΩB为配电网内全部节点集合，Ui为节点i的电压幅值，ProbUi为节点i电压越限概率，Uimax和Uimin分别为节点i电压安全运行的上下限，Pr{·}定义为事件{·}发生的概率；

[0088] 其中，所述的置信度电压波动范围分别采用在一给定的置信度下系统所有节点电压波动范围的平均值和最大值进行表示，具体表达式为：

[0089]

[0090] 式中： 和 分别为在置信度α下配电网内所有节点电压波动范围的平均值和最大值， 和 分别为在置信度α下节点i的电压的上边界和下边界，n为配电网内节点总数，ΩB为配电网内全部节点集合，；

[0091] 其中，所述的线路传输功率越限概率采用配电网内所有支路传输功率越限概率的最大值进行表示，具体表达式为：

[0092]

[0093] 式中：λpower为线路传输功率越限概率，l为配电网内支路编号，ΩS为配电网内全部支路的集合，Pl为支路l的传输功率，ProbPl为支路l传输功率越限概率，Plmax和Plmin分别为支路l允许传输功率的上下限，Pr{·}定义为事件{·}发生的概率。

[0094] (10)基于步骤(8)的概率潮流结果，计算步骤(9)所建立的运行风险关键评估指标值。

[0095] (11)输出含风储联合系统的配电网运行风险评估结果。

[0096] 本发明实施例提供的技术方案，通过采集含风储联合系统的配电网的网络拓扑结构，在该拓扑结构中对各节点进行编号并确定节点类型；输入原始数据，采集配电网中风机、光伏电源的技术参数；基于历史统计数据设置风速、光照强度以及负荷的概率分布参数；由风速与风电出力映射关系以及光照强度与光伏电源出力映射关系分别计算风电出力和光伏电源出力；设置配电网内各风电，光伏电源出力之间的相关系数矩阵；运用虚拟节点解耦法对风储联合系统的风电出力与储能设备出力进行解耦；运用基于半不变量法的含QV节点概率潮流解析算法，计算含风储联合系统的配电网概率潮流；构建含风储联合系统的配电网运行风险关键评估指标；基于概率潮流结果，计算所建立的运行风险关键评估指标值；输出含风储联合系统的配电网运行风险评估结果。总体思路为利用数据采集以及参数设置后进行含QV节点的概率潮流求解，并计算相应的评估指标数值，实现含风储联合系统的配电网运行风险动态评估。本发明能够解决含并网点电压控制的风储联合系统配电网的运行风险动态评估问题，有利于管控电网运行风险、提高系统安全运行水平。

[0097] 本发明实施例提供一种含风储联合系统的配电网运行风险动态评估方法，以IEEE33节点配电系统作为本发明的实施例进行仿真计算。

[0098] 示例性的，图3为本发明实施例提供的含风储联合系统的配电网拓扑结构示意图，其中附图中的数字为对应的节点编号。

[0099] (1)采集含风储联合系统的配电网的网络拓扑结构，在该拓扑结构中对各节点进行编号并确定节点类型。

[0100] 所述的确定节点类型是根据配电网中各节点的运行特性，对各节点在潮流计算中的节点类型进行设定，包括：在含风储联合系统的配电网的网络拓扑结构中，把并网点电压控制的风储联合系统所在的节点设定为QV节点，把负荷、以及未配有储能设备的光伏电源、未配有储能设备的风电接入的节点设定为PQ节点，把与上层电网连接的节点设定为平衡节点。

[0101] 网络拓扑结构及节点编号参见图3。在图3中，节点1为首端节点，与外部电网相连，设为平衡节点。未配有储能设备的光伏电源接入节点为节点11和节点30；未配有储能设备的风电接入节点为节点5、节点22；节点18为风储联合系统，并采用并网点电压控制运行。因此，节点1为平衡节点，节点18为QV节点，其余为PQ节点。

[0102] (2)输入原始数据。

[0103] 所述的原始数据包括：在含风储联合系统的配电网的网络拓扑结构中，相邻节点i、节点j之间线路电阻rij、线路电抗xij、线路接地电纳bi0和bj0、系统平衡节点电压幅值Ubalance、风储联合系统所在节点e的并网点控制电压Upcc,e。

[0104] 本实施例网络拓扑结构中(参见图3)，网架线路参数见表1，系统平衡节点电压幅值Ubalance设为1.05(标幺值，以下同)，节点18风储联合系统并网点控制电压Upcc,e设为1.0。

[0105] 表1实施例网架线路参数

[0106]

[0107]

[0108]

[0109] (3)采集配电网中风机、光伏电源的技术参数；所述的技术参数包括：接于节点iWT的风机的额定容量 切入风速 额定风速 和切出风速 接于节点iPV的光伏电源的额定容量 额定光照强度 最大光照强度 其中，iWT∈ΩWT，iPV∈ΩPV，ΩWT和ΩPV分别为配电网内风电、光伏电源接入节点的集合。

[0110] 本实施例中，节点5和节点22的风机的额定容量均为300kW，即节点18风机额定容量为500kW， 切入风速
额定风速 和切出风速 近似的，
ΩWT＝[5,18,22]，ΩPV＝[11,
30]。

[0111] (4)可选的，采用Weibull分布表征节点iWT的风电的风速的概率分布，采用Beta分布表征节点iPV的光伏电源的光照强度的概率分布，采用正态分布表征节点iL的负荷的概率分布，并基于历史统计数据设置风速、光照强度以及负荷的概率分布参数。

[0112] 本实施例中， 各节点 和 的取值如表2所示。

[0113] 表2各节点有功负荷和无功负荷期望值

[0114]节点名 有功负荷/kW 无功负荷/kW 节点名 有功负荷/kW 无功负荷/kW
2 100 60 18 90 40
3 90 40 19 90 40
4 120 80 20 90 40
5 60 30 21 90 40
6 60 20 22 90 40
7 200 100 23 90 50
8 200 100 24 420 200
9 60 20 25 420 200
10 60 20 26 60 25
11 45 30 27 60 25
12 60 35 28 60 20
13 60 35 29 120 70
14 120 80 30 200 600
15 60 10 31 150 70
16 60 20 32 210 100
17 60 20 33 60 40

[0115] 各节点 和 均设为期望值的10％，即

[0116] (5)由风速与风电出力映射关系计算节点iWT的风电出力，以及由光照强度与光伏电源出力映射关系计算节点iPV的光伏电源出力，具体计算表达式为：

[0117]

[0118]

[0119] 式中， 为节点iWT处的风电的输出功率， 为节点iWT处的风速， 为节点iPV处的光伏电源的输出功率， 为节点iPV处的光照强度。

[0120] (6)基于历史统计数据，设置配电网内各风电，光伏电源出力之间的相关系数矩阵R。

[0121] 所述的相关系数矩阵R可基于步骤(4)所获取的关于各风电、光伏电源接入节点的风速、光照强度的历史统计信息进行设置。

[0122] 实施例中，光伏电源和风电相关性系数为‑0.5；各节点光伏电源之间相关性系数为0.7；各节点风电之间相关性系数为0.4，由此R可设置为：

[0123] (7)运用虚拟节点解耦法对风储联合系统的风电出力与储能设备出力进行解耦。

[0124] 图4为本实施例提供的针对节点18风储联合系统所引入虚拟节点解耦示意图，参照图4，引入虚拟节点对节点18进行解耦包括：

[0125] 1)引入虚拟PQ节点(节点34)用以表示风储联合系统中风电的并网点；

[0126] 2)引入虚拟QV节点(节点35)用以表示风储联合系统中储能设备的并网点。

[0127] (8)运用基于半不变量法的含QV节点概率潮流解析算法，计算含风储联合系统的配电网概率潮流。

[0128] (9)构建含风储联合系统的配电网运行风险关键评估指标。

[0129] 其中，所述的含风储联合系统的配电网运行风险关键评估指标是用以反映和量化含风储联合系统的配电网运行风险的关键指标集合，包括系统节点电压越限概率、置信度电压波动范围，线路传输功率越限概率，具体计算方法如下：

[0130] 其中，所述的系统节点电压越限概率采用配电网内所有节点电压越限概率的最大值表示，具体表达式为：

[0131]

[0132] 式中：λvoltage为系统节点电压越限概率，ΩB为配电网内全部节点集合，Ui为节点i的电压幅值，ProbUi为节点i电压越限概率，Uimax和Uimin分别为节点i电压安全运行的上下限，Pr{·}定义为事件{·}发生的概率。

[0133] 本实施例中，Uimax和Uimin分别设置为1.07和0.93，Pr{Uimin≤Ui≤Uimax}由步骤(8)概率潮流计算结果求得。

[0134] 其中，所述的置信度电压波动范围分别采用在一给定的置信度下系统所有节点电压波动范围的平均值和最大值进行表示，具体表达式为：

[0135]

[0136] 式中： 和 分别为在置信度α下配电网内所有节点电压波动范围的平均值和最大值， 和 分别为在置信度α下节点i的电压的上边界和下边界，n为配电网内节点总数，ΩB为配电网内全部节点集合。

[0137] 本实施例中，置信度α设置为0.9，n＝33， 和 由步骤(8)概率潮流计算结果求得。

[0138] 其中，所述的线路传输功率越限概率采用配电网内所有支路传输功率越限概率的最大值进行表示，具体表达式为：

[0139]

[0140] 式中：λpower为线路传输功率越限概率，l为配电网内支路编号，ΩS为配电网内全部支路的集合，Pl为支路l的传输功率，ProbPl为支路l传输功率越限概率，Plmax和Plmin分别为支路l允许传输功率的上下限，Pr{·}定义为事件{·}发生的概率。

[0141] 本实施例中，Plmax和Plmin分别设置为3000kW和‑3000kW，Pr{Plmin≤Pl≤Plmax}由步骤(8)概率潮流计算结果求得。

[0142] (10)基于步骤(8)的概率潮流结果，计算步骤(9)所建立的运行风险关键评估指标值。

[0143] (11)输出含风储联合系统的配电网运行风险评估结果。

[0144] 将相关参数代入评估方法中，在MATLAB平台上编写基于半不变量法的含QV节点概率潮流解析算法，进而计算含风储联合系统的配电网概率潮流，基于概率潮流结果计算运行风险关键评估指标，各评估指标的计算结果如表3所示。

[0145] 表3各评估指标的计算结果

[0146]

[0147] 为了验证本发明所提的评估方法的合理性与准确性，设计了如下对比场景。

[0148] 场景1：节点18不接入储能设备。

[0149] 场景2：采用本实施例的运行方案，在节点18接入储能设备构成风储联合系统。

[0150] 两种场景下评估指标对比见表4，图5和图6分别为发明实施例提供的场景1和场景2节点电压的期望值和90％置信区间。

[0151] 表4两种场景下的评估指标对比

[0152]

[0153] 通过实施例分析表4、图5和图6，可得出如下结论：

[0154] 1)场景1没有储能设备提供功率的灵活调节，本发明所提的3个运行风险关键评估指标：系统节点电压越限概率、置信度电压波动范围、线路传输功率越限概率均明显高于场景2。其中，系统节点电压越限概率和线路传输功率越限概率分别高达19.73％和34.67％，表明此时系统存在较大的电压越限和传输功率越限的风险。图5的90％概率边界能够直观显示在可再生能源出力和负荷波动的环境下系统节点电压可能的运行范围，表明：光伏电源、风电所接入的节点会引起电压波动的增大，同时，越靠近线路末端，节点电压的波动范围也越大。节点18因为处于主干线最末端，且有风电接入，属于系统的薄弱节点，其波动范围高达0.0822p.u.，并伴有一定概率的电压越限风险。

[0155] 2)场景2在节点18处接入储能设备构建并网点电压控制的风储联合系统，评估指标中：电压越限概率、置信度电压波动范围、线路传输功率越限概率均得到大幅降低(电压越限概率仅为0.033％，线路传输功率越限概率仅为2.4％)，充分体现在薄弱节点处接入储能设备并采用并网点电压控制能够有效提升全系统电压质量水平，改善线路潮流分布。由图6可知，由于节点18为QV节点，附近的节点也得益于其并网点电压控制，相比图5电压质量得到明显提升。

[0156] 3)表4的两个场景下的评估结果对比表明，本发明所提的运行风险评估指标能够直接反映系统电压越限和波动，线路传输功率越限等关键表征运行风险的指标，评估结果的对比也表明所提指标的合理性与有效性。此外，本发明所提的评估方法所需输入的电气量为目前电力系统常用的基本电气量，无需额外增加传感器等其他量测装置，具备指标易获取的优点。

[0157] 为进一步验证本发明所提的含QV节点概率潮流解析算法的准确性和优越性，将本发明的方法与传统基于蒙特卡洛概率潮流算法进行对比。采用采样规模为10000次的蒙特卡洛法求取的概率潮流结果与本发明所提的含QV节点概率潮流解析算法的结果进行对比，以评估本发明所提评估方法的准确性。表5列出本发明所提的概率潮流解析算法与蒙特卡洛法的相对误差情况。图7为节点15的电压的概率密度函数(probability density function，PDF)和累计分布函数(cumulative distribution function，CDF)的对比。

[0158] 表5误差情况

[0159]

[0160] 分析由表5和图7可知，本发明所提的含QV节点概率潮流解析算法与传统基于蒙特卡洛法的概率潮流计算算法结果非常接近，期望值的相对误差基本在0.5％以内，标准差的相对误差基本在5％以内，满足工程应用需求。图7的PDF和CDF曲线对比可以看出，两条曲线高度重叠，表明本发明所提的含QV节点概率潮流解析算法能够获得高精度的PDF和CDF曲线。

[0161] 在计算耗时对比方面，本发明仅需1.62秒即可完成一次含风储联合系统配电网的运行风险评估，而基于蒙特卡洛法的运行风险评估则总耗时高达628.3秒，由此表明本发明提供的含风储联合系统的配电网运行风险评估方法在计算机时和效率方面具有明显优势，具备应用于在线计算与运行风险动态评估的可行性。

[0162] 通过实施例分析表明，本发明所提的一种含风储联合系统的配电网运行风险动态评估方法具有以下几个方面的优势：

[0163] 1)本发明能够解决含并网点电压控制的风储联合系统配电网运行风险评估问题，考虑了储能设备对风电出力行为的影响，丰富了配电网运行风险评估的对象；

[0164] 2)所构建的运行风险关键评估指标能够充分合理反映可再生能源与负荷不确定性对系统运行的影响，实现运行风险的量化，有利于运行人员管控电网运行风险；

[0165] 3)本发明的运行风险评估方法在计算机时与效率方面具有明显优势，具备应用于在线计算与运行风险动态评估的可行性，满足未来智能电网运行风险实时在线评估的需求。

[0166] 最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。