CN110516287A - 一种计及老化失效的供电系统停运概率计算方法 - Google Patents

Abstract

本发明公开了一种计及老化失效的供电系统停运概率计算方法，该方法同时考虑电网运行方式和环境温度、负荷的不确定性和随机性，通过电网能量管理系统EMS获取电网运行的数据，在考虑电网运行方式的不确定性和随机性时主要是引入发电机、线路、变压器等设备的不确定性运行状态，在考虑环境温度、负荷的不确定性时主要引入环境温度、负荷的不确定性状态，假设用威布尔分布来模拟电网运行方式中各设备组中元件的老化失效过程，在概率分析的基础上计算电网中元件可修复失效引起和老化失效引起的不可用率，为电网调度运行及可靠性评估提供必要的技术支撑。

Description

一种计及老化失效的供电系统停运概率计算方法

技术领域

本发明涉及电力系统及其自动化技术领域，具体涉及一种计及老化失效的供电系统停运概率计算方法。

背景技术

随着电力市场的发展，电力供应的稳定性成为了人们日益关注的问题。电力系统是由大量的发电机、架空输电线路、变压器及电力用户等元件组成，且这些元件大都是可修复元件，元件停运是系统故障停运失效的根本原因。系统运行方式、网络结构和运行条件发生变化时，引发的线路潮流、母线电压、系统频率等的变化将会直接影响到系统元件的停运概率，与此同时，天气、温度等外界因素也会影响到系统元件的停运概率；另外，负荷和电力系统的运行也都充满了不确定性，因而不可能对其准确预测。

一个电力系统是已运行多年且还将继续运行的系统，系统中的元件可能已经历了多次失效和修复。电力系统风险评估方法中通常采用元件的不可用率来描述系统停运模型，大多数情况下，只计入可修复的强迫停运，有时也对计划停运进行模拟，直到目前为止，电力系统的传统风险评估中仅考虑可修复失效，而忽略了老化失效。

作为一种常见现象的设备老化，已日益受到许多公司的重点关注，因此如果元件已接近寿命末期，风险评估中应包含元件的老化失效模式。老化失效是与历史(即元件服役年龄)有关的条件失效事件，不考虑老化失效将极有可能导致低估电力系统风险。由于一个元件发生老化失效将是永久性的，所以老化失效没有修复时间的概念。此外，老化失效的失效率随时间增长而增大，而环境温度和负荷都会引起元件的老化失效。

针对供电系统停运失效概率，以往通常采用潮流计算的方法来加以判断和确认，且建立的供电系统停运模型大都仅考虑了可修复失效，而忽略了老化失效，所以这种计算方法无法确定电网当中发生老化失效的供电系统停运概率，更不能确定发生老化失效的元件个数及损失。

发明内容

本发明的目的在于克服现有技术的缺点与不足，提供一种计及老化失效的供电系统停运概率计算方法，该方法能够确定系统中元件发生老化失效个数、损失和风险。

本发明的目的通过下述技术方案实现：

一种计及老化失效的供电系统停运概率计算方法，包括下述步骤：

S1，从电网能量管理系统EMS获取发电机运行数据，采用非序贯蒙特卡罗模拟方法确定第i台发电机G_i十年内处于强迫停运状态的次数，确定失效率修复率采用非序贯蒙特卡罗模拟方法确定第i台发电机发生老化失效概率的威布尔分布函数，确定形状参数尺度参数

S2，从电网能量管理系统EMS获取线路运行数据，采用非序贯蒙特卡罗模拟方法确定第i条线路L_i十年内处于强迫停运状态的次数，确定失效率修复率采用非序贯蒙特卡罗模拟方法确定第i条线路发生老化失效概率的威布尔分布函数，确定形状参数尺度参数

S3，从电网能量管理系统EMS获取变压器运行数据，采用非序贯蒙特卡罗模拟方法确定第i台变压器T_i十年内处于强迫停运状态的次数，确定失效率修复率采用非序贯蒙特卡罗模拟方法确定第i台变压器发生老化失效概率的威布尔分布函数，确定形状参数尺度参数

S4，计算并列供电系统发电机强迫停运同时存在n台发电机老化失效引起的联合不可用率，计算公式为：

式中为有n台发电机发生老化失效并使这些发电机退出运行时环境温度T_S和负荷功率S_D大于其它台发电机所能承受的最大环境温度和最大输送功率的概率；

