基于邻接矩阵法的变电站系统抗震可靠性分析
Seismic reliability analysis of substation system based on adjacency matrix
收稿日期: 2021-08-15
基金资助: |
|
Received: 2021-08-15
Fund supported: | 国家重点研发计划资助项目(2019YFC1509302);陕西省重点研发计划资助项目(2021ZDLSF06-10);西安市科技计划资助项目(2019113813CXSF016SF026) |
作者简介 About authors
刘晓航(1995—),男,博士生,从事生命线工程抗震研究.orcid.org/0000-0001-7988-5507.E-mail:
针对变电站系统后评估过程中的计算效率问题,提出基于邻接矩阵的可靠性评估流程用以分析变电站系统. 邻接矩阵所建立的边权模型能够直观地反映系统单元和设备间的逻辑关系,拟Warshall算法通过对邻接矩阵元素的布尔运算高效求解连通性矩阵,从而以Quasi-Monte Carlo模拟方法计算出整个系统的功能状态. 基于此评估流程研究了一个典型的6进线10出线220/110/10 kV变电站,计算其抗震可靠性并确定了抗震关键设备. 案例分析结果表明,基于邻接矩阵的系统评估思路是可行的,220 kV的电压互感器、断路器、隔离开关及110 kV的隔离开关被评定为变电站系统中最关键的功能设备,提高它们的抗震性能可以显著提高整个变电站系统的抗震可靠性.
关键词:
A reliability evaluation process based on adjacency matrix was proposed to analyze substation system, aiming at the problem of computational efficiency in the post evaluation of substation system. The edge weight model established by the adjacency matrix can intuitively reflect the logical relationship between the system units and the equipment. The quasi-Warshall algorithm efficiently solves the connectivity matrix by Boolean operations on the adjacency matrix elements, thereby calculating the functional status of the entire system by the Quasi-Monte Carlo simulation method. A typical 220/110/10 kV substation with 6 incoming lines and 10 outgoing lines was studied based on this evaluation process. Its seismic reliability was calculated and the key seismic equipment was determined. The case analysis results showed that the system evaluation idea based on adjacency matrix was feasible. 220 kV voltage transformers, circuit breakers, isolating switches and 110 kV isolating switches were assessed as the most critical functional equipment in the substation system. Improving their seismic performance can significantly improve the seismic reliability of the entire substation system.
Keywords:
本文引用格式
刘晓航, 郑山锁, 黄瑜, 董淑卿, 杨丰, 董晋琦.
LIU Xiao-hang, ZHENG Shan-suo, HUANG Yu, DONG Shu-qing, YANG Feng, DONG Jin-qi.
在变电站系统抗震可靠性分析方面,Hwang等[3]通过二分试验法将系统的每一种状态以状态向量表示,系统总状态数量为2n,计算量将随着单元数量的增加而急剧增加. Volkanovski等[4]通过故障树建立了变电站系统的基本模型,所建模型能够清楚、完整地表达系统单元之间的逻辑关系,并通过最小割法(minimum cut set,MCS)对建立的故障树进行化简[5]. 考虑到变电站中存在的冗余性,李吉超[6]提出基于概率的状态树方法,用以评估复杂系统的地震可靠性,状态树的建立过程实际上是在先评估过程中寻找所有的成功路径[7]. 综上所述,目前对于变电站系统抗震可靠性的研究多聚焦于先评估过程,即考虑如何表达设备间的串并联关系,建立合理的系统模型,而忽略后评估过程中的求解效率. 邻接矩阵法作为后评估的一种典型方法,能够有效表示系统节点间的拓扑关系,实现程序化计算. 近年来已在车辆工程[8]、电力系统、供水系统和道桥系统[9-11]等诸多领域有所应用,然而其在变电站系统领域的相关研究仍相对欠缺.
本研究提出新的基于邻接矩阵法的变电站系统抗震可靠性分析流程,邻接矩阵法在先评估过程中建立的系统模型能够充分考虑各单元间的相互作用,而所采用的拟Warshall算法则能够大幅减少后评估过程中的计算量[12]. 所建立的系统模型通过Quasi-Monte Carlo (QMC)模拟能够准确地计算出整个系统的失效概率. 本方法的优点是在建模过程中能够保留系统单元间的逻辑关系,结合“功能单元”概念与边权网络模型,综合考虑各供电设施的易损性,大大降低了系统网络基本模型的复杂度,使其在拟蒙特卡罗仿真中具有可计算性. 在评估过程中可以实现程序化模拟并高效得到系统的抗震可靠性分析结果.
1. 变电站系统抗震可靠性分析流程
1.1. 邻接矩阵法建模
变电站系统网络具有大量的供电设施,并且这些设施具有复杂的串并联关系. 为了研究变电站系统的抗震可靠性,首先将变电站系统抽象为一个由许多节点组成的几何图形,并通过相关的链路将这些节点连接起来. 节点代表这些供电设施,例如设备和架构,边代表忽略了其物理结构和电气参数差异的母线. 变电站系统被抽象为没有权重的有向拓扑图,以集合的形式表示如下:
式中:非空集合
为了实现程序化计算,采用拓扑图的邻接矩阵来描述其结构. 若有向拓扑图
式中:若
1.2. QMC模拟方法
QMC模拟方法与常规Monte Carlo (MC)模拟方法相似,但是其理论基础不同. QMC模拟的基本思路是采用具有较低离散型的低偏差序列(low discrepancy sequence, LDS)替代MC模拟方法所对应的伪随机数序列[13-15]. 在先前电力系统抗震可靠性研究中已经证明QMC模拟方法的高效性[16],因此其理论基础不再详述. 对于QMC模拟方法,LDS的构造是其核心内容. LDS的特性为构造过程简便,且序列均匀性较好,通常可以采用Sobol或Halton低偏差序列. 本研究采用Sobol序列作为QMC模拟的原始抽样序列,并基于线性矩阵加扰技术将Sobol低偏差序列进行随机化[17],保证其具有低偏差与随机性的双重特征. 随机过程分以下2步.
1)以一个加扰动的
式中:t为Sobol序列的维度变量,
2)对每个新的向量增加一个二进制相互独立的随机向量
式中:
1.3. 拟Warshall算法求解连通性矩阵
引入连通性矩阵的概念,连通性矩阵M指用矩阵形式来描述有向图中各节点之间经过一定数量边相连的状况,
1)对于
2)置新矩阵
3)定义所需输出的矩阵元素坐标,置其所在行号为s;
4)遍历s行的每一列j,若判定
5)重复检索,若存在未被选择的j,使得
2. 典型变电站案例研究
本案例研究选择了一个在中国分布广泛的典型6进线10出线220/110/10 kV变电站[18]. 其中,变电站的10 kV部分不易受到地震影响,因此假定该部分在地震后仍可正常工作. 220 kV部分中有6条进线,110 kV部分中有10条出线. 每条线路独立工作并连接至管母线,管母线用于收集和分配电能. 变电站系统各设备布局如图1所示. 图中,SF为变电站架构,用于支撑进出线路,DS、CT、CB、IS、PT、LA、BUS、TR分别对应为隔离开关、电流互感器、断路器、绝缘子、电压互感器、避雷器、管母线和变压器. 本研究的目的是验证所提出的可靠性分析流程的有效性,供电设施的初始易损性参数均来自震害统计,对震害统计处理后绘制的易损性曲线如图2所示. 图中,PGA为地震峰动值加速度,P为超越概率.
图 1
图 2
图 2 基于震害经验统计的供电设施抗震易损性曲线
Fig.2 Seismic fragility curves of power supply facilities based on empirical statistics of earthquake damage
2.1. 变电站功能定义
变电站功能Fs表示110 kV部分中10条出线的可靠性. 为了与功能单元的易损性相一致,文献[19]定义P(Fs)为地震中不超过N条出线可靠的概率. Fs的值为0~9/10,增量为1/10. 例如Fs=1/2表示10条线路中正常运转的线路为5条,Fs=9/10表示10条线路中正常运转的线路为9条.
变电站系统存在冗余性,其在正常运行条件下不会满负荷运行(通常低于75%). 在这种情况下,完全运行的2台变压器足以维持供电. 同时假定1条进线的负荷可以支撑2条出线. 因此,变电站的功能可以表达为
式中:
2.2. 变电站系统功能单元划分
为了展示各功能单元构成,将变电站划分为4部分:220 kV部分、母联部分、110 kV部分和变压器部分,如图3~6所示. 在立面图中红线为化简后的输电导线,实际上分为A、B、C三相传输. 这些设备彼此独立,但架构由不同线路的设备共享. 此外,母联部分用于组合2条不同的管母线,如图4所示给出了220 kV母联部分的示意图,110 kV母联与220 kV母联相似,只是设备间距不同. 考虑到该变电站系统中的大量设备和结构,共计664个组件. 首先采用功能单元概念进行化简,将各类结构设备依据不同功能划分为进线单元、变压器单元、出线单元、母联单元、母线单元、连接单元. 等价功能单元的破坏概率为各组成设施的联合破坏概率:
图 3
图 4
图 5
图 6
式中:
图 7
图 8
图 8 变电站系统的连接单元与母联单元
Fig.8 Connection unit and bus coupler unit of substation system
图 9
变压器单元(见图10)由3个SF、2个CB、2个CT、2个DS、2个LA和TR本身组成. 最易破坏的构件是安装在TR上部的瓷套管. 变压器连接到220 kV线路、110 kV线路和10 kV线路,因此存在3种类型的套管:高压套管(HB)、中压套管(MB)和低压套管(LB). TB为变压器的油箱.
图 10
2.3. 变电站系统的基本网络模型
逻辑关系图可以考虑组件之间的相互关系,并给出变电站系统所有可能序列的完整描述. 将变电站系统按功能单元作出逻辑关系图,如图11所示. 图中的系统包含了69个节点与若干连接线路. 为了研究变电站系统地震作用下的工作性能,必须先建立其网络基本模型.
图 11
图 11 基于功能单元的变电站系统逻辑关系图
Fig.11 Logic diagram of substation system based on functional unit
图 12
图 12 基于边权模型化简的变电站系统逻辑关系图
Fig.12 Logic diagram of substation system simplified based on edge weight model
2.4. 基于邻接矩阵法的变电站系统抗震可靠性分析
采用邻接矩阵法建立系统网络模型,以QMC法随机模拟变电站系统各功能单元在每一次模拟中的地震破坏状态,采用改进的拟Warshall算法对变电站系统进行抗震可靠性分析,基本步骤如下.
1)确定EDP范围e(本研究以PGA作为工程需求参数)、功能区间f及总模拟次数N.
2)求解在第e次循环下各功能单元对应的破坏概率Pc,ij,并生成初始邻接矩阵A.
3)以QMC法对变电站系统中所有功能单元产生低偏差数列rij,并与功能单元的破坏概率Pc,ij相比较,修正本次模拟的邻接矩阵A,其中,
4)定义进线单元、变压器单元及出线单元所对应的矩阵元素坐标,并以拟Warshall算法求解得到变电站系统的功能性矩阵M. 对所有的j列,若msj=1,则源点(控制室)s与功能单元j之间连通,若msj=0,则源点s与功能单元j之间不连通.
5)统计M矩阵中进线单元、变压器单元及出线单元对应元素msj=1的数量,并以式(7)判别统计,将统计结果累计到50×10阶破坏状态矩阵T中,矩阵T中的元素tef依次代表变电站系统在不同PGA下不同功能状态出现的次数. 在本案例研究中,PGA范围为0~0.5g,间隔为0.01g,因此e=50;变电站共划分为10个功能状态,取f=10. 矩阵T中元素表达式为
6)每个e值都重复步骤1)~5)N次,统计变电站系统各功能状态出现频率. 为了求得精确值,本研究取N=20000. 频率表达式为
7)将各功能状态出现频率逐级累加得到变电站各功能对应的超越概率P(Fs).
基于上述相关理论,变电站系统可靠性分析流程如图13所示.
图 13
图 13 基于邻接矩阵法的变电站系统仿真过程
Fig.13 Simulation process of substation system based on adjacency matrix method
3. 变电站系统抗震可靠性分析及关键设备识别
3.1. 抗震可靠性分析结果
变电站系统不同功能的易损性曲线如图14所示. 可以看出,当Fs<8/10时,相邻Fs的结果高度相似(例如,当Fs=6/10、7/10时,抗震能力中值分别为0.183 232、0.182 501,当Fs=0、1/10时,抗震能力中值皆为0.295 230). Fs=5/10和Fs=4/10的易损性曲线有较大差别,体现了变电站系统的关联性和冗余性. 理论上讲,1个进线单元可以支撑2个出线单元,因此当1个进线单元破坏应导致变电站系统功能下降2/10. 但是,进线单元中存在20%的冗余,因此在第2个进线单元发生破坏时系统功能才会受到影响. 同样,变压器单元中存在50%的冗余,除非3个变压器单元中的2个破坏,否则系统功能不会受到影响.
图 14
图 14 变电站系统不同功能的易损性
Fig.14 Vulnerability of different functions of substation system
如表1所示为变电站系统不同功能的中值能力mc和标准偏差sd. 可以看出,当PGA中值低于0.15g时,大多数的功能单元依然可以正常运行,变电站系统的功能处于较高值(不小于8/10). 此时,变电站系统是冗余的,其功能取决于出线单元,任何一个出线单元的破坏将使变电站系统的功能下降1/10. 随着PGA的增加,进线单元的冗余度降为0,一个进线单元的破坏将导致变电站功能下降2/10.
表 1 变电站系统不同功能的中值能力和标准偏差
Tab.1
Fs | mc | sd | R2 | RMSE |
0 | 0.295 230 | 0.206 9 | 0.999 4 | 0.010 61 |
1/10 | 0.295 230 | 0.206 7 | 0.999 5 | 0.010 71 |
2/10 | 0.269 012 | 0.199 2 | 0.999 5 | 0.010 09 |
3/10 | 0.269 012 | 0.199 1 | 0.999 5 | 0.010 20 |
4/10 | 0.235 746 | 0.216 9 | 0.999 2 | 0.012 38 |
5/10 | 0.208 670 | 0.224 3 | 0.999 3 | 0.011 49 |
6/10 | 0.183 232 | 0.249 3 | 0.998 9 | 0.014 53 |
7/10 | 0.182 501 | 0.244 9 | 0.998 8 | 0.015 24 |
8/10 | 0.144 136 | 0.326 9 | 0.996 1 | 0.024 97 |
9/10 | 0.135 742 | 0.321 6 | 0.996 8 | 0.022 41 |
为了更清楚地说明这一点,如图15所示为变电站系统随PGA变化的不同功能的频率f. 当PGA<0.15g时,变电站仍保持完全运行状态,而当PGA>0.23g时,变电站丧失大半功能. 可以观察到,Fs=6/10并没有随着Fs=8/10,Fs=4/10,Fs=2/10呈现逐级下降规律,Fs=6/10下的峰值概率远小于Fs=4/10,表明变压器单元的易损性过大导致其冗余性并不能有效保证变电站的抗震可靠性. 因此,这组曲线清楚地表达了组件冗余对系统功能的影响. 结果表明,系统设计中考虑的冗余度不足以抵御0.23g的最大可信地震. 同时,当PGA较小时,变电站功能取决于出线单元及110 kV连接单元,当PGA较大时,变电站的功能主要取决于进线单元、220 kV连接单元与主变压器单元. 这是由于当PGA较大时,即便进线单元与变压器单元存在一定的冗余度,但进线单元与变压器单元的易损性远大于出线单元和110 kV连接单元.
图 15
图 15 变电站系统不同功能的出现频率
Fig.15 Frequency of appearance of different functions of substation system
3.2. 关键设备识别
图 16
图 16 变电站系统抗震可靠性曲线中值能力变化(110 kV设备和支柱绝缘子)
Fig.16 Changes in median capacity of seismic vulnerability curve of substation system (110 kV equipment and pillar insulators)
图 17
图 17 变电站系统抗震可靠性曲线对数标准差变化(110 kV设备和支柱绝缘子)
Fig.17 Logarithmic standard deviation change of seismic vulnerability curve of substation system (110 kV equipment and pillar insulators)
图 18
图 18 变电站系统抗震可靠性曲线中值能力变化(220 kV设备)
Fig.18 Change in median capacity of seismic vulnerability curve of substation system (220 kV equipment)
图 19
图 19 变电站系统抗震可靠性曲线对数标准差变化(220 kV设备)
Fig.19 Changes in logarithmic standard deviation of seismic vulnerability curve of substation system (220 kV equipment)
由图16~17可知,110 kV设备和IS对变电站系统易损性曲线的中值能力和对数标准差影响不大(总体小于1.5%),因为这些设备具有相对较高的能力中值. 可以发现,除了110 kV的LA、IS,110 kV的CB、DS、CT、PT在Fs=9/10时对变电站的可靠性呈现整体上升,这是因为进线单元、220 kV连接单元和变压器单元在PGA较小时才是冗余的,而110 kV连接单元与出线单元包含1个CB、1个PT、1个CT和3个DS,因此110 kV的DS的变化率最大. 提高110 kV的DS的抗震性能可以显著提高变电站系统Fs=9/10的抗震可靠性. 此时,变电站功能状态由110 kV连接单元与出线单元决定,这与3.1节中的结果可以相互验证.
由图18可知,220 kV设备中CT对变电站系统易损性曲线的中值能力影响最大,随着Fs的减小在4%~7%变化. 随后是CB,中值能力随Fs的减小在3%~6%变化. 其余设备对中值能力的影响由大到小分别为PT、DS、TR、LA. 当Fs=9/10时220 kV设备的变化率显著减小,也验证了3.1节中的分析,即在系统抗震可靠性较高时,出线单元决定了系统的震后功能.
由图19可知,CB、CT、DS、PT对系统易损性曲线的对数标准差也有较大影响,但影响方式不同,即CB、CT减小对数标准差,并与Fs正相关,DS、PT增大对数标准差,与Fs成负相关. LA与PT、TR由于其较高的能力中值及较少的参与度,对变电站系统能力中值(<1.5%)的影响非常有限. 表明220 kV的LA与PT、TR单体设备具有较高的抗震冗余性. 因此,在进线单元、220 kV连接单元及变压器单元中,220 kV的CT、CB、DS被认为是变电站中最关键的功能设备,提高他们的抗震性能可以显著提高整个变电站系统的抗震可靠性.
4. 结 论
(1)提出一种变电站系统抗震可靠性评估流程. 该流程利用邻接矩阵对系统拓扑和功能进行分析,引入“功能单元”概念,综合考虑各供电设施的易损性,大大降低了系统网络基本模型的复杂度,使其在拟蒙特卡罗仿真中具有可计算性. 针对变电站系统功能性计算,基于Warshall算法提出适用于各类变电站系统可靠性分析的功能性算法,并应用此算法对典型220 kV变电站系统进行抗震可靠性分析.
(2)通过关键设备分析确定220 kV的CT、CB、DS及110 kV的DS对系统功能影响最大,是变电站中最关键的设备. 本研究提出的方法为程序化评估变电站系统抗震可靠性提供了相应的思路,并为变电站系统以及电力系统抗震韧性评估奠定了研究基础.
(3)可靠性评估结果的准确性取决于输入的设施易损性参数,为了更精准地获得评估结果,后续计划开展更为详细的变电设施地震易损性研究.
参考文献
Modelling catastrophic degradation of flexural-dominated RC columns at ultimate displacements based on fibre beam-column model
[J].DOI:10.1016/j.jobe.2021.103476 [本文引用: 1]
Prediction of residual behaviour for post-earthquake damaged reinforced concrete column based on damage distribution model
[J].DOI:10.1016/j.engstruct.2021.111927 [本文引用: 1]
Evaluation of seismic performance of an electric substation using event tree/fault tree technique
[J].DOI:10.1016/S0266-8920(97)00018-0 [本文引用: 1]
Application of the fault tree analysis for assessment of power system reliability
[J].DOI:10.1016/j.ress.2009.01.004 [本文引用: 1]
Evaluation of municipal water distribution system reliability using minimum cut-set method
[J].
Creating fragility functions for performance-based earthquake engineering
[J].DOI:10.1193/1.2720892 [本文引用: 1]
换挡工况下湿式换挡离合器变胞机理
[J].
Metamorphic mechanism of wet shift clutch in gear shifting process
[J].
Research on a seismic connectivity reliability model of power systems based on the quasi-Monte Carlo method
[J].
供水系统震后水力分析算法并行化研究
[J].
Research on parallel algorithm for post-earthquake hydraulic analysis of water supply system
[J].
基于图卷积神经网络的城市交通态势预测算法
[J].
Urban traffic flow prediction algorithm based on graph convolutional neural networks
[J].
变电站主接线系统地震易损性分析
[J].
Seismic vulnerability analysis of substation system
[J].
An applied complex systems framework for risk-based decision-making in infrastructure engineering
[J].DOI:10.1016/j.strusafe.2014.03.011 [本文引用: 1]
Quasi-Monte Carlo based uncertainty analysis: sampling efficiency and error estimation in engineering applications
[J].DOI:10.1016/j.ress.2019.106541
Different numerical estimators for main effect global sensitivity indices
[J].DOI:10.1016/j.ress.2017.04.003 [本文引用: 1]
基于拟蒙特卡罗法的电力系统抗震可靠性研究
[J].
Research on seismic reliability of power systems based on quasi-Monte Carlo method
[J].
Comparing halton and sobol sequences in integral evaluation
[J].
四川电网汶川地震电力设施受灾调研分析
[J].DOI:10.13335/j.1000-3673.pst.2008.11.021 [本文引用: 2]
Investigation and analysis of electric equipment damage in Sichuan power grid caused by wenchuang earthquake
[J].DOI:10.13335/j.1000-3673.pst.2008.11.021 [本文引用: 2]
/
〈 |
|
〉 |
