H2S是一种有毒有害气体,同时会对管路造成腐蚀,故天然气在使用之前必须进行脱硫处理。目前,最常用的脱硫方法为湿法脱硫[1-3],如醇胺吸收法、络合铁法[4]等。超重力技术作为一项新兴的高科技技术,与传统传质技术存在很多不同,它利用旋转产生的强大离心力和剪切力,使液体在多孔填料中高度分散,与气体充分接触,极大地强化了传质和反应[5-7],具有广泛的应用前景。
超重力场通过高速旋转填料模拟实现,其结构如图 1[8-9]所示。它由固定的圆柱形壳体、转轴、转子、圆环状的填料和液体分布器等组成。旋转床产生的强大离心加速度可高达几百甚至上千g(视旋转半径和转速而定)。
由于超重力技术具有上述特点,将其应用于天然气脱硫工艺可取得优于传统脱硫工艺的效果,因而有广泛的工业应用前景[10-11]。H2S脱除率的计算对超重力脱硫技术的研究较为重要,目前对超重力脱硫模型的研究大多局限于理论层面,很少涉及经验模型,故建立一个传质系数经验模型不仅简单可行,而且可以减少很多后续的工作量。
采用甲烷和H2S的混合气模拟天然气,在中试装置超重力机中进行了脱硫实验。利用数学软件Matlab对实验数据进行拟合处理和分析,获得了传质系数经验模型,并对其进行了可靠性验证。
目前,关于超重力模型的理论研究较多,但直接通过操作条件获得传质系数的经验模型却鲜少有人研究。本研究针对络合铁脱硫过程建立传质系数的经验模型,从而获得更加简单实用的传质模型。利用中试装置的实验数据进行分析处理,该超重力装置的旋转填料内径为40 mm,外径150 mm,填料高度40 mm。对于本研究所采用的络合铁法脱除天然气中的H2S,其传质是受气膜控制的。定义$\theta = \frac{{{y_{{\rm{in}}}}-{y_{{\rm{out}}}}}}{{{y_{{\rm{in}}}}}}$,其中θ为脱除率,yin为进气H2S体积分数,%,yout为出口H2S体积分数,%。根据文献[10-12]可以得到以气相浓度差表示的总传质系数Kya与脱除率之间的关系,见式(1)。
式中:Ky为气相总传质系数,m/h;a为填料比表面积,m2/m3;G为进气流量,m3/h;H为填料的轴向长度,m;r1、r2为填料内外半径,m。
此传质系数模型反映的是气量与超重力装置旋转填料的内外径、填料高度及脱除率之间的关系,只反映了气量这一个操作条件,没有表现出旋转过程中气液间的相互作用。并且此模型需要测得脱除率θ,使得模型在计算时更为复杂。因此,需建立一个以络合铁脱硫过程为基础、关于操作条件的传质系数模型,此模型仅需通过超重力机的操作条件就可得到传质系数,是更简单实用的经验模型。
原料天然气中H2S质量浓度为44.65 g/m3,体积分数为2.94%,将实验数据列于表 1,利用回归方法拟合出关于操作条件G、L、ω之间的经验关系式,见式(2)。
式中:α0、α1、α2、α3为系数;ω为转速,r/min;L为液相流量,m3/h。
根据回归方法对式(1)取对数,然后用MATLAB方法编程对实验数据进行拟合从而可以得到拟合出来的初步的经验模型见式(3)。
将式(3)得到的Kya代入到式(1)中,便可以得到计算的θ值。
偏差计算式见式(4)、式(5)。
式中:e为平均偏差;e′为偏差计算式;θcal,i为计算得到的脱除率;θexp,i为实验获得的脱除率。
对经验模型进行计算可以得到传质系数计算值,将计算值与实验得到的传质系数值进行对比,并计算它们的偏差,如图 2所示。
计算得到平均偏差e=0.382 7。从图 2中也可观察到偏差很大,造成这种严重偏差的原因是拟合的方法从宏观的角度出发,如旋转填料床的结构、物系性质及操作参数等因数对实验数据的影响,造成了结果的偏差,最终导致拟合结果偏差很大。因此,Kya=0.423 4G0.365 8L0.619 4ω0.001 554不能直接作为经验传质模型,需对其进行修正。
由于该模型为气膜控制过程,故气量对传质系数的影响很大,且在超重力场中,转速也存在一定程度的影响[13]。因此,α3的值应进行修正。根据朱慧铭教授[8]用水吸收空气中的氨气这一系统进行的超重力传质研究得出的KGa∝ω0.32L0.48,以及Munjal[14-15]用NaOH吸收空气中CO2系统得到的正比关系式KL∝ω0.34可知,α3应远大于0.001 554,这里假设α3与朱慧铭教授的实验数据相等,为0.32。
由此可以得到新的传质系数模型,见式(6)。
偏差如图 3所示。
平均偏差e=0.002 075,偏差已经很小,但由于本研究所用的体系是络合铁氧化还原法脱除天然气中的H2S超重力机高速旋转,液体在高分散、强混合及界面快速更新的环境下与气体充分接触,对于这一传质过程,根据实验数据可知,其气体处于湍流状态,故可以认为其传质可用表面更新理论[16-17]来描述,见式(7)。
式中:kG′为以气相浓度差表示的气相传质分系数,m/s;DG为气相扩散系数,m2/s;SG为气相的表面更新率,1/s。
很明显,气体的表面更新率和气体的湍动深度有关。在本实验中气体所受到的阻力与其速度相关。因此,可设气体的表面更新率与气速的平方成正比,见式(8)和式(9)。
式中:uG为气体的径向速度,m/s;ks为比例常数,m2/s;N为金属丝网填料的层数。
可以看出,其值和气体在旋转填料床中的运动情况有关,并随填料半径的变化而变化。因此,KGa∝r4a(r4为填料半径,m),再加上传质是气膜控制,可以得到KG=kGa∝r4a,由于指数α4在0~1之间,首先假设α4=0.5,可以得到修正后的传质系数方程,见式(10)。
对传质系数计算值与实验值偏差作图,如图 4所示。
从图 4中可以看出,e′偏差小,且计算出其平均偏差为e=0.025 812,从图 4看,拟合程度不是很好,α4还需要再进行修正。从e'的数据可以看出,模型计算值稍微偏大,所以减小半径为r的指数,就可以使模型计算值与实验值的偏差进一步减小。这里设其指数为0.3,则经验模型见式(11)。
针对以上数据作图,如图 5所示。
计算出e=0.003 378,从图 5中可以看出,模型计算值与实验值的拟合程度较好。其中最大拟合偏差为0.016 89。
用建立的模型计算Kyacal,i的值可以得到对应的脱除率θcal,i。由式(1)可知,
式中:Kya单位为h-1,根据经验模型计算得到的Kya单位为s-1。
表 3列出了传质的基础条件及脱除率实验测得θexp,i与由模型计算出的脱除率θcal,i。
根据表 3可知,θcal,i和θexp,i存在偏差,若两者之间的偏差很小,则可以认为所建立的经验模型是合理的。现计算偏差,结果如图 6所示。
从图 6可以看出,偏差很小,平均偏差为0.12%。模型计算出的传质系数值通过传质理论模型计算出的脱除率与实验得到的脱除率很接近。所以该经验模型是合理的。
采用络合铁法脱除天然气中H2S,以CH4和H2S的混合气模拟含硫天然气在超重力机中应用络合铁氧化还原脱硫实验数据,分析了超重力装置中原料气流量、脱硫溶液流量、超重机转子转速对气相传质系数的影响,建立了传质过程经验模型。该模型可以外推到其他体积与该超重力装置近似的装置,但气量应保持在~10 m3/h,液量在0.1~1.0 m3/h,转速在100~1 000 r/min。
该经验模型结合传质过程实验中的操作条件:气体流量、液体流量、转速,考虑填料不同位置对传质系数的影响。通过回归方法进行拟合,对拟合结果根据试验的方法进行修正,得到了经验模型。
该经验模型与实验值吻合程度较好,平均偏差为0.12%,与很多理论模型相比,该模型更为直观实用,从简单的操作条件出发,无需使用变化程度大的物性参数就能得到较高精度的传质过程相关数据。通过传质经验模型可以计算传质系数,从而对计算超重力脱硫技术的H2S脱除率,简化了反复进行实验测得脱除率的繁琐工作,有利于超重力脱硫技术研究工作的顺利开展。