CN103198205B - 适应于大规模数据潮流计算的初值计算方法 - Google Patents
适应于大规模数据潮流计算的初值计算方法 Download PDFInfo
- Publication number
- CN103198205B CN103198205B CN201310057906.2A CN201310057906A CN103198205B CN 103198205 B CN103198205 B CN 103198205B CN 201310057906 A CN201310057906 A CN 201310057906A CN 103198205 B CN103198205 B CN 103198205B
- Authority
- CN
- China
- Prior art keywords
- bus
- power
- default
- load flow
- equilibrium point
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 52
- 238000000034 method Methods 0.000 claims abstract description 40
- 238000000205 computational method Methods 0.000 claims description 3
- 238000009394 selective breeding Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 2
- 238000004088 simulation Methods 0.000 abstract description 2
- 238000004422 calculation algorithm Methods 0.000 description 11
- 230000006872 improvement Effects 0.000 description 7
- 238000000354 decomposition reaction Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000005303 weighing Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 235000019606 astringent taste Nutrition 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提供了一种适应于大规模数据潮流计算的初值计算方法，包括如下步骤：(1).预设平衡节点；(2).进行潮流计算仿真；(3).判断计算结果是否收敛，如果收敛，则进行步骤4，否则结束；(4).将预设的平衡节点改为原始节点类型；(5).进行潮流计算。本发明提供的适应于大规模数据潮流计算的初值计算方法，对大规模电网出现的潮流计算问题的原因进行了分析，在现有算法基础上按照实际需求情况对牛顿法进行改进，以解决大规模电网数据潮流计算不收敛问题。
Description
技术领域
本发明属于电力系统仿真分析领域，具体涉及一种适应于大规模数据潮流计算的初值计算方法。
背景技术
潮流计算是电力系统各种静态和动态分析计算的基础。成熟稳定、收敛速度快、收敛性能好的电力系统潮流算法，对电网规划部门和电网调度运行部门进行电网规划和年度运行方式分析至关重要。
潮流计算的本质是求解大规模非线性方程组，牛顿-拉夫逊法(简称牛顿法)是求解非线性方程组的经典方法。上个世纪60年代后期，随着稀疏技术和节点优化编号技术的引入，因其收敛性好、计算速度快、内存需求少的优点，牛顿拉夫逊法成为了电力系统潮流计算中一种的优秀算法，目前是潮流计算的主流算法。其后，PQ分解法根据电力系统的特点，对牛顿法进行了PQ解耦，进一步减少了计算时间，也得到了迅速的推广。
在以往的实际电网分析中，传统的潮流算法均能较好的收敛，且潮流结果正常合理。但近年来随着电网规模的不断扩大，特别是随着特高压电网的建设，传统的潮流算法面临收敛性能逐步恶化的趋势。目前电网已形成了华北和华中电网通过特高压同步联网，计算规模已经超过1万母线。母线数是衡量计算规模的重要指标，母线数超过1万，潮流不收敛的情况明显增加。按照国家电网公司发展规划，还将实现华北-华中电网与华东电网的联网，电网规模还会进一步扩大。因此，迫切需要改进现有的潮流算法，以满足大规模数据潮流计算的现实需求。
牛顿法是一种迭代算法，对初值设定很敏感，可以用图1来说明：
设潮流计算求解的非线性方程组为F(x)＝0，如图1所示，x0为初值，x为最终的真值。当初值x0取得合适时，通过x0→x1→x2...逼近真值；当初值为x0’时，由x0’→x1’→x2’...离开真值越来越远。一般来说，如果有良好的初值，牛顿法的求解是相当可靠的。
为了提高牛顿法的收敛性，现有技术常常采用PQ分解法、直流法、高斯塞德尔法等算法先迭代几次，取得潮流初值后再转为牛顿法的初值改善策略。PQ分解法等算法本身有一定限制，例如PQ分解法要求支路×＞＞R，直流法假定电压为1.0p.u.，高斯塞德尔法对PV节点处理较难且收敛性差，这些条件在大规模电网数据中往往不能满足，实际应用效果并不好。
发明内容
为克服上述缺陷，本发明提供了一种适应于大规模数据潮流计算的初值计算方法，对大规模电网出现的潮流计算问题的原因进行了分析，在现有算法基础上按照实际需求情况对牛顿法进行改进，以解决大规模电网数据潮流计算不收敛问题。
为实现上述目的，本发明提供一种适应于大规模数据潮流计算的初值计算方法，其改进之处在于，所述方法包括如下步骤：
(1).预设平衡节点；
(2).进行潮流计算仿真；
(3).判断计算结果是否收敛，如果收敛，则进行步骤4，否则结束；
(4).将预设的平衡节点改为原始节点类型；
(5).进行潮流计算。
本发明提供的优选技术方案中，在所述步骤1中，对平衡节点进行预设的方式包括人工设定和自动设定。
本发明提供的第二优选技术方案中，所述步骤1包括如下步骤：
(1-1).选择预设平衡节点的方式，如果是人工设定则进行步骤1-2，否则进行步骤1-3；
(1-2).人工选择预设平衡点，结束；
(1-3).计算功率最大误差母线并设置预设平衡点。
本发明提供的第三优选技术方案中，在所述步骤1-2中，由工程计算人员指定预设平衡点，在远离已设定平衡机的区域外指定一个大容量电厂作为预设平衡点；并且预设平衡点的数量为1-2个。
本发明提供的第四优选技术方案中，在所述步骤1-3中，功率缺额较大的点是容易导致不收敛的关键节点，将预设平衡点设置在功率误差最大的母线上，能最大程度地提供系统支撑；功率最大误差的计算方法如下：
ΔSmax＝Max(P1，P2，...，Pn，，Q1，Q2，...，Qn)
其中，P1到Pn为第1个母线到第n个母线的节点注入有功功率，n为母线总数；Q1到Qn为第1个母线到第n个母线的节点注入无功功率；ΔSmax为功率最大误差，功率最大误差所在的第i个母线，即满足条件Max(Pi，Qi)＝ΔSmax的节点，自动设置为预设平衡点。
与现有技术比，本发明提供的一种适应于大规模数据潮流计算的初值计算方法，从牛顿法初值改善策略方面提出了读预设平衡点策略，与PQ分解法、直流法、高斯塞德尔法等现有的初值改善策略相比，没有额外的前提条件，完全从实际电网潮流计算的数据特点出发，更能适应实际大规模电网潮流计算，有效提高了原有潮流算法的收敛性；再者，通过实践特高压联网数据的潮流计算也验证了所提方法的正确性和适应性。
附图说明
图1为现有的牛顿迭代算法的示意图。
图2为适应于大规模数据潮流计算的初值计算方法的流程图。
具体实施方式
如图2所示，一种适应于大规模数据潮流计算的初值计算方法，包括如下步骤：
(1).预设平衡节点；
(2).进行潮流计算仿真；
(3).判断计算结果是否收敛，如果收敛，则进行步骤4，否则结束；
(4).将预设的平衡节点改为原始节点类型；
(5).进行潮流计算。
在所述步骤1中，对平衡节点进行预设的方式包括人工设定和自动设定。
所述步骤1包括如下步骤：
(1-1).选择预设平衡节点的方式，如果是人工设定则进行步骤1-2，否则进行步骤1-3；
(1-2).人工选择预设平衡点，结束；
(1-3).计算功率最大误差母线并设置预设平衡点。
在所述步骤1-2中，由工程计算人员指定预设平衡点，在远离已设定平衡机的区域外指定一个大容量电厂作为预设平衡点；并且预设平衡点的数量为1-2个。
本发明提供的第四优选技术方案中，在所述步骤1-3中，功率缺额较大的点是容易导致不收敛的关键节点，将预设平衡点设置在功率误差最大的母线上，能最大程度地提供系统支撑；功率最大误差的计算方法如下：
ΔSmax＝Max(P1，P2，...，Pn，，Q1，Q2，...，Qn)
其中，P1到Pn为第1个母线到第n个母线的节点注入有功功率，n为母线总数；Q1到Qn为第1个母线到第n个母线的节点注入无功功率；ΔSmax为功率最大误差，功率最大误差所在的第i个母线，即满足条件Max(Pi，Qi)＝ΔSmax的节点，自动设置为预设平衡点。
通过以下实施例对适应于大规模数据潮流计算的初值计算方法做进一步描述。
1、预设平衡点的初值改善策略
平衡节点的选取对牛顿法潮流计算的收敛性影响很大。平衡节点(slack bus)是节点电压幅值和相角是给定的节点，它对系统起到功率平衡的作用，一般来说在一个交流互联电网中平衡节点只有一个。在实际大规模电网潮流计算中发现，在关键地方多设平衡点作为系统支撑，可以显著提高潮流收敛性。
因此，将潮流计算分为两个阶段：初值计算和正常计算。初值计算阶段，多设平衡节点进行牛顿法的迭代计算，得到潮流初值；正常计算阶段，将刚才设置的平衡节点还原回原始节点类型，以初值计算阶段得到的潮流初值来启动正常的迭代计算。计算过程中设置的“平衡节点”并不是真正的平衡节点，只是为了得到潮流初值辅助收敛的，可命名为“预设平衡点”，这种初值改善策略称为预设平衡点初值改善策略。
2、预设平衡点位置选择
预设平衡点设置的不合理，可能在潮流初值计算阶段就不收敛，该功能就无法起到作用，所以需合理设置预设平衡点的位置。设置预设平衡点的方法有两种：
1)人工设定。即由工程计算人员指定预设平衡点的位置，一般来说应在远离已设定平衡机的区域外指定一个大容量电厂作为预设平衡点。并且预设平衡点不宜设置过多，1-2个比较合适。
2)自动设定。根据分析，功率缺额较大的点是容易导致不收敛的关键节点，将预设平衡点设置在功率误差最大的母线上，能最大程度地提供系统支撑。功率最大误差的计算方法如下：
ΔSmax＝Max(P1，P2，...，Pn，，Q1，Q2，...，Qn)
其中P1到Pn为第1个母线到第n个母线的节点注入有功功率，n为母线总数；Q1到Qn为第1个母线到第n个母线的节点注入无功功率。ΔSmax为功率最大误差，功率最大误差所在的第i个母线(满足条件Max(Pi，Qi)＝ΔSmax)自动设置为预设平衡点。
需要声明的是，本发明内容及具体实施方式意在证明本发明所提供技术方案的实际应用，不应解释为对本发明保护范围的限定。本领域技术人员在本发明的精神和原理启发下，可作各种修改、等同替换、或改进。但这些变更或修改均在申请待批的保护范围内。
Claims (2)
1.一种适应于大规模数据潮流计算的初值计算方法，其特征在于，所述方法包括如下步骤：
(1).预设平衡节点；
(2).进行潮流计算仿真；
(3).判断计算结果是否收敛，如果收敛，则进行步骤4，否则结束；
(4).将预设的平衡节点改为原始节点类型；
(5).进行潮流计算；
在所述步骤1中，对平衡节点进行预设的方式包括人工设定和自动设定；
所述步骤1包括如下步骤：
(1-1).选择预设平衡节点的方式，如果是人工设定则进行步骤1-2，否则进行步骤1-3；
(1-2).人工选择预设平衡点，结束；
(1-3).计算功率最大误差母线并设置预设平衡点，功率缺额较大的点是容易导致不收敛的关键节点，将预设平衡点设置在功率误差最大的母线上，能最大程度地提供系统支撑；功率最大误差的计算方法如下：
ΔSmax＝Max(P1,P2,…,Pn,Q1,Q2,…,Qn)
其中，P1到Pn为第1个母线到第n个母线的节点注入有功功率，n为母线总数；Q1到Qn为第1个母线到第n个母线的节点注入无功功率；ΔSmax为功率最大误差，功率最大误差所在的第i个母线，即满足条件Max(Pi,Qi)＝ΔSmax的节点，自动设置为预设平衡点。
2.根据权利要求1所述的方法，其特征在于，在所述步骤1-2中，由工程计算人员指定预设平衡点，在远离已设定平衡机的区域外指定一个大容量电厂 作为预设平衡点；并且预设平衡点的数量为1-2个。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310057906.2A CN103198205B (zh) | 2013-02-25 | 2013-02-25 | 适应于大规模数据潮流计算的初值计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310057906.2A CN103198205B (zh) | 2013-02-25 | 2013-02-25 | 适应于大规模数据潮流计算的初值计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103198205A CN103198205A (zh) | 2013-07-10 |
| CN103198205B true CN103198205B (zh) | 2017-02-08 |
Family
ID=48720760
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310057906.2A Active CN103198205B (zh) | 2013-02-25 | 2013-02-25 | 适应于大规模数据潮流计算的初值计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103198205B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108899909B9 (zh) * | 2018-08-14 | 2020-04-14 | 清华大学 | 基于低非线性度网络模型的用于潮流计算的初值估计方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101141064A (zh) * | 2007-09-14 | 2008-03-12 | 清华大学 | 通过交换边界节点状态和网损信息完成分布式潮流分析的方法 |
| CA2661753C (en) * | 2006-08-26 | 2011-10-11 | Sureshchandra B. Patel | Method of loadflow computation for electrical power system |
| CN102708212A (zh) * | 2012-06-01 | 2012-10-03 | 广西大学 | 一种bpa双向数据接口方法 |
-
2013
- 2013-02-25 CN CN201310057906.2A patent/CN103198205B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CA2661753C (en) * | 2006-08-26 | 2011-10-11 | Sureshchandra B. Patel | Method of loadflow computation for electrical power system |
| CN101141064A (zh) * | 2007-09-14 | 2008-03-12 | 清华大学 | 通过交换边界节点状态和网损信息完成分布式潮流分析的方法 |
| CN102708212A (zh) * | 2012-06-01 | 2012-10-03 | 广西大学 | 一种bpa双向数据接口方法 |
Non-Patent Citations (1)
| Title |
|---|
| 电力系统潮流收敛性的实用性改进;于尔铿;《电网技术》;19950131;第19卷(第1期);23-26 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103198205A (zh) | 2013-07-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106451556B (zh) | 一种确定配电网馈线接纳分布式电源能力的方法和装置 | |
| CN104037793B (zh) | 一种应用于主动配电网的储能单元容量配置方法 | |
| CN104281737B (zh) | 一种海上风电场集电系统拓扑结构选型方法 | |
| CN103607466B (zh) | 一种基于云计算的广域多级分布式并行电网分析方法 | |
| CN103001214B (zh) | 一种基于中性点偏移的配电网三相潮流计算方法 | |
| CN111797510A (zh) | 一种新能源场站短路比的计算方法及系统 | |
| CN107294120A (zh) | 一种主动配电网混合储能容量优化配置方法及装置 | |
| CN103106314B (zh) | 太阳能光伏电源输出功率的时序概率建模方法 | |
| CN102856899B (zh) | 一种微电网降低网损的方法 | |
| CN103441527B (zh) | 一种基于实测数据的风电接入系统模型 | |
| CN112862317B (zh) | 一种协调5g通讯基站选址接入的配电网网架规划方法 | |
| CN106033888B (zh) | 基于速降直流灵敏度的交流线路故障安控量在线评估方法 | |
| CN112260326B (zh) | 一种新能源集群的等效短路比计算方法及系统 | |
| CN106684917A (zh) | 分布式电源的位置动态优化及准入容量计算方法及系统 | |
| CN115622053A (zh) | 一种用于考虑分布式电源的自动负荷建模方法及装置 | |
| CN112072692A (zh) | 一种新能源发电场站的阻抗等值方法及装置 | |
| CN103186366A (zh) | 基于cuda并行计算实现电力系统电磁暂态实时仿真测试方法 | |
| CN103198205B (zh) | 适应于大规模数据潮流计算的初值计算方法 | |
| CN105162119A (zh) | 一种适用于大规模配电网的局部修正潮流计算方法 | |
| CN109560568A (zh) | 基于短路电流裕度的双馈风电场最大可接入容量确定方法 | |
| de Lima Vianna et al. | A general analysis of the distributed generation impact on electrical energy planning | |
| CN109066709A (zh) | 一种基于元模型的分布式电源就地电压控制策略改进方法 | |
| CN113468729B (zh) | 无标定方式的电力厂站计量装置运行误差测算方法及系统 | |
| CN105046583A (zh) | 一种适用于分布式实时数据处理的电网模型分区方法 | |
| Xia et al. | Assessment of the carrying capacity of distributed generation in distribution network considering soft open point |
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 |