CN109886834B - 一种基于最大风险分析的发输电系统可靠性跟踪算法 - Google Patents

一种基于最大风险分析的发输电系统可靠性跟踪算法 Download PDF

Info

Publication number
CN109886834B
CN109886834B CN201910072739.6A CN201910072739A CN109886834B CN 109886834 B CN109886834 B CN 109886834B CN 201910072739 A CN201910072739 A CN 201910072739A CN 109886834 B CN109886834 B CN 109886834B
Authority
CN
China
Prior art keywords
risk
state
fault
power generation
elements
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
CN201910072739.6A
Other languages
English (en)
Other versions
CN109886834A (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.)
Chongqing University
Original Assignee
Chongqing University
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 Chongqing University filed Critical Chongqing University
Priority to CN201910072739.6A priority Critical patent/CN109886834B/zh
Publication of CN109886834A publication Critical patent/CN109886834A/zh
Application granted granted Critical
Publication of CN109886834B publication Critical patent/CN109886834B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种基于最大风险分析的发输电系统可靠性跟踪算法，包括以下步骤：生成相似故障元素集合，设定系统风险阈值；应用最大风险分析模型求解；生成求解结果，并生成新的系统状态集合；进行风险分摊，更新各系统元件的风险贡献度；判断相似故障集合的风险值是否高于所设定的系统风险阈值，输出各系统元件风险贡献度，确定系统薄弱环节。本发明的有益效果有：可以有效地从系统状态集合内辨识出主要故障事件，结合相似故障状态缩减技术，可以在计及元件切负荷影响下通过可靠性跟踪快速地量化元件的风险贡献，确定对应的薄弱环节。

Description

