利用生物质灰渣批量吸附Cr(VI)的放大研究：Pyomo与SymPy应用案例

《Scientific African》：Upscaling of batch adsorption of Cr(VI) using biomass ash: A case for the use of Pyomo and SymPy

　　

研究人员利用外来入侵物种Prosopis树作为锅炉生物质燃料。产生的生物质灰渣带来了新的处置问题。为满足循环经济需求，研究人员提出研究将生物质灰渣用作水溶液中重金属吸附剂的可能性，并考察了负载金属后的生物质再生能力。以Prosopis树生物质灰渣为吸附剂，处理含Cr(VI)离子的合成废水。采用中心复合可旋转设计（CCRD）改变吸附时间、固载量、温度和初始金属离子浓度。Cr(VI)在生物质灰渣上的吸附遵循Freundlich等温线和准二级动力学。负载金属的生物质灰渣可循环再生4次，之后解吸效率下降至59%，对应吸附率骤降至77%。通过量纲分析获得的无量纲时间π组用于放大研究。Python优化建模对象（Pyomo）预测的放大结果与实际放大实验相比准确率达98.6%。研究表明，Prosopis树灰可作为Cr(VI)吸附剂，Pyomo可通过量纲分析生成π组作为放大工具。

**研究背景与问题**
工业废水中重金属Cr(VI)的去除是环境治理的重要课题。传统吸附剂成本较高，而生物质灰渣作为一种固体废弃物，兼具资源化利用和污染控制潜力。Prosopis树作为非洲干旱地区的外来入侵物种，被用作锅炉生物质燃料，但其燃烧产生的灰渣目前多被填埋，造成二次污染。为满足循环经济需求，研究人员提出将这种灰渣转化为廉价吸附剂，用于去除水溶液中的Cr(VI)离子。然而，大多数批式吸附研究仅停留在实验室规模，缺乏可靠的放大方法，限制了其工业应用。现有放大方法如构建中试装置成本高昂，而数学模拟模型缺乏通用性。因此，本研究旨在评估Prosopis树生物质灰渣的吸附性能及再生能力，并借助量纲分析与Pyomo优化工具开发一种可推广的放大模型，促进实验室成果向实际工业转化。该论文发表在《Scientific African》。

**主要技术方法**
研究人员采用中心复合可旋转设计（CCRD）进行批式吸附实验，考察固载量、时间、温度和初始Cr(VI)浓度四个参数的影响。通过X射线荧光光谱（XRF）分析元素组成，傅里叶变换红外光谱（FTIR）鉴定化学键，扫描电子显微镜（SEM）观察表面形貌，并采用固体添加法测定零点电荷（PZC）。吸附数据拟合至Langmuir与Freundlich等温线以及准一级、准二级和颗粒内扩散动力学模型。再生实验使用0.01 M HCl解吸负载金属的灰渣。放大模型基于Buckingham π定理，利用SymPy库生成无量纲π组，并通过Pyomo调用内点优化求解器（IPOPT）进行求解，经验证后预测放大规模下的吸附性能。

**研究结果**

*Characterisation of the biomass ash（生物质灰渣表征）*
XRF分析显示Ca和Mg为主量元素（质量分数分别为46.95%和6.68%），烧失量24%。激光衍射测得比表面积为420.7 m2/kg，表面平均粒径11.5 μm，体积平均粒径34.9 μm。PZC为9.0，表明pH低于9时吸附剂表面带正电。SEM观察显示灰渣呈高度团聚的不规则形态，表明其为无定形多孔材料。FTIR谱图在3644 cm?1处出现O–H伸缩峰，1407.7 cm?1为O–H弯曲，1117.3 cm?1为C–O伸缩，875.54 cm?1为C–H弯曲。

*Effect of parameters on solid loading（参数对固载量的影响）*
通过相关性热图分析发现：残余浓度与固载量呈中等负相关（固载量增加，吸附位点增多，去除率升高）；与初始浓度呈中等正相关（高浓度下竞争加剧导致残余浓度升高）；与温度和时间近乎零相关，表明吸附速度快且时间延长可能导致已吸附Cr(VI)解吸。