S5，计算并列供电系统线路强迫停运同时存在n条线路老化失效引起的联合不可用率，计算公式为：

式中为有n条线路发生老化失效并使这些线路退出运行时环境温度T_S和负荷功率S_D大于其它条线路所能承受的最大环境温度和最大输送功率的概率；

S6，计算并列供电系统变压器强迫停运同时存在n台变压器老化失效引起的联合不可用率，计算公式为：

式中为有n台变压器发生老化失效并使这些变压器退出运行时环境温度T_S和负荷功率S_D大于其它台变压器所能承受的最大环境温度和最大输送功率的概率；

S7，计算并列供电系统计及老化失效的联合停运概率，计算公式为：

优选地，所述步骤S4中，并列供电系统发电机强迫停运同时存在n台发电机老化失效引起的联合不可用率的计算步骤为：

S4.1，从电网能量管理系统EMS获取发电机运行数据，采用非序贯蒙特卡罗模拟方法确定一台发电机发生老化失效并使该台发电机退出运行时环境温度T_S和负荷功率S_D大于其它台发电机所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中G_OUT1为N_G台发电机中一台发电机发生老化失效而永久退出运行的集合；

计算并列供电系统发电机强迫停运同时存在一台发电机老化失效引起的联合不可用率，其计算公式为：

S4.2，从电网能量管理系统EMS获取发电机运行数据，采用非序贯蒙特卡罗模拟方法确定两台发电机发生老化失效并使这两台发电机退出运行时环境温度T_S和负荷功率S_D大于其它台发电机所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中G_OUT2为N_G台发电机中两台发电机发生老化失效而永久退出运行的集合；

计算并列供电系统发电机强迫停运同时存在两台发电机老化失效引起的联合不可用率，其计算公式为：

S4.3，从电网能量管理系统EMS获取发电机运行数据，采用非序贯蒙特卡罗模拟方法确定三台发电机发生老化失效并使这三台发电机退出运行时环境温度T_S和负荷功率S_D大于其它台发电机所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中G_OUT3为N_G台发电机中三台发电机发生老化失效而永久退出运行的集合；

计算并列供电系统发电机强迫停运同时存在三台发电机老化失效引起的联合不可用率，其计算公式为：

S4.4，依照此规律，直到计算并确定了并列供电系统发电机强迫停运同时存在n台发电机老化失效引起的联合不可用率。

优选地，所述步骤S5中，并列供电系统线路强迫停运同时存在n条线路老化失效引起的联合不可用率的计算步骤为：

S5.1，从电网能量管理系统EMS获取线路运行数据，采用非序贯蒙特卡罗模拟方法确定一条线路发生老化失效并使该条线路退出运行时环境温度T_S和负荷功率S_D大于其它条线路所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中L_OUT1为N_L条线路中一条线路发生老化失效而永久退出运行的集合；

计算并列供电系统线路强迫停运同时存在一条线路老化失效引起的联合不可用率，其计算公式为：

S5.2，从电网能量管理系统EMS获取线路运行数据，采用非序贯蒙特卡罗模拟方法确定两条线路发生老化失效并使这两条线路退出运行时环境温度T_S和负荷功率S_D大于其它条线路所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中L_OUT2为N_L条线路中两条线路发生老化失效而永久退出运行的集合；

计算并列供电系统线路强迫停运同时存在两条线路老化失效引起的联合不可用率，其计算公式为：

S5.3，从电网能量管理系统EMS获取线路运行数据，采用非序贯蒙特卡罗模拟方法确定三条线路发生老化失效并使这三条线路退出运行时环境温度T_S和负荷功率S_D大于其它条线路所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中L_OUT3为N_L条线路中三条线路发生老化失效而永久退出运行的集合；

计算并列供电系统线路强迫停运同时存在三条线路老化失效引起的联合不可用率，其计算公式为：

S5.4，依照此规律，直到计算并确定了并列供电系统线路强迫停运同时存在n条线路老化失效引起的联合不可用率。

优选地，所述步骤S6中，并列供电系统变压器强迫停运同时存在n台变压器老化失效引起的联合不可用率的计算步骤为：

