1 引言
当前，一些学者对厌氧反应器流场和水力学性能进行了研究，例如，应用湍流模型模拟了侧伸搅拌式厌氧反应器流场，季军远等(2012)采取多釜串联及扩散模型模拟了分段组合式厌氧反应器行为，通过测定停留时间分布研究了厌氧反应器中不同反应区的流动特征.这些研究取得了一定成果，但由于模型相对简单，对于构造复杂的内循环厌氧反应器不能完全适用.同时，将动力学方程和流动模型结合以模拟厌氧反应器中反应过程的研究则较少有文献涉及.基于此，本文在建立厌氧反应动力学方程后，针对自主设计的内循环厌氧反应器采取组合模型进行模拟研究.厌氧反应过程中由于有气体生成，气体从产生到形成气泡，以及气泡在反应器中的向上运动必然会对整个反应器的流动型式产生重大的影响，使得反应器的流型偏离于理想状况.为此，本文通过测定实际反应器中流体的停留时间分布，根据停留时间分布的特征提出流动的组合模型，并据此研究实际的反应过程，以期为内循环厌氧反应器的进一步工业化放大提供有益的开发思路和方法.
2 实验方法
接种污泥为某啤酒厂UASB反应器中活性污泥 ;营养液参照产甲烷活性测定营养液配制方法进行配制;牲畜废水为牲畜粪便经过淘洗、过滤去除大颗粒泥沙及部分悬浮物后所得的液体，粪便取自某养猪场;示踪剂为硝酸钠(分析纯)，由杭州化学试剂有限公司提供.
主要仪器：FA1104N型电子天平、DFG-9053A型电热恒温鼓风干燥箱、TC-15型恒温电热套、SC-15型数控超级恒温槽、BT01-YZ1515型蠕动泵、GC-1690B型气相色谱仪
2.1 内循环厌氧反应器主反应区停留时间分布密度函数E(t)的实验测定流程
为建立厌氧反应器的流动模型，实验研究了反应器停留时间分布，并考察了不同气体流量的影响，获得了停留时间分布函数.内循环厌氧反应器由主、次两反应区组成，由于反应主要在主反应区4进行，因此，实验主要针对反应区4进行停留时间分布测定.实验流程如图 1所示.
图1 工艺流程示意图(1.进水罐，2.蠕动泵，3.气泵，4.主反应区，5.内反应区，6.出水口，7.气体出口，8.反应液出口，9.环隙反应区)
以水为工作介质，将水加入进水罐至一定高度，开启蠕动泵，保持水流量为15 L · h-1，流量计安装在泵出口处;同时开启气泵，气泵流量由安装在泵出口处的气体流量计检测，分别为0、90、180 mL · h-1.水从下而上流过反应器，反应器分为上下两部分，下部为主反应区，内径为110 mm，高550 mm，几何体积为5 L;上部为次反应区，又分为两部分，一是内反应区5，另一为环隙反应区9，环隙反应区几何体积VR4为2.44 L，内反应区几何体积VR5为2.86 L.在环隙反应区的上部安装有气体出口管，由于测定的是主反应区停留时间分布，故液体在主反应区上部出水口处排出，出水口处同时安装一电导仪，通过检测出口液体的电导率而推算出示踪物浓度的高低.
当气液两相流动达到稳定后，在进水罐上方快速加入一定量的示踪剂，同时计时，间隔一定时间，于反应器出口处记录下水的电导率随时间的变化.
2.2 验证实验工艺流程
验证实验工艺流程示意图如图 2所示.新鲜料液加入储槽中，料液由蠕动泵输送，进入换热器中，换热器由10根有机玻璃管组成，管外径为10 mm，长65 mm，料液走管程，与循环水换热后进入反应器，在次反应区顶部进行气液分离，液体进入废液槽，气体经过气体流量计后由缓冲瓶进入液封槽，液封槽容积为1500 mL，最后进入皂沫流量计计量放空.
图2 验证实验工艺流程示意图(1.原料液储槽，2.蠕动泵，3.流量计，4.换热器，5.反应器，6.反应废液槽，7.气体流量计，8.液封槽，9.缓冲瓶，10.皂沫流量计)
由于新鲜料液配料较为困难，当系统稳定，即气体流量变化不大时，每隔4 h分别测量气液两相组成、流量，一般测量2～3次结束实验.
2.3 内循环厌氧反应器主反应区流动模型建立
对反应器主反应区流动状况作出如下假设：示踪物A为瞬间脉冲进料;反应器中气体流量恒定不变，且均匀上升;进料清水流量亦恒定不变.现假设厌氧反应器主反应区流动型态可由一平推流区(PFR)与两个并联的全混区(CSTR)串联组合而成，模型示意图如图 3所示.
图3 中反应器流动模型组合示意图文名
设反应器主反应区总体积为VR，等于模型中各反应区体积VRi之和，设PFR的体积分率为f1，CSTR2、CSTR3的体积分率分别为f2、f3，总流量为v0，流经CSTR2的流量分率为I，某时刻流出CSTR2、CSTR3全混区的流体中示踪剂浓度分别为CA2、CA3，混合后出口总浓度为CA.当t=0时，在PFR反应器入口处，脉冲法加入质量为M的饱和示踪剂A.对于PFR：
式中，VRP为PFR的体积(L); p为流体在PFR中的平均停留时间(min).
对CSTR2进行物料衡算，求取CSTR2出口浓度CA2：①当t< p时，由于示踪剂A仍在PFR中，对CSTR2，CA2=0;②当t≥t¯t¯p时，对CSTR2作物料衡算，根据输入量=输出量+累计量，则有：
求解该方程，得
同理，对CSTR3进行物料衡算，可求取CSTR3出口浓度CA3：①当t< t¯t¯p时，由于CSTR3中无示踪剂，故CA3=0;②当t≥ t¯t¯p时，对CSTR3作物料衡算，可得：
最后，对CSTR2、CSTR3出口处作物料衡算，由于该处进行的是混合过程，则：①当t< t¯t¯p时，CA=0;②当t≥t¯t¯p时，CA=ICA2+(1-I)CA3.
停留时间分布密度函数E(t)的求取公式如下：
当t . 最后：
对模型求解转化为求解一个三参数I、f1、f2的E(t)方程，其中，t¯t¯p不是独立的. 3 实验结果及模型参数求解 3.1 主反应区停留时间分布测定 实验是在清水流量为15 L · h-1，气体平均流量分别为0、90.0、182.1 mL · h-1条件下进行测定的，示踪剂脉冲加入，部分实验数据如表 1所示. 根据
，各时间点的E(t)值如表 2所示.应用非线性最小二乘法求解，不同气体流量时的模型参数及方程如表 3所示. 表2 停留时间分布密度函数E(t)随时间变化关系 当vg=0、90.0、182.1 mL · h-1时，分别获得E(t)方程.E(t)方程获得后，将不同时刻t时的CA值，代入方程，可获得E(t)的计算值，其值列于表 2中.由表 3中数据可知，随着气体流量的增大，反应器中存在较大的返混，PFR体积所占比率f1不断减小，说明气体流量对反应器的流型有着较大的影响.
3.2 内循环厌氧反应器模拟 3.2.1 内循环厌氧反应器反应过程模型的建立及验证 1)反应器反应过程模型建立 根据以上停留时间分布的特点，厌氧反应器由主、次反应区组成，因此，可用图 4中的组合模型模拟反应器反应过程行为.图 4所示组合模型中，根据停留时间分布规律，1、2、3拟合的是厌氧反应器主反应区;次反应区由中间回流区和环隙区构成，中间回流区可作为平推流处理，环隙区由于有大量产物气体流过，可作为全混区处理，所以反应器次反应区可由4、5并联进行拟合. 图4 反应器反应过程模型组合示意图 模型中各参数的确定：针对本反应器，实际反应过程中生成的气量在90～180 mL · h-1之间，因此，模型中主反应区各参数取其平均值，具体如下：f1=0.36，f2=0.53，f3=0.11，I1=0.53. 次反应区流量分率为：
各反应区有效体积VRi的确定：总污泥加入量为5 L，上层即次反应区污泥所占分率按0.35计算，则：VR1+ VR2+VR3=3.25L，VR1=1.17L，VR2=1.7225L，VR3=0.3575L，VR4+ VR5=1.75L，VR4=V4/(V4+ V5)=0.8057L，VR5=0.9443L. 2)模型验证 为验证上述模型是否正确，进行了3组不同初始浓度操作状态的实验，当厌氧反应器操作稳定后，即在一定时间能气体流量基本不变，此时，操作参数分别为：废液体积流量v0L=15 L · h-1，温度T=35 ℃，废液初始COD CS0=3.3840、3.5326、3.8729 g · L-1，由前期实验所得动力学方程为：
式中，Cs为反应器中废液浓度(g · L-1)，-rs为废液COD消耗速率(g · min-1 · L-1)，rpCH4为甲烷生成速率(g · min-1 · L-1)，rpCO2为二氧化碳生成速率(g · min-1 · L-1)，模拟计算了反应器单程转化过程，主、次反应区出口即B、C两点处计算值如表 4所示.以基质初始浓度CS0=3.3840 g · L-1为例，实验数据及处理如表 5所示. 表4 主、次反应区出口处B、C两点状态数据 表5 验证实验各点状态数据一览 CH4体积流量vCH4实为193.81 mL · h-1，CH4质量流量FCH4实为2.045×10-3 g · min-1.同理，可计算CO2体积流量、质量流量分别为8.75 mL · h-1、2.54×10-4 g · min-1.误差按(实验值-计算值)/实验值×100%计算，可知本组B点基质浓度CSB误差最大，为7.86%. 由表中数据可见：不同时刻反应器出口C点的气体体积流量v大于模拟计算时的体积流量vgC，其原因可能是该验证实验时间不够长，系统中存在一定的N2，同时反应过程中生成其它气体所致.其余两组数据计算依此类推，最后得到不同浓度时，C点出口处CH4体积流量最大误差为9.67%，CO2体积流量最大误差为8.53%，可见所建立的组合模型能够反映反应器实际操作状况. 3.2.2 内循环厌氧反应器模拟计算 根据以上反应过程模型，模拟了初始COD、温度变化对厌氧反应器状态的影响. 1)初始COD变化对反应器操作状态影响的模拟研究 操作参数为：废液体积流量v0L=15 L · h-1，温度35 ℃，改变COD初始值，B、C两处模拟数据如表 6所示.由表 6中数据可知，当温度一定时，随着基质初始浓度的增高，B、C两处基质浓度也随之升高，产物CH4、CO2的质量流量也升高. 表6 反应器出口处B、C两点状态随基质浓度变化关系 2)反应温度变化对反应器操作状态影响模拟研究 操作参数为：废液体积流量v0L=15 L · h-1，COD初始值CS0=3.40 g · L-1，改变温度，B、C两处模拟数据如表 7所示.由表 7中数据可知，在基质初始浓度一定时，随着温度的升高，B、C两处的基质浓度不变，这主要是由于在该温度范围内，基质消耗活化能很低，以致消耗速率几乎不随温度变化而形成的结果;但产物CH4、CO2的质量流量随温度的升高而升高，同时CO2的增加速率要快于CH4，这主要是由于CO2生成活化能高于CH4，故其生成速率对温度更为敏感. 表7 反应器出口处B、C两点状态随基质温度变化关系 4 结论 1)通过对反应器主反应区停留时间分布的考察，主反应区流型可以用一平推流区与两个并联的全混区进行串联组成，各区所占体积分率随反应过程中所产气体流量不同而不同，随着气体流量的增大，平推流区所占体积分率减小. 2)次反应区由平推流区和全混区并联组合而成，整个厌氧反应器由主反应区和次反应区串联组合而成，根据实验数据可知，C点出口处CH4体积流量最大误差为9.67%，在可允许的范围内，说明该组合模型较好地模拟了反应器的实际状态. 3)采取单因素分析方法，模拟了基质初始COD和温度变化对反应器出口状态的影响，发现产物CH4、CO2的质量流量随基质浓度、温度的升高而升高;在304.15～312.15 K范围内，温度的变化对出口处基质浓度的影响不大.
表1 不同气体流量时分布电导率随时间变化关系
