机械蒸汽再压缩(mechanical vapor recompression,简称MVR)技术通过引入蒸汽压缩机回收利用二次蒸汽的潜热,提高了系统能源效率,相比传统蒸发技术可显著降低运行成本[1-3],其广泛应用于海水淡化、制盐、废水处理等领域,是国家重点推广的节能技术之一[4]。
作为国内科学研究和工业应用的重点项目,MVR技术已经进入快速发展时期[5]。Zhou等[6]以2%硫酸钠溶液为处理对象,建立起单效MVR稳态运行热力学模型,探讨了不同操作条件对系统运行参数的影响;李清方等[7]将MVR系统用于油田污水处理;Liang等[8]则对双效MVR系统进行了相关计算分析。
目前,国内对稳态运行时MVR系统的性能分析已有大量研究,而鲜有对于非稳态,即系统启动阶段的研究。MVR系统在热力学上属于朗肯循环,依靠压缩机的机械功作为循环的持续能源输入,在系统启动阶段一般需要提供启动热源以产生初始蒸汽。因此,启动阶段的非稳态运行过程对系统至关重要。本研究以西南地区某页岩气田压裂返排液为处理对象,在稳态计算的基础上建立起单效MVR系统的热力学模型,对系统非稳态运行过程进行了模拟与分析。
所用MVR系统的工艺流程图如图 1所示。
MVR系统启动时通常以外部通入的蒸汽作为热源,需要蒸汽发生器提供蒸汽。本研究旨在设计处理页岩气压裂返排液的撬装装置,为了简化流程,减少设备数量,该系统启动时以经过加热后的料液作为启动热源。流程中,蒸发器为竖直管降膜蒸发器,其基本结构参数为:换热管长度2.41 m,数量363根,正三角形排列;换热管外径25 mm,内径20 mm;折流板8块,间距263 mm,换热面积70 m2。
启动压缩机前,系统内为常压。开启进料泵与循环泵,向系统内通入料液并由电加热器加热至沸点,直至蒸发器与分离器内均存有一定量料液。启动压缩机后,蒸发器与分离器内压力下降,料液开始闪蒸,产生的蒸汽经压缩机压缩后进入蒸发器冷凝侧作为热源。当蒸发器与分离器内压力下降至稳态压力时,压力不再变化,料液停止闪蒸。
系统处理对象为经预处理后的页岩气压裂返排液,返排液采自西南地区某页岩气田某平台,其含盐质量分数约为2%,饱和浓度约为24%,可溶解固体主要成分为氯化钠,计算模型中料液参数均以相同浓度的氯化钠溶液参数替代,具体数据参照相关文献与手册[9-12]。其稳态运行条件见表 1。
目前,换热器瞬态,即非稳态计算过程多为无相变对流传热过程计算,其计算基本思想为对时间与空间进行微分,取换热管某一微元段建立起连续性方程与能量守恒方程,并将该方程与其他微元段所得方程联立,将得到的微分方程组进行离散化处理,利用迭代法计算得到任意时间与空间的各项参数的数值解。
有相变过程的计算分析常见于核反应堆热工水力模型,计算较为复杂,且其尚未涉及压力与浓度变化过程分析。本研究中非稳态计算过程主要考察在蒸发器结构参数确定的情况下压缩机最大进气量与最大功率等指标,蒸发器内各处传热过程随时间变化规律并非研究重点。故在换热器瞬态传热模型的基础上做出如下假设:①忽略蒸发器、预热器、分离器与管路热损失;②蒸汽压缩过程为绝热压缩过程;③忽略泵功率对系统的影响;④启动压缩机前, 蒸发器与管路内充满蒸汽,系统内不凝性气体被全部排出;⑤任意时刻蒸发器内浓缩液温度分布均匀,且为该时刻对应压力下的沸点温度;⑥任意时刻蒸发器内浓缩液浓度均匀分布。
换热器管内瞬态传热公式为[13]:
式中:下标1为管内,2为管外; 其他符号见文后符号说明, 下同。
消去方程两边dl,整理得:
预热器计算模型中,通过控制进料流量使经预热后的料液达到沸点,忽略系统内料液沸点温升的影响,则蒸发器为沸点进料,即在换热管内无单相对流换热区域,全部为饱和沸腾区,温度沿管长l不变,即:
由于设定温度分布均匀,在此基础上,可将微元体扩大至整根换热管,即dl=L,则M变为管内流体质量,式(3)变为关于时间t的一元微分方程。
将式(4)与式(5)线性化,可得:
在有相变过程中,M包含料液质量Mc与蒸汽质量Mv,即M=Mc+Mv, 则有:
式(7)则变为稳态下换热器热平衡方程。管外冷凝换热可参照上述简化过程。
为简化计算,模型中预热器并未进行详细设计,仅进行能量守恒计算。设定料液进料温度Ti为20 ℃,冷凝水出口温度T0为25 ℃,Tpr(i)为i时刻蒸发器内压力对应的料液沸点温度,其能量守恒公式为:
启动压缩机后,蒸发器内压强逐渐降低,压强变化规律参照真空系统抽气公式[14]。
式中:设$\frac{{V'}}{S} $值为100, 蒸发器所能达到的极限压力pw为68.65 kPa。
闪蒸放热量:
闪蒸产生的蒸汽质量:
蒸发器能量平衡公式如式(15)~式(18)所示:
蒸发器换热量可由式(19)计算:
总传热系数U(i)则由式(20)确定:
本文中冷凝对流传热系数采用式(21)计算,其中ts、twall分别为过热蒸汽温度与壁面温度。
沸腾传热系数采用Chen公式计算[15]。
压缩机压缩过程假设为绝热压缩过程,压缩机指示功采用式(22)计算:
式中:绝热指数k′取1.3。
MVR系统稳态运行时,以自身内部蒸发出的蒸汽经压缩后作为系统热源,除压缩机所需电能,无需外界额外提供能量。在非稳态运行时则需要外部提供运行所需启动热源,本研究主要目的为探究非稳态运行过程中蒸发量、压缩机功率等指标的变化规律,为实际操作提供一定参考。
在上述模型下,系统约在2 500 s后达到稳态运行,模拟结果见表 2。由表 2可知,非稳态模拟所得系统稳态运行结果与稳态计算结果相差不大。
在传热面积一定的情况下,传热温差与传热系数直接决定传热量大小,非稳态过程中二者变化规律如图 2、图 3所示。
由图 2、图 3可知,整个过程中二者均呈先上升后下降直至平稳的变化趋势。压缩机压缩比一定时,传热温差受蒸发器两侧压力与蒸发器内料液沸点温升影响,在启动初期,随蒸发侧压力降低,管内蒸发温度与管外冷凝温度差逐渐增大,传热温差逐渐增大直至闪蒸结束。之后,随运行过程中料液浓度的增大,料液沸点逐渐升高,逐渐减小了传热温差;传热系数则主要受冷凝侧对流传热系数影响,传热温差较低时,不利于蒸汽冷凝,故初期总传热系数较低,且传热系数随传热温差的增大而升高,闪蒸结束后,由于传热温差逐渐减小,传热系数也随之缓慢降低直至稳态。
蒸发量直接决定了非稳态过程中压缩机进气量与功率变化趋势,而系统内浓缩液浓度则是判别系统是否达到稳态的重要参数。
由图 4可知,蒸发量最初迅速增长,在348 s出现拐点,拐点处最大值为0.57 kg/s,后逐渐下降至稳态,在拐点之前,蒸汽包括传热产生的蒸汽与闪蒸产生的蒸汽两部分,随着蒸汽流量逐渐增加,蒸发量随之增加,当换热量达到该时刻蒸发器换热量的最大值时,传热量主要受蒸发器影响,在传热温差与传热系数均降低的情况下,蒸发量也随之下降,蒸发量最大值比稳态值高出32.56%;出料液含盐率与浓缩液含盐率相同,主要受蒸发量影响。如图 5所示,当2 000 s左右蒸发量基本不再变化时,出料液含盐率也随之基本稳定。
在压缩比一定的情况下,压缩机指示功会随进气量变化,非稳态过程中压缩机指示功变化趋势见图 6。
由图 6可知,压缩机指示功变化趋势与蒸发量变化趋势基本一致,模型中峰值功率出现在蒸发量最大时刻,约为45.90 kW,稳态时功率约为32.36 kW,峰值相比稳态值高出41.89%。实际运行过程中由于存在热损失等因素影响,峰值功率会有所下降。
(1) 建立起MVR非稳态模型并进行模拟,计算结果指示系统在约2 500 s时达到稳态,且达到稳态时的热力学性能与稳态计算结果相近。
(2) 蒸发量最大值与压缩机功率峰值均出现在非稳态运行过程中,最大功率值相比稳态值高出约41.89%,最大进气量相比稳态值高出约32.56%。
符号说明