基于同步热跟踪法的微量气液反应热测量

工程科学学报，第41卷，第3期：368-376,2019年3月 Chinese Journal of Engineering,Vol.41,No.3:368-376,March 2019 DOI:10.13374/j.issn2095-9389.2019.03.010:http://journals.ustb.edu.cn 基于同步热跟踪法的微量气液反应热测量 张宣凯”，张挥2)，李东”，刘应书”，张会元”，王润”，郭亚楼” 1)北京科技大学能源与环境工程学院，北京1000832)北京科技大学治金工业节能减排北京市重点实验室，北京100083 ☒通信作者，E-mail:zhanghui56@usth.cdu.cn 摘要利用同步热跟踪原理，提供一种测定微量气液反应热的研究方法.通过程序控制容器外壳温度与内部溶液同步升 温，减小温度梯度，形成“热屏障”，阻止溶液以热传导、对流、辐射的形式与外界环境进行热交换，获得动态绝热环境，提高微 量气液反应热直接测量的精度，减少样品用量，无需热补偿.采用MEA(乙醇胺)与MDEA(N甲基二乙醇胺)两类弱碱吸收 液，容积为15mL,分别在10%、20%、30%、40%和50%质量分数下，测定吸收C02的反应热.实验表明：同步热跟踪法测量更 为准确：随溶液浓度的增加，MEA反应热先降低后升高，MDEA反应热逐渐降低：在质量分数为20%~40%时，MEA、MDEA质 量分数对反应热的影响不显著：反应放热形成的升温曲线出现“下凹”现象 关键词同步热跟踪：微量反应热：吸收：化学反应器 分类号TK113 Measurement of gas-liquid reaction heat based on synchronous thermal tracking ZHANG Xuan-ai,ZHANG Hui,LI Dong,LIU Ying-shu,ZHANG Hui-yuan,WANG Run,GUO Ya-ou 1)School of Energy and Environmental Engineering,University of Science and Technology Beijing.Beijing 100083,China 2)Beijing Key Laboratory for Energy Saving and Emission Reduction of Metallurgical Industry,University of Science and Technology Beijing,Beijing 100083,Chia Corresponding author,E-mail:zhanghui56@ustb.edu.cn ABSTRACT In the process of CO2 capture by chemical absorption,regeneration energy consumption accounts for 70%-80%of the total energy consumption.Currently,the most critical issue is how to reduce the energy consumption of regeneration.Equipment such as micro-reaction calorimeter (Thermal Hazard Technology provides),differential reaction calorimeter and Setaram C80 thermal differ- ential calorimeter is used to compare the reference and sample solutions,which are simultaneously heated to compensate for heat loss of the sample solution during the measurement,but the heat of reaction cannot be directly measured.In this study,the reaction heats of MEA (ethanolamine)and MDEA (N-methyldiethanolamine)with CO2 at 10%,20%,30%,40%,and 50%mass fraction were measured by synchronous thermal tracing technique.By synchronously controlling the temperature of the shell of the container and the internal solution,the temperature gradient was reduced to form a "thermal barrier"to prevent the solution from exchanging heat with the external environment in the form of conduction,convection,or radiation.A dynamic adiabatic environment was obtained without thermal compensation.The accuracy of direct measurement of the trace gas-liquid reaction heat was improved to save the sample amount.The experimental results show that the simultaneous thermal tracking method is more accurate.With the increase of solution concentration,the reaction heat of MEA first decreases and then increases,and the reaction heat of MDEA decreases gradually.When the mass concentration of MEA and MDEA is between 20%and 40%,the mass concentration has no significant effect on the reaction heat.The curve of temperature rise formed by exothermic reaction appears to be concave. KEY WORDS synchronous thermal tracking:trace reaction heat:absorption:chemical reactor 收稿日期：201803-20 基金项目：北京市科委项目基金资助项目(D161100006016001)

工程科学学报，第 41 卷，第 3 期: 368--376，2019 年 3 月 Chinese Journal of Engineering，Vol． 41，No． 3: 368--376，March 2019 DOI: 10． 13374 /j． issn2095--9389． 2019． 03． 010; http: / /journals． ustb． edu． cn 基于同步热跟踪法的微量气液反应热测量 张宣凯1) ，张 辉1，2) ，李 东1) ，刘应书1) ，张会元1) ，王 润1) ，郭亚楼1) 1) 北京科技大学能源与环境工程学院，北京 100083 2) 北京科技大学冶金工业节能减排北京市重点实验室，北京 100083 通信作者，E-mail: zhanghui56@ ustb． edu． cn 摘 要 利用同步热跟踪原理，提供一种测定微量气液反应热的研究方法． 通过程序控制容器外壳温度与内部溶液同步升 温，减小温度梯度，形成“热屏障”，阻止溶液以热传导、对流、辐射的形式与外界环境进行热交换，获得动态绝热环境，提高微 量气液反应热直接测量的精度，减少样品用量，无需热补偿． 采用 MEA( 乙醇胺) 与 MDEA( N-甲基二乙醇胺) 两类弱碱吸收 液，容积为 15 mL，分别在 10% 、20% 、30% 、40% 和 50% 质量分数下，测定吸收 CO2的反应热． 实验表明: 同步热跟踪法测量更 为准确; 随溶液浓度的增加，MEA 反应热先降低后升高，MDEA 反应热逐渐降低; 在质量分数为 20% ～ 40% 时，MEA、MDEA 质 量分数对反应热的影响不显著; 反应放热形成的升温曲线出现“下凹”现象． 关键词 同步热跟踪; 微量反应热; 吸收; 化学反应器 分类号 TK11 + 3 收稿日期: 2018--03--20 基金项目: 北京市科委项目基金资助项目( D161100006016001) Measurement of gas--liquid reaction heat based on synchronous thermal tracking ZHANG Xuan-kai1) ，ZHANG Hui1，2)  ，LI Dong1) ，LIU Ying-shu1) ，ZHANG Hui-yuan1) ，WANG Ｒun1) ，GUO Ya-lou1) 1) School of Energy and Environmental Engineering，University of Science and Technology Beijing，Beijing 100083，China 2) Beijing Key Laboratory for Energy Saving and Emission Ｒeduction of Metallurgical Industry，University of Science and Technology Beijing，Beijing 100083，China Corresponding author，E-mail: zhanghui56@ ustb． edu． cn ABSTＲACT In the process of CO2 capture by chemical absorption，regeneration energy consumption accounts for 70% － 80% of the total energy consumption． Currently，the most critical issue is how to reduce the energy consumption of regeneration． Equipment such as micro-reaction calorimeter ( Thermal Hazard Technology provides) ，differential reaction calorimeter and Setaram C80 thermal differ￾ential calorimeter is used to compare the reference and sample solutions，which are simultaneously heated to compensate for heat loss of the sample solution during the measurement，but the heat of reaction cannot be directly measured． In this study，the reaction heats of MEA ( ethanolamine) and MDEA ( N-methyldiethanolamine) with CO2 at 10% ，20% ，30% ，40% ，and 50% mass fraction were measured by synchronous thermal tracing technique． By synchronously controlling the temperature of the shell of the container and the internal solution，the temperature gradient was reduced to form a“thermal barrier”to prevent the solution from exchanging heat with the external environment in the form of conduction，convection，or radiation． A dynamic adiabatic environment was obtained without thermal compensation． The accuracy of direct measurement of the trace gas--liquid reaction heat was improved to save the sample amount． The experimental results show that the simultaneous thermal tracking method is more accurate． With the increase of solution concentration，the reaction heat of MEA first decreases and then increases，and the reaction heat of MDEA decreases gradually． When the mass concentration of MEA and MDEA is between 20% and 40% ，the mass concentration has no significant effect on the reaction heat． The curve of temperature rise formed by exothermic reaction appears to be concave． KEY WOＲDS synchronous thermal tracking; trace reaction heat; absorption; chemical reactor

