CN103178518B - 根据轨迹及轨迹灵敏度定量评价电力系统暂态稳定性方法 - Google Patents

根据轨迹及轨迹灵敏度定量评价电力系统暂态稳定性方法 Download PDF

Info

Publication number
CN103178518B
CN103178518B CN201310065501.3A CN201310065501A CN103178518B CN 103178518 B CN103178518 B CN 103178518B CN 201310065501 A CN201310065501 A CN 201310065501A CN 103178518 B CN103178518 B CN 103178518B
Authority
CN
China
Prior art keywords
sensitivity
clearing time
generator
fault
fault clearing
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
Application number
CN201310065501.3A
Other languages
English (en)
Other versions
CN103178518A (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.)
Northeast Electric Power University
Original Assignee
Northeast Dianli University
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 Northeast Dianli University filed Critical Northeast Dianli University
Priority to CN201310065501.3A priority Critical patent/CN103178518B/zh
Publication of CN103178518A publication Critical patent/CN103178518A/zh
Application granted granted Critical
Publication of CN103178518B publication Critical patent/CN103178518B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明是一种根据轨迹及轨迹灵敏度定量评价电力系统暂态稳定性方法,其特点是:首先推导出基于实测轨迹计算电力系统中功角轨迹相对于故障切除时间灵敏度的计算方法,然后分析该灵敏度与系统临界切除时间的关系,通过大量仿真验证可以从此关系中得出用该灵敏度来作为系统稳定程度评价指标。该方法计算简单,为基于轨迹的电力系统动态分析提供重要的理论支持,具有计算快速、评价精度高、应用价值高、在线应用前景大等优点。

Description

根据轨迹及轨迹灵敏度定量评价电力系统暂态稳定性方法
技术领域
本发明涉及一种根据轨迹及轨迹灵敏度定量评价电力系统暂态稳定性方法。
背景技术
随着国民经济的不断发展,电力系统规模越趋庞大,电力系统面临的运行风险也越来越大,如何保证电力系统安全稳定运行已成为电力企业及科研人员关注的重点,因此定量评估电力系统暂态稳定性对电力系统安全稳定运行具有十分重要意义。
现有的评价系统稳定程度的方法有直接法和数值仿真法,其计算结果均受系统模型参数的准确度影响较大,且需计算能量函数和临界能量,而能量函数的构造和临界能量的计算比较困难,造成评价的结果不准确。因此,寻找更为精确的电力系统稳定程度评价方法已成为研究的重点。
轨迹灵敏度是定量描述动态系统的参数对其动态轨迹影响程度的数学工具。轨迹灵敏度是系统动态特性的反应,同时也反映了参数对动态特性的影响。因此,充分挖掘电力系统轨迹灵敏度的动态信息可以有效评估电力系统的暂态稳定性。
广域量测系统                                                的不断发展为基于轨迹的电力系统动态分析提供了重要数据支持,能够实时记录电力系统受扰后的动态过程,且精度较高。基于轨迹计算轨迹灵敏度可以避免电力系统数值仿真过程,计算效率高。利用基于轨迹的轨迹灵敏度计算结果评估电力系统暂态稳定性可以使暂态稳定评估不受模型参数精度的影响。迄今未见根据实测轨迹计算轨迹灵敏度,用此灵敏度来评价系统稳定程度的方法的报道和实际应用。
发明内容
本发明所要解决的技术问题是,克服现有技术存在的不足,提出一种计算精度高,根据轨迹及轨迹灵敏度定量评价电力系统暂态稳定性方法,该方法基于系统轨迹,从而避免求解复杂的微分代数方程组和能量函数的计算,提高了计算精度。
解决其技术问题采用的技术方案是:
一种根据轨迹及轨迹灵敏度定量评价电力系统暂态稳定性方法,其特征是,它包含有以下步骤:
1).基于实测轨迹计算故障切除时间灵敏度
对于单机或多机系统,在时刻发生三相短路,时刻故障被切除,其各发电机功角对故障切除时间灵敏度按(1)式计算:
               (1)
对于多机系统,采用(1)式计算获得各发电机功角轨迹相对于故障切除时间的灵敏度,计算两台发电机间的相对灵敏度将两台发电机功角对故障切除时间灵敏度相减得到;
2).评价系统稳定程度
    由不同运行方式下同一时刻切除故障后功角对故障切除时间与电网临界切除时间关系可知,在单机系统下,功角对故障切除时间灵敏度最大值倒数与临界切除时间呈现正相关关系,亦即临界切除时间越大,功角对故障切除时间灵敏度最大值倒数越大,而临界切除时间越大说明系统抗干扰能力越强,系统就越稳定,因此将功角对故障切除时间灵敏度最大值倒数作为评价系统稳定程度,灵敏度最大值倒数越大,系统越稳定,稳定指标按(2)计算:
                                               (2)
