王曉蕾 張學(xué)軍 陸軍亮 陳 瑤

(浙江大學(xué)制冷與低溫研究所 杭州 310027)

空分用立式徑向流雙層吸附床吸附性能研究

王曉蕾 張學(xué)軍 陸軍亮 陳 瑤

(浙江大學(xué)制冷與低溫研究所 杭州 310027)

闡述了徑向流吸附器在空分領(lǐng)域的應(yīng)用背景，分析了雙層吸附床的工作原理，建立了徑向流雙層吸附床數(shù)值計(jì)算模型，總結(jié)了吸附床層厚度布置對(duì)吸附性能及穿透時(shí)間的影響規(guī)律。計(jì)算結(jié)果表明，吸附反應(yīng)必須以監(jiān)控出口處空氣中二氧化碳的含量作為反應(yīng)結(jié)束的標(biāo)志;在吸附床總厚度和空氣流量不變的情況下，兩種吸附劑的厚度比為0.3時(shí)，雙層床吸附劑的利用率最高。

徑向流吸附器 雙層床 吸附平衡

1 引言

空氣在進(jìn)入精餾塔進(jìn)行分離之前，必須除去其中的水蒸氣，二氧化碳和乙炔等碳?xì)浠衔铮行Х乐乖诤罄m(xù)過(guò)程中冷凝結(jié)冰堵塞管路;乙炔等碳?xì)浠衔镌诳辗衷O(shè)備中的積聚甚至?xí)l(fā)爆炸，給工業(yè)生產(chǎn)帶來(lái)巨大的經(jīng)濟(jì)損失。目前大、中型空分設(shè)備中普遍采用分子篩吸附器來(lái)除去空氣中的雜質(zhì)氣體。傳統(tǒng)的臥式水平床吸附器占地面積大，床層壓降較大，且吸附床容易發(fā)生穿透，在大型空分設(shè)備中逐漸被徑向流吸附器所替代。法國(guó)液化空氣公司于1981年首次研發(fā)了立式徑向流吸附器并獲得實(shí)際應(yīng)用［1］。在相同吸附劑裝載負(fù)荷的前提下，與水平床相比，徑向吸附床與壓縮空氣的接觸面積大大增加而床層厚度卻成倍減小，從而大幅度降低了床層壓降，節(jié)約能耗，因此具有廣闊的發(fā)展前景。

分子篩吸附器內(nèi)發(fā)生的是變溫吸附反應(yīng)(Temperature Swing Adsorption，TSA)。在吸附階段，低溫物料不斷流過(guò)床層，雜質(zhì)氣體被吸附留在床層中，其余組分從床中流出，吸附過(guò)程持續(xù)到吸附床飽和為止。解吸過(guò)程中通過(guò)高溫再生氣，高溫下雜質(zhì)氣體的吸附量降低而從吸附床中解析出來(lái)，吸附劑再生，從而實(shí)現(xiàn)吸附器的循環(huán)利用。徑向流吸附器雖然具有諸多優(yōu)點(diǎn)，但是由于其體積龐大，結(jié)構(gòu)復(fù)雜，采用相應(yīng)的實(shí)驗(yàn)手段展開(kāi)研究具有很高的難度和局限性，數(shù)值計(jì)算融和了理論分析和科學(xué)試驗(yàn)的特點(diǎn)，關(guān)注于物質(zhì)性質(zhì)與反應(yīng)機(jī)理，具有很高的準(zhǔn)確性，被廣泛應(yīng)用于科學(xué)研究，工程與生產(chǎn)領(lǐng)域，尤其適用于大型吸附器的吸附性能研究［2-3］。本文在對(duì)該吸附器內(nèi)二氧化碳及水蒸氣的傳熱傳質(zhì)特性進(jìn)行深入分析的基礎(chǔ)上，建立其吸附過(guò)程的動(dòng)態(tài)模型，并進(jìn)行了理論模擬與計(jì)算，分析了吸附劑床層布置對(duì)吸附性能的影響，旨在為立式徑向流吸附器設(shè)計(jì)提供理論依據(jù)。

