中南大学资源与安全工程学院，湖南长沙，410083)
LIU Aihua', WU Chao′, XU Wenbin'
(1. School of Environmental Science and Engineering, Guangdong University of Technology, Guangzhou, 510006, China;
2. School of Resources and Safety Engineering, Central South University, Changsha, 410083, China)
基金项目(Foundation item): 国家自然科学基金资助项目(51508110)(Project(51508110) supported by National Natural Science Foundation of China)通信作者：徐文彬，教授，博士，从事安全评价、事故模拟和应急方面的研究；E-mail：xuwenbin@tom.com
城市生命线系统的概念是20世纪70年代DUKE等[]在考察电力、交通、供水及煤气等系统震害以及对整个城市功能的影响后提出来的，认为这些基础设施对城市的经济和社会生活具有特别重要的作用，因而称之为“城市生命线系统工程”，对生命线系统的高度依赖已经成为现代化城市的重要特征之一。随着城市化程度的不断提高、城市功能的日益复杂化和多样化，与城市生命线相伴出现的灾害隐患越来越多，而且城市生命线系统的灾害事故具有破坏严重、耦合放大性、次生灾害突出、社会影响大等特点。灾害风险研究是构建综合防灾减灾体系的基础和依据。而承灾体的脆弱性评估又是灾害风险评估的核心问题，也是国际上各种灾害风险评估方法共同探讨的热点问题[]。为了预防和控制城市生命线灾害事故，国内外许多学者对城市生命线的脆弱性也进行了较深入的研究，主要集中在2个方面：一是地震灾害作用下城市生命线系统的脆弱性分析和评估[3-5]；二是供电、交通等生命线子系统在地震、洪水等单灾种作用下脆弱性评估方法的研究[6-8]。城市生命线系统具有复杂网络结构特点，且不同类型系统在功能上具有耦联性。这种耦联性由造成结构破坏的随机变量或随机过程之间的相关性引起，也可由结构破坏的传递性引起。因此，城市生命线脆弱性评估如果把单功能子系统割裂开来进行评估，往往会因为忽略系统这种耦联性而使得评估结果出现偏颇，单一灾种条件下的脆弱性分析则会因忽略灾害的链式效应和放大效应而出现表征的局限性，不利于城市综合防灾减灾体系的构建。针对上述问题，本文作者在对脆弱性和灾损敏感性的概念进行界定的基础上，提出用均衡熵表征单功能网络结构的脆弱性，并应用脆性熵理论求解出各生命线子系统的脆性关联程度，结合城市生命线的质量灾损敏感性，构建生命线系统的灾损敏感性评估模型。
1城市生命线灾损敏感性评估指标体系的构建
受一定致灾强度的灾害打击后易于受损的程度，它反映了承灾体因自身的物理特性而对致灾因子破坏作用的抵御能力。物理暴露性可通过承灾体的管理部门获取相应的统计数据而确定，因此灾损敏感性是对各承灾体脆弱性分析的重点和难点。
在布设城市生命线时，许多敷设在地下，受地震的影响较大，因此生命线本身的抗震性能是影响其质量灾损敏感性的重要因素，另外使用年限、工程质量、以及生命线的网络结构等也是影响其质量灾损敏感性的主要因素。根据王飞[10]的研究成果，1985年以前投资建设的基础设施因使用年限较长而抗灾性能较差；而1994年以后投资建设的基础设施抗灾性能相对较好；近5年基础设施投资额年均增长率则反映了目前的城市生命线设施状况和发展水平。
城市生命线以整体网络结构形式支持着城市的系统功能，因此衡量生命线网络结构的灾损敏感性既要考虑各子系统的网络可靠性，也要考虑各子系统之间的关联性对整体网络结构所造成的影响。单功能网络结构的可靠性可从网络的拓扑结构性、系统的供给能力等多方面进行度量[9, 11]。生命线网络结构中脆性的存在，是生命线系统事故具有耦合放大特性的根本原因，因此，在对生命线系统的灾损敏感性进行评估时，必须对其脆性进行评估。由此，可建立城市生命线灾损敏感性评估的指标体系，如图1所示。
图1城市生命线灾损敏感性评估指标体系 Fig. 1 Index system of damage sensitivity evaluation on urban lifeline
许多研究把承灾体的灾损敏感性和脆弱性等同为一个概念，固然灾损敏感性和脆弱性都是相对承灾体而言的，但两者的内涵是不一样的。脆弱性是指在一定灾变条件下受灾体的抗御能力及其损毁程度。脆弱性包括承灾体的物理暴露性V、灾损敏感性V和应灾能力V.3个方面[9]。物理暴露性是指承灾体暴露在致灾因子影响范围之内的数目或价值；区域应灾能力是指区域社会为了预防某种灾害的袭击而对承灾体采取的各种防御措施的力度；灾损敏感性是指承灾体在遭
2单功能网络的灾损敏感性评估
城市生命线系统结构庞大、功能复杂。在结构上，呈现分块又分层的立体结构；在运行过程中，呈现出明显的非线性和高维特征，因而城市生命线系统是一个动态多变的复杂网络系统[12-13]。城市生命线网络已经发展成为世界上最复杂的人造网络之一，而每个单功能子系统在结构上又有自己独立的网络结构。
2.1网络节点的结构重要性
在每个单功能网络结构中，每个节点的地位是不同的。以交通网络为例，某些交通节点在发生事故时只会影响局部的交通畅通，而某些交通节点在发生事故时会导致整个交通系统的瘫痪。为了衡量这些节点在网络拓扑结构中的重要程度以及影响能力，引入中介中心性的概念[11,14]。在网络结构中任意节点k的中介中心性gk定义为：设某网络拓扑结构中的一条路径，开始于节点为i，中间经过一系列的节点和边，终止于节点j，为其中的一个中间节点[15-171，设c(i，j)为节点i到节点i所有最短路径的总数，在这些路径中，通过节点k 的路径数为chi，j)，则记gk为gii，j),，可表示为
(1)
gri，j)即表示节点k在网络中的结构重要程度。
2.2网络节点的功能重要性
在网络结构中，有些路径输送的资源量很大，该路径以及路径上的节点一旦遭受破坏，对系统的影响很大。因此衡量网络节点重要性的另一个指标就是节点所在路径传递的信息流，称之为功能重要性，其实际是指节点所在路径传递介质的能力，可用连接边强度来表征。对于供水、供气等资源供给的管网系统，连接边的连接强度就是管网中的流压。设任意2个连通节点k和1的连接边强度用uu表示，而且在节点i与之间存在连接强度大于0的一个边的序列，把i和i连接在起来[10-1]:
该序列被称为节点i与j的最短路径，记为Pr。
(2)
2.3单功能网络的灾损敏感性评估
节点的结构重要性和功能重要性综合重要性都影响着整个网络的灾损敏感性，因此，要综合考虑。引入综合重要度H予以衡量，定义为
(3)
如果网络系统中存在少数几个综合重要度极高的节点，则该网络称之为“不均衡”网络，称综合重要度极高的节点为“核心节点”。在此网络中，核心节点一旦失效，会导致系统资源的供给量远远小于正常的需求量，从而影响广泛，严重的可能导致整个网络崩溃或瘫痪。对于不存在核心节点的“均衡”网络，就不会因为某个节点的失效而影响系统的整体运行。可见网络节点分布的均衡性能有效反映单功能网络的灾损敏感性[10]。德国物理学家克劳修斯所提出的热力
学概念“熵”也被广泛应用于系统失序现象的表征。熵越大，表示对象越无序，而系统越均衡。由此可推，在网络结构中，同层节点综合重要度的熵越大，表明网络也就越均衡[11]，因此，可用“熵”来反映节点在网络结构中分布的均衡性。在此引入“均衡”的概念对生命线单功能网络系统的灾损敏感性进行评价在理论上具有可行性。
设某单功能网络系统共有n个节点，从源点到汇点的最短路径共有h条。分别统计计算每个节点的结构重要度，并把具有结构重要度相同节点组成节点集，设共有k个节点集，第i个节点集中有m个有节点，各个节点集的综合重要度分别为Hi，H2，·，Hk。首先对各节点集节点的综合重要度进行归一化处理，设归一化处理后各节点的综合重要度分别为 Hi,H2,.,H'.,H,则
(4)
其中：m;为第j个节点集中的节点数。
六时，系统埔 E 取得最大根据熵原理，当H；：
n值，即
(5)
那么网络系统实际的均衡的计算式为
(6)
如果网络系统实际均衡熵与最大摘的比值越接近1，表明网络的均衡性越好，网络系统的可靠性也就越高；反之，网络系统的可靠性越低，灾损敏感性也越大。因此，网络系统的灾损敏感性可表示为
(7)
其中：LG-vu 为单功能网络的灾损敏感性指数。
3基于脆性整体网络结构的灾损敏感性评估
3.1城市生命线系统间相互影响
城市生命线系统中的各个子系统之间存在相互作用，且作用方式各异，关系非常复杂。以各子系统的相互作用机理为基础，对相互作用类型进行科学分类是厘清其复杂关系的有效方法。因此，从作用机理入手，结合子系统间相互影响的事故案例，将生命线系
统的相互作用划分为4种类型：
1）功能型相互作用，是指一个子系统功能的正常发挥依赖于与之关联的其他子系统功能的正常发挥，也就是当与之关联的子系统功能失效时，会导致该系统的功能无法正常发挥。典型的案例有，供电中断会使得水泵无法运转，水泵的停止运转则会进一步导致供水系统的功能失效；天气寒冷导致供暖增加，从而导致燃气紧张等。
2）布设型相互作用，即因为子系统的相邻布设，导致一个系统受损时，而另一个系统也被破坏。例如，沿道路敷设的城市燃气、供水等地下管道受到破坏时，往往会引起道路系统的损坏，同样在道路遭受破坏时，也可能会导致燃气、供水等地下管道的损坏。
3）替代型相互作用，即-个系统的功能可用与之关联的其他子系统的功能进行替代，因此当系统失效时，会造成与其功能替换子系统的过量需求。例如，供气系统中断，作为可替换能源的电力的消耗将明显增加。
4）恢复型相互作用，即破坏后在功能恢复与修复进程中会彼此影响，当多种生命线系统失效时，其中某种生命线系统会影响其他生命线系统的修复进程。
例如电力系统破坏，其他系统恢复都会变得很困难。
相互作用的类型划分是相互影响机理的具体化。根据如上分类方式，对供电、供气、供水、交通、通讯五个典型生命线系统之间的相互作用关系进行分析，它们之间在功能严重失效情况下彼此间的相互作用类型及对应的影响方式如表 1 所示[18-19]。
3.2基于脆性炳的生命线系统脆性分析
脆性是指材料在外力作用下(如拉伸、冲击等)仅产生很小的变形即断裂破坏的性质[20]。后来人们发现，复杂系统也表现出这一特性，因此，逐步把脆性的定义和研究方法引申到对复杂系统的特性研究当中。其相应的定义为：复杂系统在内、外因素的扰动作用下，出现部分崩溃，进而影响到与之关联的其他部分，并产生连锁崩溃效应，最终导致整个复杂系统崩溃的特性，称之为脆性[21]。
可由上述研究可知：城市生命线系统是一个动态的复杂系统，各子系统的关联性必然会导致系统内脆性的存在。因此，在对城市生命线的灾损敏感性进行评估时，必须考虑系统的脆性，以免割裂或忽视了各系统之间的关联。根据脆性理论思想，结合生命线子系统的相互影响分析，利用熵理论对城市生命线子系
表1生命线系统功能严重失效时的相互作用特性 Table 1 Interaction characteristics of lifeline system in serious failure of function
注：·一布设型相互作用；o一功能型相互作用；△一恢复性相互作用；口一替代型相互作用。
统之间的脆性关联进行定量处理和计算，再结合单功能子系统的灾损敏感性，对城市生命线系统的整体灾损敏感性进行评价。
1）脆性联系。如果一个子系统X因为干扰而发生事故且不能正常工作，与之脆性关联的其他子系统必将会受到影响，表征其状态的变量也将随之改变。假设 Y因与 X脆性相关而受到影响，表征 Y的 n个状态向量中，变量 以(1 
(8)
(9)
脆性对立的定义式为
(10)
脆性波动熵的定义式为
(11)
实际上，子系统Y在X的关联作用下，应该是脆性同一、脆性对立和脆性波动的综合。因此子系统 X与子系统Y的脆性联系熵为
(12)
式中：のa，Wb，のc分别为脆性同一、脆性对立和脆性波动的权重系数，采用燕尾突变级数法可以确定相应的值:
(13)
在系统中，存在一定的概率分布组合，使 Hx达到最大，表明此时子系统Y受到X的影响最大，而崩溃的可能性也最大。与子系统 Y类似，也可以定义其他子系统脆性联系摘的相关概念。
2）生命线系统的相互作用。根据脆性联系摘的定义，在城市生命线系统当中，各子系统间的布设影响可表征为脆性同一，替换作用可表征为脆性对立；而
功能影响和恢复影响可表征为脆性波动。在事故状态下城市生命线子系统间的脆性关联特征可用布尔型矩阵来表示，矩阵中的0表示彼此间没有影响，1表示彼此间存在相互影响，那么不同类型影响方式的布尔型矩阵如下。
城市生命线子系统间的脆性联系熵计算过程如下。
1）首先对每个矩阵的列求和，得到它们耦合之后的综合矩阵：
2）对综合矩阵的每列再求和，就可得到各个子系统与其他子系统的影响关联：供电系统为9，供气系统为8，供水系统为7，交通系统为7，通讯系统为7。
3）最后再计算各子系统的脆性同一、脆性对立和脆性波动所对应的概率Pa，Ps和P。及其相应的熵Ha， Hp和H.等参数，结果如表2所示。
由表2可知：各生命线子系统的脆性关联从小到
表2脆性关联计算结果 Table 2 Computer results of britleness connection
大排列顺序为：供水系统，交通系统，供气系统，供电系统，通讯系统。表明与其他子系统脆性关联最大的是通讯系统，是灾后应该抢先修复的子系统，以此来降低由于通讯系统与其他系统关联而带来的连锁影响。
4城市生命线的灾损敏感性计算
4.1生命线网络结构的灾损敏感性评估
作为单功能网络系统的城市生命线子系统，网络结构本身存在一定的灾损敏感性，而脆性关联的存在会导致城市生命线系统整体的灾损敏感性进一步增加，因此对生命线系统的整体灾损敏感性进行分析时，既要考虑各个子系统的灾损敏感性，还必须考虑脆性关联所引发灾损敏感性的增量。在对整体灾损敏感性进行评估时，可把每个子系统的灾损敏感性视为整个系统灾损敏感性一种表现途径，即各个子系统类似并联的几个系统，另外，把每个子系统自身崩溃而引发的其他系统崩溃的情况也统计在内，由此生命线系统整体网络结构的灾损敏感性可表示为
(14)
其中：Lvlae)为生命线系统网络结构的灾损敏感性指数；LG-Vuli 为子系统的网络灾损敏感性指数；H；为各系统的脆性关联熵的归一化值。
(15)
由上述分析可知：随着城市化程度的提高，城市各子系统之间的功能、恢复、替换等交互影响将会更更加复杂化，因此生命线各子系统之间的脆性关联熵只会进一步增加，而很难降低。不过可以通过方案设计使布设更加合理从而降低彼此间的脆性关联度。通过网络结构优化设计，在保证功能的基础上尽可能提高网络的综合均衡性则是降低生命线整体网络结构灾损敏感性最有效的措施。
4.2城市生命线的灾损敏感性计算
综合考虑各项指标，城市生命线灾损敏感指数的计算方法如下：
V, = αquality Lulqualt) + net * Lvuline)
(16)
式中：αquality 和αnet 分别为城市生命线质量敏感性指数、生命线系统网络结构的灾损敏感性指数的权重系数，可采用层次分析法得出。