一种基于最大风险分析的发输电系统可靠性跟踪算法
技术领域
本发明涉及电力系统可靠性评估领域，具体涉及一种基于最大风险分析的发输电系统可靠性跟踪算法。
背景技术
当前世界电力系统可靠性与社会经济、生活和安全都紧密相关。而目前的可靠性评估算法大多数针对整个系统可靠性评估设计，可以确定系统整体的可靠性水平。而如何衡量单个元件对系统可靠性的贡献，相关研究鲜有涉及。衡量不同元件对系统可靠性的贡献，可以有效辨识出当前系统的薄弱环节，即风险贡献度最高的元件。而薄弱环节可为系统的规划运行提供指导方向，进而对实现建设资金的最大利用率和改善系统的可靠性有着重要意义。
已有学者尝试分别从电力系统拓扑结构和网络攻击两个角度对系统的薄弱环节进行分析，但这两种方法均存在不足。基于电力系统拓扑结构的薄弱环节分析在故障导致系统网架结构或运行模式变化后可能需要重新确定系统的薄弱环节。网络攻击模型属于人为故障，即不考虑故障发生概率对风险分布及脆弱环节辨识结果的影响。在正常情况的日常运行或规划问题中通常不具备指导意义。此外，以上两类方法均无法量化分析各个元件对系统风险的贡献度。
借鉴潮流跟踪思想及其在输电网损分摊的应用，已有学者提出了以可靠性追踪方式将系统风险分摊给内部元件，可实现元件风险贡献的量化，并对系统的薄弱环节进行辨识。但是每一次可靠性跟踪均需与可靠性评估同时进行，因此可靠性评估过程的计算复杂度决定了可靠性跟踪计算的复杂度。目前可靠性评估算法主要可以分为状态解析法和蒙特卡诺模拟法，但是由于电力系统的元件维度极高，致使系统状态发生对应的“维数灾”问题；而且状态分析所涉及的潮流计算、最优切负荷等均为非线性问题，求解复杂度极高。即使采用线性化的状态分析模型，由于系统状态数的指数增长，仍会导致计算量爆炸。因此依靠状态解析法或蒙特卡诺模拟进行可靠性跟踪从根本上也面临着计算爆炸的问题。
实际上系统可靠性水平(如失负荷概率，期望失负荷量等)仅由很小的一部分系统状态所决定。状态解析法和蒙特卡诺模拟需要通过产生大量的系统状态来识别出这部分系统状态。然而这种状态生成过程具有明显的不足，即所生成的众多系统状态中，大多数的状态对于可靠性水平没有影响或影响很小，但会导致算法计算量大幅增加。
基于上述缺陷和不足，有必要提出一种方法，来提高大电网可靠性追踪计算效率。
发明内容
针对上述现有技术中的不足之处，本发明提供一种基于最大风险分析的发输电系统可靠性跟踪算法，其计算量小，提高了大电网可靠性追踪计算效率。
为了达到上述目的本发明采用了一下技术方案：
一种基于最大风险分析的发输电系统可靠性跟踪算法，其特征在于，包括以下步骤：
S1、根据发输电系统原始数据生成相似故障元素集合，并生成初始系统状态集合，并设定系统风险阈值；
S2、应用最大风险分析模型求解初始系统状态集合中风险最高的状态；
S3、生成S2步骤求解结果的相似故障集合，并从初始系统状态集合中删去该相似故障集合，以生成新的系统状态集合；
S4、应用风险分摊准则对已求得的系统故障事件集合进行风险分摊，更新各系统元件的风险贡献度；
S5、判断步骤S3中相似故障集合的风险值是否高于所设定的系统风险阈值，若是则返回步骤S2；如否进入下一步骤；
S6、输出各系统元件风险贡献度，确定系统薄弱环节。
进一步地，所述步骤S1中所述的发输电系统原始数据包括发输电系统的可靠性参数和电气参数；其根据发输电系统原始数据生成相似故障元素集合，并生成初始系统状态集合的具体步骤为，
S11、根据系统设备种类将系统元件分类；
S12、在同一种类的系统元件内，分别确定各自的相似指标；
S13、采用枚举法枚举产生所有可能的系统状态，由此生成初始系统状态集合。
进一步地，步骤S2中求解发输电系统的最大风险分析模型为直流潮流模型，其具体数学模型为
Figure BDA0001957784620000031
其中，中risk表示某一系统状态对应的系统风险值，将风险定义为切负荷量和对应状态概率的乘积，pg表示发电机组的故障概率、pl表示线路的故障概率；xg、xl分别表示发电机组、线路的状态变量，
Figure BDA0001957784620000032
表示系统最优切负荷运行状态所对应的各节点切负荷量；Ng为发电机组总数、Nl为线路总数、Ng为削负荷节点总数。
发电机组、线路的状态变量的取值
xg,xl＝{0,1}
取0代表该元件故障，取1代表该元件正常工作；
各元件状态确定的系统状态应属于系统可选状态集S，即
[xg1,...,xgNg,xl1,...,xlNl]∈S
取值函数为
Figure BDA0001957784620000033
式中，Pg、
Figure BDA0001957784620000034
Fl、θb分别表示系统状态确定后，系统最优切负荷运行状态所对应的机组出力、切负荷量、线路功率、相角的取值；
每个母线节点应满足的功率平衡
Figure BDA0001957784620000035
其中
Figure BDA0001957784620000041
为节点b处的负荷，G(b)表示连接至节点b处的发电机组集合，而Lfr(b)、Lto(b)分别表示将节点b作为潮流发送、接受节点的线路集合；λb为该约束条件对应的对偶变量；
直流潮流下的线路潮流方程
Figure BDA0001957784620000042
采用了大M法对Fl＝xlbl(θi-θj)进行了线性化处理,其中bl为线路l的电纳，M为一个足够大的正数；
Figure BDA0001957784620000043
和
Figure BDA0001957784620000044
表示该约束条件对应的对偶变量，其中“+”表示对应每个母线节点应满足的功率平衡右边的不等式约束，而“-”表示对应每个母线节点应满足的功率平衡左边的不等式约束；
Figure BDA0001957784620000045
表示线路功率约束条件对应的对偶变量；
机组的出力范围约束
Figure BDA0001957784620000046
式中，
Figure BDA0001957784620000047
表示该约束条件对应的对偶变量；
母线上负荷切除量的取值约束
Figure BDA0001957784620000048
式中，
Figure BDA0001957784620000049
表示该约束条件对应的对偶变量；
母线的相角取值约束，并指定1节点处为参考节点
Figure BDA0001957784620000051
式中，
Figure BDA0001957784620000052
表示该约束条件对应的对偶变量；
进一步地，将所述发输电系统的最大风险分析模型的双层模型，通过对偶理论转化为单层模型求解，其单层模型的数学式为
Figure BDA0001957784620000053
其中，原优化问题中下层问题所对应的强对偶定理条件为
Figure BDA0001957784620000054
对偶问题的约束条件为
Figure BDA0001957784620000055
Figure BDA0001957784620000056
Figure BDA0001957784620000057
Figure BDA0001957784620000061
Figure BDA0001957784620000062
Figure BDA0001957784620000063
Figure BDA0001957784620000064
Figure BDA0001957784620000065
进一步地，对所述发输电系统的最大风险分析模型的求解，还包括实现单层模型的线性化，其过程包括，
第一步，对目标函数，即单层模型的数学式进行线性化处理，首先进行对数变换，
Figure BDA0001957784620000066
并采用分段函数对式中的
Figure BDA0001957784620000067
进行线性化处理，
y＝min(yi)
Figure BDA0001957784620000068
Figure BDA0001957784620000069
第二步，原优化问题中下层问题所对应的强对偶定理条件中的双变量乘积进行线性化处理，首先引入辅助变量
Figure BDA0001957784620000071
同时采用线性约束实现
Figure BDA0001957784620000072
的等效线性化
Figure BDA0001957784620000073
Figure BDA0001957784620000074
Figure BDA0001957784620000075
Figure BDA0001957784620000076
其中Y、Z均为一个足够大的正数；
对线性化处理后的模型进行求解可即可确定当前系统状态集S内的风险最高的系统故障状态。
进一步地，步骤S3中，包括根据步骤S1生成的相似元件集合生成步骤S2中求解得到的当前最大风险系统状态的相似故障集合，并生成新的系统状态集合，具体步骤为，
S31、以S1步骤得到的相似故障元素集为单位将当前故障状态分解为片段，并生成各状态片段的相似故障状态片段；
S32、对各个相似元素集合的相似故障片段进行组合拼接，最终可生成当前系统状态的相似故障状态集合；其中相似故障状态集合内的状态个数N为
Figure BDA0001957784620000077
式中的ns表示相似元素集合个数，Ni为第i个相似故障元素集合的相似片段个数；
S33、在原始系统状态集合中删去由步骤S32得到的当前系统状态的相似故障状态集合，生成新的系统状态集合。
进一步地，步骤S4中包括应用风险分摊准则对已求得的系统故障事件集合进行风险分摊；其中风险分摊准则包括如果某两个系统状态A与B，A状态中比B状态新增了m个故障元件，除这m个故障元件外其余元件的状态均相同；而这两个系统状态的负荷损失相同，即这m个元件的故障没有令系统状况恶化，则这m个元件不参与A状态风险值的分担；
基于以上不可靠性跟踪准则，系统元件对发输电系统可靠性的贡献评估如下：
假设该发输电系统由n个元件Xi(i＝1,2…,n)组成，同时采用xi(i＝1,2…,n)表示各元件所对应的性能参数，f(x1,x2,…,xn)表示由系统状态为(X1,X2,…,Xn)时的系统风险值risk；
以X1和X2两个元件故障时所对应的系统状态为例，当假设元件参数xi(i＝1,2…,n)相互独立，仅与自身相关，f(x1,x2,…,xn)可表示为
f(x1,x2,...,xn)＝f(x1)f(x2)f(x3，...,xn)
其中f(x1)仅与X1有关，f(x2)仅与X2有关，f(x3,…,xn)仅与除X1、X2以外的元件有关；
风险分摊依据
Figure BDA0001957784620000081
其中的λi、μi分别为元件的故障率和修复率，ΔPd、ΔPi d分别表示当其他故障元件状态不变，而元件Xi处于故障和正常状态时系统切负荷量，
其中X1和X2各自对应风险分摊的部分f(s→1)和f(s→2)可表示为:
Figure BDA0001957784620000082
Figure BDA0001957784620000091
进一步地，步骤S4中按照比例分摊准则确定各元件风险贡献比之前，应该首先判断各元件的风险分摊系数；通过搜索比当前状态风险高的故障状态集，即在截止当前已经确定的解集进行判断。
进一步地，S6骤中的各系统元件风险贡献度由每次循环中步骤S4得到的风险贡献度累加得到，并根据元件风险贡献度高低选取风险贡献前十的元件作为系统薄弱环节集合。
本发明的有益效果有：可以有效地从系统状态集合内辨识出主要故障事件，结合相似故障状态缩减技术，可以在计及元件切负荷影响下通过可靠性跟踪快速地量化元件的风险贡献，确定对应的薄弱环节。能够为电力系统的调度检修、建设规划提供有效信息，也为电力市场中的负荷损失赔偿问题的定责提供了参考建议。
附图说明
图1为本发明可靠性跟踪算法流程图；
图2为本发明中对数函数分段线性化示意图；
图3为本发明中相似故障状态生成示意图；
图4为本发明中实施例使用的RBTS测试系统接线图。
具体实施方式
下面结合具体实施例及附图来进一步详细说明本实用新型。
一种如图1至4所示的一种基于最大风险分析的发输电系统可靠性跟踪算法，包括以下步骤：
S1：根据发输电系统原始数据生成相似故障元素集合，并生成初始系统状态集合以及设定系统风险阈值。
S2：应用最大风险分析模型求解当前初始系统状态集合中风险最高的状态。
S3：生成S2步骤求解结果的相似故障集合，并从初始系统状态集合中删去该相似故障集合，以生成新的系统状态集合。
S4：应用风险分摊准则对已求得的系统故障事件集合进行风险分摊，更新各系统元件的风险贡献度。
S5：判断S3似故障集合的风险值是否高于所设定的系统风险阈值。若S3中相似故障集合的风险值高于所设定的系统风险阈值则转S2，否则转S6。
S6：输出各系统原件风险贡献度，确定系统薄弱环节。
通过上述方案，可以有效地从系统状态集合内辨识出主要故障事件，结合相似故障状态缩减技术，可以在计及元件切负荷影响下通过可靠性跟踪快速地量化元件的风险贡献，确定对应的薄弱环节。能够为电力系统的调度检修、建设规划提供有效信息，也为电力市场中的负荷损失赔偿问题的定责提供了参考建议。
进一步的，在S1步骤中输入的数据包括发输电系统的可靠性参数和电气参数。根据上述数据，可以根据相似指标生成系统相似元素集合并生成初始系统状态集合，具体生成步骤为：
S11：根据系统设备种类如机组，线路，变压器等将系统元件进行分类。
S12：在同一种类的系统元件内，分别确定各自的相似指标，如发电机组可将装机量、机组所在节点、故障率作为指标，而输电线路或变压器可将额定容量、始末节点、故障率作为指标。根据相似指标将所有的系统元件分别分配至各自的相似元素集合中。
S13：采用枚举法枚举产生所有可能的系统状态，生成初始系统状态集合。
作为进一步说明，步骤S2中求解发输电系统的最大风险分析模型的目的是实现从系统状态集中识别出对可靠性水平影响最高的状态。该模型采用双层优化模型结构，其上层问题为风险值最大化，下层问题为最优切负荷问题对应的系统状态分析。在已有的可靠性研究中，直流潮流模型由于其采用线性结构，计算复杂度比交流潮流模型大幅降低，而且精确度可以满足可靠性评估要求，本发明在状态分析中同样也采取了直流潮流模型。其具体数学模型如式(1)-(11)所示。
Figure BDA0001957784620000101
约束条件
xg,xl＝{0,1} (2)
[xg1,...,xgNg,xl1,...,xlNl]∈S (3)
Figure BDA0001957784620000111
Figure BDA0001957784620000112
Figure BDA0001957784620000113
Figure BDA0001957784620000114
Figure BDA0001957784620000115
Figure BDA0001957784620000116
Figure BDA0001957784620000117
θ1＝0 (11)
其中，式(1)中risk表示某一系统状态对应的系统风险值，将风险定义为切负荷量和对应状态概率的乘积。其中的p表示对应元件的故障概率；
式(2)中xg、xl分别表示发电机组、线路的状态变量，取0代表该元件故障，取1代表该元件正常工作；
式(3)则表示由各元件状态确定的系统状态应属于系统可选状态集S。
式(4)中arg为取值函数，表示当上层模型中各元件的状态变量确定，即系统状态确定后，系统最优切负荷运行状态所对应的机组出力Pg、切负荷量
Figure BDA0001957784620000125
线路功率Fl、相角θb等变量取值；
式(5)表示每个母线节点应满足的功率平衡，其中Db为节点b处的负荷，G(b)表示连接至节点b处的发电机组集合，而Lfr(b)、Lto(b)分别表示将节点b作为潮流发送、接受节点的线路集合；
式(6)-(7)表示直流潮流下的线路潮流方程，并采用了大M法对Fl＝xlbl(θi-θj)进行了线性化处理,其中bl为线路l的电纳，M为一个足够大的正数；
式(8)表示机组的出力范围约束；
式(9)表示母线上负荷切除量的取值约束；
式(10)-(11)表示母线的相角取值约束，并指定1节点处为参考节点。
此外，式(5)-(11)中“:”后的λb、
Figure BDA0001957784620000121
和
Figure BDA0001957784620000122
等变量表示每个约束条件各自对应的对偶变量，
Figure BDA0001957784620000123
中“+”表示对应式(5)右边的不等式约束，而“-”表示对应左边的不等式约束。
大多数的优化算法均针对单层的优化模型设计，因此双层优化模型一般需转化为了单层优化模型进行求解。本发明采用对偶理论将该双层模型转化单层模型求解。转化后的单层模型如式(1)、(13)-(21)所示。
Figure BDA0001957784620000124
约束条件：s.t.
如(2)-(11)式
Figure BDA0001957784620000131
Figure BDA0001957784620000132
Figure BDA0001957784620000133
Figure BDA0001957784620000134
Figure BDA0001957784620000135
Figure BDA0001957784620000136
Figure BDA0001957784620000137
Figure BDA0001957784620000138
Figure BDA0001957784620000139
其中式(13)为原优化问题中下层问题所对应的强对偶定理条件，其含义为下层问题取最优值的充要条件；而式(14)-(21)为原下层模型对偶问题的约束条件。
为了实现模型的线性化，对模型进行以下两步处理：
第一步，对目标函数式(1)进行线性化处理。首先对式(1)进行对数变换，如式(22)所示。并采用如式(23)-(25)所示的分段函数对式(22)中的
Figure BDA0001957784620000141
进行线性化处理，分段过程如图2所示。
Figure BDA0001957784620000142
其中
y＝min(yi) (23)
Figure BDA0001957784620000143
Figure BDA0001957784620000144
第二步，对模型约束式(13)中形如
Figure BDA0001957784620000145
的双变量乘积进行线性化处理。首先引入辅助变量
Figure BDA0001957784620000146
同时采用线性约束(26)-(29)实现
Figure BDA0001957784620000147
的等效线性化,其中Y、Z均为一个足够大的正数。
Figure BDA0001957784620000148
Figure BDA0001957784620000149
Figure BDA00019577846200001410
Figure BDA00019577846200001411
对线性化处理后的模型进行求解可即可确定当前系统状态集S内的风险最高的系统故障状态。
进一步的，在步骤S3中，需要根据步骤S1生成的相似元件集合生成步骤S2中求解得到的当前最大风险系统状态的相似故障集合，并生成新的系统状态集合，具体步骤为：
S31：以S1步骤得到的相似故障元素集为单位将当前故障状态分解为片段，并生成各状态片段的相似故障状态片段。对于任意一个相似故障元素集而言，只要状态片段内的故障元件个数相等，则均为相似故障状态片段。
S32：对各个相似元素集合的相似故障片段进行组合拼接，最终可生成当前系统状态的相似故障状态集合。其生成示意图如图3所示。其中相似故障状态集合内的状态个数N可由式(30)计算，式中的ns表示相似元素集合个数，Ni为第i个相似故障元素集合的相似片段个数。
Figure BDA0001957784620000151
S33：在原始系统状态集合中删去由步骤S32得到的当前系统状态的相似故障状态集合，生成新的系统状态集合。
作为进一步说明，在步骤S4中需要应用风险分摊准则对已求得的系统故障事件集合进行风险分摊。在已有的可靠性跟踪理论中，由故障元件对风险指标按比例进行分摊，系统的不可靠性由全部故障元件承担。但现有理论未计及各元件故障所造成的负荷损失贡献，即可能存在以下情况：在一个固定的系统故障状态中，某个元件的故障不会引起系统失负荷，因此该元件不应分担本次失负荷责任。
针对上述问题，本发明提出新的可靠性准则对原有准则作为补充：
如果某两个系统状态A与B，A状态中比B状态新增了m个故障元件，除这m个故障元件外其余元件的状态均相同；而这两个系统状态的负荷损失相同，即这m个元件的故障没有令系统状况恶化，则这m个元件不参与A状态风险值的分担。
基于以上不可靠性跟踪准则，系统元件对发输电系统可靠性的“贡献”评估如下：
假设该发输电系统由n个元件Xi(i＝1,2…,n)组成，同时采用xi(i＝1,2…,n)表示各元件所对应的性能参数，f(x1,x2,…,xn)表示由系统状态为(X1,X2,…,Xn)时的系统风险值risk。
以X1和X2两个元件故障时所对应的系统状态为例，当假设元件参数xi(i＝1,2…,n)相互独立，仅与自身相关，f(x1,x2,…,xn)可表示为式(31)。其中f(x1)仅与X1有关，f(x2)仅与X2有关，f(x3,…,xn)仅与除X1、X2以外的元件有关。
f(x1,x2,...,xn)＝f(x1)f(x2)f(x3，...,xn) (31)
本发明采用式(32)作为新的风险分摊依据，其中的λi、μi分别为元件的故障率和修复率，ΔPd、ΔPi d分别表示当其他故障元件状态不变，而元件Xi处于故障和正常状态时系统切负荷量。式(32)的物理含义为在一次故障事件中，故障概率高和故障前后导致负荷损失变化较多的元件应承担较多的责任。
Figure BDA0001957784620000161
由所提出的可靠性跟踪准则可知，该系统状态风险应由故障元件X1和X2分担。其中X1和X2各自对应分摊的部分f(s→1)和f(s→2)可表示为:
Figure BDA0001957784620000162
Figure BDA0001957784620000163
作为进一步说明，需要特别注意的是，步骤S4中按照比例分摊准则确定各元件风险贡献比之前，应该首先判断各元件的风险分摊系数。而本方法已经根据风险对故障状态完成了排序，只需通过搜索比当前状态风险高的故障状态集，即在截止当前已经确定的解集进行判断即可，无需重新进行状态分析。这是因为对于一个k阶故障，如要根据准则确定这k个故障元件是否参与风险分摊，至少要与其他的k个k-1阶故障进行切负荷状况对比，而如果第i个元件不应参加风险分担，则对应的k-1阶故障与k阶故障切负荷状态相同，但k-1阶故障概率大于k阶故障，因此k-1阶风险大于k阶风险。换言之，如果第i个元件不需要参与分担该风险，则对应的k-1阶故障一定存在了之前求解确定的系统状态集内。
进一步的，S6骤中的各系统元件风险贡献度由每次循环中步骤S4得到的风险贡献度累加得到。根据元件风险贡献度高低选取风险“贡献”前十的元件作为系统薄弱环节集合。
本发明的一种实施例
结合图1可以看出，一种基于最大风险分析的发输电系统可靠性跟踪算法，按照以下步骤进行：
S1：根据发输电系统原始数据生成相似故障元素集合，并生成初始系统状态集合,设定系统风险阈值。
在本实施例中，本实例使用的算例为RBTS测试系统，该测试系统接线图如图4所示，其中包含11台发电机和9条输电线路，系统总装机容量240MW，系统年最大峰荷185MW，各元件可靠性参数和电气参数详见文献“Billinton R,Kumar S,Chowdhury N,et al.Areliability test system for educational purposes-basic data[J].IEEETransactions on Power Systems,1990,4(3):1238-1244.”
根据上述数据，可以根据相似指标生成系统相似元素集合并生成初始系统状态集合，具体生成步骤为：
S11：根据系统设备种类如机组，线路，变压器等将系统元件分类。
S12：在同一种类的系统元件内，分别确定各自的相似指标，根据相似指标将所有的系统元件分别分配至各自的相似元素集合中。
S13：采用枚举法枚举产生所有可能的系统状态，生成初始系统状态集合。
S2：应用最大风险分析模型求解当前系统状态集合中风险最高的状态。
S3：生成S2步骤求解结果的相似故障集合，并从系统状态集合中删去该相似故障集合，以生成新的系统状态集合。具体步骤为：
S31：以S1步骤得到的相似故障元素集为单位将当前故障状态分解为片段，并生成各状态片段的相似故障状态片段。对于任意一个相似故障元素集而言，只要状态片段内的故障元件个数相等，则均为相似故障状态片段。
S32：对各个相似元素集合的相似故障片段进行组合拼接，最终可生成当前系统状态的相似故障状态集合。其生成示意图如图3所示。
S33：在原始系统状态集合中删去由步骤S32得到的当前系统状态的相似故障状态集合，生成新的系统状态集合。
S4：应用风险分摊准则对已求得的系统故障事件集合进行风险分摊，更新各系统元件的风险贡献度。
以X1和X2两个元件故障时所对应的系统状态为例
当假设元件参数xi(i＝1,2…,n)相互独立，仅与自身相关，f(x1,x2,…,xn)可表示为式(30)。其中f(x1)仅与X1有关，f(x2)仅与X2有关，f(x3,…,xn)仅与除X1、X2以外的元件有关。
f(x1,x2,...,xn)＝f(x1)f(x2)f(x3，...,xn) (31)
本发明采用式(31)作为新的风险分摊依据，其中的λ_i、μ_i分别为元件的故障率和修复率，、分别表示当其他故障元件状态不变，而元件Xi处于故障和正常状态时系统切负荷量。式(32)的物理含义为在一次故障事件中，故障概率高和故障前后导致负荷损失变化较多的元件应承担较多的责任。
Figure BDA0001957784620000181
由所提出的可靠性跟踪准则可知，该系统状态风险应由故障元件X1和X2分担。其中X1和X2各自对应分摊的部分f(s→1)和f(s→2)可表示为:
Figure BDA0001957784620000182
Figure BDA0001957784620000183
S5：判断S3似故障集合的风险值是否高于所设定的系统风险阈值。若S3中相似故障集合的风险值高于所设定的系统风险阈值则转S2，否则转S6。
S6：输出各系统原件风险贡献度，确定系统薄弱环节。
本发明所提出的可靠性跟踪算法对RBTS系统在峰荷185MW下进行元件风险“贡献”分析，其结果如表1所示，实施例结果表明G1、G2、L9和G7四个元件对风险的贡献度占比达到了85％，因此可将其作为系统的薄弱环节。在运行维护、规划设计的时候，都应该着重考虑这四个元件的影响。表1同时也给出了基于状态解析法和蒙特卡诺模拟法进行可靠性跟踪得到的结果。其中状态解析法分析了4阶及4阶以下的系统故障状态；蒙特卡诺模拟次数为10000次；本文方法设定的风险阈值分别为10-5MW和10-6MW。
表1峰荷185MW下RBTS系统可靠性跟踪结果
Figure BDA0001957784620000191
应当指出，从本实施例的结果可以看出本发明所采用的方法计算的可靠性跟踪结果与状态解析法的精度基本一致，验证了本发明方法的有效性。通过本实施例的对比，虽然蒙特卡诺抽样数远大于解析法状态数，其对元件的风险“贡献”的衡量相比状态解析法仍存在较大的误差。相对于状态解析法，本发明方法在较高的精确度下两种风险阈值下的计算时间分别降低了82％和52％。实施例结果对比验证了本发明方法拥有更高的实用性。
以上对本实用新型实施例所提供的技术方案进行了详细介绍，本文中应用了具体个例对本实用新型实施例的原理以及实施方式进行了阐述，以上实施例的说明只适用于帮助理解本实用新型实施例的原理；同时，对于本领域的一般技术人员，依据本实用新型实施例，在具体实施方式以及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本实用新型的限制。

