危险品具有爆炸、易燃、腐蚀性等性质，在运输、贮存和保管的过程中，极易引起人身伤亡和财产损失，因此需要特别防护。危险货物运输作为一种特种运输，需拥有能保障危险品安全运输的设施设备，并通过国家相关职能部门审核批准，才有资格进行危险品运输。据最新统计，我国年危险品运输量约10亿t，专业车辆有31万辆，从业人员达120万人。危险品运输事故一方面威胁人类健康、财产安全，另一方面会对道路、环境、交通等造成重大影响。因此，危险品运输界如何降低运输过程中风险事故发生频率是很具现实意义的问题。

针对危险品运输风险评价，国内外学者进行了相关研究。DONG等^[1]基于理论方法建立了多层级危险品运输安全评估指标体系。Yuliia等^[2]为确定铁路运输危险品期间的环境风险及其管理原则，提出了一种能够优化铁路运输环境管理的新功能策略，强调了优化环境安全监测和组织环境风险管理的重要性。HUANG等^[3]基于尖端灾难模型建立了铁路危险货物运输系统，并对铁路危险品运输系统的风险状态变化进行分析。杨婷等^[4]基于N-K模型从人-机-物-环-管5类风险因素定量分析了道路危险品运输系统耦合风险，得出耦合风险与参与耦合风险因素有关。吴金中等^[5]应用模糊综合评价法，建立了危险品道路运输风险评价模型。黄文成等^[6]利用道路危险品安全事故数据对道路危险品运输系统进行耦合协调度分析，建立了危险品道路运输系统风险评价方法。任常兴等^[7]在基于NFPA进行危险品运输风险数字分级的基础上考虑了现实风险因素对风险等级评定的影响，建立了危险品运输风险分级指数程序。陈跃等^[8]在现有风险分级指数模型的基础上提出了基于指数评价模型的集成化风险管控流程并验证了其有效性和可行性。左博睿等^[9]基于模糊推理技术的Fuzzy-FRAM模型，克服了危险品运输事故分析中无法定量分析的缺陷，得到了功能的风险指数。贺政纲等^[10]考虑了危险品运输的时变特性，构建了人口风险评估模型，得出最佳出发时间，为有效降低危险品运输人口风险提供新途径。李升朝^[11]基于AHP和熵权计算方法构建了模糊综合评价矩阵, 对危险品道路运输风险进行评价。

综上所述，目前研究成果多数集中在评价指标、评价模型的定性分析和风险影响因子重要度的定量分析，而未考虑评价指标存在的模糊性与随机性问题。基于此，本研究引用熵权法获取评价指标权重，并在云理论的基础上，构建基于熵权二维云模型的危险品运输风险评价方法，为保障危险品运输安全，降低危险品运输风险事故提供新的思路与方法。

熵权法是一种高精度、强客观性的较成熟赋权方法，熵值可以判断评价指标的变异性。若某个指标的熵值越小，则指标变异程度越大，信息量越多，权重越重要。但若某个指标的值均相等，则该指标不可作为综合评价的评价指标。评价指标的赋权步骤如下：构造设有待评价项目m个，评价指标n个的原始判断矩阵R=(r_ij)_m×n，第j个指标的熵值计算为：

(1)

式中，P_n为观权重向量；k为常数；m为评价项目数量；r_ij为指标j在项目i下的评分；e_j为指标j的熵值。

指标j的熵权计算为：

(2)

式中ω_j为指标j的熵权。

云模型的概念最早由李德毅在1995年提出，自提出至今，在数据挖掘、智能控制、决策分析等多个方面得到广泛应用^[12]。二维云则是在一维正态云的基础上提出的思想，用以表示定性与定量不确定转换的过程。本研究运用该模型处理危险品运输中风险危害和风险概率两大因素相互作用下的不确定性问题。定义二维云数字特征为：(期望值，熵，超熵)。设F为服从正态分布的二维随机函数，其中Ex和En_x为期望值，Ey和Ey_n为标准差，He为超熵；μ为确定度，故称满足式(3)的云模型为二维正态云模型^[13]。

(3)

危险品运输风险是各项风险因素综合作用导致的，为了采取有效的安全管理、控制措施应对风险，本研究根据我国道路危险品运输现行的相关标准《生产过程危险和有害因素分类与代码》(GB/T13861—2009)、《公路物流主要单证要素要》(GB/T 33458—2016)并参考对公路危险品运输风险识别及风险发生后应急处理办法具有成熟理论国家的现行标准《危险品公路运输-运输过程中的紧急处置和覆盖-特性、尺寸和装填》(ABNT NBR 7503—2017)，从人员因素、货物与设备因素、环境因素、管理因素4个方面辨识危险源。人员因素包含人员健康、人员素质和人员专业技能；危险品货物性质、运输设备、应急预防设备为货物与设备影响因子；环境影响因素包含运输途中道路特征、天气状况、交通状况和影响人员分布情况；安全管理、风险应急预案、运输时间以及运输路线是管理范畴的影响因子。危险品运输风险评价指标体系见图 1。