2 數(shù)學(xué)模型建立與求解

徑向流吸附器內(nèi)進(jìn)行的變溫吸附過(guò)程分為4個(gè)步驟:升壓，吸附，降壓，解吸。吸附床內(nèi)填充活性氧化鋁和分子篩雙層吸附劑，氧化鋁吸附水蒸氣，分子篩吸附殘余的水蒸氣和二氧化碳。吸附和解吸過(guò)程的氣流方向采用逆向流動(dòng)方式。圖1所示為立式徑向流吸附器的結(jié)構(gòu)示意圖。吸附時(shí)空氣從吸附器底部進(jìn)氣口進(jìn)入，通過(guò)分配器均勻流動(dòng)進(jìn)入最外層環(huán)形空間，由外向內(nèi)徑向依次穿過(guò)氧化鋁和分子篩雙層吸附床，干燥純凈的空氣從集流管上部流出;再生時(shí)高溫污氮?dú)獾牧鲃?dòng)方向相反。

圖1 立式徑向流吸附器結(jié)構(gòu)示意圖1.空氣進(jìn)氣管;2.空氣出口及再生進(jìn)氣口;3.填料口;4.孔板筒;5.分子篩;6.導(dǎo)氣錐;7.活性氧化鋁;8.再生氣出口管。Fig.1 Schematic diagram of vertical radial flow adsorber

式中:ci為i組分在空氣中摩爾濃度，mol/m3，i=1為水蒸氣，i=2為二氧化碳;u為空氣流速，m/s;ρp為吸附劑密度，kg/m3;ε為孔隙率;qi為i組分的吸附量，mol/kg;τ為時(shí)間，s;r為床層徑向位置，m。

(2)吸附質(zhì)氣體能量方程

固定床分子篩吸附器中的氣體吸附是非穩(wěn)態(tài)傳質(zhì)過(guò)程，此過(guò)程中床層內(nèi)氣體濃度分布隨時(shí)間和沿床層位置不斷變化，氣體的流動(dòng)速度、溫度、濃度、壓力等參數(shù)互相耦合，對(duì)吸附分離過(guò)程也產(chǎn)生復(fù)雜的影響。對(duì)于徑向流吸附器，一般認(rèn)為其床層壓降沿軸向分布是均勻的［4］，此時(shí)吸附器內(nèi)氣體流動(dòng)達(dá)到理想的均布狀態(tài)，可以簡(jiǎn)化為一維模型，并進(jìn)行假設(shè):(1)空氣可以視為理想氣體;(2)吸附劑顆?？醋鬟B續(xù)均勻介質(zhì)，物性參數(shù)不發(fā)生變化;(3)忽略吸附床沿軸向上的溫度梯度、濃度梯度;(4)忽略床層壓降。R.T.Yang論證了線(xiàn)性驅(qū)動(dòng)力模型(LDF)對(duì)雙組份吸附模型的適用性［5］。建立完整的雙組份吸附模型包括物料平衡方程，傳質(zhì)方程，吸附等溫線(xiàn)方程及相應(yīng)的邊界條件和初始條件。

(1)質(zhì)量守恒方程:

式中:cpg為空氣比定壓熱容，J/(kg·K);cpp為床層比定壓熱容，J/(kg·K);Tg為空氣溫度，K;Tp為吸附劑溫度，K;ρg為空氣密度，kg/m3;ap為吸附劑顆粒表面積與體積之比，m-1;hp為氣體與吸附顆粒表面?zhèn)鳠嵯禂?shù)，W/(m2·K)。

(3)吸附劑能量方程

吸附反應(yīng)是物理吸附，吸附時(shí)放熱導(dǎo)致床層溫度升高。利用能量守恒關(guān)系，得出床層溫度變化方程

式中:ΔHi為 i組分的吸附熱，J/mol。

(4)傳質(zhì)方程