张宣凯等：基于同步热跟踪法的微量气液反应热测量 ·369· 在化学吸收法捕集C0,的过程中，再生能耗占 散失越严重.在反应容器的外部设置同步升温壁 C0,捕集过程总能耗的70%~80%m,目前最关键 垒，与内部溶液形成小的温差，构成一道“热屏障”， 的问题是如何降低再生能耗-)，因此，气液反应热 由于温度梯度减小，热散失驱动力不足，从而减少热 的测量对筛选低能耗吸收剂具有重大意义 损失.传统的温度控制方法多采用PID控制系统根 微反应量热计(Thermal Hazard Technolog罗公司 据温度设定值进行控制6-，而热屏温度的控制是 提供)t、DRC系统量热计)、Setaram C80导热差 追踪变化的设定点，并且不能采用“波浪形”温控曲 示量热仪80等设备，均采用参考溶液与样品溶液 线，如果热屏温度高于内部，会导致向系统输入热 同时加热进行对比的方式来抵消样品溶液在测量中 量，造成反应热增大，因此，热屏温度应“尾随”并且 的热损失，但不能直接测量反应热，存在较大误差. 最大限度的“逼近”内部溶液的温度.本研究采用固 Kim和Svendsen用容积为20O0mL的机械搅 态继电器替代机械继电器缩短加热电路通断时间： 拌反应器（瑞典ChemiSens AB,型号CPA122)进行 采用上位机自由编程读取溶液温度和控制热屏温 量热，无保温措施：郭东方等回使用250mL真实热 升，从而实现热屏温度与内部溶液温度的同步，造成 流反应器（型号CPA201),将反应器置于恒温液体 “零”温度梯度。具体实验流程如图1所示，c系统与 中，器壁采用双层玻璃形成中间真空层，虽然采用了 m系线表示反应器系统各项的比热容与质量，Q为反 保温措施，但仍无法阻止热散失.上述两种方法均 应过程放热量. 采用增加实验样品量的方式，使散热损失占总放热 1.3实验装置 量的比例变小，可近似忽略，不适用于微量气液反应 所有实验仪器主要有：温度变送器，A1702,厦 热的测量 门宇电自动化科技有限公司，测量精度0.2级；单相 Svensson等2使用化学过程分析仪（型号 电参数盘面表，DSW9002C,创鸿（台湾）电子仪表 CPA202),将容量为250mL的反应器置于恒温液体 有限公司，测量精度0.4级，0.02W:数采控制器， 浴中，使其主动随着反应器升温而升温，但对于微量 JYDAM0808,北京聚英翱翔电子有限公司：进气数 样品的情况，恒温液体浴的升温过程响应驰豫时间 字式质量流量控制器，CS200,北京七星华创流量计 长，控制精度低，且反应器顶部和底部还可进行热传 有限公司，测量精度0.35级；出气质量流量传感器， 导，实验误差大 FS4001200L,美国矽翔微机电系统有限公司，测 本文提出一种同步热跟踪法测微量气液反应 量精度1.5级：磁力加热搅拌器，79」，上海梅香仪 热，程序控制外部热屏障自动与反应器温度保持同 器有限公司：集热式搅拌器，DF101S,金坛市科析 步，响应快速，无需热补偿，可以减小上述方法所带 仪器有限公司 来的误差，适用于微量气液反应热的测量，计算与操 测量质量精度为0.03%，温度变送器设置量热 作简便，准确度高.此外，其他方法将反应器置于恒 装置温度并控制在0.01K以内，气体流量控制导致 温环境中4，溶液向外散热量大，几乎所有文献 的不确定度为0.35%，外壳同步升温精度控制在 均未给出溶液发生反应的升温过程分析.本文采用 3%以内，气体传热误差为1%.考虑到质量、温度、 工业中广泛使用的MEA(伯胺)溶液进行实验，而 压力、气体热容及量热误差等方面，焓不确定度估计 MDEA(叔胺)具有黏性大、反应过程中不形成氨基 为10%， 甲酸酯等与MEA不同的性质，将两者进行对比，分 (1)液体比热容测定装置 析其吸收反应热及升温过程 如图2，反应器外面有加热丝包裹，加热丝外有 保温层，使加热丝能够快速响应升温，两者共同构成 1 方法原理与实验装置 同步热跟踪的外壳，外壳与反应器之间为真空层 1.1实验材料 温度变送器连接热电偶，将溶液及外壳温度数据传 实验材料主要有：去离子水；乙醇胺(MEA),西 输给数采控制器.上位机程序控制数采控制器， 陇化工股份有限公司，分析纯：N甲基二乙醇胺 OUT2控制反应器内加热体，使溶液升温或是在保 (MDEA),国药集团化学试剂有限公司，化学纯 温状态，功率表测得其加热功率.OUT1控制外壳对 1.2方法原理 反应器同步热跟踪.反应器内搅拌使溶液温度场 气液反应热大部分被溶液及盛装容器吸收，从 均匀 而提升系统的温度，随着温度梯度的加大，热量以传 实验过程中，将溶液注入反应器后，统一控制溶 导、对流与辐射方式散至环境中，温度梯度越大，热 液及外壳稳定在实验起始温度，实验从起始温度开

张宣凯等: 基于同步热跟踪法的微量气液反应热测量 在化学吸收法捕集 CO2的过程中，再生能耗占 CO2捕集过程总能耗的 70% ～ 80%［1］，目前最关键 的问题是如何降低再生能耗［2--7］，因此，气液反应热 的测量对筛选低能耗吸收剂具有重大意义［8--15］． 微反应量热计( Thermal Hazard Technology 公司 提供) ［16］、DＲC 系统量热计［17］、Setaram C80 导热差 示量热仪［18--20］等设备，均采用参考溶液与样品溶液 同时加热进行对比的方式来抵消样品溶液在测量中 的热损失，但不能直接测量反应热，存在较大误差． Kim 和 Svendsen［21］用容积为 2000 mL 的机械搅 拌反应器( 瑞典 ChemiSens AB，型号 CPA122) 进行 量热，无保温措施; 郭东方等［22］使用 250 mL 真实热 流反应器( 型号 CPA201) ，将反应器置于恒温液体 中，器壁采用双层玻璃形成中间真空层，虽然采用了 保温措施，但仍无法阻止热散失． 上述两种方法均 采用增加实验样品量的方式，使散热损失占总放热 量的比例变小，可近似忽略，不适用于微量气液反应 热的测量． Svensson 等［23］ 使用化学过程分析仪 ( 型 号 CPA202) ，将容量为 250 mL 的反应器置于恒温液体 浴中，使其主动随着反应器升温而升温，但对于微量 样品的情况，恒温液体浴的升温过程响应驰豫时间 长，控制精度低，且反应器顶部和底部还可进行热传 导，实验误差大． 本文提出一种同步热跟踪法测微量气液反应 热，程序控制外部热屏障自动与反应器温度保持同 步，响应快速，无需热补偿，可以减小上述方法所带 来的误差，适用于微量气液反应热的测量，计算与操 作简便，准确度高． 此外，其他方法将反应器置于恒 温环境中［24--25］，溶液向外散热量大，几乎所有文献 均未给出溶液发生反应的升温过程分析． 本文采用 工业中广泛使用的 MEA( 伯胺) 溶液进行实验，而 MDEA( 叔胺) 具有黏性大、反应过程中不形成氨基 甲酸酯等与 MEA 不同的性质，将两者进行对比，分 析其吸收反应热及升温过程． 1 方法原理与实验装置 1. 1 实验材料 实验材料主要有: 去离子水; 乙醇胺( MEA) ，西 陇化工 股 份 有 限 公 司，分 析 纯; N-甲 基 二 乙 醇 胺 ( MDEA) ，国药集团化学试剂有限公司，化学纯． 1. 2 方法原理 气液反应热大部分被溶液及盛装容器吸收，从 而提升系统的温度，随着温度梯度的加大，热量以传 导、对流与辐射方式散至环境中，温度梯度越大，热 散失越严重． 在反应容器的外部设置同步升温壁 垒，与内部溶液形成小的温差，构成一道“热屏障”， 由于温度梯度减小，热散失驱动力不足，从而减少热 损失． 传统的温度控制方法多采用 PID 控制系统根 据温度设定值进行控制［26--28］，而热屏温度的控制是 追踪变化的设定点，并且不能采用“波浪形”温控曲 线，如果热屏温度高于内部，会导致向系统输入热 量，造成反应热增大，因此，热屏温度应“尾随”并且 最大限度的“逼近”内部溶液的温度． 本研究采用固 态继电器替代机械继电器缩短加热电路通断时间; 采用上位机自由编程读取溶液温度和控制热屏温 升，从而实现热屏温度与内部溶液温度的同步，造成 “零”温度梯度． 具体实验流程如图 1 所示，c系统 与 m系统表示反应器系统各项的比热容与质量，Qabs为反 应过程放热量． 1. 3 实验装置 所有实验仪器主要有: 温度变送器，AI-702，厦 门宇电自动化科技有限公司，测量精度 0. 2 级; 单相 电参数盘面表，DSW-9002C，创鸿( 台湾) 电子仪表 有限公司，测量精度 0. 4 级，0. 02 W; 数采控制器， JY-DAM0808，北京聚英翱翔电子有限公司; 进气数 字式质量流量控制器，CS200，北京七星华创流量计 有限公司，测量精度 0. 35 级; 出气质量流量传感器， FS4001-200 mL，美国矽翔微机电系统有限公司，测 量精度 1. 5 级; 磁力加热搅拌器，79-1，上海梅香仪 器有限公司; 集热式搅拌器，DF-101S，金坛市科析 仪器有限公司． 测量质量精度为 0. 03% ，温度变送器设置量热 装置温度并控制在 0. 01 K 以内，气体流量控制导致 的不确定度为 0. 35% ，外壳同步升温精度控制在 3% 以内，气体传热误差为 1% ． 考虑到质量、温度、 压力、气体热容及量热误差等方面，焓不确定度估计 为 10% ． ( 1) 液体比热容测定装置． 如图 2，反应器外面有加热丝包裹，加热丝外有 保温层，使加热丝能够快速响应升温，两者共同构成 同步热跟踪的外壳，外壳与反应器之间为真空层． 温度变送器连接热电偶，将溶液及外壳温度数据传 输给数采控制器． 上位机程序控制数采控制器， OUT2 控制反应器内加热体，使溶液升温或是在保 温状态，功率表测得其加热功率． OUT1 控制外壳对 反应器同步热跟踪． 反应器内搅拌使溶液温度场 均匀． 实验过程中，将溶液注入反应器后，统一控制溶 液及外壳稳定在实验起始温度，实验从起始温度开 · 963 ·

