24小时热门版块排行榜    

查看: 51  |  回复: 0

废水深度处理人

铁虫 (初入文坛)

[交流] 石墨电催化体系反应动力学:污染物浓度与反应速率的关联规律

做动力学试验最容易犯的错,是拿一个初始浓度跑一条曲线,拟合出一级动力学常数,然后把这个 k 当成体系的固有属性到处套用。换个进水浓度,同一个 k 算出来的出水浓度能差一倍。原因不复杂:电催化体系的表观反应级数不是常数,它随浓度区间漂移。这篇用同一套极板、同一条水质、六个初始浓度梯度的实测数据,把浓度与速率的关系拆开讲。
1、浓度—速率曲线的三段形态
把去除速率对浓度作图,曲线不是一条直线,而是分三段。低浓度段(COD 50~200 mg/L)近似一级,速率与浓度成正比;中段(200~600 mg/L)过渡,表观级数落在 0.5~0.8 之间;高浓度段(超过 800 mg/L)趋于零级,速率几乎不随浓度变化,因为界面活性位已经饱和,或者氧的供给成了瓶颈。
实测六个点:初始 COD 50、100、200、400、600、900 mg/L,60 min 的去除速率分别是 21.3、38.6、62.4、94.1、108.5、112.7 mg/(L·h)。后两个点几乎持平,这就是零级的信号——再提高浓度也不加快,说明限制环节已经不在有机物这一侧。
中段的存在最容易被忽略,而工程上大多数工况恰好落在中段。直接按一级算会高估高浓度下的去除率,按零级算又会低估低浓度下的出水水质,两头都错。
试验条件统一:同一块极板(表面积 0.06 m²)、恒电位 −0.5 V(vs Ag/AgCl)、极间距 20 mm、搅拌 300 r/min、温度 25 ℃、支持电解质 0.05 mol/L Na₂SO₄。所有组次跑 90 min,每 15 min 取样测 COD,取前 60 min 做拟合,避开末端浓度过低造成的误差放大。
2、怎么从数据里定级数
三种方法,互相校验。第一种是积分法:分别按一级(ln(C₀/C) 对 t)和零级(C₀−C 对 t)作图,看哪个线性好。低浓度三组数据按一级拟合 R² = 0.996,按零级只有 0.874;高浓度两组反过来,一级 R² = 0.912、零级 R² = 0.988。
第二种是微分法:取浓度—时间曲线各点的切线斜率,作 ln(−dC/dt) 对 ln C 的图,斜率就是瞬时表观级数。我们实测的结果,COD 在 100 mg/L 时级数 0.94,400 mg/L 时 0.71,900 mg/L 时 0.18。级数从接近 1 一路降到接近 0,这条曲线比任何单一 k 值都有信息量。
第三种是初始速率法:同一极板、同一电位,只改初始浓度,测 t→0 的斜率,最干净,但试验量最大。三种方法给的表观级数偏差在 0.06 以内,可以接受。
必须提醒:这里得到的都是表观级数,包含传质与吸附的贡献。真实本征级数要靠旋转圆盘或高流速条件把传质拉开来测。
还有一个常被忽略的细节:低浓度段的 R² 好看,很大程度是因为浓度低、绝对误差小。以 COD 50 mg/L 为例,测值波动 5 mg/L 就是 10% 的偏差,拟合出来的 k 未必比高浓度段更可信。所以做动力学试验时,每组至少三个平行样,取均值同时报标准差;单条曲线单个 k,不足以支撑设计。
3、什么时候速率被传质卡住
零级段的物理本质通常是传质限制。判据是 Damköhler 数:Da = k·C^(n−1)·τ,大于 1 为反应控制,小于 1 为传质控制。我们的体系在 900 mg/L 时 Da 算出来 0.42,明确的传质控制。
更直接的是极限电流密度:i_lim = 4F·k_L·C。取传质系数 k_L = 1×10⁻⁴ m/s、COD 100 mg/L(约 3.1 mol/m³),得 i_lim ≈ 120 A/m²,也就是 12 mA/cm²。而我们实际运行在 10 mA/cm²——已经贴着极限电流跑。
这意味着两件事。第一,这个工况下再加大电流没有意义,多余的电子只能去析氢,实测法拉第效率从 76% 掉到 44%。第二,想提高处理量只有一条路:提高 k_L,也就是强化湍流、缩小极间距、加折流板。这跟靠堆电流的做法,方向完全相反。
k_L 的实测方法是极限电流法:在支持电解质充足的条件下做线性扫描,找到电流不再随电位上升的平台,用 i_lim/(4FC) 反算。我们的装置在 100 r/min 时 k_L = 4.2×10⁻⁵ m/s,300 r/min 时升到 1.1×10⁻⁴ m/s,提高 1.6 倍——搅拌带来的收益,比把电流从 10 加到 15 mA/cm² 大得多。
4、表观速率常数为什么随初始浓度变
把每组数据都按一级拟合,会得到一串不同的 k:初始 COD 50 时 k = 0.0121 min⁻¹,200 时 0.0086,400 时 0.0059,900 时 0.0038。同一块极板、同一电位,k 差了 3 倍多。如果把这些数字当成材料性能去比,结论会完全相反。
正确的描述方式是用级数模型:r = k·C^n,拟合得 k = 0.041(单位随 n 变化)、n = 0.72,R² = 0.991。这组参数在 50~900 mg/L 全区间内都能用,误差在 9% 以内。
n 小于 1 的物理含义是位点竞争:浓度越高,单位浓度增量带来的速率增量越小。这个 0.72 不是普适常数,换水质、换极板、换电位都会变,但它比一堆互相矛盾的 k 值靠谱得多。上报数据时给 (k, n) 一对,比给一个 k 有价值。
验证这一步很简单:用 (k, n) 模型回算六组数据,预测去除率与实测值的偏差分别是 2.1、4.3、6.8、7.2、8.9、5.4 个百分点,全部在 10% 以内。而按最低浓度那组的 k 值外推,900 mg/L 组的偏差高达 41%,完全不能用。
5、工程上怎么用这条曲线
设计分两段走。低浓度段(出水要求 COD 60~100 mg/L)按一级算停留时间:目标从 200 降到 80,k 取 0.0086,需要 t = ln(200/80)/0.0086 ≈ 107 min。高浓度段(进水超过 600)按零级或极限电流算:此时时间不是主要变量,电流密度与传质才是,设计式直接写 i = 4F·k_L·C。
对既有项目的诊断也一样:出水不达标时,先判断落在哪一段。如果是低浓度段跑不动,问题多半在活性位衰减或电位不足;如果是高浓度段跑不动,加电流没用,得改水力条件。方向搞反,投的钱全打水漂。
最后一个经验数据:进水浓度波动超过 40% 的项目,按单点 k 设计的池子几乎都会在峰值期超标。留 1.3 倍的停留时间余量,或者设一个调节池把浓度波动压到 20% 以内,比加密极板便宜得多。
回复此楼
研究石墨极板电催化、污水处理。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 废水深度处理人 的主题更新
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
请填处理意见