其中,为系统稳定指标,为单机系统下发电机功角对故障切除时间灵敏度;
在多机系统下,首先采用以发电机角速度作为源数据,获取协方差矩阵及协方差矩阵的特征根和特征相量,由此确定发电机角速度的主成分,然后通过比较各发电机对主成分的载荷系数实现对发电机的同调分群的主成成分,将系统发电机分为两群,然后计算两群中距离最远两台发电机的功角对故障切除时间灵敏度,以该灵敏度最大值倒数作为系统稳定程度评价指标,稳定指标按(3)计算:
          )                   (3)
其中,为1群中第台发电机功角对故障切除时间灵敏度,为2群中第台发电机功角对故障切除时间灵敏度,为系统稳定指标。
本发明根据轨迹及轨迹灵敏度定量评价电力系统暂态稳定方法与传统的解微分代数方程组的评价方法不同,本发明仅利用了电力系统受扰后的轨迹和当前的故障切除时间即可评价出系统稳定程度,计算快速、评价精读高、具有在线应用前景。
附图说明
图为10机39节点系统接线示意图。
具体实施方式
下面对本发明作进一步说明。
本发明根据轨迹及轨迹灵敏度定量评价电力系统暂态稳定性方法,包含有以下步骤:
1.基于实测轨迹计算故障切除时间灵敏度
对于单机或多机系统,在时刻发生三相短路,时刻故障被切除,其各发电机功角对故障切除时间灵敏度按式(1)计算:
        (1)
对于多机系统,采用式(1)获得各发电机功角轨迹相对于故障切除时间的灵敏度,计算两台发电机的相对灵敏度将两台发电机功角对故障切除时间灵敏度相减得到;
2.评价系统稳定程度
    (1)由不同运行方式下同一时刻切除故障后功角对故障切除时间与电网临界切除时间关系可知,在单机系统下,功角对故障切除时间灵敏度最大值倒数与临界切除时间呈现正相关关系,亦即临界切除时间越大,功角对故障切除时间灵敏度最大值倒数越大,而临界切除时间越大说明系统抗干扰能力越强,系统就越稳定,因此将功角对故障切除时间灵敏度最大值倒数作为评价系统稳定程度,如式(2),灵敏度最大值倒数越大,系统越稳定;
                                                (2)
其中,为系统稳定指标,为单机系统下发电机功角对故障切除时间灵敏度;
(2)在多机系统下,首先采用以发电机角速度作为源数据,获取协方差矩阵及协方差矩阵的特征根和特征相量,由此确定发电机角速度的主成分,然后通过比较各发电机对主成分的载荷系数实现对发电机的同调分群的主成成分,将系统发电机分为两群,然后计算两群中距离最远两台发电机的功角对故障切除时间灵敏度,以该灵敏度最大值倒数作为系统稳定程度评价指标,稳定指标按(3)计算:
               ()             (3)
其中,为1群中第台发电机功角对故障切除时间灵敏度,为2群中第台发电机功角对故障切除时间灵敏度,为系统稳定指标。
具体应用实例:参见图1所示的10机39节点,在Bus12母线上发生100ms三相短路故障,使用PSAT对该故障进行仿真,发电机采用2阶模型,以仿真得到的发电机的功角和角速度轨迹作为轨迹,采用式(1)计算得到功角轨迹相对于故障切除时间的灵敏度,改变系统运行方式,得到不同运行方式下,以功角对故障切除时间灵敏度最大值倒数作为评价系统稳定程度指标得到的各运行方式下系统稳定程度排序和通过计算各运行方式下临界切除时间对不同运行方式下的排序结果完全吻合。

Claims (1)

1.一种根据轨迹及轨迹灵敏度定量评价电力系统暂态稳定性方法,其特征是,它包含有以下步骤:
1).基于实测轨迹计算故障切除时间灵敏度
对于单机或多机系统,在t0时刻发生三相短路,t1时刻故障被切除,其各发电机功角对故障切除时间灵敏度按(1)式计算:
Δ δ t 1 ( t ) * [ ϵ ( t - t 0 ) - ϵ ( t - t 1 ) ] = ω 0 · Δω ( t ) * ϵ ( t - t 1 ) - - - ( 1 )
对于多机系统,采用(1)式计算获得各发电机功角轨迹相对于故障切除时间的灵敏度,计算两台发电机间的相对灵敏度,将两台发电机功角对故障切除时间灵敏度相减得到;
2).评价系统稳定程度
由不同运行方式下同一时刻切除故障后功角对故障切除时间与电网临界切除时间关系可知,在单机系统下,功角对故障切除时间灵敏度最大值倒数与临界切除时间呈现正相关关系,亦即临界切除时间越大,功角对故障切除时间灵敏度最大值倒数越大,而临界切除时间越大说明系统抗干扰能力越强,系统就越稳定,因此将功角对故障切除时间灵敏度最大值倒数作为评价系统稳定程度,灵敏度最大值倒数越大,系统越稳定,稳定指标按(2)式计算:
S = 1 | max ( η ( t ) ) | - - - ( 2 )
其中,S为系统稳定指标,η(t)为单机系统下发电机功角对故障切除时间灵敏度;
在多机系统下,首先采用以发电机角速度作为源数据,获取协方差矩阵及协方差矩阵的特征根和特征相量,由此确定发电机角速度的主成分,然后通过比较各发电机对主成分的载荷系数实现对发电机的同调分群的主成成分,将系统发电机分为两群,然后计算两群中距离最远两台发电机的功角对故障切除时间灵敏度,以该灵敏度最大值倒数作为系统稳定程度评价指标,稳定指标按(3)式计算:
S = 1 | max ( η i ( t ) - η j ( t ) ) | ( i ≠ j , i , j = 1,2 , . . n ) - - - ( 3 )
其中,ηi(t)为1群中第i台发电机功角对故障切除时间灵敏度,ηj(t)为2群中第j台发电机功角对故障切除时间灵敏度,S为系统稳定指标。
CN201310065501.3A 2013-03-01 2013-03-01 根据轨迹及轨迹灵敏度定量评价电力系统暂态稳定性方法 Expired - Fee Related CN103178518B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310065501.3A CN103178518B (zh) 2013-03-01 2013-03-01 根据轨迹及轨迹灵敏度定量评价电力系统暂态稳定性方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310065501.3A CN103178518B (zh) 2013-03-01 2013-03-01 根据轨迹及轨迹灵敏度定量评价电力系统暂态稳定性方法