S6.1，从电网能量管理系统EMS获取变压器运行数据，采用非序贯蒙特卡罗模拟方法确定一台变压器发生老化失效并使该台变压器退出运行时环境温度T_S和负荷功率S_D大于其它台变压器所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中T_OUT1为N_T台变压器中一台变压器发生老化失效而永久退出运行的集合；

计算并列供电系统变压器强迫停运同时存在一台变压器老化失效引起的联合不可用率，其计算公式为：

S6.2，从电网能量管理系统EMS获取变压器运行数据，采用非序贯蒙特卡罗模拟方法确定两台变压器发生老化失效并使这两台变压器退出运行时环境温度T_S和负荷功率S_D大于其它台变压器所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中T_OUT2为N_T台变压器中两台变压器发生老化失效而永久退出运行的集合；

计算并列供电系统变压器强迫停运同时存在两台变压器老化失效引起的联合不可用率，其计算公式为：

S6.3，从电网能量管理系统EMS获取变压器运行数据，采用非序贯蒙特卡罗模拟方法确定三台变压器发生老化失效并使这三台变压器退出运行时环境温度T_S和负荷功率S_D大于其它台变压器所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中T_OUT3为N_T台变压器中三台变压器发生老化失效而永久退出运行的集合；

计算并列供电系统变压器强迫停运同时存在三台变压器老化失效引起的联合不可用率，其计算公式为：

S6.4，依照此规律，直到计算并确定了并列供电系统变压器强迫停运同时存在n台变压器老化失效引起的联合不可用率。

本发明与现有技术相比具有以下的有益效果：

本发明可计算出在一定运行周期(1年、5年或10年等)内，多发电机、线路、变压器并列供电系统由强迫停运(可修复失效)以及在给定各元件已服役T年的条件下同时给定后续时间期间t内由老化失效引起的联合不可用率，为系统风险评估提供技术方法，为电网可靠运行提供技术支撑。

附图说明

图1为本发明所针对的老化失效模型所给定的后续时间期间t内的时间分段示意图；

图2为本发明所针对的并列供电系统组成和连接关系示意图；

图3为本发明的流程框图。

图中附图标记为：1、并列运行的第一台发电机；2、并列运行的第N_G台发电机；3、发电机机端母线；4、并列运行的第一条线路；5、并列运行的第N_L条线路；6、变压器高压母线；7、并列运行的第一台变压器；8、并列运行的第N_T台变压器；9、变压器低压母线；10、负荷。

具体实施方式

下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。

本发明提供了一种计及老化失效的供电系统停运概率计算方法，该方法的基本原理是同时考虑电网运行方式和环境温度、负荷的不确定性和随机性，通过电网能量管理系统EMS获取电网运行的数据，在考虑电网运行方式的不确定性和随机性时主要是引入发电机、线路、变压器等设备的不确定性运行状态，在考虑环境温度、负荷的不确定性时主要引入环境温度、负荷的不确定性状态，假设用威布尔分布来模拟电网运行方式中各设备组中元件的老化失效过程，在概率分析的基础上计算电网中元件可修复失效引起和老化失效引起的不可用率，为电网调度运行及可靠性评估提供必要的技术支撑。

该方法针对图2所示的多发电机、多线路、多变压器并列供电系统，并考虑发电机、输电线路和变压器中发生老化失效的元件个数，再结合并列发电机运行方式改变、输电线路运行方式改变和变压器运行方式改变的不确定性和随机性，提出并列供电系统中计及老化失效的停运概率计算方法；图2中，由N_G台发电机G₁、G₂、…、G_NG并列运行，N_L条线路L₁、L₂、L₃、…、L_NL并列运行和N_T台变压器T₁、T₂、T₃、…、T_NT并列运行组成供电系统，假设负荷功率为S_D(S_D＝P_D+jQ_D)。这是一种由多发电机、多线路、多变压器及负荷组成的并列供电系统。元件老化失效引起的不可用率被定义为，给定元件已服役T年的条件下在时间t内不可用的概率。为分析老化失效引起的不可用率，把考察的后续时间(t)按时间长度Δx等分成N个时段，分段的端点用t_j(j＝0，1，2，...，N)标示，如图1所示。

如图3所示，一种计及老化失效的供电系统停运概率计算方法，包括下述步骤：

