CN102790389B - 一种暂态稳定裕度概率分布的计算方法 - Google Patents
一种暂态稳定裕度概率分布的计算方法 Download PDFInfo
- Publication number
- CN102790389B CN102790389B CN201210263026.6A CN201210263026A CN102790389B CN 102790389 B CN102790389 B CN 102790389B CN 201210263026 A CN201210263026 A CN 201210263026A CN 102790389 B CN102790389 B CN 102790389B
- Authority
- CN
- China
- Prior art keywords
- power
- stability margin
- group
- driving power
- failure model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000001052 transient effect Effects 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 title claims abstract description 18
- 230000035945 sensitivity Effects 0.000 claims abstract description 45
- 238000004088 simulation Methods 0.000 claims description 26
- 238000004458 analytical method Methods 0.000 claims description 25
- 238000000205 computational method Methods 0.000 claims description 23
- 238000004364 calculation method Methods 0.000 claims description 11
- 230000008859 change Effects 0.000 claims description 11
- 230000001186 cumulative effect Effects 0.000 claims description 10
- 238000005315 distribution function Methods 0.000 claims description 10
- 238000006467 substitution reaction Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 230000003068 static effect Effects 0.000 claims description 2
- 230000009897 systematic effect Effects 0.000 claims description 2
- 238000012821 model calculation Methods 0.000 claims 2
- 238000013461 design Methods 0.000 abstract description 3
- 238000004445 quantitative analysis Methods 0.000 abstract description 3
- 239000000243 solution Substances 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 208000024780 Urticaria Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004064 recycling Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000000528 statistical test Methods 0.000 description 1
- 238000002948 stochastic simulation Methods 0.000 description 1
- 238000004804 winding Methods 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种暂态稳定裕度概率分布的计算方法,将分布在电力系统中的各个节点功率的概率数字特征,变换成失稳模式对应的两群驱动功率和负荷功率的概率数字特征;根据扩展等面积准则EEAC,将系统稳定裕度按Taylor级数表示成以两群驱动功率和负荷功率为变量、以稳定裕度灵敏度为系数的线性化表达式,以求取稳定裕度的概率数字特征;再利用Gram-Charlier级数计算出系统暂态稳定裕度的概率分布。可直接用于指导电力系统规划、设计、运行和控制等各个环节中安全稳定量化分析与优化决策的具体实践。
Description
技术领域
本发明属于电力系统及其自动化技术领域,具体涉及一种应用在节点功率注入不确定条件下的暂态稳定裕度概率分布的计算方法。
背景技术
电力工业是国民经济发展中最重要的基础能源产业,电力系统的安全稳定高效运行具有重大意义。大规模新能源发电的并网运行和需求侧响应机制是现代电网的重要特征,但风能、太阳能发电的随机性、间歇性和波动性特点,以及用电侧灵活多变的需求方式,却给电力系统的安全稳定特性带来了深刻的变化。
稳定分析是系统安全性评估的不可或缺的一部分。传统的稳定分析方法局限于用确定性的观点对待系统工况和扰动场景,凭经验根据时域仿真得到的系统受扰轨迹进行稳定性的定性判断,因而无法量化主要随机因素对稳定性的影响。针对传统稳定性分析准则所面临的挑战,国内外的发展趋势是采用风险的观点审视电力系统的规划、设计、运行和控制等各个环节,而概率暂态稳定分析是实现安全性风险评估技术的基础。
电力系统概率暂态稳定分析的任务是根据影响稳定性的主要随机因素的统计特性来确定系统的暂态稳定性概率指标。已有研究成果中,蒙特卡罗仿真法基于对随机变量的统计试验和随机模拟,稳定域法采用系统工况点处于稳定域外的比例来反映系统失稳概率。这些方法只给出系统稳定或失稳的概率,却无法评估系统稳定或失稳的程度,更无法提供决策所需要的灵敏度信息,而且,由于计算费时或方法复杂等原因,难以处理不确定变量个数较多的情况。
间歇性电源和用电负荷分布在电力系统广域空间中,构成了大量的节点功率注入不确定性因素。鉴于在针对确定性系统工况和扰动场景的安全稳定评估技术领域,扩展等面积准则(EEAC)是惟一得到严格证明并投入工程应用的稳定性量化分析工具,可给出稳定裕度、稳定模式、控制措施灵敏度等信息,因此,有必要将EEAC向处理不确定因素的方向扩展,提出新的节点功率注入不确定条件下暂态稳定裕度概率分布的计算方法。
发明内容
针对现有技术的不足,本发明提供在一种节点功率注入不确定条件下的暂态稳定裕度概率分布的计算方法,为含有分布在广域空间中多个随机变量的电力系统稳定裕度概率分布计算提供了高效率的解决方法。不确定是指节点功率注入量随机,但遵循概率分布。
本发明提供的一种暂态稳定裕度概率分布的计算方法,其改进之处在于,所述方法包括如下步骤:
(1)确定初始状态下受扰系统的失稳模式和稳定裕度初始值η0,并由用户确定采用的计算模式,包括基于纯解析法的暂态稳定裕度概率分布计算和综合时域仿真与解析法的暂态稳定裕度概率分布计算;
所述基于纯解析法的暂态稳定裕度概率分布计算包括如下步骤:
(2)计算所述受扰系统稳定裕度对临界群S驱动功率的灵敏度系数I:令临界群S的驱动功率Pms和余下群A的驱动功率Pma的初始值分别等于各自群内机组驱动功率均值的总和,通过静态EEAC(SEEAC)算法求取稳定裕度对临界群驱动功率的灵敏度系数I;所述灵敏度系数I包括一阶灵敏度系数和二阶灵敏度系数其中临界群和余下群中发电机机组通常有多个,临界群S的驱动功率Pms和余下群A的驱动功率Pma的初始值分别等于各自群内不同机组的驱动功率均值求和。
(3)计算临界群S驱动功率的概率数字特征,即各阶半不变量γυ(Pms):根据临界群内各发电机组驱动功率的概率分布,先求取每台发电机驱动功率Pmk的各阶半不变量γυ(Pmk),再按临界群S中发电机组个数对所述各阶半不变量γυ(Pmk)求和,求取Pms的各阶半不变量γυ(Pms),为
其中,k为机组编号,υ表示半不变量的阶数,Ns为临界群中发电机组个数。
(4)求取受扰系统稳定裕度的各阶中心矩βυ(η):利用稳定裕度η对临界群S驱动功率Pms的灵敏度系数I和临界群S驱动功率的各阶半不变量γυ(Pms),计算稳定裕度的各阶半不变量γυ(η),进而通过半不变量和中心矩的标准数学转换公式,求取稳定裕度的各阶中心矩βυ(η)。
(5)计算稳定裕度的累积分布函数,完成理想两群失稳模式的暂态稳定裕度概率分布计算。将稳定裕度的初始值η0及其各阶中心矩βυ(η)代入Gram-Charlier级数计算公式,获得受扰系统稳定裕度η的累积分布函数,完成了基于纯解析法、针对理想两群失稳模式的暂态稳定裕度概率分布的快速计算。其中,理想两群是指互补群群内各发电机转子角轨迹间的间隙在整个动态过程中保持不变。
所述综合时域仿真与解析法的暂态稳定裕度概率分布计算包括如下步骤:
(6)确定节点功率变化范围内的失稳模式变化情况:根据系统中具有功率不确定性的发电机和负荷功率变化范围的上下限,计算出对应步骤(1)初始分群的临界群驱动功率Pms、余下群驱动功率Pma、临界群负荷功率Pls、余下群负荷功率Pla变化范围的上下限,分别按照增加临界群驱动功率和余下群负荷功率,同时减少余下群驱动功率和临界群负荷功率,以及减少临界群驱动功率和余下群负荷功率,同时增加余下群驱动功率和临界群负荷功率的两个方向,依据人工设定的功率变化步长及其在各机组和负荷之间的分配规则,调整并形成一组系统稳态运行工况,再针对每一个系统稳态运行工况,进行针对设定故障下的详细时域仿真,识别失稳模式变化,并记录每一个新的失稳模式i所对应的稳定裕度η′(i),以及临界群的驱动功率Pms0(i)、临界群的负荷功率Pls0(i)、余下群驱动功率Pma0(i)和余下群的负荷功率Pla0(i)。节点功率变化范围是指所计算系统中节点功率实际变化范围,由各节点功率概率分布的上下限确定,初始分群是指系统中各节点功率等于其概率分布均值时,针对设定故障,利用详细时域仿真和扩展等面积准则(EEAC)得到的系统分群情况。
(7)针对每个失稳模式,计算稳定裕度的灵敏度系数II:对每个失稳模式,采用摄动法,通过时域仿真和EEAC准则求取稳定裕度对所述临界群驱动功率Pms的灵敏度系数临界群负荷功率Pls的灵敏度系数余下群驱动功率Pma的灵敏度系数和余下群负荷功率Pla的灵敏度系数
(8)针对每个失稳模式,求取该失稳模式下稳定裕度线性变化对应的各节点功率等于其概率分布均值时的稳定裕度初始值。若失稳模式的个数为1,则稳定裕度初始值η′0(1)=η0;若失稳模式有变化,则根据稳定裕度的线性化展开式,结合步骤(7)得到的灵敏度系数II和步骤(6)得到的失稳模式变化后稳定裕度η′(i),利用数学反推获得每个失稳模式下的稳定裕度初始值η′0(i)。
(9)针对每个失稳模式,计算临界群驱动功率Pms的各阶半不变量、临界群负荷功率Pls的各阶半不变量、余下群驱动功率Pma的各阶半不变量和余下群负荷功率Pla的各阶半不变量:先根据系统内所有机组驱动功率和负荷功率的概率分布,求取单个节点的机组驱动功率和负荷功率的各阶半不变量γυ,再根据每个失稳模式的分群情况,分别按照临界群及余下群内机组和负荷节点的个数分类求和,获得Pms、Pls、Pma和Pla的各阶半不变量。
(10)针对每个失稳模式,求取稳定裕度η的各阶中心矩:利用步骤(7)得到的稳定裕度对所述临界群驱动功率Pms的灵敏度系数II、临界群负荷功率Pls的灵敏度系数II、余下群驱动功率Pma的灵敏度系数II和余下群负荷功率Pla的灵敏度系数II,以及步骤(9)得到的所述临界群驱动功率Pms的各阶半不变量、临界群负荷功率Pls的各阶半不变量、余下群驱动功率Pma的各阶半不变量和余下群负荷功率Pla的各阶半不变量,计算稳定裕度η的各阶半不变量,进而通过半不变量和中心矩的标准数学转换公式,求取稳定裕度η的各阶中心矩βυ(η)。
(11)求取稳定裕度的累积分布函数:针对每个失稳模式,将稳定裕度的各阶中心矩βυ(η)和稳定裕度初始值η′0(i)代入Gram-Charlier级数计算公式,获得稳定裕度η的累积分布函数。若在节点功率变化范围内失稳模式发生变化,则按失稳模式进行分段,通过分段函数综合每个失稳模式对应的稳定裕度概率分布,完成了综合采用时域仿真与解析方法的暂态稳定裕度概率分布计算。
其中,所述步骤(4)中,在计算稳定裕度η的方差D(η)(即2阶中心距)时,按式(1)计入临界群驱动功率Pms方差E(Δ2Pms)的影响,而在计算稳定裕度的高阶(大于等于3阶)中心矩时则忽略Pms方差的影响;其中
其中,所述步骤(4)中,按式(2)求取稳定裕度η对临界群驱动功率Pms的一阶灵敏度系数按式(3)求取稳定裕度的各阶半不变量γυ(η);
其中,所述步骤(10)中,按式(4)根据时域仿真和EEAC获得的灵敏度系数 和求取稳定裕度η的各阶半不变量γυ(η);
其中,所述步骤(8)中,对每个失稳模式i,根据稳定裕度的线性化展开关系,按式(5)进行数学反推,求取各节点功率等于其概率分布均值时的稳定裕度初始值η′0(i)。
式中 和分别为各节点功率等于其概率分布均值时的临界群驱动功率、临界群负荷功率、余下群驱动功率和余下群负荷功率。
与现有技术比,本发明的有益效果为:
1)本发明将分布在电力系统中的各个节点功率的概率数字特征,变换成失稳模式对应的两群驱动功率和负荷功率的概率数字特征,再与系统稳定裕度进行数学关联,为含有分布在广域空间中多个随机变量的电力系统稳定裕度概率分布计算提供了高效率的解决方法。
2)本发明将基于扩展等面积准则(EEAC)的稳定裕度线性化展开式与Gram-Charlier级数有机结合,先求解系统稳定裕度的概率数字特征,再计算其概率分布,避免了大量的时域仿真计算,提高了概率暂态稳定评估的效率。
3)本发明提出的纯解析方法能够计及系统内发电机机组驱动功率的不确定性,在理想两群失稳模式下计算精度高,具有速度快的特点;综合采用时域仿真与解析方法的混合算法能计及机组驱动功率和负荷功率的不确定性,通过时域仿真求取各项灵敏度系数,适用于理想两群、非理想两群和多个失稳模式的情况,具有精度高的优点。
4)本发明为系统运行工况不确定性条件下的概率稳定裕度计算提供了解决方案,可直接用于指导电力系统规划、设计、运行和控制等各个环节中安全稳定量化分析与优化决策的具体实践。
附图说明
图1为本发明提供的一种暂态稳定裕度概率分布的计算方法流程图。
图2为本发明提供的3机9节点系统接线图。
图3为本发明提供的单一失稳模式下暂态稳定裕度概率分布计算结果示意图。
图4为本发明提供的系统稳定裕度随临界群驱动功率变化的曲线图。
图5为本发明提供的用分段函数描述的系统稳定裕度概率分布图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步的详细说明。
本发明的主要过程是:将分布在电力系统中的各个节点功率的概率数字特征,变换成失稳模式对应的两群驱动功率和负荷功率的概率数字特征;根据扩展等面积准则EEAC,将系统稳定裕度按Taylor级数表示成以两群驱动功率和负荷功率为变量、以稳定裕度灵敏度为系数的线性化表达式,以求取稳定裕度的概率数字特征;再利用Gram-Charlier级数计算出系统暂态稳定裕度的概率分布。本发明提供了两种计算模式,其中纯解析方法能够计及系统内发电机机组驱动功率的不确定性,在理想两群失稳模式下计算精度高,具有速度快的特点;综合采用时域仿真与解析方法的混合算法能计及机组驱动功率和负荷功率的不确定性,通过时域仿真求取各项灵敏度系数,并用分段函数的方法解决由失稳模式变化导致稳定裕度或灵敏度系数突变的描述问题,适用于理想两群、非理想两群和多个失稳模式的情况,具有精度高的优点。
本实施例的流程图如图1所示,本实施将本发明方法应用于IEEE标准3机9节点系统,其系统接线图如图2所示,在已知发电机G2驱动功率概率分布的前提下,计算系统在线路8-9三相短路故障下的稳定裕度概率分布。为便于说明问题,假设系统中仅发电机G2的驱动功率具有随机性。
(1)令发电机G2的驱动功率等于其驱动功率概率分布均值,运用EEAC对该稳态运行工况进行详细时域仿真,得到系统的失稳模式(发电机G2为临界群,发电机G1,G3为余下群),以及系统稳定裕度初始值η0,本实施例分别用基于纯解析法的暂态稳定裕度概率分布计算和综合时域仿真与解析法的暂态稳定裕度概率分布计算两种计算模式进行计算说明。
基于纯解析法的暂态稳定裕度概率分布计算包括如下步骤:
(2)由SEEAC中灵敏度系数求解公式,求出系统稳定裕度对临界群驱动功率的一阶灵敏度系数和二阶灵敏度系数
(3)通过发电机G2驱动功率的概率分布,求取G2驱动功率的各阶半不变量γυ(Pm2)。由于该失稳模式下临界群仅含有G2一台发电机,因此临界群驱动功率的各阶半不变量即等于G2驱动功率的各阶半不变量γυ(Pms)=γυ(Pm2);
(4)按照式(1)计算稳定裕度的方差,即2阶中心距,式中E(Δ2Pms)为临界群驱动功率的方差,在本例中即为发电机G2驱动功率的方差;概率分布的数值计算表明,在Gram-Charlier级数中,越是高阶的半不变量,对最终累积分布函数的影响就越小,为了减少计算量,忽略6阶及以上的各阶半不变量,按照式(3)计算系统稳定裕度的3,4,5阶半不变量,然后通过半不变量和中心距的数学转换公式求出稳定裕度的3,4,5阶中心距;
(5)将稳定裕度的初始值η0及其各阶中心矩βυ(η)代入Gram-Charlier级数计算公式,获得稳定裕度η的累积分布函数,完成基于纯解析法的暂态稳定裕度概率分布的快速计算。结果如附图3所示,横坐标为系统稳定裕度,纵坐标为概率。由图可知系统稳定裕度不大于横坐标所对应稳定裕度值的概率,例如A点表示系统在设定故障下稳定裕度不大于40%的概率为0.7。
综合时域仿真与解析法的暂态稳定裕度概率分布计算包括如下步骤:
(6)由于本实施例中仅有发电机G2的驱动功率具有不确定性,其余发电机驱动功率及负荷功率是确定的,因此只需要根据发电机G2驱动功率的概率分布即可计算出步骤(1)初始分群的临界群驱动功率变化范围的上下限。依据人工设定的功率步长,分别按照逐渐增加G2的驱动功率和逐渐减少G2的驱动功率这两个方向调整运行方式,形成新的运行工况,并对每一个新的运行工况进行详细时域仿真。
得到附图4所示的系统稳定裕度随临界群驱动功率的变化曲线,图中ηm和η′m均在时域仿真时获得。由图可知:当发电机G2驱动功率在0到P′ms之间时,系统失稳模式的临界群仅包含发电机G2;当发电机G2驱动功率大于P′ms时,系统失稳模式发生变化,临界群中包含发电机G2和G3。此时,系统分群情况发生变化:由临界群仅包含发电机G2、余下群包含发电机G1和发电机G3,变成了临界群中包含发电机G2和G3、余下群仅包含发电机G1。图中曲线不连续的原因是,当失稳模式发生变化后,系统的临界群驱动功率变为G2和G3驱动功率之和;
(7)分别针对上述2个失稳模式,通过时域仿真摄动求取稳定裕度对临界群驱动功率Pms的一阶灵敏度系数和由于本实施例中负荷功率及余下群驱动功率是确定的,因此不必要计算和
(8)针对失稳模式1(临界群包含发电机G2,余下群包含发电机G1和G3),根据步骤1),稳定裕度初始值η′0(1)=η0;针对失稳模式2(临界群包含发电机G2和G3,余下群包含发电机G1),根据式(5)求解稳定裕度初始值η′0(2);
(9)针对每个失稳模式,计算临界群驱动功率的各阶半不变量。第一个失稳模式下,γυ(Pms)=γυ(Pm2);第二个失稳模式下,由于G3驱动功率恒定,所以γυ(Pms)=γυ(Pm2)依旧成立;
(10)针对每个失稳模式,根据式(4)求取稳定裕度的各阶半不变量,继而根据半不变量和中心距的数学转换关系,求取稳定裕度的各阶中心距;
(11)针对每个失稳模式,将稳定裕度的各阶中心矩和稳定裕度初始值代入Gram-Charlier级数计算公式,分别得到两个失稳模式下的稳定裕度概率分布曲线。由于系统失稳模式变化,临界群包含机组不同,需要结合两个失稳模式下的稳定裕度概率分布,通过分段函数综合描述系统的稳定裕度概率分布,如图5所示。至此,完成了综合采用时域仿真与解析方法的暂态稳定裕度概率分布计算。
本实施例的字母含义说明:
η:表示受扰系统的稳定裕度,其中η0为稳定裕度初始值;
S:表示临界群;A:表示余下群;
Pms:表示临界群的驱动功率,其中Pmk表示临界群内每台发电机的驱动功率;分别表示稳定裕度对临界群驱动功率的一阶灵敏度系数和二阶灵敏度系数;
Pls:表示临界群的负荷功率;
Pma:表示余下群的驱动功率;Pla:表示余下群的负荷功率;
表示第i个失稳模式下稳定裕度对临界群驱动功率的一阶灵敏度系数;
表示第i个失稳模式下稳定裕度对临界群负荷功率的一阶灵敏度系数;
表示第i个失稳模式下稳定裕度对余下群驱动功率的一阶灵敏度系数;
表示第i个失稳模式下稳定裕度对余下群负荷功率的一阶灵敏度系数;
γυ(Pms):表示临界群驱动功率的各阶半不变量,其中γυ(Pmk)表示每台发电机驱动功率的各阶半不变量;
βυ(η):表示受扰系统的稳定裕度的各阶中心矩;
γυ(η):表示受扰系统的稳定裕度的各阶半不变量;
η′(i):表示所述步骤(6)中,随着节点功率的变化,系统由第i-1个失稳模式变化成第i个失稳模式后,新失稳模式对应运行工况下的系统稳定裕度;
η′0:表示仅有一个失稳模式时,系统中各发电机驱动功率、负荷功率等于其概率分布均值时的系统稳定裕度,亦称为该失稳模式下的稳定裕度初始值;η′0(i)表示有多个失稳模式时,按照第i个失稳模式对系统进行分群,系统中各发电机驱动功率、负荷功率等于其概率分布均值时的系统稳定裕度,在本发明中亦称为第i个失稳模式下的稳定裕度初始值;
Pms0(i):表示随着节点功率的变化,系统由第i-1个失稳模式变化为第i个失稳模式后,新失稳模式对应运行工况下的系统临界群驱动功率;
Pma0(i):表示随着节点功率的变化,系统由第i-1个失稳模式变化为第i个失稳模式后,新失稳模式对应运行工况下的系统临界群负荷功率;
Pls0(i):表示随着节点功率的变化,系统由第i-1个失稳模式变化为第i个失稳模式后,新失稳模式对应运行工况下的系统余下群驱动功率;
Pla0(i):表示随着节点功率的变化,系统由第i-1个失稳模式变化为第i个失稳模式后,新失稳模式对应运行工况下的系统余下群负荷功率;
γυ:表示第ν阶半不变量;
D(η):表示受扰系统稳定裕度的方差(2阶中心矩);
E(Δ2Pms):表示临界群驱动功率的方差;
分别表示系统稳定裕度对临界群驱动功率、临界群负荷功率、余下群驱动功率和余下群负荷功率的灵敏度系数的ν次方;
分别表示第i个失稳模式下,系统稳定裕度对临界群驱动功率、临界群负荷功率、余下群驱动功率和余下群负荷功率的灵敏度系数;
ΔPms(i)、ΔPls(i)、ΔPma(i)、ΔPla(i):表示第i个失稳模式下,Pms0(i)Pls0(i)、Pma0(i)和Pla0(i)与系统中各节点功率等于其概率分布均值时的临界群驱动功率临界群负荷功率余下群驱动功率和余下群负荷功率的偏差;
P′ms:表示发电机G2驱动功率变化导致系统失稳模式改变时驱动功率的临界值,P′ms的具体取值要根据实施例的情况确定。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求范围当中。
Claims (17)
1.一种暂态稳定裕度概率分布的计算方法,其特征在于,所述方法包括如下步骤:
(1)确定初始状态下受扰系统的失稳模式和稳定裕度初始值η0,并由用户确定采用的计算模式,包括基于纯解析法的暂态稳定裕度概率分布计算和综合时域仿真与解析法的暂态稳定裕度概率分布计算;
所述基于纯解析法的暂态稳定裕度概率分布计算包括如下步骤:
(2)计算所述受扰系统稳定裕度对临界群S驱动功率的灵敏度系数I;
(3)计算所述临界群S驱动功率的概率数字特征;
(4)计算所述受扰系统稳定裕度的各阶中心矩;
(5)计算所述受扰系统稳定裕度的累积分布函数;
所述综合时域仿真和解析法的暂态稳定裕度概率分布计算包括如下步骤:
(6)计算确定节点功率变化范围内的失稳模式变化情况;
(7)对所述各失稳模式计算稳定裕度对临界群驱动功率Pms、临界群负荷功率Pls、余下群驱动功率Pma、余下群负荷功率Pla的灵敏度系数II;
(8)对所述各失稳模式计算该失稳模式下稳定裕度线性变化对应的各节点功率等于其概率分布均值时的稳定裕度初始值η′0(i);
(9)对所述失稳模式计算临界群驱动功率的概率数字特征、临界群负荷功率的概率数字特征、余下群驱动功率的概率数字特征和余下群负荷功率的概率数字特征;
(10)对所述各失稳模式计算稳定裕度的各阶中心矩;
(11)通过分段函数求取稳定裕度的累积分布函数。
2.如权利要求1所述的计算方法,其特征在于,所述概率数字特征即各阶半不变量γυ。
3.如权利要求1所述的计算方法,其特征在于,所述步骤(1)是令所述受扰系统中各节点功率等于其概率分布均值,对所述受扰系统中的设定故障,利用详细时域仿真和扩展等面积准则,得到所述受扰系统的失稳模式和稳定裕度初始值η0。
4.如权利要求1所述的计算方法,其特征在于,所述步骤(2)是令所述临界群S的驱动功率Pms和余下群A的驱动功率Pma的初始值分别等于各自群内机组驱动功率均值的总和,通过静态扩展等面积准则算法求取所述受扰系统的稳定裕度对所述临界群S驱动功率Pms的灵敏度系数I;所述灵敏度系数I包括一阶灵敏度系数和二阶灵敏度系数
5.如权利要求1所述的计算方法,其特征在于,所述步骤(3)是根据所述临界群S内每台发电机组驱动功率Pmk的概率分布,求取每台发电机驱动功率Pmk的各阶半不变量γυ(Pmk),按所述临界群S中发电机组个数对所述各阶半不变量γυ(Pmk)求和;求取Pms的各阶半不变量γυ(Pms),为
其中,k为机组编号,υ表示半不变量的阶数,Ns为临界群机组个数。
6.如权利要求1所述的计算方法,其特征在于,所述步骤(4)是用所述灵敏度系数I和所述各阶半不变量γυ(Pms),计算所述受扰系统稳定裕度的各阶半不变量γυ(η),通过半不变量和中心矩的标准数学转换公式,求取所述受扰系统稳定裕度的各阶中心矩βυ(η)。
7.如权利要求1所述的计算方法,其特征在于,所述步骤(5)是将所述受扰系统稳定裕度的初始值η0及其各阶中心矩βυ(η)代入Gram‐Charlier级数计算公式,得到所述受扰系统稳定裕度η的累积分布函数。
8.如权利要求1所述的计算方法,其特征在于,所述步骤(6)是根据所述受扰系统中具有功率不确定性的发电机和负荷功率变化范围的上下限,计算出对应步骤(1)初始状态下受扰系统失稳模式的临界群驱动功率Pms、余下群驱动功率Pma、临界群负荷功率Pls、余下群负荷功率Pla变化范围的上下限,分别按照增加所述临界群S驱动功率和余下群A负荷功率,同时减少余下群A驱动功率和临界群S负荷功率,以及减少所述临界群S驱动功率和余下群A负荷功率,同时 增加余下群A驱动功率和临界群S负荷功率两个方向,依据设定的功率变化步长及其在各机组和负荷之间的分配规则,生成一组系统稳态运行工况,再针对每一个系统稳态运行工况,进行针对设定故障下的时域仿真,识别失稳模式变化,记录每一个新的失稳模式i所对应的稳定裕度η′(i),以及临界群S的驱动功率Pms0(i)、临界群S的负荷功率Pls0(i)、余下群A的驱动功率Pma0(i)和余下群A的负荷功率Pla0(i)。
9.如权利要求1所述的计算方法,其特征在于,所述步骤(7)是对每个失稳模式,采用摄动法,通过时域仿真和EEAC求取稳定裕度对所述临界群驱动功率Pms的灵敏度系数、临界群负荷功率Pls的灵敏度系数、余下群驱动功率Pma的灵敏度系数和余下群负荷功率Pla的灵敏度系数
10.如权利要求1所述的计算方法,其特征在于,所述步骤(8)中若所述受扰系统的失稳模式的个数为1,则稳定裕度初始值η’0(i)=η0;若所述受扰系统的失稳模式有变化,则根据稳定裕度的线性化展开式,结合步骤(7)得到的灵敏度系数II和步骤(6)得到的失稳模式变化后稳定裕度η’(i),利用数学反推获得每个失稳模式下的稳定裕度初始值η’0(i)。
11.如权利要求1所述的计算方法,其特征在于,所述步骤(9)是根据所述受扰系统内发电机组驱动功率和负荷功率的概率分布,求取单个节点的发电机组驱动功率和负荷功率的各阶半不变量γυ,再根据每个失稳模式的分群情况,分别按照所述临界群S及余下群A内机组和负荷节点的个数分类求和,对应得到所述临界群驱动功率Pms的各阶半不变量、临界群负荷功率Pls的各阶半不变量、余下群驱动功率Pma的各阶半不变量和余下群负荷功率Pla的各阶半不变量。
12.如权利要求1所述的计算方法,其特征在于,所述步骤(10)是用步骤(7)得到所述灵敏度系数II和步骤(9)得到的所述临界群驱动功率Pms的各阶半不变量、临界群负荷功率Pls的各阶半不变量、余下群驱动功率Pma的各阶半不 变量和余下群负荷功率Pla的各阶半不变量,计算所述稳定裕度η的各阶半不变量,进而通过半不变量和中心矩的标准数学转换公式,求取稳定裕度η的各阶中心矩βυ(η)。
13.如权利要求1所述的计算方法,其特征在于,所述步骤(11)是针对每个失稳模式,将其稳定裕度的各阶中心矩βυ(η)和稳定裕度初始值η′0(i)代入Gram‐Charlier级数计算公式,获得各失稳模式下稳定裕度η的累积分布函数,并利用分段函数综合描述。
14.如权利要求1所述的计算方法,其特征在于,所述步骤(4)中,在计算稳定裕度η的方差D(η)时,表达式如下:
“D(η)”为受扰系统稳定裕度的方差(2阶中心矩),为表示稳定裕度对临界群驱动功率的一阶灵敏度系数的平方,表示临界群驱动功率的方差,“η”表示受扰系统的稳定裕度,其中“η0”为稳定裕度初始值。
15.如权利要求1所述的计算方法,其特征在于,所述步骤(4)中,按式(2)求取稳定裕度η对临界群驱动功率Pms的一阶灵敏度系数按式(3)求取稳定裕度的各阶半不变量γυ(η);
表示求偏导符号,“d”表示求导符号,“P”表示功率,“Pms”表示临界群的驱动功率, 分别表示稳定裕度对临界群驱动功率的一阶灵敏度系数的平方、二阶灵敏度系数,“Pma”表示余下群的驱动功率,γυ(Pms)”表示临界群驱动功率的各阶半不变量,“γυ(Pmk)”表示每台发电机驱动功率的各阶半不变量,“γυ(η)”表示受扰系统的稳定裕 度的各阶半不变量 ,“γυ”表示第V阶半不变量。
16.如权利要求1所述的计算方法,其特征在于,所述步骤(10)中,按式(4)根据时域仿真和EEAC获得的灵敏度系数灵敏度系数灵敏度系数和灵敏度系数求取稳定裕度η的各阶半不变量γυ(η);
(4)
“Pms(i)”表示第i个失稳模式下临界群驱动功率“γυ(η)”表示受扰系统稳定裕度的各阶半不变量,分别表示系统稳定裕度对临界群驱动功率、临界群负荷功率、余下群驱动功率和余下群负荷功率的灵敏度系数的ν次方。
17.如权利要求1所述的计算方法,其特征在于,所述步骤(8)中,对每个失稳模式i,根据稳定裕度的线性化展开关系,按式(5)进行数学反推,求取各节点功率等于其概率分布均值时的稳定裕度初始值η′0(i);
式中ΔPms(i)=Pms0(i)-Pms,ΔPls(i)=Pls0(i)-Pls,ΔPma(i)=Pma0(i)-Pma,ΔPla(i)=Pla0(i)-Pla,Pms、Pls、Pma和Pla分别为各节点功率等于其概率分布均值时的临界群驱动功率、临界群负荷功率、余下群驱动功率和余下群负荷功率,“η′(i)”表示随着节点功率的变化系统由第i-1个失稳模式变化成第i个失稳模式后,新失稳模式对应运行工况下的系统稳定裕度,“Pms0(i)”、“Pma0(i)”、“Pls0(i)”和“Pla0(i)”分别表示随着节点功率的变化,系统由第i-1个失稳模式变化为第i个失稳模式后,新失稳模式对应运行工况下的系统临界群驱动功率、临界群负荷功率、余下群驱动功率和余下群负荷功率,“η′0(i)”表示有多个失稳模式时,按照第i个失稳模式对系统进行分群,系统中各发电机 驱动功率、负荷功率等于其概率分布均值时的系统稳定裕度,在本发明中亦称为第i个失稳模式下的稳定裕度初始值,分别表示在第i个失稳模式下,临界群负荷功率、余下群驱动功率和余下群负荷功率的灵敏度系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210263026.6A CN102790389B (zh) | 2012-07-27 | 2012-07-27 | 一种暂态稳定裕度概率分布的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210263026.6A CN102790389B (zh) | 2012-07-27 | 2012-07-27 | 一种暂态稳定裕度概率分布的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102790389A CN102790389A (zh) | 2012-11-21 |
CN102790389B true CN102790389B (zh) | 2014-09-10 |
Family
ID=47155716
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210263026.6A Active CN102790389B (zh) | 2012-07-27 | 2012-07-27 | 一种暂态稳定裕度概率分布的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102790389B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103515951B (zh) * | 2013-08-12 | 2015-07-08 | 南京南瑞集团公司 | 一种电力系统直流紧急功率调制和切负荷暂态功角稳定控制敏感因子计算方法 |
CN103762579B (zh) * | 2013-12-22 | 2015-11-11 | 国家电网公司 | 一种提高电网暂态功角稳定性的直流功率紧急控制性能指标计算方法 |
CN105160586B (zh) * | 2015-02-17 | 2019-01-18 | 国网电力科学研究院 | 电力系统失稳模式易变性的快速评估方法 |
CN105656038B (zh) * | 2016-03-15 | 2018-02-13 | 国电南瑞科技股份有限公司 | 一种计及电网接纳能力的概率潮流计算方法 |
CN106992516B (zh) * | 2017-04-18 | 2020-02-14 | 国网上海市电力公司 | 概率静态电压稳定裕度优化实现方法 |
CN107706948B (zh) * | 2017-11-03 | 2020-12-01 | 华北电力大学 | 多维阶数控制的多步Taylor级数暂态稳定分析方法 |
CN114253157B (zh) * | 2021-12-21 | 2024-05-14 | 华中科技大学 | 一种基于二阶灵敏度分析的电机多参数优化方法和系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1835336A (zh) * | 2006-03-10 | 2006-09-20 | 天津大学 | 基于概率的电力系统电压稳定在线监视和预防控制方法 |
-
2012
- 2012-07-27 CN CN201210263026.6A patent/CN102790389B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1835336A (zh) * | 2006-03-10 | 2006-09-20 | 天津大学 | 基于概率的电力系统电压稳定在线监视和预防控制方法 |
Non-Patent Citations (8)
Title |
---|
兰强等.基于EEAC的考虑暂态安全稳定约束的最优潮流计算.《电力系统自动化》.2010,第34卷(第8期),34-38,115. |
关于暂态稳定不确定性分析的评述;薛禹胜等;《电力系统自动化》;20070725;第31卷(第14期);1-6,75 * |
刘强等.基于稳定域及条件概率的暂态稳定不确定性分析.《电力系统自动化》.2007,第31卷(第19期),1-6. |
基于EEAC的考虑暂态安全稳定约束的最优潮流计算;兰强等;《电力系统自动化》;20100425;第34卷(第8期);34-38,115 * |
基于动态安全域与埃奇沃斯级数的电力系统暂态稳定概率评估;汪隆君;《中国电机工程学报》;20110105;第31卷(第1期);52-58 * |
基于稳定域及条件概率的暂态稳定不确定性分析;刘强等;《电力系统自动化》;20071010;第31卷(第19期);1-6 * |
汪隆君.基于动态安全域与埃奇沃斯级数的电力系统暂态稳定概率评估.《中国电机工程学报》.2011,第31卷(第1期),52-58. |
薛禹胜等.关于暂态稳定不确定性分析的评述.《电力系统自动化》.2007,第31卷(第14期),1-6,75. |
Also Published As
Publication number | Publication date |
---|---|
CN102790389A (zh) | 2012-11-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102790389B (zh) | 一种暂态稳定裕度概率分布的计算方法 | |
Wen et al. | A review on reliability assessment for wind power | |
CN103746370B (zh) | 一种风电场可靠性建模方法 | |
Feng et al. | Developing a successive linear programming model for head-sensitive hydropower system operation considering power shortage aspect | |
CN115693757A (zh) | 一种基于数字孪生技术的光伏能源优化调控方法 | |
CN103700036B (zh) | 一种适于电力系统多时间尺度的暂态稳定性投影积分方法 | |
CN102509027B (zh) | 一种基于交叉熵理论的风电功率组合预测方法 | |
CN104269867A (zh) | 一种节点扰动功率转移分布均衡度分析方法 | |
CN109636035A (zh) | 负荷预测模型创建方法及装置、电力负荷预测方法及装置 | |
CN107730111A (zh) | 一种考虑用户负荷和新能源接入的配网电压风险评估模型 | |
Fu et al. | Use of a second-order reliability method to estimate the failure probability of an integrated energy system | |
CN106655177A (zh) | 基于扩展二阶锥规划的分布式电源最大接入能力计算方法 | |
CN105226650A (zh) | 基于微燃机-储能联合运行策略的微电网可靠性计算方法 | |
CN103473393A (zh) | 一种考虑随机概率的输电裕度控制模型建模方法 | |
CN103455716A (zh) | 一种基于超短期风电功率预测的电力系统电压稳定裕度计算方法 | |
Liu et al. | A decomposition-coordination planning method for flexible generation resources in isolated microgrids | |
Dimopoulou et al. | A Markov decision process for managing a hybrid energy storage system | |
CN106786608A (zh) | 一种适用于分布式电源接入的不确定谐波潮流计算方法 | |
Deng et al. | Economic control for a residential photovoltaic-battery system by combining stochastic model predictive control and improved correction strategy | |
Zhang et al. | Frequency-constrained unit commitment for power systems with high renewable energy penetration | |
CN105633950A (zh) | 考虑风电注入不确定性的多目标随机模糊动态最优潮流 | |
CN106056305A (zh) | 一种基于状态聚类的发电系统可靠性快速评估方法 | |
CN109711612A (zh) | 一种优化回声状态网络的风力发电功率预测方法及装置 | |
Bai et al. | Probabilistic production simulation of a wind/photovoltaic/energy storage hybrid power system based on sequence operation theory | |
Zhang et al. | An efficient approach for robust SCOPF considering load and renewable power uncertainties |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |