本发明属于天然气液化领域，具体涉及一种快速确定天然气液化工艺冷剂适宜配比的方法。
背景技术：
：随着环境污染问题的日益加剧，天然气作为一种新型的优质清洁燃料，在能源、交通、化工等领域的应用越来越广泛。目前，天然气输送主要有两种方式：一种是管道输送，虽然输气量大，但输送成本高，沿线需加设加压站和调节站；另一种是液化输送，即净化后的天然气采用制冷工艺，在一个大气压下冷冻至-162℃变为液体，液化后体积仅为气态时1/600，可有效解决气体不方便运输的问题，为实现天然气的跨国贸易和回收偏、散、小的天然气资源提供了有利条件。国内外天然气液化的方法较多，主要包括阶式制冷液化流程、多级单组分制冷液化流程、混合冷剂制冷液化流程、带预冷的混合冷剂制冷液化流程和膨胀机制冷液化流程等。冷剂作为这些制冷液化流程中冷能共同的媒介和载体，其配比大小对于天然气液化工艺实际运行效果的好坏具有非常重要的作用，若从天然气液化装置的冷能平衡角度来看，当冷剂提供的冷量小于天然气液化所需的冷量时，天然气液化率不足，产品收率降低；当冷剂提供的冷量大于天然气液化所需的冷量时，冷剂所携冷量富余，装置能耗增加。因此，适宜的天然气液化工艺冷剂配比不但可以显著降低装置能耗、提高产品收率，同时还能够大幅减少生产成本和提高企业竞争力。然而，由于现场操作人员在装置实际运行过程中往往无法准确、可靠地通过调节冷剂配比来及时应对生产状况的变化，导致液化天然气吸收的冷量和冷剂提供的冷量难以形成理想匹配，故经常造成设备处理能力偏低或系统能耗增加。技术实现要素：本发明的目的是提供一种快速确定天然气液化工艺冷剂适宜配比的方法，解决了现有技术中存在的液化天然气吸收的冷量和冷剂提供的冷量难以形成理想匹配的问题。本发明是通过以下技术方案来实现：一种快速确定天然气液化工艺冷剂适宜配比的方法，包括如下步骤：第一步，以x1、x2、…、xn-1和xn表示冷剂中n种组分各占的百分比，每种组分的百分比均为非负的且各种组分的百分比总和为1；xi受到如下约束条件的限制，xi≥ai，i＝1,2,...,n-1,n第二步，对xi进行编码，将自然变量xi变成规范变量zi，使得xi转换成0≤zi≤1，其中，i＝1,2,...,n-1,n；则第三步，根据冷剂中的组分数n确定相应的基础数据组合方案，并通过流程模拟软件得到不同规范变量zi组合条件下的冷剂单位冷量第四步，按照基础数据组合方案模拟得到的冷剂单位冷量建立冷剂单位冷量与规范变量zi之间的高阶回归方程，并进一步求解得到该高阶回归方程的各项回归系数；第五步，将第四步高阶回归方程中的冷剂单位冷量替换为天然气液化所需的适宜冷剂单位冷量并将各项回归系数代入，计算得到适宜冷剂单位冷量条件下的各规范变量zi的取值；第六步，将适宜冷剂单位冷量条件下的各规范变量zi转换成自然变量xi，得到适宜的冷剂配比。优选的，ai＝0，则xi＝zi，i＝1,2,...,n-1,n。优选的，第三步中，基础数据组合方案如以下表格所示：优选的，第四步中，其中，优选的，第五步中，通过规划求解的方法对高阶回归方程进行计算，得到适宜冷剂单位冷量条件下的各规范变量zi的取值。优选的，第五步中，所述的适宜冷剂单位冷量是根据制冷液化流程的理论单位能耗，并结合制冷装置的实际天然气加工能力、实际能量利用效率、实际保冷率和实际冷剂循环量计算得到的。进一步的，适宜冷剂单位冷量计算公式为：式中，为适宜冷剂单位冷量，kj/kg冷剂；为不同制冷液化流程的理论单位能耗，kj/m3天然气；qlng为实际天然气加工能力，m3天然气/h；η为实际能量利用效率；ζ为实际保冷率；qr为实际冷剂循环量，kg冷剂/h。优选的，第六步中，将适宜冷剂单位冷量条件下的各规范变量zi转换成自然变量xi，具体公式为：x1＝[1-(a1+a2+...+an-1+an)]z1+a1x2＝[1-(a1+a2+...+an-1+an)]z2+a2...xn-1＝[1-(a1+a2+...+an-1+an)]zn-1+an-1xn＝[1-(a1+a2+...+an-1+an)]zn+an。优选的，第三步中，所述流程模拟软件为aspenplus、aspenhysys或chemicalcad。与现有技术相比，本发明具有以下有益的技术效果：本发明通过对冷剂中n种组分xi进行编码，将自然变量xi变成规范变量zi，再根据实验方案利用流程模拟软件得到不同参数条件下的冷剂单位冷量建立起冷剂单位冷量与规范变量zi之间的高阶回归方程，通过规划求解的方法得到该高阶回归方程的各项回归系数，进一步计算得到适宜冷剂单位冷量条件下各规范变量zi的取值，从而可以得到与适宜冷剂单位冷量相对应的冷剂中各组分的百分比，即得到天然气液化工艺冷剂适宜配比。装置现场验证结果表明，按照本发明方法优化冷剂配比后，天然气液化工艺的实际单位能耗大大降低，与计算所得的理论单位能耗相差不大，同时装置的能量利用率达到了最大化，节能效果十分显著。具体实施方式下面结合具体的实施例对本发明做进一步的详细说明，所述是对本发明的解释而不是限定。本发明快速确定天然气液化工艺冷剂适宜配比的方法，具体步骤如下：第一步，根据不同制冷液化流程的理论单位能耗(详见表1)，结合各装置的实际天然气加工能力、实际能量利用效率、实际保冷率和实际冷剂循环量来计算适宜冷剂单位冷量。表1不同制冷液化流程的理论单位能耗制冷液化流程理论单位能耗(kj/m3天然气)阶式制冷液化流程1152单级混合冷剂制冷液化流程1440丙烷预冷单级混合冷剂制冷液化流程1325多级混合冷剂制冷液化流程1210单级膨胀制冷液化流程2304丙烷预冷单级膨胀制冷液化流程1958两级膨胀制冷液化流程1958式中，为天然气液化所需的适宜冷剂单位冷量，kj/kg冷剂；为不同制冷液化流程的理论单位能耗，kj/m3天然气；qlng为实际天然气加工能力，m3天然气/h；η为实际能量利用效率；ζ为实际保冷率；qr为实际冷剂循环量，kg冷剂/h。第二步，设定冷剂配比的约束条件。以x1、x2、…、xn-1和xn表示冷剂中n种组分各占的百分比，要求每种组分的比例必须是非负的且它们的总和必须是1，即同时，除了上述约束条件外，若xi还要受到其他约束条件的限制，则存在xi≥ai，i＝1,2,...,n-1,nai为各自然变量xi对应的最小值(下限)。第三步，对xi(i＝1,2,...,n-1,n)进行编码，将自然变量xi变成规范变量zi，使得xi转换成0≤zi≤1，此时即对于由n种组分组成的天然气液化工艺冷剂来说，各组分所占的百分比分别为x1、x2、…、xn-1、xn，由如下约束条件x1≥a1、x2≥a2、...、xn-1≥an-1、xn≥an可知a1+a2+...+an-1+an＜x1+x2+...+xn-1+xn＝1因此，x1＝[1-(a1+a2+...+an-1+an)]z1+a1x2＝[1-(a1+a2+...+an-1+an)]z2+a2...xn-1＝[1-(a1+a2+...+an-1+an)]zn-1+an-1xn＝[1-(a1+a2+...+an-1+an)]zn+an需要说明的是，如果各组分xi的下界都是0，则xi＝zi，这时自然变量xi与规范变量zi的值相等。第四步，根据冷剂中的组分数n确定相应的基础数据组合方案，并通过流程模拟实验得到不同参数条件下的冷剂单位冷量，详见表2。所述流程模拟软件为aspenplus、aspenhysys或chemicalcad。表2不同参数条件下的冷剂单位冷量第五步，按照基础数据组合方案模拟得到冷剂单位冷量建立冷剂单位冷量与规范变量zi之间的高阶回归方程，并进一步求解该高阶回归方程的各项回归系数。其中，yji、yj和yi为试验指标。第六步，根据冷剂单位冷量与规范变量zi之间的高阶回归方程以及有关约束条件，通过规划求解的方法可预测适宜冷剂单位冷量条件下各规范变量zi的取值。第七步，将适宜冷剂单位冷量条件下的各规范变量zi转换成自然变量xi，就能够得到适宜的冷剂配比。实施例一：第一步，陕北某天然气液化装置采用多级混合冷剂制冷液化流程(理论单位能耗为1210kj/m3天然气)，其实际天然气加工能力为21316m3天然气/h，实际能量利用效率为0.25，实际保冷率为0.94，实际冷剂循环量为73079kg冷剂/h。天然气液化所需的适宜冷剂单位冷量为第二步，设定冷剂配比的约束条件。以x1、x2、x3、x4、x5分别表示冷剂中氮气、甲烷、乙烯、丙烷和异戊烷这5种组分各占的百分比，要求每种组分的比例必须是非负的且它们的总和必须是1，即同时，除了上述约束条件外，x1、x2、x3、x4、x5均不受其他约束条件的限制，即a1＝a2＝a3＝a4＝a5＝0。第三步，对x1、x2、x3、x4和x5进行编码，将自然变量x1、x2、x3、x4、x5变成规范变量z1、z2、z3、z4、z5。鉴于a1＝a2＝a3＝a4＝a5＝0，因此x1＝z1，x2＝z2，x3＝z3，x4＝z4，x5＝z5。第四步，根据冷剂中的组分数确定相应的基础数据组合方案，并通过aspenhysys流程模拟实验得到不同参数条件下的冷剂单位冷量，详见表3。表3实施例一不同参数条件下的冷剂单位冷量第五步，按照基础数据组合方案得到不同参数下的实验结果，建立冷剂单位冷量与规范变量z1、z2、z3、z4、z5之间的高阶回归方程，并进一步求解该高阶回归方程的各项回归系数。将每号实验结果代入上述方程，可得此时，第六步，根据冷剂单位冷量与规范变量z1、z2、z3、z4、z5之间的高阶回归方程以及有关约束条件来预测适宜冷剂单位冷量条件下各规范变量z1、z2、z3、z4、z5的取值。当冷剂时，1502＝732z1+1168z2+2096z3+2758z4+3241z5+136z1z2+592z1z3+484z1z4+894z1z5+112z2z3+92z2z4+94z2z5+116z3z4-2z3z5+434z4z5通过规划求解的方法对以上高阶回归方程进行计算，最终得到规范变量z1、z2、z3、z4和z5分别为第七步，将适宜冷剂单位冷量条件下的各规范变量zi转换成自然变量xi，就能够得到适宜的冷剂配比。因为x1＝z1，x2＝z2，x3＝z3，x4＝z4，x5＝z5，所以因此，当冷剂中氮气的百分比为49.51％(v/v)，甲烷的百分比为19.48％(v/v)，乙烯的百分比为13.91％(v/v)，丙烷的百分比为10.10％(v/v)且异戊烷的百分比为7.00％(v/v)时，液化天然气吸收的冷量可与冷剂提供的冷量形成理想的匹配。装置现场验证结果表明，冷剂配比优化前该装置的实际单位能耗高达1982kj/m3天然气，而冷剂配比优化后该装置的实际单位能耗仅为1475kj/m3天然气，较优化前降低了约26％且与理论单位能耗1210kj/m3天然气相差不大，节能效果十分显著。当前第1页1 2 3