Claims (5)

1.一种基于最大风险分析的发输电系统可靠性跟踪算法，其特征在于，包括以下步骤：
S1、根据发输电系统原始数据生成相似故障元素集合，并生成初始系统状态集合，并设定系统风险阈值；
S2、应用最大风险分析模型求解初始系统状态集合中风险最高的状态；求解发输电系统的最大风险分析模型为直流潮流模型，其具体数学模型为
Figure FDA0003923456770000011
其中，中risk表示某一系统状态对应的系统风险值，将风险定义为切负荷量和对应状态概率的乘积，其中的pg表示发电机组的故障概率、pl表示线路的故障概率；xg、xl分别表示发电机组、线路的状态变量，
Figure FDA0003923456770000012
表示系统最优切负荷运行状态所对应的各节点切负荷量；Ng为发电机组总数、Nl为线路总数、Ng为削负荷节点总数；
发电机组、线路的状态变量的取值
xg,xl＝{0,1}
取0代表元件故障，取1代表元件正常工作；
各元件状态确定的系统状态应属于系统可选状态集S，即
[xg1,...,xgNg,xl1,...,xlNl]∈S
取值函数为
Figure FDA0003923456770000013
式中，Pg、
Figure FDA0003923456770000014
Fl、θb分别表示系统状态确定后，系统最优切负荷运行状态所对应的机组出力、切负荷量、线路功率、相角的取值；
每个母线节点应满足的功率平衡
Figure FDA0003923456770000021
其中
Figure FDA0003923456770000022
为节点b处的负荷，G(b)表示连接至节点b处的发电机组集合，而Lfr(b)、Lto(b)分别表示将节点b作为潮流发送、接受节点的线路集合；λb为该约束条件对应的对偶变量；
直流潮流下的线路潮流方程
Figure FDA0003923456770000023
采用了大M法对Fl＝xlbl(θi-θj)进行了线性化处理,其中bl为线路l的电纳，M为一个足够大的正数；
Figure FDA0003923456770000024
和
Figure FDA0003923456770000025
表示该约束条件对应的对偶变量，其中“+”表示对应每个母线节点应满足的功率平衡右边的不等式约束，而“-”表示对应每个母线节点应满足的功率平衡左边的不等式约束；
Figure FDA0003923456770000028
表示线路功率约束条件对应的对偶变量；
机组的出力范围约束
Figure FDA0003923456770000026
式中，
Figure FDA0003923456770000027
表示该约束条件对应的对偶变量；
母线上负荷切除量的取值约束
Figure FDA0003923456770000031
式中，
Figure FDA0003923456770000032
表示该约束条件对应的对偶变量；
母线的相角取值约束，并指定1节点处为参考节点
Figure FDA0003923456770000033
式中，
Figure FDA0003923456770000034
表示该约束条件对应的对偶变量；
将所述发输电系统的最大风险分析模型的双层模型，通过对偶理论转化为单层模型求解，其单层模型的数学式为
Figure FDA0003923456770000035
其中，原优化问题中下层问题所对应的强对偶定理条件为
Figure FDA0003923456770000036
对偶问题的约束条件为
Figure FDA0003923456770000037
Figure FDA0003923456770000038
Figure FDA0003923456770000041
Figure FDA0003923456770000042
Figure FDA0003923456770000043
Figure FDA0003923456770000044
Figure FDA0003923456770000045
Figure FDA0003923456770000046
S3、生成S2步骤求解结果的相似故障集合，并从初始系统状态集合中删去该相似故障集合，以生成新的系统状态集合；对所述发输电系统的最大风险分析模型的求解，还包括实现单层模型的线性化，其过程包括，
第一步，对目标函数，即单层模型的数学式进行线性化处理，首先进行对数变换，
Figure FDA0003923456770000047
并采用分段函数对式中的
Figure FDA0003923456770000048
进行线性化处理，
y＝min(yi)
Figure FDA0003923456770000051
Figure FDA0003923456770000052
第二步，原优化问题中下层问题所对应的强对偶定理条件中的双变量乘积进行线性化处理，首先引入辅助变量
Figure FDA0003923456770000053
同时采用线性约束实现
Figure FDA0003923456770000054
的等效线性化
Figure FDA0003923456770000055
Figure FDA0003923456770000056
Figure FDA0003923456770000057
Figure FDA0003923456770000058
其中Y、Z均为一个足够大的正数；
对线性化处理后的模型进行求解可即可确定当前系统状态集S内的风险最高的系统故障状态；
S4、步骤S3中，包括根据步骤S1生成的相似元件集合生成步骤S2中求解得到的当前最大风险系统状态的相似故障集合，并生成新的系统状态集合，具体步骤为，
S31、以S1步骤得到的相似故障元素集为单位将当前故障状态分解为片段，并生成各状态片段的相似故障状态片段；
S32、对各个相似元素集合的相似故障片段进行组合拼接，最终可生成当前系统状态的相似故障状态集合；其中相似故障状态集合内的状态个数N为
Figure FDA0003923456770000061
式中的ns表示相似元素集合个数，Ni为第i个相似故障元素集合的相似片段个数；
S33、在原始系统状态集合中删去由步骤S32得到的当前系统状态的相似故障状态集合，生成新的系统状态集合；应用风险分摊准则对已求得的系统故障事件集合进行风险分摊，更新各系统元件的风险贡献度；
S5、判断步骤S3中相似故障集合的风险值是否高于所设定的系统风险阈值，若是则返回步骤S2；如否进入下一步骤；
S6、输出各系统元件风险贡献度，确定系统薄弱环节。
2.根据权利要求1所述的一种基于最大风险分析的发输电系统可靠性跟踪算法，其特征在于：所述步骤S1中所述的发输电系统原始数据包括发输电系统的可靠性参数和电气参数；其根据发输电系统原始数据生成相似故障元素集合，并生成初始系统状态集合的具体步骤为，
S11、根据系统设备种类将系统元件分类；
S12、在同一种类的系统元件内，分别确定各自的相似指标；
S13、采用枚举法枚举产生所有可能的系统状态，由此生成初始系统状态集合。
3.根据权利要求1所述的一种基于最大风险分析的发输电系统可靠性跟踪算法，其特征在于：步骤S4中包括应用风险分摊准则对已求得的系统故障事件集合进行风险分摊；其中风险分摊准则包括如果某两个系统状态A与B，A状态中比B状态新增了m个故障元件，除这m个故障元件外其余元件的状态均相同；而这两个系统状态的负荷损失相同，即这m个元件的故障没有令系统状况恶化，则这m个元件不参与A状态风险值的分担；
基于以上不可靠性跟踪准则，系统元件对发输电系统可靠性的贡献评估如下：
假设该发输电系统由n个元件Xi(i＝1,2…,n)组成，同时采用xi(i＝1,2…,n) 表示各元件所对应的性能参数，f(x1,x2,…,xn)表示由系统状态为(X1,X2,…,Xn)时的系统风险值risk；
以X1和X2两个元件故障时所对应的系统状态为例，当假设元件参数xi(i＝1,2…,n)相互独立，仅与自身相关，f(x1,x2,…,xn)可表示为
f(x1,x2,…,xn)＝f(x1)f(x2)f(x3，…,xn)
其中f(x1)仅与X1有关，f(x2)仅与X2有关，f(x3,…,xn)仅与除X1、X2以外的元件有关；
风险分摊依据
Figure FDA0003923456770000071
其中的λi、μi分别为元件的故障率和修复率，ΔPd、ΔPi d分别表示当其他故障元件状态不变，而元件Xi处于故障和正常状态时系统切负荷量，
其中X1和X2各自对应风险分摊的部分f(s→1)和f(s→2)可表示为:
Figure FDA0003923456770000072
Figure FDA0003923456770000073
4.根据权利要求3所述的一种基于最大风险分析的发输电系统可靠性跟踪算法，其特征在于：步骤S4中按照比例分摊准则确定各元件风险贡献比之前，应该首先判断各元件的风险分摊系数；通过搜索比当前状态风险高的故障状态集，即在截止当前已经确定的解集进行判断。
5.根据权利要求1所述的一种基于最大风险分析的发输电系统可靠性跟踪算法，其特征在于：S6骤中的各系统元件风险贡献度由每次循环中步骤S4得到的风险贡献度累加得到，并根据元件风险贡献度高低选取风险贡献前十的元件作为系统薄弱环节集合。
CN201910072739.6A 2019-01-25 2019-01-25 一种基于最大风险分析的发输电系统可靠性跟踪算法 Active CN109886834B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910072739.6A CN109886834B (zh) 2019-01-25 2019-01-25 一种基于最大风险分析的发输电系统可靠性跟踪算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910072739.6A CN109886834B (zh) 2019-01-25 2019-01-25 一种基于最大风险分析的发输电系统可靠性跟踪算法