图3中的步骤1描述发电机强迫停运概率的参数确定及其发生老化失效概率的威布尔分布函数确定和参数计算的过程和方法。从电网能量管理系统EMS获取并列运行发电机的运行数据，按照抽取10年(1年作为每一个时段)的数据规模进行处理、计算和分析。对于可修复失效，重点抽取发电机强迫停运次数、失效频率、修复时间及其发生的条件(温度、负载率)等数据，然后通过计算得到发电机的失效率和修复率；对于老化失效，重点抽取各台发电机投运和退役年份及其发生的条件(温度、负载率)等数据，进而确定各发电机的服役年龄及存活概率，采用概率分析方法验证这些数据是否具备威布尔分布特征，并确定其概率分布函数。具体来说，从电网能量管理系统EMS获取发电机运行数据，采用非序贯蒙特卡罗模拟方法确定第i台发电机G_i十年内处于强迫停运(可修复失效)状态的次数，确定失效率修复率采用非序贯蒙特卡罗模拟方法确定第i台发电机发生老化失效概率的威布尔分布函数，确定形状参数尺度参数

图3中的步骤2描述线路强迫停运概率的参数确定及其发生老化失效概率的威布尔分布函数确定和参数计算的过程和方法。从电网能量管理系统EMS获取并列运行线路的运行数据，按照抽取10年(1年作为每一个时段)的数据规模进行处理、计算和分析。对于可修复失效，重点抽取线路强迫停运次数、失效频率、修复时间及其发生的条件(温度、负载率)等数据，然后通过计算得到线路的失效率和修复率；对于老化失效，重点抽取各条线路投运和退役年份及其发生的条件(温度、负载率)等数据，进而确定各线路的服役年龄及存活概率，采用概率分析方法验证这些数据是否具备威布尔分布特征，并确定其概率分布函数。具体来说，采用非序贯蒙特卡罗模拟方法确定第i条线路L_i十年内处于强迫停运(可修复失效)状态的次数，确定失效率修复率采用非序贯蒙特卡罗模拟方法确定第i条线路发生老化失效概率的威布尔分布函数，确定形状参数尺度参数

图3中的步骤3描述变压器强迫停运概率的参数确定及其发生老化失效概率的威布尔分布函数确定和参数计算的过程和方法。从电网能量管理系统EMS获取并列运行变压器的运行数据，按照抽取10年(1年作为每一个时段)的数据规模进行处理、计算和分析。对于可修复失效，重点抽取变压器强迫停运次数、失效频率、修复时间及其发生的条件(温度、负载率)等数据，然后通过计算得到变压器的失效率和修复率；对于老化失效，重点抽取各台变压器投运和退役年份及其发生的条件(温度、负载率)等数据，进而确定各变压器的服役年龄及存活概率，采用概率分析方法验证这些数据是否具备威布尔分布特征，并确定其概率分布函数。具体来说，采用非序贯蒙特卡罗模拟方法确定第i台变压器T_i十年内处于强迫停运(可修复失效)状态的次数，确定失效率修复率采用非序贯蒙特卡罗模拟方法确定第i台变压器发生老化失效概率的威布尔分布函数，确定形状参数尺度参数

图3中的步骤4描述计算并列供电系统发电机强迫停运同时存在n台发电机老化失效引起的联合不可用率，计算公式为：

式中为有n台发电机发生老化失效并使这些发电机退出运行时环境温度T_S和负荷功率S_D大于其它台发电机所能承受的最大环境温度和最大输送功率的概率；

所述步骤S4中，并列供电系统发电机强迫停运同时存在n台发电机老化失效引起的联合不可用率的计算步骤为：

S4.1，从电网能量管理系统EMS获取发电机运行数据，采用非序贯蒙特卡罗模拟方法确定一台发电机发生老化失效并使该台发电机退出运行时环境温度T_S和负荷功率S_D大于其它台发电机所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中G_OUT1为N_G台发电机中一台发电机发生老化失效而永久退出运行的集合；

计算并列供电系统发电机强迫停运同时存在一台发电机老化失效引起的联合不可用率，其计算公式为：

S4.2，从电网能量管理系统EMS获取发电机运行数据，采用非序贯蒙特卡罗模拟方法确定两台发电机发生老化失效并使这两台发电机退出运行时环境温度T_S和负荷功率S_D大于其它台发电机所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中G_OUT2为N_G台发电机中两台发电机发生老化失效而永久退出运行的集合；