·370· 工程科学学报，第41卷，第3期 测定气液反应热△， 测定溶液比热容©滚 1.给定浓度配液，称量质量m液 1相同浓度，质量m黄配液 2.设置进气流量Gm 2.设置输入电压V,测量得电流/ 业 3.外壳同步热跟踪 3,外壳同步热跟踪 4,进气开始，计时，=0，起始温度T, 相同 测量出气流量G。 4加热体通电，计时=0，起始温度T 5.吸收饱和，实验结束，升温区间内温差 相同 5.实验结束，选取相同温度T,△T=T-T=-△T △T=-T-T),到达峰值温度T,对应时间= T T对应时间='，Q=PG-D=P' 6Q=6ea"et+∑egnA7 6A=+∑end Q Q Q cc点c6 得到△T内的平均c落 图1气液反应热测量实验流程图 Fig.1 Flow chart of gas-iquid reaction heat measurement 3S48 19 20 F-1015g C3200 MFC flo =【5232= OOO·o000000ooooooOO 1234 JUYING +02 ELECTRONIC cin port:Rs232Rs485RF 23 C01 MO品E UDIAN ■et export■W日■USB 6品 A-7021 p0 werRS2321BS485广口DI■AH7T口D=AI7 02 6o88868556666660o0 24 G0。山。宽界D 910a 91o1四1包间141的a ooaoy aood 3 1314■6 ▣234567□8 :*a IN:I POWER AC 90-250 D3W.9902 1一C0,气瓶：2一减压阀：3一进气质量流量传感器：4一集热式搅拌器：5一进气管：6一温度变送器：7、10一热电偶：8一球形网状材料：9一保 温材料：11一加热丝：12、14一玻璃容器：13一加热体：15一蜂窝状多孔材料：16一磁力搅拌子：17一支撑圆锥体：18一磁力搅拌器：19一出气 管：20一干燥管：21一出气质量流量传感器：22一上位计算机：23一数采控制器：24一功率表 图2气液反应热及液体比热容测量装置图 Fig.2 Schematic of gas-iquid reaction heat and liquid specific heat capacity measurement device 始加热，并记录加热功率及温度变化. 一起始温度，加热体停止加热，开始通气，反应开始 (2)气液反应热测定装置. 进行：当进气量与出气量相等时，吸收达到饱和，实 如图2，实验采用质量流量传感器控制进气流 验结束 量保持不变，压力保持恒定，Mathonat等四也验证 2数据处理 过压力与反应热的关系相互独立，因此本文中不考 虑压力的影响.进气管路经过水浴加热，使其保持 2.1数据处理案例 在反应温度 (1)液体比热容计算. 在实验过程中，统一控制溶液及外壳稳定在同 为计算吸收反应升温区间内的平均比热容

工程科学学报，第 41 卷，第 3 期 图 1 气液反应热测量实验流程图 Fig． 1 Flow chart of gas--liquid reaction heat measurement 1—CO2气瓶; 2—减压阀; 3—进气质量流量传感器; 4—集热式搅拌器; 5—进气管; 6—温度变送器; 7、10—热电偶; 8—球形网状材料; 9—保 温材料; 11—加热丝; 12、14—玻璃容器; 13—加热体; 15—蜂窝状多孔材料; 16—磁力搅拌子; 17—支撑圆锥体; 18—磁力搅拌器; 19—出气 管; 20—干燥管; 21—出气质量流量传感器; 22—上位计算机; 23—数采控制器; 24—功率表 图 2 气液反应热及液体比热容测量装置图 Fig． 2 Schematic of gas--liquid reaction heat and liquid specific heat capacity measurement device 始加热，并记录加热功率及温度变化． ( 2) 气液反应热测定装置． 如图 2，实验采用质量流量传感器控制进气流 量保持不变，压力保持恒定，Mathonat 等［29］也验证 过压力与反应热的关系相互独立，因此本文中不考 虑压力的影响． 进气管路经过水浴加热，使其保持 在反应温度． 在实验过程中，统一控制溶液及外壳稳定在同 一起始温度，加热体停止加热，开始通气，反应开始 进行; 当进气量与出气量相等时，吸收达到饱和，实 验结束． 2 数据处理 2. 1 数据处理案例 ( 1) 液体比热容计算． 为计算吸收反 应 升 温 区 间 内 的 平 均 比 热 容 · 073 ·

张宣凯等：基于同步热跟踪法的微量气液反应热测量 ·371· C漆液，选取时间一温度T曲线中从起始温度T,到升 温区间的峰值温度T,区间段，如图3(a)所示，以质 是=(m精+∑c) (1) 量分数30%的MEA溶液为例，根据吸收反应实验 其中，加热量Q只包含对溶液及反应器从T,升温到 中的升温区间（从T到T),对应到测定溶液比热容 T,的加热量，由于外壳热屏的同步热跟踪效应，与环 实验中的升温曲线，选取同样的温度区间，如图 境热交换产生的散热量可近似忽略.因此根据公式 3(b)所示，将选取的测定溶液比热容实验中的升温 (1)即可得到溶液在升温区间△T内的平均比热容 曲线（从T,到T2)对应到同时间段（从t到t)的功 c溶液，其中△T=T2一T1,m溶液为溶液质量，C系统与 率曲线，将此段功率P对时间积分得到此温度区 m系统表示反应器系统各项的比热容与质量，包括玻 间内的加热量Q 璃容器、加热体、搅拌子等 10 70 10 60 40 ■一T测比热容实验) 30 -▲T气液反应实验) 20 20 -20002004006008001000120014001600 0100200300,400500g600700 e 图315mL质量分数30%的MEA溶液的相关实验曲线.()测定溶液比热容实验中的升温曲线与气液反应过程中的升温曲线：(b)在 测定溶液比热容的实验中的升温曲线与加热功率曲线 Fig.3 Experimental curves of 15 ml30%MEA solution:(a)the rising curves of specific heat capacity and gas-iquid reaction measurement:(b) the rising temperature curve and heating power curve for the determination of solution specific heat capacity (2)气液反应热计算. 为每次测量的反应热值，x为5次测量比热容的平 根据公式(2)即可得到反应过程放热量Q,仅 均值，则置信区间为：(x-t(0.05,4)S,,x+t(0.05, 供给反应器及溶液升温，实验测得本实验过程中气 4)S,),其中S为平均值的标准误差，计算结果如 体带走的热量微小，其他形式热损失可忽略不计. 表1所示 其中△T=T,-T2,c系统与m系统与公式(1)相同. 公式(3)中Gm与Gn分别为出气与进气瞬时流 (x-2 量，反应起始时间为，峰值温度T,所对应的时间 S:=± W5×(5-1) (4) 为2，从t,到2对时间积分，得到出气流量与进气流 2.3气液反应热实验值与文献值对比 量，进气量与出气量的差值即为升温区间内的吸收 如图4所示，Mathonat等给出了反应温度T 气量，放热量Q与吸收气量的比值即为反应热 为40、80和120℃时，质量分数为30%的MEA溶液 △Hhs 吸收CO,趋近于饱和时的反应热.根据这些反应热 Q山=（C溶液m溶液+∑c系铁m系统△T (2) 值拟合得出直线y=0.263x+70.271,得到20℃时 的反应热拟合外延值为-75.49kJ·mol-1.本实验 △Hhs= Q (3) fc.-cd 在反应温度为20℃条件下的反应热测得值为 -(73.12±4.06)kJ.mol-,在实验误差范围内. 2.2MEA、MDEA气液反应热计算结果 3实验结果与分析 每种溶液实验重复测量5次，由于样本容量较 小（实验次数为5，自由度为5-1=4），误差分析选 3.1同步热跟踪与非同步热跟踪实验结果对比 用t分布.选取置信概率为95%（显著性水平α为 15mL质量分数为30%的MEA溶液，在反应温 0.05),计算置信区间.查t分布表得1(0.05,4)= 度为20℃无同步热跟踪（以下简称“非热跟踪”） 2.77,根据公式(4)得到样本均值的标准差.其中x: 时，进行吸收CO2气液反应实验，测得反应热为

张宣凯等: 基于同步热跟踪法的微量气液反应热测量 c溶液，选取时间 t--温度 T 曲线中从起始温度 T1到升 温区间的峰值温度 T2区间段，如图 3( a) 所示，以质 量分数 30% 的 MEA 溶液为例，根据吸收反应实验 中的升温区间( 从 T1到 T2 ) ，对应到测定溶液比热容 实验中的升温曲线，选取同样的温度区间，如 图 3( b) 所示，将选取的测定溶液比热容实验中的升温 曲线( 从 T1到 T2 ) 对应到同时间段( 从 t'1到 t'2 ) 的功 率曲线，将此段功率 P 对时间 t'积分得到此温度区 间内的加热量 Q． Q ΔT' = ( c溶液m溶液 + ∑ c系统m系统 ) ( 1) 其中，加热量 Q 只包含对溶液及反应器从 T1升温到 T2的加热量，由于外壳热屏的同步热跟踪效应，与环 境热交换产生的散热量可近似忽略． 因此根据公式 ( 1) 即可得到溶液在升温区间 ΔT'内的平均比热容 c溶液，其中 ΔT' = T2 － T1，m溶液 为溶液质量，c系统 与 m系统表示反应器系统各项的比热容与质量，包括玻 璃容器、加热体、搅拌子等． 图 3 15 mL 质量分数 30% 的 MEA 溶液的相关实验曲线． ( a) 测定溶液比热容实验中的升温曲线与气液反应过程中的升温曲线; ( b) 在 测定溶液比热容的实验中的升温曲线与加热功率曲线 Fig． 3 Experimental curves of 15 mL 30% MEA solution: ( a) the rising curves of specific heat capacity and gas--liquid reaction measurement; ( b) the rising temperature curve and heating power curve for the determination of solution specific heat capacity ( 2) 气液反应热计算． 根据公式( 2) 即可得到反应过程放热量 Qabs，仅 供给反应器及溶液升温，实验测得本实验过程中气 体带走的热量微小，其他形式热损失可忽略不计． 其中 ΔT = T1 － T2，c系统与 m系统与公式( 1) 相同． 公式( 3) 中 Gout与 Gin分别为出气与进气瞬时流 量，反应起始时间为 t1，峰值温度 T2 所对应的时间 为 t2，从 t1到 t2对时间积分，得到出气流量与进气流 量，进气量与出气量的差值即为升温区间内的吸收 气量，放热量 Qabs 与吸收气量的比值即为反应热 ΔHabs． Qabs = ( c溶液m溶液 + ∑ c系统m系统 )·ΔT ( 2) ΔHabs = Qabs ∫ 2t t1 Gin dt － ∫ 2t t1 Goutdt ( 3) 2. 2 MEA、MDEA 气液反应热计算结果 每种溶液实验重复测量 5 次，由于样本容量较 小( 实验次数为 5，自由度为 5 － 1 = 4) ，误差分析选 用 t 分布． 选取置信概率为 95% ( 显著性水平 α 为 0. 05) ，计算置信区间． 查 t 分布表得 t( 0. 05，4) = 2. 77，根据公式( 4) 得到样本均值的标准差． 其中 xi 为每次测量的反应热值，x 为 5 次测量比热容的平 均值，则置信区间为: ( x － t( 0. 05，4) S x － ，x + t( 0. 05， 4) S x － ) ，其中 S x － 为平均值的标准误差，计算结果如 表 1 所示． S x － = ± ∑ 5 i = 1 ( xi － x) 2 槡5 × ( 5 － 1) ( 4) 2. 3 气液反应热实验值与文献值对比 如图 4 所示，Mathonat 等［29］给出了反应温度 T 为 40、80 和 120 ℃时，质量分数为 30% 的 MEA 溶液 吸收 CO2趋近于饱和时的反应热． 根据这些反应热 值拟合得出直线 y = 0. 263x + 70. 271，得到 20 ℃ 时 的反应热拟合外延值为 － 75. 49 kJ·mol － 1 ． 本实验 在反 应 温 度 为 20 ℃ 条件下的反应热测得值为 － ( 73. 12 ± 4. 06) kJ·mol － 1，在实验误差范围内． 3 实验结果与分析 3. 1 同步热跟踪与非同步热跟踪实验结果对比 15 mL 质量分数为 30% 的 MEA 溶液，在反应温 度为 20 ℃ 无同步热跟踪( 以下简称“非热跟踪”) 时，进 行 吸 收 CO2 气液 反 应 实 验，测 得 反 应 热 为 · 173 ·

