专利摘要

本发明提供一种快速确定天然气液化工艺冷剂适宜配比的方法，通过对冷剂中n种组分xi进行编码，将自然变量xi变成规范变量zi，再根据实验方案利用流程模拟软件得到不同参数条件下的冷剂单位冷量建立起冷剂单位冷量与规范变量zi之间的高阶回归方程，通过规划求解的方法得到该高阶回归方程的各项回归系数，进一步计算得到适宜冷剂单位冷量条件下各规范变量zi的取值，从而可以得到与适宜冷剂单位冷量相对应的冷剂中各组分的百分比，即得到天然气液化工艺冷剂适宜配比。装置现场验证结果表明，按照本发明方法优化冷剂配比后，天然气液化工艺的实际单位能耗大大降低。

权利要求

1.一种快速确定天然气液化工艺冷剂适宜配比的方法，其特征在于，包括如下步骤：

第一步，以x1、x2、…、xn-1和xn表示冷剂中n种组分各占的百分比，每种组分的百分比均为非负的且各种组分的百分比总和为1；xi受到如下约束条件的限制，ai为各自然变量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，得到适宜的冷剂配比；

第三步中，基础数据组合方案如以下表格所示：

第三步中，所述流程模拟软件为ASPEN Plus、ASPEN HYSYS或CHEMICAL CAD；

第四步中，

其中，

yji、yj和yi为试验指标；bi、bji分别为高阶回归方程的一次项回归系数和交互项回归系数；

第五步中，所述的适宜冷剂单位冷量是根据制冷液化流程的理论单位能耗，并结合制冷装置的实际天然气加工能力、实际能量利用效率、实际保冷率和实际冷剂循环量计算得到的。

2.根据权利要求1所述的快速确定天然气液化工艺冷剂适宜配比的方法，其特征在于，ai＝0，则xi＝zi，i＝1,2,...,n-1,n。

3.根据权利要求1所述的快速确定天然气液化工艺冷剂适宜配比的方法，其特征在于，第五步中，通过规划求解的方法对高阶回归方程进行计算，得到适宜冷剂单位冷量条件下的各规范变量zi的取值。

4.根据权利要求1所述的快速确定天然气液化工艺冷剂适宜配比的方法，其特征在于，适宜冷剂单位冷量计算公式为：

式中，为适宜冷剂单位冷量，kJ/kg冷剂；为不同制冷液化流程的理论单位能耗，kJ/m3天然气；QLNG为实际天然气加工能力，m3天然气/h；η为实际能量利用效率；ζ为实际保冷率；QR为实际冷剂循环量，kg冷剂/h。

5.根据权利要求1所述的快速确定天然气液化工艺冷剂适宜配比的方法，其特征在于，第六步中，将适宜冷剂单位冷量条件下的各规范变量zi转换成自然变量xi，具体公式为：

x1＝[1-(a1+a2+...+an-1+an)]z1+a1

x2＝[1-(a1+a2+...+an-1+an)]z2+a2

...

xn-1＝[1-(a1+a2+...+an-1+an)]zn-1+an-1

xn＝[1-(a1+a2+...+an-1+an)]zn+an。

说明书

技术领域

本发明属于天然气液化领域，具体涉及一种快速确定天然气液化工艺冷剂适宜配比的方法。

背景技术

随着环境污染问题的日益加剧，天然气作为一种新型的优质清洁燃料，在能源、交通、化工等领域的应用越来越广泛。目前，天然气输送主要有两种方式：一种是管道输送，虽然输气量大，但输送成本高，沿线需加设加压站和调节站；另一种是液化输送，即净化后的天然气采用制冷工艺，在一个大气压下冷冻至-162℃变为液体，液化后体积仅为气态时1/600，可有效解决气体不方便运输的问题，为实现天然气的跨国贸易和回收偏、散、小的天然气资源提供了有利条件。

