CN115864393B - 一种新能源场站暂态电压跌落严重度指标计算方法 - Google Patents

一种新能源场站暂态电压跌落严重度指标计算方法 Download PDF

Info

Publication number
CN115864393B
CN115864393B CN202211707733.XA CN202211707733A CN115864393B CN 115864393 B CN115864393 B CN 115864393B CN 202211707733 A CN202211707733 A CN 202211707733A CN 115864393 B CN115864393 B CN 115864393B
Authority
CN
China
Prior art keywords
new energy
transient voltage
voltage drop
calculating
subinterval
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
CN202211707733.XA
Other languages
English (en)
Other versions
CN115864393A (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.)
China University of Mining and Technology Beijing CUMTB
Economic and Technological Research Institute of State Grid Gansu Electric Power Co Ltd
Original Assignee
China University of Mining and Technology Beijing CUMTB
Economic and Technological Research Institute of State Grid Gansu 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 China University of Mining and Technology Beijing CUMTB, Economic and Technological Research Institute of State Grid Gansu Electric Power Co Ltd filed Critical China University of Mining and Technology Beijing CUMTB
Priority to CN202211707733.XA priority Critical patent/CN115864393B/zh
Publication of CN115864393A publication Critical patent/CN115864393A/zh
Application granted granted Critical
Publication of CN115864393B publication Critical patent/CN115864393B/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
    • 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/50Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种新能源场站暂态电压跌落严重度指标计算方法,属于新能源电力系统暂态电压稳定性分析技术领域,包括:将新能源场站电压跌落过程离散划分成多个子区间;计算各子区间内暂态电压的变化速率;确定暂降—恢复阶段各子区间系数;计算新能源场站暂态电压跌落严重度初始指标;确定电压暂态跌落过程中新能源出力场景;基于新能源出力场景对初始指标进行修正,得到新能源场站暂态电压跌落严重度指标。本发明提供的新能源场站暂态电压跌落严重度指标计算方法,综合考虑了暂态电压不同变化区间和新能源出力变化场景对稳定性的影响,能够进一步提高暂态电压跌落严重度和稳定性评估的准确性。

Description

