CN102570422B - 根据实测轨迹计算电力系统极限切除时间的方法 - Google Patents
根据实测轨迹计算电力系统极限切除时间的方法 Download PDFInfo
- Publication number
- CN102570422B CN102570422B CN201210002069.9A CN201210002069A CN102570422B CN 102570422 B CN102570422 B CN 102570422B CN 201210002069 A CN201210002069 A CN 201210002069A CN 102570422 B CN102570422 B CN 102570422B
- Authority
- CN
- China
- Prior art keywords
- delta
- track
- centerdot
- clearing time
- generator
- 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
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明是一种根据实测轨迹计算电力系统极限切除时间的方法,其特点是:首先推导出电力系统中功角轨迹相对于故障切除时间灵敏度与轨迹之间的关系,基于实测轨迹获取功角轨迹相对于故障切除时间的各阶灵敏度,然后根据计算出的轨迹灵敏度预测出其他故障切除时间下的功角轨迹,最后通过将预测出的各发电机功角轨迹之间做差,判断电力系统的暂态稳定性,从而确定极限切除时间。该方法首先根据实测轨迹获取轨迹灵敏度,并以轨迹灵敏度为工具,估计出电力系统的极限切除时间,避免了现有的极限切除时间计算方法对电力系统仿真模型的依赖。该方法计算简单,为基于轨迹的电力系统动态分析提供重要的理论支持,具有灵敏度高、精度高、应用价值高等优点。
Description
技术领域
本发明涉及一种根据实测轨迹计算电力系统极限切除时间的方法。
背景技术
极限切除时间是电力系统暂态稳定分析与控制中十分重要的一个参数,它表示电力系统遭受某一故障后在不失去稳定的前提下能够承受的最长故障持续时间。根据极限切除时间和当前的故障切除时间可计算出电力系统的稳定裕度,进而可以设计合理的暂态稳定控制措施,提高电力系统运行的安全水平。
现有的极限切除时间求取方法均是依赖于电力系统的模型,其计算结果受系统模型的精度影响较大。而国内外多次电力系统事故后分析和大扰动试验结果均表明现有的电力系统模型存在一定的误差,使得基于现有模型的动态仿真结果与实测轨迹不相吻合。因此,寻找更为精确的电力系统极限切除时间方法已成为研究的重点。
广域量测系统WAMS的不断发展为基于轨迹的电力系统动态分析提供了重要数据支持,WAMS能够实时记录电力系统受扰后的动态过程,且精度较高。在现有技术中,研究人员利用广域量测系统,提出了基于轨迹与轨迹灵敏度之间的关系计算轨迹灵敏度的方法,但是,这种轨迹灵敏度的计算方法只给出基于轨迹计算一阶轨迹灵敏度,而无法计算更高阶次的轨迹灵敏度。迄今未见根据轨迹计算高阶轨迹灵敏度的方法,特别是未见根据轨迹直接预测其他故障切除时间下的动态变量轨迹,从而确定故障极限切除时间方法的报道和实际应用。
发明内容
本发明所要解决的技术问题是,克服现有技术存在的不足,提出一种计算精度高,根据实测轨迹计算电力系统极限切除时间的方法,该方法可以避免极限切除时间及轨迹灵敏度计算对电力系统模型的依赖,从而提高了计算精度。
解决其技术问题采用的技术方案是:一种根据实测轨迹计算电力系统极限切除时间的方法,其特征是:它包含有以下步骤:
1).基于实测轨迹计算故障切除时间灵敏度
对于单机或多机系统,在t0时刻发生三相短路,t1时刻故障被切除,则整个过程中系统经历了短路前、短路中和短路切除后3个阶段;
设在短路发生前发电机内电势点与无穷大点之间的阻抗为xΙ,在短路期间变成了xΙΙ,短路切除后又变为xΙ,则故障前后发电机内电势点至无穷大节点出的阻抗变化值Δx可表示为:
Δx=ΔxΙ[ε(t-t0)-ε(t-t1)] (1)
其中:
ΔxΙ=xΙ-xΙΙ (2)
将单机系统的动态响应看成是由发电机内电势点至无穷大节点之间的电抗变化这一激励引发,利用线性系统的理论获得轨迹相对于故障切除时间灵敏度和轨迹之间的关系;
设发电机为2阶模型,不考虑励磁系统模型和调速系统模型,其数学模型为:
其中Mm、Me分别为发电机的机械功率和电磁功率,Tj和D分别为发电机惯性时间常数和阻尼系数。其中电磁功率可表达为:
对其进行线性化并转换到复频域,则模型为:
ΔMe(s)=k1Δδ(s)+kxΔx(s) (7)
其中k1、kx均为常系数;
由于不考虑调速器模型,发电机的机械功率不变,即ΔMm=0,
整理式(5)~(7)得:
F1(s)·Δδ(s)=kx·Δx(s) (8)
其中:
两边同时对故障切除时间t1求导,可得:
联合式(8)~(10),可得:
进行拉普拉斯反变换,即:
式(12)两端同时对t1再次进行求导,可得:
其中 分别表示功角轨迹相对于t1的2阶、3阶灵敏度,是待求量;Δδ′′(t)、Δδ′′′(t)分别表示功角轨迹的2阶、3阶导数,根据实际功角轨迹计算出来; 分别表示功角轨迹的一阶导数相对于t1的2阶灵敏度和功角轨迹的2阶导数相对于t1的1阶灵敏度;
式(12)~(14)表明直接根据功角轨迹和阶跃函数轨迹即可计算出功角轨迹相对于故障切除时间t1的1、2、3阶灵敏度,同理也可进一步求得功角轨迹相对于tcl的其他更高阶灵敏度;
对于多机系统,可采用以上方法获得各发电机功角轨迹相对于故障切除时间的灵敏度。
2).其他故障切除时间下受扰轨迹的计算
设电力系统当前的故障切除时间为tcl,在已知该故障切除时间下的功角轨迹δ(t,tcl)后,可预测其他故障切除时间下的功角轨迹:
对于多机系统,设有N台发电机,也可根据以上的分析方法求得任意故障切除时间下每台发电机的功角轨迹;
3).对于电力系统暂态稳定性的判断
当多机系统中任意两台发电机功角之差的最大值大于180°或2π(rad)时,表明此电力系统失去功角稳定,则暂态稳定性的判别方法为
max|δij(t)|<2π (16)
4).对于实测轨迹的极限切除时间计算
(a)获取当前的故障切除时间tcl,以及当前故障切除时间下的受扰轨迹δ(t,tcl),根据式(12)~(14)计算出功角轨迹相对于故障切除时间tcl的各阶灵敏度;
(b)设某一待求的故障切除时间t0,根据式(15)计算出t0下电力系统中各台发电机的功角轨迹δi(t0,tcl+Δt);
(c)计算各台发电机功角轨迹之差Δδij(t0,tcl+Δt),并判断max|Δδij(t0,tcl+Δt)|<2π(rad),若成立,则增大t0,返回(b),直至max|Δδij(t0,tcl+Δt)|=2π,则此时的t0为极限切除时间;
(d)若max|Δδij(t0,tcl+Δt)|>2π,减小t0,返回(b),直至max|Δδij(t0,tcl+Δt)|=2π,则此时的t0为极限切除时间。
本发明的根据实测轨迹计算电力系统极限切除时间的方法与传统的基于模型的极限切除时间计算方法不同,本发明仅利用了电力系统受扰后的实测轨迹和当前的故障切除时间即可计算出极限切除时间,准确判断其暂态稳定性,从而避免了极限切除时间的计算对系统模型的依赖,具有灵敏度高、精度高、应用价值高等优点。
附图说明
图为单机无穷大系统接线示意图。
具体实施方式
下面对本发明作进一步说明。
本发明根据实测轨迹计算电力系统极限切除时间的方法,包含有以下步骤:
1.基于实测轨迹计算故障切除时间灵敏度
对于单机或多机系统,在t0时刻发生三相短路,t1时刻故障被切除,则整个过程中系统经历了短路前、短路中和短路切除后3个阶段;
设在短路发生前发电机内电势点与无穷大点之间的阻抗为xΙ,在短路期间变成了xΙΙ,短路切除后又变为xΙ,则故障前后发电机内电势点至无穷大节点出的阻抗变化值Δx可表示为:
Δx=ΔxΙ[ε(t-t0)-ε(t-t1)] (1)
其中:
ΔxΙ=xΙ-xΙΙ (2)
将单机系统的动态响应看成是由发电机内电势点至无穷大节点之间的电抗变化这一激励引发,利用线性系统的理论获得轨迹相对于故障切除时间灵敏度和轨迹之间的关系;
设发电机为2阶模型,不考虑励磁系统模型和调速系统模型,其数学模型为:
其中Mm、Me分别为发电机的机械功率和电磁功率,Tj和D分别为发电机惯性时间常数和阻尼系数。其中电磁功率可表达为:
对其进行线性化并转换到复频域,则模型为:
ΔMe(s)=k1Δδ(s)+kxΔx(s) (7)
其中k1、kx均为常系数;
由于不考虑调速器模型,发电机的机械功率不变,即ΔMm=0,
整理式(5)~(7)得:
F1(s)·Δδ(s)=kx·Δx(s) (8)
其中:
两边同时对故障切除时间t1求导,可得:
联合式(8)~(10),可得:
进行拉普拉斯反变换,即:
式(12)两端同时对t1再次进行求导,可得:
其中 分别表示功角轨迹相对于t1的2阶、3阶灵敏度,是待求量;Δδ′′(t)、Δδ′′′(t)分别表示功角轨迹的2阶、3阶导数,根据实际功角轨迹计算出来; 分别表示功角轨迹的一阶导数相对于t1的2阶灵敏度和功角轨迹的2阶导数相对于t1的1阶灵敏度;
式(12)~(14)表明直接根据功角轨迹和阶跃函数轨迹即可计算出功角轨迹相对于故障切除时间t1的1、2、3阶灵敏度,从而避免基于模型微分进行轨迹灵敏度计算,同理也可进一步求得功角轨迹相对于tcl的其他更高阶灵敏度;
对于多机系统,可采用以上方法获得各发电机功角轨迹相对于故障切除时间的灵敏度。
2.其他故障切除时间下受扰轨迹的计算
设电力系统当前的故障切除时间为tcl,在已知该故障切除时间下的功角轨迹δ(t,tcl)后,可预测其他故障切除时间下的功角轨迹:
对于多机系统,设有N台发电机,也可根据以上的分析方法求得任意故障切除时间下每台发电机的功角轨迹;
3.对于电力系统暂态稳定性的判断
当多机系统中任意两台发电机功角之差的最大值大于180°或2π(rad)时,表明此电力系统失去功角稳定,则暂态稳定性的判别方法为
max|δij(t)|<2π (16)
4.对于实测轨迹的极限切除时间计算
(a)获取当前的故障切除时间tcl,以及当前故障切除时间下的受扰轨迹δ(t,tcl),根据式(12)~(14)计算出功角轨迹相对于故障切除时间tcl的各阶灵敏度;
(b)设某一待求的故障切除时间t0,根据式(15)计算出t0下电力系统中各台发电机的功角轨迹δi(t0,tcl+Δt);
(c)计算各台发电机功角轨迹之差Δδij(t0,tcl+Δt),并判断max|Δδij(t0,tcl+Δt)|<2π(rad),若成立,则增大t0,返回(b),直至max|Δδij(t0,tcl+Δt)|=2π,则此时的t0为极限切除时间;
(d)若max|Δδij(t0,tcl+Δt)|>2π,减小t0,返回(b),直至max|Δδij(t0,tcl+Δt)|=2π,则此时的t0为极限切除时间。
具体应用实例:参见图1所示的单机无穷大系统,在Bus02母线上发生100ms三相短路故障,使用电力系统分析综合程序(PSASP)对该故障进行仿真,发电机采用2阶模型,以仿真得到的发电机G1的功角和角速度轨迹作为实测轨迹,采用卷积方法得到功角轨迹相对于故障切除时间的各阶灵敏度,并根据式(12)~(14)计算出其他故障切除时间下发电机G1的功角轨迹,最后进行比较得到的极限故障切除时间为0.342s,而根据仿真分析计算得到的极限切除时间为0.35s,可见二者是比较吻合的。
Claims (1)
1.一种根据实测轨迹计算电力系统极限切除时间的方法,其特征是:它包含有以下步骤:
1).基于实测轨迹计算故障切除时间灵敏度
对于单机或多机系统,在t0时刻发生三相短路,t1时刻故障被切除,则整个过程中系统经历了短路前、短路中和短路切除后3个阶段;
设在短路发生前发电机内电势点与无穷大点之间的阻抗为xΙ,在短路期间变成了xΙΙ,短路切除后又变为xΙ,则故障前后发电机内电势点至无穷大节点出的阻抗变化值Δx可表示为:
Δx=ΔxΙ[ε(t-t0)-ε(t-t1)] (1)
其中:
ΔxΙ=xΙ-xΙΙ (2)
将单机系统的动态响应看成是由发电机内电势点至无穷大节点之间的电抗变化这一激励引发,利用线性系统的理论获得轨迹相对于故障切除时间灵敏度和轨迹之间的关系;
设发电机为2阶模型,不考虑励磁系统模型和调速系统模型,其数学模型为:
其中Mm、Me分别为发电机的机械功率和电磁功率,Tj和D分别为发电机惯性时间常数和阻尼系数,其中电磁功率可表达为:
对其进行线性化并转换到复频域,则模型为:
ΔMe(s)=k1Δδ(s)+kxΔx(s) (7)
其中k1、kx均为常系数;
由于不考虑调速器模型,发电机的机械功率不变,即ΔMm=0,
整理式(5)~(7)得:
F1(s)·Δδ(s)=kx·Δx(s) (8)
其中:
两边同时对故障切除时间t1求导,可得:
联合式(8)~(10),可得:
进行拉普拉斯反变换,即:
式(12)两端同时对t1再次进行求导,可得:
其中分别表示功角轨迹相对于t1的2阶、3阶灵敏度,是待求量;Δδ′′(t)、Δδ′′′(t)分别表示功角轨迹的2阶、3阶导数,根据实际功角轨迹计算出来;分别表示功角轨迹的一阶导数相对于t1的2阶灵敏度和功角轨迹的2阶导数相对于t1的1阶灵敏度;
式(12)~(14)表明直接根据功角轨迹和阶跃函数轨迹即可计算出功角轨迹相对于故障切除时间t1的1、2、3阶灵敏度,同理也可进一步求得功角轨迹相对于tcl的其他更高阶灵敏度;
对于多机系统,可采用以上方法获得各发电机功角轨迹相对于故障切除时间的灵敏度;
2).其他故障切除时间下受扰轨迹的计算
设电力系统当前的故障切除时间为tcl,在已知该故障切除时间下的功角轨迹δ(t,tcl)后,可预测其他故障切除时间下的功角轨迹:
对于多机系统,设有N台发电机,也可根据以上的分析方法求得任意故障切除时间下每台发电机的功角轨迹;
3).对于电力系统暂态稳定性的判断
当多机系统中任意两台发电机功角之差的最大值大于180°或2π(rad)时,表明此电力系统失去功角稳定,则暂态稳定性的判别方法为
max|δij(t)|<2π (16)
4).对于实测轨迹的极限切除时间计算(a)
获取当前的故障切除时间tcl,以及当前故障切除时间下的受扰轨迹δ(t,tcl),根据式(12)~(14)计算出功角轨迹相对于故障切除时间tcl的各阶灵敏度;
(b)设某一待求的故障切除时间t0,根据式(15)计算出t0下电力系统中各台发电机的功角轨迹δi(t0,tcl+Δt);
(c)计算各台发电机功角轨迹之差Δδij(t0,tcl+Δt),并判断max|Δδij(t0,tcl+Δt)|<2π(rad),若成立,则增大t0,返回(b),直至max|Δδij(t0,tcl+Δt)|=2π,则此时的t0为极限切除时间;
(d)若max|Δδij(t0,tcl+Δt)|>2π,减小t0,返回(b),直至max|Δδij(t0,tcl+Δt)|=2π,则此时的t0为极限切除时间。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210002069.9A CN102570422B (zh) | 2012-01-05 | 2012-01-05 | 根据实测轨迹计算电力系统极限切除时间的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210002069.9A CN102570422B (zh) | 2012-01-05 | 2012-01-05 | 根据实测轨迹计算电力系统极限切除时间的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102570422A CN102570422A (zh) | 2012-07-11 |
CN102570422B true CN102570422B (zh) | 2014-03-12 |
Family
ID=46415157
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210002069.9A Expired - Fee Related CN102570422B (zh) | 2012-01-05 | 2012-01-05 | 根据实测轨迹计算电力系统极限切除时间的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102570422B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103178518B (zh) * | 2013-03-01 | 2014-08-13 | 东北电力大学 | 根据轨迹及轨迹灵敏度定量评价电力系统暂态稳定性方法 |
CN103280783B (zh) * | 2013-05-09 | 2015-07-22 | 哈尔滨工程大学 | 基于临界机组对的故障临界切除时间计算方法 |
CN109038542A (zh) * | 2018-06-27 | 2018-12-18 | 国网天津市电力公司电力科学研究院 | 一种基于轨迹灵敏度的在线安全稳定预防控制方法 |
CN115146233B (zh) * | 2022-07-18 | 2024-07-02 | 广西电网有限责任公司电力科学研究院 | 一种电力系统中极限切除时间的批量计算方法和相关装置 |
CN116526583B (zh) * | 2023-07-03 | 2023-09-26 | 国网北京市电力公司 | 一种光热最小装机容量配置方法、系统、设备及介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3987359A (en) * | 1974-12-23 | 1976-10-19 | Westinghouse Electric Corporation | Phase angle function generator for static VAR control |
CN1110838A (zh) * | 1995-02-25 | 1995-10-25 | 电力部电力自动化研究院 | 电力系统暂态稳定的在线监视和预防控制系统 |
EP1691462A2 (en) * | 2005-02-09 | 2006-08-16 | The Tokyo Electric Power Co., Inc. | Group-based BCU methods for on-line dynamical security assessments and energy margin calculations of practical power systems |
CN1901317A (zh) * | 2006-07-17 | 2007-01-24 | 天津大学 | 电力系统暂态稳定分析与控制方法 |
CN101635465A (zh) * | 2009-07-22 | 2010-01-27 | 东北电力大学 | 根据稳定指标灵敏度计算电力系统暂态稳定极限值的方法 |
CN101969199A (zh) * | 2010-08-26 | 2011-02-09 | 天津大学 | 用于暂态功角稳定风险评估的故障损失估计方法 |
CN102035203A (zh) * | 2010-10-20 | 2011-04-27 | 中国电力科学研究院 | 一种快速求取联络线暂态稳定极限传输功率的计算方法 |
-
2012
- 2012-01-05 CN CN201210002069.9A patent/CN102570422B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3987359A (en) * | 1974-12-23 | 1976-10-19 | Westinghouse Electric Corporation | Phase angle function generator for static VAR control |
CN1110838A (zh) * | 1995-02-25 | 1995-10-25 | 电力部电力自动化研究院 | 电力系统暂态稳定的在线监视和预防控制系统 |
EP1691462A2 (en) * | 2005-02-09 | 2006-08-16 | The Tokyo Electric Power Co., Inc. | Group-based BCU methods for on-line dynamical security assessments and energy margin calculations of practical power systems |
CN1901317A (zh) * | 2006-07-17 | 2007-01-24 | 天津大学 | 电力系统暂态稳定分析与控制方法 |
CN101635465A (zh) * | 2009-07-22 | 2010-01-27 | 东北电力大学 | 根据稳定指标灵敏度计算电力系统暂态稳定极限值的方法 |
CN101969199A (zh) * | 2010-08-26 | 2011-02-09 | 天津大学 | 用于暂态功角稳定风险评估的故障损失估计方法 |
CN102035203A (zh) * | 2010-10-20 | 2011-04-27 | 中国电力科学研究院 | 一种快速求取联络线暂态稳定极限传输功率的计算方法 |
Non-Patent Citations (4)
Title |
---|
刘洪波.基于量测轨迹及轨迹灵敏度的电力系统动态分析.《中国博士学位论文全文数据库》.2007,第16-18、37-44、61-65页. |
基于轨迹灵敏度的电力系统动态安全预防控制算法研究;孙景强等;《电网技术》;20041130;第28卷(第21期);第26-30页 * |
基于量测轨迹及轨迹灵敏度的电力系统动态分析;刘洪波;《中国博士学位论文全文数据库》;20071205;第16-18、37-44、61-65页 * |
孙景强等.基于轨迹灵敏度的电力系统动态安全预防控制算法研究.《电网技术》.2004,第28卷(第21期),第26-30页. |
Also Published As
Publication number | Publication date |
---|---|
CN102570422A (zh) | 2012-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102570422B (zh) | 根据实测轨迹计算电力系统极限切除时间的方法 | |
Zeng et al. | Online identification of inertia distribution in normal operating power system | |
CN100486073C (zh) | 考虑负荷电压特性的非线性动态状态估计算法 | |
CN102163844B (zh) | 基于相量测量装置的电力系统状态检测方法 | |
CN103645422B (zh) | 一种发电厂内部扰动引起电网强迫功率振荡在线分析方法 | |
CN103474992B (zh) | 电力系统节点电压稳定状态的实时在线辨识判据 | |
CN102901929A (zh) | 用于计算电池单元阻抗的装置及电池阻抗测量系统 | |
CN105206888A (zh) | 一种锂离子电池内部温度监测方法 | |
CN102520353A (zh) | 同步发电机模型参数分步辨识方法 | |
CN102570428B (zh) | 基于电子式互感器微分输出的故障测距与距离保护方法 | |
CN102608542A (zh) | 动力电池荷电状态估计方法 | |
CN102680795A (zh) | 一种二次电池内阻的实时在线估计方法 | |
CN102013703B (zh) | 基于量测注入的发电机模型参数校核方法 | |
CN102779238A (zh) | 一种基于自适应卡尔曼滤波的无刷直流电机系统辨识方法 | |
CN103033716A (zh) | 一种电网综合负荷模型中各负荷成分所占比例的计算方法 | |
CN103543376A (zh) | 用于小电流接地系统故障选线的径向基神经网络方法 | |
CN107565867A (zh) | 一种基于轨迹灵敏度的同步发电机参数辨识方法 | |
CN103972912A (zh) | 一种含风电电力系统频率响应的频域分析方法 | |
CN103178518B (zh) | 根据轨迹及轨迹灵敏度定量评价电力系统暂态稳定性方法 | |
CN105445614B (zh) | 一种基于小波分析的双端行波故障定位的方法和系统 | |
CN102570457B (zh) | 一种基于内网量测数据的两端口外网静态等值方法 | |
CN105244887A (zh) | 电力系统暂态功角失稳的闭环控制方法 | |
CN103036252A (zh) | 基于轨迹灵敏度法的风电场模型校核系统及方法 | |
Lin et al. | Experience with synchronous generator parameter identification using a Kalman filter | |
CN102280884A (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: 20140312 Termination date: 20220105 |
|
CF01 | Termination of patent right due to non-payment of annual fee |