电力系统暂态安全稳定输电断面关联度指标快速计算方法
技术领域
本发明属于电力系统安全稳定分析技术领域,更准确地说,本发明涉及适用于对电力系统当前运行方式下各个输电断面暂态安全稳定关联程度的计算。
背景技术
电力系统的运行方式决定了电力系统的安全稳定运行水平,为了便于监视电力系统的安全稳定运行水平,调度运行人员通过监视输电断面的功率变化来掌握电力系统安全稳定运行水平的变化。实际上,大电网中各个输电断面的安全稳定水平可能不仅与本输电断面的输送功率密切相关,而且受其它输电断面输送功率大小的影响,即各个输电断面之间是关联的,关联的程度与电力系统的运行方式紧密相关。
暂态安全稳定包括暂态功角稳定(简称为TAS)、暂态电压安全稳定和暂态频率安全3个方面,其中暂态电压安全稳定包括暂态电压跌落安全(简称为TVDS)和暂态电压稳定(简称为TVS)两个方面,暂态频率安全包括暂态频率跌落安全(简称为TFDS)和暂态频率上升安全(简称为TFUS)两个方面。本发明中的输电断面暂态安全稳定关联度指标是指一个输电断面暂态安全稳定水平的变化对另一个输电断面暂态安全稳定水平变化的影响程度,采用两个输电断面之间的TAS关联度、TVDS关联度、TVS关联度、TFDS关联度和TFUS关联度共5项指标来反映。若两个输电断面的暂态安全稳定的5项关联度指标都大于0,表示电力系统运行状态变化对两个输电断面暂态安全稳定水平影响的方向是一致的,且指标值越大,表示两者之间的关联程度越高,5项关联度指标都等于1,表示电力系统运行状态变化对两个输电断面暂态安全稳定水平的影响程度相同;若两个输电断面的暂态安全稳定的5项关联度指标都小于0,表示电 力系统运行状态变化对两个输电断面暂态安全稳定水平影响的方向是相反的,且指标值负得越多,表示两者之间的关联程度越高,5项关联度指标都等于-1,表示电力系统运行状态变化对两个输电断面暂态安全稳定水平的影响程度相同、方向相反;若两个输电断面的暂态安全稳定5项关联度指标都等于0,表示一个输电断面暂态安全稳定水平与另一个输电断面暂态安全稳定水平之间的没有关联。若两个输电断面的暂态安全稳定的5项关联度指标中有正有负,则表示电力系统运行状态变化对两个输电断面中关联度指标为正的暂态安全稳定子问题的影响方向是一致的,对两个输电断面中关联度指标为负的暂态安全稳定子问题的影响方向是相反的。
专利申请“基于并行计算模式的关联输电断面极限功率计算方法”(受理号:201110111966.9)提出考虑多种安全稳定约束的关联输电断面中单个或多个断面输送功率最大值,以及各输电断面极限功率相互影响因子的计算方法。该专利提出的各输电断面极限功率相互影响因子计算方法建立在已经得到输电断面极限功率的基础上,计算速度还比较慢,对于大电网其计算时间大约需要10分钟左右。此外,该专利申请给出各输电断面极限功率之间的相互影响因子是综合性指标,没有针对暂态安全稳定中的5类安全稳定问题,给出各输电断面之间的关联度指标。
本发明提出的快速计算方法建立在一次安全稳定量化评估的基础上,针对各个输电断面关键故障集中各个故障的TAS裕度及模式、TVDS裕度及模式、TVS裕度及模式、TFDS裕度及模式和TFUS安全裕度及模式之间的关联关系进行分析,只需增加极小的计算量即可得到各个输电断面之间的暂态安全稳定关联度,包括TAS关联度、TVDS关联度、TVS关联度、TFDS关联度和TFUS 关联度共5项指标。对于大电网采用基于大规模集群计算模式进行暂态安全稳定量化评估的计算时间只需要1分钟左右,本发明可为调度运行人员快速掌握当前运行方式下各个输电断面的安全稳定水平之间的交互影响因素及程度提供技术手段。
发明内容
本发明的目的是:针对电力系统的当前运行方式,快速计算出各个输电断面之间的暂态安全稳定关联度,包括TAS关联度、TVDS关联度、TVS关联度、TFDS关联度和TFUS关联度共5项指标,为调度运行人员快速掌握当前运行方式下各个输电断面的暂态安全稳定水平之间的交互影响因素及程度提供技术手段。
具体地说,本发明是采取以下的技术方案来实现的,包括下列步骤:
1)基于电力系统的单个运行方式,对各个输电断面的考核故障集进行暂态安全稳定量化评估,对于不同考核故障集中的相同故障只进行一次暂态安全稳定量化评估,得到该运行方式下所有输电断面的考核故障集中每个故障的TAS裕度及发电机的TAS参与因子、TVDS裕度及监视节点的TVDS参与因子、TVS裕度及监视负荷的TVS参与因子、TFDS裕度及监视节点的TFDS参与因子和监视发电机的TFDS参与因子、TFUS裕度及监视节点的TFUS参与因子和监视发电机的TFUS参与因子;
2)分别针对每个输电断面,从各输电断面的考核故障集中分别筛选出TAS裕度、TVDS裕度、TVS裕度、TFDS裕度和TFUS裕度小于相应的设定的裕度门槛值的故障,即筛选出相对较低的故障,构成每个输电断面的TAS关键故障集、TVDS关键故障集、TVS关键故障集、TFDS关键故障集和TFUS关键故障集;
3)分别针对每个输电断面,根据输电断面的TAS关键故障集中各个故障的TAS裕度及发电机的TAS参与因子,计算出发电机在输电断面TAS中的综合参 与因子;根据输电断面的TVDS关键故障集中各个故障的TVDS裕度及TVDS监视节点的TVDS参与因子,计算出TVDS监视节点在输电断面TVDS中的综合参与因子;根据输电断面的TVS关键故障集中各个故障的TVS裕度及TVS监视负荷的TVS参与因子,计算出TVS监视负荷在输电断面TVS中的综合参与因子;根据输电断面的TFDS关键故障集中各个故障的TFDS裕度及TFDS监视节点的TFDS参与因子和TFDS监视发电机的TFDS参与因子,计算出TFDS监视节点和TFDS监视发电机在输电断面TFDS中的综合参与因子;根据输电断面的TFUS关键故障集中各个故障的TFUS裕度及TFUS监视节点的TFUS参与因子和TFUS监视发电机的TFUS参与因子,计算出TFUS监视节点和TFUS监视发电机在输电断面TFUS中的综合参与因子;
4)针对所有输电断面中每两个输电断面,分别基于两个输电断面中发电机的TAS综合参与因子、TVDS监视节点的TVDS综合参与因子、TVS监视负荷的TVS综合参与因子、TFDS监视节点和TFDS监视发电机的TFDS综合参与因子、TFUS监视节点和TFUS监视发电机的TFUS综合参与因子,计算出两个输电断面之间的TAS关联度指标、TVDS关联度指标、TVS关联度指标、TFDS关联度指标和TFUS关联度指标。
本发明的进一步特征在于:所述步骤2)中的设定的裕度门槛值,是指:将[1-(1-ηa.min)k]作为从输电断面考核故障集中筛选出暂态功角稳定关键故障集的裕度门槛值,其中ηa.min为相应输电断面考核故障集中暂态功角稳定裕度的最小值;将[1-(1-ηvd.min)k]作为从输电断面考核故障集中筛选出暂态电压跌落安全关键故障集的裕度门槛值,其中ηvd.min为相应输电断面考核故障集中暂态电压跌落安全裕度的最小值;将[1-(1-ηvs.min)k]作为从输电断面考核故障集中筛选出暂态电压稳定关键故障集的裕度门槛值,其中ηvs.min为相应输电断面考核故障集中暂态电压稳定裕度的最小值;将[1-(1-ηfd.min)k]作为从输电断面考核故障集 中筛选出暂态频率跌落安全关键故障集的裕度门槛值,其中ηfd.min为相应输电断面考核故障集中暂态频率跌落安全裕度的最小值;将[1-(1-ηfu.min)k]作为从输电断面考核故障集中筛选出暂态频率上升安全关键故障集的裕度门槛值,其中ηfu.min为相应输电断面考核故障集中暂态频率上升安全裕度的最小值;以上k为设定参数,0<k<1。
本发明的进一步特征在于:在步骤3)中,采用公式(1)计算各台发电机在单个输电断面暂态功角稳定中的综合参与因子λa,采用公式(2)计算各个暂态电压跌落安全监视节点在单个输电断面暂态电压跌落安全中的综合参与因子λvd,采用公式(3)计算各个暂态电压稳定监视负荷在单个输电断面暂态电压稳定中的综合参与因子λvs,采用公式(4)计算各个暂态频率跌落安全监视节点或暂态频率跌落安全监视发电机在单个输电断面暂态频率跌落安全中的综合参与因子λfd,采用公式(5)计算各个暂态频率上升安全监视节点或暂态频率上升安全监视发电机在单个输电断面暂态频率上升安全中的综合参与因子λfu:
其中,Fa为步骤2)得到的该输电断面暂态功角稳定关键故障集中的故障数,ηa.i为其中第i个故障的暂态功角稳定裕度,λa.i为发电机在其中第i个故障中的暂态功角稳定参与因子;
Fvd为步骤2)得到的该输电断面暂态电压跌落安全关键故障集中的故障数,ηvd.i为其中第i个故障的暂态电压跌落安全裕度,λvd.i为暂态电压跌落安全监视节点在其中第i个故障中的暂态电压跌落安全参与因子;
Fvs为步骤2)得到的该输电断面暂态电压稳定关键故障集中的故障数,ηvs.i为其中第i个故障的暂态电压稳定裕度,λvs.i为暂态电压稳定监视负荷在其中第i个故障中的暂态电压稳定参与因子;
Ffd为步骤2)得到的该输电断面暂态频率跌落安全关键故障集中的故障数,ηfd.i为其中第i个故障的暂态频率跌落安全裕度,λfd.i为暂态频率跌落安全监视节点或暂态频率跌落安全监视发电机在其中第i个故障中的暂态频率跌落安全参与因子;
Ffu为步骤2)得到的该输电断面暂态频率上升安全关键故障集中的故障数,ηfu.i为其中第i个故障的暂态频率上升安全裕度,λfu.i为暂态频率上升安全监视节点或暂态频率上升安全监视发电机在第i个故障中的暂态频率上升安全参与 因子。
本发明的进一步特征在于:在步骤4)中,采用公式(6)计算两个输电断面的暂态功角稳定关联度指标Ra,采用公式(7)计算两个输电断面的暂态电压跌落安全关联度指标Rvd,采用公式(8)计算两个输电断面的暂态电压稳定关联度指标Rvs,采用公式(9)计算两个输电断面的暂态频率跌落安全关联度指标Rfd,采用公式(10)计算两个输电断面的暂态频率上升安全关联度指标Rfu:
其中,N为发电机台数,λa.i.1为其中第i台发电机在步骤3)中得到的与这两个输电断面中第1个输电断面对应的暂态功角稳定综合参与因子,λa.i.2为该 发电机在步骤3)中得到的与这两个输电断面中第2个输电断面对应的暂态功角稳定综合参与因子;
Nvd为暂态电压跌落安全监视节点的个数,λvd.i.1为其中第i个暂态电压跌落安全监视节点在步骤3)中得到的与这两个输电断面中第1个输电断面对应的暂态电压跌落安全综合参与因子,λvd.i.2为该暂态电压跌落安全监视节点在步骤3)中得到的与这两个输电断面中第2个输电断面对应的暂态电压跌落安全综合参与因子;
Nvs为暂态电压稳定监视负荷的个数,λvs.i.1为其中第i个暂态电压稳定监视负荷在步骤3)中得到的与这两个输电断面中第1个输电断面对应的暂态电压稳定综合参与因子,λvs.i.2为该暂态电压稳定监视负荷在步骤3)中得到的与这两个输电断面中第2个输电断面对应的暂态电压稳定综合参与因子;
Nfd为暂态频率跌落安全监视节点和暂态频率跌落安全监视发电机的总数,λfd.i.1为其中第i个暂态频率跌落安全监视节点或暂态频率跌落安全监视发电机在步骤3)中得到的与这两个输电断面中第1个输电断面对应的暂态频率跌落安全综合参与因子,λfd.i.2为该暂态频率跌落安全监视节点或暂态频率跌落安全监视发电机在步骤3)中得到的与这两个输电断面中第2个输电断面对应的暂态频率跌落安全综合参与因子;
Nfu为暂态频率上升安全监视节点和暂态频率上升安全监视发电机的总数,λfu.i.1为其中第i个暂态频率上升安全监视节点或暂态频率上升安全监视发电机在步骤3)中得到的与这两个输电断面中第1个输电断面对应的暂态频率上升安全综合参与因子,λfu.i.2为该暂态频率上升安全监视节点或暂态频率上升安全监视发电机在步骤3)中得到的与这两个输电断面中第2个输电断面对应的暂态频 率上升安全综合参与因子。
本发明的有益效果如下:输电断面的暂态安全稳定水平是由关键故障集决定的,关键故障集中各个故障的TAS裕度及模式、TVDS裕度及模式、TVS裕度及模式、TFDS裕度及模式和TFUS裕度及模式是输电断面的暂态安全稳定水平的量化指标,本发明针对各个输电断面关键故障集中各个故障的TAS裕度及模式、TVDS裕度及模式、TVS裕度及模式、TFDS裕度及模式和TFUS裕度及模式之间的关联关系进行分析,能够准确反映各个输电断面的暂态安全稳定关联程度。由于针对各个输电断面考核故障集的量化评估是当前运行方式暂态安全稳定计算的基础工作,本发明只增加了极小的计算量即可得到各个输电断面之间的暂态安全稳定关联程度,包括TAS关联度指标、TVDS关联度指标、TVS关联度指标、TFDS关联度指标和TFUS关联度指标。对于大电网采用基于大规模集群计算模式进行暂态安全稳定量化评估的计算时间只需要1分钟左右,因此,本发明能够为调度运行人员快速掌握当前运行方式下各个输电断面的暂态安全稳定水平之间的交互影响因素及程度提供技术手段。
附图说明
图1为本发明方法的流程图。
具体实施方式
下面参照附图并结合实例对本发明作进一步详细描述。
图1中步骤1描述的是基于电力系统的单个运行方式,对各个输电断面的考核故障集进行暂态安全稳定量化评估,对于不同考核故障集中的相同故障只进行一次暂态安全稳定量化评估,得到该运行方式下所有输电断面的考核故障集中每个故障的TAS裕度及发电机的TAS参与因子,TVDS裕度及监视节点的TVDS参与因子,TVS裕度及监视负荷的TVS参与因子,TFDS裕度及监视节点的TFDS参与因子和监视发电机的TFDS参与因子,以及TFUS裕度及监视节 点的TFUS参与因子和监视发电机的TFUS参与因子。
关于暂态安全稳定量化评估是基于对考核故障集中的故障进行一次时域仿真得到的。其中,基于扩展等面积准则(EEAC)得到TAS的裕度和模式(摆次、主导机群、余下机群)。关于各台发电机参与因子的计算采用下述方法。
若TAS裕度小于0,对于主导机群发电机,计算受扰轨迹经过动态鞍点(DSP)时主导机群中各台发电机的加速动能,以主导机群中发电机加速动能的最大值作为基准,将主导机群中各台发电机的加速动能与该基准值的比值作为各台发电机的TAS参与因子;对于余下机群发电机,计算受扰轨迹经过DSP时余下机群中各台发电机的减速动能,同样以主导机群中发电机加速动能的最大值作为基准,将余下机群中各台发电机的减速动能与该基准值的比值的相反数作为各台发电机的TAS参与因子。
若TAS裕度大于等于0,对于主导机群发电机,首先要确定受扰轨迹在稳定模式中给出的摆次中主导机群等值发电机的加速动能达到最大值的时刻,以该时刻主导机群中发电机加速动能的最大值作为基准,将主导机群中各台发电机该时刻的加速动能与该基准值的比值作为各台发电机的TAS参与因子;对于余下机群发电机,同样以该时刻主导机群中发电机加速动能的最大值作为基准,将余下机群中各台发电机该时刻的减速动能与该基准值的比值的相反数作为各台发电机的TAS参与因子。
基于暂态电压安全稳定量化分析方法得到所有监视节点的TVDS裕度和所有监视负荷的TVS裕度,并分别以所有监视节点TVDS裕度的最小值ηvd.min、所有监视负荷TVS裕度的最小值ηvs.min作为该故障下电力系统的TVDS裕度和TVS裕度。
针对所有的TVDS监视节点,以(1-ηvd.min)为基准,将(1-ηvd.i)/(1-ηvd.min)作为第i个监视节点的TVDS参与因子,其中,ηvd.i为第i个监视节点的TVDS裕度。针对所有的TVS监视负荷,以(1-ηvs.min)为基准,将(1-ηvs.i)/(1-ηvs.min) 作为第i个监视负荷的TVS参与因子,ηvs.i为第i个监视负荷的TVS裕度。
基于暂态频率安全量化分析方法得到所有监视节点、监视发电机的TFDS裕度和TFUS裕度,并分别以所有监视节点和监视发电机TFDS裕度的最小值ηfd.min、所有监视节点和监视发电机TFUS裕度的最小值ηfu.min作为该故障下电力系统的TFDS裕度和TFUS裕度。
针对所有的TFDS监视节点和监视发电机,以(1-ηfd.min)为基准,将(1-ηfd.i)/(1-ηfd.min)作为第i个监视节点或监视发电机的TFDS参与因子,其中,ηfd.i为第i个监视节点或监视发电机的TFDS裕度。针对所有的TFUS监视节点和监视发电机,以(1-ηfu.min)为基准,将(1-ηfu.i)/(1-ηfu.min)作为第i个监视节点或监视发电机的TFUS参与因子,其中,ηfu.i为第i个监视节点或监视发电机的TFUS裕度。
图1中步骤2描述的是针对每个输电断面,将输电断面考核故障集中TAS裕度小于[1-(1-ηa.min)k]的故障筛选出来构成该输电断面的TAS关键故障集,其中ηa.min为该输电断面考核故障集中TAS裕度的最小值,k为设定参数,0<k<1;将输电断面考核故障集中TVDS裕度小于[1-(1-ηvd.min)k]的故障筛选出来构成该输电断面的TVDS关键故障集,其中ηvd.min为该输电断面考核故障集中TVDS裕度的最小值;将输电断面考核故障集中TVS裕度小于[1-(1-ηvs.min)k]的故障筛选出来构成该输电断面的TVS关键故障集,其中ηvs.min为该输电断面考核故障集中TVS裕度的最小值;将输电断面考核故障集中TFDS裕度小于[1-(1-ηfd.min)k]的故障筛选出来构成该输电断面的TFDS关键故障集,其中ηfd.min为该输电断面的考核故障集中TFDS裕度的最小值;将输电断面考核故障集中TFUS裕度小于[1-(1-ηfu.min)k]的故障筛选出来构成该输电断面的TFUS关键故障集,其中ηfu.min为该输电断面的考核故障集中TFUS裕度的最小值。
图1中步骤3描述的是分别针对每个输电断面,采用公式(1)计算各台发电机在该输电断面TAS中的综合参与因子λa,其中Fa为步骤2)得到的该输电断面TAS关键故障集中的故障数,ηa.i为其中第i个故障的TAS裕度,λa.i为发电机在第i个故障中的TAS参与因子;采用公式(2)计算各个TVDS监视节点在该输电断面TVDS中的综合参与因子λvd,其中Fvd为步骤2)得到的该输电断面TVDS关键故障集中的故障数,ηvd.i为其中第i个故障的TVDS裕度,λvd.i为TVDS监视节点在第i个故障中的TVDS参与因子;采用公式(3)计算各个TVS监视负荷在该输电断面TVS中的综合参与因子λvs,其中Fvs为步骤2)得到的该输电断面TVS关键故障集中的故障数,ηvs.i为其中第i个故障的TVS裕度,λvs.i为TVS监视负荷在第i个故障中的TVS参与因子;采用公式(4)计算各个TFDS监视节点或发电机在该输电断面TFDS中的综合参与因子λfd,其中Ffd为步骤2)得到的该输电断面TFDS关键故障集中的故障数,ηfd.i为其中第i个故障的TFDS裕度,λfd.i为TFDS监视节点或发电机在第i个故障中的TFDS参与因子;采用公式(5)计算各个TFUS监视节点或发电机在该输电断面TFUS中的综合参与因子λfu,其中Ffu为步骤2)得到的该输电断面TFUS关键故障集中的故障数,ηfu.i为其中第i个故障的TFUS裕度,λfu.i为TFUS监视节点或发电机在第i个故障中的TFUS参与因子。
图1中步骤4描述的是针对所有输电断面中每两个输电断面,分别基于两个输电断面中发电机的TAS综合参与因子、监视节点的TVDS综合参与因子、监视负荷的TVS综合参与因子、监视节点和发电机的TFDS综合参与因子、以及监视节点和发电机的TFUS综合参与因子,计算出两个输电断面之间的TAS关联度指标、TVDS关联度指标、TVS关联度指标、TFDS关联度指标和TFUS关联度指标。具体计算方法如下。
采用公式(6)计算两个输电断面的TAS关联度指标Ra,其中,N为发电机台数,λa.i.1为其中第i台发电机在步骤3)中得到的与这两个输电断面中第1个输电断面对应的TAS综合参与因子,λa.i.2为该发电机在步骤3)中得到的与这两个输电断面中第2个输电断面对应的TAS综合参与因子。
采用公式(7)计算两个输电断面的TVDS关联度指标Rvd,其中,Nvd为TVDS监视节点个数,λvd.i.1为其中第i个监视节点在步骤3)中得到的与这两个输电断面中第1个输电断面对应的TVDS综合参与因子,λvd.i.2为该监视节点在步骤3)中得到的与这两个输电断面中第2个输电断面对应的TVDS综合参与因子。
采用公式(8)计算两个输电断面的TVS关联度指标Rvs,其中,Nvs为TVS 监视负荷个数,λvs.i.1为其中第i个监视负荷在步骤3)中得到的与这两个输电断面中第1个输电断面对应的TVS综合参与因子,λvs.i.2为该监视负荷在步骤3)中得到的与这两个输电断面中第2个输电断面对应的TVS综合参与因子。
采用公式(9)计算两个输电断面的TFDS关联度指标Rfd,其中,Nfd为TFDS监视节点和发电机总数,λfd.i.1为其中第i个监视节点或发电机在步骤3)中得到的与这两个输电断面中第1个输电断面对应的TFDS综合参与因子,λfd.i.2为该监视节点或发电机在步骤3)中得到的与这两个输电断面中第2个输电断面对应的TFDS综合参与因子。
采用公式(10)计算两个输电断面的TFUS关联度指标Rfu,其中,Nfu为TFUS监视节点和发电机总数,λfu.i.1为其中第i个监视节点或发电机在步骤3)中得到的与这两个输电断面中第1个输电断面对应的TFUS综合参与因子,λfu.i.2为该监视节点或发电机在步骤3)中得到的与这两个输电断面中第2个输电断面对应的TFUS综合参与因子。
虽然本发明已以较佳实施例公开如上,但实施例并不是用来限定本发明的。在不脱离本发明之精神和范围内,所做的任何等效变化或润饰,同样属于本发明之保护范围。因此本发明的保护范围应当以本申请的权利要求所界定的内容为标准。