计算并列供电系统发电机强迫停运同时存在两台发电机老化失效引起的联合不可用率，其计算公式为：

S4.3，从电网能量管理系统EMS获取发电机运行数据，采用非序贯蒙特卡罗模拟方法确定三台发电机发生老化失效并使这三台发电机退出运行时环境温度T_S和负荷功率S_D大于其它台发电机所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中G_OUT3为N_G台发电机中三台发电机发生老化失效而永久退出运行的集合；

计算并列供电系统发电机强迫停运同时存在三台发电机老化失效引起的联合不可用率，其计算公式为：

S4.4，依照此规律，直到计算并确定了并列供电系统发电机强迫停运同时存在n台发电机老化失效引起的联合不可用率。

图3中的步骤5描述计算并列供电系统线路强迫停运同时存在n条线路老化失效引起的联合不可用率，计算公式为：

式中为有n条线路发生老化失效并使这些线路退出运行时环境温度T_S和负荷功率S_D大于其它条线路所能承受的最大环境温度和最大输送功率的概率；

所述步骤S5中，并列供电系统线路强迫停运同时存在n条线路老化失效引起的联合不可用率的计算步骤为：

S5.1，从电网能量管理系统EMS获取线路运行数据，采用非序贯蒙特卡罗模拟方法确定一条线路发生老化失效并使该条线路退出运行时环境温度T_S和负荷功率S_D大于其它条线路所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中L_OUT1为N_L条线路中一条线路发生老化失效而永久退出运行的集合；

计算并列供电系统线路强迫停运同时存在一条线路老化失效引起的联合不可用率，其计算公式为：

S5.2，从电网能量管理系统EMS获取线路运行数据，采用非序贯蒙特卡罗模拟方法确定两条线路发生老化失效并使这两条线路退出运行时环境温度T_S和负荷功率S_D大于其它条线路所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中L_OUT2为N_L条线路中两条线路发生老化失效而永久退出运行的集合；

计算并列供电系统线路强迫停运同时存在两条线路老化失效引起的联合不可用率，其计算公式为：

S5.3，从电网能量管理系统EMS获取线路运行数据，采用非序贯蒙特卡罗模拟方法确定三条线路发生老化失效并使这三条线路退出运行时环境温度T_S和负荷功率S_D大于其它条线路所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中L_OUT3为N_L条线路中三条线路发生老化失效而永久退出运行的集合；

计算并列供电系统线路强迫停运同时存在三条线路老化失效引起的联合不可用率，其计算公式为：

S5.4，依照此规律，直到计算并确定了并列供电系统线路强迫停运同时存在n条线路老化失效引起的联合不可用率。

图3中的步骤6描述计算并列供电系统变压器强迫停运同时存在n台变压器老化失效引起的联合不可用率，计算公式为：

式中为有n台变压器发生老化失效并使这些变压器退出运行时环境温度T_S和负荷功率S_D大于其它台变压器所能承受的最大环境温度和最大输送功率的概率；

所述步骤S6中，并列供电系统变压器强迫停运同时存在n台变压器老化失效引起的联合不可用率的计算步骤为：

S6.1，从电网能量管理系统EMS获取变压器运行数据，采用非序贯蒙特卡罗模拟方法确定一台变压器发生老化失效并使该台变压器退出运行时环境温度T_S和负荷功率S_D大于其它台变压器所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中T_OUT1为N_T台变压器中一台变压器发生老化失效而永久退出运行的集合；

计算并列供电系统变压器强迫停运同时存在一台变压器老化失效引起的联合不可用率，其计算公式为：

S6.2，从电网能量管理系统EMS获取变压器运行数据，采用非序贯蒙特卡罗模拟方法确定两台变压器发生老化失效并使这两台变压器退出运行时环境温度T_S和负荷功率S_D大于其它台变压器所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中T_OUT2为N_T台变压器中两台变压器发生老化失效而永久退出运行的集合；

计算并列供电系统变压器强迫停运同时存在两台变压器老化失效引起的联合不可用率，其计算公式为：

S6.3，从电网能量管理系统EMS获取变压器运行数据，采用非序贯蒙特卡罗模拟方法确定三台变压器发生老化失效并使这三台变压器退出运行时环境温度T_S和负荷功率S_D大于其它台变压器所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中T_OUT3为N_T台变压器中三台变压器发生老化失效而永久退出运行的集合；

