电力系统风险评估

1、电力系统风险评估,绪论,什么叫风险？ 电力系统中存在哪些风险？ 大规模风电入网给电力系统带来哪些风险？ 教材及参考文献： 电力系统风险评估模型、方法和应用李文沅著 郭永基电力系统可靠性分析 静铁岩，大规模风电并网条件下的电力系统有功功率平衡理论研究，博士论文第一章绪论部分。 张丽英，叶廷路等，大规模风电接入电网的相关问题及措施. 中国电机工程学报. Vol.30(25)，2010,1.1 电力系统中的风险,风险和可靠性的关系 描述同一个事实的两个方面：更高的风险即更低的可靠性 本课程的目的：讨论电力系统风险评估的模型、方法及其工程应用 电力系统风险的根源：设备的随机故障，负荷的波动，市场环境下

2、多变的市场需求及市场参与者的行为，可间歇性能源入网。,风险管理包括： 实施风险的定量评估：建立表征系统风险的指标 确定降低风险的措施：意味着投资 确认可接受的风险水平：对技术、经济、社会以及环境等因素进行综合评估,电力系统风险评估适用的领域： 发、输、配系统的规划 电源规划 运行风险评估 设备维修 设备备用分析 变电站接线分析 输电服务 。,1.2 电力系统风险评估的基本概念,风险 = 事件发生的概率 X 事件发生带来的损失 系统风险评估步骤 确定元件停运模型（第二章） 选择系统状态，计算发生的概率（两种方法） 评估所选择状态的后果（可能是简单的功率平衡分析，网络结构连通性分析，优化潮流，电压

3、稳定分析） 计算风险指标（指标的概率分布，随机变量的期望值）,元件停运模型参数 由元件停运统计数据计算得到； 数据的预处理问题 数据的动态特性,评估所选系统状态的后果 由评估分类决定具体的运算 充裕性分析：表明系统设施是否能充分满足用户的负荷需求和系统运行的约束条件，只涉及到系统的稳态条件，不要求动态和暂态分析。 安全性分析：系统对动态和暂态扰动的响应能力。因此通常要进行动态、暂态或电压稳定性分析。 基于历史停运统计数据的评估,风险指标 反应停电规模，电压的变化或均值，潮流的变化或均值等 单位停电损失 基于用户损失函数的方法； 基于投资核算的方法； 基于国民生产总值的方法；,第二章 系统元件的

4、停运模型,可修复强迫失效,Component two-state model:,Failure rate:,Repair rate:,Failure rate and repair rate: from statistics of the target component. For example:,Availability A:,Mean Time To Failure (MTTF):,Mean Time To Repair (MTTR):,Unavailability U:,Failure frequency f :,Mean Time Between Failure (MTBF): MT

5、BF = MTTF+MTTR,Example:,Statistics of a 200MW generator:,The reliability of a system: depends on the reliability of the components involved, and systems topology Series connected system:,Example:,A simple power supply system: the availability of each component is: generator A1=0.990099, transformer

6、A2=0.999933, busbar A3=0.999965, circuit breaker A4=0.999833, transmission line A5=0.999334. what is the availability of the whole system?,Parallel connected system: only when all components fail, the system fails.,(a) Parallel connected system,(b) Equivalent system,Example:,Two transformers are in

7、parallel operation. Statistics are:,What is the failure rate of this parallel connected system and the MTTF of the system?,Integrated system: break down into series connected subsystem and parallel connected subsystem,Example,A: power supply bus; C: load point. Suppose busbar A is always available.

8、Question: what is the failure rate of point C and the total down time; What is the reliability of point c,Parallel connected 1,series connected 2,Parallel connected 3,Parallel connected 1: L1 / L2,Series connected 2: L12 +L3,Parallel connected 3: L123 / L4,Unavailability: Total down time at busbar C

9、: Reliability at busbar C:,作业： Matlab编程，实现可修复强迫失效相关计算。,计划停运,最简单的蒙特卡洛仿真流程,随机生产模拟：元件的投入状态，负荷的随机变化，间歇性电源出力变化,考虑连锁故障的Monte Carlo 仿真流程,对系统进行抽样,元件的失效模型参考文献,考虑保护的隐性故障 J.Chen, J.S.Thorp, I.Dobson. Cascading dynamics and mitigation assessment in power system disturbances via a hidden failure model. Internati

10、onal Journal of Electrical Power and Energy Systems. 2005, 27: 318-326 考虑天气变化对导体散热影响 M.Anghel, K.A. Werley, A.E. Motter. Stochastic model for power grid dynamics. The 40th HICSS. 2007 环境相依失效 卢锦玲，朱永利. 基于暂态能量裕度的电力系统脆弱性评估. 电工技术学报，2010，25(6)：96-103 （看2.1节）,保护的隐性故障,当有线路断开时，与此线路两端相连的所有线路的保护均存在发生隐性故障的可能。 仅