一种新能源场站暂态电压跌落严重度指标计算方法
技术领域
本发明公开了一种新能源场站暂态电压跌落严重度指标计算方法,属于新能源电力系统暂态电压稳定性分析技术领域。
背景技术
越来越多的新能源场站并入电网,导致并网点近区无功支撑能力不足,电压暂态稳定性下降,易发生暂态电压跌落现象。当跌落深度较大、跌落时间较长、跌落速度较快、新能源出力较大时,更加容易导致新能源场站脱网运行,系统电压崩溃等一系列严重后果。新能源场站与常规电源并网点的暂态电压变化特性具有显著区别,增加了稳定性评估与判定的难度。因此,亟需提出一种新能源场站暂态电压跌落严重度指标计算方法,以准确量化评估其暂态电压稳定程度,为故障情况下的暂态电压稳定控制提供参考依据。
发明内容
本发明公开了一种新能源场站暂态电压跌落严重度指标计算方法,解决现有技术中,暂态电压跌落评估精度低的问题。
一种新能源场站暂态电压跌落严重度指标计算方法,包括:
步骤1:将新能源场站电压跌落过程离散划分成多个子区间;
步骤2:计算各子区间内暂态电压的变化速率;
步骤3:确定暂降-恢复阶段各子区间系数;
步骤4:计算新能源场站暂态电压跌落严重度初始指标;
步骤5:确定电压暂态跌落过程中新能源出力场景;
步骤6:基于新能源出力场景对初始指标进行修正,得到新能源场站暂态电压跌落严重度指标。
步骤1包括:
步骤1.1:根据新能源场站提供的暂态电压跌落过程时序曲线u(t)和电压跌落门槛值Ud确定暂态电压跌落严重度计算的起止时间,当曲线u(t)上的暂态电压数值跌落到Ud时作为开始,恢复到Ud时作为结束,得到起止区间[0,t],并统计新能源场站电压额定值UN以及区间内的最低跌落值Umin
步骤1.2:将[0,t]划分为:[0,t1),[t1,t2),[t2,t3),[t3,t4),[t4,t5),[t5,t],共6个子区间,其中t1为暂态电压由Ud跌落至的时刻,t2为暂态电压由/>跌落至/>的时刻,t3为暂态电压由/>跌落至Umin的时刻,t4为暂态电压由Umin恢复至/>的时刻,t5为暂态电压由/>恢复至/>的时刻。
步骤2包括:
步骤2.1:计算电压跌落门槛值Ud和最低跌落值Umin之间的变化速率v0作为基准速率值,其中
步骤2.2:计算各子区间内暂态电压首尾时刻之间的变化速率
i=1,2,3,...,6,其中Ui为ti时刻对应的暂态电压值,t0=0;
步骤2.3:基于变化速率vi,计算得到各子区间内暂态电压变化速率的标幺值vi′=vi/v0,i=1,…,6。
步骤3包括:
步骤3.1:计算各子区间内暂态电压平均值进而得到各子区间的严重度系数/>其中第一个子区间[0,t1)的严重度系数为1;
步骤3.2:计算各子区间的速率系数li,其中暂降阶段的速率系数与速率成正关系,即li=vi′,i=1,2,3,恢复阶段速率系数与速率成反关系,即li=1/v′i,i=4,5,6。
步骤4包括:
步骤4.1:基于各子区间的暂态电压跌落平均深度和区间长度,计算得到各子区间的严重度指标
步骤4.2:综合考虑各子区间的严重度系数ki和速率系数li,计算得到电压跌落过程严重度初始指标
步骤5包括:
步骤5.1:统计既定历史时间周期T内的新能源出力数值P;
步骤5.2:将历史新能源出力划分为小发、中发和大发三个区间,即L={P|Pmin≤P<P1}、M={P|P1≤P<P2}和H={P|P2≤P<Pmax},其中Pmin和Pmax分别为历史新能源出力的最小值和最大值,历史数据中出力小于P1的频率为0.2,即f(P)=0.2|(P<P1),出力小于P2的频率为0.8,即f(P)=0.8|(P<P2)。
步骤6包括:
步骤6.1:根据步骤5.2中划分的新能源出力区间,当新能源出力位于小发区间L时,新能源场站暂态电压跌落严重度指标τ′=0.9τ;
步骤6.2:当新能源出力位于中发区间M时,新能源场站暂态电压跌落严重度指标τ′=τ;
步骤6.3:当新能源出力位于大发区间H时,新能源场站暂态电压跌落严重度指标σ′=1.1σ。
本发明有益效果为:本发明提供的新能源场站暂态电压跌落严重度指标计算方法,综合考虑了暂态电压不同变化区间和新能源出力变化场景对稳定性的影响,能够进一步提高暂态电压跌落严重度和稳定性评估的准确性。
附图说明
图1是本发明的步骤流程图;
图2是暂态电压跌落过程中电压暂降-恢复阶段示意图。
具体实施方式
下面结合具体实施例对本发明的具体实施方式做进一步说明:
一种新能源场站暂态电压跌落严重度指标计算方法,如图1,包括:
步骤1:将新能源场站电压跌落过程离散划分成多个子区间;
步骤2:计算各子区间内暂态电压的变化速率;
步骤3:确定暂降-恢复阶段各子区间系数;
步骤4:计算新能源场站暂态电压跌落严重度初始指标;
步骤5:确定电压暂态跌落过程中新能源出力场景;
步骤6:基于新能源出力场景对初始指标进行修正,得到新能源场站暂态电压跌落严重度指标。
步骤1包括:
步骤1.1:根据新能源场站提供的暂态电压跌落过程时序曲线u(t)和电压跌落门槛值Ud确定暂态电压跌落严重度计算的起止时间,当曲线u(t)上的暂态电压数值跌落到Ud时作为开始,恢复到Ud时作为结束,得到起止区间[0,t],并统计新能源场站电压额定值UN以及区间内的最低跌落值Umin
步骤1.2:将[0,t]划分为:[0,t1),[t1,t2),[t2,t3),[t3,t4),[t4,t5),[t5,t],共6个子区间,其中t1为暂态电压由Ud跌落至的时刻,t2为暂态电压由/>跌落至/>的时刻,t3为暂态电压由/>跌落至Umin的时刻,t4为暂态电压由Umin恢复至/>的时刻,t5为暂态电压由/>恢复至/>的时刻。
步骤2包括:
步骤2.1:计算电压跌落门槛值Ud和最低跌落值Umin之间的变化速率v0作为基准速率值,其中
步骤2.2:计算各子区间内暂态电压首尾时刻之间的变化速率
i=1,2,3,...,6,其中Ui为ti时刻对应的暂态电压值,t0=0;
步骤2.3:基于变化速率vi,计算得到各子区间内暂态电压变化速率的标幺值vi′=vi/v0,i=1,…,6。
步骤3包括:
步骤3.1:计算各子区间内暂态电压平均值进而得到各子区间的严重度系数/>其中第一个子区间[0,t1)的严重度系数为1;
步骤3.2:计算各子区间的速率系数li,其中暂降阶段的速率系数与速率成正关系,即li=vi′,i=1,2,3,恢复阶段速率系数与速率成反关系,即li=1/v′i,i=4,5,6。
步骤4包括:
步骤4.1:基于各子区间的暂态电压跌落平均深度和区间长度,计算得到各子区间的严重度指标
步骤4.2:综合考虑各子区间的严重度系数ki和速率系数li,计算得到电压跌落过程严重度初始指标
步骤5包括:
步骤5.1:统计既定历史时间周期T内的新能源出力数值P;
步骤5.2:将历史新能源出力划分为小发、中发和大发三个区间,即L={P|Pmin≤P<P1}、M={P|P1≤P<P2}和H={P|P2≤P<Pmax},其中Pmin和Pmax分别为历史新能源出力的最小值和最大值,历史数据中出力小于P1的频率为0.2,即f(P)=0.2|(P<P1),出力小于P2的频率为0.8,即f(P)=0.8|(P<P2)。
步骤6包括:
步骤6.1:根据步骤5.2中划分的新能源出力区间,当新能源出力位于小发区间L时,新能源场站暂态电压跌落严重度指标τ′=0.9τ;
步骤6.2:当新能源出力位于中发区间M时,新能源场站暂态电压跌落严重度指标τ′=τ;
步骤6.3:当新能源出力位于大发区间H时,新能源场站暂态电压跌落严重度指标σ′=1.1σ。
实施例中,步骤1:暂态电压跌落曲线u(t)如图2所示,取Ud=0.8,UN=1.0,起止区间[0,t]=[0s,0.038s],暂态电压跌落门槛值Umin;6个子区间分别为
[0s,0.005s),[0.005s,0.009s),[0.009s,0.017s),[0.017s,0.022s),[0.022s,0.025s),[0.025s,0.038s]。
步骤2:计算基准速率值v0=11.76;计算子区间的电压变化速率v1=13.4,v2=16.5,v3=8.38,v4=13.4,v5=22,v6=5.15;计算子区间的电压变化速率标幺值v1′=1.14,v2′=1.40,v3′=0.71,v4′=1.14,v5′=1.87,v6′=0.44。
步骤3:计算各子区间内暂态电压平均值 计算各子区间的严重度系数k1=1,k2=1.1,k3=1.2,k4=1.2,k5=1.1,k6=1;计算各子区间的速率系数l1=1.14,l2=1.40,l3=0.71,l4=0.88,l5=0.53,l6=2.27。
步骤4:计算得到各子区间的严重度指标σ1=0.00115,σ2=0.0012,σ3=0.00288,σ4=0.0018,σ5=0.0009,σ6=0.00299;计算得到电压跌落过程严重度初始指标σ=0.01483。
步骤5:统计一天内96个新能源出力数值P,Pmin=312MW,Pmax=2568MW,P1=566MW,P2=1985MW。
步骤6:若电压跌落期间有功出力属于场景L={P|312MW≤P<566MW},则σ′=0.01334;若电压跌落期间有功出力属于场景M={P|566MW≤P<1985MW},则τ′=0.01483;若电压跌落期间有功出力属于场景H={P|1985MW≤P<2568MW},则:σ′=0.01631。
当然,上述说明并非是对本发明的限制,本发明也并不仅限于上述举例,本技术领域的技术人员在本发明的实质范围内所做出的变化、改型、添加或替换,也应属于本发明的保护范围。