·372 工程科学学报，第41卷，第3期 表1不同质量分数醇胺溶液比热容及反应热计算结果 Table 1 Values of specific heat capacity and reaction heat for different concentrations of alcohol amine solution 溶液 质量 温度 比热容/ 比热容置信区间/ -△He/ 反应热置信区间/ 种类 分数/% 区间/℃ (kJkg1℃1) (kJkg.℃) (kJ-mol-1) (kJ-mol-1) 10 19.5-20.5 3.33 (3.33-0.09,3.33+0.09) 90.46 (90.46-13.45,90.46+13.45) 20 19.5-20.5 3.61 (3.61-0.10,3.61+0.10) 80.16 (80.26-5,91,80.26+5.91) MEA 30 19.5~20.5 3.72 (3.72-0.10,3.72+0.10) 73.12 (73.12-4.06,73.12+4.06) 40 19.5-20.5 3.87 (3.87-0.10,3.87+0.10) 73.42 (73.42-8.18,73.42+8.18) 50 19.520.5 3.95 (3.95-0.14,3.95+0.14) 77.60 (77.60-9.23,77.60+9.23) 10 19.5-20.5 3.76 (3.76-0.16,3.76+0.16) 40.62 (40.62-12.53,40.62-12.53) 20 19.5-20.5 3.81 (3.81-0.16,3.81+0.16) 30.45 (30.45-8.94,30.45+8.94) MDEA 30 19.5-20.5 3.60 (3.60-0.13,3.60+0.13) 29.30 (29.30-8.65,29.30+8.65) 40 19.5-20.5 3.34 (3.34-0.18,3.34+0.18) 26.59 (26.59-9.17,26.59+9.17) 50 19.5~20.5 3.38 (3.38-0.20,3.38+0.20) 14.92 (14.92-11.28,14.92+11.28) 环境散失的热量越多，反应放热量也越大，对应需 105 ▲文献值 ●拟合值 要更多的CO,参加反应，与温度最高点对应的出气 100 ★实验值 一拟合直线 流量曲线中C02出气量越少，表现为温度曲线最高 95 点左移.热跟踪时，反应放热大部分用于溶液的升 90 温，仅有微小部分散失到环境中，因此只需少量的 85 反应热，此时对应吸收剂耗尽状态，CO2出气量接 80 近最大值.因此热跟踪温度曲线的最高点更接近 75 于出气流量曲线中C0,的最大值，表现为温度曲线 最高点向右移动.由于本实验装置没有完全实现 20 0 6080100120 反应器的“立体”热屏障，存在小部分散热，因此温 度最高点出现在吸收接近饱和的位置. 图4质量分数30%的MEA溶液吸收C02时，不同温度T下的 非热跟踪实验中，当时间超过温度最大值对应 反应热文献值、拟合值与实验值对比 Fig.4 Literature value,fitted value,and experimental value of dif- 的时刻时，反应放热不足以维持溶液热损失，此时吸 ferent temperatures of 30%MEA solution with absorbed CO, 收剂逐渐耗尽，但仍有部分未参与反应，可继续吸收 C02,此时吸收的C02未能考虑到吸收热的计算中， -42.25 kJmol-1.其他条件相同，有同步热跟踪（以 造成测量误差.热跟踪时，溶液的温度最高点升高， 下简称“热跟踪”)时测得反应热为-(73.12± 升温时间变长，升温阶段的吸收气量增多，即 4.06)kJ.mol-1,绝对值远大于非热跟踪时的实验结 图5(b)中的虚线区域面积与图5(a)相比变大.无 果，更接近Mathonat等网给出的文献值 论有无同步热跟踪，耗尽相同量吸收剂所需的C0, 如图5(a)所示，非热跟踪实验中，当反应温度 气量应该相同，热跟踪过程中，由于出气流量曲线中 T达到最大值时，参加反应的C0,开始出现剩余， CO2开始逸出点右移，出气流量曲线平台段对应的 表现为出气流量曲线呈现“抬头”趋势：但热跟踪 C0,吸收气量增大，因此在出气流量开始上升后的 测定时，如图5(b)所示，反应温度T最大值对应 阶段，C0,吸收量则会减小.所以，热跟踪时出气流 溶液吸收近饱和状态，此时，参加反应的C02剩余 量上升阶段曲线的上升速率更快，使温度达到最高 量迅速增加，出气流量曲线陡然上扬.两者对比， 点后没有计算在内的剩余吸收气量更小，可近似忽 当C0,开始出现过剩时（对应出气流量曲线“抬 略.综上所述，热跟踪时，升温区间内的吸收气量更 头”点)，热跟踪时温度曲线最高点向右发生了移 接近吸收饱和时的气量，此时吸收气量的计算准确 动.其主要原因如下：非热跟踪时，反应产生的热 度更高 一部分用于溶液升温，另一部分散失到环境中.当 如图5(C)所示，热跟踪时，温度T曲线的峰值 反应热与热损失相同时温度曲线出现最高点，向 远高于非热跟踪，达到温度最高点所需的时间更长

