CN114840982A - 台风灾害下考虑连锁故障演化的电力系统弹性评估方法 - Google Patents
台风灾害下考虑连锁故障演化的电力系统弹性评估方法 Download PDFInfo
- Publication number
- CN114840982A CN114840982A CN202210397176.XA CN202210397176A CN114840982A CN 114840982 A CN114840982 A CN 114840982A CN 202210397176 A CN202210397176 A CN 202210397176A CN 114840982 A CN114840982 A CN 114840982A
- Authority
- CN
- China
- Prior art keywords
- typhoon
- fault
- line
- power
- transmission line
- 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.)
- Pending
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 30
- 238000000034 method Methods 0.000 claims abstract description 90
- 230000005540 biological transmission Effects 0.000 claims abstract description 72
- 230000008569 process Effects 0.000 claims abstract description 47
- 238000011084 recovery Methods 0.000 claims abstract description 32
- 230000001052 transient effect Effects 0.000 claims abstract description 15
- 238000004458 analytical method Methods 0.000 claims abstract description 9
- 238000000342 Monte Carlo simulation Methods 0.000 claims abstract description 4
- 238000004088 simulation Methods 0.000 claims description 54
- 230000004224 protection Effects 0.000 claims description 28
- 230000009471 action Effects 0.000 claims description 11
- 238000009826 distribution Methods 0.000 claims description 9
- 230000009467 reduction Effects 0.000 claims description 8
- 238000005520 cutting process Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 230000007423 decrease Effects 0.000 claims description 5
- 230000005251 gamma ray Effects 0.000 claims description 3
- 238000010248 power generation Methods 0.000 claims description 3
- 230000006641 stabilisation Effects 0.000 claims description 3
- 238000011105 stabilization Methods 0.000 claims description 3
- 238000007619 statistical method Methods 0.000 claims description 3
- 238000011161 development Methods 0.000 abstract description 7
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000004364 calculation method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000007123 defense Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 229910052731 fluorine Inorganic materials 0.000 description 1
- 125000001153 fluoro group Chemical group F* 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/04—Power grid distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS 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/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Marketing (AREA)
- Educational Administration (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Game Theory and Decision Science (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Health & Medical Sciences (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种台风灾害下考虑连锁故障演化的电力系统弹性评估方法,涉及电力系统分析领域。该方法通过台风灾害下输电线路故障模型生成故障场景;利用基于短、长两时间尺度的连锁故障演化模型模拟电力系统因暂态稳定和热稳定问题引起的连锁故障及其演化过程;结合电力系统恢复模型,采用序贯蒙特卡洛模拟法,根据系统性能情况评估电力系统弹性指标值。相对于未考虑连锁故障演化的电力系统弹性评估方法而言,本发明能够有效模拟台风灾害导致的初始故障情况及连锁故障演化过程,对于故障场景的刻画更全面客观,更贴近于台风灾害下电力系统实际故障发展情况,从而使弹性评估结果更准确,有利于电网运行人员准确评估电力系统应对台风灾害的能力,制定有针对性的恢复策略。
Description
技术领域
本发明涉及电力系统分析领域,具体为一种台风灾害下考虑连锁故障演化的电力系统弹性评估方法。
背景技术
全球气温升高将导致台风灾害发生频率提高,显著增加了电力系统发生大面积停电的风险,严重影响社会经济健康发展和人民群众正常生产生活。因此,为准确评估电力系统应对日益频发台风灾害的能力,迫切需要深化电力系统弹性评估研究。基于现有电力系统弹性评估研究对连锁故障演化过程考虑不够全面的问题,可以看出,在电力系统弹性评估过程中考虑到连锁故障演化过程,对准确地感知电力系统应对台风灾害的能力具有重要的意义。
专利202011599832.1“一种高分辨度电网弹性评估方法”及专利202110659995.2“一种极端灾害条件下的电力系统弹性评估方法和系统”在弹性评估方法中未考虑可能发生的连锁故障及其演化过程;专利201910295707.2“台风天气条件下电力系统的多时间尺度连锁故障预测方法”未考虑故障恢复过程,未对电力系统弹性进行评估。基于上述评估方法或系统中所存在的缺陷和不足,需要改进在台风灾害下对电力系统的弹性评估方法。
发明内容
本发明为了解决现有电力系统弹性评估方法对连锁故障演化过程考虑不足,得出的弹性评估结果偏乐观的问题,本发明在电力系统弹性评估过程中计及可能发生的连锁故障演化过程,提供了一种台风灾害下考虑连锁故障演化的电力系统弹性评估方法。
图1为台风灾害下电力系统性能变化情况,T0为台风开始侵袭系统时刻,T1为系统性能开始下降时刻,T2为系统降额运行时刻,T3为系统性能开始恢复时刻,T4为台风灾害离境时刻,T5为系统性能恢复到正常状态时刻,按照台风灾害发展过程可将电力系统状态分为灾前预防、灾时抵御及降额运行和灾后恢复四个阶段。SI(t)为不发生故障的理想系统性能曲线,S(t)为台风灾害下系统实际性能曲线,阴影区域为系统性能损失情况。本发明采用负荷供应量来描述系统性能。
本发明是通过如下技术方案来实现的:一种台风灾害下考虑连锁故障演化的电力系统弹性评估方法,主要是通过台风灾害下输电线路故障模型生成故障场景;通过短、长两时间尺度的连锁故障演化模型分别模拟因暂态稳定和热稳定问题引起的连锁故障及其演化过程;短时间指的通常是故障刚发生的机电暂态过程,几秒或十几秒达到稳态,与电力系统稳定性有关;长时间通常是0.5~1小时,本发明中取0.5小时;结合电力系统恢复模型,采用序贯蒙特卡洛模拟法,根据系统性能情况评估电力系统弹性指标值。具体包括如下步骤:
步骤1、输入电力系统数据和台风灾害预报信息,设置模拟台风轮次n=1,初始仿真时刻tk为台风开始影响电力系统时刻T0,T0为台风开始侵袭系统时刻,仿真时间间隔为Δt;
步骤2、更新tk时刻台风灾害情况和负荷水平,根据台风灾害下输电线路故障模型计算求取受灾输电线路在时间段(tk,tk+Δt)内的故障概率,通过抽样法得到输电线路故障集,将其作为tk时刻新增故障;比如:某条线路故障概率为0.1,生成一个0到1的随机数,若随机数小于0.1即认为该线路故障,若大于0.1则该线路不故障;
具体的,台风灾害下输电线路故障模型为:
(1)以台风与受灾输电线路位置关系为基础,r1为台风影响区域半径,r2为台风最大风速区域半径,vt(t)为台风移动速度,受灾输电线路故障概率随台风移动发生变化;为得到尽可能准确的输电线路故障概率,以0.5km为平均档距对各条输电线路分段,认为每一段内气象情况相同;第l条输电线路第n段Ll,n在t时刻的有效风速为:
式中:(σx(t),σy(t))为台风中心坐标;(x,y)为输电线路Ll,n的中心坐标;rl,n(t)为台风中心到输电线路Ll,n的距离;α1和α2分别为台风的两个最大风速参数;γ1和γ2为台风衰减参数;β(t)为输电线路Ll,n与台风环绕方向的夹角;
(2)台风灾害下输电线路Ll,n的有效风速和故障概率的近似关系为:
式中:pl,n(t)为单位长度输电线路的故障概率,单位为1/(50km·h);vd为输电线路设计风速;a和b为模型参数,通过统计分析历史数据获得;
(3)第l条输电线路在t时刻的每小时故障概率为:
式中:L′l,n为输电线路Ll,n的长度,单位为km;nn为输电线路l的分段数;
(4)输电线路l在时间段(tk,tk+Δt)内的故障概率为:
步骤3、判断系统中是否存在故障,若存在则根据电力系统连锁故障演化模型进行连锁故障演化分析以检测是否发生连锁故障,可得到系统性能开始下降时刻T1和系统降额运行时刻T2,根据故障情况更新系统状态,若不存在则直接转入步骤5;
为准确模拟由电稳定问题引起的连锁故障快速发展过程和由热稳定问题引起的连锁故障缓慢发展过程,采取暂态和稳态交替进行的连锁故障仿真模型模拟两种时间尺度的连锁故障演化过程。
为了检测tk时刻初始故障后系统可能发生的连锁故障及其演化过程,本发明将时间段(tk,tk+Δt)内的连锁故障演化过程划分为短、长两种时间尺度,先对短时间尺度连锁故障过程进行模拟,再模拟长时间尺度连锁故障过程,长时间尺度连锁故障过程期间对线路过载情况进行判断并采取必要的过载控制措施。
具体的,电力系统连锁故障演化模型过程如下:
(1)判断tk时刻是否有新增故障,若有则进行步骤(2),若没有则转入步骤(3),用于分析负荷变化后是否有新增的过载线路;
(2)进行短时间尺度连锁故障过程模拟过程如下:
a、在模拟由电稳定问题引起的短时间尺度连锁故障过程中,先判断故障线路是否被正确切除以得到停运线路情况,若故障线路保护因存在隐性故障而发生拒动、失灵情况,通过停运相邻线路来切除故障,并假设故障线路保护存在隐性故障发生拒动的概率为ph,其他保护均正确动作;
b、然后,将系统故障情况输入暂态仿真模块,暂态仿真模块将模拟系统运行情况,并判定保护动作情况,若系统失稳则采取紧急控制措施,考虑到的保护为发电机电压越限和频率越限的保护,紧急控制措施为低压减载和低频减载的系统稳定措施,若系统稳定则结束模拟;发电机电压和频率保护:当发电机的电压或频率偏移过大时,需将发电机切机以保护设备安全;系统低压减载和低频减载保护:当系统中发生故障或扰动引起电压或频率超出稳定运行边界时,低压减载和低频减载保护会切除部分负荷,以维护系统稳定。
(3)根据故障情况对系统状态进行判别,通过识别系统连通性得到子系统情况,并进行功率平衡调整来保证各子系统功率平衡,当出现功率不平衡时,通过发电机出力调整、切机或切负荷办法使功率平衡;
(4)计算各个子系统的潮流,将各条线路潮流和潮流正常值上限进行对比,检查是否存在过载线路,若存在则通过线路过载控制模型进行过载控制,若不存在则记录系统状态后结束模拟;
线路过载控制模型为:
通过发电机调度和负荷削减的办法控制线路过载情况,采用基于直流潮流的最优潮流模型计算最小负荷削减量,目标函数为:
约束条件为:
式中:Fl、Fmax,l分别为第i台发电机的有功下限和上限、第l条线路有功潮流和线路的有功潮流限值;η为电网运行人员对于线路过载的控制力度;Pi G、分别为第i台发电机有功出力、第j个节点处的负荷需求量;为第j个节点处的负荷削减量;为第j个节点处的负荷需求量;plf为第l条线路停运概率;xl、θfl、θtl分别为线路l的电抗、线路l首端和末端节点的电压相位角;μL,l为输电线路状态变量,μL,l=1为运行,μL,l=0为停运;μG,i为发电机状态变量,μG,i=1为运行,μG,i=0为停运;NG、ND、NL分别为发电机、负荷节点以及输电线路的集合;
(5)根据线路l基于潮流的停运概率公式计算过载线路故障概率,通过抽样法获取故障集,若集合为非空集则将该故障集作为新增故障进行短时间尺度连锁故障过程模拟,若集合为空集则记录系统状态后结束模拟;
忽略线路自身故障率,线路l基于潮流的停运概率公式为:
步骤4、然后根据电力系统故障恢复模型得到故障和负荷的恢复情况,根据故障情况更新系统状态,可得到T3,T3为系统性能开始恢复时刻,设定因紧急控制措施切除的负荷可在下个仿真时刻恢复供应,因元件停运损失的负荷在停运元件恢复后的下个仿真时刻恢复供应;设定在台风不再影响停运线路时刻起对停运线路进行修复,并假设恢复资源充裕,可同时开展对多个停运元件恢复工作;假设发电机具备快速启停功能,设定其停运后重新并网时间为4h;
输电线路非瞬时性停运后的恢复时间服从对数正态分布,选择对数正态分布作为停运线路的恢复模型,其概率密度函数为:
式中:μ和σ分别为对数正态分布的均值和方差;
步骤5、进入下一个仿真时刻(台风是一个过程,需要一段一段的仿真。举个例子:台风20h,可以按照0.5h,0.5h的时间段进行仿真),重复步骤2~4直到台风离境后电力系统性能恢复正常,可得到台风灾害离境时刻T4,系统性能恢复到正常状态时刻T5;SI(t)为不发生故障的理想系统性能曲线,S(t)为台风灾害下系统实际性能曲线,阴影区域为系统性能损失情况,采用负荷供应量来描述系统性能,第k轮台风灾害下弹性指标rk为:
式中:rk∈[0,1],rk越大表示系统应对第k轮台风灾害的能力越强;
步骤6、进行下一轮台风模拟(一个台风过程具有偶然性,需要多次进行重复试验),重复步骤2~5直到达到预设的台风模拟轮次;
步骤7、计算电力系统弹性指标值:通过对台风过程进行多次模拟的方法评估电力系统弹性,根据电力系统在第k轮台风灾害下弹性指标rk,考虑灾害持续时间的电力系统弹性指标为:
式中:NMCS为模拟台风总轮次;R∈[0,1],R越大表示系统应对台风灾害的能力越强。
与现有技术相比本发明具有以下有益效果:本发明所提供的一种台风灾害下考虑连锁故障演化的电力系统弹性评估方法,能够有效模拟台风灾害导致的初始故障情况,通过短、长两时间尺度的连锁故障演化模型分别模拟因暂态稳定和热稳定问题引起的连锁故障及其演化过程,对于故障场景的刻画更为全面和客观,更贴近于台风灾害下电力系统实际故障发展情况,从而使弹性评估结果更加准确,有利于电网运行人员准确评估电力系统应对台风灾害的能力,制定更有针对性的恢复策略。
附图说明
图1为本发明台风灾害下电力系统状态示意图。
图2为本发明台风灾害下考虑连锁故障演化的电力系统弹性评估流程。
图3为本发明台风与输电线路位置示意图。
图4为本发明短时间尺度连锁故障过程模拟。
图5为本发明长时间尺度连锁故障过程模拟。
图6为具体实施例所涉及的电力系统地理拓扑图。
图7为具体实施例所涉及的IEEE 39节点系统线路编号。
图8为具体实施例所涉及IEEE 39节点系统的阶梯式日负荷波动图。
图9为具体实施例所涉及的台风期间部分输电线路故障概率。
图10为具体实施例所涉及第k轮台风场景下系统负荷曲线。
图11为具体实施例所涉及弹性指标值与台风轮次关系图。
具体实施方式
以下结合具体实施例对本发明作进一步说明。
本实施例以IEEE39系统为例进行仿真分析,假设其位于我国华南沿海地区。电力系统地理拓扑图如图6所示,其中数字为节点编号,虚线为根据气象部门历史台风数据设定的台风路径;各节点的经纬度设定见表1,线路编号见附图7。
表1节点经纬度表
图1为台风灾害下电力系统性能变化情况,T0为台风开始侵袭系统时刻,T1为系统性能开始下降时刻,T2为系统降额运行时刻,T3为系统性能开始恢复时刻,T4为台风灾害离境时刻,T5为系统性能恢复到正常状态时刻,按照台风灾害发展过程可将电力系统状态分为灾前预防、灾时抵御及降额运行和灾后恢复四个阶段。SI(t)为不发生故障的理想系统性能曲线,S(t)为台风灾害下系统实际性能曲线,阴影区域为系统性能损失情况。本实施例采用负荷供应量来描述系统性能。
本实施例的一种台风灾害下考虑连锁故障演化的电力系统弹性评估方法,流程如图2所示,主要是通过台风灾害下输电线路故障模型生成故障场景;利用基于短、长两时间尺度的连锁故障演化模型分析电力系统可能发生的连锁故障情况,短时间指的通常是故障刚发生的机电暂态过程,几秒或十几秒达到稳态,与电力系统稳定性有关;长时间通常是0.5~1小时,本实施例中取0.5小时;结合电力系统恢复模型,采用序贯蒙特卡洛模拟法,根据系统性能情况评估电力系统弹性指标值。具体包括如下步骤:
步骤1、输入电力系统数据和台风灾害预报信息,设置模拟台风轮次n=1,初始仿真时刻tk为台风开始影响电力系统时刻T0,T0为台风开始侵袭系统时刻,仿真时间间隔为Δt;
步骤2、更新tk时刻台风灾害情况和负荷水平,根据台风灾害下输电线路故障模型计算求取受灾输电线路在时间段(tk,tk+Δt)内的故障概率,通过抽样法得到输电线路故障集,将其作为tk时刻新增故障;
具体的,台风灾害下输电线路故障模型为:
(1)以台风与受灾输电线路位置关系为基础,如图3所示,r1为台风影响区域半径,r2为台风最大风速区域半径,vt(t)为台风移动速度,受灾输电线路故障概率随台风移动发生变化;为得到尽可能准确的输电线路故障概率,以0.5km为平均档距对各条输电线路分段,认为每一段内气象情况相同;第l条输电线路第n段Ll,n在t时刻的有效风速为:
式中:(σx(t),σy(t))为台风中心坐标;(x,y)为输电线路Ll,n的中心坐标;rl,n(t)为台风中心到输电线路Ll,n的距离;α1和α2分别为台风的两个最大风速参数;γ1和γ2为台风衰减参数;β(t)为输电线路Ll,n与台风环绕方向的夹角;
(2)台风灾害下输电线路Ll,n的有效风速和故障概率的近似关系为:
式中:pl,n(t)为单位长度输电线路的故障概率,单位为1/(50km·h);vd为输电线路设计风速;a和b为模型参数,通过统计分析历史数据获得;
(3)第l条输电线路在t时刻的每小时故障概率为:
式中:L′l,n为输电线路Ll,n的长度,单位为km;nn为输电线路l的分段数;
(4)输电线路l在时间段(tk,tk+Δt)内的故障概率为:
步骤3、判断系统中是否存在故障,若存在则根据电力系统连锁故障演化模型进行连锁故障演化分析以检测是否发生连锁故障,可得到系统性能开始下降时刻T1和系统降额运行时刻T2,根据故障情况更新系统状态,若不存在则直接转入步骤5;
具体的,电力系统连锁故障演化模型过程如下,如图4所示:
(1)判断tk时刻是否有新增故障,若有则进行步骤(2),若没有则转入步骤(3),用于分析负荷变化后是否有新增的过载线路;
(2)进行短时间尺度连锁故障过程模拟具体过程如下:
a、在模拟由电稳定问题引起的短时间尺度连锁故障过程中,先判断故障线路是否被正确切除以得到停运线路情况,若故障线路保护因存在隐性故障而发生拒动、失灵情况,通过停运相邻线路来切除故障,并假设故障线路保护存在隐性故障发生拒动的概率为ph,其他保护均正确动作;
b、然后,将系统故障情况输入暂态仿真模块,暂态仿真模块将模拟系统运行情况,并判定保护动作情况,若系统失稳则采取紧急控制措施,考虑到的保护为发电机电压越限和频率越限的保护,紧急控制措施为低压减载和低频减载的系统稳定措施,若系统稳定则结束模拟;
(3)根据故障情况对系统状态进行判别,通过识别系统连通性得到子系统情况,并进行功率平衡调整来保证各子系统功率平衡,当出现功率不平衡时,通过发电机出力调整、切机或切负荷办法使功率平衡;
(4)计算各个子系统的潮流,将各条线路潮流和潮流正常值上限进行对比,检查是否存在过载线路,若存在则通过线路进行过载控制模型对线路进行过载控制,若不存在则记录系统状态后结束模拟:
线路过载控制模型为:
通过发电机调度和负荷削减的办法控制线路过载情况,采用基于直流潮流的最优潮流模型计算最小负荷削减量,目标函数为:
约束条件为:
式中:Fl、Fmax,l分别为第i台发电机的有功下限和上限、第l条线路有功潮流和线路的有功潮流限值;η为电网运行人员对于线路过载的控制力度;Pi G、分别为第i台发电机有功出力、第j个节点处的负荷需求量;为第j个节点处的负荷削减量;为第j个节点处的负荷需求量;plf为第l条线路停运概率;xl、θfl、θtl分别为线路l的电抗、线路l首端和末端节点的电压相位角;μL,l为输电线路状态变量,μL,l=1为运行,μL,l=0为停运;μG,i为发电机状态变量,μG,i=1为运行,μG,i=0为停运;NG、ND、NL分别为发电机、负荷节点以及输电线路的集合;
(5)根据线路l基于潮流的停运概率公式计算过载线路故障概率,通过抽样法获取故障集,若集合为非空集则将该故障集作为新增故障进行短时间尺度连锁故障过程模拟,若集合为空集则记录系统状态后结束模拟:
忽略线路自身故障率,线路l基于潮流的停运概率公式为:
步骤4、然后根据电力系统故障恢复模型得到故障和负荷的恢复情况,根据故障情况更新系统状态,可得到T3,T3为系统性能开始恢复时刻,设定因紧急控制措施切除的负荷可在下个仿真时刻恢复供应,因元件停运损失的负荷在停运元件恢复后的下个仿真时刻恢复供应;设定在台风不再影响停运线路时刻起对停运线路进行修复,并假设恢复资源充裕,可同时开展对多个停运元件恢复工作;假设发电机具备快速启停功能,设定其停运后重新并网时间为4h;
输电线路非瞬时性停运后的恢复时间服从对数正态分布,选择对数正态分布作为停运线路的恢复模型,其概率密度函数为:
式中:μ和σ分别为对数正态分布的均值和方差;
步骤5、进入下一个仿真时刻,重复步骤2~4直到台风离境后电力系统性能恢复正常,可得到,可得到台风灾害离境时刻T4,系统性能恢复到正常状态时刻T5;SI(t)为不发生故障的理想系统性能曲线,S(t)为台风灾害下系统实际性能曲线,阴影区域为系统性能损失情况,采用负荷供应量来描述系统性能,第k轮台风灾害下弹性指标rk为:
式中:rk∈[0,1],rk越大表示系统应对第k轮台风灾害的能力越强;
步骤6、进行下一轮台风模拟,重复步骤2~5直到达到预设的台风模拟轮次;
步骤7、计算电力系统弹性指标值:通过对台风过程进行多次模拟的方法评估电力系统弹性,根据电力系统在第k轮台风灾害下弹性指标rk,考虑灾害持续时间的电力系统弹性指标为:
式中:NMCS为模拟台风总轮次;R∈[0,1],R越大表示系统应对台风灾害的能力越强。
本实施例中,通过权衡计算精度和效率,仿真时间间隔设置为0.5h,阶梯式日负荷波动图见图8。其他参数设定见表2所示。考虑到不同停运线路恢复时间差异,将因台风导致的停运线路参数设定为μ=3.7和σ=2.4,将因连锁故障而停运线路的参数设定为μ=3.2和σ=1.4。本实施例使用电网和保护系统动态仿真模型进行暂态仿真分析,负荷采用恒阻抗模型,并使用Matpower软件进行稳态仿真分析。
表2参数设定
暂态仿真考虑到的发电机电压和频率保护,系统低压减载和低频减载保护,具体为:
本实施例设定的过电压和低电压保护动作临界值分别为1.3p.u.和0.8p.u.,动作时限均为1.0s;设定的高频和低频保护动作临界值分别为51Hz和48Hz,动作时限均为10.0s。对于每个负荷节点,本实施例设定的低压减载的动作临界值为0.87p.u.,动作时限为0.5s,削减负荷量为25%;设定的低频减载动作临界值为48Hz.,动作时限为0.5s,削减负荷量为25%。
假设台风从第一天0时开始侵袭系统,路径如图6所示,通过计算得到台风期间不同受灾线路的故障概率,受灾最严重的五条输电线路故障概率如图7所示。
进行第k轮台风模拟时,线路L41、L7和L25受台风影响分别在不同时刻发生故障。从不考虑和考虑可能发生连锁故障的角度生成故障场景,并分别设置为场景1和场景2,具体故障情况见表3和表4。
表3场景1电力系统故障情况
表4场景2电力系统故障情况
场景1与场景2的负荷变化情况如图8所示,通过分析可以看出,两个场景负荷供应情况相似,仅在时间段(18.5h,30h)期间不同,虚线所围成的面积为连锁故障造成的负荷缺额。由此可见,本实施例所提方法不仅能够模拟台风灾害下的初始故障场景,而且能够模拟后续可能发生的连锁故障演化过程,对故障场景的刻画更贴近于实际情况,验证了本实施例所提方法的可行性与客观性。
通过对多轮台风过程进行模拟,并根据传统未考虑连锁故障演化过程的弹性评估方法(以下简称传统评估方法)和本发明评估方法分别计算弹性指标值,计算结果如图9所示,弹性评估结果如表5所示。通过图9与表5可以看出,传统评估方法的弹性高于本发明评估方法的弹性评估结果,这是由于传统评估方法在刻画故障场景时,仅考虑了初始故障,没有考虑后续可能发生的连锁故障,使弹性评估结果偏于乐观,本发明评估方法的评估结果更为客观,验证了本发明所提方法的有效性。
表5电力系统弹性评估结果
本发明要求保护的范围不限于以上具体实施方式,而且对于本领域技术人员而言,本发明可以有多种变形和更改,凡在本发明的构思与原则之内所作的任何修改、改进和等同替换都应包含在本发明的保护范围之内。
Claims (2)
1.一种台风灾害下考虑连锁故障演化的电力系统弹性评估方法,其特征在于:包括如下步骤:通过台风灾害下输电线路故障模型生成故障场景;利用基于短、长两时间尺度的连锁故障演化模型模拟电力系统可能发生的连锁故障及其演化过程情况;结合电力系统恢复模型,采用序贯蒙特卡洛模拟法,根据系统性能情况评估电力系统弹性指标值。
2.根据权利要求1所述的一种台风灾害下考虑连锁故障演化的电力系统弹性评估方法,其特征在于:具体包括如下步骤:
步骤1、输入电力系统数据和台风灾害预报信息,设置模拟台风轮次n=1,初始仿真时刻tk为台风开始影响电力系统时刻T0,T0为台风开始侵袭系统时刻,仿真时间间隔为Δt;
步骤2、更新tk时刻台风灾害情况和负荷水平,根据台风灾害下输电线路故障模型计算求取受灾输电线路在时间段(tk,tk+Δt)内的故障概率,通过抽样法得到输电线路故障集,将其作为tk时刻新增故障;
具体的,台风灾害下输电线路故障模型为:
(1)以台风与受灾输电线路位置关系为基础,r1为台风影响区域半径,r2为台风最大风速区域半径,vt(t)为台风移动速度,受灾输电线路故障概率随台风移动发生变化;0.5km为平均档距对各条输电线路分段,第l条输电线路第n段Ll,n在t时刻的有效风速为:
式中:(σx(t),σy(t))为台风中心坐标;(x,y)为输电线路Ll,n的中心坐标;rl,n(t)为台风中心到输电线路Ll,n的距离;α1和α2分别为台风的两个最大风速参数;γ1和γ2为台风衰减参数;β(t)为输电线路Ll,n与台风环绕方向的夹角;
(2)台风灾害下输电线路Ll,n的有效风速和故障概率的近似关系为:
式中:pl,n(t)为单位长度输电线路的故障概率,单位为1/(50km·h);vd为输电线路设计风速;a和b为模型参数,通过统计分析历史数据获得;
(3)第l条输电线路在t时刻的每小时故障概率为:
式中:Ll′,n为输电线路Ll,n的长度,单位为km;nn为输电线路l的分段数;
(4)输电线路l在时间段(tk,tk+Δt)内的故障概率为:
步骤3、判断系统中是否存在故障,若存在则根据电力系统连锁故障演化模型进行连锁故障演化分析以检测是否发生连锁故障,可得到系统性能开始下降时刻T1,系统降额运行时刻T2,根据故障情况更新系统状态,若不存在则直接转入步骤5;
具体的,电力系统连锁故障演化模型过程如下:
(1)判断tk时刻是否有新增故障,若有则进行步骤(2),若没有则转入步骤(3),用于分析负荷变化后是否有新增的过载线路;
(2)进行短时间尺度连锁故障过程模拟具体过程如下:
a、在模拟由电稳定问题引起的短时间尺度连锁故障过程中,先判断故障线路是否被正确切除以得到停运线路情况,若故障线路保护因存在隐性故障而发生拒动、失灵情况,通过停运相邻线路来切除故障,并假设故障线路保护存在隐性故障发生拒动的概率为ph,其他保护均正确动作;
b、然后,将系统故障情况输入暂态仿真模块,暂态仿真模块将模拟系统运行情况,并判定保护动作情况,若系统失稳则采取紧急控制措施,考虑到的保护为发电机电压越限和频率越限的保护,紧急控制措施为低压减载和低频减载的系统稳定措施,若系统稳定则结束模拟;
(3)根据故障情况对系统状态进行判别,通过识别系统连通性得到子系统情况,并进行功率平衡调整来保证各子系统功率平衡,当出现功率不平衡时,通过发电机出力调整、切机或切负荷办法使功率平衡;
(4)计算各个子系统的潮流,将各条线路潮流和潮流正常值上限进行对比,检查是否存在过载线路,若存在则通过线路过载控制模型进行过载控制,若不存在则记录系统状态后结束模拟;
线路过载控制模型为:
通过发电机调度和负荷削减的办法控制线路过载情况,采用基于直流潮流的最优潮流模型计算最小负荷削减量,目标函数为:
约束条件为:
式中:Fl、Fmax,l分别为第i台发电机的有功下限和上限、第l条线路有功潮流和线路的有功潮流限值;η为电网运行人员对于线路过载的控制力度;Pi G、分别为第i台发电机有功出力、第j个节点处的负荷需求量;为第j个节点处的负荷削减量;Pj D为第j个节点处的负荷需求量;plf为第l条线路停运概率;xl、θfl、θtl分别为线路l的电抗、线路l首端和末端节点的电压相位角;μL,l为输电线路状态变量,μL,l=1为运行,μL,l=0为停运;μG,i为发电机状态变量,μG,i=1为运行,μG,i=0为停运;NG、ND、NL分别为发电机、负荷节点以及输电线路的集合;
(5)根据线路l基于潮流的停运概率公式计算过载线路故障概率,通过抽样法获取故障集,若集合为非空集则将该故障集作为新增故障进行短时间尺度连锁故障过程模拟,若集合为空集则记录系统状态后结束模拟;
忽略线路自身故障率,线路l基于潮流的停运概率公式为:
步骤4、然后根据电力系统故障恢复模型得到故障和负荷的恢复情况,根据故障情况更新系统状态,可得到T3,T3为系统性能开始恢复时刻,设定因紧急控制措施切除的负荷可在下个仿真时刻恢复供应,因元件停运损失的负荷在停运元件恢复后的下个仿真时刻恢复供应;设定在台风不再影响停运线路时刻起对停运线路进行修复,并假设恢复资源充裕,可同时开展对多个停运元件恢复工作;假设发电机具备快速启停功能,设定其停运后重新并网时间为4h;
输电线路非瞬时性停运后的恢复时间服从对数正态分布,选择对数正态分布作为停运线路的恢复模型,其概率密度函数为:
式中:μ和σ分别为对数正态分布的均值和方差;
步骤5、进入下一个仿真时刻,重复步骤2~4直到台风离境后电力系统性能恢复正常,可得到台风灾害离境时刻T4,系统性能恢复到正常状态时刻T5;SI(t)为不发生故障的理想系统性能曲线,S(t)为台风灾害下系统实际性能曲线,阴影区域为系统性能损失情况,采用负荷供应量来描述系统性能,第k轮台风灾害下弹性指标rk为:
式中:rk∈[0,1],rk越大表示系统应对第k轮台风灾害的能力越强;
步骤6、进行下一轮台风模拟,重复步骤2~5直到达到预设的台风模拟轮次;
步骤7、计算电力系统弹性指标值:通过对台风过程进行多次模拟的方法评估电力系统弹性,根据电力系统在第k轮台风灾害下弹性指标rk,考虑灾害持续时间的电力系统弹性指标为:
式中:NMCS为模拟台风总轮次;R∈[0,1],R越大表示系统应对台风灾害的能力越强。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210397176.XA CN114840982A (zh) | 2022-04-15 | 2022-04-15 | 台风灾害下考虑连锁故障演化的电力系统弹性评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210397176.XA CN114840982A (zh) | 2022-04-15 | 2022-04-15 | 台风灾害下考虑连锁故障演化的电力系统弹性评估方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114840982A true CN114840982A (zh) | 2022-08-02 |
Family
ID=82566665
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210397176.XA Pending CN114840982A (zh) | 2022-04-15 | 2022-04-15 | 台风灾害下考虑连锁故障演化的电力系统弹性评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114840982A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115374654A (zh) * | 2022-10-24 | 2022-11-22 | 国网四川省电力公司电力科学研究院 | 一种评估电力系统弹性的方法和系统 |
CN117422438A (zh) * | 2023-10-07 | 2024-01-19 | 国家电网有限公司华东分部 | 一种输电线路加固方案的确定方法及装置 |
-
2022
- 2022-04-15 CN CN202210397176.XA patent/CN114840982A/zh active Pending
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115374654A (zh) * | 2022-10-24 | 2022-11-22 | 国网四川省电力公司电力科学研究院 | 一种评估电力系统弹性的方法和系统 |
CN117422438A (zh) * | 2023-10-07 | 2024-01-19 | 国家电网有限公司华东分部 | 一种输电线路加固方案的确定方法及装置 |
CN117422438B (zh) * | 2023-10-07 | 2024-03-29 | 国家电网有限公司华东分部 | 一种输电线路加固方案的确定方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109103874B (zh) | 考虑负荷特性与分布式电源接入的配电网可靠性评估方法 | |
CN114840982A (zh) | 台风灾害下考虑连锁故障演化的电力系统弹性评估方法 | |
US9727036B2 (en) | Operation plan creating method, computer product, and operation plan creating apparatus | |
Dong et al. | Frequency prediction of power systems in FNET based on state-space approach and uncertain basis functions | |
CN102623990B (zh) | 自适应电网运行方式的低频减载参数在线校核方法 | |
CN109447441B (zh) | 一种考虑新能源机组不确定性的暂态稳定风险评估方法 | |
EP2381094A1 (en) | Energy network and control thereof | |
CN109034461A (zh) | 一种基于实际电网监测信息的电压暂降随机预估方法 | |
CN113077075B (zh) | 新能源不确定性电力系统安全风险预防控制方法及装置 | |
CN104574202A (zh) | 一种风电场出力特性分析方法 | |
CN103246939A (zh) | 基于安全稳定裕度的电网运行安全风险事件在线辨识方法 | |
CN117557065B (zh) | 一种基于bim技术的建筑工程施工进度监管系统 | |
CN111209535B (zh) | 一种电力设备相继故障风险识别方法及系统 | |
CN113609637B (zh) | 一种考虑故障连锁的多灾害配电网弹性评估方法 | |
CN114418138A (zh) | 一种多设备联合自检的智能电网运维方法及系统 | |
JP7074932B1 (ja) | エネルギー運用システム、エネルギー運用方法、および記憶媒体 | |
CN116502771B (zh) | 一种基于电力物资预测的配电方法及系统 | |
Shokrollahi et al. | Reliability assessment of distribution system using fuzzy logic for modelling of transformer and line uncertainties | |
Liu | Short-term operational reliability evaluation for power systems under extreme weather conditions | |
CN111146779B (zh) | 一种电力设备相继故障的大电网柔性安全控制方法及系统 | |
Begovic et al. | Four-year performance assessment of the 342 kW PV system at Georgia Tech | |
CN110365047B (zh) | 一种含分布式光伏发电系统的电网短路电流概率评估方法 | |
CN112670991A (zh) | 一种基于敏感用户暂降损失的电网防雷治理优化方法 | |
Wei et al. | Resilience Assessment Using Simulation System of Distribution Network under Extreme Weather | |
Abdelaziz et al. | Future Climate Change Impact On Wind Gust For Offshore Wind Turbines Cut-Out Within The UK Exclusive Economic Zone |
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 |