CN108281967A - 一种改善电力系统潮流计算收敛性的算法启动方法 - Google Patents
一种改善电力系统潮流计算收敛性的算法启动方法 Download PDFInfo
- Publication number
- CN108281967A CN108281967A CN201810009362.5A CN201810009362A CN108281967A CN 108281967 A CN108281967 A CN 108281967A CN 201810009362 A CN201810009362 A CN 201810009362A CN 108281967 A CN108281967 A CN 108281967A
- Authority
- CN
- China
- Prior art keywords
- node
- voltage
- nodes
- electric power
- power system
- 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
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Abstract
本发明公开了一种改善电力系统潮流计算收敛性的算法启动方法。目前的平滑启动是电力系统潮流计算中最常采用的算法启动方法,其应用于病态电力系统时,容易出现初始解偏差大,潮流计算收敛困难的情况。本发明根据各PV节点及平衡节点的电压幅值取为相应设定值,电压相位取为0,且各PQ节点的功率注入相量取为0,构造初始潮流解,用于电力系统潮流计算的算法启动。本发明能够有效减小初始潮流解的偏差,为病态电力系统潮流计算提供较好的初始潮流解,改善电力系统潮流计算的收敛性。
Description
技术领域
本发明属于电力系统潮流计算领域,特别涉及一种改善电力系统潮流计算收敛性的算法启动方法。
背景技术
潮流计算是电力系统常用的分析计算功能,其数学本质是非线性方程组的数值求解。复杂电力系统潮流计算的计算机算法均为迭代算法,数值求解时均需要提供潮流初始解,用于潮流计算的算法启动。初始潮流解的质量是影响潮流计算收敛性的重要因素。常用的电力系统潮流计算方法(如牛顿拉夫逊法、PQ分解法和最优乘子法)均对潮流计算的初值有较高的要求。
平滑启动是电力系统潮流计算中最常采用的算法启动方法,能够为良态电力系统提供较为合理的初始值,在实际电力系统中有着广泛的应用,但将其应用于病态电力系统时,容易出现初始解偏差大,潮流计算收敛困难的情况。
三绕组变压器中存在的小的负电抗支路,以及用于模拟闭合联络开关的小阻抗支路是电力系统潮流计算出现病态的重要因素。随着电力负荷的快速增长及西电东送、区域电网互联战略的实施,一个装机容量数亿千瓦、横跨数千公里的广域电力系统已经在我国形成。大规模电力系统中三绕组变压器多,小阻抗支路多,其潮流计算经常存在严重病态的情况,导致潮流计算收敛困难。为潮流计算提供偏差更小的初始潮流解是改善潮流计算收敛性的重要途径。
发明内容
本发明的目的是提供一种改善电力系统潮流计算收敛性的算法启动方法,以解决病态电力系统潮流计算的算法启动问题,有效避免因潮流初始解不合理引起的潮流计算收敛困难问题。
为实现上述目的,本发明采用以下技术方案:一种改善电力系统潮流计算收敛性的算法启动方法,包括如下步骤:
步骤1),读入电网模型及参数,进行网络拓扑分析,并生成导纳矩阵;
步骤2),各PV节点及平衡节点的电压幅值取其设定值,电压相位取为0;
步骤3),各PQ节点的功率注入相量取为0;
步骤4),根据网络方程计算各PQ节点的电压相量;
步骤5),采用前述给出的各节点电压相量作为潮流计算的初始解,用于潮流计算的算法启动。
作为上述技术方案的补充,步骤4)中,采用如下方法确定各PQ节点的电压相量:
交流电力系统的网络方程为:
其中,表示节点导纳矩阵,表示节点电压复向量,表示节点注入电流复向量;
根据节点类型对电压复向量及电流复向量进行分组,
其中,下标v表示PV节点及平衡节点对应的部分,下标q表示PQ节点对应的部分;
由上式得
根据步骤2),为给定的常量,根据步骤3),故有
求解上式的复线性方程组,即给出各PQ节点的电压相量,能够用于潮流计算的算法启动。
作为上述技术方案的补充,步骤4)中,忽略导纳矩阵的实数部分,所述的复线性方程组替换为
其中,Bqq表示的虚数部分,Bqv表示的虚数部分;
由步骤2)可知,各PV节点及平衡节点的电压相量均为实数,即为实数向量,由上式可知,也为实数向量,故有
BqqVq=-BqvVv
其中,Vq为相位取0时的实数描述,表示各PQ节点电压幅值构成的向量;Vv为相位取0时的实数描述,表示各PV节点及平衡节点的电压幅值构成的向量;
求解上式的实线性方程组获得各PQ节点的电压幅值,取各节点的电压相位为0,用于潮流计算的算法启动。
本发明所述的算法启动方法能够有效减小初始潮流解的偏差,为病态电力系统潮流计算提供较好的初始潮流解,改善电力系统潮流计算的收敛性。
附图说明
图1为本发明的流程图;
图2为本发明实施例中采用两种不同启动方法下最优乘子法潮流计算过程中误差随迭代次数的变化情况图。
具体实施方式
下面结合说明书附图和具体实施方式对本发明作进一步说明。
如图1所示的一种改善电力系统潮流计算收敛性的算法启动方法,包括如下步骤:
步骤1:读入电网模型及参数,进行网络拓扑分析,并生成导纳矩阵;
步骤2:各PV节点及平衡节点的电压幅值取其设定值,电压相位取为0;
步骤3:各PQ节点的功率注入相量取为0;
步骤4:根据网络方程计算各PQ节点的电压相量;
步骤5:采用步骤2~4给出的各节点电压相量作为潮流计算的初始解,用于潮流计算的算法启动。
步骤4中,采用如下方法确定各PQ节点的电压相量。
交流电力系统的网络方程可描述为:
其中,表示节点导纳矩阵,表示节点电压复向量,表示节点注入电流复向量。
根据节点类型对电压复向量及电流复向量进行分组,
其中,下标v表示PV节点及平衡节点对应的部分,下标q表示PQ节点对应的部分。
由式(2)可得
根据步骤2,为给定的常量,根据步骤3,故有
求解式(4)的复线性方程组,即可给出各PQ节点的电压相量,可以用于潮流计算的算法启动。
实用中,复线性方程组的计算量相对较大,忽略导纳矩阵的实数部分,式(4)可替换为
其中,Bqq表示的虚数部分,Bqv表示的虚数部分。
由步骤2可知,各PV节点及平衡节点的电压相量均为实数,即事实上为实数向量,由式(5)可知,也为实数向量,故有
BqqVq=-BqvVv(6)
其中,Vq为相位取0时的实数描述,表示各PQ节点电压幅值构成的向量;Vv为相位取0时的实数描述,表示各PV节点及平衡节点的电压幅值构成的向量;
求解式(6)的实线性方程组可获得各PQ节点的电压幅值,取各节点的电压相位为0,即可用于潮流计算的算法启动。
本发明应用时,以某34964节点实际电网为例,对本发明提出的启动法与电力系统常用的平滑启动法进行比较,图2给出了2种算法启动方法下最优乘子法潮流计算过程中误差随迭代次数的变化情况。本发明给出的启动法(即零注入启动)能够为潮流计算提供更好的初始解(初始误差为249.835),潮流计算的收敛性更好,经8次迭代即获得收敛潮流解,而传统平滑启动方法所提供的初始解较差(初始误差为1268.667),需要经过13次迭代才能收敛。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (3)
1.一种改善电力系统潮流计算收敛性的算法启动方法,其特征在于,包括如下步骤:
步骤1),读入电网模型及参数,进行网络拓扑分析,并生成导纳矩阵;
步骤2),各PV节点及平衡节点的电压幅值取其设定值,电压相位取为0;
步骤3),各PQ节点的功率注入相量取为0;
步骤4),根据网络方程计算各PQ节点的电压相量;
步骤5),采用前述给出的各节点电压相量作为潮流计算的初始解,用于潮流计算的算法启动。
2.根据权利要求1所述的改善电力系统潮流计算收敛性的算法启动方法,其特征在于,步骤4)中,采用如下方法确定各PQ节点的电压相量:
交流电力系统的网络方程为:
其中,表示节点导纳矩阵,表示节点电压复向量,表示节点注入电流复向量;
根据节点类型对电压复向量及电流复向量进行分组,
其中,下标v表示PV节点及平衡节点对应的部分,下标q表示PQ节点对应的部分;
由上式得
根据步骤2),为给定的常量,根据步骤3),故有
求解上式的复线性方程组,即给出各PQ节点的电压相量,用于潮流计算的算法启动。
3.根据权利要求2所述的改善电力系统潮流计算收敛性的算法启动方法,其特征在于,步骤4)中,忽略导纳矩阵的实数部分,所述的复线性方程组替换为
其中,Bqq表示的虚数部分,Bqv表示的虚数部分;
由步骤2)可知,各PV节点及平衡节点的电压相量均为实数,即为实数向量,由上式可知,也为实数向量,故有
BqqVq=-BqvVv
其中,Vq为相位取0时的实数描述,表示各PQ节点电压幅值构成的向量;Vv为相位取0时的实数描述,表示各PV节点及平衡节点的电压幅值构成的向量;
求解上式的实线性方程组获得各PQ节点的电压幅值,取各节点的电压相位为0,用于潮流计算的算法启动。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810009362.5A CN108281967B (zh) | 2018-01-04 | 2018-01-04 | 一种改善电力系统潮流计算收敛性的算法启动方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810009362.5A CN108281967B (zh) | 2018-01-04 | 2018-01-04 | 一种改善电力系统潮流计算收敛性的算法启动方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108281967A true CN108281967A (zh) | 2018-07-13 |
CN108281967B CN108281967B (zh) | 2020-06-26 |
Family
ID=62803124
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810009362.5A Active CN108281967B (zh) | 2018-01-04 | 2018-01-04 | 一种改善电力系统潮流计算收敛性的算法启动方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108281967B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112260283A (zh) * | 2020-10-22 | 2021-01-22 | 云南电网有限责任公司 | 一种电力系统输电断面潮流调整方法及装置 |
CN117828243A (zh) * | 2024-03-06 | 2024-04-05 | 国网上海能源互联网研究院有限公司 | 一种fpga潮流并行计算系统和方法 |
CN117828243B (zh) * | 2024-03-06 | 2024-05-14 | 国网上海能源互联网研究院有限公司 | 一种fpga潮流并行计算系统和方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130226482A1 (en) * | 2012-02-24 | 2013-08-29 | Hongbo Sun | Decoupled Three-Phase Power Flow Analysis Method for Unbalanced Power Distribution Systems |
CN106356860A (zh) * | 2016-10-17 | 2017-01-25 | 大连海事大学 | 一种配电系统三相潮流计算的电压初值设置方法 |
CN106374497A (zh) * | 2016-10-20 | 2017-02-01 | 天津大学 | 基于ami的中低压混合配电网一体化状态估计方法 |
CN107016489A (zh) * | 2017-03-09 | 2017-08-04 | 中国电力科学研究院 | 一种电力系统抗差状态估计方法和装置 |
CN107359616A (zh) * | 2017-07-18 | 2017-11-17 | 国电南瑞科技股份有限公司 | 一种解决大规模电网方式调整后潮流计算不收敛的方法 |
-
2018
- 2018-01-04 CN CN201810009362.5A patent/CN108281967B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130226482A1 (en) * | 2012-02-24 | 2013-08-29 | Hongbo Sun | Decoupled Three-Phase Power Flow Analysis Method for Unbalanced Power Distribution Systems |
CN106356860A (zh) * | 2016-10-17 | 2017-01-25 | 大连海事大学 | 一种配电系统三相潮流计算的电压初值设置方法 |
CN106374497A (zh) * | 2016-10-20 | 2017-02-01 | 天津大学 | 基于ami的中低压混合配电网一体化状态估计方法 |
CN107016489A (zh) * | 2017-03-09 | 2017-08-04 | 中国电力科学研究院 | 一种电力系统抗差状态估计方法和装置 |
CN107359616A (zh) * | 2017-07-18 | 2017-11-17 | 国电南瑞科技股份有限公司 | 一种解决大规模电网方式调整后潮流计算不收敛的方法 |
Non-Patent Citations (4)
Title |
---|
FENG XU, HAIHONG YU, YI LU: "A Shunt Type Power Flow Controller for Meshed DC Grids", 《12TH IET INTERNATIONAL CONFERENCE ON AC AND DC POWER TRANSMISSION (ACDC 2016)》 * |
李智欢: "几种提高牛顿法潮流收敛性的初值给定方法研究", 《电工电气》 * |
范翔: "提高电力系统潮流计算收敛性方法的比较分析", 《电力系统及其自动化学报》 * |
赵晋泉: "潮流计算中PV-PQ 节点转换逻辑的研究", 《中国电机工程学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112260283A (zh) * | 2020-10-22 | 2021-01-22 | 云南电网有限责任公司 | 一种电力系统输电断面潮流调整方法及装置 |
CN112260283B (zh) * | 2020-10-22 | 2023-08-22 | 云南电网有限责任公司 | 一种电力系统输电断面潮流调整方法及装置 |
CN117828243A (zh) * | 2024-03-06 | 2024-04-05 | 国网上海能源互联网研究院有限公司 | 一种fpga潮流并行计算系统和方法 |
CN117828243B (zh) * | 2024-03-06 | 2024-05-14 | 国网上海能源互联网研究院有限公司 | 一种fpga潮流并行计算系统和方法 |
Also Published As
Publication number | Publication date |
---|---|
CN108281967B (zh) | 2020-06-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107171341B (zh) | 一种基于分布式计算的输配电网一体化无功优化方法 | |
WO2018103317A1 (zh) | 一种含upfc电力系统的通用潮流计算方法 | |
CN108683191B (zh) | 一种下垂控制型孤岛微电网的三相潮流分析方法 | |
CN110797874B (zh) | 含电力电子变压器的交直流混合配电网状态估计方法 | |
CN102545252B (zh) | 基于三阶收敛牛顿法vsc-hvdc潮流计算方法 | |
CN109638839B (zh) | 一种双极柔性直流输电系统潮流计算方法 | |
CN108429294B (zh) | 一种含能量路由器的交直流网络潮流模型及求解方法 | |
CN110086200A (zh) | 一种孤岛模式下混合串并联微电网的协调控制方法 | |
CN106159955B (zh) | 基于连续惩罚对偶分解的电力系统分布式最优潮流方法 | |
CN103490428A (zh) | 微电网无功补偿容量配置方法及系统 | |
CN108808681A (zh) | 基于混合注入模型的光伏并网潮流计算方法 | |
CN107196321A (zh) | 一种提高电力弹簧稳态运行范围的方法 | |
CN106058898A (zh) | 一种配电网三相不对称直流潮流计算方法 | |
Miao et al. | Generalized steady-state model for energy router with applications in power flow calculation | |
CN108281967A (zh) | 一种改善电力系统潮流计算收敛性的算法启动方法 | |
CN105896558B (zh) | 一种基于vsc的upfc机电暂态模块化建模方法 | |
Chitsazan et al. | State estimation for large-scale power systems and facts devices based on spanning tree maximum exponential absolute value | |
Yang et al. | Implementation of a novel unified power flow controller into Newton-Raphson load flow | |
Tong et al. | A novel method of power flow analysis with UPFC considering limit violations of variable | |
CN107508318B (zh) | 一种基于电压灵敏度分区的有功控制方法及系统 | |
CN116304487A (zh) | 一种电力系统潮流反函数计算方法 | |
CN113270868B (zh) | 用于交流电气化铁路列车供电系统的动态潮流计算方法 | |
CN106877337B (zh) | 一种基于线性负荷模型的配电网单相潮流计算方法 | |
Jianming et al. | An improved power flow algorithm using equation changing method for AC/DC power system with VSC-HVDC | |
CN107800137A (zh) | 一种基于mmc的upfc机电暂态仿真方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |