上升流式生物流化床反应器膨胀行为的数学建模与分析.docVIP

上升流式生物流化床反应器膨胀行为的数学建模与分析 姜 虹，李 峰 (西安工业大学 计算机科学与工程学院，西安 7100) 摘要：本文基于Richardson-Zaki液固流态化方程，(方法)运用数值模拟分析的方法针对上升流式生物流化床反应器的膨胀行为展开了研究，并建立了相应数学模型。同时以厌氧膨胀颗粒污泥床反应器为例进行了数值模拟研究，(结果)数值模拟研究结果很好反映了该反应器的特点且相关实验结果的平均相对误差分别为4.74% 和8.3%。(结论)关键数值模拟；Richardson–Zaki方程；膨胀行为 中图分类号:TP302　　文献标识码:以其表面积巨大、传质高效、污泥负荷和容积负荷高、抗冲击能力强、生物活性好、占地面积少等优点[1]。一般来讲，生物流化床反应器中通过固定化与自固定化两种形式所形成的生物颗粒（或生物膜颗粒）具有良好的沉降性能以及很大的比表面积，从而能有效将微生物保留在反应器内，大大提高了单位容积的微生物持有量，进而极大地提高了反应器的水力负荷与处理能力。 根据载体流态化原理[2]，废水及气体自下而上流过生物颗粒床时, 生物颗粒与液体间将出现不同的相对运动，呈现出不同的工作状态。过度膨胀导致颗粒的大量流失，降低反应器处理能力与效率甚至导致反应器运行的失败。膨胀颗粒污泥床反应器ρ的废水以Vup的上流速度自下而上流动，废水从底部进入，从上部离开，如图1所示。 图1 上升流生物反应器示意图 Fig.2 Diagram of UBFR 实际情况下生物颗粒的沉降性能受到颗粒球形度、壁面效应等多种因素的影响。通过分析，假设颗粒的沉降过程符合以下条件：① 生物颗粒为球形；② 颗粒分散状态良好，颗粒之间互不干扰，为自由沉降过程。在此基础上根据经典的颗粒沉降理论，生物颗粒的沉降速度ut可按下式计算[5]： (1) 式中dp—生物颗粒直径(或当量球径)，m ρp—生物颗粒的密度 ，kg/m3 (—生物颗粒与流体的相对运动阻力系数 阻力系数(是流体相对于颗粒运动时的颗粒雷诺数Ret的函数，即 (2) 由于生物颗粒的沉降过程属过渡流状态（1＜Ret＜500)，故( 按下式计算[5]： (3) 将(3)代入(1)式中可得生物颗粒的沉降速度可按下式计算： (4) 式中μ—废水粘度，Pa·s；ρ—废水密度，kg/m3。 生物流化床的膨胀率与床层流体物性、流速、温度、通过床层的固体颗粒大小等因素有关。严格意义上讲，生物流化床属于气－液－固三相流化床系统，气体会对床层膨胀性能有着一定程度的影响。但由于气体（其密度远小于废水和生物颗粒的密度）对颗粒沉降和反应器膨胀性能影响很小，因此在研究其膨胀行为时可将其简化为液－固两相系统（生物颗粒－废水）来研究。 此时可运用Richardson－Zaki方程来确定废水流速和床层空隙率的关系。尽管该式是一个经验关联式，至今还没有为这一公式找到理论依据，但多年的实践证明，用该公式来描述液－固流化床的膨胀行为是十分准确的[6]。该方程具体表达如下： (5) 式中 ε—床层空隙率 u—颗粒相对于液体的运动速度，即液体表面上升流速Vup，m/h n—n只与Ret有关，两者之间存在以下关系[7]： 由于生物颗粒的沉降过程属过渡区域，故将n = 4.4Ret-0.1代入（5）式可得： (6) 膨胀率与空隙率之间的关系可如下表达[8]： (7) 其中ε0为床层的初始空隙率。对于近似球形颗粒，ε0一般取为0.4~0.45，本文ε0 取0.4。 将(6)式代入(7)式中得生物流化床膨胀率计算公式如下： (8) 根据(8)式可对不同液体表面上升流速Vup条件下不同生物流化床反应器的膨胀行为进行研究。 2 生物流化床膨胀过程的模拟分析 以EGSB反应器为例，运用上述所建立的数学模型进行数值模拟分析。由于工业条件下常见厌氧颗粒污泥的平均直径