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

3 蒙特卡洛法在电力系统可靠性评估中的应用3.1电力系统可靠性评估的内容与意义可靠性指的是处于某种运行条件下的元件、设备或系统在规定时间内完成预定功能的概率。

电力系统可靠性是指电网在各种运行条件下，向用户持续提供符合一定质量要求的电能的能力。

电力系统可靠性包括充裕度(Adequacy)和安全性(seeurity)两个方面。

充裕度是指在考虑电力元件计划与非计划停运以及负荷波动的静态条件下，电力系统维持连续供应电能的能力，因此又被称为静态可靠性。

安全性指的是电力系统能够承受如突然短路或未预料的失去元件等事件引起的扰动并不间断供应电能的能力，安全性又被称为动态可靠性。

目前国内外学者对充裕度评估的算法和应用关注较多，且在理论和实践中取得了大量的研究成果，但随着研究的深入也出现了很多函待解决的新课题。

电力系统的安全性评估以系统暂态稳定性的概率分析为基础，在原理、建模、算法和应用等方面都处于起步和探索阶段。

由于电力系统的规模很大，通常根据功能特点将其分为不同层次的子系统，如发电、输电、发输电组合、配电等子系统，对电力系统的可靠性评估通常也是对上述子系统单独进行。

不同层次的子系统的可靠性评估的任务、模型、算法都有较大区别。

电力系统在正常运行情况下，系统能够正常供电，不会出现切负荷的事件。

如果系统受到某些偶发事件的扰动，如元件停运(包括机组、线路、变压器等电力元件的计划停运与故障停运)、负荷水平变化等，可能会引起系统功率失衡、线路潮流越限和节点电压越限等故障状态，进而导致切负荷。

电力系统可靠性研究的主要内容是基于系统偶发故障的概率分布及其后果分析，对系统持续供电能力进行快速和准确的评价，并找出影响系统可靠性水平的薄弱环节以寻求改善可靠性水平的措施，为电力系统规划和运行提供决策支持。

3.2电力系统可靠性评估的基本方法电力系统可靠性评估方法可分为确定性方法和概率性方法两类。

确定性方法主要是对几种确定的运行方式和故障状态进行分析，校验系统的可靠性水平。

在电源规划中，典型的确定性的可靠性判据有百分备用指标和最大机组备用指标;电网规划中，确定性的可靠性判据主要是校验负荷的最小供电回路数。

电力系统是一个具有随机特性的系统，负荷水平的波动、元件故障等都具有随机性，确定性方法难以考虑各种状态的概率分布特性，评估结果存在较大偏差，因此概率性方法在电力系统的可靠性评估领域得到更加广泛应用，并在理论和实践方面取得很大的进展。

概率性可靠性评估方法主要有解析法和模拟法两大类，后者一般又被称作蒙特卡洛法。

两者的共同点是都以系统随机状态发生的概率对随机状态的后果(切负荷功率)进行加权，即不仅考虑故障的严重性，同时考虑其概率性，且对随机状态的分析方法是一致的。

两者的根本区别在于获取系统随机状态及其概率值的方法不同，解析法通过故障枚举来获得系统随机状态，通过解析计算获得系统随机状态发生的概率;蒙特卡洛法通过随机抽样的方法获得系统随机状态，采用统计的方法以随机状态的频率来估算概率。

解析法的数学模型精确，得到的可靠性指标计算精度高，但该方法的缺点也非常突出。

首先，采用解析法要分析的系统状态数目随着系统元件数目的增长呈指数规律增长，因此难以应用于大规模电力系统可靠性评估的场合。

采用忽略多重故障状态的“故障筛选技术”来解决这一问题，但显然会在一定程度上削弱解析法在计算精度方面的优势。

其次，采用解析法难以获得频率和持续时间指标，而这些又是非常重要的可靠性信息。

最后，解析法难以处理系统中随机因素的影响，如负荷的波动、水库水位的变化等，也不易模拟运行人员对系统的控制措施及其后果，因此影响到了计算结果的可信度。

