CN105140909A - 基于启发式潮流计算的发电机出力灵敏度计算方法 - Google Patents
基于启发式潮流计算的发电机出力灵敏度计算方法 Download PDFInfo
- Publication number
- CN105140909A CN105140909A CN201510458158.8A CN201510458158A CN105140909A CN 105140909 A CN105140909 A CN 105140909A CN 201510458158 A CN201510458158 A CN 201510458158A CN 105140909 A CN105140909 A CN 105140909A
- Authority
- CN
- China
- Prior art keywords
- idle
- meritorious
- sensitivity
- node
- generating set
- 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.)
- Granted
Links
- 230000035945 sensitivity Effects 0.000 title claims abstract description 53
- 238000004364 calculation method Methods 0.000 title claims abstract description 35
- 238000000034 method Methods 0.000 claims abstract description 8
- 230000001105 regulatory effect Effects 0.000 claims description 9
- 150000001875 compounds Chemical class 0.000 claims description 6
- 238000000205 computational method Methods 0.000 claims description 4
- 238000009795 derivation Methods 0.000 abstract description 2
- 238000005457 optimization Methods 0.000 description 6
- 238000004422 calculation algorithm Methods 0.000 description 3
- 230000001276 controlling effect Effects 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 239000003990 capacitor Substances 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000005070 ripening Effects 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于启发式潮流计算的发电机出力灵敏度计算方法。该方法通过将全部PV节点转化为PQ节点,令被研究发电机组单独增加或减少适当出力,然后计算两个状态断面的潮流,对关注的响应量求取两个潮流断面的结果的差值,由此可计算出,某台发电机有功或无功出力对于电网中其他线路潮流、线损、网损、节点电压的灵敏度系数。解决了常规灵敏度系数计算方法中,公式推导复杂,计算量大,容易出错,编程求解开发工作量大的问题。
Description
技术领域
本发明属于电力系统计算领域,特别涉及一种发电机出力灵敏度的计算方法。
背景技术
发电机有功和无功出力变化对电力系统中各支路潮流、节点电压以及网损变化的影响是进行电力系统有功和无功优化和控制的依据。传统的方法是建立全网的潮流方程,并将潮流方程组针对某台发电机的有功出力P或无功出力Q求导数,然后通过求解非线性方程,求解出某一响应量对发电机P或Q的导数。然而,实际上由于线路、变压器数目庞大,其响应量受全网所有发电机影响,而发电机数目往往也很多,因此采用上述传统方法要针对不同发电机的有功和无功建立大量的导数方程,非常容易出错,并且需要针对每台发电机的每个变量求解涉及全网的非线性方程,其计算求解量也很大。
此外,对于全网网损、区域网损等响应量,需要在常规潮流方程的计算上,进一步求出各支路的网损,并对各支路网损进行总加得出全网网损或区域网损公式,然后对总加后的网损求导,由于不是显示解,并受全网所有机组的影响,需要将网损的求导方程再与潮流方程组的求导方程进行联立求解。上述,求导工作量很大,并且极易出错,相应的非线性方程组求解的程序开发工作量也很大,因此发电机出力的灵敏度系数虽然在优化和控制中应该能起很重要的作用,但是实际上由于程序开发实现的难度,很少见有实际的应用。
本申请提出的方法,利用常规成熟的,并且已模块化的潮流计算程序,依据灵敏度的原始定义,通过单独改变一个变量,并保持其他注入量的不变,从而实现发电机灵敏度的方便快速计算。从而为发电机出力灵敏度在优化和控制中的广泛应用提供了基础。
发明内容
为了解决常规的发电机出力灵敏度系数计算方法中公式推导复杂,计算量大,容易出错,编程求解开发工作量大的问题。本申请提出了基于启发式潮流计算的发电机出力灵敏度计算方法。
本发明具体采用以下技术方案:
一种基于启发式潮流计算的发电机出力灵敏度计算方法,其特征在于:该方法通过将全部PV节点(即有功出力和节点电压固定的发电机节点)转化为PQ节点(即有功出力和无功出力固定的节点),令待计算灵敏度的机组单独增加或减少适当出力,然后计算改变出力前后两个状态断面的潮流,对关注的响应量求取两个潮流断面的结果的差值,由此可计算出,某台发电机有功或无功出力对于电网中其他线路潮流、线损、网损、节点电压的灵敏度系数。
基于启发式潮流计算的发电机出力灵敏度计算方法,其特征在于,所述计算方法包括以下步骤:
步骤1:获取需要计算灵敏度的发电机组所处的当前潮流断面参数和状态数据,以及该机组的有功或无功出力变量及其当前值;
步骤2:在需要计算灵敏度的发电机组之外,选取容量大于本电网参与自动发电控制的各发电机组容量的外网等值机组做为松弛节点,或者选取本电网中参与自动发电控制的发电机组中容量最大的1个作为松弛节点,其余发电机节点作为PV节点,负荷节点作为PQ节点,计算对应的全网潮流值;
步骤3:以步骤2中的潮流计算结果为基础,将全部PV节点转化为PQ节点,将需要计算灵敏度的发电机组i的有功pi或无功qi变化小的调节量适当值Δpi或Δqi,保持步骤2中的松弛节点的电压幅值、相角数值和PQ节点的有功和无功数值不变,计算对应的全网潮流值;
其中Δpi或Δqi的选取,应以不引起本电网中各支路、变压器以及发电机组有功、无功和电压越限为前提;
步骤4:计算各响应量即各支路有功、各支路无功、各支路损耗、组合断面有功、组合断面无功、总网损、区域网损、节点电压在前后两个潮流断面,即步骤2和步骤3两个潮流断面的差值,即变化量;
步骤5:用Δpi或Δqi除各响应量在两个断面间的变化量,由此得到机组i的有功或无功对各响应量的灵敏度系数。
本申请进一步优选以下方案:
在步骤3中,Δpi或Δqi选取为待计算灵敏度的发电机组额定容量的1%。
将需要计算灵敏度的发电机组i的有功pi或无功qi变化小的调节量适当值Δpi或Δqi时,优先采用增加待计算灵敏度的发电机组有功或无功出力的变化方式,若增加的有功Δpi或无功Δqi导致电网中任何一个设备的有功、无功或电压的越限,则采用减少Δpi或Δqi的变化方式;若增减Δpi或Δqi均导致前述越限情况的发生,则将调节量适当值减半,再继续按先增加后减少的方式寻找变化方向,直到找到的调节量适当值不引起本电网中各支路、变压器以及发电机组的有功、无功和电压越限为止。
上述方法,基于常规、成熟的潮流计算程序模块,依据灵敏度的原始定义,实现了发电机出力灵敏度的方便快速计算,降低了开发的工作量,从而为发电机出力灵敏度在优化和控制中的广泛应用提供了基础。
附图说明
图1基于启发式潮流计算的发电机出力灵敏度计算方法
具体实施方式
下面结合说明书附图和具体实施例对本发明的技术方案做进一步详细介绍。
如图1所示为基于启发式潮流计算的发电机出力灵敏度计算方法流程图。依据本发明所提出的基于启发式潮流计算的发电机出力灵敏度计算方法所开发的灵敏度计算模块运行于电力调度控制中心的SCADA/EMS平台,并可被其他优化控制应用的算法所调用,例如有功优化和无功优化应用。对于同一运行断面的不同发电机出力需要分别调用本算法,无功和有功出力的灵敏度也应分别计算。当运行断面改变时,通常应该重新计算上述灵敏度,但是在实际应用中,当网络结构没有变化,并且运行点即各机组和负荷的功率变化不大时,可以仍然使用之前计算出的灵敏度来做优化或控制决策,从而减少计算量。本发明所提出的优化算法的具体实施步骤如下:
步骤1:获取被研究机组所处的当前潮流断面参数(例如:线路、变压器、电容器、电抗器阻抗,变压器分接头档位设置、发电机出力上下限)和状态数据(PV发电机有功和电压幅值、Vθ发电机电压幅值和相角、负荷有功和无功),以及需要计算灵敏度的机组和相应的有功或无功变量;
步骤2:在需要计算灵敏度的发电机组之外,首先选取容量大于本电网参与自动发电控制的各发电机组容量的外网等值机组作为松弛节点,若无外网等值机组,则从网内参与自动发电控制的发电机组中选取容量最大的1个作为松弛节点,其余发电机节点作为PV节点,负荷节点作为PQ节点,针对这一状态断面计算对应的全网潮流值,确定各母线电压幅值和相角、各支路有功无功和线损,以及全网的网损;
步骤3:以步骤2中的潮流计算结果为基础,将全部PV节点转化为PQ节点,将待计算灵敏度的发电机组i的有功pi或无功qi变化小的适当值Δpi或Δqi,保持步骤2中的松弛节点的电压幅值、相角数值和PQ节点的有功和无功数值不变,计算该断面对应的全网潮流值;
其中Δpi或Δqi的选取,应以不引起电网中各线路、变压器以及发电机有功、无功和电压越限为限;
步骤4:计算各响应量即各支路有功、各支路无功、各支路损耗、组合断面有功、组合断面无功、总网损、区域网损、节点电压在前后两个潮流断面,即步骤2和步骤3两个潮流断面的差值,即变化量;
步骤5:用Δpi或Δqi除各响应量在两个断面间的变化量,由此得到机组i的有功或无功对各响应量的灵敏度系数。
其中,在步骤3中,Δpi或Δqi选取为待计算灵敏度的发电机组额定容量的1%。
将待计算灵敏度的发电机组的有功pi或无功qi变化小的适当值Δpi或Δqi时,优先采用增加待计算灵敏度的发电机组有功或无功出力的变化方式,若增加的有功Δpi或无功Δqi导致电网中任何一个设备的有功、无功或电压的越限,则采用减少Δpi或Δqi的变化方式;若增减Δpi或Δqi均导致前述越限情况的发生,则将调节量Δpi或Δqi减半,再继续按前述方法寻找变化方向,直到找到的Δpi或Δqi不引起电网中任何一个设备的有功、无功和电压越限为止。
申请人结合说明书附图对本发明的实施例做了详细的说明与描述,但是本领域技术人员应该理解,以上实施例仅为本发明的优选实施方案,详尽的说明只是为了帮助读者更好地理解本发明精神,而并非对本发明保护范围的限制,相反,任何基于本发明的发明精神所作的任何改进或修饰都应当落在本发明的保护范围之内。
Claims (4)
1.一种基于启发式潮流计算的发电机出力灵敏度计算方法,其特征在于:
该方法通过将全部PV节点即有功出力和节点电压固定的发电机节点转化为PQ节点即有功出力和无功出力固定的节点,令待计算灵敏度的机组单独增加或减少适当出力,然后计算改变出力前后两个状态断面的潮流,对关注的响应量求取两个潮流断面的结果的差值,由此可计算出,某台发电机有功或无功出力对于电网中其他线路潮流、线损、网损、节点电压的灵敏度系数。
2.根据权利要求1所述的基于启发式潮流计算的发电机出力灵敏度计算方法,其特征在于,所述计算方法包括以下步骤:
步骤1:获取需要计算灵敏度的发电机组所处的当前潮流断面参数和状态数据,以及该机组的有功或无功出力变量及其当前值;
步骤2:在需要计算灵敏度的发电机组之外,选取容量大于本电网参与自动发电控制的各发电机组容量的外网等值机组做为松弛节点,或者选取本电网中参与自动发电控制的发电机组中容量最大的1个作为松弛节点,其余发电机节点作为PV节点,负荷节点作为PQ节点,计算对应的全网潮流值;
步骤3:以步骤2中的潮流计算结果为基础,将全部PV节点转化为PQ节点,将需要计算灵敏度的发电机组i的有功pi或无功qi变化小的调节量适当值Δpi或Δqi,保持步骤2中的松弛节点的电压幅值、相角数值和PQ节点的有功和无功数值不变,计算对应的全网潮流值;
其中Δpi或Δqi的选取,应以不引起本电网中各支路、变压器以及发电机组有功、无功和电压越限为前提;
步骤4:计算各响应量即各支路有功、各支路无功、各支路损耗、组合断面有功、组合断面无功、总网损、区域网损、节点电压在前后两个潮流断面,即步骤2和步骤3两个潮流断面的差值,即变化量;
步骤5:用Δpi或Δqi除各响应量在两个断面间的变化量,由此得到机组i的有功或无功对各响应量的灵敏度系数。
3.根据权利要求2所述的基于启发式潮流计算的发电机出力灵敏度计算方法,其特征在于:
在步骤3中,Δpi或Δqi选取为待计算灵敏度的发电机组额定容量的1%。
4.根据权利要求3所述的基于启发式潮流计算的发电机出力灵敏度计算方法,其特征在于:
将需要计算灵敏度的发电机组i的有功pi或无功qi变化小的调节量适当值Δpi或Δqi时,优先采用增加待计算灵敏度的发电机组有功或无功出力的变化方式,若增加的有功Δpi或无功Δqi导致电网中任何一个设备的有功、无功或电压的越限,则采用减少Δpi或Δqi的变化方式;若增减Δpi或Δqi均导致前述越限情况的发生,则将调节量适当值减半,再继续按先增加后减少的方式寻找变化方向,直到找到的调节量适当值不引起本电网中各支路、变压器以及发电机组的有功、无功和电压越限为止。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510458158.8A CN105140909B (zh) | 2015-07-30 | 2015-07-30 | 基于启发式潮流计算的发电机出力灵敏度计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510458158.8A CN105140909B (zh) | 2015-07-30 | 2015-07-30 | 基于启发式潮流计算的发电机出力灵敏度计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105140909A true CN105140909A (zh) | 2015-12-09 |
CN105140909B CN105140909B (zh) | 2017-07-18 |
Family
ID=54726167
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510458158.8A Active CN105140909B (zh) | 2015-07-30 | 2015-07-30 | 基于启发式潮流计算的发电机出力灵敏度计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105140909B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106410784A (zh) * | 2016-06-02 | 2017-02-15 | 国网江西省电力公司赣东北供电分公司 | 变电站有功负荷对地区电网有功网损的灵敏度计算方法 |
CN106712009A (zh) * | 2017-01-20 | 2017-05-24 | 清华大学 | 基于分布式光储的主动配电网安全运行优化方法 |
CN107181284A (zh) * | 2017-04-24 | 2017-09-19 | 中国电力科学研究院 | 一种线路越限电量的调整方法和装置 |
CN109586346A (zh) * | 2018-12-04 | 2019-04-05 | 深圳供电局有限公司 | 一种交直流大电网动态等值的平衡机出力调整方法和系统 |
CN112436526A (zh) * | 2020-08-21 | 2021-03-02 | 广西大学 | 一种基于矢量化技术的特征值解析灵敏度计算方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101771276A (zh) * | 2009-01-04 | 2010-07-07 | 华东电网有限公司 | 一种多缓冲网自动识别方法 |
CN104134994A (zh) * | 2014-07-08 | 2014-11-05 | 国电南瑞科技股份有限公司 | 一种利用agc进行稳定断面潮流越限的校正控制方法 |
-
2015
- 2015-07-30 CN CN201510458158.8A patent/CN105140909B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101771276A (zh) * | 2009-01-04 | 2010-07-07 | 华东电网有限公司 | 一种多缓冲网自动识别方法 |
CN104134994A (zh) * | 2014-07-08 | 2014-11-05 | 国电南瑞科技股份有限公司 | 一种利用agc进行稳定断面潮流越限的校正控制方法 |
Non-Patent Citations (2)
Title |
---|
江岳文等: "基于灵敏度计算的发电机无功资源价值分析", 《继电器》 * |
肖晓春: "组合系统可靠性价值评估中有功校正策略的研究", 《中国优秀博硕士学位论文全文数据库(硕士)(工程科技II辑)》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106410784A (zh) * | 2016-06-02 | 2017-02-15 | 国网江西省电力公司赣东北供电分公司 | 变电站有功负荷对地区电网有功网损的灵敏度计算方法 |
CN106410784B (zh) * | 2016-06-02 | 2020-03-17 | 国网江西省电力公司赣东北供电分公司 | 变电站有功负荷对地区电网有功网损的灵敏度计算方法 |
CN106712009A (zh) * | 2017-01-20 | 2017-05-24 | 清华大学 | 基于分布式光储的主动配电网安全运行优化方法 |
CN106712009B (zh) * | 2017-01-20 | 2019-06-07 | 清华大学 | 基于分布式光储的主动配电网安全运行优化方法 |
CN107181284A (zh) * | 2017-04-24 | 2017-09-19 | 中国电力科学研究院 | 一种线路越限电量的调整方法和装置 |
CN107181284B (zh) * | 2017-04-24 | 2021-04-16 | 中国电力科学研究院 | 一种线路越限电量的调整方法和装置 |
CN109586346A (zh) * | 2018-12-04 | 2019-04-05 | 深圳供电局有限公司 | 一种交直流大电网动态等值的平衡机出力调整方法和系统 |
CN112436526A (zh) * | 2020-08-21 | 2021-03-02 | 广西大学 | 一种基于矢量化技术的特征值解析灵敏度计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105140909B (zh) | 2017-07-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Wang et al. | Inverter-based voltage control of distribution networks: A three-level coordinated method and power hardware-in-the-loop validation | |
CN105140909A (zh) | 基于启发式潮流计算的发电机出力灵敏度计算方法 | |
Ranamuka et al. | Flexible AC power flow control in distribution systems by coordinated control of distributed solar-PV and battery energy storage units | |
CN108134401B (zh) | 交直流混合系统多目标潮流优化及控制方法 | |
CN107017661B (zh) | 一种模态依赖的微网多态运行切换控制方法 | |
CN103020853B (zh) | 一种短期交易计划安全校核的方法 | |
CN109524979B (zh) | 一种含vsc-mtdc的交直流互联电网连续潮流模型 | |
CN102780226B (zh) | 基于斩控均压的链式statcom直流侧电压控制方法及控制电路 | |
CN103280821A (zh) | 智能配电系统多时段动态无功优化方法 | |
CN109713694A (zh) | 一种三相不对称配电网的网络动态重构方法 | |
CN103701132B (zh) | 一种基于有功-无功协调的控制主动配电网过电压的方法 | |
CN105322535A (zh) | 含统一潮流控制器的电力系统两阶段最优潮流计算方法 | |
CN101976837A (zh) | 一种电力系统潮流无解调整方法 | |
CN108493985A (zh) | 含分布式电源的配电网电压越限薄弱环节的识别方法 | |
CN105162144A (zh) | 一种配电网电压无功优化的线性逼近求解方法 | |
CN102930141A (zh) | 一种电网复杂监控断面潮流控制的快速计算方法 | |
CN108306334A (zh) | 基于粒子群优化算法的风电场内部无功优化策略 | |
CN104779609A (zh) | 一种用于互联电网的潮流协同控制方法 | |
CN115481856A (zh) | 考虑综合需求响应的综合能源系统多尺度调度方法及系统 | |
CN114899829A (zh) | 一种电压稳定性的评估方法及装置 | |
CN110350538B (zh) | 一种基于主动需求侧响应的微电网协调控制方法 | |
CN103138595A (zh) | 中点箝位式变流器的控制系统及控制方法 | |
Murugan et al. | A new approach for voltage control of IPFC and UPFC for power flow management | |
CN113241753B (zh) | 一种用于直流微电网的改进虚拟发电机控制方法 | |
CN111181189A (zh) | 一种用于多端直流输电系统的非线性下垂控制方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |