1、城市潜在地震次生火灾风险属性区间识别理论 刘晓然 王威 北京建筑大学 北京工业大学抗震减灾研究所 摘 要: 针对城市地震潜在次生火灾风险分析影响因素的复杂性, 首先构建了基于属性区间识别理论的城市地震潜在次生火灾风险评价模型, 并采用综合权重法确定权重向量, 以最大程度地消除各种不确定性;其次从风险的角度综合考虑地震潜在次生火灾系统中致灾因子和孕灾环境的危险性, 承灾体的易损性和暴露性, 及城市地震潜在次生火灾风险的抗灾因子, 构建了城市地震潜在次生火灾害风险评价指标体系, 并制定了相应评价指标的等级标准;最后, 应用属性识别理论对某城市地震潜在次生火灾系统中危险性、易损性、暴露性, 及抗灾能
2、力进行了计算, 并应用证据理论对区域火灾风险性进行了验证分析。结果表明:当城市发生、度地震时, 中等风险以上的区域面积分别达到了50.20%、61.00%、64.90%, 说明针对城市的风险级别较高的区域单元应重点防范并采取相应的措施以降低风险。文章证实了属性区间识别理论及证据理论在城市地震潜在次生火灾风险评估中的应用, 提高了风险评价结果的稳定性和可靠性, 为城市地震潜在次生火灾风险评估提供一种易操作且可行的方法。关键词: 城市; 地震; 次生火灾; 属性区间识别; 证据理论; 最大熵原理; 风险分析; 作者简介:刘晓然 (1988-) , 女, 山东聊城人, 博士, 讲师, 研究方向为地震
3、灾害系统仿真与灾害模拟技术等方面的研究.E-mail:收稿日期:2017-02-14基金:国家自然科学基金项目 (51678017) Potential Risk Analysis of Urban Post-earthquake Fire Based on Attribute Interval Recognition TheoryLIU Xiaoran WANG We Beijing University of Civil Engineering and Architecture; Institute of Earthquake Resistance and Disaster Reducti
4、on, Beijing University of Technology; Abstract: For the influence factors complexity of the urban post-earthquake fire, theory about the attribute interval recognition is proposed for the potential risk analysis of urban post-earthquake fire. Firstly, an attribute interval recognition model is estab
5、lished to eliminate all kinds of uncertainty to maximize based on the principle of maximum entropy; Secondly, the evaluation index system for the urban post-earthquake fire was presented in terms of urban disaster hazard, exposure and vulnerability of disaster bearing body and disaster prevention an
6、d mitigation capability. Meanwhile, the grading standards of evaluation index are also given; finally, potential risk analysis for a citys post-earthquake fire is provided to validate and prove its efficiency and feasibility by the evidence theory and attribute interval recognition theory. The resul
7、ts show that the middle and higher risk area reaches 50. 20%, 61. 00% and 64. 90% respectively when , , degree earthquake occurs. The application of attribute interval recognition and evidence theory for potential risk assessment of urban post-earthquake fire were verified, which improved the stabil
8、ity and reliability of the risk assessment results, to provide a measurable and achievable approach for the potential risk assessment of urban post-earthquake fire.Keyword: urban post-earthquake fire; attribute interval recognition; evidence theory; principle of maximum entropy; risk analysis; Recei
9、ved: 2017-02-14我国地震活动频繁、破坏性大。破坏性地震发生后, 极易引发次生灾害。现代城市是区域政治、经济和文化中心, 燃气管道、压力管道、城市电网等各种管道错综复杂。这都将使人口集中、建筑物集中、生产经营集中、财富集中的现代化城市更容易在震后集中发生特大火灾和燃爆事故, 有时甚至超过地震直接损失。城市地震潜在次生火灾是地震次生灾害最重要的方面, 因此应当引起高度重视。为了有效地减轻地震次生火灾灾害, 减少火灾损失, 有必要进行地震次生火灾风险性分析及风险性区划的基础性研究工作1。目前, 国内外关于城市火灾风险的研究已经取得了一些成果2-5, 日本的“城市等级法”, 美国的“灭火
10、分级制”等, 都是根据区域建筑物、环境密度等来确定区域风险等级, 进而设置消防设施。我国近年才展开了对城市火灾风险评估的研究, 主要方法有统计分析法、指数法、专家调查法、层次分析法等, 但以上这些理论与方法大多是针对城市平时整体火灾风险或危险评估, 对地震次生火灾不是十分适用。向琴琴6、张志华7建立了城市小区地震次生火灾危险性评估模型, 但对于如何考虑地震破坏的影响没有具体可行的好办法。城市地震火灾风险分析是对致灾因子和孕灾环境的危险性、承灾体的易损性和暴露性、城市抗灾防灾能力的分析。以往研究的主要特点是, 确定区域灾害风险评价的指标, 通过数学处理手段进行综合得到风险评价结果。然而, 如何解
11、决灾害风险指标的不确定性及其间的相互影响仍值得进一步研究, 也是关系到风险评价结果准确性的一个重要问题。属性区间识别理论模型8-9是一种以最小代价原则、最大测度原则、置信度准则和评分准则为基础的综合评价方法, 同时最大熵原理可排除人为因素、风险因素等的干扰, 最大化反映评价对象的客观信息。目前以城市区域为研究对象, 以区域地震次生火灾风险评估为研究内容, 建立针对区域地震次生火灾风险评估模型, 还没有完整的、现成的方法可供使用。1 基于属性区间识别理论模型基于以上所述, 本文提出了基于属性区间识别理论及证据理论对城市地震潜在次生火灾风险进行评价与分析, 详细技术路线如图 1 所示。图 1 城市
12、地震火灾风险评价技术路线 下载原图城市地震潜在次生火灾风险评价中, 已知 n 个评价单元的 m 个指标 X= (x1, x2, xn) , 根据 k 个风险等级构造区间矩阵 Iab= (a, bjk) , 式中 j=1, 2, , m;k=1, 2, K。j 和 k 分别是地震潜在次生火灾风险评价指标数目和风险等级评价级别数目9-12。上届风险矩阵为 A=ajk, 下届风险矩阵为 B=bjk。步骤 1:建立上下届标准隶属度矩阵 S=sjk (mK) 和样本值相对于上下届的隶属度矩阵 F=fjk (mn) 。具体计算如下:正指标:负指标:式中:i=1, 2, n;j=1, 2, m, k=1,
13、2, K。步骤二:采用层次分析法和熵权法相结合的综合权重法确定危险、暴露、易损、抗灾防灾能力等级权重向量。式中: j为各项评价指标的权重; j、 j分为第个评价指标的主、客观权重;01 为各评价指标主、客观权重的偏好系数, 步骤三:属性测度区间计算。为消除随机性和模糊性, 引入最大熵原理计算 -ik和 ik:式中:B 一般取 10。步骤四:综合属性测度计算。采用均化系数对属性测度区间进行转换, 得到样本xi属于第 k 类的综合属性测度。式中: 为属性测度区间均化系数。步骤五:评价样本类别识别。风险的危险性、暴露性、危险性、易损性和抗震防灾能力的评价等级确定方法为:另一种算法是采用置信度准则评价
14、样本等级:式中: 为置信度, 值越大越趋向于保守。步骤六:风险等级确定。为有效获得样本的风险评价等级, 采用证据理论 Dempster 合成法则对步骤四计算的综合属性测度进行耦合13:其中:式中:A 为可能的结果的集合;m (A) 为命题 A 的基本信度;k 为归一化系数, 反映证据间冲突的大小。信任函数:式中:B 表示 A 的子集, 则该函数为 上的信任函数, 表示 A 的全部子集所对应的基本信任分配之和。似然函数:表示 A 不否定的信任度, 是所有 A 相交子集的基本信任分配之和。式中:珔 A 表示集合 A 的补集, B 表示和 A 相交不为空集的集合。似然函数与信任函数之差为不确定性的大
15、小。2 城市地震潜在次生火灾风险指标体系城市地震潜在次生灾害风险评价12涉及致灾因子、孕灾环境及承灾体等众多因素。且由于许多因子无统一定量标准, 使得评价指标体系很复杂, 难以操作。从灾害风险生成的动力学角度看, 防灾减灾能力与灾害风险生成的作用方向是相反的, 即特定地区的防灾减灾能力越强, 灾害危险性、易损性和暴露性生成灾害风险的作用力就越会受到限制。本文基于灾害系统理论, 从致灾因子、孕灾环境、承灾体、抗灾能力出发, 遵循系统性、易量化、可操作、普适性的原则, 构建了城市地震潜在次生火灾风险综合评价指标体系 (图 2) 。图 2 城市地震风险评价指标体系 下载原图(1) 危险性指标。地震发
16、生造成的建筑物倒塌、危险品泄露是引发次生火灾的主要原因, 而交通系统和供水系统的震害程度影响到次生火灾救援的便利性, 震害指数越大, 救援的不利因素越多, 火灾蔓延速度越大, 故选取危险源密度、建筑物震害指数和主要生命线系统震害指数作为危险性评价指标。建筑物或主要生命线震害指数, 指一定地区范围内所有建筑物或生命线工程的震害指数的平均值, 即各破坏等级的建筑物或生命线工程所占的比率与其相应的震害程度量化值乘积之和。危险源的级别不同造成灾害程度不同, 重大危险源主要是指有危险品仓库的化工企业、加油站等, 一般危险源主要是指民用危险源。(2) 暴露性指标。城市地震火灾承灾体的暴露性主要体现在人和财
17、富上, 所以选取人口密度和建筑物密度为暴露性因子。(3) 易损性指标。易损性主要是体现承灾体的易损程度, 同暴露性, 承灾体是人和建筑物。易损性人口因子选取人口老少比例 (年龄中大于 65, 小于 14 人口所占总人口比例) , 建筑物因子选取建筑结构类型指数、薄弱区比例 (多为13 层的老旧民宅) 、3、4 级或无耐火建筑比例为易损性指标。其中建筑结构类型指数是指各类建筑结构类型指数乘以相应建筑面积所占权重所得。综合考虑建筑物的耐火等级及抗倒塌能力所得, 各类建筑结构类型指数如表 1 所示。表 1 各类建筑结构类型指数 下载原表 (4) 抗灾能力指标。绿地对次生火灾的蔓延具有阻力作用, 道路
18、的等级及覆盖度对震后火灾救援具有重要意义。消防能力的大小取决于消防有效覆盖度及消防人员的数量, 根据易量化原则, 选取人均绿地面积、道路覆盖度、消防覆盖度 (消防服务面积所占区域面积比) 、消防员人均保护数量为抗灾防灾因子, 其中道路覆盖度取决于区域单元内道路等级、道路宽度、道路冗余度。城市地震潜在次生火灾风险评价指标尚无明确的统一标准, 在参考国家、行业及地方规定或颁布的有关标准基础上, 引用相关文献中确定标准的方法, 结合研究区域经济发展及次生火灾源的实际情况和特点, 将地震次生火灾危险、暴露、易损、风险、抗灾能力等级均分为极低、低、中等、高、极高 5 个等级, 具体等级区间值如表 2 所
19、示13-17。3 实例应用以某城市为研究对象, 为综合考虑城市区域自然环境和建设环境特点、行政管理权限, 从抗震防灾的角度, 将该城市划分为 110 个疏散分区 (图 3) 。因为划分单元数量过多, 根据建设用地功能的不同, 选取不同功能用地作为示例单元进行计算说明, 8-4 工业区、9-6 居住区、12-13 风景旅游区、13-3 商贸区。示例单元及所取值如表 3 所示。表 2 城市地震潜在次生火灾风险指标分区标准 下载原表 表 3 示例单元地震潜在次生火灾风险指标取值 下载原表 表 4 不同烈度下危险性指标权重 下载原表 表 5 基于属性区间识别模型危险性评价结果 下载原表 表 6 不同烈
20、度下城市地震潜在次生火灾危险等级分区数对比 下载原表 图 3 城市疏散分区图 下载原图图 4 城市火灾危险性分区图 下载原图3.1 地震次生火灾危险性评价按照前述属性区间识别模型的计算步骤首先进行危险性计算。首先可判断评价样本数目 n=110, 指标数目 m=6 以及评价等级 K=5。因为不同的地震烈度导致不同的建筑震害指数、交通震害指数、供水震害指数, 按照层次分析法和熵权法组合权重法, 参考文献18确定 =0.6, 可以得到不同地震烈度下危险性指标权重如表 4 所示。采用式 (4) 式 (8) 计算各评价样本的 ik和 ik, 取均化系数为 0.5, 得到综合属性测度, 并根据式 (8)
21、和式 (9) 确定样本所属危险等级, 其中取置信度为 0.5, 同理可得暴露性、易损性、抗震防灾能力等级, 详细计算结果见表 5。由结果绘制地震烈度、度下该城市地震潜在次生火灾的危险性专题图, 见图 4。根据地震烈度、度下城市地震火灾危险等级图的对比, 可得到各烈度下各危险等级的单元个数, 见表 6。极高风险等级单元个数均为 0, 在此不列出。由表 6 可看出地震烈度度时, 区域灾害危险性主要集中在低危险等级以下, 地震烈度、度时, 危险等级高的单元个数及面积明显增加。其中 4-9、8-1、9-5、9-6、9-8、10-1、10-2、10-3、10-6、10-8、11-2、11-8、11-9、
22、12-4、12-8 均为老城区, 老旧房屋、道路失修所占比例多, 导致建筑、道路震害指数大且一般危险源数目较多导致危险等级较高, 属于重点防范区。3.2 地震潜在次生火灾暴露性、易损性、抗震防灾能力评价同理可得基于属性区间识别模型的暴露性、易损性、抗震防灾能力的计算结果如表 7 所示。由结果绘制该城市地震潜在次生火灾的暴露性、易损性及抗震防灾能力专题图层 (图 5、图 6、图 7) 。从图 5 中可以看出, 由于老城区的人口和建筑物密集, 导致老城区分区单元的暴露性较高。由图 6 可以看出, 1-9、4-3、4-4、7-1、7-2、7-5、8-5、9-1、12-6 易损等级较高, 主要是因为城
23、中的较为发达, 导致大量青年进城工作, 城市外围的老少人口较多, 且较不发达地区的房屋多为砖混结构, 建筑结构类型指数高、耐火等级低导致的这些区域单元易损性高。由图 7 可以看出, 2-3、4-3、4-8、7-6、9-4 抗灾防灾能力较差, 主要是道路等级低导致的有效覆盖度不够级消防设置不完善导致。图 5 城市地震潜在次生火灾承灾体暴露性等级图 下载原图3.3 基于证据理论的城市地震潜在次生火灾风险评价城市地震潜在次生火灾风险等于区域危险性与承灾体暴露性、易损性及抗灾防灾能力的耦合, 且风险值随危险性、暴露性、易损性的增大而增大, 随抗灾防灾能力的增大而减小。第一步首先采用危险性、暴露性、易损
24、性、抗灾能力的属性测度分别作为证据 1、证据 2、证据 3 和证据 4, 应用式 (11) 计算冲突系数 k。第二步多证据组合。以以下片区为例, 应用式 (10) 、 (12) 及 (13) 计算各片区在地震烈度分为度、度、度下的质量函数、信任函数和似然函数 (表 8) 。根据烈度、度下城市地震潜在次生火灾风险等级对比, 可得到各烈度下各风险等级的单元个数及各风险等级面积所占比 (表 9) 。图 6 城市地震潜在次生火灾承灾体易损性等级图 下载原图图 7 城市地震潜在次生火灾抗灾防灾能力示意图 下载原图如表 9 所示, 随着地震烈度的增加, 高风险等级的单元个数及所占面积比增加, 当地震烈度为
25、时, 中等及以上风险等级面积已经达到 60%。该城市地震潜在次生火灾风险等级较高, 经过结果对比, 判定 1-9、2-1、4-9、8-8、10-2、10-3、10-11、11-5、11-6、12-8、12-12 单元为高风险区, 其中 4-9、10-2、10-3、10-11、11-5、11-6、12-8 属于老城, 人口密集、老旧民房较多, 一般危险源数目多, 危险性、暴露性、易损性等级高, 消防配置不够导致风险等级高。1-9、2-1、8-8 为工业区或新建区, 重大危险源数目较多, 主要是危险性等级较高, 消防设施不完善导致的风险等级高。根据分区图采取相应措施降低灾害的危险性、承灾体的暴露性
26、与易损性, 提高城市的抗震防灾能力, 进而降低城市地震潜在次生火灾风险。对于中等风险区, 也应采取相应的措施进一步防范。表 7 暴露性、易损性、抗震防灾能力评价结果 下载原表 表 8 基于证据理论模型风险性评价结果 下载原表 表 9 不同烈度下城市地震潜在次生火灾风险等级对比 下载原表 4 结论探讨了属性区间识别理论及证据理论在城市地震潜在次生火灾风险评估中的应用。在深入分析城市地震潜在次生火灾风险成因的基础上, 从城市地震发生后的致灾因子、孕灾环境的灾害状况和承灾体即抗灾防灾能力四个方面构建了城市地震潜在次生火灾风险评价指标体系。以某城市为例, 运用基于最大熵理论的属性区间识别理论模型对城市
27、分区单元不同地震烈度下危险等级、暴露性等级、易损等级、抗灾防灾能力等级做出科学合理的计算, 并由风险等级分区矩阵得到各评价单元的风险等级, 评价结果比其他方法更加客观。根据结果找出区域潜在风险, 采取相应措施, 与实际情况也较为吻合, 提高了风险评价结果的稳定性和可靠性, 为城市地震潜在次生火灾风险评估提供一种易操作且可行的方法。参考文献1郭惠.城市地震火灾风险评估指标体系研究J.中国安全生产科学技术, 2013, 9 (3) :56-59. 2USFA.Information on the Risk, Hazard and Value EvaluationR.1999. 3Michael S
28、 Wright.Dwelling Risk Assessment ToolkitR.1999. 4杜霞, 张欣, 刘庭全, 等.国外区域火灾风险评估技术及应用现状J.消防科学与技术, 2004, 23 (2) :137-139 5连旦军, 董希琳, 吴立志.城市区域火灾风险评估综述J.消防科学与技术, 2004, 23 (3) :240-242. 6向琴琴.基于 GIS 的地震次生火灾危险性评估及动态蔓延模拟D.北京:北京工业大学, 2005. 7张志华.城市区域地震次生火灾危险性评估J.消防管理研究, 2008, 27 (8) :602-605. 8刘松涛, 刘文利, 孙旋, 等.基于能力和脆弱性分析的火灾风险综合评价体系研究J.科技导报 2009, 27 (22) :76-80. 9李群, 宁利.属性区间识别理论模型研究及其应用J.数学的实践与认识, 2002, 32 (1) :50-54. 10邹强, 周建中, 周超, 等.基于最大熵原理和属性区间识别理论的洪水灾害风险分析J.水科学进展, 2012, 23 (12) :323-333. 11任政, 郝振纯.基于属性区间识别理论的水资源可持续利用评价J.水电能源科学, 2009, 27 (3) :22-24