Publications (2)

Publication Number Publication Date
CN109886834A CN109886834A (zh) 2019-06-14
CN109886834B true CN109886834B (zh) 2023-01-10

Family

ID=66926907

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910072739.6A Active CN109886834B (zh) 2019-01-25 2019-01-25 一种基于最大风险分析的发输电系统可靠性跟踪算法

Country Status (1)

Country Link
CN (1) CN109886834B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110969355B (zh) * 2019-12-03 2023-06-09 重庆大学 增量风险事件的筛选方法、装置及计算机可读介质
CN111651889B (zh) * 2020-06-04 2024-04-26 重庆大学 输电系统高风险事件筛选方法、装置、设备及存储介质
CN111881588B (zh) * 2020-07-30 2023-10-03 重庆大学 基于最小点火能的爆炸风险量化方法
CN116565872A (zh) * 2023-02-06 2023-08-08 天津大学 基于状态相似性的交流电网可靠性快速控制方法与装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102737286A (zh) * 2012-04-23 2012-10-17 中国电力科学研究院 一种地区电网在线风险分析系统和方法
CN103426056A (zh) * 2013-07-18 2013-12-04 清华大学 基于风险评估的电力系统薄弱环节辨识方法
CN103985065A (zh) * 2014-05-20 2014-08-13 天津大学 一种基于故障预扫描的电力系统风险评估方法
CN104331788A (zh) * 2014-12-03 2015-02-04 太原理工大学 发电系统可靠性评估的严重故障集快速筛选方法
CN106503923A (zh) * 2016-11-15 2017-03-15 国网四川省电力公司电力科学研究院 一种电力系统连锁故障风险评估的方法及系统

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8417360B2 (en) * 2001-08-10 2013-04-09 Rockwell Automation Technologies, Inc. System and method for dynamic multi-objective optimization of machine selection, integration and utilization
CN102708411A (zh) * 2012-04-23 2012-10-03 中国电力科学研究院 一种地区电网在线风险评估方法
CN103793854B (zh) * 2014-01-21 2015-09-30 中国南方电网有限责任公司超高压输电公司检修试验中心 多重组合优化的架空输电线路运行风险信息化评估方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102737286A (zh) * 2012-04-23 2012-10-17 中国电力科学研究院 一种地区电网在线风险分析系统和方法
CN103426056A (zh) * 2013-07-18 2013-12-04 清华大学 基于风险评估的电力系统薄弱环节辨识方法
CN103985065A (zh) * 2014-05-20 2014-08-13 天津大学 一种基于故障预扫描的电力系统风险评估方法
CN104331788A (zh) * 2014-12-03 2015-02-04 太原理工大学 发电系统可靠性评估的严重故障集快速筛选方法
CN106503923A (zh) * 2016-11-15 2017-03-15 国网四川省电力公司电力科学研究院 一种电力系统连锁故障风险评估的方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于模糊理论的地区电网在线供电风险分析方法研究;罗雅迪;《电工电能新技术》;20141231;第33卷(第12期);1-6 *

