本发明提供一种基于静态可靠性概率指标的电网可靠性水平判定方法，本发明考虑了元件可能停运的概率、运行环境中的不确定性和多重约束导致的风险增加等因素，能计算元件、负荷点和系统的可靠性指标，记录具体的切负荷节点和切负荷量，并给出了判定电网可靠性水平的一种评判标准，能对电网的静态可靠性进行较为全面和客观的评价，对电网的可靠性水平进行判定。因此，基于静态可靠性概率指标的电网可靠性水平判定方法具有较好的实际指导意义和应用价值。
1.一种基于静态可靠性概率指标的电网可靠性水平判定方法，其特征在于，所述方法包括如下步骤：(1)输入元件的可靠性参数；
(2)搭建电力系统，形成初始状态的潮流计算模型，进行潮流计算；
(4)判断所述电力系统是否解列，若是则根据解列信息形成子网络集，否则在判定电力系统电力供给充足后转入步骤(6)；
(5)选择一个子网络，判断该子网络是否采取了有功校正措施；
(6)在初始状态的潮流计算模型中删除故障元件，形成故障后的潮流计算数据文件，再次进行潮流计算；
(7)根据步骤(6)计算结果进行节点电压越限和线路过载检测判断；
(8)输出节点电压越限信息、线路过载信息和节点负荷削减量信息；
(9)对所有的子网络和所有的状态进行上述步骤分析，并计算系统、负荷点和元件的可靠性指标；
(10)结合电网规模和可靠性历史数据，确定评价电网静态可靠性水平的基准值，判定电网可靠性水平；
步骤(1)中，所述可靠性参数包括元件故障率、故障平均修复时间、计划检修率和计划检修时间；
步骤(9)中，所述计算系统和负荷点可靠性指标的方法，包括如下步骤：步骤9-1、计算系统和负荷节点的可靠性水平指标，公式如下：式中，PSLOLP为停电概率，当不考虑动态切负荷时，认为该指标即为系统停电概率；P(Xi)为系统状态的概率分布函数；P(Xi,k)为元件的概率分布函数；
式中，FSLOLF为停电频率，单位：次/年；当不考虑动态切负荷时，认为该指标即为系统停电频率；SLC为存在切负荷的系统状态集合；Mdown为系统中停运元件集合；Mup为系统中运行元件集合；μk为修复率；λk为故障率；
式中，PSEPNS为停电功率期望值，单位：MW；当不考虑动态切负荷时，认为该指标即为系统停电功率期望值；PS,i为第i个系统状态的停电电量，单位为：MWh；
式中，ESEENS为停电电量期望值，单位：MWh/年；当不考虑动态切负荷时，认为该指标即为系统停电电量期望值；T为研究时段小时数，一年取8760小时；
式中，HSLOLE为年停电小时数，单位为：h/年；
式中，DSLOLD均次停电时间，单位：h/次；当不考虑动态切负荷时，认为该指标即为系统停电时间期望值；
式中，ESAEC为均次停电电量，单位：MWh/次；
PSMPC＝max(PS,i)，i＝1,2,...,NLC式中，PSMPC为最大停电功率，单位：MW；NLC为存在切负荷的系统状态的总个数；
ESMEC＝max(ES,i)，i＝1,2,...,NLC式中，ESMEC为最大停电电量，单位：MWh；ES,i为第i个系统状态的停电电量，单位：MWh；
式中，ISISI为停电严重性指标，单位：分钟；Lmax为系统年峰荷；
步骤9-2、计算负荷节点的可靠性水平的指标，公式如下：式中，PVIO,l为负荷点l存在电压越界的概率，SVIO,l为负荷点l存在电压越界的系统状态集合；
步骤9-3、计算系统可靠性指标在负荷节点间分配情况的平均值指标，公式如下：式中，FSYAIF为负荷点平均停电频率，单位：次/(年.负荷点)；NL为系统中负荷点个数；
式中，HSYAIH为负荷点年平均停电时数，单位：h/(年.负荷点)；
式中，DSYAID为负荷点平均停电持续时间，单位：h/(次.负荷点)；
式中，PSYAPC为负荷点平均停电功率，单位：MW/(年.负荷点)；
式中，ESYAEC为负荷点平均停电电量，单位：MWh/(年.负荷点)；
步骤(9)中，计算元件可靠性指标，公式如下：式中，POL,k为线路k的过负荷概率；SOL,k为线路k存在过载的系统状态集合；
ΔPmax,k＝max(ΔPk,j)，(j＝1,2,...,NOL,k)式中，ΔPk,j为线路k过负荷量，NOL,k为线路k存在过载的系统状态的总个数，ΔPmax,k为线路k的有功最大过负荷量；
步骤(10)中，所述判定电网可靠性水平包括：将可靠性最低的元件不可用率作为阈值，将步骤(9)中所计算的可靠性指标与所述阈值做比较，若低于阈值，则认为可靠性满足要求，否则不满足要求。
2.根据权利要求1所述判定方法，其特征在于，步骤(2)中，所述电力系统包括发电机、负荷、平衡节点、变压器和线路。
3.根据权利要求1所述判定方法，其特征在于，步骤(4)中，所述解列信息为所述电力系统的一部分和所述电力系统其它部分的联系切断，分成相互独立、互不联系的部分。
4.根据权利要求1所述判定方法，其特征在于，步骤(5)中，若采取了有功校正措施，则转入步骤(8)输出节点负荷削减量信息；若未采取有功校正措施，则继续判断该子网络是否包含了负荷、发电机节点和平衡节点信息；如果子网络包含了负荷、发电机节点和平衡节点，则转入步骤(6)；否则，将采取校正措施直至子网络包含完整的负荷、发电机节点和平衡节点信息。
5.根据权利要求1所述判定方法，其特征在于，步骤(7)中，所述进行节点电压越限和线路过载检测判断包括，若存在节点电压越限或线路过载问题，则采取校正措施，并跳转至步骤(8)。
一种基于静态可靠性概率指标的电网可靠性水平判定方法
技术领域
[0001] 本发明涉及一种可靠性水平判定方法，具体涉及一种基于静态可靠性概率指标的电网可靠性水平判定方法。
背景技术
[0002] 随着电网元件的数量、种类不断增多，电力系统自动化程度不断提高，以及用户对电能质量要求的提高，使得电力系统的规划、运行、维修等工作面临极大的挑战。近年来，国内外频发的大面积停电事故，带来了巨大的经济损失和沉重的社会影响，因此保持电力系统的安全可靠运行，不仅仅可以为社会生产带来效益保障，更可以保证全社会健康稳定的发展。传统的可靠性评估往往仅考虑单一的或少量的几种不确定因素，仅计算常规的电力或电量损失指标，且未能给出判定电网可靠性水平的评判标准。
[0003] 一种基于静态可靠性概率指标的电网可靠性水平判定方法，考虑了元件可能停运的概率、运行环境中的不确定性和多重约束导致的风险增加等因素，能计算元件、负荷点和系统的可靠性指标，记录具体的切负荷节点和切负荷量，并给出了判定电网可靠性水平的一种评判标准，能对电网的静态可靠性进行较为全面和客观的评价，对电网的可靠性水平进行判定。因此，基于静态可靠性概率指标的电网可靠性水平判定方法具有较好的实际指导意义和应用价值。
[0004] 在静态可靠性概率指标计算中，需要采用概率可靠性评估模型。
发明内容
[0005] 为了克服上述现有技术的不足，本发明提供一种基于静态可靠性概率指标的电网可靠性水平判定方法。本发明基于静态可靠性概率指标的电网可靠性水平判定方法能对电网的静态可靠性进行较为全面和客观的评价，对电网的可靠性水平进行判定。
[0006] 为了实现上述发明目的，本发明采取如下技术方案：
[0007] 一种基于静态可靠性概率指标的电网可靠性水平判定方法，所述方法包括如下步骤：
[0008] (1)输入元件的可靠性参数；
[0009] (2)搭建电力系统，形成初始状态的潮流计算模型，进行潮流计算；
[0010] (3)选择一个故障事件，计算其发生的概率；
[0011] (4)判断所述电力系统是否解列，若是则根据解列信息形成子网络集，否则在判定电力系统电力供给充足后转入步骤(6)；
[0012] (5)选择一个子网络，判断该子网络是否采取了有功校正措施；
[0013] (6)在初始状态的潮流计算模型中删除故障元件，形成故障后的潮流计算数据文件，再次进行潮流计算；
[0014] (7)根据步骤(6)计算结果进行节点电压越限和线路过载检测判断；
[0015] (8)输出节点电压越限信息、线路过载信息和节点负荷削减量信息；
[0016] (9)对所有的子网络和所有的状态进行上述步骤分析，并计算系统、负荷点和元件的可靠性指标；
[0017] (10)结合电网规模和可靠性历史数据，确定评价电网静态可靠性水平的基准值，判定电网可靠性水平。
[0018] 优选的，所述步骤(1)中，所述可靠性参数包括元件故障率、故障平均修复时间、计划检修率和计划检修时间。
[0019] 优选的，所述步骤(2)中，所述电力系统数据包括发电机、负荷、平衡节点、变压器和线路。
[0020] 优选的，所述步骤(4)中，所述解列信息为所述电力系统的一部分和所述电力系统其它部分的联系切断，分成相互独立、互不联系的部分。
[0021] 优选的，所述步骤(5)中，若采取了有功校正措施，则转入步骤(8)输出节点负荷削减量信息；若未采取有功校正措施，则继续判断该子网络是否包含了负荷、发电机节点和平衡节点信息；如果子网络包含了负荷、发电机节点和平衡节点，则转入步骤(6)；否则，将采取校正措施直至子网络包含完整的负荷、发电机节点和平衡节点信息。
[0022] 优选的，所述步骤(7)中，所述进行节点电压越限和线路过载检测判断包括，若存在节点电压越限或线路过载问题，则采取校正措施，并跳转至步骤(8)。
[0023] 优选的，所述步骤(9)中，所述计算系统和负荷点可靠性指标的方法，包括如下步骤：
[0024] 步骤9-1、计算系统和负荷节点的可靠性水平指标，公式如下：
[0025]
[0026] 式中，PSLOLP为停电概率，当不考虑动态切负荷时，认为该指标即为系统停电概率；P(Xi)为系统状态的概率分布函数；P(Xi,k)为元件的概率分布函数；
[0027]
[0028] 式中，FSLOLF为停电频率，单位：次/年；当不考虑动态切负荷时，认为该指标即为系统停电频率；SLC为存在切负荷的系统状态集合；Mdown为系统中停运元件集合；Mup为系统中运行元件集合；μk为修复率；λk为故障率；
[0029]
[0030] 式中，PSEPNS为停电功率期望值，单位：MW；当不考虑动态切负荷时，认为该指标即为系统停电功率期望值；PS,i为第i个系统状态的停电电量，单位为：MWh；
[0031] ESEENS＝PSEPNS·T
[0032] 式中，ESEENS为停电电量期望值，单位：MWh/年；当不考虑动态切负荷时，认为该指标即为系统停电电量期望值；T为研究时段小时数，一年取8760小时；
[0033] HSLOLE＝PSLOLP·T
[0034] 式中，HSLOLE为年停电小时数，单位为：h/年；
[0035] DSLOLD＝HSLOLE/FSLOLF
[0036] 式中，DSLOLD均次停电时间，单位：h/次；当不考虑动态切负荷时，认为该指标即为系统停电时间期望值；
[0037] PSAPC＝PSEPNS/FSLOLF
[0038] 式中，PSAPC为均次停电功率，单位：MW/次；
[0039] ESAEC＝ESEENS/FSLOLF
[0040] 式中，ESAEC为均次停电电量，单位：MWh/次；
[0041] PSMPC＝max(PS,i)，i＝1,2,...,NLC
[0042] 式中，PSMPC为最大停电功率，单位：MW；NLC为存在切负荷的系统状态的总个数；
[0043] ESMEC＝max(ES,i)，i＝1,2,...,NLC
[0044] 式中，ESMEC为最大停电电量，单位：MWh；ES,i为第i个系统状态的停电电量，单位：MWh；
[0045] ISISI＝ESEENS/Lmax×60
[0046] 式中，ISISI为停电严重性指标，单位：分钟；Lmax为系统年峰荷；
[0047] 步骤9-2、计算负荷节点的可靠性水平的指标，公式如下：
[0048]
[0049] 式中，PVIO,l为负荷点l存在电压越界的概率，SVIO,l为负荷点l存在电压越界的系统状态集合；
[0050] 步骤9-3、计算系统可靠性指标在负荷节点间分配情况的平均值指标，公式如下：
[0051]
[0052] 式中，FSYAIF为负荷点平均停电频率，单位：次/(年.负荷点)；NL为系统中负荷点个数；
[0053]
[0054] 式中，HSYAIH为负荷点年平均停电时数，单位：h/(年.负荷点)；
[0055]
[0056] 式中，DSYAID为负荷点平均停电持续时间，单位：h/(次.负荷点)；
[0057] PSYAPC＝PSEPNS/NL
[0058] 式中，PSYAPC为负荷点平均停电功率，单位：MW/(年.负荷点)；
[0059] ESYAEC＝ESEENS/NL
[0060] 式中，ESYAEC为负荷点平均停电电量，单位：MWh/(年.负荷点)；
[0061] 优选的，所述步骤(9)中，计算元件可靠性指标，公式如下：
[0062]
[0063] 式中，POL,k为线路k的过负荷概率；SOL,k为线路k存在过载的系统状态集合；
[0064] ΔPmax,k＝max(ΔPk,j)，(j＝1,2,...,NOL,k)
[0065] 式中，ΔPk,j为线路k过负荷量，NOL,k为线路k存在过载的系统状态的总个数，ΔPmax,k为线路k的有功最大过负荷量。
[0066] 优选的，所述步骤(10)中，所述判定电网可靠性水平包括：将可靠性最低的元件不可用率作为阈值，将步骤(9)中所计算的可靠性指标与所述阈值做比较，若低于阈值，则认为可靠性满足要求，否则不满足要求。
[0067] 与现有技术相比，本发明的有益效果在于：
[0068] 本发明考虑了元件可能停运的概率、运行环境中的不确定性和多重约束导致的风险增加等因素，能计算元件、负荷点和系统的可靠性指标，能对电网的静态可靠性进行较为全面和客观的评价，记录具体的切负荷节点和切负荷量、描述系统在不同区域的失负荷情况，解决了电网静态可靠性评估缺乏统一标准的问题。
附图说明
[0069] 图1是本发明提供的一种基于静态可靠性概率指标的电网可靠性水平判定方法[0070] 图2是元件二状态模型
[0071] 图3是某地区电网系统局部接线图
具体实施方式
[0072] 下面结合附图对本发明作进一步详细说明。
[0073] 如图1所示，一种基于静态可靠性概率指标的电网可靠性水平判定方法，该方法包括如下步骤：
[0074] (1)根据计算需求搭建数据；输入元件可靠性参数，包括：元件的故障率(次/年)、故障平均修复时间(小时/次)、计划检修率(次/年)、计划检修时间(小时/次)等。
[0075] (2)搭建发电机、负荷、平衡节点、变压器和线路数据，形成初始状态的潮流计算模型，进行潮流计算。
[0076] (3)选择一个故障事件，计算其发生的概率。
[0077] (4)判断所述电力系统是否解列，若是则根据解列信息形成子网络集，否则在判定电力系统电力供给充足后转入步骤(6)。
[0078] (5)选择一个子网络，判断该子网络是否采取了有功校正措施。若采取了有功校正措施，则转入(8)输出节点负荷削减量信息；若未有功校正措施，则继续判断该子网络是否包含了负荷、发电机节点和平衡节点信息。如果子网络包含了负荷、发电机节点和平衡节点，则转入(6)；否则，将采取校正措施直至子网络包含完整的负荷、发电机节点和平衡节点信息。
[0079] (6)形成BPA潮流数据文件，进行故障状态下的潮流计算。
[0080] (7)根据(6)的计算结果进行节点电压越限和线路过载检查。若存在节点电压越限和(或)线路过载问题，则采取校正措施，并转入(8)输出节点电压越限、线路过载信息以及相应的负荷削减量信息。
[0081] (8)输出节点电压越限信息、线路过载信息、节点负荷削减量信息。
[0082] (9)对所有的子网络和所有的状态进行上述分析，并计算系统、负荷点以及元件的可靠性指标。
[0083] (10)结合电网规模、可靠性历史数据，确定评价电网静态可靠性水平的基准值，判[0084] 定电网是否满足安全性要求。
[0085] 所述步骤(9)中，所述计算系统和负荷点可靠性指标的方法，包括如下步骤：
[0086] 步骤4-1、计算系统和负荷节点的可靠性水平指标，公式如下：
[0087]
[0088] 式中，PSLOLP为停电概率，当不考虑动态切负荷时，认为该指标即为系统停电概率；P(Xi)为系统状态的概率分布函数；P(Xi,k)为元件的概率分布函数；
[0089]
[0090] 式中，FSLOLF为停电频率，单位：次/年；当不考虑动态切负荷时，认为该指标即为系统停电频率；SLC为存在切负荷的系统状态集合；Mdown为系统中停运元件集合；Mup为系统中运行元件集合；μk为修复率；λk为故障率；
[0091]
[0092] 式中，PSEPNS为停电功率期望值，单位：MW；当不考虑动态切负荷时，认为该指标即为系统停电功率期望值；PS,i为第i个系统状态的停电电量，单位为：MWh；
[0093] ESEENS＝PSEPNS·T
[0094] 式中，ESEENS为停电电量期望值，单位：MWh/年；当不考虑动态切负荷时，认为该指标即为系统停电电量期望值；T为研究时段小时数，一年取8760小时；
[0095] HSLOLE＝PSLOLP·T
[0096] 式中，HSLOLE为年停电小时数，单位为：h/年；
[0097] DSLOLD＝HSLOLE/FSLOLF
[0098] 式中，DSLOLD均次停电时间，单位：h/次；当不考虑动态切负荷时，认为该指标即为系统停电时间期望值；
[0099] PSAPC＝PSEPNS/FSLOLF
[0100] 式中，PSAPC为均次停电功率，单位：MW/次；
[0101] ESAEC＝ESEENS/FSLOLF
[0102] 式中，ESAEC为均次停电电量，单位：MWh/次；
[0103] PSMPC＝max(PS,i)，i＝1,2,...,NLC
[0104] 式中，PSMPC为最大停电功率，单位：MW；NLC为存在切负荷的系统状态的总个数；
[0105] ESMEC＝max(ES,i)，i＝1,2,...,NLC
[0106] 式中，ESMEC为最大停电电量，单位：MWh；ES,i为第i个系统状态的停电电量，单位：MWh；
[0107] ISISI＝ESEENS/Lmax×60
[0108] 式中，ISISI为停电严重性指标，单位：分钟；Lmax为系统年峰荷；
[0109] 步骤4-2、计算负荷节点的可靠性水平的指标，公式如下：
[0110]
[0111] 式中，PVIO,l为负荷点l存在电压越界的概率，SVIO,l为负荷点l存在电压越界的系统状态集合；
[0112] 步骤4-3、计算系统可靠性指标在负荷节点间分配情况的平均值指标，公式如下：
[0113]
[0114] 式中，FSYAIF为负荷点平均停电频率，单位：次/(年.负荷点)；NL为系统中负荷点个数；
[0115]
[0116] 式中，HSYAIH为负荷点年平均停电时数，单位：h/(年.负荷点)；
[0117]
[0118] 式中，DSYAID为负荷点平均停电持续时间，单位：h/(次.负荷点)；
[0119] PSYAPC＝PSEPNS/NL
[0120] 式中，PSYAPC为负荷点平均停电功率，单位：MW/(年.负荷点)；
[0121] ESYAEC＝ESEENS/NL
[0122] 式中，ESYAEC为负荷点平均停电电量，单位：MWh/(年.负荷点)；
[0123] 优选的，所述步骤(9)中，计算元件可靠性指标，公式如下：
[0124]
[0125] 式中，POL,k为线路k的过负荷概率；SOL,k为线路k存在过载的系统状态集合；
[0126] ΔPmax,k＝max(ΔPk,j)，(j＝1,2,...,NOL,k)
[0127] 式中，ΔPk,j为线路k过负荷量，NOL,k为线路k存在过载的系统状态的总个数，ΔPmax,k为线路k的有功最大过负荷量。
[0128] 如图3所示，某地区电网包含24台发电机，140条输电线路，104台变压器(其中三绕组变有37台)，53个负荷节点，包含有220kV、110kV这两层电压等级。根据基于静态可靠性概率指标的电网可靠性水平判定方法计算得到各可靠性如下表所示。
[0129] 表1 系统可靠性指标
[0130]PSLOLP PSEPNS(MW) ESEENS(MWh/y)
0.00239169 0.03809812 333.7395
HSLOLE(h/y) POL,k PSMPC(MW)
[0131]20.95128 0.00136305 270
[0132] 表2 元件可靠性指标(按线路过负荷概率由大到小)
[0133]线路名称 POL,k ΔPmax,k(MW)
BUS33–BUS27 4.662072e-004 5.450000e+001
BUS8–BUS30 1.351481e-004 9.710000e+001
BUS36–BUS19 2.106149e-007 4.520000e+001
BUS33–BUS36 2.009073e-007 6.660000e+001
BUS30–BUS32 1.100947e-007 1.140000e+001
[0134] 表3 负荷点可靠性指标
[0135]负荷点年平均停电时数 负荷点平均停电功率 负荷点平均停电电量
1.812997e-001 1.788644e-004 1.566852
[0136] 表4 负荷节点可靠性指标(按负荷点停电概率由大到小)
[0137]节点名称 PSLOLP PSEPNS(MW)
BUS27 5.765937e-004 1.042017e-003
BUS9 5.393999e-004 1.519624e-002
BUS15 1.351174e-004 2.187986e-003
BUS38 1.351080e-004 1.418368e-004
BUS1 1.351080e-004 1.620992e-004
[0138] 以单一元件故障(N-1)是否造成系统失稳作为衡量系统稳定性的标准，阈值为1e-003，上述元件的可靠性参数均小于阈值，由此可见系统符合安全标准。
[0139] 基于静态可靠性概率指标的电网可靠性水平判定方法计算比较方便，概念清晰，在实际电网中得到了应用和验证。
[0140] 最后应当说明的是：以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者等同替换，而未脱离本发明精神和范围的任何修改或者等同替换，其均应涵盖在本发明的权利要求范围当中。
