摘要：在对燃气管道进行地震反应判别的基础上，计算管道的抗震可靠度，并结合图论将城市燃气管网简化为源点、汇点和边线，建立管网简化模型。采用蒙特卡罗法对某六节点网络模型进行抗震连通可靠性分析，结果表明：随着迭代次数的增加，近似解以更小的误差接近精确解的概率更大；迭代次数取为5000时，能够较好地满足精度要求。通过对连通性矩阵进行分析，探讨了设定源点提高目标汇点连通可靠度的可行性。最后，以保定市区的大型燃气管网为例，计算了Ⅶ～Ⅸ度地震烈度下管网的连通可靠度，对比了设定源点对汇点连通可靠度的影响，为该市燃气管网的震害预测和抗震优化提供参考。
关键词：蒙特卡罗法；城市燃气管网；设定源点；汇点；抗震连通可靠性度
中图分类号：TV352 文献标识码：A 文章编号：1000-0666（2015）02-0292-05
0 引言
城市燃气管网是城市生命线工程的重要组成部分，也是地震次生灾害源之一（Michael，Ron-ald， 1996；周伟国等，2009）。燃气管网的抗震研究主要包括管道单体震害机理研究和网络可靠性分析两个方面（Li，He，2002），其中，网络可靠性分析是研究热点。
目前，网络可靠性分析方法主要有解析算法和蒙特卡罗法（刘威，李杰，2002）。解析算法可以获得网络连通可靠度的精确解，但对于大型网络（节点超过100个）的求解，会因非多项式增长问题而导致计算网难。蒙特卡罗法是一种随机模拟方法，能够给出解的近似值（徐钟济，1985），在网络连通性的计算分析中不考虑拓扑结构的复杂性，因而计算效率较高，随着迭代次数的增加，精度在一定范围内也会有不断提高的趋势，但难以预测精度是否满足要求。
1 燃气管道抗震可靠度计算
1.1 管道地震反应判别
城市燃气管网主要由埋地管道构成，其震害程度与地震动强度、场地土特质、管道特征等诸多因素有关（汤爱平，欧进萍，2002）。埋地管道在地震波作用下的可靠性分析主要有两种判别指标：对于接口式管道，以接口轴向形变作为判别指标；对于连续焊接管道，以应力作为判别指标。管道地震反应计算方法见郭恩栋等（1996）文献。
1.2 管道抗震可靠度预测模型
在地震作用下管道的状态功能函数设为
2 燃气管网抗震连通可靠性分析
2.1 基于图论的管网简化模型
图包括点和连接各点的矢量线，它是顶点集合和顶点的有序偶集合（E・米涅卡，1984），图论正是研究图形而发展起来的理论。燃气管网通过管线和节点的相瓦连接构建成网络，将燃气储配站作为源点，管线交汇点、用户接人点作为汇点，管线作为连接各点的边，便可简化为图。通过分析计算不同地震烈度影响下的管段损坏的概率，对边赋予相应的权值，可以对源点与汇点的连通可靠度进行分析，进而得出网络的连通可靠度。
可用邻接矩阵A的形式表示图，便于进行数值分析。
2.2 基于蒙特卡罗法的网络连通性分析
（1）用线性同余法产生0到1之间最大周期的伪随机数rij（吴新瞻，吴新垣，1990），与计算所得的管线抗震可靠度pij比较，得到所有节点的邻接矩阵A，矩阵中元素aij的取值规则为
近精确解的概率更大。
采用迭代5000次的蒙特卡罗法对图1中每个汇点各重复计算6000次，获得6000组连通可靠度数据，用显著水平为0.000063的置信区间表示近似解，并与全概率分解法（梅启智等，1992）计算的精确解对比，可以看出蒙特卡罗方法在大量数据支撑下，能够无限逼近精确解（表1）。由于实际应用中一般只计算一组连通可靠度数据，即单次模拟，再用这6000组数据的均值与标准差来表示单次模拟所求近似解的精度。
从表1和图2可以看出，利用蒙特卡罗法迭代5000次，能够将近似解的精度控制在一定的范围，根据管线的抗震可靠度计算精度的实际需求，可以直接应用单次模拟所求近似解。
网络连通性分析中，在不改变拓扑结构的情况下，汇点的连通可靠度与源点的相对位置有一定的关系，这在连通性矩阵Pij中有所表现，第i行中的每一个元素代表汇点（j≠1）与设定源点的连通性。图1的连通性矩阵中，元素与设定源点、汇点的对应关系见表2。通过连通性矩阵Pij可以看出，选取不同的节点作为源点，对诸多节点的连通可靠度有一定程度影响。
3 工程实例
3.1 工程概况
保定市城区燃气管网总长约229 km，设计压力分级为中压A （0.4 MPa）、中压B（0.03～0.07 MPa）和低压（1.1～2.2 kPa）。由于中压A管网以燃气储配站为供气源，独立成系统，可单独进行分析。中压A管网于2002年初投入使用，主要采用无缝钢管和Q235B直焊缝钢管，公称商径为250～500 mm。市区场地土以Ⅲ类场地为主，部分地区为Ⅱ类场地。
3.2 管网抗震连通性分析
现以燃气储配站为源点，将中压A管网简化为113个节点、141条边线的火型网络图（图3）。以应力为判别指标，计算各管线在Ⅶ度、Ⅷ度和Ⅸ度地震烈度下的抗震可靠度，并赋权给相应的边线，再用蒙特卡罗法迭代5000次，求得不同汇点与源点的连通可靠度。
在不同地震烈度下，选取不同节点作为设定源点会对汇点抗震连通可靠度造成不同程度的影响。表3为在Ⅶ度、Ⅷ度和Ⅸ度地震烈度情况下，选取39号节点为设定源点时，与1号节点（燃气储配站）为源点时相比，部分汇点抗震连通可靠度变化幅度。从表中可以看出，在不同地震烈度下，设定源点对某些汇点连通可靠度的影响存在一定差别。
通过分析连通性矩阵可知，设定源点以改善局部重要节点或人部分节点的抗震连通可靠度是可行的（表4）。设定较优的源点，需要通过分析连通性矩阵，综合评定。如选取39号节点为设定源点，与1号节点为源点时相比，所有汇点抗震连通可靠度平均提高幅度为：Ⅶ度情况下，排名第1位，为10.22%：Ⅷ度情况下，排名第8位，为36.46%，但与排名首位的节点仪相差1.26%：Ⅸ度时，排名仪第25位，为203.55%，但此时汇点的连通可靠度已经普遍较低，排名高位也难以在实质上提高汇点的供气能力。此外，Ⅶ度、Ⅷ度和Ⅸ度情况下，连通可靠度≥0. 800 0的汇点数目分别为
8
1、30和0，与排名首位的节点相比，相差分别仪为
1、6和4。因此，39号节点为较优的设定源点。
较优的设定源点多分布在管线密集的市中心，因此要考虑实际情况，在市区边缘选择合适的节点作为源点，而同时也需要着重对源点周边的网络拓扑结构进行优化。
4 结语
在城市燃气管网抗震连通可靠性分析中利用蒙特卡罗法，迭代5000次即能满足相应的精度要求。
每个汇点与源点的连通可靠性与管道单体的抗震可靠性有关，但单纯依靠增大管径或采用抗震新型管材，经济成本较大。管网各汇点的抗震可靠性与源点的相对位置密切相关，需要综合实际情况来选取地理位置合适的源点。此外，增设新的源点或改善管网的拓扑结构来提高整个系统抗震可靠性也是很有成效的途径，但如何在一定的时间跨度内经济合理地实现这一目的，有待进一步研究。