Publications (2)

Publication Number Publication Date
CN103178518A CN103178518A (zh) 2013-06-26
CN103178518B true CN103178518B (zh) 2014-08-13

Family

ID=48638202

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310065501.3A Expired - Fee Related CN103178518B (zh) 2013-03-01 2013-03-01 根据轨迹及轨迹灵敏度定量评价电力系统暂态稳定性方法

Country Status (1)

Country Link
CN (1) CN103178518B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104578053B (zh) * 2015-01-09 2017-01-04 北京交通大学 基于受扰电压轨迹簇特征的电力系统暂态稳定性预测方法
CN105071380A (zh) * 2015-07-28 2015-11-18 天津大学 一种电力系统广域同调辨识方法及其装置
CN105762796B (zh) * 2016-04-25 2018-06-19 武汉大学 一种解决区域电网故障后正确选切发电机和负荷的方法
CN107069697B (zh) * 2016-10-27 2019-06-25 西南交通大学 基于轨迹特征根的系统暂态稳定判断和扰动类型筛选方法
CN109038542A (zh) * 2018-06-27 2018-12-18 国网天津市电力公司电力科学研究院 一种基于轨迹灵敏度的在线安全稳定预防控制方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN100403618C (zh) * 2006-07-17 2008-07-16 天津大学 电力系统暂态稳定分析与控制方法
CN101635465A (zh) * 2009-07-22 2010-01-27 东北电力大学 根据稳定指标灵敏度计算电力系统暂态稳定极限值的方法
CN102570422B (zh) * 2012-01-05 2014-03-12 东北电力大学 根据实测轨迹计算电力系统极限切除时间的方法

Also Published As

Publication number Publication date
CN103178518A (zh) 2013-06-26

Similar Documents

Publication Publication Date Title
CN103178518B (zh) 根据轨迹及轨迹灵敏度定量评价电力系统暂态稳定性方法
CN104239681B (zh) 基于脉冲激励响应信号互功率谱函数的主轴系统运行模态分析方法
CN101282040B (zh) 非侵入式电力负荷实时分解方法
CN103033716B (zh) 一种电网综合负荷模型中各负荷成分所占比例的计算方法
CN101741086B (zh) 基于pmu数据的稳定计算模型的综合评估方法
CN105098763B (zh) 一种基于wams和ems的暂态电压稳定在线评估方法
CN102163844B (zh) 基于相量测量装置的电力系统状态检测方法
CN103592575A (zh) 一种基于多传感器系统的自适应加权数据融合故障测距方法
CN106443246A (zh) 基于pmu量测数据的小干扰稳定参数的在线辨识方法
CN103401238A (zh) 一种基于总体测辨法的电力负荷建模方法
CN103324841B (zh) 基于故障录波器数据的在线动态负荷建模方法
CN102664409B (zh) 一种基于实测数据风电场风电功率实时预测的计算方法
CN102364490B (zh) 基于层次分析模型的自动同调识别方法
CN101833602A (zh) 一种硅太阳能电池工程用数学模型的建模方法
CN106054023A (zh) 一种输电线路单端测距中估计两侧系统阻抗的方法
CN104142219A (zh) 一种基于多点脉冲激励的主轴系统运行模态分析方法
CN201600446U (zh) 采用瓦秒法检测谐波条件下电能表误差的装置
CN105045973A (zh) 一种变弧长自适应采样方法
CN101599644A (zh) 以电网正常或非正常运行信息辨识等值综合负荷模型
CN102510072B (zh) 一种电网系统暂态失稳判别方法
CN102842906B (zh) 一种综合负荷模型中电动机功率比例的计算方法
CN102280884B (zh) 一种电网等值方法
CN102570422B (zh) 根据实测轨迹计算电力系统极限切除时间的方法
CN101762801B (zh) 一种采用瓦秒法检测谐波条件下电能表误差的装置的检测方法
CN110098610B (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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140813

Termination date: 20200301