石油与天然气化工  2017, Vol. 46 Issue (2): 53-58
基于层次分析法的冷箱冻堵调整
杨远 1, 唐文洁 2, 张春光 2, 肖春生 2, 马玉华 2     
1. 油气资源与勘探技术教育部重点实验室(长江大学);
2. 新疆油田公司石西油田作业区
摘要:为了研究冷箱冻堵原因,量化分析问题,找出最佳调整方案, 结合某厂液化装置冷箱冻堵实例,收集冻堵前与紧急处理后数据,运用层次分析法结合专家经验建立数学模型。利用该模型,仿真模拟了在当前成本与风险下的最佳处理方案和应急处理方案, 为冷箱冻堵调整提供了量化依据。研究结果表明,当前情况下可以均衡降量,同时关小LNG线、冷剂线截流阀进行应急处理,但不能根治问题。由此推测冷箱盘管有异物堵塞,需停产并拆卸冷箱对盘管、过滤器等附件进行检修。停产检修拆解情况与仿真模拟吻合, 彻底解决了问题,说明该方法可行,可在其他装置上分析和处理类似问题。但层次分析法只能核算出各方案层权重值,为方案决策做出参考,暂不能得出可调区间及调整程度,下步仍需改进。
关键词冷箱    冻堵    层次分析法    模拟仿真    应急处理    
Adjustment of frozen blockage of cold box based on the analytic hierarchy process
Yang Yuan1 , Tang Wenjie2 , Zhang Chunguang2 , Xiao Chunsheng2 , Ma Yuhua2     
1. Key Laboratory of Exploration Technologies for Oil and Gas Resources (Yangtze University), Wuhan, Hubei, China;
2. Oil Field Operation Area, Xinjiang Oil Field Co., Karamay, Xinjiang, China
Abstract: In order to study the reason of frozen blockage, analyze problem quantitatively, and find out the best adjustment scheme, combined with the example of frozen blockage of cold box in one liquefying plant, the data before frozen blockage and after the emergency treatment are collected, then the mathematical model is established by using analytic hierarchy process (AHP) combined with expert experience. By using the model, the best treatment scheme and the emergency treatment plan under current cost and risk are simulated, and quantitative basis for adjustment of cold box frozen blockage can be provided. The simulating result shows that the quantity can be lowered equilibrium currently, at the same time, the intercept valves of LNG line and cooling line can be shut off for emergency treatment, but the problem can not be completely solved. So it is suspected that there is a blockage in cold box coil, shutdown must be done and the cold box must be disassembled to overhaul coil, filters and other accessories. The overhaul situation is in accordance with the simulation result. The problem has been solved completely, and it shows that the method is feasible and can be used to analyze and treat the same problem in other devices. But the AHP can be only used for accounting layer weight value of the program, and to make a reference for scheme decision-making. The adjustable range and degree can not be obtained at present, and the simulation model needs to be improved next step.

层次分析法(Analytic Hierarchy Process,简称AHP)是将与决策有关的元素分解成目标、准则、方案等层次,并在此基础之上进行定性和定量分析的决策方法[1-2]。该法最早运用在电力分配问题上,是一种简洁、实用、系统性的分析方法。

冷箱作为一种紧凑、易维护的板翅式换热器,广泛运用于LNG行业[3]。但因盘管材质、设计精度、环境等问题,导致冷箱不能进行DCS自动控制。手动控制情况下还会发生因操作员不同而引发换热效率低,导致产品不达标,甚至盘管堵塞、联锁停机等情况。据统计,以苯及其衍生物为主的环状芳香烃凝固导致的冷箱天然气管线冻堵是LNG行业中的常见现象。当前的研究方向主要集中在优化冷箱设计和完善自控仪表方面[4]。对冻堵情况的分析也仅仅是现场实例叙述,缺乏调整分析。本研究基于层次分析法,对使用混合冷剂换热的冷剂线冻堵实例进行了分析,并优选了调整方案。

1 工艺简介

某厂液化装置采用多级单组分制冷液化工艺,该工艺总体上为传统阶式制冷工艺,最后一级改进为混合冷剂制冷,3种冷剂分别为丙烯、乙烯和甲烷。甲烷制冷循环的制冷剂为以甲烷为主的配方冷剂,包含了氮气、甲烷和乙烯3种组分,其中氮气体积分数为5%~10%,乙烯15%~30%,其余为甲烷。氮气由工厂氮气站供给,同时将已处理的原料气用于甲烷补偿。乙烯、丙烯从各自罐车供应。该装置具体采用了3个制冷循环[5-7]。第一级丙烯制冷循环为天然气、制冷剂乙烯和制冷剂甲烷提供冷量;第二级乙烯制冷循环为天然气和制冷剂甲烷提供冷量;第三级甲烷混合冷剂制冷剂循环为天然气及自身提供冷量。通过6个换热器和一个板翅式换热器(冷箱)冷却,天然气的温度逐渐降低,直至液化。而最后一级冷箱的调控以及混合冷剂的配比是影响LNG产品的关键[8]。冷箱的不稳定会导致冷箱冻堵、LNG产品温度不达标、能耗过大等多方面不良后果。本文抛开前端复杂工艺,着重对冷箱运行情况进行分析并记录调整过程。