工程科学学报，第 41 卷，第 3 期 表 1 不同质量分数醇胺溶液比热容及反应热计算结果 Table 1 Values of specific heat capacity and reaction heat for different concentrations of alcohol amine solution 溶液 种类 质量 分数/% 温度 区间/℃ 比热容/ ( kJ·kg － 1·℃ － 1 ) 比热容置信区间/ ( kJ·kg － 1·℃ － 1 ) － ΔHabs / ( kJ·mol － 1 ) 反应热置信区间/ ( kJ·mol － 1 ) 10 19. 5 ～ 20. 5 3. 33 ( 3. 33 － 0. 09，3. 33 + 0. 09) 90. 46 ( 90. 46 － 13. 45，90. 46 + 13. 45) 20 19. 5 ～ 20. 5 3. 61 ( 3. 61 － 0. 10，3. 61 + 0. 10) 80. 16 ( 80. 26 － 5，91，80. 26 + 5. 91) MEA 30 19. 5 ～ 20. 5 3. 72 ( 3. 72 － 0. 10，3. 72 + 0. 10) 73. 12 ( 73. 12 － 4. 06，73. 12 + 4. 06) 40 19. 5 ～ 20. 5 3. 87 ( 3. 87 － 0. 10，3. 87 + 0. 10) 73. 42 ( 73. 42 － 8. 18，73. 42 + 8. 18) 50 19. 5 ～ 20. 5 3. 95 ( 3. 95 － 0. 14，3. 95 + 0. 14) 77. 60 ( 77. 60 － 9. 23，77. 60 + 9. 23) 10 19. 5 ～ 20. 5 3. 76 ( 3. 76 － 0. 16，3. 76 + 0. 16) 40. 62 ( 40. 62 － 12. 53，40. 62 － 12. 53) 20 19. 5 ～ 20. 5 3. 81 ( 3. 81 － 0. 16，3. 81 + 0. 16) 30. 45 ( 30. 45 － 8. 94，30. 45 + 8. 94) MDEA 30 19. 5 ～ 20. 5 3. 60 ( 3. 60 － 0. 13，3. 60 + 0. 13) 29. 30 ( 29. 30 － 8. 65，29. 30 + 8. 65) 40 19. 5 ～ 20. 5 3. 34 ( 3. 34 － 0. 18，3. 34 + 0. 18) 26. 59 ( 26. 59 － 9. 17，26. 59 + 9. 17) 50 19. 5 ～ 20. 5 3. 38 ( 3. 38 － 0. 20，3. 38 + 0. 20) 14. 92 ( 14. 92 － 11. 28，14. 92 + 11. 28) 图 4 质量分数 30% 的 MEA 溶液吸收 CO2 时，不同温度 T 下的 反应热文献值、拟合值与实验值对比 Fig． 4 Literature value，fitted value，and experimental value of dif￾ferent temperatures of 30% MEA solution with absorbed CO2 － 42. 25 kJ·mol － 1 ． 其他条件相同，有同步热跟踪( 以 下简 称“热 跟 踪”) 时测得反应热为 － ( 73. 12 ± 4. 06) kJ·mol － 1，绝对值远大于非热跟踪时的实验结 果，更接近 Mathonat 等［29］给出的文献值． 如图 5( a) 所示，非热跟踪实验中，当反应温度 T 达到最大值时，参加反应的 CO2开始出现剩余， 表现为出气流量曲线呈现“抬头”趋势; 但热跟踪 测定时，如图 5( b) 所示，反应温度 T 最大值对应 溶液吸收近饱和状态，此时，参加反应的 CO2剩余 量迅速增加，出气流量曲线陡然上扬． 两者对比， 当 CO2开始出 现 过 剩 时( 对应出气流量曲线“抬 头”点) ，热跟踪时温度曲线最高点向右发生了移 动． 其主要原因如下: 非热跟踪时，反应产生的热 一部分用于溶液升温，另一部分散失到环境中． 当 反应热与热损失相同时温度曲线出现最高点，向 环境散失的热量越多，反应放热量也越大，对应需 要更多的 CO2参加反应，与温度最高点对应的出气 流量曲线中 CO2出气量越少，表现为温度曲线最高 点左移． 热跟踪时，反应放热大部分用于溶液的升 温，仅有微小部分散失到环境中，因此只需少量的 反应热，此时对应吸收剂耗尽状态，CO2 出气量接 近最大值． 因此热跟踪温度曲线的最高点更接近 于出气流量曲线中 CO2的最大值，表现为温度曲线 最高点向右移动． 由于本实验装置没有完全实现 反应器的“立体”热屏障，存在小部分散热，因此温 度最高点出现在吸收接近饱和的位置． 非热跟踪实验中，当时间超过温度最大值对应 的时刻时，反应放热不足以维持溶液热损失，此时吸 收剂逐渐耗尽，但仍有部分未参与反应，可继续吸收 CO2，此时吸收的 CO2未能考虑到吸收热的计算中， 造成测量误差． 热跟踪时，溶液的温度最高点升高， 升温 时 间 变 长，升温阶段的吸收气量增多，即 图 5( b) 中的虚线区域面积与图 5( a) 相比变大． 无 论有无同步热跟踪，耗尽相同量吸收剂所需的 CO2 气量应该相同，热跟踪过程中，由于出气流量曲线中 CO2开始逸出点右移，出气流量曲线平台段对应的 CO2吸收气量增大，因此在出气流量开始上升后的 阶段，CO2吸收量则会减小． 所以，热跟踪时出气流 量上升阶段曲线的上升速率更快，使温度达到最高 点后没有计算在内的剩余吸收气量更小，可近似忽 略． 综上所述，热跟踪时，升温区间内的吸收气量更 接近吸收饱和时的气量，此时吸收气量的计算准确 度更高． 如图 5( c) 所示，热跟踪时，温度 T 曲线的峰值 远高于非热跟踪，达到温度最高点所需的时间更长， · 273 ·

张宣凯等：基于同步热跟踪法的微量气液反应热测量 ·373· 温度上升阶段曲线形状也更趋近于无散热损失情况 算的升温区间更接近完全绝热时的升温区间，利用 下的直线。综上表明，热跟踪能有效阻止溶液向环 公式(2)计算的反应放热量精度更高，使反应热测 境散热，反应热更多为溶液及反应器所吸收，此时计 量值更为准确 36@ 60 40 b▲ 32 30 30 28 20 20 10 10 溶液温度 鱼一进气流量 9 -●一出气流量 0 -200020040060080010001200140016001800 -200020040060080010001200140016001800 1/s 60m 55 50 35 30 25 ■一同步热跟踪 20 ▲一非同步热跟踪 500 1000 1500 2000 图515mL质量分数30%的MEA溶液温度T曲线及进气、出气流量曲线.(a)非同步热跟踪：(b)同步热跟踪：()温度T曲线对比 Fig.5 Temperature curves,inlet and outlet flow curves for 15 ml30%MEA solution:(a)non-synchronous thermal tracking:(b)synchronous thermal tracking:(e)comparison of the temperature curves 3.2MEA、MDEA不同浓度反应热对比 Kim和Svendsen L研究了不同种类胺溶液在 100 最一MEA MDEA 不同浓度下吸收CO2的反应热，指出胺浓度对反应 热的影响不显著，但并未给出具体范围及反应热随 浓度变化的关系.本实验针对伯胺(MEA)与叔胺 (MDEA)两种胺溶液浓度与反应热之间的关系进行 HV- 40 研究. 如图6所示，MEA与MDEA两种溶液在质量分 数为20%~40%时，反应热变化较缓，认为反应热 随胺溶液浓度变化不显著.吸收热包含物理溶解 20 20 40 热、化学反应热两部分，Svensson等四指出物理溶 质量分数/% 图6不同质量分数MEA、MDEA溶液反应热对比 解的CO,不与醇胺溶液反应，采用亨利定律计算的 Fig.6 Reaction heats at different mass fractions of MEA and MDEA CO,物理溶解量低于总吸收量的2%，可忽略不计. solutions 式(5)、(6)分别为MEA(其分子式简单表示为 RRNH)、MDEA溶液吸收CO,的反应方程式，对于 等B利用vant Hoff公式(7)计算得出MEA溶液反 MEA溶液，化学反应热主要包括质子化作用热、氨 应的温度效应主要由氨基甲酸酯形成热决定，其中 基甲酸酯形成热等，而MDEA没有氨基甲酸酯的形 K是化学反应平衡常数，R是普适气体常数，T是反 成，因此其反应焓不包括氨基甲酸酯形成热.Kim 应温度，下角标p表示压力，△H是反应热

张宣凯等: 基于同步热跟踪法的微量气液反应热测量 温度上升阶段曲线形状也更趋近于无散热损失情况 下的直线． 综上表明，热跟踪能有效阻止溶液向环 境散热，反应热更多为溶液及反应器所吸收，此时计 算的升温区间更接近完全绝热时的升温区间，利用 公式( 2) 计算的反应放热量精度更高，使反应热测 量值更为准确． 图 5 15 mL 质量分数 30% 的 MEA 溶液温度 T 曲线及进气、出气流量曲线． ( a) 非同步热跟踪; ( b) 同步热跟踪; ( c) 温度 T 曲线对比 Fig． 5 Temperature curves，inlet and outlet flow curves for 15 mL 30% MEA solution: ( a) non-synchronous thermal tracking; ( b) synchronous thermal tracking; ( c) comparison of the temperature curves 3. 2 MEA、MDEA 不同浓度反应热对比 Kim 和 Svendsen［30］研究了不同种类胺溶液在 不同浓度下吸收 CO2的反应热，指出胺浓度对反应 热的影响不显著，但并未给出具体范围及反应热随 浓度变化的关系． 本实验针对伯胺( MEA) 与叔胺 ( MDEA) 两种胺溶液浓度与反应热之间的关系进行 研究． 如图 6 所示，MEA 与 MDEA 两种溶液在质量分 数为 20% ～ 40% 时，反应热变化较缓，认为反应热 随胺溶液浓度变化不显著． 吸收热包含物理溶解 热、化学反应热两部分，Svensson 等［23］指出物理溶 解的 CO2不与醇胺溶液反应，采用亨利定律计算的 CO2物理溶解量低于总吸收量的 2% ，可忽略不计． 式( 5) 、( 6) 分别为 MEA( 其分子式简单表示为 ＲＲ'NH) 、MDEA 溶液吸收 CO2的反应方程式，对于 MEA 溶液，化学反应热主要包括质子化作用热、氨 基甲酸酯形成热等，而 MDEA 没有氨基甲酸酯的形 成，因此其反应焓不包括氨基甲酸酯形成热． Kim 图 6 不同质量分数 MEA、MDEA 溶液反应热对比 Fig． 6 Ｒeaction heats at different mass fractions of MEA and MDEA solutions 等［31］利用 van't Hoff 公式( 7) 计算得出 MEA 溶液反 应的温度效应主要由氨基甲酸酯形成热决定，其中 K 是化学反应平衡常数，Ｒ 是普适气体常数，T 是反 应温 度，下 角 标 p 表 示 压 力，ΔHabs 是 反 应 热． · 373 ·

·374 工程科学学报，第41卷，第3期 Mathonat等也得出结论，形成氨基甲酸酯的伯胺 3.3MEA、MDEA反应过程升温曲线对比 比不形成氨基甲酸酯的叔胺产生的吸收热多，但只 15mL质量分数为30%的MEA与MDEA溶液 对比了质量分数为30%时MEA与MDEA的反应 进行吸收CO2反应过程中反应温度T曲线如图7所 热.图6为反应温度为20℃时，MEA、MDEA吸收 示.由于MDEA的反应热小于MEA,MDEA溶液的 CO2的质量分数与反应热的关系曲线图.从图中可 温升也远小于MEA.MDEA吸收CO2的反应更缓 以看出，MEA在各个质量分数对应的反应热均比 慢，到达温度最高点需要更长的时间，主要原因是 MDEA大，进一步验证了Mathonat等的结论在不同 MDEA溶液黏度远大于MEA溶液，没有反应的 质量分数下均成立 MDEA与CO2接触困难，反应的传质性能变差，不利 2RRNH +CO,RR'NH,+RR'NCOO (5) 于反应的进行，反应热释放缓慢 MDEA +CO2 +H2O=MDEAH*+HCO;(6) 55 △Hs (), (7) 50 MDEA只有质子化的过程，反应的进行受溶液 中H20含量的影响较大，因此随着浓度的升高，水 ,…直线轴助线 一质量分数30%MEA 含量减少，阻碍了MDEA质子化的进程，使反应更 35 0 400800.120016002000 难发生，反应热呈现降低的趋势 30 C02进入溶液后，首先，如式(8)所示，伯胺直接 一一直线辅助线 25 与C0,反应生成两性离子；两性离子一部分与胺分 量-质量分数30%MEA ▲质量分数30%MDEA 子反应生成氨基甲酸酯，如式(9)所示，另一部分与 水发生水解反应生成氨基甲酸酯，如式(10)所示. 100020003000400050006000 t的 反应热主要来自于上述中的式(9)与式(10)，而这 图715mL质量分数30%的MEA与MDEA溶液吸收C02反应 两种反应中前者没有水参加反应，后者有水参加反 过程中的升温曲线对比 应，因此，反应热的大小不仅取决于水含量的多少， Fig.7 Comparison of the heating curves of 15 mL.30%MEA/MDEA 即溶液浓度的高低。 solution in gas-liquid reaction RR'NH +CO,RR'NH*COO (8) 图7中右上角小图为MEA溶液的反应温度曲 RR'NH*COO+RR'NH 线，在低温区升温曲线接近红色虚线表示的直线辅 RRNH,+RRNCOO (9) 助线，溶液温度与环境的温差较小，向外部环境散热 RR'NH*COO+H,O-H,O*+RR'NCOO- 的驱动力不足，表现为40℃以下，散热很少，可忽略 (10) 不计：当温度大于40℃时，曲线开始偏离直线趋势， 如图6所示，MEA在30%左右呈现两种不同的 温度越高，偏离越大，直到温度达到最高点，这部分 趋势.低于30%时，MEA主要与CO2发生式(8)所 曲线表现为明显的“下凹”现象.在Mathonat等四 示反应，生成两性离子，并且MEA几乎不与两性离 给出的文献中，吸收热随负载量变化的曲线所得出 子反应，因此，生成的两性离子只能与H20反应生 的规律也是反应开始时为直线，反应一段时间后，直 成氨基甲酸酯，随着EA质量分数的增加，H,O含 线开始弯曲，而后反应热逐渐达到最高值，Km等在 量减少，式(9)生成氨基甲酸酯的放热反应被抑制， 文献21]中也提到同样的规律.由于实验装置没有 反应热降低：在30%~50%范围内，随着MEA的比 达到完全绝热，系统会向外部环境散热，因此，当温 例增加，过量的MEA与两性离子反应生成氨基甲酸 度到达最高点后，两者都没有保持最高温度，而是有 酯，此时存在两种放热反应，一种是两性离子与水的 不同程度的下降趋势. 反应，另一种是两性离子与胺分子的反应，MEA质 影响“下凹”现象的主要因素包括以下几点： 量分数越大，两性离子与胺分子的反应热也越大，两 (I)溶液温度T越高、与环境温差越大，散热量越 种反应的总体热效应比低于30%的放热大，表现为 多，而由于采用了热屏蔽技术，溶液与热屏间的温差 曲线上升.由于胺溶液浓度的增加降低了溶液中水 接近于零，从图7可以看出，40℃以下升温曲线基 的含量，导致式(10)对应的反应减弱，放热量有所 本维持直线状态，说明实验中可以很好的阻止溶液 降低，因此曲线上升较缓 热量的散失，因此，温差散热量很小；(2)随着溶液