由于解析法存在上述难以克服的缺点，在大型电力系统可靠性评估的场合应用较少，而蒙特卡洛法则得到了广泛的应用。

蒙特卡洛方法（又被称作统计试验方法)或随机抽样技术，其提出可以追溯到19世纪末期，20世纪40年代中期之后随着科学技术的发展和电子计算机的发明，该方法得到了快速的发展和应用。

几十年来，随着计算技术的迅速发展，蒙特卡洛方法的应用范围日趋广阔。

目前它已经被广泛应用到包括电力系统可靠性分析在内的各类科学研究与工程设计领域中，成为计算数学的一个重要分支。

采用蒙特卡洛方法评估电力系统可靠性，存在着明显的优势。

第一，在一定的精度要求下，蒙特卡洛方法的抽样次数与系统的规模无关，因此特别适用于大型电力系统的评估计算。

第二，采用蒙特卡洛方法评估可靠性，不但能够获得概率性指标，而且能够得到频率和持续时间指标，得到的可靠性信息更加丰富、实用。

第三，基于蒙特卡洛方法的程序数学模型相对简单，且容易模拟负荷变化等随机因素和系统的校正控制措施，因此计算结果更加符合工程实际。

电力系统规模日趋扩大、元件众多、控制策略复杂，因此蒙特卡洛法在其可靠性评估中获得了日益广泛的应用。

3.3蒙塔卡洛法的基本内容3.3.1基本参数介绍电力系统元件众多，在可靠性评估中可根据计算需要对发电机组、输电线路、变压器、电抗器、电容器、保护元件、自动重合闸装置、母线等可修复元件进行状态模拟。