2 问题描述

正常生产时,过冷至-155 ℃的甲烷混合冷剂进入冷箱,再经J-T阀节流至表压140 kPa、-158.9 ℃后返回冷箱,为LNG、甲烷冷剂自身提供冷量,实现甲烷制冷系统的冷剂循环。由于生产需要,LNG产量由冷箱处产品线控制阀提量(300×104 m3/d提至350×104 m3/d),同时加大冷量供给,保持换热平衡。提量后,由于冷剂线供给增大导致差压报警。紧急降量,冷箱工艺流程图见图 1,各测点及相关运行数据时间轴见表 1

图 1     冷箱工艺流程图 Figure 1     Process flow chart of cold box

表 1    各测点及相关运行数据时间轴 Table 1    Time axis of each measurement point and the related running data

图 2~图 5可以看出,紧急降量后冷剂线回温响应快,天然气线无明显响应。天然气线堵塞程度比冷剂线严重。分析原因后怀疑是甲烷水露点不合格引起天然气线形成水合物而导致冰堵,同时怀疑天然气线、冷剂线滤网破裂导致机械杂质进入引起盘管堵塞[9]

图 2     冷剂线温度堆积折线图 Figure 2     Temperature stacking line chart of refrigerant line

图 3     天然气线温度堆积折线图 Figure 3     Temperature stacking line chart of natural gas line

图 4     差压堆积折线图 Figure 4     Pressure difference stacking line chart

图 5     LNG产量堆积折线图 Figure 5     LNG yield stacking line diagram

针对以上分析,现需要一套科学可行的调整方案。在此结合操作经验和专家建议,选取5种调整方案,运用层次分析法针对该厂工况进行科学分析。

3 模拟仿真
3.1 计算原理

层次分析法数学模拟首先建立有支配关系的分组型层次结构。然后通过矩阵构造,综合比较各因素间相互影响关系,计算其在系统中的权重,最终确定各因素的相对重要性,达到分析或决策的目的[10]。其计算步骤如下:

(1) 确定问题,构建物理模型并收集数据。

(2) 依次构建最高层(目标层)、中间层(准则层)以及最低层(方案层)层次结构模型。

(3) 两两比较打分,确定下层对上层的分数。

准则层中的各准则在目标衡量中所占的比重并不一定相同,在决策者的心目中,它们各占有一定的比例。引用数字1~9及其倒数作为标度来定义判断矩阵A=(aij)nxn (见表 2)。

表 2    判断矩阵标度定义 Table 2    Scale definition of judgment matrix

(4) 层次合成计算与一致性检验。

1) 计算一致性指标CI

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

式中:λmax为判断矩阵的最大特征值。

2) 查找一致性指标RI(见表 3)。

表 3    平均随机一致性指标 Table 3    Consistency index of average random

3) 计算一致性比例CR

$CR = \frac{{CI}}{{RI}}$ (2)

CR<0.10时,认为通过了一次性检验,否则应作适当修正。

(5) 计算权重向量W。层次分析法有4种计算方法求权重:算术平均法、几何平均法、特征向量法。在此,本文选用几何平均法(方根法)。

几何平均法(方根法):

${W_i} = \frac{{{{\left( {\prod _{j = 1}^n{a_{ij}}} \right)}^{\frac{1}{n}}}}}{{\sum\nolimits_{i = 1}^n {{{\left( {\prod _{j = 1}^n{a_{ij}}} \right)}^{\frac{1}{n}}}} }},i = 1,2, \cdots ,n$ (3)

计算步骤:

1) A的元素按行相乘得一新向量;

2) 将新向量的每个分量开n次方;

3) 将所得向量归一化即为权重向量。

3.2 运算过程

基于现有工况、处理手段、考虑细节,以冷箱冻堵调整为目的,指导下步工作。首先根据调整原则确定建模中间层(准则层)。然后结合操作经验和专家建议,选取5种调整方案作为最低层(方案层)构建层次结构模型(见图 6)。

图 6     冷箱冻堵调整结构模型图 Figure 6     Structure model of frozen blockage adjustment of cold box

冷箱冻堵调整各判断矩阵如表 4所示。

表 4    冷箱冻堵调整各判断矩阵(4a~4g) Table 4    Judgment matrices of frozen blockage adjustment of cold box(4a-4g)