工程科学学报，第 41 卷，第 3 期 Mathonat等［29］也得出结论，形成氨基甲酸酯的伯胺 比不形成氨基甲酸酯的叔胺产生的吸收热多，但只 对比了质量分数为 30% 时 MEA 与 MDEA 的反应 热． 图 6 为反应温度为 20 ℃ 时，MEA、MDEA 吸收 CO2的质量分数与反应热的关系曲线图． 从图中可 以看出，MEA 在各个质量分数对应的反应热均比 MDEA 大，进一步验证了 Mathonat 等的结论在不同 质量分数下均成立． 2ＲＲ'NH + CO2 幑幐ＲＲ'NH + 2 + ＲＲ'NCOO － ( 5) MDEA + CO2 + H2O MDEAH 幑幐 + + HCO － 3 ( 6 ( ) dln K d( 1 /T ) ) p = ΔHabs Ｒ ( 7) MDEA 只有质子化的过程，反应的进行受溶液 中 H2O 含量的影响较大，因此随着浓度的升高，水 含量减少，阻碍了 MDEA 质子化的进程，使反应更 难发生，反应热呈现降低的趋势． CO2进入溶液后，首先，如式( 8) 所示，伯胺直接 与 CO2反应生成两性离子; 两性离子一部分与胺分 子反应生成氨基甲酸酯，如式( 9) 所示，另一部分与 水发生水解反应生成氨基甲酸酯，如式( 10) 所示． 反应热主要来自于上述中的式( 9) 与式( 10) ，而这 两种反应中前者没有水参加反应，后者有水参加反 应，因此，反应热的大小不仅取决于水含量的多少， 即溶液浓度的高低． ＲＲ'NH + CO2 幑幐ＲＲ'NH + COO － ( 8) ＲＲ'NH + COO － + ＲＲ'NH 幑幐 ＲＲ'NH + 2 + ＲＲ'NCOO － ( 9) ＲＲ'NH + COO － + H2O H 幑幐 3O + + ＲＲ'NCOO － ( 10) 如图 6 所示，MEA 在 30% 左右呈现两种不同的 趋势． 低于 30% 时，MEA 主要与 CO2发生式( 8) 所 示反应，生成两性离子，并且 MEA 几乎不与两性离 子反应，因此，生成的两性离子只能与 H2O 反应生 成氨基甲酸酯，随着 MEA 质量分数的增加，H2O 含 量减少，式( 9) 生成氨基甲酸酯的放热反应被抑制， 反应热降低; 在 30% ～ 50% 范围内，随着 MEA 的比 例增加，过量的 MEA 与两性离子反应生成氨基甲酸 酯，此时存在两种放热反应，一种是两性离子与水的 反应，另一种是两性离子与胺分子的反应，MEA 质 量分数越大，两性离子与胺分子的反应热也越大，两 种反应的总体热效应比低于 30% 的放热大，表现为 曲线上升． 由于胺溶液浓度的增加降低了溶液中水 的含量，导致式( 10) 对应的反应减弱，放热量有所 降低，因此曲线上升较缓． 3. 3 MEA、MDEA 反应过程升温曲线对比 15 mL 质量分数为 30% 的 MEA 与 MDEA 溶液 进行吸收 CO2反应过程中反应温度 T 曲线如图 7 所 示． 由于 MDEA 的反应热小于 MEA，MDEA 溶液的 温升也远小于 MEA． MDEA 吸收 CO2 的反应更缓 慢，到达温度最高点需要更长的时间，主要原因是 MDEA 溶 液 黏 度 远 大 于 MEA 溶 液，没 有 反 应 的 MDEA 与 CO2接触困难，反应的传质性能变差，不利 于反应的进行，反应热释放缓慢． 图 7 15 mL 质量分数 30% 的 MEA 与 MDEA 溶液吸收 CO2 反应 过程中的升温曲线对比 Fig． 7 Comparison of the heating curves of 15 mL 30% MEA /MDEA solution in gas--liquid reaction 图 7 中右上角小图为 MEA 溶液的反应温度曲 线，在低温区升温曲线接近红色虚线表示的直线辅 助线，溶液温度与环境的温差较小，向外部环境散热 的驱动力不足，表现为 40 ℃以下，散热很少，可忽略 不计; 当温度大于 40 ℃时，曲线开始偏离直线趋势， 温度越高，偏离越大，直到温度达到最高点，这部分 曲线表现为明显的“下凹”现象． 在 Mathonat 等［29］ 给出的文献中，吸收热随负载量变化的曲线所得出 的规律也是反应开始时为直线，反应一段时间后，直 线开始弯曲，而后反应热逐渐达到最高值，Kim 等在 文献［21］中也提到同样的规律． 由于实验装置没有 达到完全绝热，系统会向外部环境散热，因此，当温 度到达最高点后，两者都没有保持最高温度，而是有 不同程度的下降趋势． 影响“下凹”现象的主要因素包括以下几点: ( 1) 溶液温度 T 越高、与环境温差越大，散热量越 多，而由于采用了热屏蔽技术，溶液与热屏间的温差 接近于零，从图 7 可以看出，40 ℃ 以下升温曲线基 本维持直线状态，说明实验中可以很好的阻止溶液 热量的散失，因此，温差散热量很小; ( 2) 随着溶液 · 473 ·

