一种估算CO2和SO2在多元胺水溶液体系中溶解度的数学模型VIP

第 3 O卷 第 5期 2 0 1 3年 5月 2 8日 计算机与应用化学 Co mp u t e r s a n d Ap p l i e d Ch e mi s t r y Vo1 30， No 5 M a y 28 ，2 01 3 一 种估算 C O2 和 S O2 在多元胺水溶液体系中溶解度的数学模型 王端阳 ，李炼 ，李娜 ，周屈兰 ，荆强征 ，沈吉兆 ( 1 西安交通大学动力工程多相流国家重点实验室，陕西，西安，7 1 0 0 4 9 ； 2 西安特种设备检验检测研究院，陕西，西安，7 1 0 0 6 8 ； 3 江苏扬州晨光特种 设备有限公司，江 苏，扬 州，2 2 5 8 0 0 ) 摘要 ：通 过 离子反应 式 和亨利 定律 建立 数学模 型 ，在对模 型进 行适 当简化 的基础 上再 借助 编程手 段可 以得 到 C Oz - S O ：在 H2 O MDE A 活化胺复杂 多元溶液体系 中溶解度 的数值解 。研究表 明，在 MDE A 质量分率 3 O的水溶液 中 S O2 比 C 02 优先溶 解 ，溶液只有在基本 上完全 吸收 S O2 之后才会 吸收 C O2 ；而且其溶解度比较大，常温 常压 下 l mo l 的 MDE A可吸收超过 1 4 mo l 的 S O 2 气体 。在气 体分压 为( 1 1 0 0 ) k P a范围内，S O2 的溶解度 受压 力参数 的影响较 大：而在此压 力范 围之外 ，温度参数则成 为 主要影响因素。在 MDE A水溶液 中添加 ME A或 DE A活化剂，随着活化剂质量分率的增加，S O2 的溶解度有较 大程度 的下 降； 若同时添加 MEA DE A2种活化剂 ，S O2 的溶解度会随着 ME A质量 分率 的增加而减小 ，但减少 的幅度并不大 。 关键词：气体；数学模型 ；溶解度：MDE A：活化剂 中图分类号 ：T Q0 1 5 9 ； T P 3 9 1 9 ； 06 3 9 文献标识码：A 文章编号 ：1 0 0 1 - 4 1 6 0 ( 2 0 1 3 ) 0 5 4 9 2 4 9 6 D0I ：1 O 1 1 7 1 9 c o m a p pc h e m2 01 3 05 1 0 1 引言 脱除S O ： C O ： 提供更为多样可行的方案。 N 甲基二 乙醇胺( MDE A) 水溶液作为化学吸收剂被 广泛用于选择性脱除天然气等工业气体中的酸性气体如 二氧化碳和硫 化氢 ，其主要优点是脱碳和脱硫 的能力 高，腐蚀性小，降解程度低 J 。 MD E A 是叔 胺 ，不 能直接和 C O2反 应 ，只 能 由 MDE A的弱碱性来促进 CO 2 在水中的溶解 ，因而 MDE A 水溶液吸收 C O 2 的速率很小【 2 】 。如果在溶液中加入少量 活化剂可 以大大提高 MD E A水溶液吸收C O 2 的速率，常 见 的活 化 剂 包 括 一 乙醇 胺 ( ME A) 和 二 乙醇 胺 ( D E A) 等 。 近 年来 ，关于酸 性气 体尤其是 C O 2 或 H 2 S在 以 MDE A为基础 的选择性脱除溶剂中气液平衡的模型和实 验的文献在国内外屡见不鲜 。国外，从早期 Au s t g e n D M【 4 】 等获得 H2 S C O2 在胺水溶液中气液平衡的模型和数 据到后来 Ma mu m s 【 5 等在更加宽泛的温度压力范围精 确地测量出 C O： 在 MDE A 中的溶解度数据，再到后来 Al i H【 6 】 等采用电解质 UNI QUAC - NR F模型进一步优化了 这一吸收过程 ；国内，徐 国文【 等最早测量 出 C O 2 在 MDE A水溶液 中的溶解度，陈健【 】 】 等又把 H 2 S加入气液 平衡体系进行了研究分析，而陆建刚【 8 】 、徐宏建【 9 】 等人 分别将溶液体系以MDE A为基础进行了不同的扩展，并 提 出了相关的模型 。 S O 作为 电厂烟气中主要排放成分之一，到 目前为 止关于其在胺水溶液中溶解度 的报道却很少。因此，有 必要对 S O2 以及 S O2 C O2 同时在胺水溶液尤其是 MDE A 溶液中溶解度性质开展基础工作的研究，为工业 中联合 2 气液平衡模型的描述 2 1 基本模型的建立 选择性脱除剂胺水溶液吸收 C O ： S O 气体溶液体系 中发生的离子反应如下： MDE A H 山H + MDE A ( 1 ) ME A山H + ME A ( 2 ) DE A 蔓 H + D E A ( 3 ) H2 o 生 _ _ H oH一 ( 4 ) C O 2 ( a q ) + H 2 o H + H C O 3 一 ( 5 ) S O 2 ( a q ) + H 2 o 山H + H S O 3 一 ( 6 ) HC Oj 一 L H + C O 一 ( 7 ) HS O3 一七 H + S O3 一 ( 8 ) 其中， 岛 ( ， = 1 ， ， 8 ) 为反应平衡常数 ，此平衡常数 已考虑离 子、分子的活度 。 为作简化，忽略胺水溶液中的 C O3 2 - 与 S O3 2 离子 ， 因 为这些 离子 与其 他类 型 的分 子离子 相 比浓度 低 得 多 依据电荷守恒和质量守恒可得： A mi n e Hi + + H = HC O3 -+ HS O3 一 ( 9 ) 一 m i =l A m i n e H l + A m i n e ( 1 o ) m a 。= H C O 3 - + c o ( a q ) ( 1 1 ) m as= H S 0 3 - + s o ： ( a q ) ( 1 2 ) m = ( 1 3 ) 收稿 日期：2 0 1 3 - 0 2 - 2 0 ；修回日期：2 0 1 3 0 4 0 8 基金项 目：国家 自然科学基金资助项 ( 5 1 0 7 6 1 2 7 ) ；国家级大学生创新训练项 目( 2 0 1 2 1 0 6 9 8 0 1 4 ) 作者简介：王端( 1 9 9 o _ _ ) ，男，河南人，硕士研究生，E ma i l ： w d y 5 5 s t u x j tu e d u c n 联系人：周屈兰( 1 9 7 4 一) ，男，四川人，博士，教授，博士生导师，E ma i l ： q l z h o u m a i l x j t I1 e d u c a 2 0 1 3 ， 3 0 ( 5 ) 王端阳， 等： 一种估算C O 2 和 S O 2 在多元胺水溶液体系中溶解度的数学模型4 9 3 其 中，下标 C 、S 分别指 C O 2 与 S O 2 气体。 依据亨利定律可得： P c o ： = H c L C O z ( a q ) j ( 1 4 ) o： l s o ： ( a q ) J ( 1 5 ) 在一定温度和压力下 日 与 K， 都可以看作为定 值，而 m也是一确定的常量，本模型的关键在于通过确 定氢离子浓度从而得到酸性气体 C O ： 与 S O ： 溶解度值。 因此上述的反应和等式可以继续进行如下的变形 。 由式( 5 ) 、式( 1 4 ) 可得 ： m w + + c , ) H P c o 2 ( 1 6 ) 同理 由式( 6 ) 、式( 1 5 ) 可得： n w + + c 6 ) H e s o ( 1 7 ) 由式( 9 ) 、式( 1 0 ) 、式( 1 1 ) 与式( 1 2 ) 可得 ： 呻 = 南+ + H 】 两边 同时乘 以H 可得： + 南 c o ： (aq ) s o ： (aq ) (19) 在一定的温度、压力下等式右边的值为常数，而等 式左边只有 H 未知，下面根据溶剂 的不 同分别进行讨 论。 2 2 MDEA水溶液吸收 C O2 S O 模型 若溶液中只有 MD E A，可将式( 1 9 ) 进行化简变形得 如下方程： H j + ( + 毛 ) H 2 - A H A = 0( 2 0 ) 式中： A = c o z ( a q ) + k 6 S O ( a q ) ( 2 1 ) 式( 2 0 ) 中 H 3 可以忽略， H 的系数( + 毛) 中， k 1 可以忽略 8 】 ，于是继续将式( 2 1 ) 整理如下： 可以求出： H + = A + x A2 + 4 Ak l ml ( 2 3 ) 2 将 求出的 H 】 值代入式( 1 6 ) 、式( 1 7 ) 解 出 与 的值。 2 3 CO2 与 S O2 在混合溶剂中溶解度的求取模型 对于 2种及 2种以上溶剂的胺水溶液 中酸性气体的 吸收，不再像单种 MD E A溶液吸收模型一样得到比较精 确的解析解 。但是借助计算机编程，同样可 以得到十分 精确 的数值解，甚至优于单种吸收模型的解析解 。 继 续分 析式 ( 1 9 ) ，等式 中只 有 H 为变 量 ，利用 MA T L AB编程，采取循环逐步逼近的方法，将误差控制 在一定的范围，从而可 以得到 H 的值。 P_ l ! I x l 0 0 (2 4 ) l A l 其 中 ， 为第 步 迭代 等 式 左 边 的值 ；e 为 MA T L A B迭代控制误差精度 ，单位为。 规定 e 0 0 1 时， H 的值作为方程的准确解。 2 4 模型 中的参数估算 对于本模型计算过程中所需的参数 ，有 k 1 以及亨 利常数 风和 。对于这些参数的取值可 以作如下处 理 。 所涉及的反应平衡系数 k 的计算采用温度因变函 数，其中 k i - k 5 可以在文献【 】 中查到，而 可以借助化工 软件 As p e n的数据库查得 ，数据列表见表 l 。 l n k =- C 1 + C T+ c 3 l n T + c 4 T ( 2 5 ) 其 中，C i ( i = l ， 2 ， 3 ， 4 ) 为温度 因变平衡常数表达式系数 。 其 中所涉 及 的亨利系数 ，凰可 以从文献 4 中查 得，而 凰 可以从 As p e n软件 自带的数据库中查得 。 1 n 日c=1 7 0 7 8 4 7 7 7 一2 1 9 5 7 4 3 1 n T+0 0 0 5 7 8 1 T( 2 6 ) l n nS =8 3 9 6 0 6 5 5 7 8 8 一8 7 6 1 5 2 T ( 2 7 ) 式( 2 6 ) 、式( 2 7 ) 中 的单位为 K。 m l H 】 - A H 】 A= 0 ( 2 2 ) 表 l 温度因变反应平衡常数表达式 中相关系数 Ta b l e 1 T e mp e r a t u r e d e o e n d e n c e o fc o u i l i b r i u m c o n s t a n t s f o r r e a c t i o n s I - 6 3 模型结果及讨论 3 1 模型的验证与分析 本文首先采用所建立的模型计算 了 C O2 气体在 MD E A水溶液中随温度、压力变化 的溶解度 曲线，并与 文献值进行了比较 。 F a n g Y u a n J o u I n 等人测得的有关 CO2 溶解度的实验 4 9 4 计算机与应用化学 2 0 1 3 ， 3 0 ( 5 ) 数据与 S h e n和 L 】 测得的数据 原本就存在一定的差 异。而从图 1中可以看出，模拟的结果又在一定程度上 偏离 了实验数据，存在误差的范围主要在低温 区域与 中 压区域 。模型中完全忽略 C O3 2 - 离子、推导过程中进行 相关简化、离子平衡常数 以及亨利常数关系式只适用于 特定温度、压力范围等因素是造成模拟数据偏小的主要 原因。一方面，本文 旨在快速研究、预测 C O2 S O 2 在活 化 MD E A水溶液中溶解变化规律 ，对其溶解度的精确值 要求并不高，而且未来工作 中会加入实验部分以对溶解 度进行精确测量 ；另一方面 ，通过本模型模拟得到结果 的误差不大 ，在误差允许的范 围之 内。因此，没有必要 再对数学模型进行修正，本文正是基于此模型展开相关 的研究工作。 1 O l 0 重m 。 1 0 。 1 0- 1 0 2 0 4 0 4 0 6 0 8 1 0 1 - 2 1 4 1 6 c o2 l o a d in g ( m o l o f C O 2 to o l o f MDE A) C O2 ( m o l o f C O m o l o f MDE A) F i g 1 Co mp a r i s o n o f s o l u bi l i t y o f CO2 i n 3 0 wt M DEA s o l u t i o n 图 1 C O 2 在质量分数 3 0 的MD E A水溶液中溶解度的对 比 3 2 S O2 一 C O2 在 MDE A水溶液中溶解度 S O 2 C O 2 气体的分压、温度 以及胺水溶液有机胺的 溶度是影 响溶解度的主要因素，分别对这些影响因素进 行分析可 以获得溶解度变化的规律。 当 MD E A 的质量浓度为 3 0，分别考虑温度和压 力对溶解度的影响。经过计算并使用 O r i g i n软件绘制得 到图 2与图 3 。 0 2 0 4 0 6 0 8 0 1 0 0 l 2 O 1 4 0 1 6 0 T o c F i g 2 I n f l u e n c e o f t e mp e mt u r e o n t h e s o l u b i l i ty o f C O2 一 S O2 图 2 温度对 C O： S O 溶解度 的影 响 图 2是在气体分压为 1 0 0 k P a时，单独 CO 2 或 S O2 气体 以及混合 的CO 2 与 S O2 气体在 MDE A水溶液中随温 度变化的溶解度 曲线 。从图 1中可以知，在条件相同的 情况下，S O 的溶解度明显高于 C O2 的溶解度，而随着 温度 的升高，气 体在 溶液中的溶解 度均逐渐降低 。在 2种气体混合的时候，S O2的溶解度 曲线几乎与单独气 体的溶解度 曲线重合 ，而 C O2 的溶解度近乎为 0 ，这就 说明MD E A水溶液对于 S O 有优先选择吸收的性质，对 于 C O2 的溶解吸收则建立在 S O 吸收完毕的基础上 。 l 0 1 0 1 0 4 1 05 p P a F i g 3 I n fl u e n c e o f p r e s s u r e o n the s o l u b i l i ty o f C O2 一 S 02 图 3 压力对 C O2 - S O 2 溶解度的影响 图 3是在温度为 2 5时，单独 C O2 或 S O2 气体 以 及混合的 C O2 与 S O 2 气体在 MDE A水溶液中随着压力变 化 的溶解度 曲线。从图 2中可以看 出，单种气体时，在 低压区间 S O 2 的溶解度 明显高于 C O 的溶解度 ，而随着 压力的增 加，二者的差距逐渐减小 。对于单种气体 ， CO 2的溶解度几乎呈线性上升趋势 ，而 S O： 的溶解度斜 率不大，在超过 1 1 K P a的时候才有 明显的增长。同样 地，2种气体混合下，CO2的溶解度趋于 0 ，也反映出 MDE A水溶液优先吸收 S O， 。 3 3 S O2 在 MDE A- ME A水溶液 中溶解度 图 4是 S O 2 在温度为 2 5，质量分数为 2 4的 MDE A 和 6 ME A 的溶液中溶解度随气体分压的变 化 。从图中可也看出压力在( 1 1 0 )k P a的范围内，气体 溶 解度 曲线 的斜 率较 大 ，而且 摩 尔分 率在 不 同温 度 ( 2 0 1 0 0 )时差别不大，这就说明了在此区间内温度不 是影响溶解度 的主要 因素，压力 的影响占据主导。而在 此压力范围之外 ，溶解度差异开始变大 ，说明了温度因 素开始发挥作用。 0 0 0 5 1 0 I 5 2 0 2 5 3 0 3 5 4 0 S O， mol e f r a c t i o n mo l mo l S O 摩尔分率， m0 l mo l。 Fi g 4 I n f l u e nc e o f p r e s s u r e a n d t e mp e r a t u r e o n the s o l u b i l i t y o fS O2 i n M DE A M EA a q u e o u s 图4 压力温度对 S O2 在 MDE A ME A水溶液中溶解度的影响 4 2 0 8 6 4 2 O 1 1 0 0 O O O 一0 暑 10 暑， u 0 I 】u 鐾 0 吕 p 1 0 、 ， 一0 gl Io m 料 赶 l。 m m 一 0 昌 一 0 暑 0 一 趸 白 0 g 。 Ip l 。 一 _1 0 g_I o u v 求 掣 2 0 1 3 ， 3 0 ( 5 ) 王端阳， 等：一种估算C O 2 和S O 2 在多元胺水溶液体系中溶解度的数学模型4 9 5 图 5是研究 S O 2 在不同胺浓度下溶解度的差异性 。 从图中可以看到，随着 MDE A浓度 的降低，S 02 的溶解 度逐渐减小。这说明虽然加入的活化剂 ME A 能适当的 增加整个过程 的吸收速率 ，可也在一定程度上 降低 了 S O2的溶解度。于此 同时，在高温 区间( 6 0 1 0 0 ) S O 2 溶解度 曲线差异较小，这表 明超过一定温度范围，溶解 度主要与溶液的成分有关系。因此选择合适浓度的胺溶 液 以及适当的温度区间对于 S O 气体的吸收和溶解具有 十分重要的工程意义。 ：垂 20 4 0 60 8O 1 OO MDEA ma s s ffa c i t i o n i n mi x t u r e s o fMEA +MDEA s o l u t i o r d MD E A 在混 合ME A 与MD E A溶液中质量分数 Fi g5 S O2 s o l u b i l i t y c u r v e s i n di ffe r e n t ma s s r a t i o s o l u t i o ns o f M DEA M EA 图 5 S O2 在不同质量配 比下MD E A ME A水溶液中溶解度曲线 3 4 S O2 在 MDE A ME A- DEA水溶液 中溶解度 S O 2 在分压为 1 0 0 k P a ，溶液中MDE A的质量浓度为 2 4，H2 O的质量浓度为 7 0，研究在ME A DE A不同 质量配比下 S O 2 溶解度变化 的性质 。经过计算，可以得 到图 6 。 从图 6 可以看出，随着活化剂 ME A浓度的增加 即活 化剂 DE A 浓度 的减小，S O2 在溶液 中的溶解度逐渐减 小 。但在相 同温度下，S O 的摩尔分率值变化并不大 。 这就说明添加的这 2种活化剂对 以MDE A为基础的脱除 剂中 S O 2 的溶解度的影响并不大。因此，挑选适当的活 化剂则主要从活化剂 自身的价格 以及提升反应速率的性 能方面考虑 。 1 2 = _趸 1 1 至 10 篓 毂毛 0 8 2 0 4 0 6 0 8 0 l 0 O T o C F i g 6 S O2 s o l u b i l i ty c u r v e s i n d i f f e r e n t ma s s r a t i o s o l u t i o n s o f MD EA M EA DEA 图 6 S O 2 在不同质量配比下 MD E A ME A DE A水溶液中溶解度曲线 4 结果与讨论 本文利用物理模型和编程技术，比较系统地计算了 温度( 2 0 1 0 0 )，压力( 0 0 1 1 0 0 0 ) k P a 范围内，酸性气 体 C O2 S O2 在 以MDE A为基础的选择性脱碳脱硫剂中溶 解度 曲线。具体结论如下： ( 1 )虽然 C O 2 在胺 水溶液中的溶解度并不小 ，可 C O2 的溶解能力却远低于 S O2 ，选择性吸收剂 MDE A溶 液往往先吸收 S O 2 之后才吸收 C O2 ； ( 2 )加入 的活化剂 ME A 和 DE A 虽然能够 增强 MD E A 水溶液溶解吸收 S O 2 的能力，但会降低 S O 2 的溶 解度 。因此 ，选择合适 的胺水溶液 的浓度对于实际工程 的经济性具有十分重要的意义。 ( 3 )在 MDE A水溶液中同时添加多种活化剂，不 同 活化剂对于 S O 在以MDE A为基础的脱硫溶液中溶解度 的影响并不大。 ( 4 )本文只针对气液平衡状态下 的溶解度进行 了模 拟与分析，并未对实际溶液吸收有害气体的速率规律展 开研究，这也是本作者今后将致力的方向。 综上所述 ，本研究为电厂采用循环再生有机溶液联 合脱除 C O S O ： 提供了较为详细、系统的理论数据 ，也 为快速研究、预测其他气体在混合溶液中溶解度提供了 较为简便易行 的方法 。 符号说明 q =l ， ， 8 ) 反应平衡常数 不同种类胺的物质的量，too l m 胺的总物质的量，mo l ： 胺水溶液中 C O 2 的亨利系数，mo l mo 1 日s 胺水溶液中 S O 2 的亨利系数，mo Y mo l 0 2 C O 2 平衡分压 ，P a P s 0 2 S O 2 平衡分压，P a 溶解的 C O 2 与胺的物质的量之比，mo l mo l 溶解的 s O 2 与胺的物质的量之比，mo l mo l Re f e r e nc e s ： 1 Che n J i a n ， M i J i a n g u o a n d Li u J i n c h e nCa l c u l a t i o n of g a s s o l u bi l i t y i n t he s y s t e m M DEA H2 O CO2 一 H2 S Na t u r a l Ga s C h e mi c a l I n d u s t r y , 2 0 0 1 2 6 ( 3 ) ： 5 7 6 1 2 Xu eQu a n mi n ， Z h o uY a p i n g a n d S uWe i S tu d y o n a b s o r p t i o n o f CO2 i n t o a q u e o us N me t h y l d i e tha n o l a m i n e Ch e mi c a l En g i n e e r i n g ( Ch i n a ) ， 2 0 0 9 ， 3 7 ( 9 ) ： 1 4 3 M a nd a l B Gu h a M ，Bi s wa s A K a n d Ban dy o pa d h y a y S S Re mo v a l o f c a r b o n d i o x i d e b y a b s o rpt i o n i n mi x e d a mi n e ： mo d e l i n g o f a b s o r p t i o n i n a q u e o u s M DEA MEA an d AMP MEA Ch e m i c a l Eng i n e e rin g S c i e nc e ， 2 0 0 1 ， 5 6 ： 6 21 7 - 6 2 2 4 4 Aus e n D M an d Ro c h e l l e G M o d e l of v a p o r - l i q u i d e q u i l i b ria f o r a q ue o u s a c i d g a s a l k an o l a mi n e I n d En g Ch e m Re s ，1 9 91 ， 3 0 ( 3 ) ： 5 4 3 5 5 5 5 M a mu m SNi l s e n R a n d S ve n d s e n H F S o l u b i l i t y o f c a r b o n d i o x i d e in 3 0 ma s s mo n o e t h a n o l a mi n e an d 5 0 ma s s me thyl d i e t h an o l a m i n e s o l u t i o n s J Ch e m En g Da t a ， 2 0 0 5， 5 0 ： 6 3 0 - 6 3 4 6 Al i H，Mo h a mm an d D El e c t r o l y t e UNI OUI C NRF mo d e l t o s t u d y t h e s o l u bi l i t y o f a c i d g a s e s i n a l k a n o t a mi n e s I nd En g Ch e m Re s ， 2 0 0 7 ， 4 6 ( 1 8 ) ： 6 0 5 3 6 0 6 0 6 4 2 O 8 6 4 2 0 l 1 0 O O O O 4 9 6 计算机与应用化学 2 0 1 3 ， 3 0 ( 5 ) 7 Xu Gu o we n Z h a n g Ch e n g f a n g a n d Qi n S h u j u n S o hb i l i t y o f C02 ， H2 ， N2 i n a q u e o us M DEA s o l u t i o n s J o u m a l o f Ch e mi c a l E n g i n e e r i n g o f C h i n e s e Un i v e r s i ti e s ， 1 9 9 6 ， 1 0 ( 4 ) ：4 0 6 - 4 1 0 8 Lu J i a n g a n g， Li uCo n g ， F an Fan， J i Yan a n dZh a n gHl l i S o l u b i l i t y o f c a r b o n d i o x i d e i n mi xe d s o l u t i o n o f MDEA AM J o u r n a l o f F u e l Ch e mi s t r y and T e c lm o l o g y , 2 0 1 0 ， 3 8 ( 6 ) ： 7 3 0 7 3 2 9 Xu H o n g i i an ， Z h an g Hu i x i n and Z h an g C h e n g f ang ， S o l u b i l i t y o f h y d r o g e n s ul fi de an d c a r b o n d i o x i d e i n a s o l u t i o n o f me t h y l di e t han o l a mi n e mi x e d wi t h e t h y l e n e g l y c o 1 Ch e roi c a l Re a c t i o n E n g ine e n n g and r e c h n o l o g y , 2 0 1 0 。 2 6 ( 3 ) ： 2 6 9 2 7 4 1 0 Oi u S h a n u a n g L i u Y 0 u z h i and B a i Me i A s i mp l e mo d e l f o r s o l u b i l i t y o f H 2 S and C O2 i n MDE A i n a i r t i pJ a t c o mp a r t me n t Ch e ro i c a l I nt e r m e d i a t e s 2 0 1 1 。 3 ： 3 0 3 2 1 l J o u F Y Ot t o F D a n d M a th e r A E Va p o r - Li q u i d e q u i l i b riu m o f c a r b o n d i o x i d e in a q u e o u s mi x t u r e s o f mo n o e t h a no l a mi ne an d me thyl d i e t h a no l a mi ne Ind En g Ch e m Re s 1 9 9 4 。 3 3 ： 2 0 0 2 2 0 0 5 1 2 Li M H an d S h e n K P _De ns i t i e s an d s o l u b i l i t i e s o f s o l u t i o n s of c arbo n d i o x i d e in wa t e r + m o no e t h an o l a mi n e + N- m e t h yl d i e t h a n o l a mi n e J Ch e m En g Da t a ，1 9 9 2 ， 3 7 ： 2 8 8 2 9 0 中文参考文献 l 陈建，密建 国，刘金晨 MDE A H2 O C O2 一 H2 S体系 的气体溶解 度计算 J 天然气化工， 2 0 0 1 ， 2 6 ( 3 ) ： 5 7 6 1 2 薛全 民，周 亚平，苏伟 甲基 二乙醇胺水溶液 吸收 C O2的研 究 J 化学工程， 2 0 0 9 ， 3 7 ( 9 ) ： 1 - 4 7 徐国文，张成芳，钦淑 均 C O2 、H 2 、N2 在 MDE A 水溶液中的 溶解度 J 高校化学工程学报， 1 9 9 6 ， 1 0 ( 4 ) ： 4 0 6 -41 0 8 陆 建 刚 ，刘 聪 ，樊 瑶 ，嵇 艳 ，张 慧 C O：在 复 合 溶 液 MDE A A MP 中 的 溶 解 度 J 燃 料 化 学 学 报 ， 2 0 1 0 ， 3 8 ( 6 ) ： 7 3 0 - 7 3 2 9 徐宏 建，张 慧鑫 ，张成 芳甲基 二乙醇胺 乙二胺 水混 合溶剂 中 C O2 和 H2 S的溶解度 J 化学反应工程与工艺， 2 0 1 0 ， 2 6 ( 3 ) ： 2 6 9 2 7 4 1 O 邱 尚煌，刘有智，白梅密 闭空 间中 H2 S ， CO 2 在 MDE A中的理 论溶解度模 型 J 化工中间体， 2 0 1 1 ， 3 ： 3 0 3 2 A s i mpl e m a t he m a t i c a l m o de l f o r e s t i ma t i n g t he s o l ub i l i t y o f CO2 a n d 802 i n t he pl o y - a mi ne s a que o us s o l u t i o ns Wa n g Du a n y a n g ， L i L i a n ， L i Na ， Z h o u Qu l a n ， J i n Qi a n g z h e n g 2 a n d Sh e n J i z h a o ( 1 S t a t e Ke y L a b Mu l t i p h a s e F l o w, Xi A n J i a o T o n g U n i v e r s i t y , X i a n 7 1 0 0 4 9 ， S h a a n x i ， C h i n a ) ( 2 Xi a n S p e c i a l E q u i p me n t I n s p e c t i o n a n d T e s t i n g I n s t i t u t e ， X i an 7 1 0 0 6 8 ， S h a a n x i ， C h i n a ) ( 3 Y ang z h o u C h e n G u a n g S p e c i a l E q u i p me n t C o L T D， Y ang z h o u 2 2 5 8 0 0 ， J i a n g s u ， C h ina ) Ab s t r a c t ： A r e l a t i v e l y a c c u r a t e n u me ri c a l s o l uti o n t o the s o hb i l i ty o f C O2 - S O 2 in H2 O- MDE A DE A M匝A c o u l d b e 0 b t a i n e d wi th the h e l D o f s p e c i fic p r o g r a mmi n g me tho d and a s i mp l i fie d ma the ma t i c a l mo de l b