11、考虑单重隐性故障。,考虑天气变化对导体散热影响,随机线路故障模型 线路l上发生的随机故障由泊松过程描述。即在时间段t中发生的故障的次数服从泊松分布，均值为,风功率波动建模,Wind Turbine Generator Modeling,P(SWt): wind turbine output SWt : Wind speed Vci: cut-in wind speed Vr : rated speed Vco: cut-out speed,P-W characteristic established for 1.5-MW doubly-fed (continuous) and 1.5-MW f

12、ixed speed asynchronous wind generators (dotted).,Wind speed modeling: Drawing of an uniformly distributed number u on the interval 0,1; Application of that drawn random number u to the Weibull Cumulative Distribution Function F(W,A,B) in order to determine the associated wind speed W.,Note: Cumulat

13、ive Distribution Function: p ( x X ),序贯蒙特卡洛仿真,状态时间抽样法 第1步：指定所有元件的初始状态（一般假定均处于运行状态） 第2步：对每一元件停留在当前状态的持续时间进行抽样。 第3步：在所研究的时间跨度内重复第2步，并记录所有元件的每一状态持续时间的抽样值，获得给定时间跨度内每一元件的随机状态转移过程。 第4步：组合所有元件的状态转移过程，建立系统时序状态转移循环过程。 第5步：对每一个不同系统状态进行系统分析，计算风险指标函数。,负荷的处理,负荷保持不变 年负荷曲线,年负荷曲线，并考虑随机变化,负荷分区、增长及随机波动 假设在蒙特卡罗仿真的每一个循

14、环内，负荷为该循环的平均负荷与某一随机变化量的叠加，且在该循环中保持不变。所有负荷节点按其所处位置划分为NF个区域，处于同一区域的负荷变化幅度相同，不同区域之间负荷的变化幅度不同。 代表平均负荷的增长速度，略大于1，可设定。随机因子r在区间 上变化。则对于给定的 ，所引入的负荷波动系数为 。,第t+1循环的平均负荷,缓慢的负荷增长,P98 负荷曲线模型,多级负荷水平模型 将各个母线符合按其遵循的不同负荷曲线分类为相同的母线组。每个母线组有其使用的负荷持续曲线。 使用聚类技术 选择聚类均值Mij的初值，i代表聚类，j代表曲线j 利用下式计算每小时负荷点至每个聚类均值的欧拉距离,负荷点分配到最近的

15、聚类，重新编组，并利用下式计算新的聚类均值 重复第2步和第3步，直到全部聚类均值在迭代中保持不变为止。 使用收敛后的聚类均值Mij作为多级负荷模型中每一曲线每一聚类的负荷水平。,故障分析,基于交流潮流的最优潮流模型 目标函数：社会效益最大化（等效成发电费用最小），切负荷量最小，网损最小。 约束条件： 等式约束：功率平衡 不等式约束：运行极限,期望却供电量 (EENS),NL : various load levels Ti : the no. of hours that system has the ith load level C(s) : load shed of state s n(s)

16、 : the no. of runs that system has a load shed of C(s) Ni : the total no. of runs in the simulation Fi : under the ith load level, the set of runs that have non-zero load shed.,EENS： 其中 S 导致切负荷的系统状态集合, Ci 是状态 i下，系统的切负荷量, ti 是系统状态i的蒙特卡洛循环次数, T是蒙特卡洛总的循环次数.,负荷削减频率,n(s)：抽样中S状态的发生数； Fi是多级负荷模型中，第i级负荷水平下系统

17、失效状态的集合； Ti是第i负荷水平的时间长度h； T是负荷曲线的总时间区间h； NL是负荷水平分级数,Average line flow limit PerMW served : it is used to show the efficiency of the investment. where Pd(t) is the average load demand of the test system at the tth simulation run. NLine is the no. of transmission lines in the system. Fjmax(t) is the t

18、ransmission capacity limit of line j at the tth simulation run. The average line flow limit per MW served is used as a measure of grid utilization. Since the grid investment is related to the maximum line flow limits and the societal benefit is related to the power served, the average lineflow limit

19、 per MW served is one way to indicate the ratio of societal benefit to the grid investment.,参考文献 R. Billionton, and W. Wangdee, “Reliability-based Transmission Reinforcement Planning Associated with Large-scale Wind Farms”, IEEE Trans. on Power Systems, Vol.22, No.1, Feb. 2007. 6.4例1 6.5例2,指标简介,Basi

20、c adequacy indices LOLP: 失负荷概率 ( of an electric power system), the probability of all the outage events leaving the system with an available capacity lower than expected load Lo.,LOLE: the loss of load expectation. Widely accepted LOLE risk criterion: “one day in ten years”,Ok: the magnitude of the

21、k-th outage in the system pk: the probability of a capacity outage of magnitude Ok. tk: the number of days that an outage of magnitude Ok would cause a loss of load in the system,daily peak load curve,Loss of Energy Expectation: the ratio of expected non-served energy to total energy demand over a period of time.,Ek is the energy not supplied due to a capacity outage Ok; E is the total energy during the period of time.,期望负荷削减频率,Probability of load curtailment PLC P(s): probability of state s Fi: the set of all states under load level i with non-zero load shed. T: total hours of the sim