图 1 危险品运输风险评价指标体系 Fig. 1 Dangerous goods transport risk evaluation indicator system

图选项

邀请专家采用打分法对评价对象的各项指标进行危害等级和概率等级评定，分值区间为[0, 10]，分值精确到小数点后1位。通过计算得到危害等级和概率等级的风险云，即二级风险云，风险云数字特征计算为：

(4)

式中，q为样本数；Ex为样本期望；S²为样本方差；x_k为第k位专家的评分；En为熵；He为超熵。

依据从严从高的风险等级判定原则将[0, 10]分成4个子区间，与之相对应的将评价指标的风险等级、危害等级、概率等级划分为Ⅰ~Ⅳ级，标准云数字特征计算为：

(5)

式中，Ex，En，He分别为标准云的期望、熵和超熵；S_j^min和S_j^max分别为区间j的上、下限值。

为评价危险品运输的风险等级，将二级风险云矩阵与二级权重矩阵合成运算，得到一级风险云矩阵，再将其与一级权重矩阵合成运算，最终得到综合风险云。合成运算为：

(6)

式中，Ex′为一级风险云的期望；En′为一级风险云的熵；He′为一级综合风险云的超熵。风险等级数字特征见表 1。

运用MATLAB编程实现正向云发生器并生成风险云图，得到危险品运输风险等级，同时得到各项评价指标风险等级^[14]。

为实际、客观的判定风险等级，提出一种数学计算方法表示实际、标准风险云的贴近程度。其中，相近度越大，评价结果越相似。运算公式为：

(7)

式中，L为风险等级相近度；Ex，Ex′分别为标准云的危害、概率等级期望；Ex，Ex′分别为实际风险云的危害、概率等级期望^[15-18]。

某危险品运输公司创立于2010年，是一家以危险品运输为特色集运输、仓储、配送、货运代理为一体的国家3A级综合物流企业。该企业先后通过了ISO9001国际质量管理体系认证、ISO14001环境管理体系认证、OHSAS18001职业健康安全管理体系认证、和RSQAS道路安全管理认证，被授予“中国危险品物流百强企业”称号。企业针对化危险品运输现有化工专家2名，国家安全注册工程师3名，专业化工操作人员100余名，危险品运输车100余辆，主要支持第2，3，4，5，6类1项，8，9类危险品运输。同时企业建有先进的GPS监控中心，可实现对城市危险品运输车辆的实时监控，以及规范危险化学品仓储和运输的安全管理。为确定企业危险品运输风险等级，邀请了6位资深专家对危险品运输风险评价体系中各指标的重要程度、危害等级和概率等级进行独立性评分，分值区间为[0, 10]。

6名专家打分结果如表 2所示。

运用式(1)、(2)，对各项二级指标进行熵值计算，根据式(3)、(4)、(5)、对二级指标风险云数值进行计算，最后运用式(6)进行合成运算，得到综和风险云数字特征。

运用MATLAB编程得到综合风险云与标准风险云对比云图如图 2所示。

图 2 综合风险云 Fig. 2 Comprehensive risk cloud

图选项

根据式(6)和表 2、表 3得出二级指标的数字特征可知综合风险云的数字特征为(6.86，0.15，0.13)，(3.11，0.15，0.13)。最终在MATLAB中利用正向云发生器实现二维云成像。根据图 2及式(7)的计算结果，得出本企业危险品运输风险等级为Ⅱ级，即该企业在危险品道路运输过程中发生风险的概率较大、风险发生后危害较大、企业需采取相应管理措施规避风险。风险评价结果与该企业所获安全生产标准化达标建设等级一致。

(1) 基于二维云模型，综合考虑危险品运输评价指标的随机性和模糊性等特点，提出危险品运输风险等级评价方法，实现了定性、定量概念不确定转换。通过输出云图，得到直观、方便的风险等级评价结果。

(2) 使用熵权法全面获取了各评价指标权重系数，减少了系统风险评价的不确定性，使评价结果更准确可靠。

(3) 使用数学贴近度计算公式，真实的反映了危险品运输风险等级，进一提升评价结果可信度使得评价结果更具说服性。

(4) 基于熵权法二维云模型对危险品道路运输风险影响因素进行综合评价，评价结果与传统指数模型评价结果一致且与某危险品运输企业安全考评标准结果具有良好吻合性，验证了本评价模型的科学性和实际应用性。