由于濃度差的存在，氣體中各組份將向固體吸附顆粒表面擴(kuò)散。對(duì)于多組份氣體，一般采用LDF線(xiàn)性驅(qū)動(dòng)力方程［6］

式中:kci為i組分傳質(zhì)系數(shù)，可以結(jié)合Perry在化工手冊(cè)中總結(jié)的系列公式［7］進(jìn)行求解;qi*為 i組分飽和吸附量，mol/kg。

(5)吸附平衡方程

Yao［8］對(duì)擴(kuò)展 Langmuir方程的計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行比較，結(jié)果顯示擴(kuò)展Langmuir方程與實(shí)驗(yàn)結(jié)果有更好地吻合性，且計(jì)算較簡(jiǎn)單，適用于描述雙層床氣固吸附平衡

由Clausius-Clapeyron方程可以推出能量控制方程中的等量吸附熱ΔH的解析形式為［9］:

(6)邊界條件和初始條件

初始條件對(duì)吸附循環(huán)達(dá)到穩(wěn)定狀態(tài)沒(méi)有影響，故選擇干凈吸附床作為初始條件。c(r，τ)r=0= cinlet，T(r，τ)r=0=Tinlet，即吸附器入口處氣流溫度濃度保持不變，即吸附器出口處氣體溫度和濃度變化率為零。c(r，τ)τ=0=0，q(r，τ)τ=0=0，吸附開(kāi)始之前，吸附氣體尚未進(jìn)入吸附床中，因而吸附器內(nèi)固相與氣相中吸附質(zhì)濃度都為零。

式(1)—(6)為一組偏微分方程組，是時(shí)間和空間的函數(shù)，必須進(jìn)行數(shù)值求解。采用Talor展開(kāi)法進(jìn)行時(shí)間和空間上的離散求解。為求解該6元非線(xiàn)性方程組，采用牛頓法構(gòu)造線(xiàn)性迭代矩陣，從而把一個(gè)非線(xiàn)性問(wèn)題轉(zhuǎn)化為求解若干次線(xiàn)性方程組的迭代過(guò)程。圖2所示為徑向流吸附器的計(jì)算單元和網(wǎng)格劃分。表1所示為計(jì)算過(guò)程的物性參數(shù)、結(jié)構(gòu)參數(shù)和運(yùn)行參數(shù)等。

圖2 計(jì)算單元和離散網(wǎng)格劃分Fig.2 Computing unit and discrete mesh

表1 計(jì)算參數(shù)Table 1 Parameters of calculation

3 討論與分析

在徑向流吸附器運(yùn)行過(guò)程中，低溫壓縮空氣徑向流過(guò)吸附床，流動(dòng)速率隨床層半徑的減小而增大，如圖3所示。氣體流動(dòng)速度越快，單位長(zhǎng)度吸附床達(dá)到吸附飽和狀態(tài)的時(shí)間越短，氣體濃度波面前進(jìn)越快，參數(shù)方程中的傳熱傳質(zhì)系數(shù)也會(huì)隨之發(fā)生相應(yīng)的變化，從而對(duì)吸附器的動(dòng)態(tài)吸附效果和穿透時(shí)間產(chǎn)生影響［10］。

圖3 吸附床內(nèi)部氣體流動(dòng)速度分布Fig.3 Velocity distribution of feed gas through adsorption bed

吸附床內(nèi)部溫度分布如圖4所示。圖中3條曲線(xiàn)分別代表吸附床不同半徑處在整個(gè)吸附過(guò)程中的溫度變化。r*為所在位置床層半徑比氧化鋁床層厚度的無(wú)量綱數(shù)，r*=0.8表示位于氧化鋁床層內(nèi)，r*=1表示位于雙層床的交界處，r*=1.2表示位于分子篩床層內(nèi)。吸附床入口處溫度迅速上升，形成陡峭的溫度波面并逐漸向床層內(nèi)部推移。計(jì)算結(jié)果顯示床層不同位置處的溫度峰值并不相同，隨著氣體流動(dòng)峰值降低，這是由于水蒸氣被吸附而放出大量的吸附熱，且高于二氧化碳的吸附熱，故氧化鋁床層內(nèi)溫度上升幅度高于分子篩床層。吸附床吸附飽和時(shí)達(dá)到溫度峰值，之后無(wú)吸附熱產(chǎn)生，只存在與空氣的換熱，故溫度下降。