Claims (4)

1.一种新能源场站暂态电压跌落严重度指标计算方法,其特征在于,包括:
步骤1:将新能源场站电压跌落过程离散划分成多个子区间;
步骤2:计算各子区间内暂态电压的变化速率;
步骤2.1:计算电压跌落门槛值Ud和最低跌落值Umin之间的变化速率v0作为基准速率值,其中
步骤2.2:计算各子区间内暂态电压首尾时刻之间的变化速率i=1,2,3,...,6,其中Ui为ti时刻对应的暂态电压值,t0=0;
步骤2.3:基于变化速率vi,计算得到各子区间内暂态电压变化速率的标幺值vi′=vi/v0,i=1,…,6;
步骤3:确定暂降—恢复阶段各子区间系数;
步骤3.1:计算各子区间内暂态电压平均值进而得到各子区间的严重度系数/>其中第一个子区间[0,t1)的严重度系数为1;
步骤3.2:计算各子区间的速率系数li,其中暂降阶段的速率系数与速率成正关系,即li=vi′,i=1,2,3,恢复阶段速率系数与速率成反关系,即li=1/v′i,i=4,5,6;
步骤4:计算新能源场站暂态电压跌落严重度初始指标;
步骤4.1:基于各子区间的暂态电压跌落平均深度和区间长度,计算得到各子区间的严重度指标
步骤4.2:综合考虑各子区间的严重度系数ki和速率系数li,计算得到电压跌落过程严重度初始指标
步骤5:确定电压暂态跌落过程中新能源出力场景;
步骤6:基于新能源出力场景对初始指标进行修正,得到新能源场站暂态电压跌落严重度指标。
2.根据权利要求1所述的一种新能源场站暂态电压跌落严重度指标计算方法,其特征在于,步骤1包括:
步骤1.1:根据新能源场站提供的暂态电压跌落过程时序曲线u(t)和电压跌落门槛值Ud确定暂态电压跌落严重度计算的起止时间,当曲线u(t)上的暂态电压数值跌落到Ud时作为开始,恢复到Ud时作为结束,得到起止区间[0,t],并统计新能源场站电压额定值UN以及区间内的最低跌落值Umin
步骤1.2:将[0,t]划分为:[0,t1),[t1,t2),[t2,t3),[t3,t4),[t4,t5),[t5,t],共6个子区间,其中t1为暂态电压由Ud跌落至的时刻,t2为暂态电压由/>跌落至/>的时刻,t3为暂态电压由/>跌落至Umin的时刻,t4为暂态电压由Umin恢复至/>的时刻,t5为暂态电压由/>恢复至/>的时刻。
3.根据权利要求1所述的一种新能源场站暂态电压跌落严重度指标计算方法,其特征在于,步骤5包括:
步骤5.1:统计既定历史时间周期T内的新能源出力数值P;
步骤5.2:将历史新能源出力划分为小发、中发和大发三个区间,即L={P|Pmin≤P<P1}、M={P|P1≤P<P2}和H={P|P2≤P<Pmax},其中Pmin和Pmax分别为历史新能源出力的最小值和最大值,历史数据中出力小于P1的频率为0.2,即f(P)=0.2|(P<P1),出力小于P2的频率为0.8,即f(P)=0.8|(P<P2)。
4.根据权利要求3所述的一种新能源场站暂态电压跌落严重度指标计算方法,其特征在于,步骤6包括:
步骤6.1:根据步骤5.2中划分的新能源出力区间,当新能源出力位于小发区间L时,新能源场站暂态电压跌落严重度指标σ′=0.9σ;
步骤6.2:当新能源出力位于中发区间M时,新能源场站暂态电压跌落严重度指标σ′=σ;
步骤6.3:当新能源出力位于大发区间H时,新能源场站暂态电压跌落严重度指标σ′=1.1σ。
CN202211707733.XA 2022-12-29 2022-12-29 一种新能源场站暂态电压跌落严重度指标计算方法 Active CN115864393B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211707733.XA CN115864393B (zh) 2022-12-29 2022-12-29 一种新能源场站暂态电压跌落严重度指标计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211707733.XA CN115864393B (zh) 2022-12-29 2022-12-29 一种新能源场站暂态电压跌落严重度指标计算方法

