扎根理论是质性研究方法的一种，其基本思想是围绕研究问题，以经验事实为依据对原始资料中的概念进行归纳，并进一步提炼相关概念和分析其相互关系，自下而上建立一套实质理论 [13] 。扎根理论的研究方法包括文献分析、典型案例分析、实地调研访谈等。扎根步骤包括：开放式编码、主轴式编码、选择式编码和理论饱和度检验，在对新的资料编码后不再出现新范畴时认为理论达到饱和，最终构建以核心范畴为核心、主范畴为维度、初始范畴为指标的理论框架。文本利用典型案例分析的方式构建危化品仓储风险因素体系，收集到的事故报告数量多、描述详尽，可以达到理论饱和。

本文选取了2004年~2024年间的108起典型危化品仓储事故，收集详细事故调查报告资料，报告中须包含事故发生时间地点、伤亡情况、事故经过、详细事故原因分析(直接、间接)等。本文是以危化品企业管理视角探讨其仓储风险因素，因此仅选取事故报告中对于企业自身相关原因的文本描述，政府和其他相关机构管理因素、地形气候等环境因素等暂不做探讨。

本文案例数据来源主要有：中华人民共和国应急管理部、安全管理网、国家危险化学品安全公共服务互联网平台(NRCC)、中国化学品安全协会官方网站、以及各地应急管理部通报、网页新闻报道等。

表1 是部分事故调查报告相关数据节选。

为了确保扎根理论所构建的风险因素理论框架科学完整，需要对编码结果进行饱和度检验 [14] 。对108份事故报告按照事故类型进行分类，并按一定比例选取编码案例与饱和度检验案例，最终选取88份案例报告作为编码资料，其余20份案例报告作为理论饱和度检验资料。事故调查报告按事故类型分类如 表2 所示。

扎根理论编码过程包含三级编码，首先进行一级编码，即开放式编码(Open Coding)，需要从原始资料中对原始语句进行类属识别并对概念进行命名。以按比例选取后的88份危化品仓储事故调查报告为编码资料一一提取出其中描写事故原因的原始文本，分别将每个案例事故原因概念化，接着将近似的概念进行合并，最终提炼出52个初始概念，即初始范畴。部分编码过程如 表3 所示。

接着进行第二级编码，即主轴式编码(Axial Coding)，需要找出初始范畴风险因素之间的联系，对其进行概念类属的归类，在归类时要充分考虑研究对象的特点和原始资料语境。

对52个初始范畴因素反复比较，发现它们并不是独立的，分析它们之间的联系，根据逻辑关系、从属关系等归纳成更高级别的21个相对独立的事故影响因素，即21个主范畴因素。

最后进行第三级编码，即选择式编码(Selective Coding)，需要对已经归类好的概念类属进行进一步地核心归类，核心类属应当是一个比较宽泛的理论范围，就有集中性，将最大量的概念类属囊括在内，在核心类属归纳完成之后，理论框架也已初步形成。

对21个主范畴因素进一步提炼和分类归纳，得出最终的5个核心范畴，分别为设施设备管理、人员安全行为管理、企业日常安全管理、作业现场管理和应急救援管理。

通过所选案例中预留的20份调查报告作为新增原始资料对得出的编码结果进行饱和度检验。按照上述步骤，继续对这20份资料文本一一进行概念化和逐级编码，发现并没有出现新的范畴与关系，也就是说这20份案例对应的编码结果从属于前文88份案例编码结果。因此，本文对危化品仓储事故影响因素的归纳已达到理论饱和。

最后构建包含5个核心范畴，21个主范畴和52个初始范畴的危化品仓储风险因素理论体系如 图1 所示。

依据前文所得扎根理论模型三级编码结果，构建危化品仓储风险因素递阶层次结构模型。以“危化品仓储风险因素”作为顶层目标层，并包含5个一级因素指标，即五大核心范畴因素，和21个二级因素指标，即主范畴因素。根据从属关系构建的危化品仓储风险因素递阶层次结构模型如 图2 所示。

以问卷调查的形式，成立专家组对各层因素重要度进行打分。为了尽可能地减少不同类型因素之间相互比较的难度，采用因素间两两比较的相对尺度以提高准确度。即通过对不同因素两两比较其重要性打分，打分准则选择1~9分标度法。记 $a_{ij}$为因素i与因素j重要度之比，n为单一准则下因素个数，即判断矩阵阶数。则单一准则下构造的判断矩阵为：

$A={\left(a_{ij}\right)}_{n\times n},\left(i=1,2,\cdots ,n;j=1,2\cdots ,n\right)$(1)

本次接受问卷调查的对象共100名，包括高校相关研究领域教授与H危化品公司操作人员、安全员、安全负责人等，共发放100份问卷，收回100份问卷，其中有效问卷94份，问卷有效率达94%。共得到6个判断矩阵如下：