圖4 吸附床不同位置溫度變化曲線(xiàn)Fig.4 Temperature distribution along bed

圖5為氧化鋁床層與分子篩床層的穿透曲線(xiàn)。縱坐標(biāo)分別為吸附床出口處水蒸氣及二氧化碳體積濃度的無(wú)量綱數(shù)。對(duì)比發(fā)現(xiàn)，吸附床內(nèi)二氧化碳首先發(fā)生了穿透，水蒸氣的穿透時(shí)間則遠(yuǎn)遠(yuǎn)超過(guò)二氧化碳，故而在實(shí)際生產(chǎn)中以監(jiān)控出口處二氧化碳的濃度作為判斷吸附循環(huán)時(shí)間的標(biāo)準(zhǔn)是合理的。吸附反應(yīng)結(jié)束后，氧化鋁床層遠(yuǎn)沒(méi)有達(dá)到飽和狀態(tài)，此時(shí)氧化鋁床層沒(méi)有得到充分的利用，造成系統(tǒng)與設(shè)備的浪費(fèi)。故在雙層吸附床的設(shè)計(jì)中，兩種吸附劑的厚度比理論上存在一個(gè)理想的最優(yōu)值，此時(shí)水蒸氣和二氧化碳分別在氧化鋁與分子篩床層同時(shí)發(fā)生穿透，吸附劑達(dá)到最大利用率。定義氧化鋁與分子篩床層的厚度比Ψ，整個(gè)吸附過(guò)程相應(yīng)的的原料空氣凈化量Qc，在原料空氣流量和吸附床總厚度相同的情況下，Ψ與Qc的關(guān)系如圖6所示。當(dāng)Ψ減小，即氧化鋁床層厚度減小，分子篩床層厚度增加時(shí)，空氣凈化量Qc先增大后減小。Ψ＜0.3時(shí)，水蒸氣在氧化鋁床層上吸附飽和后，又進(jìn)入到分子篩床層被吸附，只有部分分子篩吸附二氧化碳飽和，降低了吸附器對(duì)二氧化碳雜質(zhì)的處理量，在實(shí)際生產(chǎn)中會(huì)影響分子篩的使用壽命;Ψ＞0.3時(shí)如圖5所示，二氧化碳先發(fā)生穿透，則吸附過(guò)程結(jié)束，此時(shí)氧化鋁床層尚未飽和，吸附劑沒(méi)有得到充分利用;在Ψ=0.30，即氧化鋁床床層厚度為20 mm，分子篩床層厚度為66.5 mm時(shí)，空氣凈化量達(dá)到最大值，吸附床的利用率最高。

圖5 吸附床穿透曲線(xiàn)Fig.5 Adsorption breakthrough curve

圖6 吸附劑厚度比對(duì)吸附量影響Fig.6 Influence of dual-bed adsorbent thickness ratio

4 結(jié)論

建立了空分用立式徑向流雙層床分子篩吸附器的非等溫非平衡吸附模型，得到了吸附床內(nèi)部的溫度分布，揭示了雙層吸附床的變溫吸附規(guī)律，得到如下結(jié)論:

(1)在徑向流吸附器中，氣體穿越吸附床時(shí)流動(dòng)速率隨床層半徑的減小而增大，使得單位長(zhǎng)度吸附床達(dá)到吸附飽和狀態(tài)的時(shí)間縮短。

(2)吸附床入口處溫度迅速上升，形成陡峭的溫度波面，在吸附飽和時(shí)達(dá)到溫度峰值，之后與空氣換熱溫度降低;水蒸氣的吸附熱高于二氧化碳，故氧化鋁床層內(nèi)溫度高于分子篩床層。