*Variation in pH with time（pH随时间变化）*
pH随时间增加（从约5升至7），归因于Ca和Mg的溶出，但始终低于PZC（9.0），有利于Cr(VI)的吸附。

*Regression modeling（回归建模）*
多元线性回归获得显著模型（Eq.9），仅保留固载量、初始浓度、固载量平方项及固载量与浓度交互项。模型相关系数R2=0.96，调整R2=0.936，F检验显著（p=1.37×10??），Durbin-Watson统计量2.3（无自相关），多重共线性测试值6.06（无共线性），均方根误差5.12。系数置信区间均不含零，验证模型有效性。

*Response surface（响应面）*
三维响应面显示：增加固载量显著降低残余浓度，但存在边际递减效应；温度影响微弱；初始浓度升高导致残余浓度急剧上升，尤其在低固载量时。最优固载量为0.8 g/L，处理100 ppm Cr(VI)溶液可实现98%以上去除率。

*Adsorption isotherms（吸附等温线）*
Freundlich等温线拟合度（R2）高于Langmuir等温线，且Langmuir常数较小表明吸附较弱。Freundlich参数K_f=77.05（mg/g）（L/mg）^1/n，n=6.93，1/n<1，表明吸附容量大且为优惠吸附，吸附表面非均质。

*Adsorption kinetics（吸附动力学）*
颗粒内扩散模型R2低（0.69）且直线不通过原点，不适用。准一级和准二级模型R2较高，但准二级模型计算的平衡吸附容量q_e与实验值绝对误差仅0.05，而准一级误差大于5，故Cr(VI)吸附遵循准二级动力学，表明化学吸附占主导。FTIR谱图对比显示峰位移动（>10–30 cm?1）及新峰出现（<800 cm?1），证实了表面配位键（M-O键）的形成，进一步支持化学吸附机制。

*Reusability（可重复使用性）*
生物质灰渣在3次循环再生中吸附容量无明显下降；第4次循环后解吸效率降至59%，对应吸附率降至77%。酸再生导致Ca含量从46.95%降至12.58%，Ca2?在溶液pH调节中消耗，是吸附容量下降的主因。解吸效率比吸附效率对老化更敏感，可作为吸附剂失效的预警指标。

*Upscaling（放大）*
通过Buckingham π定理和SymPy得到四个无量纲π组：π?=C?×K_L（亲和力负载），π?=q_e/q_max（材料常数），π?=S/C?（固液比），π?=k?t/q_max（无量纲时间）。归一化吸附容量与无量纲时间图可用于预测不同规模下的吸附性能。放大时保持π?和π?恒定，以维持吸附驱动力和动力学机制一致。Pyomo优化输出不同规模下的吸附剂质量（如处理25 ppm、10 L废水需4.16 g）。验证实验（10 L溶液，25–100 ppm Cr(VI)）的实际去除率（96.41%–91.23%）与预测（95%）高度吻合，t统计量（1.26）小于临界值（3.18），p=0.29，表明预测与实验结果无显著差异，放大模型准确率98.6%。

**总结与结论**
研究证实，Prosopis树生物质灰渣对Cr(VI)的吸附遵循Freundlich等温线和准二级动力学模型，非均质表面化学吸附为主导机制。通过量纲分析和Buckingham π定理可生成无量纲π组，结合Pyomo优化可实现批式吸附过程的高精度放大（准确率98.6%），为实验室研究向工业规模转化提供了经济有效的工具。研究局限性在于使用合成废水而非实际工业废水，且放大验证仅基于有限条件，未来需在多种工业条件下进一步测试。研究结论如下：本研究强化了吸附理论，证明Cr(VI)去除符合Freundlich等温线（非均质表面吸附）和准二级动力学（化学吸附主导机制）。量纲分析与π组的成功应用进一步证实其连接实验室结果与工业规模的相关性。Pyomo建模的整合凸显了计算优化在过程放大中的重要性。从实践角度看，Prosopis生物质灰渣是一种可行、低成本且可持续的重金属废水吸附剂；可多次再生的特性增强了其经济性，而高精度的放大模型展示了其在真实工程系统中的适用性。关键启示包括：入侵生物质（如Prosopis）可价值化为符合循环经济原则的功能材料；再生对可持续性至关重要，但设计时需考虑性能衰减。