蒙特卡洛法在电力系统可靠性评估中的应用

《蒙特卡洛法在电力系统可靠性评估中的应用》由会员分享，可在线阅读，更多相关《蒙特卡洛法在电力系统可靠性评估中的应用（25页珍藏版）》请在金锄头文库上搜索。

1、蒙特卡洛法在电力系统可靠性评估中的应用3.1电力系统可靠T生评估的容与意义可靠性指的是处于某种运行条件下的元件、设备或系统在规定时间完成预定功能的概率。电力系统可 靠性是才旨电网在各种运行条件下，向用户持续提供符合一定质量要求的电能的能 力。电力系统可靠性包括 充裕度（Adequacy）和安全性（security）两个方面。充裕度是指在考虑电力元件计划与非计划停运以及 负荷波动的静态条件下，电力系统维持连续供应电能的能力，因此又被称为静态可靠性。安全性指的是电力 系统能够承受如突然短路或未预料的失去元件等事件引起的扰动并不间断供应电能的能力，安全性又被称为 动态可靠性。目前国外学者对充裕度评估

2、的算法和应用矢注较多，且在理论和实践中取得了大量的研究成果， 但随着研究的深入也出现了很多函待解决的新课题。电力系统的安全性评估以系统暂态稳定性的概率分析为 基础，在原理、建模、算法和应用等方面都处于起步和探索阶段。由于电力系统的规模很大，通常根据功能 特点将其分为不同层次的子系统，如发电、输电、发输电组合、配电等子系统，对电力系统的可靠性评估通常 也是对上述子系统单独进行。不同层次的子系统的可靠性评估的任务、模型、算法都有较大区别。电力系统在 正常运行情况下，系统能够正常供电，不会出现切负荷的事件。如果系统受到某些偶发事件的扰动，如元件停 运（包括机组、线路、变压器等电力元件的计划停运与故障

3、停运）、负荷水平变化等，可能会引起系统功率失 衡、线路潮流越限和节点电压越限等故障状态，进而导致切负荷。电力系统可箕性研究的主要容是基于系统 偶发故障的概率分布及其后果分析，对系统持续供电能力进行快速和准确的评价，并找出影响系统可靠性水 平的薄弱环节以寻求改善可靠性水平的措施，为电力系统规划和运行提供决策支持。一 七 AJZTV匀日巨I工 l/TILJLIll空平/J/Za电力系统可靠性评估方法可分为确定性方法和概率性方法两类。确定性方法主要是对 几种确定的 运行方式和故障状态进行分析，校验系统的可靠f生水平。在电源规划中，典型的确定性的可靠性判据有百 分备用指标和最大机组备用指标；电网规划中

4、，确定性的可靠性判据主要是校验负荷的最小供电回路数。电 力系统是一个具有随机特T生的系统，负荷水平的波动、元件故障等都具有随机性，确定性方法难以考虑各 种状态的概率分布特性，评估结果存在较大偏 差，因此概率性方法在电力系统的可靠性评估领域得到更加广 泛应用，并在理论和实践方面取得很大的进展。概率性可靠性评估方法主要有解析法和模拟法两大类，后者 一般又被称作蒙特 卡洛法。两者的共同点是都以系统随机状态发生的概率对随机状态的后果（切负荷功率） 进行 加权，即不仅 考虑故障的严重性，同时考虑其概率性，且对随机状态的分析方法是一致的。两 者的根 本区别在于获取系统随机状态及其概率值的方法不同，解析法通

5、过故障枚举来获得系 统随机状态，通过解析 计算获得系统随机状态发生的概率；蒙特卡洛法通过随机抽样的方法获得系统随机状态，采用统计的方法以 随机状态的频率来估算概率。解析法的数学模型精确，得到的可靠性指标计算精度高，但该方法的缺点也非常突出。首先，采用解析 法要分析的系统状态数目随着系统元件数目的增长呈指数规律增长，因此难以应 用于大规模电力系统可靠性 评估的场合。采用忽略多重故障状态的“故障筛选技术”来解决这一问题，但显然会在一定程度上削弱解析 法在计算精度方 面的优势。其次，采用解析法难以获 得频率和持续时间指标，而这些又是非常重要的 可靠 性信息。最后，解析法难以处理系统中随 机因素的影响