( 3) 二氧化碳只被分子篩吸附，而水蒸氣既能被氧化鋁吸附，又能被分子篩吸附，故二氧化碳先于水蒸氣發(fā)生穿透吸附床，在實(shí)際生產(chǎn)中必須選擇監(jiān)控出口處二氧化碳的濃度作為判斷吸附床吸附飽和的標(biāo)準(zhǔn)，且吸附循環(huán)時(shí)間不能超過(guò)二氧化碳穿透所需的時(shí)間。

( 4) 當(dāng)雙層吸附劑厚度分別為20 mm 和66. 5mm 時(shí)，既保證雙層吸附床同時(shí)達(dá)到最高的利用效率，又能夠有效延長(zhǎng)設(shè)備的使用壽命，充分發(fā)揮雙層吸附床的優(yōu)勢(shì)。

1 楊學(xué)申. 立式雙層床徑向流吸附器［J］. 深冷技術(shù)，1999( 2) : 16-19.

2 H Ahn，Lee C H. Adsorption dynamics of water in layered bed for airdryingTSA process［J］. Separations，2003，49( 6) : 1601-1609.

3 Schoofs G R. Repressurization of adsorption purifiers for cryogenic airseparation［J］. Chem. Eng. Data，1993，48( 4) : 753-760.

4 Kim J H，Lee C H，Kim W S，et al. Adsorption equilibria of watervapor on alumina，zeolite 13X，and a zeolite X/activated carbon composite［J］. Chem. Eng. Data，2003，48( 1) : 137-141.

5 Kapoor A，Yang R T. Kinetic separation of methane-carbon dioxidemixture by adsorption on molecular sieve carbon［J］. Chem. Eng.Sci. ，1989，44( 8) : 1723.

6 Glueckauf E，Coates J J. Theory of chromatography. Part IV. The influencesof incomplete equilibrium on the front boundary of chromatogramsand on the electiveness of separation［J］. Journal of ChemicalSociety，1947，1315-1321.

7 Perry R H，Green D W，Maloney J O. Perry’s chemical engineershandbook( 7th Ed.) ［M］. America: McGraw-Hill，1999.

8 Yao C C. Extended and improved Langmuir equation for correlatingadsorption equilibrium data. Separation and Purification Technology［J］. 2000( 19) : 237-242.

9 Ko D，Moon I，Choi D K. Analysis of the contact time in a cyclicthermal swing adsorption process［J］. Ind. Eng. Chem. Res. ，2002，

41( 6) : 1603-1615.

10 Malek A，F(xiàn)arooq S，Rathor M N， et al. Effect of velocity variation dueto adsorption-desorption on equilibrium data from breakthrough experiments［J］，Chem. Eng. Sci. ，1995，50( 4) : 737-740.

Study of adsorption performance of vertical radial flow dual-bed adsorber used in air sparation unit

Wang Xiaolei Zhang Xuejun Lu Junliang Chen Yao

(Institute of Refrigeration and Cryogenics，Zhejiang University，Hangzhou 310027，China)

The application background of vertical radial flow adsorber used in cryogenic air separation unit was discussed and working principle of dual-bed adsorption was analyzed.A numerical model to describe the heat and mass transfer process in the packed bed was established，revealing the influence law of dual-bed thickness arrangement on adsorption properties and penetration time.The result shows that carbon dioxide at the outlet must be monitored as the end of the reaction.Besides，with the total bed thickness as well as flow rate of the feed air remains constant，the dual-bed adsorber reaches its best performance when the thickness ratio of the two adsorbents is 0.3.

radial flow adsorber;dual-bed;adsorption equilibrium

TB66

A

1000-6516(2013)01-0019-04

2012-12-10;

2013-02-01

重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973項(xiàng)目)(No.2011CB706501)和國(guó)家自然科學(xué)基金(No.51176164)資助。

王曉蕾，女，25歲，碩士研究生。