国内外天然气液化的方法较多，主要包括阶式制冷液化流程、多级单组分制冷液化流程、混合冷剂制冷液化流程、带预冷的混合冷剂制冷液化流程和膨胀机制冷液化流程等。冷剂作为这些制冷液化流程中冷能共同的媒介和载体，其配比大小对于天然气液化工艺实际运行效果的好坏具有非常重要的作用，若从天然气液化装置的冷能平衡角度来看，当冷剂提供的冷量小于天然气液化所需的冷量时，天然气液化率不足，产品收率降低；当冷剂提供的冷量大于天然气液化所需的冷量时，冷剂所携冷量富余，装置能耗增加。因此，适宜的天然气液化工艺冷剂配比不但可以显著降低装置能耗、提高产品收率，同时还能够大幅减少生产成本和提高企业竞争力。

然而，由于现场操作人员在装置实际运行过程中往往无法准确、可靠地通过调节冷剂配比来及时应对生产状况的变化，导致液化天然气吸收的冷量和冷剂提供的冷量难以形成理想匹配，故经常造成设备处理能力偏低或系统能耗增加。

发明内容

本发明的目的是提供一种快速确定天然气液化工艺冷剂适宜配比的方法，解决了现有技术中存在的液化天然气吸收的冷量和冷剂提供的冷量难以形成理想匹配的问题。

本发明是通过以下技术方案来实现：

一种快速确定天然气液化工艺冷剂适宜配比的方法，包括如下步骤：

第一步，以x₁、x₂、…、x_n-1和x_n表示冷剂中n种组分各占的百分比，每种组分的百分比均为非负的且各种组分的百分比总和为1；x_i受到如下约束条件的限制，

x_i≥a_i，i＝1,2,...,n-1,n

第二步，对x_i进行编码，将自然变量x_i变成规范变量z_i，使得x_i转换成0≤z_i≤1，其中，i＝1,2,...,n-1,n；则

第三步，根据冷剂中的组分数n确定相应的基础数据组合方案，并通过流程模拟软件得到不同规范变量z_i组合条件下的冷剂单位冷量

第四步，按照基础数据组合方案模拟得到的冷剂单位冷量 建立冷剂单位冷量 与规范变量z_i之间的高阶回归方程，并进一步求解得到该高阶回归方程的各项回归系数；

第五步，将第四步高阶回归方程中的冷剂单位冷量 替换为天然气液化所需的适宜冷剂单位冷量 并将各项回归系数代入，计算得到适宜冷剂单位冷量条件下的各规范变量z_i的取值；

第六步，将适宜冷剂单位冷量 条件下的各规范变量z_i转换成自然变量x_i，得到适宜的冷剂配比。

优选的，a_i＝0，则x_i＝z_i，i＝1,2,...,n-1,n。

优选的，第三步中，基础数据组合方案如以下表格所示：

优选的，第四步中，

其中，

优选的，第五步中，通过规划求解的方法对高阶回归方程进行计算，得到适宜冷剂单位冷量条件下的各规范变量z_i的取值。

优选的，第五步中，所述的适宜冷剂单位冷量是根据制冷液化流程的理论单位能耗，并结合制冷装置的实际天然气加工能力、实际能量利用效率、实际保冷率和实际冷剂循环量计算得到的。

进一步的，适宜冷剂单位冷量计算公式为：

式中， 为适宜冷剂单位冷量，kJ/kg冷剂； 为不同制冷液化流程的理论单位能耗，kJ/m³天然气；Q_LNG为实际天然气加工能力，m³天然气/h；η为实际能量利用效率；ζ为实际保冷率；Q_R为实际冷剂循环量，kg冷剂/h。

优选的，第六步中，将适宜冷剂单位冷量 条件下的各规范变量z_i转换成自然变量x_i，具体公式为：

x₁＝[1-(a₁+a₂+...+a_n-1+a_n)]z₁+a₁

x₂＝[1-(a₁+a₂+...+a_n-1+a_n)]z₂+a₂

...

x_n-1＝[1-(a₁+a₂+...+a_n-1+a_n)]z_n-1+a_n-1

x_n＝[1-(a₁+a₂+...+a_n-1+a_n)]z_n+a_n。

优选的，第三步中，所述流程模拟软件为ASPEN Plus、ASPEN HYSYS或CHEMICALCAD。

与现有技术相比，本发明具有以下有益的技术效果：

