CN106410789B - 一种针对孤网的电力系统黑启动分区方法 - Google Patents
一种针对孤网的电力系统黑启动分区方法 Download PDFInfo
- Publication number
- CN106410789B CN106410789B CN201610896601.4A CN201610896601A CN106410789B CN 106410789 B CN106410789 B CN 106410789B CN 201610896601 A CN201610896601 A CN 201610896601A CN 106410789 B CN106410789 B CN 106410789B
- Authority
- CN
- China
- Prior art keywords
- unit
- subregion
- launched
- black starting
- starting
- 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.)
- Expired - Fee Related
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
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种针对孤网的电力系统黑启动分区方法,包括以下步骤:S1:确定分区个数:确定孤网中黑启动机组、待启动机组和负荷这三类节点的数量,再根据黑启动机组和待启动机组的数量关系确定分区个数;S2:机组分区:结合孤网实际拓扑结构,对黑启动机组和待启动机组按照以下三个原则进行划分:(1)以最短时间恢复电力系统中的重要输电线路;(2)以最短时间恢复电网的发电能力;(3)以最短时间建立恢复网架;S3:负荷分区:计算每个待启动机组分区的有功和无功容量,对各个待启动机组分区划入负荷节点,以平衡待启动机组分区的功率。本发明实现了孤网拓扑的简化,构建出最佳恢复网架,增强了孤网恢复的快速性和稳定性。
Description
技术领域
本发明涉及电力系统故障恢复,特别是涉及一种针对孤网的电力系统黑启动分区方法。
背景技术
随着经济的快速发展,现代电力系统日趋复杂化,导致电网抵抗大扰动能力却在不断被弱化。大规模电力系统中,当局部电网发生严重故障时,一旦调度人员操作不当或保护装置不正确动作,将会引发电力系统的连锁反应,扩大故障影响,导致大范围停电甚至是全网的崩溃。如果能在大停电事故时,采取合理快捷的黑启动措施,及时排除故障,减小对经济、社会、民生的不良影响,具有非常重要的意义。但是常规黑启动方案只研究单个黑启动机组的恢复方法,无法同时启动多台机组,浪费了启动时间。除此以外,近年来不断引入智能算法求解电力系统的恢复网架以及分区方案,存在计算量较大,结果稳定性较差等问题,难以应用于工程实践。为了能实现复数台机组同时启动,降低分区方法的计算量并提高结果稳定性,需要研究一种新的启动方法。
发明内容
发明目的:本发明的目的是提供一种能够解决现有技术中存在的缺陷的一种针对孤网的电力系统黑启动分区方法。
技术方案:为达到此目的,本发明采用以下技术方案:
本发明所述的针对孤网的电力系统黑启动分区方法,其特征在于:包括以下步骤:
S1:确定分区个数:确定孤网中黑启动机组、待启动机组和负荷这三类节点的数量,再根据黑启动机组和待启动机组的数量关系确定分区个数;
S2:机组分区:结合孤网实际拓扑结构,对黑启动机组和待启动机组按照以下三个原则进行划分:(1)以最短时间恢复电力系统中的重要输电线路;(2)以最短时间恢复电网的发电能力;(3)以最短时间建立恢复网架;
S3:负荷分区:计算每个分区的有功和无功容量,对各个分区划入负荷节点,以平衡单个分区的功率。
2、根据权利要求1所述的针对孤网的电力系统黑启动分区方法,其特征在于:所述步骤S1中的分区个数按照如下的方法来确定:如果黑启动机组个数小于待启动机组个数,则分区数等于黑启动机组个数;如果黑启动机组个数大于待启动机组个数,则分区数等于待启动机组个数。
3、根据权利要求1所述的针对孤网的电力系统黑启动分区方法,其特征在于:所述步骤S2包括以下步骤:
S2.1:定义分区目标函数:
s.t.s1g+s2g+....+sng=1 (1)
式(1)中,m为待启动机组数量;n为黑启动机组数量;[s1g,s2g,....,sng]T表示g号待启动机组是否属于黑启动机组的布尔变量:sig=0表明g号待启动机组不属于i号黑启动机组所在的分区,sig=1表明g号待启动机组属于i号黑启动机组所在的分区,1≤i≤n;[b1g,b2g,....,bng]T中的big表示g号待启动机组到i号黑启动机组所需的最短启动时间,1≤i≤n;
S2.2:定义约束条件:
S2.3:待启动机组分区:
按照式(3)对sig进行求解:
4、根据权利要求1所述的针对孤网的电力系统黑启动分区方法,其特征在于:所述步骤S2.1中big按照如下的步骤计算得到:
S2.11:初始化:选取g号待启动机组为起始点s,集合P为空,并标记起始点s,记k=s,孤网中剩余节点设为未标记点;
S2.12:寻找一个不在P中的节点j,检验从所有已标记的点k到与其直接连接的未标记的点j的距离,并设置
dj=min[dj,dk+lkj] (4)
式(4)中,lkj为从点k到j的直接连接距离,dj为起始点s到点j的距离;
S2.13:选取下一个点:从所有未标记的点中,选取到起始点s距离最小的点,并将该点进行标记;
S2.14:检查是否所有的点都已经标记:如果都标记,则最后一个未标记节点到起始点s的距离即为big;否则,返回步骤S2.12。
5、根据权利要求1所述的针对孤网的电力系统黑启动分区方法,其特征在于:所述步骤S3包括以下步骤:
S3.1:定义分区目标函数:
s.t.c1k+c2k+....+cnk=1
(5)
式(5)中,l为负荷节点数量;n为黑启动机组数量;Yk=[y1k,y2k,....,ynk]T表示k号负荷节点是否属于黑启动机组的布尔变量:yik=0表明k号负荷节点不属于i号黑启动机组所在的分区,yik=1表明k号负荷节点属于i号黑启动机组所在的分区,1≤i≤n;[c1k,c2k,....,cnk]T中的cik表示k号负荷到i号黑启动机组所需的最短启动时间,1≤i≤n;
S3.2:定义约束条件:
有功功率平衡约束如式(6)所示:
式(6)中,PGa为a号黑启动机组所在的分区中机组的有功容量之和,pk为k号负荷节点消耗的有功功率,α1和β1为有功功率平衡约束系数,yak表示k号负荷节点与a号黑启动机组的从属关系;
无功功率平衡约束如式(6)所示:
式(7)中,QGa为a号黑启动机组所在的分区中机组的无功容量之和,qk为k号负荷节点消耗的无功功率,α2和β2为无功功率平衡约束系数;
S3.3:负荷分区:
按照式(8)对yik进行求解:
有益效果:本发明针对含有复数台黑启动机组的孤网分区问题进行了较为细致的研究,并建立了一个完整的分区模型。本发明依次通过确定分区数量、机组分区和负荷分区这三个步骤,将复杂的多目标非线性规划问题化简为单目标多约束的线性规划问题。本发明实现了孤网拓扑的简化,构建出最佳恢复网架,增强了孤网恢复的快速性和稳定性,为黑启动路径的分析工作奠定了良好的基础。
附图说明
图1为本发明具体实施方式的总体流程图;
图2为本发明具体实施方式中最短启动时间算法的流程图;
图3为本发明具体实施方式中的新英格兰10机39节点系统网络的分区结果图。
具体实施方式
本具体实施方式公开了一种针对孤网的电力系统黑启动分区方法,如图1所示,包括以下步骤:
S1:确定分区个数:确定孤网中黑启动机组的数量n、待启动机组的数量m和负荷的数量l,再根据黑启动机组和待启动机组的数量关系确定分区个数;
S2:机组分区:结合孤网实际拓扑结构,对黑启动机组和待启动机组按照以下三个原则进行划分:(1)以最短时间恢复电力系统中的重要输电线路;(2)以最短时间恢复电网的发电能力;(3)以最短时间建立恢复网架;因此,以机组恢复时间最短为目标函数,黑启动机组平衡为约束构造单目标多约束的目标函数,对分区结果进行求解;
S3:负荷分区:计算每个待启动机组分区的有功和无功容量,对各个待启动机组分区划入负荷节点,以平衡待启动机组分区的功率。因此,以负荷恢复时间最短为目标函数,分区有功、无功平衡为约束构造单目标多约束的目标函数,对分区结果进行求解。
下面以新英格兰10机39节点系统为例进行计算,来证明本方法的可行性。
表1新英格兰10机39节点系统中线路的恢复时间
输电线路的恢复时间如表1所示,假设节点30、31和34为黑启动机组,32、33、35、36、37、38、39为待启动机组,剩下的节点即为负荷节点。
S1:确定分区个数:
比较黑启动机组和待启动机组的个数,取两者中的最小值3为分区个数。
S2:机组分区:
表2新英格兰10机39节点系统中待恢复机组的恢复矩阵
以节点30、31、34分别为分区1,分区2,分区3的源节点。针对7台待启动机组构造7个恢复时间矩阵。分别以节点32、33、35、36、37、38、39为末节点进行最短启动时间计算,最短启动时间的算法流程图如图2所示,求解的数值放入7个3维恢复时间矩阵B1、B2、…、B7中。对于系统中任意一个待启动节点,如节点37,其从3个分区的黑启动机组获取启动功率的时间分别为6min、20min、19min,这样节点37所对应的恢复时间矩阵为[6,20,19]T,待启动机组的恢复时间向量如表2所示。
根据求解得到的恢复时间矩阵,可以得到分区机组矩阵为:
由求解得到的分区机组矩阵可得,节点37、38、39与黑启动节点30同属于分区1;节点33、35、36与黑启动节点34同属于分区2;节点32与黑启动节点31同属于分区3。
S3:负荷分区:
在得到待启动机组分区后,可以继续进行负荷分区的计算。负荷分区计算与待启动机组分区相类似,取约束系数α1、β1为0.8,α2、β2为1.2,分别以3个黑启动节点为源点,29个负荷节点为末端进行最短启动时间计算,构造29个3维恢复时间矩阵,得到的分区结果如表3和图3所示。
表3新英格兰10机39节点负荷节点分区结果
为证明本发明的有效性,将本发明与现有文献所求得的分区恢复网架进行比分析。现有文献将39号待启动机组划入黑启动机组31所在的分区,启动所需时间为15min,而本发明将39号机组划划入黑启动机组30所在的分区,启动所需时间为10min。除此以外,文献所求得的分区1中负荷消耗有功功率超过机组发出有功17%,分区2机组发出的无功功率较负荷高出13%;本发明采用了相关的约束条件进行计算,功率不平衡量可以保证在5%以内。因此,采用本发明的分区方法和恢复网架更有利于孤网黑启动。详细的对比参数如表4所示。
表4分区参数对比
Claims (4)
1.一种针对孤网的电力系统黑启动分区方法,其特征在于:包括以下步骤:
S1:确定分区个数:确定孤网中黑启动机组、待启动机组和负荷这三类节点的数量,再根据黑启动机组和待启动机组的数量关系确定分区个数;
S2:机组分区:结合孤网实际拓扑结构,对黑启动机组和待启动机组按照以下三个原则进行划分:(1)以最短时间恢复电力系统中的重要输电线路;(2)以最短时间恢复电网的发电能力;(3)以最短时间建立恢复网架;
S3:负荷分区:计算每个分区的有功和无功容量,对各个分区划入负荷节点,以平衡单个分区的功率;
所述步骤S3包括以下步骤:
S3.1:定义分区目标函数:
s.t. y1k+y2k+....+ynk=1
(5)
式(5)中,l为负荷节点数量;n为黑启动机组数量;Yk=[y1k,y2k,....,ynk]T表示k号负荷节点是否属于黑启动机组的布尔变量:yik=0表明k号负荷节点不属于i号黑启动机组所在的分区,yik=1表明k号负荷节点属于i号黑启动机组所在的分区,1≤i≤n;[c1k,c2k,....,cnk]T中的cik表示k号负荷到i号黑启动机组所需的最短启动时间,1≤i≤n;
S3.2:定义约束条件:
有功功率平衡约束如式(6)所示:
式(6)中,PGa为a号黑启动机组所在的分区中机组的有功容量之和,pk为k号负荷节点消耗的有功功率,α1和β1为有功功率平衡约束系数,yak表示k号负荷节点与a号黑启动机组的从属关系;
无功功率平衡约束如式(7)所示:
式(7)中,QGa为a号黑启动机组所在的分区中机组的无功容量之和,qk为k号负荷节点消耗的无功功率,α2和β2为无功功率平衡约束系数;
S3.3:负荷分区:
按照式(8)对yik进行求解:
2.根据权利要求1所述的针对孤网的电力系统黑启动分区方法,其特征在于:所述步骤S1中的分区个数按照如下的方法来确定:如果黑启动机组个数小于待启动机组个数,则分区数等于黑启动机组个数;如果黑启动机组个数大于待启动机组个数,则分区数等于待启动机组个数。
3.根据权利要求1所述的针对孤网的电力系统黑启动分区方法,其特征在于:所述步骤S2包括以下步骤:
S2.1:定义分区目标函数:
s.t. s1g+s2g+....+sng=1 (1)
式(1)中,m为待启动机组数量;n为黑启动机组数量;[s1g,s2g,....,sng]T表示g号待启动机组是否属于黑启动机组的布尔变量:sig=0表明g号待启动机组不属于i号黑启动机组所在的分区,sig=1表明g号待启动机组属于i号黑启动机组所在的分区,1≤i≤n;[b1g,b2g,....,bng]T中的big表示g号待启动机组到i号黑启动机组所需的最短启动时间,1≤i≤n;
S2.2:定义约束条件:
S2.3:待启动机组分区:
按照式(3)对sig进行求解:
4.根据权利要求3所述的针对孤网的电力系统黑启动分区方法,其特征在于:所述步骤S2.1中big按照如下的步骤计算得到:
S2.11:初始化:选取g号待启动机组为起始点s,集合P为空,并标记起始点s,记k=s,孤网中剩余节点设为未标记点;
S2.12:寻找一个不在P中的节点j,检验从所有已标记的点k到与其直接连接的未标记的点j的距离,并设置
dj=min[dj,dk+lkj] (4)
式(4)中,lkj为从点k到j的直接连接距离,dj为起始点s到点j的距离;
S2.13:选取下一个点:从所有未标记的点中,选取到起始点s距离最小的点,并将该点进行标记;
S2.14:检查是否所有的点都已经标记:如果都标记,则最后一个未标记节点到起始点s的距离即为big;否则,返回步骤S2.12。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610896601.4A CN106410789B (zh) | 2016-10-13 | 2016-10-13 | 一种针对孤网的电力系统黑启动分区方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610896601.4A CN106410789B (zh) | 2016-10-13 | 2016-10-13 | 一种针对孤网的电力系统黑启动分区方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106410789A CN106410789A (zh) | 2017-02-15 |
CN106410789B true CN106410789B (zh) | 2018-12-28 |
Family
ID=59228568
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610896601.4A Expired - Fee Related CN106410789B (zh) | 2016-10-13 | 2016-10-13 | 一种针对孤网的电力系统黑启动分区方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106410789B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107368930B (zh) * | 2017-08-07 | 2020-11-03 | 国网(天津)综合能源服务有限公司 | 一种适用于电力系统分区恢复模式的黑启动电源选址方法 |
CN110676869B (zh) * | 2018-07-02 | 2022-10-11 | 国网山东省电力公司潍坊供电公司 | 一种确定新增黑启动机组布点的方法 |
CN109713669B (zh) * | 2019-01-28 | 2022-01-21 | 广东电网有限责任公司 | 一种抗灾保底电网的规划方法、装置和设备 |
CN111313461B (zh) * | 2020-01-20 | 2022-07-05 | 国网福建省电力有限公司 | 一种基于分步求解的黑启动分区方法 |
CN112288237A (zh) * | 2020-10-12 | 2021-01-29 | 国网江苏省电力有限公司镇江供电分公司 | 一种基于协同度优化的黑启动子区域划分方法 |
CN112884608B (zh) * | 2021-03-11 | 2023-07-25 | 国网陕西省电力公司电力科学研究院 | 一种供电网黑启动分区方案计算方法及计算系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103579990A (zh) * | 2013-10-17 | 2014-02-12 | 华南理工大学 | 基于混合决策图的含分布式电源配电系统故障恢复方法 |
CN103792924A (zh) * | 2014-02-14 | 2014-05-14 | 暨南大学 | 一种含有微网的电力系统扩展黑启动方法 |
CN105470986A (zh) * | 2015-12-17 | 2016-04-06 | 国家电网公司 | 一种电力系统分区方法 |
-
2016
- 2016-10-13 CN CN201610896601.4A patent/CN106410789B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103579990A (zh) * | 2013-10-17 | 2014-02-12 | 华南理工大学 | 基于混合决策图的含分布式电源配电系统故障恢复方法 |
CN103792924A (zh) * | 2014-02-14 | 2014-05-14 | 暨南大学 | 一种含有微网的电力系统扩展黑启动方法 |
CN105470986A (zh) * | 2015-12-17 | 2016-04-06 | 国家电网公司 | 一种电力系统分区方法 |
Non-Patent Citations (2)
Title |
---|
A Sectionalizing Method in Power System Restoration Based on WAMS;S. Arash Nezam Sarmadi 等;《IEEE TRANSACTIONS ON SMART GRID》;20110331;第2卷(第1期);190-197 * |
大停电后电力系统黑启动分区的两步策略;孙磊 等;《电力自动化设备》;20150930;第35卷(第9期);14-21 * |
Also Published As
Publication number | Publication date |
---|---|
CN106410789A (zh) | 2017-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106410789B (zh) | 一种针对孤网的电力系统黑启动分区方法 | |
CN103296677B (zh) | 一种在线大电网恢复辅助决策系统 | |
CN101719182B (zh) | 一种交直流电力系统分割并行电磁暂态数字仿真方法 | |
Zhu et al. | Adaptive power flow method for distribution systems with dispersed generation | |
CN105260952B (zh) | 基于马尔科夫链蒙特卡洛方法的光伏电站可靠性评估方法 | |
CN102799950B (zh) | 基于粒子群算法的船舶电网重构优化方法 | |
CN109217295B (zh) | 预防系统过载的潮流灵敏度计算方法和计算机装置 | |
CN107292481B (zh) | 一种基于节点重要度的电网关键节点评估方法 | |
CN106058876A (zh) | 考虑暂态电压稳定性的动态无功规划选址分析方法及系统 | |
CN113300367B (zh) | 电网关键断面在线快速识别方法、输电限额制定方法及装置 | |
CN106972517B (zh) | 基于双极对称特征的特高压直流输电系统可靠性计算方法 | |
CN107622360A (zh) | 一种综合考虑主客观因素的关键线路识别方法 | |
CN108830451B (zh) | 一种用户侧分布式储能的汇聚潜力评估方法及系统 | |
CN103824234A (zh) | 一种基于分块和层级结构的配电系统可靠性评估方法 | |
CN104505821A (zh) | 一种控制短路电流水平的电网运行方式优化方法 | |
CN105552899B (zh) | 一种计算大停电后电网恢复力的方法 | |
CN105656036A (zh) | 考虑潮流和灵敏度一致性等值的概率静态安全分析方法 | |
CN104504524A (zh) | 一种应用于有源配网的负荷削减法及可靠性评估方法 | |
CN108964151B (zh) | 一种微电网网架结构设计方法及系统 | |
CN103296678B (zh) | 一种在线大电网恢复辅助决策方法 | |
CN106655167B (zh) | 电网连锁故障的仿真方法及装置 | |
CN104283208B (zh) | 大型电网概率性可用输电能力的分解协调计算方法 | |
CN104578057B (zh) | 一种考虑源流路径的含风电电网电压主导节点识别方法 | |
TWI459677B (zh) | Analysis Method of High Performance Micro - grid Isolated Operation Fault | |
CN106058856B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181228 Termination date: 20191013 |
|
CF01 | Termination of patent right due to non-payment of annual fee |