$A=\left[\begin{array}{ccccc}1& 1/2& 1/3& 1/2& 2\\ 2& 1& 1/2& 1& 3\\ 3& 2& 1& 2& 4\\ 2& 1& 1/2& 1& 2\\ 1/2& 1/3& 1/4& 1/2& 1\end{array}\right]$ $A_1=\left[\begin{array}{ccc}1& 3& 2\\ 1/3& 1& 1\\ 1/2& 1& 1\end{array}\right]$ $A_2=\left[\begin{array}{ccccc}1& 5& 1& 3& 7\\ 1/5& 1& 1/4& 1& 3\\ 1& 4& 1& 3& 7\\ 1/3& 1& 1/3& 1& 3\\ 1/7& 1/3& 1/7& 1/3& 1\end{array}\right]$

$A_3=\left[\begin{array}{cccccc}1& 5& 2& 1& 3& 5\\ 1/5& 1& 1/2& 1/3& 1& 2\\ 1/2& 2& 1& 1/2& 2& 3\\ 1& 3& 2& 1& 4& 6\\ 1/3& 1& 1/2& 1/4& 1& 1\\ 1/5& 1/2& 1/3& 1/6& 1& 1\end{array}\right]$ $A_4=\left[\begin{array}{cccc}1& 2& 4& 3\\ 1/2& 1& 3& 3\\ 1/4& 1/3& 1& 1\\ 1/3& 1/3& 1& 1\end{array}\right]$ $A_5=\left[\begin{array}{ccc}1& 1/2& 3\\ 2& 1& 4\\ 1/3& 1/4& 1\end{array}\right]$

因素权重是指每个因素对于上一级因素的相对重要性，本文利用算术平均法计算各级因素权重，先对判断矩阵按列归一化处理，公式如下：

$a_{ij}=\frac{a_{ij}}{{\displaystyle {\sum }_{i=1}^n{a}_{ij}}},\left(i=1,2,\cdots ,n;j=1,2,\cdots ,n\right)$(2)

再计算每一行的平均值，即可得到每个因素的权重值，公式如下：

$w_i=\frac{1}{n}{\displaystyle \underset{j=1}{\overset{n}{\sum }}{a}_{ij}},\left(i=1,2,\cdots ,n;j=1,2,\cdots ,n\right)$(3)

计算得到一个一级指标判断矩阵权重向量，和五个二级指标判断矩阵权重向量，汇总结果见 表4 。

为了检验打分结果是否科学可靠，判断打分结果是否存在相互矛盾，需要进行一致性检验。需引入最大特征值定义： $n$阶判断矩阵 $A$与特征向量的乘积等于矩阵乘以某个值 ${\lambda }_{\max }$，即 $AW=A{\lambda }_{\max }$，该值即为最大特征值。最大特征值的近似计算公式如下：

${\lambda }_{\max }={\displaystyle \underset{i=1}{\overset{n}{\sum }}\frac{{\left(Aw\right)}_i}{n{w}_i}}$(4)

矩阵一致性检验系数CR计算公式为：

$\text{CR}=\frac{\text{CI}}{\text{RI}}$(5)

式中，CI为一致性指标，计算公式如下：

$\text{CI}=\frac{{\lambda }_{\max }-n}{n-1}$(6)

表5. 随机一致性指标标准值

若CR < 0.1，则一致性检验通过，判断矩阵构造合理，所求权重值也视为合理。

计算得到一致性检验结果汇总表格如 表6 所示。

由 表6 可知CR值均小于0.1，则一致性检验通过， 表4 即为最终权重值。

模糊综合评价法(FCE)是以模糊数学为理论基础的综合评价方法，用于解决受到多种因素制约的非确定性现实问题，在层次分析法确定因素权重的基础上，依据模糊数学中的“最大隶属度”原则对难以量化的问题和模糊概念进行定量评价 [17] 。模糊综合评价法既依托于人的专业知识与主观经验判断，同时又运用科学严谨的数学理论，二者结合使得评价结果清晰、系统性强，应用非常广泛。

H危化品公司是一家省重点国有物流企业，该公司整合优势物流资源，拥有完善的仓储设施。港口目前共拥有各类危化品储罐21个，包括6个万方储罐，储存能力达10万立方。主要储存化学品包括甲醇、乙醇、醋酸、液体二氧化碳、液态硫磺等，其中，甲醇、乙醇、醋酸属重大危险源。本文以H危化品公司仓储风险水平评价为实证研究目标，具有典型性，能够为众多危化品企业提高风险识别能力、降低仓储事故发生概率提供理论支持。

危化品仓储风险水平评价需要考虑的风险因素较多，若要对众多元素分配权重，则每个风险因素的权值都很小，为了解决这些难题，可将因素划分为两个层级，即采用二级模糊综合评价的方法进行本研究综合评价。为了获取H公司仓储风险现状相关资料，研究者深入现场调研，并成立专家小组开展访谈，并发放问卷获取数据开展H公司仓储风险等级评价。具体评价步骤如下文。

基于扎根理论编码结果和层次分析构建的递阶层次结构模型，将5个一级评价因素(核心范畴概念)细分为包含21个二级评价因素(主范畴概念)的因素层，这些二级因素能够直观地描述风险环节，具有可操作性。