Publications (2)

Publication Number Publication Date
CN115864393A CN115864393A (zh) 2023-03-28
CN115864393B true CN115864393B (zh) 2023-08-15

Family

ID=85655914

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211707733.XA Active CN115864393B (zh) 2022-12-29 2022-12-29 一种新能源场站暂态电压跌落严重度指标计算方法

Country Status (1)

Country Link
CN (1) CN115864393B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106291244A (zh) * 2016-08-24 2017-01-04 浙江群力电气有限公司 一种多次电压暂降严重度评估方法及系统
CN111537839A (zh) * 2020-06-10 2020-08-14 河海大学 适用于高比例新能源送端电网的暂态电压安全性评估方法及安全性评估系统
CN112510688A (zh) * 2020-11-20 2021-03-16 东北电力大学 基于相位校正李雅普诺夫指数的暂态电压稳定监测方法
CN112952810A (zh) * 2021-02-24 2021-06-11 中国电力科学研究院有限公司 一种风电暂态过程划分方法及装置
CN113949069A (zh) * 2021-12-20 2022-01-18 中国电力科学研究院有限公司 确定高比例新能源电力系统暂态电压稳定性的方法及系统

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10054647B2 (en) * 2014-09-15 2018-08-21 Atmel Corporation Fault detection

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106291244A (zh) * 2016-08-24 2017-01-04 浙江群力电气有限公司 一种多次电压暂降严重度评估方法及系统
CN111537839A (zh) * 2020-06-10 2020-08-14 河海大学 适用于高比例新能源送端电网的暂态电压安全性评估方法及安全性评估系统
CN112510688A (zh) * 2020-11-20 2021-03-16 东北电力大学 基于相位校正李雅普诺夫指数的暂态电压稳定监测方法
CN112952810A (zh) * 2021-02-24 2021-06-11 中国电力科学研究院有限公司 一种风电暂态过程划分方法及装置
CN113949069A (zh) * 2021-12-20 2022-01-18 中国电力科学研究院有限公司 确定高比例新能源电力系统暂态电压稳定性的方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
薛安成 等.采用多二元表判据的实用暂态电压稳定裕度指标研究.《中国电机工程学报》.2018,第14卷(第38期),第4117-4125+4317页. *