本发明通过对冷剂中n种组分x_i进行编码，将自然变量x_i变成规范变量z_i，再根据实验方案利用流程模拟软件得到不同参数条件下的冷剂单位冷量 建立起冷剂单位冷量 与规范变量z_i之间的高阶回归方程，通过规划求解的方法得到该高阶回归方程的各项回归系数，进一步计算得到适宜冷剂单位冷量 条件下各规范变量z_i的取值，从而可以得到与适宜冷剂单位冷量相对应的冷剂中各组分的百分比，即得到天然气液化工艺冷剂适宜配比。装置现场验证结果表明，按照本发明方法优化冷剂配比后，天然气液化工艺的实际单位能耗大大降低，与计算所得的理论单位能耗相差不大，同时装置的能量利用率达到了最大化，节能效果十分显著。

具体实施方式

下面结合具体的实施例对本发明做进一步的详细说明，所述是对本发明的解释而不是限定。

本发明快速确定天然气液化工艺冷剂适宜配比的方法，具体步骤如下：

第一步，根据不同制冷液化流程的理论单位能耗(详见表1)，结合各装置的实际天然气加工能力、实际能量利用效率、实际保冷率和实际冷剂循环量来计算适宜冷剂单位冷量。

表1不同制冷液化流程的理论单位能耗

式中， 为天然气液化所需的适宜冷剂单位冷量，kJ/kg冷剂； 为不同制冷液化流程的理论单位能耗，kJ/m³天然气；Q_LNG为实际天然气加工能力，m³天然气/h；η为实际能量利用效率；ζ为实际保冷率；Q_R为实际冷剂循环量，kg冷剂/h。

第二步，设定冷剂配比的约束条件。以x₁、x₂、…、x_n-1和x_n表示冷剂中n种组分各占的百分比，要求每种组分的比例必须是非负的且它们的总和必须是1，即

同时，除了上述约束条件外，若x_i还要受到其他约束条件的限制，则存在

x_i≥a_i，i＝1,2,...,n-1,n

a_i为各自然变量x_i对应的最小值(下限)。

第三步，对x_i(i＝1,2,...,n-1,n)进行编码，将自然变量x_i变成规范变量z_i，使得x_i转换成0≤z_i≤1，此时

即

对于由n种组分组成的天然气液化工艺冷剂来说，各组分所占的百分比分别为x₁、x₂、…、x_n-1、x_n，由如下约束条件x₁≥a₁、x₂≥a₂、...、x_n-1≥a_n-1、x_n≥a_n可知

a₁+a₂+...+a_n-1+a_n＜x₁+x₂+...+x_n-1+x_n＝1

因此，

x₁＝[1-(a₁+a₂+...+a_n-1+a_n)]z₁+a₁

x₂＝[1-(a₁+a₂+...+a_n-1+a_n)]z₂+a₂

...

x_n-1＝[1-(a₁+a₂+...+a_n-1+a_n)]z_n-1+a_n-1

x_n＝[1-(a₁+a₂+...+a_n-1+a_n)]z_n+a_n

需要说明的是，如果各组分x_i的下界都是0，则x_i＝z_i，这时自然变量x_i与规范变量z_i的值相等。

第四步，根据冷剂中的组分数n确定相应的基础数据组合方案，并通过流程模拟实验得到不同参数条件下的冷剂单位冷量，详见表2。所述流程模拟软件为ASPEN Plus、ASPENHYSYS或CHEMICAL CAD。

表2不同参数条件下的冷剂单位冷量

第五步，按照基础数据组合方案模拟得到冷剂单位冷量 建立冷剂单位冷量 与规范变量z_i之间的高阶回归方程，并进一步求解该高阶回归方程的各项回归系数。

其中，

y_ji、y_j和y_i为试验指标。

第六步，根据冷剂单位冷量 与规范变量z_i之间的高阶回归方程以及有关约束条件，通过规划求解的方法可预测适宜冷剂单位冷量 条件下各规范变量z_i的取值。

第七步，将适宜冷剂单位冷量 条件下的各规范变量z_i转换成自然变量x_i，就能够得到适宜的冷剂配比。

实施例一：

第一步，陕北某天然气液化装置采用多级混合冷剂制冷液化流程(理论单位能耗为1210kJ/m³天然气)，其实际天然气加工能力为21316m³天然气/h，实际能量利用效率为0.25，实际保冷率为0.94，实际冷剂循环量为73079kg冷剂/h。