张宣凯等：基于同步热跟踪法的微量气液反应热测量 ·375· 温度升高，溶液的比热容发生变化，但溶液比热容变 研究.化工学报，2010,61(5)：1233) 化很小，不致造成大程度的弯曲效果；(3)Kim等 B]Zhang K F,Liu Z L,Wang YY,et al.Analysis and calculation 提出，随着温度的升高，C02的溶解越来越困难，更 of the desorption energy consumption of CO2 capture process by chemical absorption method.Chem Ind Eng Prog,2013,32 利于逆向的解吸反应进行，反应进行一段时间后，正 (12):3008 向吸收反应产物累积越来越多，由于化学平衡，抑制 (张克舫，刘中良，王远亚，等.化学吸收法CO2捕集解吸能 了正向吸收反应的发生，吸收反应逐渐变慢，放热量 耗的分析计算.化工进展，2013,32(12)：3008) 逐渐减少，升温速率越来越小，产生“下凹”现象，直 4 Ahmady A,Hashim MA,Aroua M K.Absorption of carbon diox- 至溶液吸收达到饱和，不再产生放热量. ide in the aqueous mixtures of methyldiethanolamine with three types of imidazolium-based ionic liquids.Fluid Phase Equilib 对比MEA与MDEA升温曲线，由于热跟踪系 2011,309(1):76 统制造了一个近似绝热的环境，在40℃以下，MEA [5]Zhang F,Ma J W,Zhou Z,et al.Study on the absorption of car- 溶液升温曲线呈直线上升，没有明显的散热，绝热效 bon dioxide in high concentrated MDEA and ILs solutions.Chem 果很好，MDEA升温曲线应如图7中虚线所示直线 Emg,2012,181H82:222 上升，在温度达到最高点后保持不变，但MDEA实 [6]Arcis H,Rodier L,Ballerat-Busserolles K,et al.Modeling of 际升温曲线在30℃左右仍有“下凹”现象，说明“下 (vapor liquid)equilibrium and enthalpy of solution of carbon di- oxide (CO2)in aqueous methyldiethanolamine (MDEA)solu- 凹”现象不是温差热损失引起的，而是由于反应达 tions.J Chem Thermodyn,2009,41(6):783 到近平衡状态时正向吸收反应放热量减少所致，因 ] Ci LE,Lundberg J,Pedersen M,et al.Measurements of CO2 ab- 此，MDEA升温曲线出现的“下凹”现象仍然归因于 sorption and heat consumption in laboratory rig.Energy Procedia, 自身产热量减少，在计算反应热时不应忽略，从而提 2014,63:1569 高反应热测量值的准确度. [8]QiG J.Wang S J,Yu H,et al.Prediction model of absorption heat for CO2 capture using aqueous ammonia.C/ESC J,2013,64 4结论 (9):3079 (齐国杰，王淑娟，YuHa,等.氨水吸收CO2的吸收热预测 (1)同步热跟踪法测微量气液反应热与文献中 模型.化工学报，2013,64(9)：3079) 数据对比，实验结果一致：与非同步热跟踪法相比， ] Liu B C,Shi C H,Li Q L,et al.Experimental investigation of 无论是吸收气量、升温区间的计算，还是最终得到的 performances and renewable energy consumption of absorbent DEA 反应热值，结果都更为准确，实验精度高 in CO2 capture system.Eng J Wuhan Unir,2012,45(6):817 (刘炳成，史澄辉，李庆领，等.电厂CO2捕集系统DEA化学 (2)随溶液质量分数的增加，MEA的反应热值 吸收剂性能与能耗实验研究.武汉大学学报（工学版），2012， 呈现先降低后升高的趋势，在质量分数为30%左右 45(6):81) 达到最低值，而MDEA的反应热值则一直处于下降 [10]Chen J,Luo W L,Li H.A review for research on thermodynam- 趋势.MEA、MDEA溶液质量分数在20%~40%时， ics and kinetics of carbon dioxide absorption with organic amines. 对反应热的影响不显著 CIESC J,2014,65(1):12 (陈健，罗伟亮，李晗.有机胺吸收二氧化碳的热力学和动 (3)胺溶液吸收CO2反应对应的升温曲线在开 力学研究进展.化工学报，2014,65(1)：12) 始阶段一直呈直线上升趋势，反应进行一段时间后， [11]Gupta M,Svendsen H F.Temperature dependent enthalpy of 反应产物不断累积，开始抑制正向吸收反应的进行， CO,absorption for amines and amino acids from theoretical calcu- 使曲线出现“下凹”，将该部分反应热考虑到计算中 lations at infinite dilution.Energy Procedia,2014,63:1106 可以提高结果的准确度. [12]Van Nierop E A,Hormoz S,House K Z,et al.Effect of absorp- tion enthalpy on temperature-swing CO2 separation process per- 参考文献 formance.Energy Procedia,2011,4:1783 03] Liu JZ,Wang SJ,Zhao B,et al.Absorption of carbon dioxide [Du M.Zhang L,Feng B.Effect of pH controlling method on ener- in aqueous ammonia.Energy Procedia,2009,1(1):933 gy consumption of CO2 desorption from rsch-solvent.J Chongqing [14]Liu JZ,Wang S J,Qi G J,et al.Kinetics and mass transfer of 0nim,2010,33(8):123 carbon dioxide absorption into aqucous ammonia.Energy Proce- (杜敏，张力，冯波.pH值调节剂对CO2吸收富液解吸能耗 dia,2011,4:525 的影响.重庆大学学报，2010,33(8)：123) [15]Srisang W,Pouryousefi F,Osei P A,et al.Evaluation of the Qin F,WangSJ,Svendsen HF,et al.Research on heat require- heat duty of catalyst-aided amine-based post combustion CO2 cap- ment of aqua ammonia regeneration for CO2 capture.C/ESC J, ture.Chem Eng Sci,2017,170:48 2010,61(5):1233 06]Abdulkadir A,Rayer A V,Quang D V,et al.Heat of absorption (秦锋，王淑娟，Svendsen H F,等.氨法脱碳系统再生能耗的 and specific heat of carbon dioxide in aqueous solutions of mono-

张宣凯等: 基于同步热跟踪法的微量气液反应热测量 温度升高，溶液的比热容发生变化，但溶液比热容变 化很小，不致造成大程度的弯曲效果; ( 3) Kim 等［24］ 提出，随着温度的升高，CO2 的溶解越来越困难，更 利于逆向的解吸反应进行，反应进行一段时间后，正 向吸收反应产物累积越来越多，由于化学平衡，抑制 了正向吸收反应的发生，吸收反应逐渐变慢，放热量 逐渐减少，升温速率越来越小，产生“下凹”现象，直 至溶液吸收达到饱和，不再产生放热量． 对比 MEA 与 MDEA 升温曲线，由于热跟踪系 统制造了一个近似绝热的环境，在 40 ℃ 以下，MEA 溶液升温曲线呈直线上升，没有明显的散热，绝热效 果很好，MDEA 升温曲线应如图 7 中虚线所示直线 上升，在温度达到最高点后保持不变，但 MDEA 实 际升温曲线在 30 ℃左右仍有“下凹”现象，说明“下 凹”现象不是温差热损失引起的，而是由于反应达 到近平衡状态时正向吸收反应放热量减少所致，因 此，MDEA 升温曲线出现的“下凹”现象仍然归因于 自身产热量减少，在计算反应热时不应忽略，从而提 高反应热测量值的准确度． 4 结论 ( 1) 同步热跟踪法测微量气液反应热与文献中 数据对比，实验结果一致; 与非同步热跟踪法相比， 无论是吸收气量、升温区间的计算，还是最终得到的 反应热值，结果都更为准确，实验精度高． ( 2) 随溶液质量分数的增加，MEA 的反应热值 呈现先降低后升高的趋势，在质量分数为 30% 左右 达到最低值，而 MDEA 的反应热值则一直处于下降 趋势． MEA、MDEA 溶液质量分数在 20% ～ 40% 时， 对反应热的影响不显著． ( 3) 胺溶液吸收 CO2反应对应的升温曲线在开 始阶段一直呈直线上升趋势，反应进行一段时间后， 反应产物不断累积，开始抑制正向吸收反应的进行， 使曲线出现“下凹”，将该部分反应热考虑到计算中 可以提高结果的准确度． 参 考 文 献 ［1］ Du M，Zhang L，Feng B． Effect of pH controlling method on ener￾gy consumption of CO2 desorption from rsch-solvent． J Chongqing Univ，2010，33( 8) : 123 ( 杜敏，张力，冯波． pH 值调节剂对 CO2 吸收富液解吸能耗 的影响． 重庆大学学报，2010，33( 8) : 123) ［2］ Qin F，Wang S J，Svendsen H F，et al． Ｒesearch on heat require￾ment of aqua ammonia regeneration for CO2 capture． CIESC J， 2010，61( 5) : 1233 ( 秦锋，王淑娟，Svendsen H F，等． 氨法脱碳系统再生能耗的 研究． 化工学报，2010，61( 5) : 1233) ［3］ Zhang K F，Liu Z L，Wang Y Y，et al． Analysis and calculation of the desorption energy consumption of CO2 capture process by chemical absorption method． Chem Ind Eng Prog，2013，32 ( 12) : 3008 ( 张克舫，刘中良，王远亚，等． 化学吸收法 CO2 捕集解吸能 耗的分析计算． 化工进展，2013，32( 12) : 3008) ［4］ Ahmady A，Hashim M A，Aroua M K． Absorption of carbon diox￾ide in the aqueous mixtures of methyldiethanolamine with three types of imidazolium-based ionic liquids． Fluid Phase Equilib， 2011，309( 1) : 76 ［5］ Zhang F，Ma J W，Zhou Z，et al． Study on the absorption of car￾bon dioxide in high concentrated MDEA and ILs solutions． Chem Eng J，2012，181-182: 222 ［6］ Arcis H，Ｒodier L，Ballerat-Busserolles K，et al． Modeling of ( vapor + liquid) equilibrium and enthalpy of solution of carbon di￾oxide ( CO2 ) in aqueous methyldiethanolamine ( MDEA) solu￾tions． J Chem Thermodyn，2009，41( 6) : 783 ［7］ i L E，Lundberg J，Pedersen M，et al． Measurements of CO2 ab￾sorption and heat consumption in laboratory rig． Energy Procedia， 2014，63: 1569 ［8］ Qi G J，Wang S J，Yu H，et al． Prediction model of absorption heat for CO2 capture using aqueous ammonia． CIESC J，2013，64 ( 9) : 3079 ( 齐国杰，王淑娟，Yu Hai，等． 氨水吸收 CO2 的吸收热预测 模型． 化工学报，2013，64( 9) : 3079) ［9］ Liu B C，Shi C H，Li Q L，et al． Experimental investigation of performances and renewable energy consumption of absorbent DEA in CO2 capture system． Eng J Wuhan Univ，2012，45( 6) : 817 ( 刘炳成，史澄辉，李庆领，等． 电厂 CO2 捕集系统 DEA 化学 吸收剂性能与能耗实验研究． 武汉大学学报( 工学版) ，2012， 45( 6) : 81) ［10］ Chen J，Luo W L，Li H． A review for research on thermodynam￾ics and kinetics of carbon dioxide absorption with organic amines． CIESC J，2014，65( 1) : 12 ( 陈健，罗伟亮，李晗． 有机胺吸收二氧化碳的热力学和动 力学研究进展． 化工学报，2014，65( 1) : 12) ［11］ Gupta M，Svendsen H F． Temperature dependent enthalpy of CO2 absorption for amines and amino acids from theoretical calcu￾lations at infinite dilution． Energy Procedia，2014，63: 1106 ［12］ Van Nierop E A，Hormoz S，House K Z，et al． Effect of absorp￾tion enthalpy on temperature-swing CO2 separation process per￾formance． Energy Procedia，2011，4: 1783 ［13］ Liu J Z，Wang S J，Zhao B，et al． Absorption of carbon dioxide in aqueous ammonia． Energy Procedia，2009，1( 1) : 933 ［14］ Liu J Z，Wang S J，Qi G J，et al． Kinetics and mass transfer of carbon dioxide absorption into aqueous ammonia． Energy Proce￾dia，2011，4: 525 ［15］ Srisang W，Pouryousefi F，Osei P A，et al． Evaluation of the heat duty of catalyst-aided amine-based post combustion CO2 cap￾ture． Chem Eng Sci，2017，170: 48 ［16］ Abdulkadir A，Ｒayer A V，Quang D V，et al． Heat of absorption and specific heat of carbon dioxide in aqueous solutions of mono- · 573 ·

·376 工程科学学报，第41卷，第3期 ethanolamine,3-piperidinemethanol and their blends.Energy CO2 with alkanolamine solutions predicted from reaction equilibri- Procedia,2014,63:2070 um constants.Chem Eng Sci,2009.64(9)2027 [17]Xie Q,Aroonwilas A,Veawab A.Measurement of heat of CO2 25] Arcis H,Ballerat-Busserolles K,Rodier L,et al.Measurement absorption into 2-mino2-methyl--propanol (AMP)/piperazine and modeling of enthalpy of solution of carbon dioxide in aqueous (PZ)blends using differential reaction calorimeter.Energy Pro- solutions of diethanolamine at temperatures of (322.5 and cedia,2013,37:826 372.9)K and pressures up to 3 MPa.J Chem Eng Data,2012, [18]Arcis H,Rodier L,Ballerat-Busserolles K,et al.Enthalpy of so- 57(3):840 lution of CO in aqueous solutions of methyldiethanolamine at T= 6]Wang X N.Design and implementation of constant temperature 372.9K and pressures up to 5 MPa.J Chem Thermodyn,2009, control software based on improved PID.Comput Simul,2015, 41(7):836 32(4):371 [19]Arcis H,Rodier L,Ballerat-Busserolles K,et al.Enthalpy of so- (王晓娜.基于改进PD的恒温控制软件设计与实现.计算 lution of CO,in aqueous solutions of methyldiethanolamine at T= 机仿真，2015,32(4)：371) 322.5K and pressure up to 5 MPa.J Chem Thermodyn,2008, 7]Li J,Gu JJ.Development of tuning methods for PID parameters. 40(6):1022 Inform Electr Pouer,2001 (3):11 0]Arcis H.Rodier L,Coxam J Y.Enthalpy of solution of CO2 in (李剑，谷俊杰.PD参数整定方法进展.电力情报，2001 aqucous solutions of 2-amino2-methyl--propanol.Chem Ther- (3):11) mon,2007,39(6):878 [28]Yang Z Y,Guan Y C,Zhao S W.Temperature controller design 21]Kim I,Svendsen H F.Heat of absorption of carbon dioxide based on self-uning PID.Shandong Ind Technol,2017(3):3 (CO2)in monoethanolamine (MEA)and 2-(Aminoethyl)etha- (杨振元，关艳翠，赵硕伟.基于自整定P⑩的温度控制器 nolamine (AEEA)solutions.Ind Eng Chem Res,2007,46 设计.山东工业技术，2017(3)：3) (17):5803 129]Mathonat C,Majer V,Mather A E,et al.Use of flow calorimetry 22]Guo D F,Gao S W,Luo W L,et al.Effect of sulfolane on CO2 for determining enthalpies of absorption and the solubility of CO absorption and desorption by monethanolamine aqueous solution. in aqueous monoethanolamine solutions.Ind Eng Chem Res, CIESC J,2016,67(12):5244 1998,37(10):4136 (郭东方，部时旺，罗伟亮，等.环丁砜对乙醇胺溶液吸收和 B0]Kim I,Svendsen H F.Comparative study of the heats of absorp- 解吸C02的影响.化工学报，2016,67(12)：5244) tion of post-combustion CO2 absorbents.Int J Greenhouse Gas 23]Svensson H,Hulteberg C,Karlsson H T.Heat of absorption of Control,2011,5(3):390 CO2 in aqueous solutions of N-methyldiethanolamine and pipera- 31]Kim I,Hoff K A.Mejdell T.Heat of absorption of CO2 with zine.Int J Greenhouse Gas Control,2013,17:89 aqueous solutions of MEA:new experimental data.Energy Proce- 24]Kim I,Hoff K A,Hessen E T,et al.Enthalpy of absorption of dia,2014,63:1446