计算并列供电系统变压器强迫停运同时存在三台变压器老化失效引起的联合不可用率，其计算公式为：

S6.4，依照此规律，直到计算并确定了并列供电系统变压器强迫停运同时存在n台变压器老化失效引起的联合不可用率。

图3中的步骤7描述计算并列供电系统计及老化失效的联合停运概率，计算公式为：

计及老化失效的并列供电系统停运概率，与一定运行周期内发电机、线路和变压器发生可修复失效引起的不可用率有关系；与一定运行周期内发电机、线路和变压器发生老化失效引起的不可用率有关系；与一定运行周期内发电机、线路和变压器所处环境温度、负荷有关系。

本发明可计算出在一定运行周期(1年、5年或10年等)内，多发电机、线路、变压器并列供电系统由强迫停运(可修复失效)以及在给定各元件已服役T年的条件下同时给定后续时间期间t内由老化失效引起的联合不可用率，为系统风险评估提供技术方法，为电网可靠运行提供技术支撑。

上述为本发明较佳的实施方式，但本发明的实施方式并不受上述内容的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

Claims (4)

1.一种计及老化失效的供电系统停运概率计算方法，其特征在于，包括下述步骤：

S1，从电网能量管理系统EMS获取发电机运行数据，采用非序贯蒙特卡罗模拟方法确定第i台发电机G_i十年内处于强迫停运状态的次数，确定失效率修复率采用非序贯蒙特卡罗模拟方法确定第i台发电机发生老化失效概率的威布尔分布函数，确定形状参数尺度参数

S2，从电网能量管理系统EMS获取线路运行数据，采用非序贯蒙特卡罗模拟方法确定第i条线路L_i十年内处于强迫停运状态的次数，确定失效率修复率采用非序贯蒙特卡罗模拟方法确定第i条线路发生老化失效概率的威布尔分布函数，确定形状参数尺度参数

S3，从电网能量管理系统EMS获取变压器运行数据，采用非序贯蒙特卡罗模拟方法确定第i台变压器T_i十年内处于强迫停运状态的次数，确定失效率修复率采用非序贯蒙特卡罗模拟方法确定第i台变压器发生老化失效概率的威布尔分布函数，确定形状参数尺度参数

S4，计算并列供电系统发电机强迫停运同时存在n台发电机老化失效引起的联合不可用率，计算公式为：

式中为有n台发电机发生老化失效并使这些发电机退出运行时环境温度T_S和负荷功率S_D大于其它台发电机所能承受的最大环境温度和最大输送功率的概率；

S5，计算并列供电系统线路强迫停运同时存在n条线路老化失效引起的联合不可用率，计算公式为：

式中为有n条线路发生老化失效并使这些线路退出运行时环境温度T_S和负荷功率S_D大于其它条线路所能承受的最大环境温度和最大输送功率的概率；

S6，计算并列供电系统变压器强迫停运同时存在n台变压器老化失效引起的联合不可用率，计算公式为：

式中为有n台变压器发生老化失效并使这些变压器退出运行时环境温度T_S和负荷功率S_D大于其它台变压器所能承受的最大环境温度和最大输送功率的概率；

S7，计算并列供电系统计及老化失效的联合停运概率，计算公式为：

2.根据权利要求1所述的计及老化失效的供电系统停运概率计算方法，其特征在于，所述步骤S4中，并列供电系统发电机强迫停运同时存在n台发电机老化失效引起的联合不可用率的计算步骤为：

S4.1，从电网能量管理系统EMS获取发电机运行数据，采用非序贯蒙特卡罗模拟方法确定一台发电机发生老化失效并使该台发电机退出运行时环境温度T_S和负荷功率S_D大于其它台发电机所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中G_OUT1为N_G台发电机中一台发电机发生老化失效而永久退出运行的集合；

计算并列供电系统发电机强迫停运同时存在一台发电机老化失效引起的联合不可用率，其计算公式为：

S4.2，从电网能量管理系统EMS获取发电机运行数据，采用非序贯蒙特卡罗模拟方法确定两台发电机发生老化失效并使这两台发电机退出运行时环境温度T_S和负荷功率S_D大于其它台发电机所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中G_OUT2为N_G台发电机中两台发电机发生老化失效而永久退出运行的集合；