将不同方案各准则要素的权重矩阵(表 4a)Wi与各准则要素的相对权重矩阵(见表 4b~表 4g)Wi相乘,得到各方案层要素权重并排序(见表 5)。经计算还能得到中间层要素权重并排序(见表 6)。

表 5    方案层权重排序 Table 5    Ranking of scheme layer weights

表 6    中间层权重排序 Table 6    Ranking of the middle layer weights

3.3 结果与验证

经以上模拟分析,由表 5得到冷箱冻堵调整排序为: C5>C3>C4>C1>C2。由此可知,针对某厂冷箱冻堵事件,停产检修,拆卸冷箱对盘管、过滤器等附件检修(C5)是解决问题的最佳方案。通过模拟怀疑冷箱有过滤器损坏,盘管异物堵塞等问题。若执行此项方案必须停机,有损耗放空和其他成本损失。而均衡降量,同时关小LNG线、冷剂线截流阀(C3)能进行应急处理,但不能彻底解决问题。调节冷剂组分配比(C4)也能进行应急处理,但须专业技术人员指导,并实时取样化验,较为繁琐。单纯降产量,关小LNG线截流阀(C1)与单纯降冷量,关小冷剂线截流阀(C2)两种方案权重最小,不推荐。

通过模拟分析,某厂进行了均衡降量,同时关小LNG线、冷剂线截流阀应急处理,冻堵情况有所缓解。但后续无论如何调节冷剂组分配比和冷量换热调整都无法根治。在停工检修时拆解冷箱天然气与冷剂管线过滤器时发现,过滤器150 μm纱网均不同程度损坏,冷剂线过滤器损坏尤其突出。检修结果与模拟相吻合。再次开工复产后,提量稳产至400×104 m3/d。从4月开工至7月运行数据(见表 7)可以看出,冻堵问题已根治。

表 7    冷箱复产提量后稳产数据 Table 7    Stable output data after reproduction of cold box

4 结论

(1) 冷箱冻堵问题复杂。原因多,可能性复杂。最常见的有苯及苯衍生物为主的环状芳香烃凝固导致的冷箱天然气管线冻堵,也有因为特殊温度、压力条件导致天然气水合物的形成。冷剂线冻堵较为少见。而管道及设备安装时都会进行吹扫干燥。异物、水分引起的堵塞更是稀少。但也不排除因过滤器损坏而导致的管道上游填料、粉尘进入冷箱盘管导致阻塞的情况。也有因盘管中有异物而形成凝结核而变相导致的复合性水合物冻堵。在不能随意停工检修的情况下只能依靠经验进行判断。

(2) 层次分析法能将冷箱冻堵复杂问题分解化,并根据当前情况和专家经验分析进行模拟。最后通过量化,算出最佳调整办法。根据现场检修情况与调整后数据得出,层次分析法模拟结果与现场检修拆解情况吻合,此方法可行。下步可运用在其他装置的问题分析与决策中。

(3) 层次分析法暂时只能核算出各方案层权重值,为方案决策作出参考,而不能提出可调区间及调整程度,该法仍需进一步改进。

参考文献
[1]
邓雪, 李家铭, 曾浩健, 等. 层次分析法权重计算方法分析及其应用研究[J]. 数学的实践与认识, 2012, 42(7): 93-100.
[2]
王衍东, 顾洁, 胡斌. 基于AHP的实用中长期电力负荷预测综合模型[J]. 华东电力, 2005, 33(1): 28-31.
[3]
王世杰, 薛维海, 孙开友, 等. 天然气深冷冷箱存在问题分析[J]. 油田节能, 2003(2): 33-34.
[4]
王坤, 徐风雨, 李红艳, 等. 小型MRC天然气液化装置中板翅式换热器动态特性仿真研究[J]. 低温工程, 2007(3): 44-49.
[5]
王治红, 吴明鸥, 王小强, 等. 富含CO2天然气低温分离防冻堵工艺研究[J]. 天然气与石油, 2012, 30(4): 26-29.
[6]
钟志良, 汪宏伟, 白宇恒, 等. LNG工厂冷箱堵塞问题解决方案研究[J]. 天然气与石油, 2013, 31(5): 26-29.
[7]
刘盛鹏. 广安LNG装置冷箱积液处理[J]. 石油与天然气化工, 2013, 42(2): 154-155.
[8]
朱玉琴, 缪斌. 并联管组流动特性的数值模拟[J]. 西安石油大学学报(自然科学版), 2011, 26(3): 94-96.
[9]
孙恒, 余霆, 舒丹. LNG液化装置冷箱降温过程的动态模拟[J]. 低温与超导, 2009, 37(7): 14-15.
[10]
洪志国, 李焱, 范植华, 等. 层次分析法中高阶平均随机一致性指标(RI)的计算[J]. 计算机工程与应用, 2002, 38(12): 45-47. DOI:10.3321/j.issn:1002-8331.2002.12.017