工程科学学报，第 41 卷，第 3 期 ethanolamine，3-piperidinemethanol and their blends． Energy Procedia，2014，63: 2070 ［17］ Xie Q，Aroonwilas A，Veawab A． Measurement of heat of CO2 absorption into 2-Amino-2-methyl-1-propanol ( AMP) / piperazine ( PZ) blends using differential reaction calorimeter． Energy Pro￾cedia，2013，37: 826 ［18］ Arcis H，Ｒodier L，Ballerat-Busserolles K，et al． Enthalpy of so￾lution of CO2 in aqueous solutions of methyldiethanolamine at T = 372. 9 K and pressures up to 5 MPa． J Chem Thermodyn，2009， 41( 7) : 836 ［19］ Arcis H，Ｒodier L，Ballerat-Busserolles K，et al． Enthalpy of so￾lution of CO2 in aqueous solutions of methyldiethanolamine at T = 322. 5 K and pressure up to 5 MPa． J Chem Thermodyn，2008， 40( 6) : 1022 ［20］ Arcis H，Ｒodier L，Coxam J Y． Enthalpy of solution of CO2 in aqueous solutions of 2-amino-2-methyl-1-propanol． J Chem Ther￾modyn，2007，39( 6) : 878 ［21］ Kim I，Svendsen H F． Heat of absorption of carbon dioxide ( CO2 ) in monoethanolamine ( MEA) and 2-( Aminoethyl) etha￾nolamine ( AEEA) solutions． Ind Eng Chem Ｒes，2007，46 ( 17) : 5803 ［22］ Guo D F，Gao S W，Luo W L，et al． Effect of sulfolane on CO2 absorption and desorption by monethanolamine aqueous solution． CIESC J，2016，67( 12) : 5244 ( 郭东方，郜时旺，罗伟亮，等． 环丁砜对乙醇胺溶液吸收和 解吸 CO2 的影响． 化工学报，2016，67( 12) : 5244) ［23］ Svensson H，Hulteberg C，Karlsson H T． Heat of absorption of CO2 in aqueous solutions of N-methyldiethanolamine and pipera￾zine． Int J Greenhouse Gas Control，2013，17: 89 ［24］ Kim I，Hoff K A，Hessen E T，et al． Enthalpy of absorption of CO2 with alkanolamine solutions predicted from reaction equilibri￾um constants． Chem Eng Sci，2009，64( 9) : 2027 ［25］ Arcis H，Ballerat-Busserolles K，Ｒodier L，et al． Measurement and modeling of enthalpy of solution of carbon dioxide in aqueous solutions of diethanolamine at temperatures of ( 322. 5 and 372. 9) K and pressures up to 3 MPa． J Chem Eng Data，2012， 57( 3) : 840 ［26］ Wang X N． Design and implementation of constant temperature control software based on improved PID． Comput Simul，2015， 32( 4) : 371 ( 王晓娜． 基于改进 PID 的恒温控制软件设计与实现． 计算 机仿真，2015，32( 4) : 371) ［27］ Li J，Gu J J． Development of tuning methods for PID parameters． Inform Electr Power，2001( 3) : 11 ( 李剑，谷俊杰． PID 参数整定方法进展． 电力情报，2001 ( 3) : 11) ［28］ Yang Z Y，Guan Y C，Zhao S W． Temperature controller design based on self-tuning PID． Shandong Ind Technol，2017( 3) : 3 ( 杨振元，关艳翠，赵硕伟． 基于自整定 PID 的温度控制器 设计． 山东工业技术，2017( 3) : 3) ［29］ Mathonat C，Majer V，Mather A E，et al． Use of flow calorimetry for determining enthalpies of absorption and the solubility of CO2 in aqueous monoethanolamine solutions． Ind Eng Chem Ｒes， 1998，37( 10) : 4136 ［30］ Kim I，Svendsen H F． Comparative study of the heats of absorp￾tion of post-combustion CO2 absorbents． Int J Greenhouse Gas Control，2011，5( 3) : 390 ［31］ Kim I，Hoff K A，Mejdell T． Heat of absorption of CO2 with aqueous solutions of MEA: new experimental data． Energy Proce￾dia，2014，63: 1446 · 673 ·