天然气液化混合冷剂配方优化研究.pdf

第 4 4卷第 3期 石 油 与天 然 气 化工 CHEMI CAL E NGI NE ERI NG OF OI L 8 L GAS 6 5 天然气液化 混合冷剂 配方优化研 究 孟毅 明 王 东军 陈 博 1 ．中国石 油大学 北京 2 ．中国石油天然 气管道局 天津设计 院 摘 要 混合冷 剂制 冷是 目前天 然 气 低 温 液化 的 重要 手 段 ， 由 于操 作 工 况 复 杂 、 冷 荆介 质 多相 共 存 ， 制 冷剂 的选择 及 配 方优 化 已成 为天 然 气混合 冷 剂液 化技 术 面临 的主 要 问题 。本 文 采 用均 匀设 计 实 验方法对混合冷剂配方进行研 究， 首先依据不同制冷剂组分在不 同温 区制冷的原理 ， 初 步选定基本冷剂 组 分 ； 然后 通过 HYS YS模 拟 不 同组成 混合 冷剂 的制 冷效 果 ， 通 过 实验 结果 分析 混合 冷剂 中各组 分 的作 用； 在保持操作压力及处理量不变的条件下， 以冷剂的最小循环量为 目标参数进行优化分析， 最终形成 混 合 冷剂优 化 配方 ， 从 而避免 了采用 复杂 原理探 讨 和试 算的传 统 方法 。 关键词 天 然气液 化 混合 冷 剂 流程模 拟 配 方 均 匀设 计 优化 中图 分类 号 TE 6 2 4 ． 9 文 献标 志码 A D o I 1 0 ． 3 9 6 o ／ j ． i s s n ． 1 0 0 7 3 4 2 6 ． 2 0 1 5 ． O 3 ． O 1 4 S t u d y o n mi x e d r e f r i g e r a n t f o r mu l a t i o n o p t i mi z a t i o n f o r n a t u r a l g a s l i q u e f a c t i o n Me n g Yi mi n g ，Wa n g Do n g j a n ，Ch e n B o 1 ．C h i n a U n i v e r s i t y o f Pe t r o l e u m B e i j i n g ， B e i j i n g 1 0 2 2 4 9 ，C h i n a 2 ．C h i n a Pe t r o l e u m Pi p e l i n e Bu r e a u Ti a n j i n De s i g n I n s t i t u t e ，Ti a n j i n 1 2 0 1 0 7，Ch i n a Ab s t r a c t M i x e d r e f r i ge r a nt m e t ho d i s a n i m po r t a n t m e a ns t o l i q ue f y n a t u r a l ga s a t l o w t e mpe r a t u r e ． Due t o t he c ompl e x o pe r a t i ng c o nd i t i o n a nd c o o l i ng m e di u m m uhi pha s e c o e x i s t e n c e， s e l e c t i o n a n d f o r mu l a t i o n o p t i mi z a t i o n o f mi x e d r e f r i g e r a n t b e c o me t h e ma j o r c h a l l e n g e f o r t h e n a t u r a l g a s l i q u e f a c t i on t e c hn i qu e ．Thi s a r t i c l e us e d uni f o r m de s i g n e x pe r i m e n t s t o s t ud y m i xe d r e f r i ge r a nt f o r mul a t i on ． Th a t i s，p r e l i m i n a r i l y s e l e c t t he ba s i c c om p o ne n t s by t he pr i n c i p l e t ha t d i f f e r e n t r e f r i g e r a nt ha s a n e f f e c t i n i t s p a r t i c ul a r r e f r i ge r a t i ng r a ng e；t he n s i mul a t e t he r e f r i ge r a t i ng e f f e c t o f d i f f e r e nt c o m p o n e nt s b y us i ng p r o c e s s s i m u l a t i o n s o f t wa r e ASPEN HYSYS，a nd a na l y z e t he f unc t i on o f e a c h c o m po n e nt wi t h l e s s e xpe r i m e n t s；und e r t he c o ns t a nt pr e s s ur e a n d t r e a t m e nt c a p a c i t y，t a ki n g t he mi ni m u m c i r c u l a t i o n qu a n t i t y o f m i x e d r e f r i ge r a n t a s t a r ge t pa r a me t e r t o ma ke op t i ma l a na l y s i s，f i na l l y o bt a i n t h e o p t i mi z a t i o n f o r mu l a t i o n o f mi x e d r e f r i g e r a n t ，wh i c h a v o i d a d o p t i n g t h e c o n v e n t i o n a l a p p r o a c h w i t h c o m p l e x p r i nc i p l e a n a l ys i s a n d t r i a l c a l c ul a t i on ． Ke y wo r d sna t ur a l ga s l i q u e f a c t i on，m i x e d r e f r i ge r a nt ，p r oc e s s s i mul a t i on，f o r m ul a t i o n，un i f or m de s i gn，o pt i m i z a t i o n 随着 环境 污染 问题 的 日益加 剧 ， 天然 气 作 为一 种 优 质清 洁燃 料 ， 在 能源 、 交 通 、 化 工 等领 域 的应 用 越 来 越 广泛 。液化 天然 气便 于运 输 和使用 ， 因此 ， 开展天 然 气液化技术研究对扩大天然气的应用范围具有十分重 要 的意 义 。 目前 ， 天然 气液 化 主要通 过制 冷 实现 ， 大致 分为 级 联式、 膨胀式和混合冷剂 3种方式 。与前两者相 比， 混 合 冷剂 技术 具有 流 程 简单 、 机 组设 备 少 、 投 资少 、 能耗 低 等 特点 。 目前 世 界 上 8 O 的基 本 负 荷 型 天 然 气 液 化装置都采用了混合冷剂液化流程n ] 。 由于操 作工 况 复杂 、 冷剂介 质 多相共 存 ， 天然 气混 合 冷剂 液化 技术 面 临 的主要挑 战是 制冷 剂 的组成 确定 作者 简介 孟毅 明 1 9 9 1 一 ， 男 ， 天津 市人 ， 中国石油大学 北京 化学工程学 院 ， 硕士研究生 。E - ma i l me n g y m7 1 6 3 ． c o m 学兔兔 w w w .x u e t u t u .c o m 6 6 孟毅明 等 天然气液 化混 合冷剂配方优化研究 及配方优化 ] 。鉴于混合冷剂制冷过程中各单纯制冷 工质问复杂的热力学交互作用 ， 采用高效的均匀设计 确定 初步 试验 方案 ， 通 过 HYS YS流 程模 拟 软 件 对 混 合制冷过程进行模拟研究 ， 归纳预测混合冷剂优化组 成 并 进一 步模 拟验 证 。 1 混合冷剂天 然气 液化流程 混合冷剂天然气液化流程包括制冷循环和天然气 液化 ， 具体 工 艺流 程见 图 1 [ 3 _ 。 1 ． 1 制 冷循 环 低压混合冷剂气体首先进入一级压缩机 C 1 0 1压 缩 ， 经出口冷却器冷却降温， 使部分高沸点组分凝结后 进入 气液 分 离器 V l 。分 离器 底 部 液 相经 预 冷 换 热 器 换 热并通 过 节 流 阀 VL V一 1节 流 后 ， 与来 自主 换 热 器 L NG 一 2的逆 流低 压 气体汇 合 ， 作 为预 冷 换 热器 L NG 一 1 的 冷流 ； 分离 出的气 相 混 合 冷 剂 ， 经 高 压 压 缩 机 C l O 2 压缩 ， 由换 热器E 一 2 降 温后进入 分 离器V一 2 。 分 离器 图1 混合冷剂天然气液化流程图 Figu r e 1 Pr o c es s fl o w d i a gr a m o f na t ur al g as l i q u ef ac t i o n b y mix ed r e f r iger a nt 底 部 液相经 预 冷 换 热器 后 进 入 主换 热 器 L NG - 2中冷 却 ， 经节流阀 VL V一 3 节流后与来 自过冷换热器 L NG 一 3的低压 冷 剂 汇 合 ， 并 返 回 主 换 热 器 L NG 一 2提 供 冷 量 ； 气 相 先 后 进 入 主 换 热 器 L NG 一 2和 过 冷 换 热 器 L NG一 3完成 冷 凝 和 过 冷 ， 混 合 冷 剂 经 节 流 阀 VL V一 4 节流降温后为天然气和冷剂 自身提供冷量。 1 ． 2 天 然气 液化 净 化后 的天 然气 先后 经 预 冷换 热 器 L NG 一 1 、 主换 热器 L NG 一 2逐 步冷 却 ， 然后 进入 气液 分离 器 V一 4脱 除 气体 中少量重组分。重组分返 回换热器 E 一 3为冷剂预 冷提供部分冷量 ， 天然气中的轻组分则继续进入过冷 换 热 器 L NG 一 3冷 凝 、 液 化 ， 液 化 天 然 气 经 节 流 阀 VL V一 5节流降压分离 出不凝气后注人 L NG储罐 。 2 天然气液化混合冷 剂配方研究 混 合 冷剂 组 成的确 定 目前 主要采 用优 化 模型计 算 法 和原 理优 化 法 ] 。计 算 法是 通 过 选 取 实 验 影 响 因素 编写 优化 方程 ， 利 用热 力学 公式 优化 冷 剂组 成 ， 该 方法需要对制冷系统建立复杂的热力 学模 型， 优化计 算 困难 ； 原理优化法是结合流程模拟和原理判断 ， 通过 提高 或 降低混 合 物 中某 一 种组 分 的含量 来 最终 确定 混 合物 组 成 ， 该 方法 需 要大 量 的试 算 且 不 易 寻 求 最 优 结 果 。鉴 于混合 冷 剂制 冷过 程是 复 杂 的热 力 学效 应交 互 作用 的结果 ， 本文采用原理初选结合均匀实验设计 的 混合冷剂配方优化方法， 通过 HYS YS软件对混合制 冷 过程 进行 流程 模 拟 ， 利 用 少 量具 有 代 表 性 的实 验 结 果 建立 考察 指标 与 冷剂 组 成 问定 量 关 系 ， 进 而 预 测 混 合 冷剂 的优 化 组成 ， 并对 预测 结 果进 行模 拟验 证 。 2 ． 1混合 冷剂 初选 混合 冷剂 制冷 主要 是通 过蒸 发 吸热 及复 杂 的换 热 过程实现的， 这与阶式制冷的原理相类似 ] ， 即高沸点 冷剂 组分 主要 在 高温 区发 挥 制 冷 效果 ， 低 沸 点 冷 剂 组 分 主要在 低温 区发 挥制 冷效 果 。不 同 的 是 ， 混 合 冷 剂 中各 组 分可 共 同在单 压缩 循 环下逐 级 蒸发 、 冷 凝 ， 实现 学兔兔 w w w .x u e t u t u .c o m 第 4 4卷第 3期 石油与天然气化工 CHEMI CAL ENGI NE ERI NG OF OI L 8 L GAS 不 同沸 点组 分 在不 同温 区制 冷 。 各个组 分 都有 其 特 定 的 制冷 温 区范 围 ， 制 冷 温 区 位 于 一1 8 0 2 0 0℃的多种 冷 剂组分 在 0 ． 3 MP a及 4 ． 0 MP a压力 下 的温焓 特 性如 图 2所示 。 鲁 虫 避一 1 一’ 温度 ／ ℃ 图2 混合冷剂组分在O ． 3 MP a 、4 ． 0 MP a 下的 等压温焓特性曲线 Figur e 2 T- H c h ar ac t e ds t i c c ur v es o f diff er e n t r e f r i g er an t s u nd er t h e c o ns t an t pr es s ur e o f 0．3 MPa a nd 4． 0 MPa 由 图 2可 知 ， 各 制 冷 剂 组 分 均存 在 特定 的制 冷 温 区 ， 且 随着 冷剂 组分 的相 对分 子质 量 的增加 ， 冷 剂组 分 的沸 点逐 渐升 高 ， 相 应 的制 冷温 区也 向较 高温 区移 动 ， 气 化 潜热 增大 ， 液 化 所 需 的 功耗 较 大 。相邻 组 分 间 的 制冷 温 区有部 分 重叠 现象 。相 比而 言 ， 在 同一 制 冷 温 区， 为达到相同制冷效果， 相对分子质量高的组分功耗 小但 制 冷效果 不 足 ； 相对 分 子 质 量 低 的组 分 制冷 效 果 较好 ， 但 功耗 大 。 因此 ， 需 均 衡 选 择 适 宜 的 冷剂 组 分 ， 温区过多重叠的组分制冷效果是相似的， 在冷剂筛选 中应 尽 量避 免温 区重 叠 。 为 满 足 常规 天 然 气液 化 要 求 ， 需 选择 一 1 8 0 ～2 0 0 。 C的温 区范 围 。对 于 N 及 C ～C 制 冷体 系 ， 首 先 选 择 N 、 C H 组分 作 为低 温 区制冷 组分 。鉴 于 C H 与 CH 及 c 。 H 的制冷温 区重叠较小 ， 因此保 留 C 。 H 组 分 ， 剔 除 同温 区的 C H。 组分 。鉴 于 C 。 H。 及 C H 。 能 够连 续 覆 盖 中高 温 区 ， 因此 ， 剔 除 中 间温 区 重 叠 的 C H 。 组 分 。综 上 所 述 ， 最 终 选 择 N 、 CH 、 C 。 H 、 C 。 H 、 c H 5种组分作为混合冷剂的基本组成 。 2 ． 2混 合冷 剂配 方均 匀 实验研 究 均匀 实验 设计 的实 验 点 具 有代 表 性 强 、 试 验次 数 少 的特 点 。通 过 数学 拟 合 ， 可 建 立 考 察 指 标 与 冷剂 组 成的定量关系， 进而实现定量分析 ] 。 本 文首 先确 定实 验 影 响 因素 及 优 化指 标 ， 再结 合 流程模 拟 对不 同组 成 的 配方 进 行 均 匀 实 验 ， 并 对 实 验 结果 进行 优化 分析 及 验 证 ， 最 终 确 定适 宜 的混 合 冷 剂 配方 。 流程 模拟 参考 图 1中的混 合 冷 剂 液 化 流 程 ， 固定 操作压力 、 处理气量等条件 ， 考察不同混合冷剂配 比对 制 冷效果 的影 响 ， 并 以 最低 冷 剂 循 环 量 作 为考 察 不 同 配 方混 合冷 剂效 果 的指标 ， 进 行均 匀实 验 。 根 据某 工程 实 际 ， 选 取 原料 气 的组 成 及操 作 条件 分 别 如表 1和 表 2 所 示 。 表 1 原料天 然气 组成 皇 堡 呈 呈 旦 竺 塑 坌 旦 旦 y ／ o ． 0 2 o ． 6 4 o ． 0 0 1 9 9 6 ． 0 5 2 ． 4 6 0 ． 5 2 0 ． 1 0 ． 1 0 ． 0 5 0 ． 0 2 0 ． 0 3 表 2 混合冷剂天然气液化工艺操作参数 Table 2 Ope r at in g pa r amet er o f n a t u r a l ga s l i q u ef a ct ion p r oc e s s by m i x ed r e f r i g er an t 现利 用初 步选 择 的 5种 冷剂 组分进 行 配方评 选 实 验 。参考 AP C I 。 提 出 的 混 合 制 冷 剂 构 成 的 一 般 原 则， 确定其摩尔分数及变化范围如表 3所示 。 表 3 实验 因素 及 变 化 范 围 Ta bl e 3 Exp er i men t al f a c t or s an d t h ei r v ar i a t i on r a n ge s 在表 3规定的冷剂组成范围内， 各因素取 1 O个水 平 ， 并通过表 4所示 的均匀设计表 U 进行混合冷剂 的制 冷实验 l 1 ] 。考虑 各组 分摩 尔含 量涉 及归 一化 的 问题 ， 本文 采用 超 低温 区液化 的 He 作 为 稀 释剂 ， 保 证 混合冷剂摩尔组成加和为一， 待求得优化配 比再移除 稀释 组分 。后 续 实 验 分 析 结 果 表 明 ， He的加 入 量 对 优化 指标 的影 响不 显著 ， 因此 ， 可 作 为冷剂 配方 实验 的 稀释 组分 。 通过对组成如表 4所示 的混 合冷剂进行模 拟实 验 ， 调节 冷剂循 环 量 以满足 天然 气液 化所 需冷 量 ， 同时 学兔兔 w w w .x u e t u t u .c o m 6 8 孟 毅明 等 天然气液化混合冷剂 配方 优化研究 表 4均 匀 实验 设 计 表 Table 4 Unif o r m e xp er i men t a l de si g n t a bl e 保证 HYS YS工艺模型运行正常。将不同组成混合冷 剂 实验 得到 的最 低 冷剂 循 环 量 分 别 记 录 ， 实 验 结 果 如 图 3所 示 。 1 2 3 ’4 5 6 7 8 9 1 0 实验序号 图3 不同混合冷剂配方时天然气液化混合冷剂最小循环量 Figu r e 3 Mi n i mu m c i r c ulat i on qu a nt i t y o f mix ed r e f r iger an t f or n a t u r al ga s l i q ue f a c t i on b y diffe r e n t r e f r ige r a n t f or mulat i o n 由图 3 可 知 ， 不 同配方 混合 冷剂 的制 冷效 率 不 同 ， 其中 1 O 配方混合冷剂的循环量最低 ， 为 2 1 9 7 k mo l ／ h 。为便 于 优化 分析 ， 采 用 式 1 的计 算 值 作 为 优化 指 标 ， 均匀实验结果及对应的优化指标如表 5 所示 。 表 5 液化 天然气混合冷剂配方均 匀设计 实验结果 Tab l e 5 Un i f or m d es i g n ex per i men t al r es ul t s by mi x ed r e f r i g er a n t f or mul a t i o n y M 1 L 式中， y 为 优 化 指 标 ； M 为 转 换 因 子 ， 取 1 0 0 0 0 0 k mo l ／ h ； L 为循环 量 ， k mo l ／ h 。 2 ． 3实 验结 果优 化 对表 5中的实 验 数 据 进 行 二 次 多项 式 回归 处 理 ， 得到 的 回归方 程如 式 2 ] Y 一 一 5 7． 5 O1 4 23 2 41 2 ． 1 03 6 1 9 5 X2 4- 1 65 3 ． 77 3 1 4 5 5 X4 21 2 ． 60 4 3 8 8 8 6 X 3 X 3 4 37 0 ． 9 07 5 8 4 X 4 X4 4 - 1 5 8 0． 0 76 1 5 7 8 X 1 X2 2 9 2 6 ． 1 9 7 6 9 3 9 Xl X4 1 1 9 0 ． 2 5 5 9 1 8 7 X2 X4 2 3 5 ． 3 8 0 4 0 8 9 4 X3 X 2 对 该方 程进 行显 著性 检 验得 出 其相 关 系 数 R 为 0 ． 9 9 9 9 ， 显著性水平 为 0 ． 0 0 6 5 ， 剩余标准差 s为 0 ． 0 4 2 。因此， 该 回归方程在 实验范 围内， 回归方程能 较 好地 描述 混合 冷 剂 组成 与 最 小 冷 剂 循 环 量 的关 系 。 同时 ， X 未 引 入 方 程 ， 表 明 在 约 定 范 围 内 ， 相 应 的 c H 含 量 对 结 果 影 响 相 对 不 明显 ， 可 用 于 组 分 归 一 化 调整 。 以各 组分 合 理 范 围为 约 束 条 件 ， 以优 化 指 标 y 最 大化 为优 化 目标 ， 对 式 2 进行 最 优 化 求 解 ， 得 出冷 剂 组分 的优化 配方 见表 6 。优 化 计算 得 到 的 各 组分 摩 尔分数 加 和小 于 1 ， 为 保证 N 、 CH 、 C 。 H 、 C 。 H 的组 成为优化值 ， 需在 0 ． 6 ～1 范围内适当调整 x 含 量以使得组成归一化 。如表 6所示 ， 优化计算得到的 最小 冷 剂循 环量 为 2 0 1 4 ． 3 7 k mo l ／ h 。对 优 化 配 方 进 行模 拟 实 验 验 证 ， 得 到 最 小 冷 剂 循 环 量 为 1 9 7 6 k r n o l ／ h 。结果表明， 实验值 与优化计算值较 为接近 ， 且 低 于均匀 设计 的 全部 实 验 点 ， 从 而得 到约 束 条 件 下 天然气 液 化混合 冷 剂 的 优化 配 方 。通 过 拟合 方 程 2 可对该制冷系统不同冷剂配方下的最小冷剂循环量进 行 预测 ， 同时也可对不同条件的混合冷剂配 比进行优化。 表 6 混合冷剂制冷液化天 然气 均匀设计实验优化 结果 Table 6 Unlf o r m d es i g n ex per i men t a l op t i mi z a t i on r es u l t s b y mi x ed r e f r i g er an t f o r mul a t i on O 0 O 0 0 如 仰 如 ∞ 3 3 2，‘ 1 一 一 。 要 一 、 蚺 _軎黎 学兔兔 w w w .x u e t u t u .c o m 第 4 4卷第 3期 石油与天然气化工 CHEMI CAL ENGI NEERI NG OF OI L ＆ GAS 6 9 2 ． 4 因素 分析 根据得到的拟合方程 2 ， 将其他组分 固定为最优 配 比的条 件下 ， 分 析各 冷 剂 组 分 含 量 的变 化 对 制 冷 效 果 的影响 ， 结果 如 图 4 所 示 。 雌 颦 0 ． 1 2 0 ． 1 4 O ． 1 6 0 ． 1 8 O ． 2 0 0 ． 2 2 0 ． 2 4 O ． 2 6 0 ． 2 8 0 ． 3 0 0 ． 3 2 0 ． 3 4 Y 冷 剂组 成 ／ ％ 图4 混合冷剂各组分含量与冷剂循环量的关系曲线 F i gu r e 4 Rel a t ion sh i p b etwe en mol e c on t en t o f d i ff er en t c o mp on en t s a nd c i r c u l a t ion q ua n t it y of r e f r i g er an t 由图 4可 知 ， C 。 H。因 处 于 过 渡 温 区 ， 适 度 增 加 c 。 H 含量能够降低冷剂循环量 ， 提高制冷效率 ； 但 当 其增加至一定值后 最高点对应的冷剂组成 ， 过高 的 c 。 H 含量则会降低 总体 的制冷效率 。在 目前冷剂组 成 范 围 内 ， N 、 C H 、 C H 含 量 的 增 加 均 有 利 于 降低 冷剂 循环 量 ， 增强 冷剂 制冷 效果 。 表 7中的显 著性 检验 结果 也佐 证 了上述 结论 。 表 7实验因素显著性检验 Table 7 Sign i f i c a n ce t e s t o f exp er i men t a l f a c t o r s 由显著 性检 验结 果可 知 ， x 项 未 被选 入 方程 。因 此 ， 在加量范围内 C H 。 的加入 对天然气 液化 的影响 并不显著。根据显著性水平 P的大小可知 ， 各因素及 其 交互 作用 对冷 剂循 环量 的影 响 由大 到 小 的顺 序 为 X 2 X X 1 X2 X1 X X 2 X X3 X3 X 4 ， 其 中 的乘 积项表明， N z 、 C H 、 C H 、 C 。 H。等冷剂组分间存在一 定 的交 互作 用 。 3 结 论 本 文 通过对 天 然气 液化 混合冷 剂进 行理 论 和模 拟 实验研 究 ， 形成 了原 理初 选 、 模拟 实验 与均匀 设计 优 化 分析相结合 的天然气 液化混合冷剂配方高效评选 方 法 ， 主要研 究结 论 如下 1 不 同制 冷 剂 组 分 在 不 同 温 区制 冷 ， 应 根 据 天 然气 组 成确定 制 冷温 区 ， 从 而 选 择 对应 温 区 的混 合 冷 剂组 分 。 2 在 同一 制冷 温 区 ， 为 达到 相 同制 冷 效果 ， 分 子 量高的组分功耗小但制冷效果不足； 分子量低 的组分 制 冷效 果较 好 ， 但 功 耗大 。 因此 ， 需 均衡 选择 适宜 的冷 剂 组分 ， 同时在冷 剂 组 分筛 选 中应 尽 量 避 免 冷 剂组 分 的温 区重叠 。 3 提 出采 用原 理初 选 、 HY S YS模 拟实 验 与均匀 设 计优 化 分析相 结 合 的天然 气液化 混合 冷剂 配方 高效 评 选方 法 ， 即通 过均 匀 设 计选 取少 量 具 有 代 表 性 的模 拟实验 ， 回归分析混合冷剂中各组分的作用 ， 并通过优 化计算 确 定混 合冷 剂 的 优化 配 方 ， 对 优 化 配 方 模 拟 实 验验证 表 明 ， 此 方 法是 可行 的 。 参 考 文 献 E 1 ]赵敏 ，厉彦忠．C3 ／ MRC液化 流程 中原 料气成 分及制 冷剂组 分 匹 配[ J ] ．化工学报 ，2 0 0 9 ， 6 0 增刊 1 5 0 5 7 ． [ 2 ]曹文胜 ，吴集迎 ， 鲁雪生 ，等．撬装 型混合制冷剂液化 天然气流程 的热力学分析[ J ] ．化工学报 ，2 0 0 8 ，5 9 增刊 2 5 3 5 9 ． [ 3 ]王勇 ，张玉玺 ，白剑锋．L NG制冷 HYS YS计算模 型研究 E J ] ．天 然气 与石油 ，2 0 1 2 ，3 0 4 3 0 3 2 ． [ 4 ]尹全森．混合制冷剂循环优化设 计和动 态特性研究 E D] ．哈尔滨 哈尔滨工业 大学 ，2 O 1 O ． E 5 ]GE OF F RE Y G H， CHE N J ．Un i v e r s i t y o f L e e d s ，L e e d s L S 2 9 J T， U K A c o m p u t e r s i mu l a t i o n p r o g r a m f o r m i x e d r e { r i g e r a nt a i r c o n d i t i o n i n g [ J ] ．En e r g y ，2 0 1 3 4 6 1 5 5 ． [ 6 ]蒋浩 ， 高永和 ，王文军 ，等． MR C流程混合制冷剂的选择[ J ] ．煤气 与热 力，2 0 1 2 ，3 2 2 4 - 6 ． [ 7 ]张镨 ，鹿来运 ， 郭 开华．一种混合工 质遴选 、 配 比方法探 讨及实验 验证 口] ．低温工程 ，2 0 1 3 6 2 6 3 0 ． [ 8 ]毛东森．均匀设计及其在化学化 工中的应用 [ J ] ．化工进 展，1 9 9 7 5 1 6 - 1 9 ． [ 9 ]牛亚楠．多元混合制冷剂小型天然气液化装嚣的模拟研 究[ D] ．上 海 同济大学 ，2 0 0 7 ． [ 1 o ]王柱 ，方开泰．带有定 性因素均匀设计 的均匀性度量准则f- j ] ．数 据统计 与管理 ，2 0 0 0 ，1 9 3 2 8 3 2 ． [ 1 1 ]高辉 ，胡 良平 ，郭晋 ，等．如何正 确处理正交设计 和均匀设计 定 量资料 E J ] ．中西 医结合学报 ，2 0 0 8 ， 6 8 8 7 3 8 8 1 ． [ 1 2 ]姚钟尧．回归分析法在 均匀设计数据分析 中的地位[ J ] ．特种橡 胶 制 品 ，1 9 9 8 ．1 9 2 3 5 - 4 1 ． 收稿日期 2 0 1 5 0 4 2 4 ； 编辑 康莉 学兔兔 w w w .x u e t u t u .c o m