假定某可修复元件的故障率和修复率分别为λ、μ，平均无故障工作时间和平均维修时间分别为MTTF、MTTR，则存在以下重要关系式)1.3(/1λ=MTTF )2.3(/1μ=MTTR可修复强迫失效可以通过“运行-停运-运行”的循环过程来模拟，如图一所示：图3.1 可修复元件运行和停运循环过程平均不可用率，其数学形式可由下列三个定义 之一来表达：)3.3(8760MTTR f MTTR MTTF MTTR U ⨯=+=+=μλλλ为失效率（失效次数/年）；μ为修复率（修复次数/年）；MTTR 为平均修复时间（小时）；MTTF 为失效前平均时间（小时）；f 为平均失效频率（失效次数/年）。

d = MTTF/8760及r = MTTR/8760，则d 和r 是以年为单位计的MTTF 和MTTR 。

λ、μ是蒙特卡洛算法中模拟元件持续时间与状态转移特性的基本参数。

其反映的元件状态转移特性如图3一2所示，其数值可通过对元件长期运行的寿命过程和随机状态信息统计得到。

图3.2 可修复元件状态空间图3.3.2非序贯蒙特卡洛模拟法非序贯蒙特卡洛模拟法常常被称为状态抽样法，它被广泛用在电力系统风险评估中。

这个方法的依据是：一个系统状态是所有元件状态的组合，且每一元件状态可由对元件出现在该状态的概率进行抽样来确定。

每一元件可用一个在[0，1]区间的均匀分布来模拟。

假设每一元件有失效和工作两个状态，且元件失效是相互独立的。

令s i 代表元件i 的状态，Q i 代表其失效概率，则对元件i 产生一个在[0，1]区间均匀分布的随机数R i , 使)4.3(0)(1)(0⎩⎨⎧≤≤>=i i ii i Q R if Q R if s 失效状态工作状态具有N 个元件的系统状态由矢量s 表示：S=(s 1，…，s i ，…s N ) (3.5)一个系统状态在抽样中被选定后，即进行系统分析以判断其是否是失效状态，如果是，则对该状态的风险指标函数进行估计。

当抽样的数量足够大时，系统状态s 的抽样频率可作为其概率的无偏估计，即)6.3()s ()(M m s P =式中：M 是抽样数；m (s )是在抽样中系统状态s 出现的次数。

当每一个系统状态的概率通过抽样估计以后，就可计算系统失效概率、系统失效频率、系统失效平均持续时间、以及系统其它风险指标。

非序贯蒙特卡洛法和状态枚举法之间明显的区别在于：如何选择系统状态和如何计算单个系统状态的概率。

在实际应用中，应注意以下几个方面：(1)必要的一步是产生每一个元件的随机数序列，这些随机数必须满足三个基本条件：均匀性、独立性和足够长的重复周期。

(2)蒙特卡洛法是一个波动收敛过程，因此估计出的风险指标总是有一个相应的置信范围。

不能保证增加少量的样本就一定会减少误差，但置信范围确实会随样本数的增加而变窄。

(3)适当的收敛判据是确保蒙特卡洛模拟法精度的关键之一。

方差系数常被用作为终止抽样的判据。

在电力系统风险评估中，不同的风险指标有不同的收敛速度。

已经发现，期望缺供电量(EENS)指标的方差系数收敛速率最低，因此应作为多个指标研究时的收敛判据。

另一种方法是用预定的最大抽样数作为终止抽样的判据。

当模拟过程结束时，校验方差系数是否足够小，如果否，则需要增加样本数再进行新的抽样。

当用户并不知道需要用多少计算时间才能达到足够小的方差系数时，可用这种替代方法。

(4)非序贯模拟过程仅需要元件的失效概率作为抽样过程的输入数据。

这个特点使我们能够容易地同时模拟可修复和老化失效引起的不可用率。

为元件建立两个独立的随机数，一个是模拟可修复失效引起的不可用率，另一个是模拟老化失效引起的不可用率。

(5)状态抽样的理念不仅适用于元件失效事件，而且也可推广应用到电力系统风险评估中其它参数的状态抽样，例如：负荷水平、水文和气候状态等。

而且，这个方法并不局限于年度为基础的模拟，还可很方便地用于进行任意时间段（周、月、季或年度）的模拟。

(6)较之状态枚举法，状态抽样法更适用于规模较大的系统或具有较高元件失效概率的系统评估。

在这些情况下，为获得相同的精度，状态枚举法需要大得多的CPU 时间。

(7)与状态枚举法相似，非序贯蒙特卡洛模拟法不能计及时间相关事件的时序信息，因而得出的系统失效频率和平均失效持续时间乃是近似估计。

3.3.3 序贯蒙特卡洛模拟法序贯蒙特卡洛法是按照时序，在一个时间跨度上进行的模拟。

其中对建立虚拟系统状态转移循环过程有不同的方法。

最通用的是在这里讨论的所谓状态持续时间抽样法。

状态持续时间抽样法是基于对元件状态持续时间的概率分布进行抽样，它分为以下几步：第1步：指定所有元件的初始状态，通常是假设所有元件开始处于运行状态。

第2步：对每一元件停留在当前状态的持续时间进行抽样。

应当设定状态持续时间的概率分布。

对不同的状态，如运行或修复过程，可以假设有不同的状态持续时间概率分布。

例如，下式给出指数分布的状态持续时间的抽样值：)7.3(ln 1i ii R D λ= 式中，R i 是对应于第i 个元件在[0，1]区间均匀分布的随机数。

如果当前的状态是运行状态，则λi 是第i 个元件的失效率；而如果当前的状态是停运状态，则λi 是第i 个元件的修复率。

服从不同概率分布的随机变量的产生方法是不同的。

相关书籍有详细论述。

第3步：在所研究的时间跨度（大量的抽样年）内重复第2步，并记录所有元件的每一状态持续时间的抽样值,即可获得给定时间跨度内每一元件的时序状态转移过程，如图3.3所示。

图3.3 元件时序状态转移过程第4步：组合所有元件的状态转移过程，以建立系统时序状态转移循环过程，如图3.4所示。

图3.4 系统时序状态转移过程第5步：通过对每一个不同系统状态的系统分析，计算风险指标函数。