天然气液化所需的适宜冷剂单位冷量 为

第二步，设定冷剂配比的约束条件。以x₁、x₂、x₃、x₄、x₅分别表示冷剂中氮气、甲烷、乙烯、丙烷和异戊烷这5种组分各占的百分比，要求每种组分的比例必须是非负的且它们的总和必须是1，即

同时，除了上述约束条件外，x₁、x₂、x₃、x₄、x₅均不受其他约束条件的限制，即

a₁＝a₂＝a₃＝a₄＝a₅＝0。

第三步，对x₁、x₂、x₃、x₄和x₅进行编码，将自然变量x₁、x₂、x₃、x₄、x₅变成规范变量z₁、z₂、z₃、z₄、z₅。鉴于a₁＝a₂＝a₃＝a₄＝a₅＝0，因此x₁＝z₁，x₂＝z₂，x₃＝z₃，x₄＝z₄，x₅＝z₅。

第四步，根据冷剂中的组分数确定相应的基础数据组合方案，并通过Aspen HYSYS流程模拟实验得到不同参数条件下的冷剂单位冷量，详见表3。

表3实施例一不同参数条件下的冷剂单位冷量

第五步，按照基础数据组合方案得到不同参数下的实验结果，建立冷剂单位冷量 与规范变量z₁、z₂、z₃、z₄、z₅之间的高阶回归方程，并进一步求解该高阶回归方程的各项回归系数。

将每号实验结果代入上述方程，可得

此时，

第六步，根据冷剂单位冷量 与规范变量z₁、z₂、z₃、z₄、z₅之间的高阶回归方程以及有关约束条件来预测适宜冷剂单位冷量 条件下各规范变量z₁、z₂、z₃、z₄、z₅的取值。

当 冷剂时，

1502＝732z₁+1168z₂+2096z₃+2758z₄+3241z₅+

136z₁z₂+592z₁z₃+484z₁z₄+894z₁z₅+112z₂z₃+92z₂z₄+94z₂z₅+116z₃z₄-2z₃z₅+434z₄z₅

通过规划求解的方法对以上高阶回归方程进行计算，最终得到规范变量z₁、z₂、z₃、z₄和z₅分别为

第七步，将适宜冷剂单位冷量 条件下的各规范变量z_i转换成自然变量x_i，就能够得到适宜的冷剂配比。

因为x₁＝z₁，x₂＝z₂，x₃＝z₃，x₄＝z₄，x₅＝z₅，所以

因此，当冷剂中氮气的百分比为49.51％(V/V)，甲烷的百分比为19.48％(V/V)，乙烯的百分比为13.91％(V/V)，丙烷的百分比为10.10％(V/V)且异戊烷的百分比为7.00％(V/V)时，液化天然气吸收的冷量可与冷剂提供的冷量形成理想的匹配。

装置现场验证结果表明，冷剂配比优化前该装置的实际单位能耗高达1982kJ/m³天然气，而冷剂配比优化后该装置的实际单位能耗仅为1475kJ/m³天然气，较优化前降低了约26％且与理论单位能耗1210kJ/m³天然气相差不大，节能效果十分显著。

一种快速确定天然气液化工艺冷剂适宜配比的方法专利购买费用说明

Q：办理专利转让的流程及所需资料

A：专利权人变更需要办理著录项目变更手续，有代理机构的，变更手续应当由代理机构办理。

1：专利变更应当使用专利局统一制作的“著录项目变更申报书”提出。

2：按规定缴纳著录项目变更手续费。

3：同时提交相关证明文件原件。

4：专利权转移的，变更后的专利权人委托新专利代理机构的，应当提交变更后的全体专利申请人签字或者盖章的委托书。

Q:专利著录项目变更费用如何缴交

A：（1）直接到国家知识产权局受理大厅收费窗口缴纳,（2）通过代办处缴纳，（3）通过邮局或者银行汇款,更多缴纳方式

Q：专利转让变更，多久能出结果

A：著录项目变更请求书递交后，一般1-2个月左右就会收到通知，国家知识产权局会下达《转让手续合格通知书》。