Also Published As

Publication number Publication date
CN115864393A (zh) 2023-03-28

Similar Documents

Publication Publication Date Title
CN107885959B (zh) 一种基于置信等效功率曲线带的风电建模及性能评估方法
CN107918103B (zh) 一种基于灰色粒子滤波的锂离子电池剩余寿命预测方法
CN111965559B (zh) 一种锂离子电池soh在线估计方法
US20170356419A1 (en) Method, device and system for compensating output power of wind turbine generator set
CN110826791A (zh) 一种基于长短时记忆神经网络的混合风功率预测方法
CN111428201B (zh) 基于经验模态分解和前馈神经网络对时序数据的预测方法
CN106570790B (zh) 一种计及风速数据分段特性的风电场出力数据修复方法
CN110289631B (zh) 一种风电场储能装置容量的计算方法及系统
CN111293686A (zh) 基于armax系统辨识的电力系统惯量实时评估方法
CN107947228B (zh) 基于Markov理论的含风电电力系统随机稳定性分析方法
CN107221933B (zh) 一种概率潮流计算方法
CN115864393B (zh) 一种新能源场站暂态电压跌落严重度指标计算方法
CN113791351B (zh) 基于迁移学习和差值概率分布的锂电池寿命预测方法
CN111030109A (zh) 一种网络约束的电网备用容量减扣计算方法及相关装置
CN113505909A (zh) 一种面向短期风电功率趋势预测的误差补偿方法
CN101529347A (zh) 识别有延时的调节对象的方法和装置、调节装置和计算机程序产品
CN111044907B (zh) 一种基于微片段数据和电压滤波的soh统计方法
CN112763918A (zh) 电动汽车动力电池放电功率控制方法
CN111931967A (zh) 一种风电场短期功率预测方法
CN111756039B (zh) 基于概率统计的新能源电力系统惯量估计方法
CN114285091B (zh) 一种包含多个光伏发电的区域电网数据采集异常检测方法
CN112257997B (zh) 一种基于神经网络的pcs转化率计算方法及系统
CN113489005B (zh) 配电网潮流计算的配变负荷估算方法及系统
CN109378817B (zh) 一种基于一次调频稳定性评价方法
CN114880790A (zh) 一种汽动引风机变工况机理建模方法、系统及存储介质

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant