CN104332974B - 考虑电网调度员实际操作特性的连锁故障模拟方法 - Google Patents

考虑电网调度员实际操作特性的连锁故障模拟方法 Download PDF

Info

Publication number
CN104332974B
CN104332974B CN201410654244.1A CN201410654244A CN104332974B CN 104332974 B CN104332974 B CN 104332974B CN 201410654244 A CN201410654244 A CN 201410654244A CN 104332974 B CN104332974 B CN 104332974B
Authority
CN
China
Prior art keywords
cascading failure
node
load
branch road
branch
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
Application number
CN201410654244.1A
Other languages
English (en)
Other versions
CN104332974A (zh
Inventor
余晓鹏
黄少伟
赵阳
姚锐
高昆
梅生伟
李晓萌
刘锋
张雪敏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tsinghua University
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Henan Electric Power Co Ltd
Original Assignee
Tsinghua University
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Henan Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tsinghua University, State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Henan Electric Power Co Ltd filed Critical Tsinghua University
Priority to CN201410654244.1A priority Critical patent/CN104332974B/zh
Publication of CN104332974A publication Critical patent/CN104332974A/zh
Application granted granted Critical
Publication of CN104332974B publication Critical patent/CN104332974B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04SSYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00Systems supporting electrical power generation, transmission or distribution
    • Y04S10/20Systems supporting electrical power generation, transmission or distribution using protection elements, arrangements or systems

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种考虑电网调度员实际操作特性的连锁故障模型,包括以下步骤:(1)、设定连锁故障模拟天数上限,支路重载负载率阈值;设系统各个节点的交流潮流基准状态为:PQ、PV节点的有功注入、无功注入,和电压,其中为节点编号;(2)、初始化连锁故障天数;(3)、模拟第天连锁故障过程;(4)、;(5)、若,则转到步骤(3),否则结束。本发明通过在连锁故障模型中引入考虑实际电网调度员调整模拟,使连锁故障模拟体现了电网连锁故障过程与调度员的交互特征,更准确地反映连锁故障的发展特征,提高了连锁故障模型的实用性。

Description

考虑电网调度员实际操作特性的连锁故障模拟方法
技术领域
本发明涉及电力系统连锁故障模拟与风险评估领域,具体涉及一种考虑电网调度员实际操作特性的连锁故障模拟方法。
背景技术
电力系统是一个动态复杂的大规模非线性系统,随着电力系统规模日趋庞大,负荷需求不断增加,保障电网安全稳定运行的难度日趋增加。近年来,世界范围内发生多起由于自身事故、外部因素扰动等原因引发的电力系统连锁性故障,并导致大规模停电,造成了严重的社会经济损失。
为了对电力系统连锁故障风险进行评估,并研究连锁故障的发展机理,进行合理、高效的电力系统连锁故障模拟必不可少。目前,基于准稳态潮流的电力系统连锁故障模拟模型广泛应用于连锁故障研究与评估工作中。
连锁故障模型应包括连锁故障过程中的关键参与因素,并反映连锁故障发展的基本特征。电力系统作为在人和控制器操纵下的大型复杂系统,其调度操作是连锁故障发展的重要影响因素。由于连锁故障过程中,调度员会对故障状态下的系统进行控制操作,以降低系统的风险并保障最大程度的供电。传统的连锁故障模型往往将其建立为一个最优潮流模型,但在实际系统中,调度员很难按照最优潮流模型的方式消除所有过载,往往是通过经验或者通过局部逐渐调整的方式完成对系统状态的修正。实际调度员操作会在调整速度和精度方面存在局限性,从而可能使连锁故障持续发展,而采用最优潮流模型来模拟调度员操作可能会高估调度员的调整能力,从而使连锁故障风险评估结果出现偏差。另外在交流潮流模型下的最优潮流可能会因收敛性、非线性、非凸性等问题出现计算不收敛或者陷入局部最优解等情况造成由于计算原因而出现偏差,影响连锁故障模拟的合理性。同时调度员调整可能因为各种原因而未成功,在实际调度操作中还会有备用调度控制方案,通过紧急切负荷等操作缓解风险。
因此,传统连锁故障模型中基于最优潮流的调度员模拟在表征实际系统中调度员操作的特性方面存在明显欠缺,需要在连锁故障模型中改进调度员操作模拟,以更好地模拟实际系统连锁故障中调度员的行为,从而更为合理地模拟连锁故障的发展特征,更加准确地评估连锁故障风险。
发明内容
本发明的目的在于提供一种考虑电网调度员实际操作特性的连锁故障模拟方法,实现了对调度员在系统故障状况下操作特征的合理模拟,提高连锁故障模型的合理性,从而实现更为准确的连锁故障模拟。
本发明采用的技术方案是:
一种考虑电网调度员实际操作特性的连锁故障模拟方法,包括以下步骤:
(1)、设定连锁故障模拟天数上限kmax,支路重载负载率阈值λheavy;设系统各个节点的交流潮流基准状态为:PQ、PV节点的有功注入无功注入和电压其中j为节点编号;
(2)、初始化连锁故障天数k=1;
(3)、模拟第k天连锁故障过程;
步骤(3)具体包括:
(3.1)、随机选定第k天连锁故障模拟初始时负荷水平,将系统基准运行方式下PQ节点的有功、无功注入,以及PV节点的有功注入乘以一个[0.7,1.4]区间内平均分布的随机数,作为第k天连锁故障模拟开始时的系统状态;
(3.2)、随机模拟初始故障,随机选1条支路开断;
(3.3)、检查交流潮流是否收敛,若不收敛,表明系统发生了电压崩溃,转到步骤(3.8);若收敛,继续步骤(3.4);
(3.4)、计算所有支路的支路负载率λi,若存在支路负载率λi≥λheavy,则进入步骤(3.5),否则转到步骤(3.8);
(3.5)、模拟调度员操作的启发式模型;
(3.6)、模拟支路跳闸随机事件,计算支路负载率,并按照跳闸概率与支路负载率的关系计算各支路的跳闸概率,并通过随机采样模拟支路跳闸;
(3.7)、若有支路跳闸,则搜索总负荷量最大的电气岛继续进行连锁故障模拟,转到步骤(3.3);若没有支路跳闸,则转到步骤(3.8);
(3.8)、统计当天负荷损失、跳闸记录数据;进入步骤(4);
(4)、k=k+1;
(5)、若k≤kmax,则转到步骤(3),否则结束。
所述的步骤(3.5)具体包括以下步骤:
(3.5.1)、设调度员调整节点切有功负荷量为ΔPDj,切无功负荷量为ΔQDj,初始发电机有功功率为PGj,发电机有功出力调整量为ΔPGj,求解公式(3-1)所示的优化模型:
公式(3-1)中,CD、CG分别为目标函数中切负荷和调整发电机的权重因子,可取CD=100,CG=1;为支路i最大潮流,Pli为支路i当前潮流,λli为调度员支路i潮流风险偏好因子,为不大于1的正值,ΔPli(ΔPG,ΔPD)、ΔQli(ΔPG,ΔPD)为支路i潮流增量,分别由后面(3-10)、(3-11)定义,ΔPG,ΔPD分别为节点负荷发电调整量ΔPDj、ΔPGj的向量形式;分别为节点j发电机的出力上下限,PDj为当前负荷;
假设公式(3-1)中负荷节点的功率因数不变,即
ΔQDj=QDj/PDj·ΔPDj(3-2)
其中QDj为初始负荷无功功率;
设式(3-1)中支路i的潮流视在功率限值为支路i当前有功、无功潮流为Pli和Qli,将支路i的潮流视在功率限值转换为有功、无功限值
其中Sli为支路视在功率潮流,ΔPli(ΔPG,ΔPD)和ΔQli(ΔPG,ΔPD)分别是用节点注入变化量表示的支路有功、无功潮流变化量,由潮流方程线性化后推导而来,推导过程如下:
根据求导的链式法则,有:
其中Pl,Ql表示支路潮流,PS,QS表示节点注入,V,θ分别表示节点电压幅值和相角;
由(3-6)得到灵敏度矩阵S的计算式:
若不计支路并联支路,支路潮流有功Puv、无功Quv用状态变量表示为
其中Guv,Buv分别为支路的电导和电纳,Vu,Vv分别为节点u,v的电压幅值,θuv为节点u,v的相角差;
从而容易得到支路潮流对状态变量的偏导数,即(3-7)第一项的矩阵;而第二个矩阵实际上就是潮流方程Jacobian矩阵;这样即可求得给定支路潮流——节点注入灵敏度矩阵S;
将灵敏度矩阵S按照节点类型和变量类型进行分类,得到矩阵分别对应发电调整、负荷调整的子块:
表1灵敏度矩阵子块
这样,在调整了发电ΔPG和负荷ΔPD之后,支路潮流的近似增量可表示为
ΔPl=(Sppd+SpqdR)ΔPD+SppgΔPG (3-10)
ΔQl=(Sqpd+SqqdR)ΔPD+SqpgΔPG (3-11)
即ΔPli(ΔPG,ΔPD)和ΔQli(ΔPG,ΔPD);
(3.5.2)、若求解优化模型(3-1)不收敛,则认为调度员调整失败,转至步骤(3.5.3);若求解优化模型(3-1)收敛,则用求得的ΔPDj、ΔQDj和ΔPGj修正当前系统状态量,并转到步骤(3.5.4);
(3.5.3)、调度员应当采用次优的紧急切负荷措施,具体模拟方法为:每个负荷节点切除比例为αD的负荷,同时各个发电节点降低出力以保证功率平衡,以此达到缓解过载的目的;
(3.5.4)、若调度员连续调整次数m达到预设上限mmax,则模拟调度员未在一定时间内完成完全消除过载的目的,跳闸有可能发生,转到步骤(3.6);若m未达到预设上限mmax,则转到步骤(3.4)。
本发明的有益效果在于:
本发明提出了考虑电网调度员实际操作特性的连锁故障模拟方法,将传统的最优潮流模型改为符合实际电网调度员操作特点的启发式模型,该模型中包括基于交流系统灵敏度的调度员调整模型,以及调整不成功情况下的紧急切负荷调整模拟,体现了调度员在连锁故障中的操作特点。
本发明通过在连锁故障模拟中引入考虑实际电网调度员调整模拟,使连锁故障模拟体现了电网连锁故障过程与调度员的交互特征,更准确地反映连锁故障的发展特征,提高了连锁故障模型的实用性。在电网实际应用中,该模型用于研究调度员操作阈值与最大操作次数对风险的影响,对电网培训调度员和进行系统风险定量评估提供了实际参考。
附图说明:
图1是本发明的流程图;
图2是跳闸概率与支路负载率的关系示意图;
图3是连锁故障模拟负荷损失分布示意图。
具体实施方式:
本发明公开了一种考虑电网调度员实际操作特性的连锁故障模拟方法,以IEEE-30节点系统为例:IEEE-30节点系统包含30个节点,41条支路和6台发电机。在本算例中,将支路负载率重载阈值设为λheavy=0.8,触发保护负载率值为λprotection=1.5,保护误动概率P0=0.005;模拟天数kmax=500。对IEEE-30节点系统进行连锁故障模拟,采用考虑电网调度员实际操作特性的连锁故障模型,如图1所示,具体包括以下步骤:
(1)、设定连锁故障模拟天数上限kmax,支路重载负载率阈值λheavy;设系统各个节点的交流潮流基准状态为:PQ、PV节点的有功注入无功注入和电压其中j为节点编号;
(2)、初始化连锁故障天数k=1;
(3)、模拟第k天连锁故障过程;
(4)、统计当天损失负荷量、线路开断数等指标;
(5)、k=k+1;
(6)、若k≤kmax,则转到步骤(3),否则结束。
所述的步骤(3)具体包括以下步骤:
(3.1)、随机选定第k天连锁故障模拟初始时负荷水平,将系统基准运行方式下PQ节点的有功、无功注入,以及PV节点的有功注入乘以一个[0.7,1.4]区间内平均分布的随机数,作为第k天连锁故障模拟开始时的系统状态;
(3.2)、随机模拟初始故障,随机选1条支路开断;
(3.3)、检查交流潮流是否收敛,若不收敛,表明系统发生了电压崩溃,转到步骤(3.8);若收敛,继续步骤(3.4);
(3.4)、计算所有支路的支路负载率λi,若存在支路负载率λi≥λheavy,则进入步骤(3.5),否则转到步骤(3.8);
(3.5)、模拟调度员操作的启发式模型;
所述的步骤(3.5)具体包括以下步骤:
(3.5.1)、设调度员调整节点切有功负荷量为ΔPDj,切无功负荷量为ΔQDj,初始发电机有功功率为PGj,发电机有功出力调整量为ΔPGj,求解公式(3-1)所示的优化模型:
公式(3-1)中,CD、CG分别为目标函数中切负荷和调整发电机的权重因子,可取CD=100,CG=1;为支路i最大潮流,Pli为支路i当前潮流,λli为调度员支路i潮流风险偏好因子,为不大于1的正值,ΔPli(ΔPG,ΔPD)、ΔQli(ΔPG,ΔPD)为支路i潮流增量,分别由后面(3-10)、(3-11)定义,ΔPG,ΔPD分别为节点负荷发电调整量ΔPDj、ΔPGj的向量形式;分别为节点j发电机的出力上下限,PDj为当前负荷;
假设公式(3-1)中负荷节点的功率因数不变,即
ΔQDj=QDj/PDj·ΔPDj (3-2)
设式(3-1)中支路i的潮流视在功率限值为支路i当前有功、无功潮流为Pli和Qli,将支路i的潮流视在功率限值转换为有功、无功限值
其中Sli为线路视在功率潮流,ΔPli(ΔPG,ΔPD)和ΔQli(ΔPG,ΔPD)分别是用节点注入变化量表示的线路有功、无功潮流变化量,由潮流方程线性化后推导而来,推导过程如下:
根据求导的链式法则,有:
其中Pl,Ql表示支路潮流,PS,QS表示节点注入,V,θ分别表示节点电压幅值和相角。
由(3-6)得到灵敏度矩阵S的计算式:
若不计支路并联支路,支路潮流有功Puv、无功Quv用状态变量表示为
其中Guv,Buv分别为支路的电导和电纳,Vu,Vv分别为节点u,v的电压幅值,θuv为节点u,v的相角差。
从而容易得到支路潮流对状态变量的偏导数,即(3-7)第一项的矩阵;而第二个矩阵实际上就是潮流方程Jacobian矩阵;这样即可求得给定支路潮流——节点注入灵敏度矩阵S;
将灵敏度矩阵S按照节点类型和变量类型进行分类,得到矩阵分别对应发电调整、负荷调整的子块:
表1灵敏度矩阵子块
这样,在调整了发电ΔPG和负荷ΔPD之后,线路潮流的近似增量可表示为
ΔPl=(Sppd+SpqdR)ΔPD+SppgΔPG (3-10)
ΔQl=(Sqpd+SqqdR)ΔPD+SqpgΔPG (3-11)
即ΔPli(ΔPG,ΔPD)和ΔQli(ΔPG,ΔPD);
(3.5.2)、若求解优化模型(3-1)不收敛,则认为调度员调整失败,转至步骤(3.5.3);若求解优化模型(3-1)收敛,则用求得的ΔPDj、ΔQDj和ΔPGj修正当前系统状态量,并转到步骤(3.5.4);
(3.5.3)、调度员应当采用次优的紧急切负荷措施,具体模拟方法为:每个负荷节点切除比例为αD的负荷,同时各个发电节点降低出力以保证功率平衡,以此达到缓解过载的目的;
(3.5.4)、若调度员连续调整次数m达到预设上限mmax,则模拟调度员未在一定时间内完成完全消除过载的目的,跳闸有可能发生,转到步骤(3.6);若m未达到预设上限mmax,则转到步骤(3.4)。
(3.6)、模拟线路跳闸随机事件,计算支路负载率,并按照图2所示跳闸概率与支路负载率的关系计算各支路的跳闸概率,并通过随机采样模拟支路跳闸;图2中纵坐标为线路跳闸概率P,横坐标为线路负载率λ
(3.7)、若有支路跳闸,则搜索总负荷量最大的电气岛继续进行连锁故障模拟,转到步骤(3.3);若没有支路跳闸,则转到步骤(3.8);
(3.8)、统计当天负荷损失、跳闸记录数据。
对IEEE-30节点系统执行步骤(1)-(5),并与传统连锁故障模型(模拟调度员操作为最优潮流模型)进行对比,负荷损失分布如图3所示,图3中横坐标为连锁故障损失负荷量(标幺值),纵坐标为连锁故障损失小于对应损失负荷量的概率,该图体现了系统连锁故障的整体风险分布。
传统连锁故障模型中,交流最优潮流总共被调用1255次,其中收敛884次,未收敛371次,将近1/3的交流最优潮流调用未能正常收敛,而实际上大部分未收敛的情况是由于数值问题导致的,这样的结果会严重影响评估的准确度。由于考虑电网调度员实际操作特性的连锁故障模型中的调度员模拟启发式模型是一个线性规划问题,不存在可解性和收敛性问题,因而可以在计算中保证调度员模拟效果,提高了连锁故障模型的合理性和实用性。

Claims (2)

1.一种考虑电网调度员实际操作特性的连锁故障模拟方法,其特征在于:包括以下步骤:
(1)、设定连锁故障模拟天数上限kmax,支路重载负载率阈值λheavy;设系统各个节点的交流潮流基准状态为:PQ、PV节点的有功注入无功注入和电压其中j为节点编号;
(2)、初始化连锁故障天数k=1;
(3)、模拟第k天连锁故障过程;
步骤(3)具体包括:
(3.1)、随机选定第k天连锁故障模拟初始时负荷水平,将系统基准运行方式下PQ节点的有功、无功注入,以及PV节点的有功注入乘以一个[0.7,1.4]区间内平均分布的随机数,作为第k天连锁故障模拟开始时的系统状态;
(3.2)、随机模拟初始故障,随机选1条支路开断;
(3.3)、检查交流潮流是否收敛,若不收敛,表明系统发生了电压崩溃,转到步骤(3.8);若收敛,继续步骤(3.4);
(3.4)、计算所有支路的支路负载率λi,若存在支路负载率λi≥λheavy,则进入步骤(3.5),否则转到步骤(3.8);
(3.5)、模拟调度员操作的启发式模型;
(3.6)、模拟支路跳闸随机事件,计算支路负载率,并按照跳闸概率与支路负载率的关系计算各支路的跳闸概率,并通过随机采样模拟支路跳闸;
(3.7)、若有支路跳闸,则搜索总负荷量最大的电气岛继续进行连锁故障模拟,转到步骤(3.3);若没有支路跳闸,则转到步骤(3.8);
(3.8)、统计当天负荷损失、跳闸记录数据;进入步骤(4);
(4)、k=k+1;
(5)、若k≤kmax,则转到步骤(3),否则结束。
2.根据权利要求1所述的考虑电网调度员实际操作特性的连锁故障模拟方法,其特征在于:所述的步骤(3.5)具体包括以下步骤:
(3.5.1)、设调度员调整节点切有功负荷量为ΔPDj,切无功负荷量为ΔQDj,初始发电机有功功率为PGj,发电机有功出力调整量为ΔPGj,求解公式(3-1)所示的优化模型:
min C D ΔP D j + C G Σ | ΔP G j | s . t . - λ l i P l i max - P l i ≤ ΔP l i ( ΔP G , ΔP D ) ≤ λ l i P l i max - P l i - λ l i Q l i max - Q l i ≤ ΔQ l i ( ΔP G , ΔP D ) ≤ λ l i Q l i max - Q l i P G j min - P G j ≤ ΔP G j ≤ P G j max - P G j 0 ≤ ΔP D j ≤ P D j - - - ( 3 - 1 )
公式(3-1)中,CD、CG分别为目标函数中切负荷和调整发电机的权重因子,可取CD=100,CG=1;为支路i最大潮流,Pli为支路i当前潮流,λli为调度员支路i潮流风险偏好因子,为不大于1的正值,ΔPli(ΔPG,ΔPD)、ΔQli(ΔPG,ΔPD)为支路i潮流增量,分别由后面(3-10)、(3-11)定义,ΔPG,ΔPD分别为节点负荷发电调整量ΔPDj、ΔPGj的向量形式;分别为节点j发电机的出力上下限,PDj为当前负荷;
假设公式(3-1)中负荷节点的功率因数不变,即
ΔQDj=QDj/PDj·ΔPDj (3-2)
其中QDj为初始负荷无功功率;
设式(3-1)中支路i的潮流视在功率限值为支路i当前有功、无功潮流为Pli和Qli,将支路i的潮流视在功率限值转换为有功、无功限值
S l i = P l i 2 + Q l i 2 - - - ( 3 - 3 )
P l i max = | P l i | S l i S l i max - - - ( 3 - 4 )
Q l i max = | Q l i | S l i S l i max - - - ( 3 - 5 )
其中Sli为支路视在功率潮流,ΔPli(ΔPG,ΔPD)和ΔQli(ΔPG,ΔPD)分别是用节点注入变化量表示的支路有功、无功潮流变化量,由潮流方程线性化后推导而来,推导过程如下:
根据求导的链式法则,有:
∂ P l ∂ V ∂ P l ∂ θ ∂ Q l ∂ V ∂ Q l ∂ θ = ∂ P l ∂ P S ∂ P l ∂ Q S ∂ Q l ∂ P S ∂ Q l ∂ Q S ∂ P S ∂ V ∂ P S ∂ θ ∂ Q S ∂ V ∂ P S ∂ θ - - - ( 3 - 6 )
其中Pl,Ql表示支路潮流,PS,QS表示节点注入,V,θ分别表示节点电压幅值和相角;
由(3-6)得到灵敏度矩阵S的计算式:
S = ∂ P l ∂ P S ∂ P l ∂ Q S ∂ Q l ∂ P S ∂ Q l ∂ Q S = ∂ P l ∂ V ∂ P l ∂ θ ∂ Q l ∂ V ∂ Q l ∂ θ ∂ P S ∂ V ∂ P S ∂ θ ∂ Q S ∂ V ∂ P S ∂ θ - 1 - - - ( 3 - 7 )
若不计支路并联支路,支路潮流有功Puv、无功Quv用状态变量表示为
P u v = - G u v V u 2 + G u v V u V v cosθ u v + B u v V u V v sinθ u v - - - ( 3 - 8 )
Q u v = B u v V u 2 - B u v V u V v cosθ u v + G u v V u V v sinθ u v - - - ( 3 - 9 )
其中Guv,Buv分别为支路的电导和电纳,Vu,Vv分别为节点u,v的电压幅值,θuv为节点u,v的相角差;
从而容易得到支路潮流对状态变量的偏导数,即(3-7)第一项的矩阵;而第二个矩阵实际上就是潮流方程Jacobian矩阵;这样即可求得给定支路潮流——节点注入灵敏度矩阵S;
将灵敏度矩阵S按照节点类型和变量类型进行分类,得到矩阵分别对应发电调整、负荷调整的子块:
表1 灵敏度矩阵子块
这样,在调整了发电ΔPG和负荷ΔPD之后,支路潮流的近似增量可表示为
ΔPl=(Sppd+SpqdR)ΔPD+SppgΔPG (3-10)
ΔQl=(Sqpd+SqqdR)ΔPD+SqpgΔPG (3-11)
即ΔPli(ΔPG,ΔPD)和ΔQli(ΔPG,ΔPD);
(3.5.2)、若求解优化模型(3-1)不收敛,则认为调度员调整失败,转至步骤(3.5.3);若求解优化模型(3-1)收敛,则用求得的ΔPDj、ΔQDj和ΔPGj修正当前系统状态量,并转到步骤(3.5.4);
(3.5.3)、调度员应当采用次优的紧急切负荷措施,具体模拟方法为:每个负荷节点切除比例为αD的负荷,同时各个发电节点降低出力以保证功率平衡,以此达到缓解过载的目的;
(3.5.4)、若调度员连续调整次数m达到预设上限mmax,则模拟调度员未在一定时间内完成完全消除过载的目的,跳闸有可能发生,转到步骤(3.6);若m未达到预设上限mmax,则转到步骤(3.4)。
CN201410654244.1A 2014-11-17 2014-11-17 考虑电网调度员实际操作特性的连锁故障模拟方法 Active CN104332974B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410654244.1A CN104332974B (zh) 2014-11-17 2014-11-17 考虑电网调度员实际操作特性的连锁故障模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410654244.1A CN104332974B (zh) 2014-11-17 2014-11-17 考虑电网调度员实际操作特性的连锁故障模拟方法

Publications (2)

Publication Number Publication Date
CN104332974A CN104332974A (zh) 2015-02-04
CN104332974B true CN104332974B (zh) 2017-04-19

Family

ID=52407651

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410654244.1A Active CN104332974B (zh) 2014-11-17 2014-11-17 考虑电网调度员实际操作特性的连锁故障模拟方法

Country Status (1)

Country Link
CN (1) CN104332974B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106295094B (zh) * 2015-05-15 2019-01-11 中国电力科学研究院 一种考虑继电保护隐性故障的电网薄弱环节识别仿真方法
CN110689178B (zh) * 2019-09-18 2023-07-18 国网湖南省电力有限公司 电网的输电线路的山火连锁故障风险的分析方法及系统
CN112564113B (zh) * 2020-12-15 2023-05-12 深圳供电局有限公司 一种基于近似化潮流电压稳定的收敛性调整方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11206019A (ja) * 1998-01-19 1999-07-30 Toshiba Corp 電力系統安定化機器の制御装置及び制御方法並びに電力系統安定化機器制御プログラムを記録した記録媒体
CN101232182A (zh) * 2008-01-18 2008-07-30 清华大学 三维协调的电网能量管理系统及对电网的评估及控制方法
CN101944766A (zh) * 2010-06-01 2011-01-12 中国电力科学研究院 一种电力系统运行及安全监控仿真系统
CN102737286A (zh) * 2012-04-23 2012-10-17 中国电力科学研究院 一种地区电网在线风险分析系统和方法
CN103326356A (zh) * 2013-06-07 2013-09-25 广东电网公司电力调度控制中心 考虑人工调度知识的电网关键断面极限传输容量计算方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11206019A (ja) * 1998-01-19 1999-07-30 Toshiba Corp 電力系統安定化機器の制御装置及び制御方法並びに電力系統安定化機器制御プログラムを記録した記録媒体
CN101232182A (zh) * 2008-01-18 2008-07-30 清华大学 三维协调的电网能量管理系统及对电网的评估及控制方法
CN101944766A (zh) * 2010-06-01 2011-01-12 中国电力科学研究院 一种电力系统运行及安全监控仿真系统
CN102737286A (zh) * 2012-04-23 2012-10-17 中国电力科学研究院 一种地区电网在线风险分析系统和方法
CN103326356A (zh) * 2013-06-07 2013-09-25 广东电网公司电力调度控制中心 考虑人工调度知识的电网关键断面极限传输容量计算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"交直流系统连锁故障模型及停电风险分析";余晓鹏等;《电力系统自动化》;20141010;第38卷(第19期);第33-39页 *
"适用于连锁故障的交直流电网静态等值方法";张振安等;《电工电能新技术》;20140331;第33卷(第3期);第1-6页 *

Also Published As

Publication number Publication date
CN104332974A (zh) 2015-02-04

Similar Documents

Publication Publication Date Title
Xu et al. Robust transient stability-constrained optimal power flow with uncertain dynamic loads
Leon et al. Hierarchical wide-area control of power systems including wind farms and FACTS for short-term frequency regulation
CN104901306B (zh) 一种考虑连锁故障的电网运行安全裕度计算方法
Wang et al. Probabilistic reliability evaluation including adequacy and dynamic security assessment
CN105958538B (zh) 基于蒙特卡洛法的配电网孤岛划分方法
CN104638646B (zh) 基于复杂网络理论计及校正控制的电网连锁故障模拟方法
Yan et al. Adaptive load shedding method based on power imbalance estimated by ANN
CN106291255A (zh) 一种含逆变型分布式电源的配电网故障计算通用方法
CN104332974B (zh) 考虑电网调度员实际操作特性的连锁故障模拟方法
Souxes et al. Effect of stochasticity on voltage stability support provided by wind farms: Application to the Hellenic interconnected system
CN104537428A (zh) 一种计及风电接入不确定性的经济运行评估方法
CN104037790A (zh) 一种基于时序蒙特卡罗模拟的新能源接纳能力评估方法
Mahmoudi et al. Application of distributed control to mitigate disturbance propagations in large power networks
Moya et al. A hierarchical framework for demand-side frequency control
CN106099991B (zh) 一种电网风电消纳能力评估方法
CN102664417A (zh) 二级电压控制方法和装置
CN106655201A (zh) 一种基于安全域的电力系统热稳定安全优化控制方法
CN105224812B (zh) 一种负荷模型中的静态负荷频率因子聚合方法
CN102638036A (zh) 一种用于计划安全校核的快速支路过载辅助决策方法
Li et al. A neural power system stabilizer of DFIGs for power system stability support
CN112560227B (zh) 一种电网调度日前计划仿真预演方法
CN105281371A (zh) 一种考虑风力发电的可伸缩有功静态安全域
Weng et al. Fuzzy power flow solution considering wind power variability and uncertainty
CN103078347B (zh) 降低电压越限风险的冲击负荷消纳方法
CN110071499B (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
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20170320

Address after: 450052 Songshan, Zhengzhou, Henan District No. 27 South Road, No. 85

Applicant after: Electric Power Research Institute, State Grid Henan Electric Power Company

Applicant after: Tsinghua University

Applicant after: State Grid Corporation of China

Address before: 100031 Xicheng District West Chang'an Avenue, No. 86, Beijing

Applicant before: State Grid Corporation of China

Applicant before: Electric Power Research Institute, State Grid Henan Electric Power Company

Applicant before: Tsinghua University

GR01 Patent grant
GR01 Patent grant