CN103311926A - 基于统一潮流控制器的电力系统连锁故障模拟方法 - Google Patents
基于统一潮流控制器的电力系统连锁故障模拟方法 Download PDFInfo
- Publication number
- CN103311926A CN103311926A CN2013102822006A CN201310282200A CN103311926A CN 103311926 A CN103311926 A CN 103311926A CN 2013102822006 A CN2013102822006 A CN 2013102822006A CN 201310282200 A CN201310282200 A CN 201310282200A CN 103311926 A CN103311926 A CN 103311926A
- Authority
- CN
- China
- Prior art keywords
- node
- power system
- subsystem
- electric power
- jie
- 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
Images
Classifications
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/80—Technologies aiming to reduce greenhouse gasses emissions common to all road transportation technologies
- Y02T10/82—Elements for improving aerodynamics
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
一种基于统一潮流控制器的电力系统连锁故障模拟方法,属于电力系统安全可靠性防护技术领域。本发明利用计算机,通过程序,先确定电力系统中统一潮流控制器的安装位置、再确定电力系统各节点的有功潮流介数阈值、然后进行基于统一潮流控制器的电力系统连锁故障模拟、最后计算反映连锁故障严重程度的衡量指标。本发明考虑了潮流的方向性和系统的运行方式,优化了装设统一潮流控制器的安装位置,在满足电力系统安全运行约束下有效进行线路参数的调整,能够更加真实地反映电力系统连锁故障的传播过程。本发明可广泛应用于电力系统的连锁故障模拟,特别适用于大型复杂电力系统的连锁故障模拟,为预防电力系统连锁故障的发生提供科学依据。
Description
技术领域
本发明属于电力系统安全可靠性防护技术领域,具体涉及一种基于统一潮流控制器的电力系统连锁故障模拟方法。
背景技术
随着经济的发展和电力需求的增长,电力系统的规模不断扩大和结构日趋复杂,电力系统运行的安全可靠性越来越受到人们的高度重视。近年来,电力系统由局部故障传播造成的大停电事故时有发生,通过对这些电力系统大停电事故的统计分析表明,大停电事故多由连锁故障引起。电力系统连锁故障是指,电力系统的部分设备发生偶然故障后,电力系统的潮流发生转移,潮流在电力系统局部范围内重新分配时可能导致其它电力设备故障退出运行,最终导致故障大范围传播甚至波及整个电力系统。电力系统连锁故障严重影响电力系统的安全运行和可靠供电,造成的经济损失和社会影响巨大。因此,构建合理的模拟电力系统连锁故障的模型,研究电力系统连锁故障的模拟方法,为降低连锁故障发生的概率及影响范围提供可靠依据具有十分重要的意义。
现有电力系统连锁故障模拟方法,如2010年第30卷第13期《中国电机工程学报》中“基于电气介数的电网连锁故障传播机制与积极防御”一文中,公开的方法是基于电气介数的电力系统连锁故障模拟方法。其具体方法是:首先,通过计算故障前电力系统各节点的电气介数确定该节点介数的阈值(运行阈值和极限阈值);其次,随机或有目的地移除系统中的一个节点来模拟初始故障的发生,计算故障后电力系统各节点的电气介数。然后,将各节点故障后的电气介数与其介数阈值进行比较来判断电力系统是否发生连锁故障,当某节点故障后的电气介数大于其极限阈值时,则将该节点和与其相连的所有线路移除来模拟连锁故障的传播;当电气介数处于运行阈值和极限阈值之间时,则按比例调整与该节点相连的线路参数来降低该节点的电气介数,以缓解连锁故障的进一步传播;当所有节点的电气介数均小于运行阈值时,连锁故障模拟结束。最后,通过计算故障衡量指标来评估连锁故障的严重程度。该方法的主要缺点是:①用电气介数判断节点的重要性,忽略了电力系统中实际电源功率的大小、方向和无功电压等因素,不能体现电力系统的实际运行情况;②在连锁故障模拟过程中,对缓解连锁故障进一步发生的线路参数的调整方式,采用的是每回线路参数均可按比例进行调整的假设,与实际系统运行情况不符;③虽然文中提到可采用灵活交流输电(FACTS)元件调整相关线路的参数,只是理论设想,没有具体实施过程,并且在电力系统的每回输电线路上均安装灵活交流输电(FACTS)元件,实际中难以实现;④调整线路参数时没有考虑系统安全约束,不能反映电力系统故障后的实际运行状况。因此,用该方法模拟电力系统连锁故障不能真实反映连锁故障的传播过程,从而不能为电力系统的安全可靠性防护提供依据,不能有效地避免电力系统大面积停电的发生。
发明内容
本发明的目的是针对现有电力系统连锁故障模拟方法的不足,提出一种基于统一潮流控制器的电力系统连锁故障模拟方法,具有在连锁故障模拟过程中,考虑了电力系统潮流的大小及方向和电力系统的实际运行方式,优化了装设统一潮流控制器(UPFC)的安装位置,在满足电力系统安全运行约束要求下有效实现线路参数的调整,能更好地反映电力系统连锁故障的传播过程等特点,从而为降低连锁故障发生的概率、控制连锁故障对电力系统影响的范围和有效预防电力系统连锁故障的发生提供科学依据,能有效地避免电力系统大面积停电的发生。
实现本发明目的之技术方案是:一种基于统一潮流控制器的电力系统连锁故障模拟方法,利用计算机,通过程序,先确定电力系统中UPFC的装设位置、再确定装设有UPFC的电力系统各节点的有功潮流介数阈值、然后对装设有UPFC的电力系统进行连锁故障模拟、最后计算反映连锁故障严重程度的衡量指标。所述方法的具体步骤如下:
(1)确定电力系统中UPFC的装设位置
1)输入基本参数
首先输入电力系统的基本参数和统一潮流控制器(UPFC)的基本参数。所述的电力系统的基本参数包括节点编号、节点类型、节点对应电压等级、各节点的有功负荷(PL)和无功负荷(QL)、与发电机连接的节点的编号、各发电机输出的有功功率(PG)和无功功率(QG)、各线路首末端节点编号、线路电阻(R)、线路电抗(X)和线路电纳(B)、线路的额定电压(UB)、基准功率(SB)、节点的运行系数(α)、节点的极限系数(β);所述的UPFC控制器的基本参数包括在系统中安装UPFC的台数(m)、每台UPFC的幅值控制参数(kse)、相角控制参数和无功控制参数(Qsh)。
2)求解电力系统实际运行的交流潮流
第(1)-1)步完成后,根据第(1)-1)步输入的电力系统的基本参数,用常规的交流潮流计算方法,求解电力系统实际运行的交流潮流,确定电力系统中各发电机的输出有功功率、各节点的电压、各条线路的交流潮流的大小及其方向、功率损耗及充电功率。
3)确定电力系统等效无损网络中各输电路径的有功功率
第(1)-2)步完成后,利用本申请人申请的申请号为201310213277.8的“基于有功潮流介数的电力系统关键节点辨识方法”的第(2)步,确定所述的等效无损网络中各输电路径Bm的有功功率。即:首先将电力系统等效成无损网络,然后确定所述的等效无损网络中的输电路径,再计算所述的等效无损网络中各节点的有功功率分配比例因子,最后确定所述的等效无损网络中各输电路径的有功功率。
4)计算电力系统等效无损网络中各线路的有功潮流介数
第(1)-3)步完成后,将第(1)-3)步得到的通过线路l的所有输电路径的有功功率加权求和来计算该线路的有功潮流介数,用其线路的有功潮流介数来反映调整该线路参数对潮流影响的程度,即用线路的有功潮流介数来量化“发电-负荷”节点对之间的线路参数对潮流控制的影响。由于当通过所述的等效无损网络中各线路的输电路径越多和输电路径的有功功率越大时,调整该线路参数对潮流的影响越大,因此线路的有功潮流介数越大,表明该线路对潮流传播的贡献越大,在该线路装设统一潮流控制器(UPFC)后对电力系统潮流的调节能力也越大。计算线路l的有功潮流介数的具体公式如下:
式中:Bf(l)为线路l的有功潮流介数;G为发电机节点集合;L为负荷节点集合;B为输电路径集合;l、y和z分别为输电路径Bm的线路、发电机节点和负荷节点;Wy为发电机节点y的权重,Wy的值为发电机节点y的有功出力值;Wz为负荷节点z的权重,Wz的值为负荷节点z的有功负荷;Pm(l),y,z为输电路径Bm经过线路l的输电路径有功功率。
5)确定电力系统中UPFC的装设位置
第(1)-4)步完成后,先设定电力系统中装设UPFC的线路数量的范围为1~M;然后将第(1)-4)步得到的电力系统等效无损网络中各线路的有功潮流介数从大到小排序,并确定排在范围1~M内的前m条线路的首端装设UPFC,以保证电力系统运行的经济性和可行性。M和m的取值根据电力系统的运行成本和实际操作情况确定。
(2)确定装设有UPFC的电力系统各节点的有功潮流介数阈值
1)计算装设有UPFC线路的等效节点有功功率注入增量和无功功率注入增量
第(1)-5)步完成后,利用装设有UPFC的线路i-j两端节点的电压计算其节点等效有功功率注入增量和节点等效无功功率注入增量,具体计算公式为:
式中:ΔPi和ΔQi分别为装设有UPFC的线路i-j在节点i的等效有功功率注入增量和等效无功功率注入增量;ΔPj和ΔQj分别为装设有UPFC的线路i-j在节点j的等效有功功率注入增量和等效无功功率注入增量;gij和bij分别为线路i-j的电导和电纳;Bc为线路i-j的对地电纳;θij为节点i和节点j的之间的相角差,θij=θi-θj;Ui和Uj分别为节点i和节点j的电压幅值;kse为UPFC的幅值控制参数;为UPFC的相角控制参数;Qsh为UPFC的无功控制参数。
2)求解装设UPFC后的电力系统的交流潮流
第(2)-1)步完成后,首先,将按第(2)-1)步得到的装有UPFC的线路i-j两端节点的节点等效有功功率注入增量ΔPi、ΔPj和节点等效无功功率注入增量ΔQi、ΔQj叠加到相应节点的节点功率方程中。具体计算公式如下:
式中:Pgi和Qgi分别为节点i的发电机有功出力和无功出力;Pdi和Qdi分别为节点i的有功负荷和无功负荷;Gij和Bij分别为节点i和节点j的之间的互电导和互电纳;j∈i表示∑号后的标号为j的节点必须直接和节点i相连,并包括j=i的情况;Lu为所有与节点i相连且装有UPFC的线路集合;其他参数的含义与公式(2)相同。
然后,用常规的交流潮流计算方法,求解电力系统实际运行的交流潮流,确定电力系统中各发电机输出有功功率、各节点的电压、各条线路的交流潮流大小及其方向、功率损耗及充电功率。
3)计算装设有UPFC的电力系统等效无损网络中各节点的有功潮流介数
第(2)-2)步完成后,根据本申请人申请的申请号为201310213277.8的“基于有功潮流介数的电力系统关键节点辨识方法”的第(2)步和第(3)-1)步,先将所述的电力系统等效成无损网络,然后确定所述的等效无损网络中的输电路径,再计算所述的等效无损网络中各节点的有功功率分配比例因子,然后确定所述的等效无损网络中各输电路径的有功功率,最后计算所述的等效无损网络各节点的有功潮流介数。
4)确定装设有UPFC的电力系统各节点的有功潮流介数阈值
第(2)-3)步完成后,确定电力系统中各节点的有功潮流介数运行阈值和有功潮流介数极限阈值,计算公式为:
式中:Sn为节点n的有功潮流介数运行阈值;Γn为节点n的有功潮流介数极限阈值;Bf(n)为故障前节点n的有功潮流介数。
(3)对装设有UPFC的电力系统进行连锁故障模拟
第(2)步完成后,先通过移除所述的电力系统中的任意节点和与其相连的所有线路来模拟实际电力系统中发电厂或变电站发生故障。再通过计算所述的故障电力系统各节点的有功潮流介数来判断系统是否发生连锁故障,包括形成所述的故障电力系统的子系统及个数、确定所述的故障电力系统的各子系统等效无损网络的节点有功潮流介数、运用得到的节点有功潮流介数判断系统是否发生连锁故障,来模拟连锁故障的传播和终止。具体步骤如下:
1)运用广度优先搜索方法形成所述的故障电力系统的子系统及个数
首先从电力系统中的任一节点n出发,搜索其相邻的节点构成集合A0,将节点n和集合A0中的元素放入集合B;再搜索与A0中各节点相邻且不属于集合B的节点构成集合A1,将集合A1中的元素放入集合B中;继续搜索与A1中节点相邻且不属于集合B的节点集合A2;如此搜索下去,直到不能再搜索到新的节点集合为止;本次搜索到的所有节点集合B形成子系统I1。然后在本次未搜索到的节点中,重复上述搜索过程,依次形成子系统I2,I3,……,IK,直到不存在未被搜索到的节点为止。K为形成的子系统个数。
2)确定所述的故障电力系统的各子系统等效无损网络的节点有功潮流介数
①第(3)-1)步完成后,先将第(3)-1)步得到的子系统I1,I2,……,IK中装设有UPFC装置的子系统构成集合Du,没有装设UPFC装置的子系统构成集合D;
②第(3)-2)-①步完成后,对集合Du和集合D中有发电机节点但没有平衡节点的子系统,将该子系统中发电机装机容量最大的节点作为平衡节点;
③第(3)-2)-②步完成后,对所述集合Du的各子系统,为减小故障电力系统的各子系统中线路潮流的越限值,用计及UPFC的交流最优潮流方法确定各子系统的各发电机输出有功功率、各节点的电压、各条线路的交流潮流大小及其方向、功率损耗及充电功率。具体步骤如下:
首先,对第(3)-2)-①步中所述的集合Du中的各子系统,按照公式(2)计算装设有UPFC的线路i-j两端节点的节点等效有功功率注入增量和节点等效无功功率注入增量。
然后,运用常规内点法求解计及UPFC的交流最优潮流模型。该模型的目标函数为所有线路潮流越限值之和最小,约束条件包括潮流方程,发电机有功和无功出力、节点电压模值、UPFC幅值及相角控制参数和无功控制参数可调范围限制等约束。具体公式如下:
Pgi,min≤Pgi≤Pgi,max (8)
Qgi,min≤Qgi≤Qgi,max (9)
Ui,min≤Ui≤Ui,max (10)
0≤kse≤kse,max (11)
-Qsh,max≤Qsh≤Qsh,max (13)
式中:T为所有线路潮流越限值之和;N为子系统中所有节点的集合;Pij和Qij分别为流过线路i-j首端的有功功率和无功功率,Tij,max为线路i-j的额定视在载流容量;Lu为所有与节点i相连且装有UPFC的线路集合;Pgi,min、Qgi,min分别为节点i的发电机有功出力和无功出力下限;Pgi,max和Qgi,max分别为节点i的发电机有功出力和无功出力上限;Ui,min和Ui,max分别为节点i的电压幅值下限和上限;kse和kse,max分别为UPFC的幅值控制参数及最大值;为UPFC的相角控制参数;Qsh和Qshmax分别为UPFC的无功控制参数及最大值;其他参数的含义与公式(3)相同。
④第(3)-2)-③步完成后,对所述的集合D的各子系统用常规的交流潮流计算方法,求解系统实际运行的交流潮流,确定各子系统中各发电机输出有功功率、各节点的电压、各条线路的交流潮流大小及其方向、功率损耗及充电功率。
⑤第(3)-2)-④步完成后,转到第(2)-3)步,计算所述的故障电力系统的各子系统等效无损网络的节点有功潮流介数。
3)模拟连锁故障的传播和终止
第(3)-2)-⑤步完成后,将第(3)-2)-⑤步得到的故障电力系统的各子系统中各节点的有功潮流介数与其有功潮流介数阈值进行比较,来模拟连锁故障的传播和终止。
①对于第(3)-2)-①步所述的集合Du中的各子系统,模拟连锁故障的步骤如下:
(a)对每一条节点设置一个标志,并置为0;
(b)第(3)-3)-①-(a)步完成后,对该子系统中的每个节点进行判断,当该节点的有功潮流介数处于其有功潮流介数运行阈值和有功潮流介数极限阈值之间时,则将标志置为1;当该节点的有功潮流介数大于其有功潮流介数极限阈值时,则移除该节点和与其相连的所有线路,将标志置为2;直到判断完该子系统中所有节点为止;
(c)第(3)-3)-①-(b)步完成后,当所有节点的标志均为0时,则该子系统的连锁故障模拟过程结束;当所有节点的标志均不为2,且有1时,则重复第(3)-2)步至(3)-3)步;当所有节点的标志中有一个为2时,则重复第(3)-1)步至(3)-3)步;直到该子系统的各节点的有功潮流介数均小于其有功潮流介数运行阈值为止。
②对于第(3)-2)-①步所述的集合D中的各子系统,模拟连锁故障的步骤如下:
(a)当该子系统中只包含一个节点,或该子系统各节点的有功潮流介数均未超过其有功潮流介数极限阈值时,则该子系统的连锁故障模拟过程结束,否则转到第(3)-3)-②-(b)步;
(b)第(3)-3)-②-(a)步完成后,当该子系统中有节点的有功潮流介数大于其有功潮流介数极限阈值的节点时,则移除该节点和与其相连的所有线路;
(c)第(3)-3)-②-(b)步完成后,首先用第(3)-1)步的方法确定该子系统的新形成的各子系统;其次用第(3)-2)-②步确定新的各子系统中的平衡节点,再用第(3)-2)-④确定新的各子系统中各发电机输出有功功率、各节点的电压、各条线路的交流潮流大小及其方向、功率损耗及充电功率;然后用第(3)-2)-⑤步的方法计算新的各子系统等效无损网络中各节点的有功潮流介数;最后转至第(3)-3)-②-(a)步,直到各子系统等效无损网络中各节点的有功潮流介数均未超过其有功潮流介数极限阈值为止。
(4)计算反映连锁故障严重程度的衡量指标
第(3)-3)步完成后,对第(3)步所述的故障电力系统的各子系统,先根据电力系统在发生故障后应该满足运行的安全约束条件,当其不能满足运行的安全约束条件时,应该通过调整发电机出力和统一潮流控制器的控制参数使故障电力系统的运行尽量满足运行的安全约束条件,当调整控制后仍不能使故障电力系统运行在安全约束范围内时,应该使其失负荷最小的原则,计算所述的故障电力系统中各子系统的最小有功切负荷量。再通过计算满足安全约束条件下的最小失负荷百分比来反映连锁故障严重程度。具体步骤如下:
1)计算所述的故障电力系统中各子系统的最小有功切负荷量
运用常规的内点法求解计及UPFC的交流最优负荷削减模型来计算所述的故障电力系统的每一个子系统满足运行安全约束条件的最小有功切负荷量。该模型的目标函数为每一个子系统削减的负荷总量最小,约束条件包括潮流方程,发电机有功和无功出力、节点电压模值、UPFC幅值及相角控制参数和无功控制参数可调范围限制、削减负荷前后功率因数保持恒定、节点削减负荷量和线路潮流等约束,具体计算公式为:
Pgi,min≤Pgi≤Pgi,max (17)
Qgi,min≤Qgi≤Qgi,max (18)
Ui,min≤Ui≤Ui,max (19)
0≤kse≤kse,max (20)
-Qsh,max≤Qsh≤Qsh,max (22)
0≤Pci≤Pdi (24)
0≤Qci≤Qdi (25)
式中,C为子系统中各节点的有功切负荷总量;Pci为节点i的有功负荷削减量;Qci为节点i的无功负荷削减量;L为子系统中所有节点的集合;其他参数的含义与公式(5)到公式(13)相同。
2)计算反映连锁故障严重程度的衡量指标
第(4)-1)步完成后,计算衡量连锁故障严重程度的最小失负荷百分比指标,其具体计算公式为:
式中:Lmin为最小失负荷百分比;S为故障后系统的子系统个数;Cs为子系统s削减有功负荷的总量;L0为电力系统故障前的有功负荷总和。
根据不同节点遭受攻击后电力系统的最小失负荷百分比,就得到连锁故障进行风险评估的结果,进而为制定降低连锁故障风险的预防策略提供了科学依据,能有效地避免电力系统大面积停电的发生。
本发明采用上述技术方案后,主要有以下效果:
1.本发明的模拟方法是基于电力系统实际运行的交流潮流分布,能同时考虑潮流的方向性和电力系统的运行方式,与现有电气介数模拟法仅从网络拓扑角度分析电力系统中连锁故障的传播机理相比,更加符合电力系统的实际运行情况,实用性强。
2.本发明通过在实际电力系统中安装UPFC来调节线路的潮流,克服了现有电气介数模拟法仅从理论上设想通过调整相关线路参数来实现线路潮流调节的缺点,才能反映电力系统的实际运行情况,因而本发明更具有实用价值,进一步提高了实用性。
3.本发明根据输电线路的有功潮流介数确定UPFC的优化安装位置,克服了现有电气介数模拟法在每回输电线路都安装FACTS元件难以实现的缺点,因而才具有可操作性。
4.本发明的模拟方法通过UPFC调节线路参数和潮流时,考虑了电力系统的安全约束,与现有电气介数模拟法按照固定比例调整相关线路参数相比,能更加真实地反映电力系统连锁故障的传播过程。
本发明可广泛应用于电力系统的连锁故障模拟,特别适用于大型复杂电力系统的连锁故障模拟。为降低连锁故障发生的概率、控制连锁故障对电力系统影响的范围和有效预防电力系统连锁故障的发生提供科学依据。
附图说明
图1为本发明的程序流程图;
图2为实施例的IEEE39节点电力系统接线原理图。
图中:G为发电机节点,1~39为节点。
具体实施方式
下面结合具体实施方式进一步说明本发明。
实施例
如图1、2所示,一种基于统一潮流控制器的电力系统连锁故障模拟方法的具体步骤如下:
(1)确定电力系统中UPFC的装设位置
1)输入基本参数
首先输入电力系统的基本参数和统一潮流控制器(UPFC)的基本参数。所述的电力系统的基本参数包括节点编号(1,2,……,39)、节点类型、节点对应电压等级、各节点的有功负荷(PL)和无功负荷(QL)、与发电机连接的节点的编号、各发电机输出的有功功率(PG)和无功功率(QG)、各线路首末端节点编号、线路电阻(R)、线路电抗(X)和线路电纳(B)、线路的额定电压(UB)、基准功率(SB)、节点的运行系数(α=1.5)、节点的极限系数(β=3);所述的UPFC控制器的基本参数包括在系统中安装UPFC的台数(m=5)、每台UPFC的幅值控制参数(kse=0.4)、相角控制参数和无功控制参数(Qsh=0.1)。
2)求解电力系统实际运行的交流潮流
第(1)-1)步完成后,根据第(1)-1)步输入的电力系统的基本参数,用常规的交流潮流计算方法,求解电力系统实际运行的交流潮流,确定电力系统中各发电机的输出有功功率、各节点的电压、各条线路的交流潮流的大小及其方向、功率损耗及充电功率。
以图2的电力系统为例,说明图中部分节点的发电机输出有功功率,如表1所示。
表1部分节点的发电机输出有功功率
节点编号 | 节点类型 | 输出有功功率 |
31 | 平衡节点 | 6.8142 |
32 | PV节点 | 6.50 |
33 | PV节点 | 6.32 |
3)确定电力系统等效无损网络中各输电路径的有功功率
第(1)-2)步完成后,利用本申请人申请的申请号为201310213277.8的“基于有功潮流介数的电力系统关键节点辨识方法”的第(2)步,确定所述的等效无损网络中各输电路径Bm的有功功率。即:首先将电力系统等效成无损网络,然后确定所述的等效无损网络中的输电路径,再计算所述的等效无损网络中各节点的有功功率分配比例因子,最后确定所述的等效无损网络中各输电路径的有功功率。
以图2的电力系统为例,说明图中发电机节点31和发电机节点38到各负荷节点的输电路径及有功功率,如表2所示。节点31的发电机输出的有功功率标幺值为6.78,节点38的发电机输出的有功功率标幺值为8.30。
表2发电机节点31、38到各负荷节点的输电路径的有功功率
4)计算电力系统等效无损网络中各线路的有功潮流介数
第(1)-3)步完成后,将第(1)-3)步得到的通过线路l的所有输电路径的有功功率加权求和来计算该线路的有功潮流介数,用其线路的有功潮流介数来反映调整该线路参数对潮流影响的程度,即用线路的有功潮流介数来量化“发电-负荷”节点对之间的线路参数对潮流控制的影响。由于当通过所述的等效无损网络中各线路的输电路径越多和输电路径的有功功率越大时,调整该线路参数对潮流的影响越大,因此线路的有功潮流介数越大,表明该线路对潮流传播的贡献越大,在该线路装设统一潮流控制器(UPFC)后对电力系统潮流的调节能力也越大。计算线路l的有功潮流介数的具体公式如下:
式中:Bf(l)为线路l的有功潮流介数;G为发电机节点集合;L为负荷节点集合;B为输电路径集合;l、y和z分别为输电路径Bm的线路、发电机节点和负荷节点;Wy为发电机节点y的权重,Wy的值为发电机节点y的有功出力值;Wz为负荷节点z的权重,Wz的值为负荷节点z的有功负荷;Pm(l),y,z为输电路径Bm经过线路l的输电路径有功功率。
以图2的电力系统为例,说明图中部分线路的有功潮流介数,如表3所示。
表3部分线路的有功潮流介数
线路编号 | 线路首-末端节点 | 线路有功潮流介数 |
1 | 1-2 | 5.127 |
14 | 6-31 | 19.805 |
18 | 10-11 | 9.645 |
25 | 15-16 | 4.751 |
33 | 19-33 | 15.060 |
37 | 22-35 | 14.606 |
44 | 26-29 | 3.319 |
45 | 28-29 | 3.491 |
由表3可见,线路6-31的有功潮流介数最大,为19.805。这是因为表2中的发电机节点31到负荷节点有11条输电路径,其中有10条输电路径通过线路6-31,最大输电路径有功功率为输电路径11,其标幺值为2.6315,因此线路6-31有功潮流介数较大。在表3中,线路26-29的有功潮流介数为3.319,由表2可见,通过线路26-29的输电路径只有两条,分别为{38,29,26}和{38,29,26,27},对应的输电路径有功功率标幺值分别为0.8961和1.0123。由于通过线路26-29的输电路径较少,其输电路径有功功率较小,因此线路26-29有功潮流介数较小。
5)确定电力系统中UPFC的装设位置
第(1)-4)步完成后,先设定电力系统中装设UPFC的线路数量的范围为1~M(本实施例中M取为10);然后将第(1)-4)步得到的电力系统等效无损网络中各线路的有功潮流介数从大到小排序,并确定排在范围1~M内的前m(本实施例中m取为5)条线路的首端装设UPFC,以保证电力系统运行的经济性和可行性。M和m的取值根据电力系统的运行成本和实际操作情况确定。
以图2的电力系统为例,说明图中线路的有功潮流介数排序结果,如表4所示。
表4线路的有功潮流介数排序结果
排序 | 首-末端点 | 有功潮流介数 |
1 | 6-31 | 19.805 |
2 | 5-6 | 19.568 |
3 | 10-32 | 19.328 |
4 | 19-33 | 15.060 |
5 | 22-35 | 14.606 |
32 | 26-29 | 3.319 |
(2)确定装设有UPFC的电力系统各节点的有功潮流介数阈值
1)计算装设有UPFC线路的等效节点有功功率注入增量和无功功率注入增量
第(1)-5)步完成后,利用装设有UPFC的线路i-j两端节点的电压计算其节点等效有功功率注入增量和节点等效无功功率注入增量,具体计算公式为:
式中:ΔPi和ΔQi分别为装设有UPFC的线路i-j在节点i的等效有功功率注入增量和等效无功功率注入增量;ΔPj和ΔQj分别为装设有UPFC的线路i-j在节点j的等效有功功率注入增量和等效无功功率注入增量;gij和bij分别为线路i-j的电导和电纳;Bc为线路i-j的对地电纳;θij为节点i和节点j的之间的相角差,θij=θi-θj;Ui和Uj分别为节点i和节点j的电压幅值;kse为UPFC的幅值控制参数;为UPFC的相角控制参数;Qsh为UPFC的无功控制参数。
2)求解装设UPFC后的电力系统的交流潮流
第(2)-1)步完成后,首先,将按第(2)-1)步得到的装有UPFC的线路i-j两端节点的节点等效有功功率注入增量ΔPi、ΔPj和节点等效无功功率注入增量ΔQi、ΔQj叠加到相应节点的节点功率方程中。具体计算公式如下:
式中:Pgi和Qgi分别为节点i的发电机有功出力和无功出力;Pdi和Qdi分别为节点i的有功负荷和无功负荷;Gij和Bij分别为节点i和节点j的之间的互电导和互电纳;j∈i表示∑号后的标号为j的节点必须直接和节点i相连,并包括j=i的情况;Lu为所有与节点i相连且装有UPFC的线路集合;其他参数的含义与公式(2)相同。
然后,用常规的交流潮流计算方法,求解电力系统实际运行的交流潮流,确定电力系统中各发电机输出有功功率、各节点的电压、各条线路的交流潮流大小及其方向、功率损耗及充电功率。
以图2的电力系统为例,说明图中部分节点的发电机输出有功功率,如表5所示。
表5部分节点的发电机输出有功功率
节点编号 | 节点类型 | 输出有功功率 |
31 | 平衡节点 | 5.1387 |
32 | PV节点 | 6.50 |
33 | PV节点 | 6.32 |
3)计算装设有UPFC的电力系统等效无损网络中各节点的有功潮流介数
第(2)-2)步完成后,根据本申请人申请的申请号为201310213277.8的“基于有功潮流介数的电力系统关键节点辨识方法”的第(2)步和第(3)-1)步,先将所述的电力系统等效成无损网络,然后确定所述的等效无损网络中的输电路径,再计算所述的等效无损网络中各节点的有功功率分配比例因子,然后确定所述的等效无损网络中各输电路径的有功功率,最后计算所述的等效无损网络各节点的有功潮流介数。
4)确定装设有UPFC的电力系统各节点的有功潮流介数阈值
第(2)-3)步完成后,确定电力系统中各节点的有功潮流介数运行阈值和有功潮流介数极限阈值,计算公式为:
式中:Sn为节点n的有功潮流介数运行阈值;Γn为节点n的有功潮流介数极限阈值;Bf(n)为故障前节点n的有功潮流介数。
以图2的电力系统为例,说明图中部分节点的运行阈值和极限阈值,如表6所示。
表6部分节点的运行阈值和极限阈值
节点编号 | 节点运行阈值 | 节点极限阈值 |
32 | 75.855 | 151.71 |
5 | 69.03 | 138.06 |
31 | 57.645 | 115.29 |
8 | 51.765 | 103.53 |
16 | 49.44 | 98.88 |
(3)对装设有UPFC的电力系统进行连锁故障模拟
第(2)步完成后,先通过移除所述的电力系统中的任意节点和与其相连的所有线路来模拟实际电力系统中发电厂或变电站发生故障。再通过计算所述的故障电力系统各节点的有功潮流介数来判断系统是否发生连锁故障,包括形成所述的故障电力系统的子系统及个数、确定所述的故障电力系统的各子系统等效无损网络的节点有功潮流介数、运用得到的节点有功潮流介数判断系统是否发生连锁故障,来模拟连锁故障的传播和终止。具体步骤如下:
1)运用广度优先搜索方法形成所述的故障电力系统的子系统及个数
首先从电力系统中的任一节点n出发,搜索其相邻的节点构成集合A0,将节点n和集合A0中的元素放入集合B;再搜索与A0中各节点相邻且不属于集合B的节点构成集合A1,将集合A1中的元素放入集合B中;继续搜索与A1中节点相邻且不属于集合B的节点集合A2;如此搜索下去,直到不能再搜索到新的节点集合为止;本次搜索到的所有节点集合B形成子系统I1。然后在本次未搜索到的节点中,重复上述搜索过程,依次形成子系统I2,I3,……,IK,直到不存在未被搜索到的节点为止。K为形成的子系统个数。
以图2的电力系统中节点16遭受攻击后形成的子系统为例,运用广度优先搜索方法形成的子系统为I1={16},I2={19,20,33,34},I3={1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,17,18,25,26,27,28,29,30,31,32,37,38,39},I4={21,22,23,24,35,36},子系统个数K=4。
2)确定所述的故障电力系统的各子系统等效无损网络的节点有功潮流介数
①第(3)-1)步完成后,先将第(3)-1)步得到的子系统I1,I2,……,IK中装设有UPFC装置的子系统构成集合Du,没有装设UPFC装置的子系统构成集合D;
以图2的电力系统为例,说明图中16节点移除后子系统分为D和Du两类,如表7所示。
表7子系统分类
②第(3)-2)-①步完成后,对集合Du和集合D中有发电机节点但没有平衡节点的子系统,将该子系统中发电机装机容量最大的节点作为平衡节点;
③第(3)-2)-②步完成后,对所述集合Du的各子系统,为减小故障电力系统的各子系统中线路潮流的越限值,用计及UPFC的交流最优潮流方法确定各子系统的各发电机输出有功功率、各节点的电压、各条线路的交流潮流大小及其方向、功率损耗及充电功率。具体步骤如下:
首先,对第(3)-2)-①步中所述的集合Du中的各子系统,按照公式(2)计算装设有UPFC的线路i-j两端节点的节点等效有功功率注入增量和节点等效无功功率注入增量。
然后,运用常规内点法求解计及UPFC的交流最优潮流模型。该模型的目标函数为所有线路潮流越限值之和最小,约束条件包括潮流方程,发电机有功和无功出力、节点电压模值、UPFC幅值及相角控制参数和无功控制参数可调范围限制等约束。具体公式如下:
Pgi,min≤Pgi≤Pgi,max (8)
Qgi,min≤Qgi≤Qgi,max (9)
Ui,min≤Ui≤Ui,max (10)
0≤kse≤kse,max (11)
-Qsh,max≤Qsh≤Qsh,max (13)
式中:T为所有线路潮流越限值之和;N为子系统中所有节点的集合;Pij和Qij分别为流过线路i-j首端的有功功率和无功功率,Tij,max为线路i-j的额定视在载流容量;Lu为所有与节点i相连且装有UPFC的线路集合;Pgi,min、Qgi,min分别为节点i的发电机有功出力和无功出力下限;Pgi,max和Qgi,max分别为节点i的发电机有功出力和无功出力上限;Ui,min和Ui,max分别为节点i的电压幅值下限和上限;kse和kse,max分别为UPFC的幅值控制参数及最大值;为UPFC的相角控制参数;Qsh和Qshmax分别为UPFC的无功控制参数及最大值;其他参数的含义与公式(3)相同。
④第(3)-2)-③步完成后,对所述的集合D的各子系统用常规的交流潮流计算方法,求解系统实际运行的交流潮流,确定各子系统中各发电机输出有功功率、各节点的电压、各条线路的交流潮流大小及其方向、功率损耗及充电功率。
⑤第(3)-2)-④步完成后,转到第(2)-3)步,计算所述的故障电力系统的各子系统等效无损网络的节点有功潮流介数。
以图2的电力系统为例,说明图中16节点移除后,集合D中子系统I1和集合Du中子系统I4中各节点的有功潮流介数。因子系统I1中只包含节点16,其有功潮流介数为0,子系统I4各节点的有功潮流如表8所示。
表8子系统I4中各节点的有功潮流介数
节点编号 | 有功潮流介数 |
21 | 7.324012 |
22 | 16.001 |
23 | 13.33378 |
24 | 8.436428 |
35 | 16.001 |
36 | 4.6616 |
3)模拟连锁故障的传播和终止
第(3)-2)-⑤步完成后,将第(3)-2)-⑤步得到的故障电力系统的各子系统中各节点的有功潮流介数与其有功潮流介数阈值进行比较,来模拟连锁故障的传播和终止。
①对于第(3)-2)-①步所述的集合Du中的各子系统,模拟连锁故障的步骤如下:
(a)对每一条节点设置一个标志,并置为0;
(b)第(3)-3)-①-(a)步完成后,对该子系统中的每个节点进行判断,当该节点的有功潮流介数处于其有功潮流介数运行阈值和有功潮流介数极限阈值之间时,则将标志置为1;当该节点的有功潮流介数大于其有功潮流介数极限阈值时,则移除该节点和与其相连的所有线路,将标志置为2;直到判断完该子系统中所有节点为止;
(c)第(3)-3)-①-(b)步完成后,当所有节点的标志均为0时,则该子系统的连锁故障模拟过程结束;当所有节点的标志均不为2,且有1时,则重复第(3)-2)步至(3)-3)步;当所有节点的标志中有一个为2时,则重复第(3)-1)步至(3)-3)步;直到该子系统的各节点的有功潮流介数均小于其有功潮流介数运行阈值为止。
②对于第(3)-2)-①步所述的集合D中的各子系统,模拟连锁故障的步骤如下:
(a)当该子系统中只包含一个节点,或该子系统各节点的有功潮流介数均未超过其有功潮流介数极限阈值时,则该子系统的连锁故障模拟过程结束,否则转到第(3)-3)-②-(b)步;
(b)第(3)-3)-②-(a)步完成后,当该子系统中有节点的有功潮流介数大于其有功潮流介数极限阈值的节点时,则移除该节点和与其相连的所有线路;
(c)第(3)-3)-②-(b)步完成后,首先用第(3)-1)步的方法确定该子系统的新形成的各子系统;其次用第(3)-2)-②步确定新的各子系统中的平衡节点,再用第(3)-2)-④确定新的各子系统中各发电机输出有功功率、各节点的电压、各条线路的交流潮流大小及其方向、功率损耗及充电功率;然后用第(3)-2)-⑤步的方法计算新的各子系统等效无损网络中各节点的有功潮流介数;最后转至第(3)-3)-②-(a)步,直到各子系统等效无损网络中各节点的有功潮流介数均未超过其有功潮流介数极限阈值为止。
(4)计算反映连锁故障严重程度的衡量指标
第(3)-3)步完成后,对第(3)步所述的故障电力系统的各子系统,先根据电力系统在发生故障后应该满足运行的安全约束条件,当其不能满足运行的安全约束条件时,应该通过调整发电机出力和统一潮流控制器的控制参数使故障电力系统的运行尽量满足运行的安全约束条件,当调整控制后仍不能使故障电力系统运行在安全约束范围内时,应该使其失负荷最小的原则,计算所述的故障电力系统中各子系统的最小有功切负荷量。再通过计算满足安全约束条件下的最小失负荷百分比来反映连锁故障严重程度。具体步骤如下:
1)计算所述的故障电力系统中各子系统的最小有功切负荷量
运用常规的内点法求解计及UPFC的交流最优负荷削减模型来计算所述的故障电力系统的每一个子系统满足运行安全约束条件的最小有功切负荷量。该模型的目标函数为每一个子系统削减的负荷总量最小,约束条件包括潮流方程,发电机有功和无功出力、节点电压模值、UPFC幅值及相角控制参数和无功控制参数可调范围限制、削减负荷前后功率因数保持恒定、节点削减负荷量和线路潮流等约束,具体计算公式为:
Pgi,min≤Pgi≤Pgi,max (17)
Qgi,min≤Qgi≤Qgi,max (18)
Ui,min≤Ui≤Ui,max (19)
0≤kse≤kse,max (20)
-Qsh,max≤Qsh≤Qsh,max (22)
0≤Pci≤Pdi (24)
0≤Qci≤Qdi (25)
式中,C为子系统中各节点的有功切负荷总量;Pci为节点i的有功负荷削减量;Qci为节点i的无功负荷削减量;L为子系统中所有节点的集合;其他参数的含义与公式(5)到公式(13)相同。
2)计算反映连锁故障严重程度的衡量指标
第(4)-1)步完成后,计算衡量连锁故障严重程度的最小失负荷百分比指标,其具体计算公式为:
式中:Lmin为最小失负荷百分比;S为故障后系统的子系统个数;Cs为子系统s削减有功负荷的总量;L0为电力系统故障前的有功负荷总和。
以图2的电力系统为例,说明系统部分节点遭受攻击后的最小失负荷百分比,如表9所示。
表9 系统部分节点遭受攻击后的最小失负荷百分比
节点编号 | 最小失负荷百分比 |
3 | 3.220 |
6 | 6.800 |
15 | 3.200 |
16 | 19.040 |
24 | 3.086 |
36 | 0.000 |
根据不同节点遭受攻击后电力系统的最小失负荷百分比,就得到连锁故障进行风险评估的结果,进而为制定降低连锁故障风险的预防策略提供了科学依据,能有效地避免电力系统大面积停电的发生。
从上述实施例的结果知:
1.利用线路的有功潮流介数确定UPFC的优化装设位置,具有可操作性,更加符合电力系统的实际运行情况。以表3中线路6-31和线路26-29为例进行说明。从图2电力系统接线原理图知,线路6-31的主要作用是将发电机节点31的出力输送到网络中;由表2的结果知,除发电机节点31到该负荷节点的输电路径外,从发电机节点31出发的其余10条输电路径都经过线路6-31,该线路中流过的所有输电路径的有功功率总和为663.16MW,占发电机节点31的发电机有功出力的97.4%,占该系统总发电机出力的10.54%,因而线路6-31对潮流传播的贡献很大;由表4知线路6-31的有功潮流介数为19.805,其排序为1。由表2的结果知,经过线路26-29的输电路径只有2条,该线路中流过的所有输电路径的有功功率总和为190.84MW,占发电机节点38有功出力的22.99%,占系统总发电机出力的3.03%,因而线路26-29对潮流传播的贡献较小,由表4知线路26-29的有功潮流介数为3.319,排序为32。
由此可见,线路的有功潮流介数量化了“发电-负荷”节点对之间的线路参数对潮流控制的影响,线路的有功潮流介数越大,表明该线路对潮流传播的贡献越大,在该线路装设统一潮流控制器(UPFC)后对电力系统潮流的调节能力也越大。
2.与现有电气介数模拟法按照固定比例调整相关线路参数相比,本发明在通过UPFC调节线路参数和潮流时,如第(3)-2)-③步中的公式(5)~(13)所示,考虑到了保证电力系统安全运行的约束条件,包括潮流方程,发电机有功和无功出力、节点电压模值、UPFC幅值及相角控制参数和无功控制参数可调范围限制等;考虑了电力系统安全运行要满足的所有线路潮流越限值之和最小的目标函数,更加符合电力系统的运行实际,因而本发明的连锁故障模拟方法能更加真实地反映电力系统连锁故障的传播过程,能有效地避免电力系统大面积停电的发生。
Claims (1)
1.一种基于统一潮流控制器的电力系统连锁故障模拟方法,利用计算机,通过程序进行计算,在电力系统中装设统一潮流控制器后对其进行连锁故障模拟其特征在于所述方法的具体步骤如下:
(1)确定电力系统中UPFC的装设位置
1)输入基本参数
首先输入电力系统的基本参数和统一潮流控制器的基本参数,所述的电力系统的基本参数包括节点编号、节点类型、节点对应电压等级、各节点的有功负荷(PL)和无功负荷(QL)、与发电机连接的节点的编号、各发电机输出的有功功率(PG)和无功功率(QG)、各线路首末端节点编号、线路电阻(R)、线路电抗(X)和线路电纳(B)、线路的额定电压(UB)、基准功率(SB)、节点的运行系数(α)、节点的极限系数(β);所述的UPFC控制器的基本参数包括在系统中安装UPFC的台数(m)、每台UPFC的幅值控制参数(kse)、相角控制参数和无功控制参数(Qsh);
2)求解电力系统实际运行的交流潮流
第(1)-1)步完成后,根据第(1)-1)步输入的电力系统的基本参数,用常规的交流潮流计算方法,求解电力系统实际运行的交流潮流,确定电力系统中各发电机的输出有功功率、各节点的电压、各条线路的交流潮流的大小及其方向、功率损耗及充电功率;
3)确定电力系统等效无损网络中各输电路径的有功功率
第(1)-2)步完成后,利用本申请人申请的申请号为201310213277.8的“基于有功潮流介数的电力系统关键节点辨识方法”的第(2)步,确定所述的等效无损网络中各输电路径Bm的有功功率,即:首先将电力系统等效成无损网络,然后确定所述的等效无损网络中的输电路径,再计算所述的等效无损网络中各节点的有功功率分配比例因子,最后确定所述的等效无损网络中各输电路径的有功功率;
4)计算电力系统等效无损网络中各线路的有功潮流介数
第(1)-3)步完成后,将第(1)-3)步得到的通过线路l的所有输电路径的有功功率加权求和来计算该线路的有功潮流介数,用其线路的有功潮流介数来反映调整该线路参数对潮流影响的程度,即用线路的有功潮流介数来量化“发电-负荷”节点对之间的线路参数对潮流控制的影响,由于当通过所述的等效无损网络中各线路的输电路径越多和输电路径的有功功率越大时,调整该线路参数对潮流的影响越大,因此线路的有功潮流介数越大,表明该线路对潮流传播的贡献越大,在该线路装设统一潮流控制器(UPFC)后对电力系统潮流的调节能力也越大,计算线路l的有功潮流介数的具体公式如下:
式中:Bf(l)为线路l的有功潮流介数;G为发电机节点集合;L为负荷节点集合;B为输电路径集合;l、y和z分别为输电路径Bm的线路、发电机节点和负荷节点;Wy为发电机节点y的权重,Wy的值为发电机节点y的有功出力值;Wz为负荷节点z的权重,Wz的值为负荷节点z的有功负荷;Pm(l),y,z为输电路径Bm经过线路l的输电路径有功功率;
5)确定电力系统中UPFC的装设位置
第(1)-4)步完成后,先设定电力系统中装设UPFC的线路数量的范围为1~M;然后将第(1)-4)步得到的电力系统等效无损网络中各线路的有功潮流介数从大到小排序,并确定排在范围1~M内的前m条线路的首端装设UPFC,以保证电力系统运行的经济性和可行性,M和m的取值根据电力系统的运行成本和实际操作情况确定;
(2)确定装设有UPFC的电力系统各节点的有功潮流介数阈值
1)计算装设有UPFC线路的等效节点有功功率注入增量和无功功率注入增量
第(1)-5)步完成后,利用装设有UPFC的线路i-j两端节点的电压计算其节点等效有功功率注入增量和节点等效无功功率注入增量,具体计算公式为:
式中:ΔPi和ΔQi分别为装设有UPFC的线路i-j在节点i的等效有功功率注入增量和等效无功功率注入增量;ΔPj和ΔQj分别为装设有UPFC的线路i-j在节点j的等效有功功率注入增量和等效无功功率注入增量;gij和bij分别为线路i-j的电导和电纳;Bc为线路i-j的对地电纳;θij为节点i和节点j的之间的相角差,θij=θi-θj;Ui和Uj分别为节点i和节点j的电压幅值;kse为UPFC的幅值控制参数;为UPFC的相角控制参数;Qsh为UPFC的无功控制参数;
2)求解装设UPFC后的电力系统的交流潮流
第(2)-1)步完成后,首先,将按第(2)-1)步得到的装有UPFC的线路i-j两端节点的节点等效有功功率注入增量ΔPi、ΔPj和节点等效无功功率注入增量ΔQi、ΔQj叠加到相应节点的节点功率方程中,具体计算公式如下:
式中:Pgi和Qgi分别为节点i的发电机有功出力和无功出力;Pdi和Qdi分别为节点i的有功负荷和无功负荷;Gij和Bij分别为节点i和节点j的之间的互电导和互电纳;j∈i表示∑号后的标号为j的节点必须直接和节点i相连,并包括j=i的情况;Lu为所有与节点i相连且装有UPFC的线路集合;其他参数的含义与公式(2)相同;
然后,用常规的交流潮流计算方法,求解电力系统实际运行的交流潮流,确定电力系统中各发电机输出有功功率、各节点的电压、各条线路的交流潮流大小及其方向、功率损耗及充电功率;
3)计算装设有UPFC的电力系统等效无损网络中各节点的有功潮流介数
第(2)-2)步完成后,根据本申请人申请的申请号为201310213277.8的“基于有功潮流介数的电力系统关键节点辨识方法”的第(2)步和第(3)-1)步,先将所述的电力系统等效成无损网络,然后确定所述的等效无损网络中的输电路径,再计算所述的等效无损网络中各节点的有功功率分配比例因子,然后确定所述的等效无损网络中各输电路径的有功功率,最后计算所述的等效无损网络各节点的有功潮流介数;
4)确定装设有UPFC的电力系统各节点的有功潮流介数阈值
第(2)-3)步完成后,确定电力系统中各节点的有功潮流介数运行阈值和有功潮流介数极限阈值,计算公式为:
式中:Sn为节点n的有功潮流介数运行阈值;Γn为节点n的有功潮流介数极限阈值;Bf(n)为故障前节点n的有功潮流介数;
(3)对装设有UPFC的电力系统进行连锁故障模拟
第(2)步完成后,先通过移除所述的电力系统中的任意节点和与其相连的所有线路来模拟实际电力系统中发电厂或变电站发生故障,再通过计算所述的故障电力系统各节点的有功潮流介数来判断系统是否发生连锁故障,包括形成所述的故障电力系统的子系统及个数、确定所述的故障电力系统的各子系统等效无损网络的节点有功潮流介数、运用得到的节点有功潮流介数判断系统是否发生连锁故障,来模拟连锁故障的传播和终止,具体步骤如下:
1)运用广度优先搜索方法形成所述的故障电力系统的子系统及个数
首先从电力系统中的任一节点n出发,搜索其相邻的节点构成集合A0,将节点n和集合A0中的元素放入集合B;再搜索与A0中各节点相邻且不属于集合B的节点构成集合A1,将集合A1中的元素放入集合B中;继续搜索与A1中节点相邻且不属于集合B的节点集合A2;如此搜索下去,直到不能再搜索到新的节点集合为止;本次搜索到的所有节点集合B形成子系统I1,然后在本次未搜索到的节点中,重复上述搜索过程,依次形成子系统I2,I3,……,IK,直到不存在未被搜索到的节点为止,K为形成的子系统个数;
2)确定所述的故障电力系统的各子系统等效无损网络的节点有功潮流介数
①第(3)-1)步完成后,先将第(3)-1)步得到的子系统I1,I2,……,IK中装设有UPFC装置的子系统构成集合Du,没有装设UPFC装置的子系统构成集合D;
②第(3)-2)-①步完成后,对集合Du和集合D中有发电机节点但没有平衡节点的子系统,将该子系统中发电机装机容量最大的节点作为平衡节点;
③第(3)-2)-②步完成后,对所述集合Du的各子系统,为减小故障电力系统的各子系统中线路潮流的越限值,用计及UPFC的交流最优潮流方法确定各子系统的各发电机输出有功功率、各节点的电压、各条线路的交流潮流大小及其方向、功率损耗及充电功率,具体步骤如下:
首先,对第(3)-2)-①步中所述的集合Du中的各子系统,按照公式(2)计算装设有UPFC的线路i-j两端节点的节点等效有功功率注入增量和节点等效无功功率注入增量;
然后,运用常规内点法求解计及UPFC的交流最优潮流模型,该模型的目标函数为所有线路潮流越限值之和最小,约束条件包括潮流方程,发电机有功和无功出力、节点电压模值、UPFC幅值及相角控制参数和无功控制参数可调范围限制等约束,具体公式如下:
Pgi,min≤Pgi≤Pgi,max (8)
Qgi,min≤Qgi≤Qgi,max (9)
Ui,min≤Ui≤Ui,max (10)
0≤kse≤kse,max (11)
-Qsh,max≤Qsh≤Qsh,max (13)
式中:T为所有线路潮流越限值之和;N为子系统中所有节点的集合;Pij和Qij分别为流过线路i-j首端的有功功率和无功功率,Tij,max为线路i-j的额定视在载流容量;Lu为所有与节点i相连且装有UPFC的线路集合;Pgi,min、Qgi,min分别为节点i的发电机有功出力和无功出力下限;Pgi,max和Qgi,max分别为节点i的发电机有功出力和无功出力上限;Ui,min和Ui,max分别为节点i的电压幅值下限和上限;kse和kse,max分别为UPFC的幅值控制参数及最大值;为UPFC的相角控制参数;Qsh和Qshmax分别为UPFC的无功控制参数及最大值;其他参数的含义与公式(3)相同;
④第(3)-2)-③步完成后,对所述的集合D的各子系统用常规的交流潮流计算方法,求解系统实际运行的交流潮流,确定各子系统中各发电机输出有功功率、各节点的电压、各条线路的交流潮流大小及其方向、功率损耗及充电功率;
⑤第(3)-2)-④步完成后,转到第(2)-3)步,计算所述的故障电力系统的各子系统等效无损网络的节点有功潮流介数;
3)模拟连锁故障的传播和终止
第(3)-2)-⑤步完成后,将第(3)-2)-⑤步得到的故障电力系统的各子系统中各节点的有功潮流介数与其有功潮流介数阈值进行比较,来模拟连锁故障的传播和终止;
①对于第(3)-2)-①步所述的集合Du中的各子系统,模拟连锁故障的步骤如下:
(a)对每一条节点设置一个标志,并置为0;
(b)第(3)-3)-①-(a)步完成后,对该子系统中的每个节点进行判断,当该节点的有功潮流介数处于其有功潮流介数运行阈值和有功潮流介数极限阈值之间时,则将标志置为1;当该节点的有功潮流介数大于其有功潮流介数极限阈值时,则移除该节点和与其相连的所有线路,将标志置为2;直到判断完该子系统中所有节点为止;
(c)第(3)-3)-①-(b)步完成后,当所有节点的标志均为0时,则该子系统的连锁故障模拟过程结束;当所有节点的标志均不为2,且有1时,则重复第(3)-2)步至(3)-3)步;当所有节点的标志中有一个为2时,则重复第(3)-1)步至(3)-3)步;直到该子系统的各节点的有功潮流介数均小于其有功潮流介数运行阈值为止;
②对于第(3)-2)-①步所述的集合D中的各子系统,模拟连锁故障的步骤如下:
(a)当该子系统中只包含一个节点,或该子系统各节点的有功潮流介数均未超过其有功潮流介数极限阈值时,则该子系统的连锁故障模拟过程结束,否则转到第(3)-3)-②-(b)步;
(b)第(3)-3)-②-(a)步完成后,当该子系统中有节点的有功潮流介数大于其有功潮流介数极限阈值的节点时,则移除该节点和与其相连的所有线路;
(c)第(3)-3)-②-(b)步完成后,首先用第(3)-1)步的方法确定该子系统的新形成的各子系统;其次用第(3)-2)-②步确定新的各子系统中的平衡节点,再用第(3)-2)-④确定新的各子系统中各发电机输出有功功率、各节点的电压、各条线路的交流潮流大小及其方向、功率损耗及充电功率;然后用第(3)-2)-⑤步的方法计算新的各子系统等效无损网络中各节点的有功潮流介数;最后转至第(3)-3)-②-(a)步,直到各子系统等效无损网络中各节点的有功潮流介数均未超过其有功潮流介数极限阈值为止;
(4)计算反映连锁故障严重程度的衡量指标
第(3)-3)步完成后,对第(3)步所述的故障电力系统的各子系统,先根据电力系统在发生故障后应该满足运行的安全约束条件,当其不能满足运行的安全约束条件时,应该通过调整发电机出力和统一潮流控制器的控制参数使故障电力系统的运行尽量满足运行的安全约束条件,当调整控制后仍不能使故障电力系统运行在安全约束范围内时,应该使其失负荷最小的原则,计算所述的故障电力系统中各子系统的最小有功切负荷量,再通过计算满足安全约束条件下的最小失负荷百分比来反映连锁故障严重程度,具体步骤如下:
1)计算所述的故障电力系统中各子系统的最小有功切负荷量
运用常规的内点法求解计及UPFC的交流最优负荷削减模型来计算所述的故障电力系统的每一个子系统满足运行安全约束条件的最小有功切负荷量,该模型的目标函数为每一个子系统削减的负荷总量最小,约束条件包括潮流方程,发电机有功和无功出力、节点电压模值、UPFC幅值及相角控制参数和无功控制参数可调范围限制、削减负荷前后功率因数保持恒定、节点削减负荷量和线路潮流等约束,具体计算公式为:
Pgi,min≤Pgi≤Pgi,max (17)
Qgi,min≤Qgi≤Qgi,max (18)
Ui,min≤Ui≤Ui,max (19)
0≤kse≤kse,max (20)
-Qsh,max≤Qsh≤Qsh,max (22)
0≤Pci≤Pdi (24)
0≤Qci≤Qdi (25)
式中,C为子系统中各节点的有功切负荷总量;Pci为节点i的有功负荷削减量;Qci为节点i的无功负荷削减量;L为子系统中所有节点的集合;其他参数的含义与公式(5)到公式(13)相同;
2)计算反映连锁故障严重程度的衡量指标
第(4)-1)步完成后,计算衡量连锁故障严重程度的最小失负荷百分比指标,其具体计算公式为:
式中:Lmin为最小失负荷百分比;S为故障后系统的子系统个数;Cs为子系统s削减有功负荷的总量;L0为电力系统故障前的有功负荷总和;
根据不同节点遭受攻击后电力系统的最小失负荷百分比,就得到连锁故障进行风险评估的结果,进而为制定降低连锁故障风险的预防策略提供了科学依据,能有效地避免电力系统大面积停电的发生。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310282200.6A CN103311926B (zh) | 2013-07-05 | 2013-07-05 | 基于统一潮流控制器的电力系统连锁故障模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310282200.6A CN103311926B (zh) | 2013-07-05 | 2013-07-05 | 基于统一潮流控制器的电力系统连锁故障模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103311926A true CN103311926A (zh) | 2013-09-18 |
CN103311926B CN103311926B (zh) | 2015-04-15 |
Family
ID=49136844
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310282200.6A Active CN103311926B (zh) | 2013-07-05 | 2013-07-05 | 基于统一潮流控制器的电力系统连锁故障模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103311926B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104008455A (zh) * | 2014-05-30 | 2014-08-27 | 重庆大学 | 一种复杂辐射形配电网可靠性评估的前推回代方法 |
CN104112238A (zh) * | 2014-06-30 | 2014-10-22 | 南方电网科学研究院有限责任公司 | 基于概率指标和严重度指标的连锁故障搜索综合评估方法 |
CN104638646A (zh) * | 2015-02-17 | 2015-05-20 | 重庆大学 | 基于复杂网络理论计及校正控制的电网连锁故障模拟方法 |
CN105322535A (zh) * | 2015-10-13 | 2016-02-10 | 河海大学 | 含统一潮流控制器的电力系统两阶段最优潮流计算方法 |
CN105591393A (zh) * | 2016-02-04 | 2016-05-18 | 重庆大学 | 基于支路无功潮流介数的电力系统无功耦合度辨识方法 |
CN106850270A (zh) * | 2016-12-30 | 2017-06-13 | 国家电网公司 | 一种iec61850关键逻辑节点辨识方法 |
CN107703913A (zh) * | 2017-09-15 | 2018-02-16 | 东南大学 | 一种upfc故障诊断方法 |
CN109038640A (zh) * | 2018-08-13 | 2018-12-18 | 深圳供电局有限公司 | 一种基于相量图解法的upfc串联侧换流器定容方法 |
CN109470990A (zh) * | 2018-10-25 | 2019-03-15 | 南京南瑞继保电气有限公司 | 一种适应upfc接入的线路变化量故障方向判断方法和装置 |
CN109615248A (zh) * | 2018-12-14 | 2019-04-12 | 浙江大学 | 一种考虑天然气系统影响的电力系统连锁故障估计方法 |
CN113723821A (zh) * | 2021-08-31 | 2021-11-30 | 广东电网有限责任公司 | 一种基于潮流介数的电网故障预警方法及装置 |
CN115051416A (zh) * | 2022-08-16 | 2022-09-13 | 阿里巴巴(中国)有限公司 | 数据处理方法、发电方法、装置和云设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008252999A (ja) * | 2007-03-29 | 2008-10-16 | Chubu Electric Power Co Inc | 電力系統の線路損失低減装置、電力系統及び電力系統構築方法 |
CN103151774A (zh) * | 2013-01-30 | 2013-06-12 | 武汉大学 | 一种小世界电网连锁故障抑制方法 |
-
2013
- 2013-07-05 CN CN201310282200.6A patent/CN103311926B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008252999A (ja) * | 2007-03-29 | 2008-10-16 | Chubu Electric Power Co Inc | 電力系統の線路損失低減装置、電力系統及び電力系統構築方法 |
CN103151774A (zh) * | 2013-01-30 | 2013-06-12 | 武汉大学 | 一种小世界电网连锁故障抑制方法 |
Non-Patent Citations (1)
Title |
---|
王韶等: "基于最优负荷削减的小世界电网连锁故障模拟", 《电网技术》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104008455A (zh) * | 2014-05-30 | 2014-08-27 | 重庆大学 | 一种复杂辐射形配电网可靠性评估的前推回代方法 |
CN104008455B (zh) * | 2014-05-30 | 2017-04-12 | 重庆大学 | 一种复杂辐射形配电网可靠性评估的前推回代方法 |
CN104112238A (zh) * | 2014-06-30 | 2014-10-22 | 南方电网科学研究院有限责任公司 | 基于概率指标和严重度指标的连锁故障搜索综合评估方法 |
CN104112238B (zh) * | 2014-06-30 | 2017-08-25 | 南方电网科学研究院有限责任公司 | 基于概率指标和严重度指标的连锁故障搜索综合评估方法 |
CN104638646A (zh) * | 2015-02-17 | 2015-05-20 | 重庆大学 | 基于复杂网络理论计及校正控制的电网连锁故障模拟方法 |
CN105322535B (zh) * | 2015-10-13 | 2017-07-11 | 河海大学 | 含统一潮流控制器的电力系统两阶段最优潮流计算方法 |
CN105322535A (zh) * | 2015-10-13 | 2016-02-10 | 河海大学 | 含统一潮流控制器的电力系统两阶段最优潮流计算方法 |
CN105591393B (zh) * | 2016-02-04 | 2018-06-22 | 重庆大学 | 基于支路无功潮流介数的电力系统无功耦合度辨识方法 |
CN105591393A (zh) * | 2016-02-04 | 2016-05-18 | 重庆大学 | 基于支路无功潮流介数的电力系统无功耦合度辨识方法 |
CN106850270A (zh) * | 2016-12-30 | 2017-06-13 | 国家电网公司 | 一种iec61850关键逻辑节点辨识方法 |
CN106850270B (zh) * | 2016-12-30 | 2020-05-05 | 国家电网公司 | 一种iec61850关键逻辑节点辨识方法 |
CN107703913B (zh) * | 2017-09-15 | 2019-09-06 | 东南大学 | 一种upfc故障诊断方法 |
CN107703913A (zh) * | 2017-09-15 | 2018-02-16 | 东南大学 | 一种upfc故障诊断方法 |
CN109038640A (zh) * | 2018-08-13 | 2018-12-18 | 深圳供电局有限公司 | 一种基于相量图解法的upfc串联侧换流器定容方法 |
CN109038640B (zh) * | 2018-08-13 | 2021-12-03 | 深圳供电局有限公司 | 一种基于相量图解法的upfc串联侧换流器定容方法 |
CN109470990A (zh) * | 2018-10-25 | 2019-03-15 | 南京南瑞继保电气有限公司 | 一种适应upfc接入的线路变化量故障方向判断方法和装置 |
CN109470990B (zh) * | 2018-10-25 | 2020-12-18 | 南京南瑞继保电气有限公司 | 一种适应upfc接入的线路变化量故障方向判断方法和装置 |
CN109615248A (zh) * | 2018-12-14 | 2019-04-12 | 浙江大学 | 一种考虑天然气系统影响的电力系统连锁故障估计方法 |
CN109615248B (zh) * | 2018-12-14 | 2023-02-10 | 浙江大学 | 一种考虑天然气系统影响的电力系统连锁故障估计方法 |
CN113723821A (zh) * | 2021-08-31 | 2021-11-30 | 广东电网有限责任公司 | 一种基于潮流介数的电网故障预警方法及装置 |
CN113723821B (zh) * | 2021-08-31 | 2024-04-26 | 广东电网有限责任公司 | 一种基于潮流介数的电网故障预警方法及装置 |
CN115051416A (zh) * | 2022-08-16 | 2022-09-13 | 阿里巴巴(中国)有限公司 | 数据处理方法、发电方法、装置和云设备 |
Also Published As
Publication number | Publication date |
---|---|
CN103311926B (zh) | 2015-04-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103311926A (zh) | 基于统一潮流控制器的电力系统连锁故障模拟方法 | |
CN104638646B (zh) | 基于复杂网络理论计及校正控制的电网连锁故障模拟方法 | |
CN102255307B (zh) | 一种配电网可靠性评估的分层等值方法 | |
CN114665498A (zh) | 一种计及新能源影响的主动配电网脆弱节点辨识方法 | |
CN101592700B (zh) | 基于事故链的大电网连锁故障分析方法 | |
CN108899904A (zh) | 一种交直流大电网连锁故障快速搜索方法及系统 | |
CN103296685A (zh) | 一种svc补偿策略最优化方法 | |
CN103310390A (zh) | 一种电网安全综合评价方法 | |
CN103795058A (zh) | 电力系统的静态电压稳定裕度分析及系统故障排序方法 | |
CN104538953B (zh) | 一种基于概率潮流控制的tcsc优化配置方法 | |
CN104037790B (zh) | 一种基于时序蒙特卡罗模拟的新能源接纳能力评估方法 | |
CN104716670B (zh) | 光伏并网下基于网络安全约束的机组组合方法 | |
CN107069705A (zh) | 一种电网元件连锁故障模拟方法 | |
CN104901306A (zh) | 一种考虑连锁故障的电网运行安全裕度计算方法 | |
CN104485660A (zh) | 一种基于有向关系图的配电网可靠性评估方法 | |
CN107622360A (zh) | 一种综合考虑主客观因素的关键线路识别方法 | |
CN104505821A (zh) | 一种控制短路电流水平的电网运行方式优化方法 | |
CN106655292A (zh) | 一种特高压等级交直流混合电网稳定性分析方法 | |
CN105305488A (zh) | 一种考虑新能源并网对输电网利用率影响的评价方法 | |
CN105375453A (zh) | 一种基于连锁故障破坏程度指标的紧急控制方法 | |
CN109921459B (zh) | 一种光伏高占比送端电网直流闭锁后的有功频率控制方法 | |
CN103560497A (zh) | 一种基于电网拓扑调整的短路电流限制方法 | |
Liu et al. | A comprehensive decision-making method for wind power integration projects based on improved fuzzy AHP | |
CN107067126A (zh) | 一种基于潮流转移比的热稳定关键输电通道识别方法 | |
CN105552880A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |