2. 西北核技术研究所, 西安 710024
2. Northwest Institute of Nuclear Technology, Xi'an 710024, China
放射性氙同位素(131Xe、133Xe和135Xe)的裂变产额大,较易释放到大气中,且半衰期适中,是全面禁止核实验条约的现场视察(OSI)核素监测中重要的放射性气体核素[1-2]。放射性氙同位素的其他来源包括核电站及医用核设施,对这些设施释放出的放射性氙同位素进行监测能够为设施的安全和防护性能评估提供数据支持。但是,空气中氙本底的体积分数仅为0.087×10-6,因此在放射性氙的监测或制备过程中通常采用物理吸附的方法对氙进行富集。由于空气中CO2等杂质气体的存在,会造成吸附剂对氙的吸附效能下降或在超低温环境下造成吸附柱堵塞,因此需要在氙吸附前进行CO2去除[3]。CO2的去除通常采用膜分离技术和沸石分子筛的物理吸附方法[4]。膜分离技术利用混合气体中不同气体组分在膜内溶解、扩散性质不同,根据渗透率大小达到气体分离的目的,其性能受气体分子的动力学直径、体积、压力、温度以及膜材料影响,其操作较为简便,但对CO2的去除不彻底[4]。沸石分子筛为结晶型且具有多孔结构,其晶格中有许多大小相同的空穴可容纳被吸附的分子,拥有规整且孔径均一的孔道结构,属于强极性吸附材料。由于CO2含有高极性四极矩,能够与沸石中的阳离子反应,从而被沸石分子筛较好地吸附[5-7]。而且,炭基吸附剂如碳分子筛对于非极性分子Xe占有优势,而对极性分子CO2的吸附能力相对较差。另外,沸石分子筛的吸附性能与炭基吸附剂相反,且成本低廉、操作简便。因此,在放射性氙吸附分离过程中采用沸石分子筛去除CO2,将有助于减小氙的损失,是较为优越的方法。本文研究了不同种类的沸石分子筛对CO2的除杂性能以及操作条件对除杂性能的影响,能够为放射性氙吸附分离系统的CO2除杂柱设计提供参考依据,同时也能为其他CO2的去除应用提供参考。
1 测量原理本文主要通过测试沸石分子筛对CO2的动态穿透曲线来表征除杂性能。吸附装置的气态流出物浓度曲线称为穿透曲线,反映了吸附质在吸附过程中的组分浓度变化,曲线呈S型。当体积分数为C0的吸附质以一定的流速通过吸附装置,可在不同的时刻监测气态流出物的体积分数Cx,将Cx/C0—t作图,即可得到穿透曲线,如图 1所示[8]。Wood[9-10]将穿透曲线划分为穿透点(5%)、 中点(50%)和饱和点(95%),对应时间t0.05、 t0.5、 t0.95。 穿透曲线的横坐标可用时间t或者比吸附体积表示。比吸附体积(specific adsorption volume)计算公式如下:
$比吸附体积=\frac{气态流出物体积}{吸附柱体积}=\frac{F\cdot t}{{{V}_{c}}}.$ | (1) |
其中:F为吸附质的体积流速,单位mL/min;t为吸附过程的时间,单位 min;Vc为吸附柱的容积,单位cm3。
根据流量率、吸附剂质量等参数,能够计算出动态吸附系数,用于吸附柱设计和优化。动态吸附系数计算公式如下[11-13]:
$~{{k}_{dB}}=\frac{F\cdot {{t}_{0.05}}}{m},$ | (2) |
${{k}_{d}}=\frac{F\cdot {{t}_{0.5}}}{m},$ | (3) |
$~{{k}_{dS}}=\frac{F\cdot {{t}_{0.95}}}{m}.$ | (4) |
其中:kdB、 kd、kdS分别为穿透动态吸附系数、动态吸附系数和饱和动态吸附系数,表示在穿透时刻t0.05、 t0.5、 t0.95时单位质量吸附剂所吸附的气体体积,单位L/g;F为气体流出吸附柱的皂沫流量计流量,单位mL/min;m为吸附剂的质量,单位g。
由于t0.5对应的C0.5/C0变化最明显,通常用kd表示吸附能力的大小。根据实验测得的kdB、 kd、 kdS,结合气体的采样体积,可以计算出吸附柱内吸附剂的填充质量,进而利用填充密度算出吸附柱的柱内体积,从而给出吸附柱尺寸的设计参数。为了获得更高的气体回收率,通常在气体采集的吸附柱设计中采用kdB,而在气体进一步纯化流程中减小吸附柱体积而采用kd或者kdS设计吸附 柱。
2 实验设计沸石分子筛的CO2去除性能实验平台如图 2所示,CO2体积分数通过N2载气调节,质量流量控制器用于气体流量的控制,调节阀用于控制气路的压力,分别使用皂沫流量计、压力传感器、露点仪测定气体的流量、压力和湿度,使用HP6890气相色谱仪测量气体中CO2体积分数,采用恒温槽控制吸附柱的温度并使其保持恒定。吸附柱分别填充不同类型沸石分子筛作为吸附剂,其中N2为非极性分子,且沸点极低(77 K),因此很难被沸石分子筛吸附,可用作载气。
HP6890的CO2分析条件为:炉温80℃,TCD检测器温度205℃;采用PORAPAK Q色谱柱(填料粒径0.18~0.25 mm、 柱长2 m、 外径 3.18 mm);柱流量10 mL/min,参比流量26 mL/min,尾吹流量2 mL/min。在22×10-6~2 178×10-6的体积分数范围内,CO2体积分数测量的标准曲线如图 3所示。
吸附剂选用上海恒业化工生产的5A、 13XAPG和13XHP分子筛 (粒径1.4~2.36 mm),3种吸附剂的结构表征数据见表 1。
材料 | BET比 表面积/(m2·g-1) | 总孔容/(cm3·g-1) | 总孔径/nm | 微孔 体积/(cc·g-1) | 特征 吸附能/(kJ·mol-1) | 微孔 孔径/nm |
5A | 622 | 0.34 | 2.17 | 0.24 | 34.3 | 0.47 |
13XAPG | 690 | 0.37 | 2.14 | 0.26 | 34.6 | 0.47 |
13XHP | 516 | 0.37 | 2.85 | 0.2 | 27.5 | 0.67 |
3 结果与讨论
1) CO2吸附剂的选取
测试了5A、13XAPG和13XHP分子筛对CO2的动态吸附性能。吸附柱尺寸为Φ1cm×40 cm,温度为25℃。吸附柱装填质量及测试结果见表 2,穿透曲线如图 4所示。
材料 | 质量/g | C0/10-6 | F/(mL·min-1) | kdB/(L·g-1) | kd/(L·g-1) |
13XAPG | 9.5 | 116.0 | 564.3 | 20.6 | 25.1 |
13XHP | 8.7 | 115.4 | 564.3 | 25.7 | 29.6 |
5A | 10.1 | 118.4 | 560.2 | 7.7 | 11.0 |
如表 2和图 4所示,5A、13XAPG和13XHP分子筛对CO2的动态穿透曲线呈S型,分子筛对CO2吸附能力按由弱到强排序为5A—13XAPG—13XHP。13XAPG和13XHP分子筛的穿透动态吸附系数和动态吸附系数均比5A分子筛大,其中13XHP分子筛的穿透动态吸附系数和动态吸附系数最大,分别达到25.7和29.6 L/g,由此可见,相比其他2种吸附剂,13XHP去除气体中CO2具有优势。这是因为13XHP的总孔容和总孔径比5A和13XAPG要大,且微孔孔径(0.67 nm)大于2倍的CO2的空气动力学分子直径(0.33 nm),因此更有利于CO2的吸附。周崇阳[4]测试了25℃时,13XHP对氙的动态吸附性能,结果表明穿透动态吸附系数仅为0.015 L/g,饱和动态吸附系数仅为0.041 L/g,同时,13XHP分子筛的特征吸附能为 27.5 kJ/mol,说明13XHP对非极性分子氙的吸附能力很弱,对极性分子吸附能力强。因此,在放射性氙测量过程中选择13XHP分子筛去除CO2是合适的。
2) CO2的体积分数影响
测试了13XHP对不同体积分数的CO2的动态吸附性能,测试结果如图 5所示。图 6为不同CO2体积分数的kdB和kd的变化曲线。
对图 6中的ln C0~ln kdB和ln C0~ln kd分别作曲线拟合,得到如下:
$\begin{align} & ~ln~{{k}_{dB}}=-0.06{{\left( ln~{{C}_{0}} \right)}^{2}}+0.486ln~{{C}_{0}} \\ & +2.309\left( {{R}^{2}}=1 \right), \\ \end{align}$ | (5) |
$\begin{align} & ln~{{k}_{d}}=-0.099{{\left( ln~{{C}_{0}} \right)}^{2}}+0.944ln~{{C}_{0}}+ \\ & 1.148\left( {{R}^{2}}=0.999 \right). \\ \end{align}$ | (6) |
如图 5和6所示,随着CO2体积分数的增加,13XHP对CO2的动态吸附能力呈下降趋势,穿透动态吸附系数和动态系数也呈下降的趋势。为了减小除杂柱的体积,提升除杂效率,有必要在除杂柱前端添加膜分离装置去除气体杂质组分,然后再使用分子筛进行深度除杂。式(5)和(6)可用于不同CO2体积分数情况下除杂柱尺寸设计。
3) 气体流量影响
测试了不同气体流量下13XHP对CO2的动态吸附性能影响,测试结果如图 7和表 3所示。
F/(mL·min-1) | kdB /(L·g-1) | kd /(L·g-1) |
564 | 25.7 | 29.6 |
909 | 26.3 | 33.1 |
1322 | 25.9 | 34.3 |
1637 | 25.5 | 34.0 |
2015 | 25.6 | 34.9 |
如图 7和表 3所示,随流量的增加,13XHP的穿透动态吸附系数在平均值25.8 L/g附近波动;动态吸附系数呈现轻微升高的趋势,但变化幅度不大。因此,在气体流量564~2015 mL/min范围内,气体流量的变化对13XHP的动态吸附性能影响不大。
4) 温度影响
测试了不同温度下,13XHP分子筛对CO2的动态吸附性能,结果如图 8和9所示。
对图 9中的ln T~ln kdB和ln C0~ln kd分别作曲线拟合,得到如下:
$\begin{align} & ln~{{k}_{dB}}=-0.408{{(ln~T)}^{2}}+1.256ln~T+ \\ & 3.433\left( {{R}^{2}}=1 \right), \\ \end{align}$ | (7) |
$\begin{align} & ln~{{k}_{d}}=-0.403{{\left( ln~T \right)}^{2}}+1.267ln~T+ \\ & 3.512\left( {{R}^{2}}=0.999 \right). \\ \end{align}$ | (8) |
如图 8和9所示,随着温度的升高,13XHP分子筛对CO2的动态吸附性能下降,5 ℃时13XHP的动态吸附系数远大于其他温度。因为随着温度的升高,气体分子的热力学运动加快,这不利于CO2的吸附。但是,采用5 ℃甚至更低的温度进行除杂,需附加低温致冷装置,将会大大增大系统体积,因此在除杂柱设计时考虑选择25℃的动态吸附系数作为设计参数。
5) 柱前压力影响
测试了不同柱前压力下,13XHP分子筛对CO2的动态吸附性能,结果如图 10和表 4所示。
如图 10和表 4所示,在柱前压力103~353kPa 范围内,随着柱前压力的增加穿透曲线无明显变化,穿透动态系数和动态吸附系数变化不大,因此在该压力变化范围内13XHP对CO2的动态吸附性能基本不受影响。
4 测量结果不确定度分析沸石分子筛对CO2动态吸附系数测量结果的不确定度分项包括:CO2标准气体体积分数的不确定度小于2%;流量控制器的流量读数的不确定度小于2%;压力传感器读数的不确定度小于2%;气相色谱测量的峰面积读数的不确定度小于3%;温度传感器读数的不确定度小于1%。以上各分项相互独立,合成不确定度小于10% (k=2)。
5 结 论本文测试了不同类型分子筛在常温下对CO2的动态吸附性能,以及操作条件对动态吸附性能的影响。结果表明:13XHP沸石分子筛对CO2的去除性能较好,温度和CO2的体积分数的影响最为显著,气体流量和柱前压力分别在564~2015 mL/min和103~353 kPa范围内对吸附剂性能影响不大。该结果能够为放射性氙吸附分离系统的CO2除杂柱设计提供参考依据,同时也能为CO2去除的其他应用提供参考。
[1] | 龙斌, 吴睿, 成智威, 等. 放射性氙测量中的天然源γ谱仪校准技术[J]. 原子能科学技术 , 2014, 48 (1) : 137–141. LONG Bin, WU Rui, CHENG Zhiwei, et al. Natural source-based calibration technique for measurement of radio-xenon with γ spectrometry[J]. Atomic Energy Science and Technology , 2014, 48 (1) : 137–141. (in Chinese) |
[2] | ZHOU Chongyang, ZHOU Guoqing, FENG Shujuan, et al. The determination of trace amount of natural 133Xe in gas under high radon activity concentration levels[J]. J Environ Radioactivity , 2013, 122 : 9–15. DOI:10.1016/j.jenvrad.2013.02.006 |
[3] | 冯淑娟, 周崇阳, 周国庆, 等. 氙在活性炭和碳分子筛上的动态吸附性能[J]. 核化学与放射化学 , 2010, 32 (5) : 274–279. FENG Shujuan, ZHOU Chongyang, ZHOU Guoqing, et al. Dynamic adsorption property of xenon on activated carbon and carbon molecular sieves[J]. Journal of Nuclear and Radiochemistry , 2010, 32 (5) : 274–279. (in Chinese) |
[4] | 周崇阳.大气样品中氙的富集分离技术研究[D].北京:清华大学,2011. ZHOU Chongyang.Study on the Technology of Xenon Enrichment and Separation from Atmosphere and its Application in Nuclear Inspection[D].Beijing:Tsinghua University,2011.(in Chinese) |
[5] | 王璐.沸石13XAPG吸附分离CO2-N2混合气过程研究及其应用[D].上海:华东理工大学,2013. WANG Lu.Investigation and Application of CO2/N2 Separation by Adsorption Process using Zeolite 13XAPG[D].Shanghai:East China University of Science and Technology,2013.(in Chinese) http://cdmd.cnki.com.cn/Article/CDMD-10251-1013211535.htm |
[6] | LIU Zhen, Grande C A, LI Ping, et al. Adsorption and desorption of carbon dioxide and nitrogen on zeolite 5A[J]. Sep Sci Technol , 2011, 46 (3) : 434–451. DOI:10.1080/01496395.2010.513360 |
[7] | Cavenati S, Grande C A, Rodrigues A E. Adsorption equilibrium of methane,carbon dioxide,and nitrogen on zeolite 13X at high pressure[J]. J Chem Eng Data , 2004, 49 (4) : 1095–1101. DOI:10.1021/je0498917 |
[8] | ZHOU Chongyang, FENG Shujuan, ZHOU Guoqing, et al. The behavior of xenon dynamic adsorption on granular activated carbon packed bed adsorber[J]. J Radional Nucl Chem , 2011, 287 : 609–616. DOI:10.1007/s10967-010-0878-9 |
[9] | Wood G O. A review of the effects of covapors on adsorption rate coefficients of organic vapors adsorbed onto activated carbon from flowing gases[J]. Carbon , 2002, 40 : 685–694. DOI:10.1016/S0008-6223(01)00185-3 |
[10] | Wood G O. Quantification and application of skew of breakthrough curves for gases and vapors eluting from activated carbon beds[J]. Carbon , 2002, 40 : 1883–1890. DOI:10.1016/S0008-6223(02)00031-3 |
[11] | Nakayama Y, Nagao H, Mochida I, et al. Adsorption of radon on active-carbon[J]. Carbon , 1994, 32 : 1544–1547. DOI:10.1016/0008-6223(94)90154-6 |
[12] | Siegworth D P,Neulander C K,Pao R T,et al.Measurement of dynamic adsorption coefficients for noble gases on activated carbon[C]//Proceedings of the 12th AEC Air Cleaning Conference.San Jose,CA,USA:National Technical Information Service,1972:28-47. |
[13] | Adams R E, Browning W E, Ackley R D, et al. Containment of radioactive fission gases by dynamic adsorption[J]. San Jose,CA,USA:National Technical Information Service , 1959, 51 : 1467–1470. |