发明名称 一种石油化工企业承灾体脆弱性评估方法
摘要 本发明提供了一种石化企业承灾体脆弱性评估方法,所构建的石化企业承灾体脆弱性评估模式可指导石化企业火灾、爆炸及有毒物质泄漏事故的控制与预防,降低石化企业事故风险。该方法构建的评估模式考虑了主客观因素的影响,评估结果更为客观合理,提供的决策建议具有较大的实际应用价值。
申请公布号 CN105825320A 申请公布日期 2016.08.03
申请号 CN201610133779.3 申请日期 2016.03.07
申请人 天津工业大学;天津市华腾理工安全科技咨询有限公司 发明人 卢丹丹;宋文华;张桂钏;程宝东
分类号 G06Q10/06(2012.01)I;G06Q50/06(2012.01)I 主分类号 G06Q10/06(2012.01)I
代理机构 代理人
主权项 一种石化企业承灾体脆弱性评估方法,其中所述方法的步骤主要包括:①确定研究区域按照重大危险源的划分方式,进行研究区域的划分;②资料的收集与统计应收集的关键信息:研究区域工艺情况、研究区域内人口、装置、建(构)筑物分布情况、研究区域内原材物料情况、研究区域内设备、设施情况、研究区域内消防设施配备情况、研究区域内安全防护措施采取情况、企业安全管理状况等;③危险源和承灾体类型辨识与统计分析a.辨识研究区域内原材物料的品种、数量、储量及运行过程中的存在部位,依据标准《危险化学品重大危险源辨识》(GB18218‑2009)分析研究区域内物质的固有危险性;b.分析生产、储存、运输过程中可能发生的火灾、爆炸、有毒有害事故及事故后果影响,依据《危险化学品重大危险源监督管理暂行规定》对研究区域进行重大危险源分级;c.辨识研究区域内可能遭受火灾、爆炸、有毒物质泄漏事故的人员、装置及建(构)筑物等各类承灾目标的数量及分布,分析事故发生时承灾目标的受伤害程度;④指标体系的构建及权重的确定考虑到石化企业内人员、装置及建(构)筑物为主要承灾目标,可构建目标层为承灾体脆弱性总体指标、准则层为人员脆弱性指标、装置脆弱性指标、建(构)筑物脆弱性指标的指标体系,各因素层指标可结合承灾目标的脆弱性影响因素分析和研究区域的实际情况进行筛选;a.层次分析法确定权重对专家的评分结果进行加权计算,并采用1~9及其倒数的标度方法<sup>[12]</sup>进行两两比较,得到判断矩阵:矩阵正规化:<maths num="0001"><math><![CDATA[<mrow><mover><msub><mi>a</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><msub><mi>a</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>a</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></mrow></mfrac><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000127502110000011.GIF" wi="1498" he="169" /></maths>求和向量:<maths num="0002"><math><![CDATA[<mrow><msub><mi>W</mi><mi>i</mi></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mover><msub><mi>a</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>&OverBar;</mo></mover><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000127502110000012.GIF" wi="1568" he="116" /></maths>求权向量:<maths num="0003"><math><![CDATA[<mrow><mover><msub><mi>W</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><msub><mi>W</mi><mi>i</mi></msub><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>W</mi><mi>i</mi></msub></mrow></mfrac><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000127502110000013.GIF" wi="1543" he="142" /></maths>求最大特征根:<maths num="0004"><math><![CDATA[<mrow><msub><mi>&lambda;</mi><mi>max</mi></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mi>i</mi><mi>n</mi></munderover><mfrac><msub><mrow><mo>&lsqb;</mo><mi>A</mi><mover><msub><mi>W</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>&rsqb;</mo></mrow><mi>i</mi></msub><mrow><mi>n</mi><mrow><mo>(</mo><mover><msub><mi>W</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000127502110000014.GIF" wi="1475" he="125" /></maths>检验公式:<maths num="0005"><math><![CDATA[<mrow><mi>C</mi><mi>I</mi><mo>=</mo><mfrac><mrow><msub><mi>&lambda;</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>-</mo><mi>n</mi></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000127502110000015.GIF" wi="1588" he="102" /></maths><maths num="0006"><math><![CDATA[<mrow><mi>C</mi><mi>R</mi><mo>=</mo><mfrac><mrow><mi>C</mi><mi>I</mi></mrow><mrow><mi>R</mi><mi>I</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000127502110000021.GIF" wi="1558" he="153" /></maths>b.熵权法确定权重由已获得的m个评价对象的n个评价指标初始数据构建判断矩阵R,判断矩阵如下:<maths num="0007"><math><![CDATA[<mrow><mi>R</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>R</mi><mn>11</mn></msub></mtd><mtd><msub><mi>R</mi><mn>12</mn></msub></mtd><mtd><mn>...</mn></mtd><mtd><msub><mi>R</mi><mrow><mn>1</mn><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>R</mi><mn>21</mn></msub></mtd><mtd><msub><mi>R</mi><mn>22</mn></msub></mtd><mtd><mn>...</mn></mtd><mtd><msub><mi>R</mi><mrow><mn>2</mn><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>R</mi><mrow><mi>n</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>R</mi><mrow><mi>n</mi><mn>2</mn></mrow></msub></mtd><mtd><mn>...</mn></mtd><mtd><msub><mi>R</mi><mrow><mi>n</mi><mi>m</mi></mrow></msub></mtd></mtr></mtable></mfenced><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2...</mn><mo>,</mo><mi>n</mi><mo>;</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>m</mi><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000127502110000022.GIF" wi="943" he="256" /></maths>权重的计算矩阵正规化:<maths num="0008"><math><![CDATA[<mrow><msub><mi>Y</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>r</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>r</mi><mrow><mi>m</mi><mi>i</mi><mi>n</mi></mrow></msub></mrow><mrow><msub><mi>r</mi><mi>max</mi></msub><mo>-</mo><msub><mi>r</mi><mi>min</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000127502110000023.GIF" wi="1564" he="118" /></maths>信息熵的计算:<maths num="0009"><math><![CDATA[<mrow><msub><mi>f</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn>1</mn><mo>+</mo><msub><mi>Y</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub></mrow><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mrow><mo>(</mo><mn>1</mn><mo>+</mo><msub><mi>Y</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000127502110000024.GIF" wi="1465" he="178" /></maths><maths num="0010"><math><![CDATA[<mrow><msub><mi>E</mi><mi>j</mi></msub><mo>=</mo><mo>-</mo><mfrac><mn>1</mn><mrow><mi>ln</mi><mi> </mi><mi>n</mi></mrow></mfrac><mrow><mo>(</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>f</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mi>ln</mi><mi> </mi><msub><mi>f</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000127502110000025.GIF" wi="1527" he="119" /></maths>熵权系数即权重的计算:<maths num="0011"><math><![CDATA[<mrow><msub><mi>u</mi><mi>j</mi></msub><mo>=</mo><mfrac><mrow><mn>1</mn><mo>-</mo><msub><mi>E</mi><mi>j</mi></msub></mrow><mrow><mi>m</mi><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>E</mi><mi>j</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000127502110000026.GIF" wi="1525" he="158" /></maths>c.组合权重的计算组合权重的计算即对层次分析法和熵权法得到的权重进行复合,组合权重计算公式为:<maths num="0012"><math><![CDATA[<mrow><msub><mi>T</mi><mi>i</mi></msub><mo>=</mo><mn>0.5</mn><msub><mi>u</mi><mi>i</mi></msub><mo>+</mo><mn>0.5</mn><mover><msub><mi>W</mi><mi>i</mi></msub><mo>&OverBar;</mo></mover><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000127502110000027.GIF" wi="1551" he="63" /></maths>⑤模糊综合评价a.确定因素集根据已构建的指标体系可建立脆弱性评估因素集u={u1,u2,u3};b.建立评价等级集依据脆弱性指数划分为等级,不脆弱等级表示优秀,不太脆弱等级表示良好,一般脆弱等级表示系统处于临界状态,应对系统薄弱部位加强管控,比较脆弱表示承灾体脆弱性程度较高,应立即采取措施加以改进,非常脆弱表示承灾体脆弱性程度极高,应立即采取措施并进行重点管控;c.构建评语集,设计专家问卷调查,确立模糊评价矩阵;根据因素等级划分以及承灾体脆弱性影响因素分析,设计各因素的评语集和专家问卷调查表,统计专家对各因素的评判意见,确立模糊评价矩阵R:<maths num="0013"><math><![CDATA[<mrow><msub><mi>R</mi><mrow><mn>1</mn><mi>i</mi></mrow></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>r</mi><mn>11</mn></msub></mtd><mtd><msub><mi>r</mi><mn>12</mn></msub></mtd><mtd><mn>...</mn></mtd><mtd><msub><mi>r</mi><mrow><mn>1</mn><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>r</mi><mn>21</mn></msub></mtd><mtd><msub><mi>r</mi><mn>22</mn></msub></mtd><mtd><mn>...</mn></mtd><mtd><msub><mi>r</mi><mrow><mn>2</mn><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>r</mi><mn>31</mn></msub></mtd><mtd><msub><mi>r</mi><mn>32</mn></msub></mtd><mtd><mn>...</mn></mtd><mtd><msub><mi>r</mi><mrow><mn>3</mn><mi>m</mi></mrow></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>r</mi><mrow><mi>n</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>r</mi><mrow><mi>n</mi><mn>2</mn></mrow></msub></mtd><mtd><mn>...</mn></mtd><mtd><msub><mi>r</mi><mrow><mi>n</mi><mi>m</mi></mrow></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FSA0000127502110000031.GIF" wi="508" he="322" /></maths>其中rij(j=1,2,...,m)表示第Ui个因素对应评语集V的模糊判断;d.多级模糊综合评价脆弱性评估模式为多层次模型,需构建多级模糊综合评价,评价过程为:求二级模糊评价矩阵:R<sub>1</sub>=T<sub>i</sub>оR<sub>1i</sub>            (12)准则层脆弱性指数:<maths num="0014"><math><![CDATA[<mrow><msub><mi>f</mi><mi>i</mi></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mi>i</mi><mi>n</mi></munderover><msub><mi>R</mi><mi>i</mi></msub><msub><mi>V</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000127502110000032.GIF" wi="1390" he="98" /></maths>一级模糊评价矩阵确定:D=TоR<sub>1</sub>=Tо(T<sub>i</sub>оR<sub>1i</sub>)            (14)目标层脆弱性指数:<maths num="0015"><math><![CDATA[<mfenced open='' close=''><mtable><mtr><mtd><msubsup><mi>D</mi><mi>i</mi><mi>l</mi></msubsup><mo>=</mo><msub><mi>D</mi><mi>i</mi></msub><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>D</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>f</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msubsup><mi>D</mi><mi>i</mi><mi>l</mi></msubsup><msub><mi>V</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>16</mn><mo>)</mo></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FSA0000127502110000033.GIF" wi="1412" he="208" /></maths>
地址 300387 天津市西青区宾水西道399号