Also Published As

Publication number Publication date
CN109886834A (zh) 2019-06-14

Similar Documents

Publication Publication Date Title
CN109886834B (zh) 一种基于最大风险分析的发输电系统可靠性跟踪算法
CN110969347B (zh) 一种输电网结构形态评估方法
CN105976257A (zh) 基于隶属度函数的模糊综合评价法的电网脆弱性评估方法
CN103296685A (zh) 一种svc补偿策略最优化方法
CN113705085B (zh) 一种智能电网多层次结构建模及风险评估方法
Kamalasadan et al. A new intelligent algorithm for online voltage stability assessment and monitoring
Contreras-Jiménez et al. Multi-state system reliability analysis of HVDC transmission systems using matrix-based system reliability method
CN108985602B (zh) 一种计及风险的电网分类项目的投入产出评价方法及系统
Shahmohammadi et al. Proper sizing and placement of distributed power generation aids the intentional islanding process
CN105160459A (zh) 一种电力系统的稳定运行状况评估方法
Souza et al. Immune system memetic algorithm for power distribution network design with load evolution uncertainty
US9261864B2 (en) Method and device for producing a state signal
Venkatesh et al. Machine Learning for Hybrid Line Stability Ranking Index in Polynomial Load Modeling under Contingency Conditions.
Haidar et al. Artificial Intelligence application to Malaysian electrical powersystem
Gholami et al. Cyber-physical power system reliability assessment considering multi-state independent components
CN106875026B (zh) 一种电力市场环境下中长期输电网拓展规划的混合性规划方法
Tudose et al. Power systems security assessment based on artificial neural networks
CN113177717A (zh) 一种基于影响增量灵敏度的输电系统韧性快速评估方法
Acharya Weather effect considerations in reliability evaluation of electrical transmission and distribution systems
CN108470256A (zh) 跨省跨区加强输电通道电网工程运行效益评价方法及系统
Gharebaghi et al. Solving the divergence problem in AC-QSS cascading failure model by introducing the effect of a realistic UVLS scheme
CN110601189B (zh) 一种考虑连锁跳闸的电网安全性评估方法
CN116231646B (zh) 基于电力系统脆弱和经济性的pmu最优配置方法及系统
US11983472B2 (en) Method for identifying fragile lines in power grids based on electrical betweenness
CN118643750B (zh) 牵引供电系统负荷场景削减方法

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