6、，如负荷的波动、水库水位的变化等，也不易模拟运 行人员对系统的控制措施及其后果，因此影响到了计算结果的可信度。由于解析法存在上述难以克服的缺点， 在大型电力系统可靠性评估的场 合应用较少，而蒙特卡洛法则得到了广泛的应用。蒙特卡洛方法（又被称作统计试验方法）或随机抽样技术，其提出可以追溯到19世纪末期，20世纪40年代中期之后随着科学技术的发展和电子计算机的发明，该方法得到了快速的发展和应用。 几十年来，随着计算技术的迅速发展，蒙特卡洛方法的应用围日趋广阔。目前它分支。采用蒙特卡洛方法评估电力系统可靠性，存在着明显的优势。第一，在一定的精度要求下，蒙特卡洛 方法的抽样次数与系统的规模无矢，因此特

7、别适用于大型电力系统 的评估计算。第二，采用蒙特卡洛方法评 估可靠性，不但能够获得概率性指标，而且能够得到 频率和持续时间指标，得到的可靠性信息更加丰富、实 用。第三，基于蒙特卡洛方法的程序数学模型相对简单，且容易模拟负荷变化等随机因素和系统的校正控制 措施，因此计算结果更加 符合H程实际。电力系统规模日趋扩大、元件众多、控制策略复杂，因此蒙特卡洛 法在其可靠性评估中获得了日益广泛的应用。3.3蒙塔卡洛法的基本容3.3.1基本参数介绍电力系统元件众多，在可靠性评估中可根据计算需要对发电机组、输电线路、变压器、电抗器、电容 器、保护元件、自动重合闸装置、母线等可修复元件进行状态模拟。假定某可修复