第一级因素集如下：

U = {设施设备管理U₁，人员安全行为管理U₂，企业日常安全管理U₃，作业现场管理U₄，应急救援管理U₅}。

第二级因素集如下：

设施设备管理U₁= {储罐故障u₁₁，罐区设备故障u₁₂，安防设施u₁₃}；

人员安全行为管理U₂= {安全素质u₂₁，人员资质u₂₂，人员违规作业u₂₃，人员失误u₂₄，其他人为因素u₂₅}；

企业日常安全管理U₃= {风险防治工作u₃₁，规范化管理u₃₂，安全培训u₃₃，安全管理制度u₃₄，企业违规建设生产u₃₅，历史问题整改u₃₆}；

作业现场管理U₄= {作业准备工作u₄₁，现场安全监管u₄₂，储罐安全处置u₄₃，危化品处置u₄₄}；

应急救援管理U₅= {应急响应u₅₁，应急准备工作u₅₂，救援工作u₅₃}。

根据危化品仓储风险因素的特点，本研究选择李克特五级量表法来衡量各因素的风险大小，将危化品仓储风险水平划分为5个等级：

$V=\left\{v_1,v_2,v_3,v_4,v_5\right\}=\left\{低，较低，中等，较高，高\right\}$

本研究共邀请了15位来自H危化品公司港口罐区的相关从业人员成立专家小组，包括1位安全总监、1位安全部长、1位副部长、4位安全员以及8位罐区操作人员。以问卷调查的形式请专家们依据对H公司危化品仓储风险管理现状的了解如实对每个因素的风险等级进行判断，对五个评价等级进行勾选。整理结果后得到每个二级风险因素对应每个评价等级的选择人数，对其归一化处理后得到总体模糊判断矩阵R，并根据因素集划分得到单因素评价矩阵R₁~R₅。

$R_1=\left[\begin{array}{ccccc}0.47& 0.46& 0.07& 0.00& 0.00\\ 0.47& 0.33& 0.20& 0.00& 0.00\\ 0.33& 0.33& 0.27& 0.07& 0.00\end{array}\right]$ $R_2\text{=}\left[\begin{array}{ccccc}0.33& 0.53& 0.07& 0.07& 0.00\\ 0.53& 0.40& 0.07& 0.00& 0.00\\ 0.47& 0.26& 0.20& 0.07& 0.00\\ 0.40& 0.53& 0.07& 0.00& 0.00\\ 0.60& 0.40& 0.00& 0.00& 0.00\end{array}\right]$ $R_3\text{=}\left[\begin{array}{ccccc}0.67& 0.20& 0.13& 0.00& 0.00\\ 0.73& 0.27& 0.00& 0.00& 0.00\\ 0.33& 0.60& 0.07& 0.00& 0.00\\ 0.53& 0.40& 0.07& 0.00& 0.00\\ 0.73& 0.20& 0.07& 0.00& 0.00\\ 0.20& 0.40& 0.27& 0.13& 0.00\end{array}\right]$ $R_4\text{=}\left[\begin{array}{ccccc}0.53& 0.40& 0.07& 0.00& 0.00\\ 0.26& 0.60& 0.07& 0.07& 0.00\\ 0.33& 0.47& 0.13& 0.07& 0.00\\ 0.27& 0.33& 0.27& 0.13& 0.00\end{array}\right]$ $R_5\text{=}\left[\begin{array}{ccccc}0.53& 0.40& 0.07& 0.00& 0.00\\ 0.60& 0.33& 0.07& 0.00& 0.00\\ 0.33& 0.60& 0.07& 0.00& 0.00\end{array}\right]$

利用层次分析得到的风险因素权重与评价矩阵作合成运算，进行单因素模糊综合评价，公式如下：

$B_i=W_i\cdot R_i\left(i=1,2,\cdots ,5\right)$(7)

计算得到各二级风险因素评判向量，综合单因素评价结果得到总评价矩阵如下：

$R=\left(\begin{array}{c}{B}_1\\ B_2\\ B_3\\ B_4\\ B_5\end{array}\right)\text{=}\left[\begin{array}{ccccc}0.436& 0.401& 0.146& 0.017& 0.000\\ 0.421& 0.415& 0.113& 0.051& 0.000\\ 0.556& 0.343& 0.094& 0.007& 0.000\\ 0.394& 0.460& 0.101& 0.045& 0.000\\ 0.544& 0.386& 0.070& 0.000& 0.000\end{array}\right]$

前文已求出一级风险因素权重为 $W=\left(0.123,0.217,0.377,0.200,0.082\right)$，计算得到综合评价结果：

$B=W\cdot R\text{=}\left(\begin{array}{ccccc}0.478& 0.393& 0.104& 0.025& 0.000\end{array}\right)$

依据最大隶属度原则，一级风险因素对于评价等级“低”的隶属度最大，为0.478。由此可知，H公司危化品仓储风险等级为“低”。这表示该公司危化品仓储环节处于较安全的状态，暂时不需要采取风险控制措施。