计算并列供电系统发电机强迫停运同时存在两台发电机老化失效引起的联合不可用率，其计算公式为：

S4.3，从电网能量管理系统EMS获取发电机运行数据，采用非序贯蒙特卡罗模拟方法确定三台发电机发生老化失效并使这三台发电机退出运行时环境温度T_S和负荷功率S_D大于其它台发电机所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中G_OUT3为N_G台发电机中三台发电机发生老化失效而永久退出运行的集合；

计算并列供电系统发电机强迫停运同时存在三台发电机老化失效引起的联合不可用率，其计算公式为：

S4.4，依照此规律，直到计算并确定了并列供电系统发电机强迫停运同时存在n台发电机老化失效引起的联合不可用率。

3.根据权利要求1所述的计及老化失效的供电系统停运概率计算方法，其特征在于，所述步骤S5中，并列供电系统线路强迫停运同时存在n条线路老化失效引起的联合不可用率的计算步骤为：

S5.1，从电网能量管理系统EMS获取线路运行数据，采用非序贯蒙特卡罗模拟方法确定一条线路发生老化失效并使该条线路退出运行时环境温度T_S和负荷功率S_D大于其它条线路所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中L_OUT1为N_L条线路中一条线路发生老化失效而永久退出运行的集合；

计算并列供电系统线路强迫停运同时存在一条线路老化失效引起的联合不可用率，其计算公式为：

S5.2，从电网能量管理系统EMS获取线路运行数据，采用非序贯蒙特卡罗模拟方法确定两条线路发生老化失效并使这两条线路退出运行时环境温度T_S和负荷功率S_D大于其它条线路所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中L_OUT2为N_L条线路中两条线路发生老化失效而永久退出运行的集合；

计算并列供电系统线路强迫停运同时存在两条线路老化失效引起的联合不可用率，其计算公式为：

S5.3，从电网能量管理系统EMS获取线路运行数据，采用非序贯蒙特卡罗模拟方法确定三条线路发生老化失效并使这三条线路退出运行时环境温度T_S和负荷功率S_D大于其它条线路所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中L_OUT3为N_L条线路中三条线路发生老化失效而永久退出运行的集合；

计算并列供电系统线路强迫停运同时存在三条线路老化失效引起的联合不可用率，其计算公式为：

S5.4，依照此规律，直到计算并确定了并列供电系统线路强迫停运同时存在n条线路老化失效引起的联合不可用率。

4.根据权利要求1所述的计及老化失效的供电系统停运概率计算方法，其特征在于，所述步骤S6中，并列供电系统变压器强迫停运同时存在n台变压器老化失效引起的联合不可用率的计算步骤为：

S6.1，从电网能量管理系统EMS获取变压器运行数据，采用非序贯蒙特卡罗模拟方法确定一台变压器发生老化失效并使该台变压器退出运行时环境温度T_S和负荷功率S_D大于其它台变压器所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中T_OUT1为N_T台变压器中一台变压器发生老化失效而永久退出运行的集合；

计算并列供电系统变压器强迫停运同时存在一台变压器老化失效引起的联合不可用率，其计算公式为：

S6.2，从电网能量管理系统EMS获取变压器运行数据，采用非序贯蒙特卡罗模拟方法确定两台变压器发生老化失效并使这两台变压器退出运行时环境温度T_S和负荷功率S_D大于其它台变压器所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中T_OUT2为N_T台变压器中两台变压器发生老化失效而永久退出运行的集合；

计算并列供电系统变压器强迫停运同时存在两台变压器老化失效引起的联合不可用率，其计算公式为：

S6.3，从电网能量管理系统EMS获取变压器运行数据，采用非序贯蒙特卡罗模拟方法确定三台变压器发生老化失效并使这三台变压器退出运行时环境温度T_S和负荷功率S_D大于其它台变压器所能承受的最大环境温度和最大输送功率的概率其计算公式为：

式中T_OUT3为N_T台变压器中三台变压器发生老化失效而永久退出运行的集合；

计算并列供电系统变压器强迫停运同时存在三台变压器老化失效引起的联合不可用率，其计算公式为：

S6.4，依照此规律，直到计算并确定了并列供电系统变压器强迫停运同时存在n台变压器老化失效引起的联合不可用率。