8、元件的故障率和修复率分 别为兄、气平均无故障工作时间和平均维修时间分别为M7TF、MTTR，则存在以下重要矢系式I。）MTTR=1/jl /(3.2)可修复强迫失效可以通过“运行停运运行的循坏过程来模拟，如图一所示图3可修复元件运行和停运循环过程平均不可用率，其数学形式可由下列三个定义之一来表达：U 2 +/*伽 mQ.3） MTTF + MTTR兄为失效率（失效次数/年）；为修复率（修复次数/年）；MTTR为平均修复时间（小时）:mttf为失效前平均时间（小时）；f为平均失效频率（失效次数/年）二MTTF/8760及r = MTTR/8760则d和r是以年为单位计的MTTF和MTTR。A、是

9、蒙特卡洛算法中模拟元件持续时间与状态转移特性的基本参数。其反映f.停延状态 .状毒图32可修复元件状态空间图的元件状态转移特性如图3-2所示，其数值可通过对元件长期运行的寿命过程和随机状态信息统计得到03.3.2非序贯蒙特卡洛模拟法非序贯蒙特卡洛模拟法常常被称为状态抽样法，它被广泛用在电力系统风险评估中。这个方法的依据 是：一个系统状态是所有元件状态的组合，且每一元件状态可由对元件出现在该状态的概率进行抽样来确 定。每一元件可用一个在, 1区间的均匀分布来模拟。假设每一元件有失效和工作两个状态，且元件 失效是相互独立的。令S代表元件；的状态，Q/代表其失效概率，则对元件，产生一个在0, 口区间

10、均匀 分布的随机数使 (工作状态)ifRi QiI (失效状(Z、OV/?zV &具有W个元件的系统状态由矢量s表示：(3.5)一个系统状态在抽样中被选定后，即进行系统分析以判新其是否是失效状态，如果是，则对该状态的 风险指标函数进行估计。当抽样的数量足够大时，系统状态S的抽样频率可作为其申既率的无偏估计 即P (沪者0.6)式中：“是抽样数；来S)是在抽样中系统状态S出现的次数。当每一个系统状态的概率通过抽样估计以后，就可计算系统失效概率、系统失效频率、系统失效平均 持续时间、以及系统其它风险指标。非序贯蒙特卡洛法和状态枚举法之间明显的区别在于：如何选择系统状态和如何计算单个系统状态的 概率

11、。在实际应用中，应注意以下几个方面：(1) 必要的一步是产生每一个元件的随机数序列，这些随机数必须满足三个基本条件：均匀性、独立性和足够长的重复周期0(2) 蒙特卡洛法是一个波动收敛过程，因此估计出的风险指标总是有一个相应的置信围。不能保证增加少量的样本就一定会减少误差，但置信围确实会随样本数的增加而变窄。(3) 适当的收敛判据是确保蒙特卡洛模拟法精度的矢键之一。方差系数常被用作为终止抽样的判据。在 电力系统风险评估中，不同的风险指标有不同的收敛速度。已经发现，期望缺供电t(EENS)指标的方差系数 收敛速率最低，因此应作为多个指标研究时的收敛判据。另一种方法是用预定的最大抽样数作为终止抽样的

12、判 据。当模拟过程结束时，校验方差系数是否足够小，如果否，则需要增加样本数再进行新的抽样。当用户并不 知道需要用多少计算时间才能达到足够小的方差系数时，可用这种替代方法。(4) 非序贯模拟过程仅需要元件的失效概率作为抽样过程的输入数据。这个特点使我们能够容易地同时 模拟可修复和老化失效引起的不可用率。为元件建立两个独立的随机数，一个是模拟可修复失效引起的不可 用率，另一个是模拟老化失效引起的不可用率。(5) 状态抽样的理念不仅适用于元件失效事件，而且也可推广应用到电力系统风险评估中其它参数的状 态抽样，例如：负荷水平、水文和气候状态等。而且，这个方法并不局限于年度为基础的模拟，还可很方便 地用

13、于进行任意时间段(周、月、季或年度)的模拟。(6) 较之状态枚举法，状态抽样法更适用于规模较大的系统或具有较高元件失效概率的系统评估。在这 些情况下，为获得相同的精度，状态枚举法需要大得多的CPU时间。(7) 与状态枚举法相似，非序贯蒙特卡洛模拟法不能计及时间相矢事件的时序信心，心1山I寸山口亦勿人X乂必人寸衣人匕KLSLJLd。3- 3.3序贯蒙特卡洛模拟法序贯蒙特卡洛法是按照时序，在一个时间跨度上进行的模拟。其中对建立虚拟系统状态转移循环过程 有不同的方法。最通用的是在这里讨论的所谓状态持续时间抽样法。状态持续时间抽样法是基于对元件状态持续时间的概率分布进行抽样，它分为以下几步：第1步：指

14、定所有元件的初始状态，通常是假设所有元件开始处于运行状态。第2步：对每一元件停留在当前状态的持续时间进行抽样。应当设定状态持续时间的概率分布。对不同 的状态，如运行或修复过程，可以假设有不同的状态持续时间概率分布。例如 下式给出指数分布的状态持 续时间的抽样值：WRj(3.7)式中，人是对应于第i个元件在】0, 口区间均匀分布的随机数。如果当前的状态是运 行状态，则儿 是第i个元件的失效率；而如果当前的状态是停运状态，则儿是第i个元件的修复率。服从不同概率分布的随 机变量的产生方法是不同的。相矢书籍有详细论述。第3步：在所研究的时间跨度(大量的抽样年)重复第2步，并记录所有元件的每一状态持续时

15、间的抽样值，即可获得给定时间跨度每一元件的时序状态转移过程，1 0 101如囹JJ所示.图33 元件时序状态转移过程第4步：组合所有元件的状态转移过程，以建立系统时序状态转移循环过程，如图34所示.10 11 01 11 10 11 01 00 01 11 亲删 S 合 I II状态仿具图34系统时序状态转移过程第5步：通过对每一个不同系统状态的系统分析，计算风险指标函数由于系统失效状态的发生 它们 的持续时间、以及后果都能被清楚地确定并记录在系统状态转移循坏过程中，因此系统风险指标的计算简单 直观。式38至3. 10是三个风险指标的通用公式。(3.8)(3.9)Ddk(3.10)D, =a -式中：巧和0,分别为系统失效概率、频率和平均持续时间；仇是第k个停运状态的持续时间；化是第）个运行状态的持续时间；女任，和冷分别为在模拟时间跨度系统失效和运行状态出现的次数。除 非失效或运行状态在抽样跨度末被截尾，否则这两个被抽取的状态数一般是相同的。可见，序贯蒙特卡洛法的矢键在于系统状态转移过程的生成，一旦这一步骤完成， 指标计算则较简单。该方法的本质是建立一个虚拟的系统运行和失效的转移循环过程。在实际应用中，重要的是认识以下几点：(1) 序贯蒙特卡洛法中至矢重要的一步是计算服从某个概率分布的状态持续时间随机变量的抽样值， 其基础是在Q 1区间均匀分布的随机数的生成。