CN110570079B - 一种并列供电系统多种运行状态联合概率的计算方法 - Google Patents

一种并列供电系统多种运行状态联合概率的计算方法 Download PDF

Info

Publication number
CN110570079B
CN110570079B CN201910660579.7A CN201910660579A CN110570079B CN 110570079 B CN110570079 B CN 110570079B CN 201910660579 A CN201910660579 A CN 201910660579A CN 110570079 B CN110570079 B CN 110570079B
Authority
CN
China
Prior art keywords
state
shutdown
probability
rate
outage
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910660579.7A
Other languages
English (en)
Other versions
CN110570079A (zh
Inventor
余方明
吴杰康
何家裕
庄仲
梁继深
杨金文
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong University of Technology
Original Assignee
Guangdong University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN201910660579.7A priority Critical patent/CN110570079B/zh
Publication of CN110570079A publication Critical patent/CN110570079A/zh
Application granted granted Critical
Publication of CN110570079B publication Critical patent/CN110570079B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
    • Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Educational Administration (AREA)
  • Marketing (AREA)
  • Development Economics (AREA)
  • Health & Medical Sciences (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种并列供电系统多种运行状态联合概率的计算方法，包括下述步骤：从电网能量管理系统EMS获取发电机运行数据，采用蒙特卡罗模拟方法确定一定运行周期内发电机G的运行状态；从电网能量管理系统EMS获取线路运行数据，采用蒙特卡罗模拟方法确定一定运行周期内第i条线路Li的运行状态以及整个n条线路的运行状态，包括全部运行或全部停运两种状态；本发明可以计算出在一定运行周期内并列供电系统发电机所处状态的概率，并列运行的线路所处状态的概率，并列运行的变压器所处状态的概率，再利用蒙特卡罗模拟法来确定系统各元件的状态，为确定供电系统状态提供技术方法，为电网可靠运行提供技术支撑。

Description

一种并列供电系统多种运行状态联合概率的计算方法
技术领域
本发明涉及电力系统及其自动化技术领域，具体涉及一种并列供电系统多种运行状态联合概率的计算方法。
背景技术
电力系统是一个动态的时变系统，其状态在各个运行时刻都在变化，而且电力系统的运行充满了不确定性。电力系统的状态表示系统在特定的时间里所处的特定状况，这些状况包括运行、维修、停运或其他状况。系统中的每一个元件可以取多种状态，只有确定了系统中每一个元件的状态，才能确定整个系统的状态。
在电力系统中，元件的停运通常可分为独立停运和相关停运两类。在考虑发电机和线路的运行状态时，主要涵盖：强迫停运状态、计划停运状态、降额状态。强迫停运随机发生，而非人力所能控制；计划停运(例如元件的维修或更换)则是人为安排，而不是由失效引起的；失效不是很严重时，一些元件(如：发电机组和高压直流输电线路)还可在降额状态下运行，这被称为部分失效模式。实际上还有一些停运介于二者之间，如变压器严重漏油不需立刻强迫停运，但必须在有限的时间内退出运行，这类停运称为半强迫停运。因而在考虑变压器的运行状态时，主要涵盖：强迫停运状态、计划停运状态、半强迫停运状态、降额状态。半强迫停运是指系统元件故障引起的一种可延迟停运，这种停运与强迫因素相关，而不可能预安排。此外，任何恶劣气候条件或暴风雨灾害都将大大增加所有暴露在外的元件的失效的可能性。
由于一个停运状态包含有一个以上元件的失效，涉及多个元件的相关停运事件的概率远大于这些相同元件同时独立停运的概率。因而，在以上各元件状态模型中引入一个共因停运状态。共因停运是指由于同一原因引起的多个元件的同时停运。例如变电站终端设备的任何失效都可能导致多条线路停运。
针对电力系统多种运行状态联合概率的计算方法，以往通常采用潮流计算的方法来加以判断和确认，近年来又引入了一种马尔可夫状态评估方法，但是这些方法一般是使用简单的两状态模型(运行和强迫停运)，无法计及电网中各元件的多种运行状态，不能确定供电系统中多种运行状态的联合概率，更不能确定供电系统中多种运行状态的类型、元件个数及损失。
发明内容
本发明的目的在于克服现有技术的缺点与不足，提供一种并列供电系统多种运行状态联合概率的计算方法，该方法能够为确定供电系统状态提供技术方法，为电网可靠运行提供技术支撑。
本发明的目的通过下述技术方案实现：
一种并列供电系统多种运行状态联合概率的计算方法，包括下述步骤：
S1，从电网能量管理系统EMS获取发电机运行数据，采用蒙特卡罗模拟方法确定一定运行周期内发电机G的运行状态，确定参数
Figure BDA0002138390920000021
Figure BDA0002138390920000031
S2，从电网能量管理系统EMS获取线路运行数据，采用蒙特卡罗模拟方法确定一定运行周期内第i条线路Li的运行状态以及整个n条线路的运行状态，包括全部运行或全部停运两种状态，1≤n≤NL，确定参数
Figure BDA0002138390920000032
S3，从电网能量管理系统EMS获取变压器运行数据，采用蒙特卡罗模拟方法确定一定运行周期内第i台变压器Ti的运行状态以及整个n台变压器的运行状态，包括全部运行或全部停运两种状态，1≤n≤NT，确定参数
Figure BDA0002138390920000033
采用序贯蒙特卡罗模拟法确定一定运行周期内第i台变压器Ti处于半强迫停运状态的次数，确定参数
Figure BDA0002138390920000034
S4，计算并列供电系统发电机多种运行状态的概率计算方法，计算公式为：
Figure BDA0002138390920000035
式中
Figure BDA0002138390920000036
为发电机处于运行状态的概率，
Figure BDA0002138390920000037
为发电机处于停运状态的概率；k是只有0或1的二值变量，如果发电机处于运行状态，其值为1，反之其值为0；
S5，计算并列供电系统n条线路其它多种运行状态和共因停运的组合状态所处的概率，计算公式为：
Figure BDA0002138390920000038
式中LSU和LSD分别为这n条线路组合状态中线路处于运行状态和线路处于停运状态的元件的两个集合，其总数加起来等于n；
Figure BDA0002138390920000039
为第i条线路处于运行状态的概率，
Figure BDA0002138390920000041
为第i条线路处于停运状态的概率；k是只有0或1的二值变量，如果这n条线路发生共因停运，其值为1，反之其值为0；pLcD(n)和pLcU(n)分别表示这n条线路共因停运发生和不发生的概率；
Figure BDA0002138390920000042
为有n条线路发生共因停运并使这些线路全部退出运行时共因停运失效率
Figure BDA0002138390920000043
大于其它多种运行状态下失效率
Figure BDA0002138390920000044
之和的概率；
S6，计算并列供电系统n台变压器其它多种运行状态和共因停运的组合状态所处的概率，计算公式为：
Figure BDA0002138390920000045
式中TSU、TSD、TBU、TBO分别为这n台变压器组合状态中处于运行状态、处于停运状态、处于半强迫停运状态中的运行和预停运状态和处于半强迫停运状态中的停运状态的元件的四个集合，其总数加起来等于n；
Figure BDA0002138390920000046
为第i台变压器处于运行状态的概率，
Figure BDA0002138390920000047
为第i台变压器处于停运状态的概率，
Figure BDA0002138390920000048
和
Figure BDA0002138390920000049
为第i台变压器处于半强迫停运状态中的运行和预停运状态的概率，
Figure BDA00021383909200000410
为第i台变压器处于半强迫停运状态中的停运状态的概率；k是只有0或1的二值变量，如果这n台变压器发生共因停运，其值为1，反之其值为0；pTcD(n)和pTcU(n)分别表示这n台变压器共因停运发生和不发生的概率；
Figure BDA00021383909200000411
为有n台变压器发生共因停运并使这些变压器全部退出运行时共因停运失效率
Figure BDA00021383909200000412
大于其它多种运行状态和半强迫停运下失效率
Figure BDA00021383909200000413
之和的概率；
S7，计算并列供电系统多种运行状态联合概率，计算公式为：
Figure BDA0002138390920000051
优选地，所述S4中发电机处于运行状态和处于停运状态的概率
Figure BDA0002138390920000052
计算所需的各运行状态的失效率、修复率的参数从电网能量管理系统EMS获取其实时值，其计算公式分别为：
Figure BDA0002138390920000053
Figure BDA0002138390920000054
优选地，所述S5的具体步骤如下：
S5.1，第i条线路处于运行状态和处于停运状态的概率
Figure BDA0002138390920000055
以及这n条线路共因停运发生和不发生的概率pLcD(n)和pLcU(n)计算所需的各运行状态的失效率、修复率的参数，从电网能量管理系统EMS获取其实时值，其计算公式分别为：
Figure BDA0002138390920000056
Figure BDA0002138390920000057
Figure BDA0002138390920000058
S5.2，从电网能量管理系统EMS获取线路运行数据，采用蒙特卡罗模拟方法确定n条线路发生共因停运并使这些线路全部退出运行时共因停运失效率
Figure BDA0002138390920000059
大于其它多种运行状态下失效率
Figure BDA00021383909200000510
之和的概率
Figure BDA00021383909200000511
其计算公式为：
Figure BDA0002138390920000061
S5.3，计算并列供电系统n条线路其它多种运行状态和共因停运的组合状态所处的概率，其计算公式为：
Figure BDA0002138390920000062
优选地，所述S6的具体步骤如下：
S6.1，第i台变压器处于运行状态和处于停运状态以及处于半强迫停运状态中的运行和预停运状态和处于半强迫停运状态中的停运状态的概率
Figure BDA0002138390920000063
再加上这n台变压器共因停运发生和不发生的概率pTcD(n)和pTcU(n)，计算所需的各运行状态的失效率、修复率的参数从电网能量管理系统EMS获取其实时值，其计算公式分别为：
Figure BDA0002138390920000064
Figure BDA0002138390920000065
Figure BDA0002138390920000066
Figure BDA0002138390920000067
Figure BDA0002138390920000068
S6.2，从电网能量管理系统EMS获取变压器运行数据，采用蒙特卡罗模拟方法确定n台变压器发生共因停运并使这些变压器全部退出运行时，共因停运失效率
Figure BDA0002138390920000071
大于其它多种运行状态和半强迫停运下失效率
Figure BDA0002138390920000072
之和的概率
Figure BDA0002138390920000073
其计算公式为：
Figure BDA0002138390920000074
S6.3，计算并列供电系统n台变压器其它多种运行状态和共因停运的组合状态所处的概率，其计算公式为：
Figure BDA0002138390920000075
本发明与现有技术相比具有以下的有益效果：
本发明可以计算出在一定运行周期内并列供电系统发电机所处状态(考虑强迫停运、计划停运、部分失效状态)的概率，并列运行的线路所处状态(考虑强迫停运、计划停运、部分失效状态以及共因停运)的概率，并列运行的变压器所处状态(考虑强迫停运、半强迫停运、计划停运、部分失效状态以及共因停运)的概率，再利用蒙特卡罗模拟法来确定系统各元件的状态，为确定供电系统状态提供技术方法，为电网可靠运行提供技术支撑。
附图说明
图1为本发明所针对的并列供电系统组成和连接关系示意图；
图2为本发明所针对的发电机、线路、变压器强迫失效、计划停运和部分失效的状态空间图；
图3为本发明所针对的变压器半强迫停运的状态空间图；
图4为本发明所针对的并列供电系统n个元件共因停运的分离模型图；
图5为本发明的流程框图。
图中附图标记为：1、发电机；2、发电机机端母线；3、并列运行的第一条线路；4、并列运行的第NL条线路；5、变压器高压母线；6、并列运行的第一台变压器；7、并列运行的第NT台变压器；8、变压器低压母线；9、负荷。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。
本发明为确定供电系统发生多种运行状态类型、元件个数、损失和风险提供技术方法，其基本原理是同时考虑发电机、线路、变压器的随机故障，负荷为确定性负荷，所以只需对发电机、线路和变压器的状态进行抽样。通过电网能量管理系统EMS获取电网运行的数据，在考虑电网运行方式的不确定性时主要是引入发电机、变压器、线路等设备的不确定性运行状态。在考虑计划停运时，假设计划停运和恢复时间服从某一给定的分布，该分布的参数可从计划停运事件的统计数据进行估计；在考虑变压器半强迫停运时，假设对系统进行短期运行规划；在考虑共因停运时，假设其它多种运行状态和共因停运之间是独立的，且不互斥。最后根据给定的电力系统各元件的失效率和修复率等参数，可以计算出各元件的状态概率，再利用蒙特卡罗模拟法来确定系统的状态，为确定电网系统状态提供必要的技术支撑。
本发明方法针对图1所示的并列供电系统，并考虑发电机运行方式改变、并列线路运行方式改变、并列变压器运行方式改变的不确定性和随机性，提出并列供电系统多种运行状态联合概率的计算方法。
如图1所示，由1台发电机G1，NL条线路L1、L2、L3、…、LNL并列运行和NT台变压器T1、T2、T3、…、TNT并列运行组成供电系统，假设负荷功率为SD(SD＝PD+jQD)。这是一种由发电机、多线路、多变压器及负荷组成的并列供电系统。
如图2所示，给出发电机、线路、变压器强迫失效、计划停运和部分失效的状态空间图。图中λ、λp、λd分别为强迫停运、计划停运和降额状态的转移率(失效率)；μ、μp、μd分别为强迫停运、计划停运和降额状态的修复率(修复次数/年)。
变压器的半强迫停运是指变压器故障引起的一种可延迟停运，这种停运与强迫因素相关，而不可能预安排。由于变压器漏油一般不会导致立即失效，但需要在一定的时间内停运。停运的延迟时间一般取决于如下几个因素：漏油的严重程度，可用检修能力，通过运行措施减少漏油的可能性，以及漏油对系统安全和对环境的影响等。
如图3所示，给出了变压器半强迫停运的状态空间图，假设每台变压器的半强迫停运状态和图2中其它多种运行状态之间是独立的，且不互斥。
图3中，λs为从运行向预停运状态的转移率(失效率)；μs为修复率，即修复时间的倒数；μso为从预停运状态向停运状态的转移率，即可延迟时间的倒数。
在此并列供电系统中，可能涉及多个元件的相关停运事件，线路和变压器的共因停运是指由于同一原因引起的多条线路或多台变压器的同时停运。在考虑共因停运时，假设每条线路或每台变压器的多种运行状态(包括变压器的半强迫停运状态)和共因停运之间是独立的，且不互斥。故应用分离模型来定义线路和变压器的多种运行状态和共因停运的组合状态。
如图4所示，给出了共因停运的分离模型，图2和图3相当于是每条线路或每台变压器的多种运行状态和每台变压器半强迫停运状态的分离模型。λc为全部运行状态向全部停运状态的转移率(失效率)，μc为共因停运的修复率。
图5中的步骤1描述发电机多种运行状态的参数确定及其参数计算的过程和方法。从电网能量管理系统EMS获取发电机的运行数据，按照抽取1年(1个月、1个季度作为每一个时段)的数据规模进行处理、计算和分析。重点抽取发电机强迫停运次数、计划停运次数、部分失效频率和各修复时间及其发生的条件等数据，然后通过计算得到发电机多种运行状态的各失效率和修复率。具体来说，采用蒙特卡罗模拟方法确定一定运行周期内发电机G的运行状态，确定参数
Figure BDA0002138390920000111
Figure BDA0002138390920000112
图5中的步骤2描述线路多种运行状态的参数确定及其参数计算的过程和方法。从电网能量管理系统EMS获取并列运行线路的运行数据，按照抽取1年(1个月、1个季度作为每一个时段)的数据规模进行处理、计算和分析。重点抽取线路强迫停运次数、计划停运次数、部分失效频率以及共因停运频率和各修复时间及其发生的条件等数据，然后通过计算得到线路多种运行状态(包括共因停运状态)的各失效率和修复率。具体来说，采用蒙特卡罗模拟方法确定一定运行周期内第i条线路Li的运行状态以及整个n条线路的运行状态，包括全部运行或全部停运两种状态，1＜n≤NL，确定参数
Figure BDA0002138390920000113
Figure BDA0002138390920000114
图5中的步骤3描述变压器多种运行状态的参数确定及其参数计算的过程和方法。从电网能量管理系统EMS获取并列运行变压器的运行数据，按照抽取1年(1个月、1个季度作为每一个时段)的数据规模进行处理、计算和分析。重点抽取变压器强迫停运次数、半强迫停运失效频率、计划停运次数、部分失效频率以及共因停运频率和各修复时间、半强迫停运可延迟时间及其发生的条件等数据，然后通过计算得到变压器多种运行状态(包括半强迫停运状态、共因停运状态)的各失效率和修复率。具体来说，采用蒙特卡罗模拟方法确定一定运行周期内第i台变压器Ti的运行状态以及整个n台变压器的运行状态，包括全部运行或全部停运两种状态，1＜n≤NT，确定参数
Figure BDA0002138390920000115
Figure BDA0002138390920000121
采用序贯蒙特卡罗模拟法确定一定运行周期内第i台变压器Ti处于半强迫停运状态(包括可延迟时间长度的各种因素)的次数，确定参数
Figure BDA0002138390920000122
图5中的步骤4描述并列供电系统发电机多种运行状态的概率的计算过程和方法。首先确定发电机处于运行状态(考虑部分失效模式)和处于停运状态(考虑强迫停运、计划停运)的概率，发电机处于运行状态(考虑部分失效模式)和处于停运状态(考虑强迫停运、计划停运)的概率
Figure BDA0002138390920000123
计算所需的各运行状态的失效率、修复率等参数从电网能量管理系统EMS获取其实时值，其计算公式分别为：
Figure BDA0002138390920000124
Figure BDA0002138390920000125
在发电机处于运行状态(考虑部分失效模式)和处于停运状态(考虑强迫停运、计划停运)的概率
Figure BDA0002138390920000126
的基础上，并列供电系统发电机所处状态的概率按照如下公式计算：
Figure BDA0002138390920000127
图5中的步骤5描述并列供电系统n条线路其它多种运行状态和共因停运的组合状态所处概率的计算过程和方法。具体步骤如下：
(1)首先确定第i条线路处于运行状态(考虑部分失效模式)和处于停运状态(考虑强迫停运、计划停运)的概率以及这n条线路共因停运发生和不发生的概率，第i条线路处于运行状态(考虑部分失效模式)和处于停运状态(考虑强迫停运、计划停运)的概率
Figure BDA0002138390920000131
以及这n条线路共因停运发生和不发生的概率pLcD(n)和pLcU(n)计算所需的各运行状态的失效率、修复率等参数从电网能量管理系统EMS获取其实时值，其计算公式分别为：
Figure BDA0002138390920000132
Figure BDA0002138390920000133
Figure BDA0002138390920000134
(2)从电网能量管理系统EMS获取线路运行数据，采用蒙特卡罗模拟方法确定n条线路发生共因停运并使这些线路全部退出运行时共因停运失效率
Figure BDA0002138390920000135
大于其它多种运行状态下失效率
Figure BDA0002138390920000136
之和的概率
Figure BDA0002138390920000137
其计算公式为：
Figure BDA0002138390920000138
(3)计算并列供电系统n条线路其它多种运行状态和共因停运的组合状态所处的概率，其计算公式为：
Figure BDA0002138390920000139
图5中的步骤6描述并列供电系统n台变压器其它多种运行状态(含半强迫停运)和共因停运的组合状态所处概率的计算过程和方法。具体步骤如下：
(1)首先确定第i台变压器处于运行状态(考虑部分失效模式)和处于停运状态(考虑强迫停运、计划停运)以及处于半强迫停运状态中的运行和预停运状态和处于半强迫停运状态中的停运状态的概率再加上这n台变压器共因停运发生和不发生的概率，第i台变压器处于运行状态(考虑部分失效模式)和处于停运状态(考虑强迫停运、计划停运)以及处于半强迫停运状态中的运行和预停运状态和处于半强迫停运状态中的停运状态的概率
Figure BDA0002138390920000141
再加上这n台变压器共因停运发生和不发生的概率pTcD(n)和pTcU(n)计算所需的各运行状态的失效率、修复率等参数从电网能量管理系统EMS获取其实时值，其计算公式分别为：
Figure BDA0002138390920000142
Figure BDA0002138390920000143
Figure BDA0002138390920000144
Figure BDA0002138390920000145
Figure BDA0002138390920000146
(2)从电网能量管理系统EMS获取变压器运行数据，采用蒙特卡罗模拟方法确定n台变压器发生共因停运并使这些变压器全部退出运行时共因停运失效率
Figure BDA0002138390920000147
大于其它多种运行状态和半强迫停运下失效率
Figure BDA0002138390920000148
之和的概率
Figure BDA0002138390920000149
其计算公式为：
Figure BDA0002138390920000151
(3)计算并列供电系统n台变压器其它多种运行状态(含半强迫停运)和共因停运的组合状态所处的概率，其计算公式为：
Figure BDA0002138390920000152
图5中的步骤7描述并列供电系统多种运行状态联合概率的计算过程和方法。并列供电系统多种运行状态联合概率的计算公式为：
Figure BDA0002138390920000153
并列供电系统多种运行状态联合概率的计算方法，与一定运行周期内发电机、线路和变压器运行、强迫停运、计划停运、部分失效状态有关系；与一定运行周期内变压器漏油导致的半强迫停运状态有关系；与一定运行周期内线路、变压器发生共因停运的元件个数有关系。
本发明可以计算出在一定运行周期内并列供电系统发电机所处状态(考虑强迫停运、计划停运、部分失效状态)的概率，并列运行的线路所处状态(考虑强迫停运、计划停运、部分失效状态以及共因停运)的概率，并列运行的变压器所处状态(考虑强迫停运、半强迫停运、计划停运、部分失效状态以及共因停运)的概率，再利用蒙特卡罗模拟法来确定系统各元件的状态，为确定供电系统状态提供技术方法，为电网可靠运行提供技术支撑。
上述为本发明较佳的实施方式，但本发明的实施方式并不受上述内容的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

Claims (4)

1.一种并列供电系统多种运行状态联合概率的计算方法，其特征在于，包括下述步骤：
S1，从电网能量管理系统EMS获取发电机运行数据，采用蒙特卡罗模拟方法确定一定运行周期内发电机G的运行状态，确定参数
Figure FDA0003952664510000011
Figure FDA0003952664510000012
S2，从电网能量管理系统EMS获取线路运行数据，采用蒙特卡罗模拟方法确定一定运行周期内第i条线路Li的运行状态以及整个n条线路的运行状态，包括全部运行或全部停运两种状态，1≤n≤NL，确定参数
Figure FDA0003952664510000013
S3，从电网能量管理系统EMS获取变压器运行数据，采用蒙特卡罗模拟方法确定一定运行周期内第i台变压器Ti的运行状态以及整个n台变压器的运行状态，包括全部运行或全部停运两种状态，1≤n≤NT，确定参数
Figure FDA0003952664510000014
采用序贯蒙特卡罗模拟法确定一定运行周期内第i台变压器Ti处于半强迫停运状态的次数，确定参数
Figure FDA0003952664510000015
S4，计算并列供电系统发电机多种运行状态的概率计算方法，计算公式为：
Figure FDA0003952664510000016
式中
Figure FDA0003952664510000017
为发电机处于运行状态的概率，
Figure FDA0003952664510000018
为发电机处于停运状态的概率；k是只有0或1的二值变量，如果发电机处于运行状态，其值为1，反之其值为0；
S5，计算并列供电系统n条线路其它多种运行状态和共因停运的组合状态所处的概率，计算公式为：
Figure FDA0003952664510000021
式中LSU和LSD分别为这n条线路组合状态中线路处于运行状态和线路处于停运状态的元件的两个集合，其总数加起来等于n；
Figure FDA0003952664510000022
为第i条线路处于运行状态的概率，
Figure FDA0003952664510000023
为第i条线路处于停运状态的概率；k是只有0或1的二值变量，如果这n条线路发生共因停运，其值为1，反之其值为0；pLcD(n)和pLcU(n)分别表示这n条线路共因停运发生和不发生的概率；
Figure FDA0003952664510000024
为有n条线路发生共因停运并使这些线路全部退出运行时共因停运失效率
Figure FDA0003952664510000025
大于其它多种运行状态下失效率
Figure FDA0003952664510000026
之和的概率；
S6，计算并列供电系统n台变压器其它多种运行状态和共因停运的组合状态所处的概率，计算公式为：
Figure FDA0003952664510000027
式中TSU、TSD、TBU、TBO分别为这n台变压器组合状态中处于运行状态、处于停运状态、处于半强迫停运状态中的运行和预停运状态和处于半强迫停运状态中的停运状态的元件的四个集合，其总数加起来等于n；
Figure FDA0003952664510000028
为第i台变压器处于运行状态的概率，
Figure FDA0003952664510000029
为第i台变压器处于停运状态的概率，
Figure FDA00039526645100000210
和
Figure FDA00039526645100000211
为第i台变压器处于半强迫停运状态中的运行和预停运状态的概率，
Figure FDA00039526645100000212
为第i台变压器处于半强迫停运状态中的停运状态的概率；k是只有0或1的二值变量，如果这n台变压器发生共因停运，其值为1，反之其值为0；pTcD(n)和pTcU(n)分别表示这n台变压器共因停运发生和不发生的概率；
Figure FDA00039526645100000213
为有n台变压器发生共因停运并使这些变压器全部退出运行时共因停运失效率
Figure FDA0003952664510000031
大于其它多种运行状态和半强迫停运下失效率
Figure FDA0003952664510000032
之和的概率；
S7，计算并列供电系统多种运行状态联合概率，计算公式为：
Figure FDA0003952664510000033
其中，以下参数定义为：
Figure FDA0003952664510000034
一定运行周期内发电机G强迫停运的转移率，失效率；
Figure FDA0003952664510000035
一定运行周期内发电机G计划停运的转移率，失效率；
Figure FDA0003952664510000036
一定运行周期内发电机G降额状态的转移率，失效率；
Figure FDA0003952664510000037
一定运行周期内发电机G强迫停运的修复率；
Figure FDA0003952664510000038
一定运行周期内发电机G计划停运的修复率；
Figure FDA0003952664510000039
一定运行周期内发电机G降额状态的修复率；
Figure FDA00039526645100000310
一定运行周期内第i条线路Li强迫停运的转移率，失效率；
Figure FDA00039526645100000311
一定运行周期内第i条线路Li计划停运的转移率，失效率；
Figure FDA00039526645100000312
一定运行周期内第i条线路Li降额状态的转移率，失效率；
Figure FDA00039526645100000313
一定运行周期内第i条线路Li强迫停运的修复率；
Figure FDA00039526645100000314
一定运行周期内第i条线路Li计划停运的修复率；
Figure FDA00039526645100000315
一定运行周期内第i条线路Li降额状态的修复率；
Figure FDA00039526645100000316
一定运行周期内整个n条线路全部运行状态向全部停运状态的转移率，失效率；
Figure FDA00039526645100000317
一定运行周期内整个n条线路共因停运的修复率；
Figure FDA00039526645100000318
一定运行周期内第i台变压器Ti强迫停运的转移率，失效率；
Figure FDA00039526645100000319
一定运行周期内第i台变压器Ti计划停运的转移率，失效率；
Figure FDA00039526645100000320
一定运行周期内第i台变压器Ti降额状态的转移率，失效率；
Figure FDA0003952664510000041
一定运行周期内第i台变压器Ti强迫停运的修复率；
Figure FDA0003952664510000042
一定运行周期内第i台变压器Ti计划停运的修复率；
Figure FDA0003952664510000043
一定运行周期内第i台变压器Ti降额状态的修复率；
Figure FDA0003952664510000044
一定运行周期内整个n台变压器全部运行状态向全部停运状态的转移率，失效率；
Figure FDA0003952664510000045
一定运行周期内整个n台变压器共因停运的修复率；
Figure FDA0003952664510000046
一定运行周期内第i台变压器Ti从运行向预停运状态的转移率，失效率；
Figure FDA0003952664510000047
一定运行周期内第i台变压器Ti的修复率，即修复时间的倒数；
Figure FDA0003952664510000048
一定运行周期内第i台变压器Ti从预停运状态向停运状态的转移率，失效率，即可延迟时间的倒数。
2.根据权利要求1所述的并列供电系统多种运行状态联合概率的计算方法，其特征在于，所述S4中发电机处于运行状态和处于停运状态的概率
Figure FDA0003952664510000049
计算所需的各运行状态的失效率、修复率的参数从电网能量管理系统EMS获取其实时值，其计算公式分别为：
Figure FDA00039526645100000410
Figure FDA00039526645100000411
3.根据权利要求1所述的并列供电系统多种运行状态联合概率的计算方法，其特征在于，所述S5的具体步骤如下：
S5.1，第i条线路处于运行状态和处于停运状态的概率
Figure FDA00039526645100000412
以及这n条线路共因停运发生和不发生的概率pLcD(n)和pLcU(n)计算所需的各运行状态的失效率、修复率的参数，从电网能量管理系统EMS获取其实时值，其计算公式分别为：
Figure FDA0003952664510000051
Figure FDA0003952664510000052
Figure FDA0003952664510000053
S5.2，从电网能量管理系统EMS获取线路运行数据，采用蒙特卡罗模拟方法确定n条线路发生共因停运并使这些线路全部退出运行时共因停运失效率
Figure FDA0003952664510000054
大于其它多种运行状态下失效率
Figure FDA0003952664510000055
之和的概率
Figure FDA0003952664510000056
其计算公式为：
Figure FDA0003952664510000057
S5.3，计算并列供电系统n条线路其它多种运行状态和共因停运的组合状态所处的概率，其计算公式为：
Figure FDA0003952664510000058
4.根据权利要求1所述的并列供电系统多种运行状态联合概率的计算方法，其特征在于，所述S6的具体步骤如下：
S6.1，第i台变压器处于运行状态和处于停运状态以及处于半强迫停运状态中的运行和预停运状态和处于半强迫停运状态中的停运状态的概率
Figure FDA0003952664510000061
再加上这n台变压器共因停运发生和不发生的概率pTcD(n)和pTcU(n)，计算所需的各运行状态的失效率、修复率的参数从电网能量管理系统EMS获取其实时值，其计算公式分别为：
Figure FDA0003952664510000062
Figure FDA0003952664510000063
Figure FDA0003952664510000064
Figure FDA0003952664510000065
Figure FDA0003952664510000066
S6.2，从电网能量管理系统EMS获取变压器运行数据，采用蒙特卡罗模拟方法确定n台变压器发生共因停运并使这些变压器全部退出运行时，共因停运失效率
Figure FDA0003952664510000067
大于其它多种运行状态和半强迫停运下失效率
Figure FDA0003952664510000068
之和的概率
Figure FDA0003952664510000069
其计算公式为：
Figure FDA00039526645100000610
S6.3，计算并列供电系统n台变压器其它多种运行状态和共因停运的组合状态所处的概率，其计算公式为：
Figure FDA00039526645100000611
Figure FDA0003952664510000071
CN201910660579.7A 2019-07-22 2019-07-22 一种并列供电系统多种运行状态联合概率的计算方法 Active CN110570079B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910660579.7A CN110570079B (zh) 2019-07-22 2019-07-22 一种并列供电系统多种运行状态联合概率的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910660579.7A CN110570079B (zh) 2019-07-22 2019-07-22 一种并列供电系统多种运行状态联合概率的计算方法

Publications (2)

Publication Number Publication Date
CN110570079A CN110570079A (zh) 2019-12-13
CN110570079B true CN110570079B (zh) 2023-02-03

Family

ID=68773192

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910660579.7A Active CN110570079B (zh) 2019-07-22 2019-07-22 一种并列供电系统多种运行状态联合概率的计算方法

Country Status (1)

Country Link
CN (1) CN110570079B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103646350A (zh) * 2013-10-25 2014-03-19 南昌大学 一种基于风险理论的电网运行状态评估方法
CN103793853A (zh) * 2014-01-21 2014-05-14 中国南方电网有限责任公司超高压输电公司检修试验中心 基于双向贝叶斯网络的架空输电线路运行状态评估方法
CN105096039A (zh) * 2015-07-17 2015-11-25 广东电网有限责任公司东莞供电局 一种多线路并列供电系统过负荷停运概率计算方法
CN107086582A (zh) * 2017-05-18 2017-08-22 国网江苏省电力公司经济技术研究院 一种电网接入统一潮流控制器的可靠性分析方法
CN108512226A (zh) * 2018-05-21 2018-09-07 西安交通大学 一种灾害下电力系统恢复力评估的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103646350A (zh) * 2013-10-25 2014-03-19 南昌大学 一种基于风险理论的电网运行状态评估方法
CN103793853A (zh) * 2014-01-21 2014-05-14 中国南方电网有限责任公司超高压输电公司检修试验中心 基于双向贝叶斯网络的架空输电线路运行状态评估方法
CN105096039A (zh) * 2015-07-17 2015-11-25 广东电网有限责任公司东莞供电局 一种多线路并列供电系统过负荷停运概率计算方法
CN107086582A (zh) * 2017-05-18 2017-08-22 国网江苏省电力公司经济技术研究院 一种电网接入统一潮流控制器的可靠性分析方法
CN108512226A (zh) * 2018-05-21 2018-09-07 西安交通大学 一种灾害下电力系统恢复力评估的方法

Also Published As

Publication number Publication date
CN110570079A (zh) 2019-12-13

Similar Documents

Publication Publication Date Title
Balijepalli et al. Modeling and analysis of distribution reliability indices
Rios et al. Value of security: modeling time-dependent phenomena and weather conditions
Milanovic et al. Probabilistic assessment of financial losses due to interruptions and voltage sags-part I: the methodology
Fotuhi-Firuzabad et al. An analytical method to consider DG impacts on distribution system reliability
Ge et al. Reliability and maintainability improvement of substations with aging infrastructure
Sang et al. An integrated preventive operation framework for power systems during hurricanes
Carpaneto et al. Evaluation of the probability density functions of distribution system reliability indices with a characteristic functions-based approach
Sun et al. Probabilistic reliability evaluation for distribution systems with DER and microgrids
CN109447418B (zh) 一种输变电两级多元件并列的电网检修停运概率计算方法
CN110516287B (zh) 一种计及老化失效的供电系统停运概率计算方法
Papic et al. Comparing a transmission planning study of cascading with historical line outage data
Li et al. A probabilistic analysis approach to making decision on retirement of aged equipment in transmission systems
Jayaweera et al. Security of energy supply with change in weather conditions and dynamic thermal limits
CN114140266A (zh) 一种低压电网故障智能研判算法
Roslan et al. Sequential and nonsequential monte carlo in assessing reliability performance of distribution network
Anders et al. A comprehensive study of outage rates of air blast breakers
CN110570079B (zh) 一种并列供电系统多种运行状态联合概率的计算方法
Li et al. A risk evaluation based approach to replacement strategy of aged HVDC components and its application at BCTC
González-Fernández et al. Generating capacity reliability assessment of the Itaipu hydroelectric plant via sequential Monte Carlo simulation
Ndawula et al. A novel formulation of lv distribution network equivalents for reliability analysis
Enkhmunkh et al. A parameter estimation of Weibull distribution for reliability assessment with limited failure data
Watts et al. Classification and discussion on methods for cascading failure analysis in transmission system
Sailaja et al. Evaluation of reliability indices using FMEA technique
Yue et al. Enhanced probabilistic contingency analysis considering fast ramping events of renewable generation
Ekisheva et al. Analysis of CDM Events in NERC and WECC Power Systems Using TADS Outage Statistics

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant