基于改进ISM-FRAM模型的天然气管道风险事故致因分析.pdf

能源化工Energy Chemical Industry第 44 卷第 3 期2023 年 6 月Vol.44 No.3Jun.,2023基于改进ISM-FRAM模型的天然气管道风险事故致因分析唐丽玲(杭州市城乡建设设计院股份有限公司，浙江杭州 310004)摘要：针对天然气管道风险事故的原因，提出了一种基于改进 ISM-FRAM 模型的天然气管道风险事故致因分析方法，介绍了解释结构模型(ISM)和功能共振分析模型(FRAM)，对天然气管道事故致因分析的结构层次进行划分，分别计算了网络层级重要度及功能共振模型中的结构重要度。研究结果表明，使用 ISM 可以对 FRAM 模型的功能网络进行合理划分，便于工作人员对事故致因进行合理的流程分析，探究天然气管道风险事故出现的表象原因以及根本原因，以便制定合理措施，有助于全面提高天然气管道运行的安全度。关键词：天然气管道；解释结构模型；功能共振分析模型；风险事故中图分类号：TE88 文献标志码：A 文章编号：2095-9834(2023)03-0069-07引文格式：唐丽玲.基于改进 ISM-FRAM 模型的天然气管道风险事故致因分析 J.能源化工,2023,44(3):69-75.Causal analysis of natural gas pipeline risk accidents based on improved ISM-FRAM modelTANG Liling(Hangzhou Urban&Rural Construction Design Institute Co.,Ltd.,Hangzhou 310004,Zhejiang,China)Abstract:A method for analyzing the causes of natural gas pipeline risk accidents based on the improved ISM-FRAM model is proposed.The interpretative structural model(ISM)and functional resonance analysis model(FRAM)are introduced,and the structural hierarchy of natural gas pipeline accident cause analysis is divided.The importance of the network hierarchy and the structural importance in the FRAM are calculated respectively.The research results indicate that the use of ISM can reasonably partition the functional network of the FRAM,facilitate staff to conduct reasonable process analysis of accident causes,explore the apparent causes and root causes of natural gas pipeline risk accidents,and develop reasonable measures to comprehensively improve the safety of natural gas pipeline operation.Key words:natural gas pipeline;ISM;FRAM;risk accident收稿日期：2022-07-07。作者简介：唐丽玲(1976)，女，贵州遵义人，本科，工程师，现主要从事电气与仪表设计工作。E-mail：。天然气长输管道系统由大量的子系统共同耦合构成，该系统的复杂度相对较高，且我国天然气管道网络较为庞大，在管道运行的过程中，不确定因素的数量也相对较多，各个作业流程之间也具有较强的耦合性，若出现安全风险问题，会造成不良的社会影响以及严重的经济损失1-3。对管道网络系统进行系统分析，并对风险事故的致因进行深入研究，可以制定出能全面提高管道系统安全性的有效预防措施4。有学者对天然气管道的风险问题进行了深入研究。Ariany 等5针对海底天然气管道的运行安全问题，基于管道所处环境的最坏假设，对海底管道运行过程中的风险因素进行了全面识别，研究发现海上船只作业以及海底土壤的沉降问题是影响海底管安全、节能与环保70能源化工2023 年 6 月道运行安全的重要因素。贺焕婷等6针对气田集输管道的安全问题，提出了一种半定量评价方法，在进行研究的过程中，将影响天然气管道安全的因素划分为 11 大类、59 小类，利用层次分析法确定了失效后果评价因子权重，以涩北气田的集输管道为例，对半定量评价方法进行了验证。都心爽等7针对工业管道的安全问题，采用三角直觉模糊故障树评价方法对某天然气净化厂管道进行了评价分析。王武昌等8将风险理论与脆弱性分析方法相互结合，进而提出了一种基于风险-脆弱度的天然气管道风险评价方法，在使用该种方法的过程中需要依靠风险偏好型效用函数，对后果的严重度进行计算，使用该类型的方法可以对管道运行过程中的风险进行全面识别。笔者主要从风险事故致因的角度出发，在对天然气管道运行风险因素进行初步识别的基础上，引入改进 ISM-FRAM 模型：首先对该种类型模型进行简单介绍，提出使用该模型对管道运行事故致因进行分析的应用过程；结合我国某天然气管道，对该模型应用的可行性以及先进性进行验证；通过对管道风险事故的致因进行全面分析，有助于提前制定合理措施以保障管道长期运行安全。1天然气管道风险事故致因分析模型1.1ISM模型ISM 模型是一种应用于系统工程的分析法，它可以把系统中各变量间的复杂关系转化为结构关系，并用有向结构图来表示9。1.2FRAM模型FRAM 模型是一种用于描述系统的振动特性的数学模型，该模型通过对系统的物理特性进行建模，可以预测系统的共振频率、振幅、模态等振动特性，从而为设计优化和故障排除提供重要的参考8-10。FRAM模型通常基于有限元方法或边界元方法建立，有限元方法是将结构分解成有限数量的小单元，通过数值求解有限元方程组得到结构的振动响应，边界元方法则是通过对结构的边界进行建模，求解边界元方程组得到结构的振动响应，对于复杂的结构系统，功能共振分析模型可以将结构分解成若干个子系统进行建模，然后将这些子系统的振动特性组合起来得到整体系统的振动特性，同时，该模型也可以考虑系统的材料特性、非线性特性、温度变化等因素对振动特性的影响11-12。在使用 FRAM 模型进行事故致因分析的过程中，主要可以分为以下 6 个步骤。1)确定系统和事件范围：首先需要确定分析的系统范围和事件范围，以便于分析人员理解分析的具体场景。2)描述系统功能结构：使用 FRAM 模型，绘制出系统的功能结构图，以便于理解系统中不同功能之间的相互关系和依赖关系。系统功能可根据六角形模型进行描述，六角形模型主要分为时间(T)、控制(C)、输出(O)、资源(R)、前提(P)和输入(I)等 6 个部分13。3)描述功能失效模式：针对每个系统功能，识别可能的失效模式，如功能缺失、错误的输入或输出、功能执行错误等。4)描述功能失效的后果：对于每种失效模式，识别可能的后果或影响，如人员伤亡、物质损失、环境破坏等。5)分析功能相互作用和共振：使用FRAM模型，分析不同功能之间的相互作用和共振效应，以识别可能导致事故的系统复杂性和非线性因素。6)确定安全控制措施：根据分析结果，确定安全控制措施，如人员培训、程序规范、设备检修等13。1.3FRAM模型的优化改进在使用 FRAM 模型进行天然气管道安全研究时，可以增加环境功能模块以更全面地考虑环境因素对管道安全的影响，通过增加环境功能模块，可以更全面、深入地研究天然气管道系统的安全性，并制定相应的管道安全管理策略和措施，从而提高管道系统的安全性14-15。具体而言，可以按照以下步骤来增加环境功能模块。1)确定环境因素：首先需要确定管道系统所处的环境因素，如温度条件、湿度条件等。2)定义环境功能模块：将环境因素作为输入变量，并定义环境功能模块，环境功能模块可以包括对环境因素的评估和对环境因素对管道系统的影响的建模。3)建立环境-管道系统关系模型：在环境功能模块中，需要建立环境因素与管道系统之间的关系模型，可以通过统计分析、数据挖掘和模型建立等技术手段来实现。71第 44 卷第 3 期4)进行系统安全分析：将环境功能模块与FRAM 模型的其他功能模块整合起来，进行系统安全分析，在分析过程中，需要考虑环境因素对管道系统安全的影响，并综合考虑各个功能模块之间的相互作用，以得出管道系统的整体安全评估结果16。在使用环境功能模块对天然气管道安全进行研究的过程中，研究内容可分为时间以及精度两个方面。在时间方面，时间方面的研究内容主要涉及管道系统与环境因素之间的关系随时间的变化，研究管道系统在不同时间段内与环境因素的交互作用如何影响管道系统的安全性，研究环境因素的变化如何对管道系统的安全性产生影响；在精度方面，研究内容主要涉及对环境因素的精确评估以及对管道系统与环境因素之间的相互作用的准确建模17-18。2改进ISM-FRAM模型2.1模型改进改进 ISM-FRAM 模型的应用具体可以分为 4个步骤。1)对天然气管道系统的功能进行分析，进而确定 FRAM 模型中的功能模块，将 FRAM 模型中的功能模块看作影响天然气管道系统安全的主要因素，各功能模块编号为 Si，Si S1,S2,Si,SN。2)通过使用 FRAM 模型的功能网络结构，进而建立 ISM 模型的关系有向图，使用 N 阶邻接矩阵 A=aij 对功能模块的相互关系进行表示，其中，N 表示为功能模块的数量。若上游某个功能模块 Si的输出可能会对下游功能模块 Sj产生影响，此时A=1；如果不产生影响，则 A=0。3)首先对邻接矩阵 A 进行布尔运算，就可以得到可达矩阵 M，M 矩阵计算公式见式(1)。M=(A+I)n=(A+I)n+1(A+I)n-1(1)式中：I表示单位矩阵。4)通过对可达矩阵的结构进行全面分析，就可依次得到可达集合 R(Si)、先行集合 Q(Si)以及共同集合 T(Si)。在 3 种类型的集合中，可达集合 R(Si)表示为可达矩阵 M 中要素 Si对应行中含有 1 的元素对应列的功能集合，先行集合 Q(Si)表示对应行的功能集合，共同集合 T(Si)表示可达集合 R(Si)以及先行集合 Q(Si)的交集19；其次，根据条件关系T(Si)=R(Si)，可以进行重复的层级抽取，最终得到改进 ISM-FRAM 网络模型的分层结构。2.2ISM-FRAM网络功能结构重要度计算针对改进 ISM-FRAM 模型网络，若该系统网络可以分为 N 层结构，各层网络使用 L 表示，L 1,2,3,N，该模型网络中各层级的权重为ML，则权重可由式(2)计算。ML=1/(N+1)-L1N(1/L)(2)将节点出入度(节点的出入度分为入度和出度两部分)引入到功能模块出入度中，可以对网络模型中各层结构的结构重要度进行计算，针对入度功能而言，其将会指向待计算的功能；针对待计算的功能，其将会指向出度功能20-21。同时，针对待计算的功能模块，其非输出端和输出端的结构重要度相同，网络模型中各层结构的结构重要度可根据式(3)计算。Ri=MLi(j MLj i+kMLi k)(3)式中：Ri为i功能模块的结构重要度；MLi 为待计算功能模块i的层级权重；MLji为指向i功能模块的j功能模块的层级权重；MLik为指向i功能模块的k功能模块的层级权重。2.3功能模块可变性定量方法在对功能模块输出结果的时间以及精度进行评分的过程中，可以根据可变性变化的时间将其划分为 4 种类型的情况，根据可变性变化的精度将其划分为 3 种类型的情况。功能模块 i 的时间变化可使用 STi,t进行评分，其中 t (1,2,3,4)分别表示为适时、过早、过迟以及未发生；功能模块的精度变化可使用 SPi,p进行评分，其中 p (1,2,3)分别表示为精确、可接受以及不精确，功能模块的变化评分见表 1。表1功能变化评分变化范围 变化程度 评分说明时间变化适时1功能正常，输出不存在偏差问题过早2输出存在一定的偏差过迟2未发生3无法完成正常输出精度变化精确1功能正常，输出不存在偏差问题可接受2输出准确率相对较低不精确3输出结果与正常结果的偏差较大唐丽玲.基于改进ISM-FRAM模型的天然气管道风险事故致因分析72能源化工2023 年 6 月研究中引入了改进区间值犹豫模糊熵的方法，进而对各个功能模块变化程度概率进行合理地计算，计算过程可以分为 4 个步骤。1)如 果 a=a-,a+=xa-x a+，其 中a-,a+R，则可以将看作为 1 个区间中的数值，如果 a-=a+，则将会 a逐渐退化为实数。天然气管道安全领域中的专家可根据区间数的方式对多个功能模块的时间变化范围进行可能性评分，其评分的范围为 0.00,1.00；在对功能模块精度变化范围进行可能性评分的过程中，其评分的范围也为0.00,1.00，区间数的波动范围不高于 0.20。2)对于专家而言，在对某一功能模块进行评分的过程中，可得到该功能模块时间变化范围以及精度变化范围的模糊元，模糊元可由式(4)和式(5)计算。h TZ(i,t)=hTZ-(K)(i,t),hTZ+(K)(i,t)(4)h PZ(i,p)=hPZ-(K)(i,p),hPZ+(K)(i,p)(5)式中：hTZ-(K)(i,t)为专家对功能模块i时间变化范围左侧的评分；hPZ-(K)(i,p)为专家对功能模块i精度变化范围左侧的评分；hTZ+(K)(i,t)为专家对功能模块i时间变化范围右侧的评分；hPZ+(K)(i,p)为专家对功能模i块精度变化范围右侧的评分。所有专家对功能模块的可变性评分可构成犹豫模糊集合hZ(i)。3)对犹豫模糊熵计算方法进行改进，对变化程度基本一致的打分值进行合理的融合，进而达到模糊化处理的效果。对于改进的函数模型，其将会呈现出单调递减的趋势，在将犹豫模糊元逐渐向 y 轴映射以后，可以更加完整地保留专家打分的相关信息。如果专家的数量为 Y，在对式(4)进行充分利用的基础上，可对某个功能模块时间变化区间数值的改进犹豫模糊熵 HTi,t 进行计算，计算公式见式(6)。HTi,t=1(2-1)YYK=1sin-14hTZ-(K)(i,t)+hTZ+(Y-K+1)(i,t)cos+114hTZ-(K)(i,t)+hTZ+(Y-K+1)(i,t)(6)在对式(5)进行充分利用的基础上，对某个功能模块精度变化区间的改进犹豫模糊熵进行计算，计算过程见式(7)。HPi,p=1(2-1)YYK=1sin-14hPZ-(K)(i,p)+hPZ+(Y-K+1)(i,t)cos+114hPZ-(K)(i,p)+hPZ+(Y-K+1)(i,t)(7)4)根据功能模块时间或者精度变化区间的改进犹豫模糊熵数值，可对时间变化概率占比和精度变化概率占比进行计算，计算公式见式(8)和式(9)。wTi,t=HTi,t/t HTi,t(8)wPi,p=HPi,p/t HPi,p(9)对于某一功能模块的可变性数值，可由式(10)进行计算。Ai=tSTi,t wTi,t +p SPi,p wPi,p (10)各个功能模块的耦合损失度主要是指各个功能模块在出现耦合现象时造成的损失程度，结合改进ISM-FRAM 组合模型的结构重要度以及功能可变性，可对耦合损失度进行计算，计算公式见式(11)。Bi=Ri tSTi,t wTi,t +p SPi,p wPi,p (11)3应用案例研究3.1改进ISM-FRAM模型分析天然气资源调度管理系统属于天然气长输管道的重要基础设施，主要由调度控制中心、应急中心、首站、末站、分输站以及管道沿线的阀室等部分组成。基于天然气资源调度管理系统的基本结构，并结合天然气管道风险事故案例，对基于改进 ISM-FRAM 模型的天然气管道风险事故致因分析中的功能模块进行完善。功能模块 F1F21 分别代表介质腐蚀性、管道所处环境、管道巡检情况、定期检测情况、管道周围人员活动频繁度、管道维护计划、介质调度计划、运行误操作情况、管道及设备运行情况把控、应急管理措施、管道信息获取便利性、管道信息利用率、管道及沿线设备使用时间、管道沿线阀门联锁情况、安全预警机制完善性、人员培训情况、信息通讯通畅性、管道及沿线设备运行参数设置、管理误操作情况、安全影响因素控制情况、管道保护情况等。对 FRAM 模型的可视化工具 FRAM Model Visualiser 进行充分利用，进而将FRAM 功能网络转化为有向图，并使用 ISM 模型对有向图的结构进行合理的层次划分，最终建立基于73第 44 卷第 3 期改进 ISM-FRAM 模型的天然气管道风险事故致因分析网络结构，网络结构见图 1。首先由天然气管道安全领域内的 5 名专家对天然气长输管道各个功能模块的变化程度进行打分，在对打分结果进行全面整理以后，将其代入到式(6)和式(7)中，最终得到各个功能模块的时间和精度变化区间值犹豫模糊熵和，然后对各个功能模块的时间以及精度变化概率占比和进行计算，最终得到计算结果见表 2。在结合各个功能模块时间变化和精度变化评分数值和的基础上，分别将其和代入到式(10)中，计算得到功能模块的可变性，计算结果见表 3。表2功能模块变化概率模块时间变化概率精度变化概率适时过早过迟未发生精确可接受不精确F10.680.190.070.030.480.360.13F20.340.320.290.060.290.560.14F30.520.320.120.040.420.420.19F40.520.280.120.050.470.460.08F50.730.140.070.020.520.420.07F60.690.230.080.020.520.410.09F70.760.130.110.010.610.380.08F80.750.140.090.010.620.340.05F90.640.220.140.010.640.330.03F100.760.140.090.010.630.320.04F110.690.110.210.010.580.350.08F120.610.520.110.060.520.340.11F130.720.110.130.010.470.410.11F140.660.240.060.020.520.390.12F150.690.230.070.020.520.410.11F160.450.340.110.050.360.480.16F170.670.030.260.010.520.410.08F180.810.080.070.020.560.380.09F190.690.030.260.010.520.410.08F200.810.080.070.030.540.380.09F210.720.120.110.040.520.360.14表3功能模块可变性功能模块可变性功能模块可变性功能模块可变性F12.969 3F82.674 2F152.919 3F23.616 2F92.741 2F163.354 2F33.242 9F100.674 6F172.845 2F43.153 1F110.815 2F182.774 1F52.843 2F123.089 2F192.846 9F62.941 6F132.941 8F202.764 1F72.691 2F142.946 8F212.965 2由表 3 可见：管道所处环境模块 F2、管道巡检情况模块 F3、定期检测情况模块 F4、管道信息利用率模块 F12、人员培训情况模块 F16 的可变性程度相对较大，对于管道所处环境模块 F2 而言，其变化情况难以直观察觉，输出精度的差异性相对较大，其功能可变性最高，管道巡检情况模块 F3、定期检测情况模块 F4 的可变性也相对较高，因此，需要进一步规范管道巡检制度，并定期开展管道检测工作，人员培训情况模块 F16 的可变性高，这说明通过提高人员技能水平和安全意识，可以有效提高管道运行的安全性。在天然气管道系统运行的过程中，众多的功能模块之间将会产生功能共振现象，功能共振可能会导致功能链接断裂，最终导致出现风险事故问题，通过对风险事故出现时功能模块的变化情况进行分析，可以制定合理的保障措施，确保管道长期处于安全状态。因此，将改进 ISM-FRAM 模型的结构分层代入式(2)得到各个网络层级的重要度数值，计算结果见表 4。将层级重要度的计算结果和各个功能模块节点的出入度分别代入式(3)得到功能模块结构重要度数值，计算结果见表 5。通过对图 1 进行分析可以发现，在事故根源层中，介质腐蚀性、管道所处环境以及人员培训情况图1改进ISM-FRAM模型事故分析网络分析唐丽玲.基于改进ISM-FRAM模型的天然气管道风险事故致因分析74能源化工2023 年 6 月表4网络层级重要度结构分层层级权重事故表象层0.038 4第2层0.045 6第3层0.051 6第4层0.059 7第5层0.071 2第6层0.088 7第7层0.127 1第8层0.192 3事故根源层0.363 8表5功能模块结构重要度功能模块结构重要度功能模块结构重要度功能模块结构重要度F10.063 1F80.005 6F150.026 3F20.234 1F90.004 8F160.185 2F30.149 6F100.038 4F170.026 4F40.145 8F110.012 6F180.018 6F50.083 4F120.035 2F190.016 8F60.041 6F130.134 2F200.010 6F70.027 3F140.026 7F210.018 2等因素将会对天然气管道风险事故的出现产生根源性影响，在事故网络模型的第8层中，管道巡检情况、定期检测情况属于保障天然气管道安全运行的基础条件，如果管道巡检不到位和没有对管道进行定期检测，则无法及时发现管道存在的风险问题，事故预防过程中的屏障将会失效，同时，第8层环境检查领域的相关信息在向第7层网络传递的过程中，功能链接将会断裂；在事故网络的第5层中，如果管道信息的利用率低，通过对管道运行数据进行分析无法及时发现管道存在的风险问题，或者管道出现了超压运行的状况，此时出现风险事故的概率必然会提升；在事故网络的第4层中，在管道信息利用率低的前提下，将无法对管道及相关设备的运行参数进行合理的设定，这也是风险事故出现的重要原因；在事故网络的第3层中，运行误操作和管理误操作等行为都将会对管道安全产生严重影响，在出现两种类型的行为以后，会导致事故网络第2层中管道及设备运行情况把控以及安全影响因素控制减弱，在事故表象层中，由于企业对管道及设备运行情况把控以及安全影响因素控制严重不足，最终会导致管道保护力度降低，管道风险事故出现的概率提升。通过进行全面的事故致因分析可以发现，通过确定各个功能模块功能共振的方式，可以得到事故网络中各个功能模块的共振影响，同时，还可以确定功能链接的失效情况。天然气管道事故致因网络结构中功能共振影响及功能链接失效表见表6。将表2和表4的计算结果代入到式(11)以后，对天然气管道事故致因网络结构的功能耦合损失度进行计算，计算结果见表7。表6功能共振影响及功能链接失效功能共振模块 功能共振影响因素失效功能链接F3F13(P)F3(O)-F13(P)F4F13(I)F4(O)-F13(I)F8F9(I)F8(O)-F9(I)F12F18(C)F12(O)-F18(C)F12F20(C)F12(O)-F20(C)F9F21(P)F9(O)-F21(P)表7功能耦合损失度功能模块耦合损失度功能模块耦合损失度功能模块耦合损失度F10.188 2F80.013 5F150.075 2F20.846 2F90.012 6F160.616 8F30.487 5F100.100 9F170.081 2F40.468 5F110.039 4F180.042 6F50.243 3F120.098 2F190.048 2F60.123 8F130.375 2F200.029 7F70.084 2F140.076 2F210.052 3由表 7 可见：F2、F3、F4、F13 以及 F16 等功能模块可能出现的耦合损失相对较大，说明这些功能模块出现耦合的前提下，天然气管道出现风险事故的概率相对较大。3.2有效功能屏障设置通过对功能模块的可变性以及各个功能模块的耦合损失度进行分析，可以筛选出天然气管道风险事故致因网络结构中需要着重关注功能模块 F2、F3、F4、F13 以及 F16，并需要根据这些功能模块制定合理的功能屏障措施，以此降低风险事故出现的概率。首先，在进行管道设计建设的过程中，需要着重关注管道沿线的环境情况，如果管道沿线出现地质灾害以及自然灾害的概率相对较大，则需要提前做好管道安全保障措施；其次，需要强化管道巡检工作，通过管道巡检的方式可以发现管道沿线存在的风险问题，并及时对其进行治理，定期对管道进行检测，及时发现管道存在的防腐层破损问题，对防腐层破损位置进行及时修复，及时掌握管道的腐蚀情况，对腐蚀较为严重的区域进行补强；最后，75第 44 卷第 3 期严格把控管道及沿线设备的使用时间，在管道及沿线设备达到使用寿命以后需要对其进行及时的更换，强化人员培训工作，人员因素属于影响管道安全的重要因素，强化人员培训可以使得人员的安全意识以及设备操作技能得到全面提升。4结论为了保障天然气管道的安全运行，通过引入基于改进 ISM-FRAM 模型的方式，从天然气管道风险事故致因的角度出发，对影响管道安全运行的因素进行了综合分析，对影响运行安全的主要因素进行了判断，并以此制定有效的功能屏障，降低风险出现的概率，得出以下结论。1)FRAM 模型可以用于对人员、技术以及组织等为了完成自身任务而进行的行为进行描述，在使用改进 FRAM 模型对天然气管道风险事故致因进行分析的过程中，可以直观地了解管道沿线环境变化与其他功能模块变化之间的耦合作用。2)使用 ISM 模型对改进 FRAM 模型的功能网络进行合理的划分以后，可以清晰地了解天然气管道风险事故出现的流程变化情况，分析风险事故出现的表象原因以及根本原因，其中，介质腐蚀性、管道所处环境以及人员培训情况等因素将会对天然气管道风险事故的出现产生根源性影响。3)通过使用改进 ISM-FRAM 模型对天然气管道风险事故致因进行案例分析发现，管道所处环境、管道巡检情况、定期检测情况、管道及沿线设备使用时间以及人员培训情况等功能模块对于风险事故的影响相对较大，因此，需要针对这些功能模块制定合理的措施，以此降低风险事故出现的概率。参考文献：1 程玉峰.保障中俄东线天然气管道长期安全运行的若干技术思考J.油气储运,2020,39(1):1-8.2 李玉星,杨昊,王武昌,等.基于属性区间识别模型的天然气管道事故后果评价技术研究J.安全与环境学报,2021,21(3):943-949.3 宋珂,武玮,马泽宇,等.天然气净化厂工业管道安全屏障的综合评价方法J.西北大学学报(自然科学版),2021,51(4):621-631.4 黄雪松,张佳,周兆明,等.输气管道第三方破坏数值分析及安全评价方法J.科学技术与工程,2020,20(1):320-326.5 ARIANY Z,SANTOSO B,SUHARTO.Risk assessment of subsea pipeline of natural gas in the cross area:case study in west java watersJ.IOP Conference Series Earth and Environmental Science,2020,557(1):1-10.6 贺焕婷,惠贤斌,吴东容,等.涩北气田集输管道风险评价方法修正J.油气储运,2020(8):892-897.7 都心爽,武玮,张金阳,等.基于三角直觉模糊故障树的工业管道安全评价方法与应用J.西北大学学报(自然科学版),2021,51(04):632-642.8 王武昌,章艺,李玉星,等.基于风险理论的天然气管网脆弱性分析方法J.油气储运,2022,41(1):42-47.9 汪应洛.系统工程理论方法与应用M.北京:高等教育出版社,1992.10 乔万冠,李新春,刘全龙.基于改进FRAM模型的煤矿重大事故致因分析J.煤矿安全,2019,50(2):249-252,256.11 尹德志,帅斌,黄文成,等.Tropos-FRAM法在道路客运事故分析中的应用J.中国安全科学学报,2020,30(8):151-157.12 左博睿,帅斌,黄文成.基于模糊推理的FRAM事故定量分析方法研究J.中国安全生产科学技术,2020,16(7):11-17.13 HIROSE T,SAWARAGI T.Development of FRAM Model based on structure of complex adaptive systems to visualize safety of socio-technical systemsJ.IFAC-PapersOnLine,2019,52(19):13-18.14 GEYER F,SAGEN H,CORNUELLE B,et al.Using a regional ocean model to understand the structure and variability of acoustic arrivals in Fram StraitJ.The Journal of the Acoustical Society of America,2020,147(2):1042-1053.15 RUTKOWSKA P,KRZYANOWSKI M.FRAM modelling of the transfer of control over aircraftJ.Scientific Journal of Silesian University of Technology-Series Transport,2018,101(8):159-166.16 HUANG W,SHUAI B,ZUO B,et al.A systematic railway dangerous goods transportation system risk analysis approach:the 24 modelJ.Journal of Loss Prevention in the Process Industries,2019,61(12):94-103.17 YOUSEFI A,HERNANDEZ M R,PEA V L.Systemic accident analysis models:a comparison study between AcciMap,FRAM,and STAMPJ.Process Safety Progress,2019,38(2):126-137.18 徐逸飞,帅斌,张玥,等.基于改进ISM-FRAM铁路运输事故致因分析J.中国安全生产科学技术,2021,17(5):161-168.19 王雪丽,程鑫,朱光华,等.基于ISM-ANP的堤防风险灰色评价模型及应用J.水电能源科学,2022,40(4):145-149.20 葛春景.基于ISM模型的国际航空物流中心发展影响因素分析J.科技和产业,2018,18(9):24-30.21 贾宝惠,史思杨,王玉鑫.基于改进DEMATEL-ISM模型的机轮刹车系统风险因素分析J.安全与环境学报,2021,21(2):506-512.唐丽玲.基于改进ISM-FRAM模型的天然气管道风险事故致因分析