基于羟基自由基的高级氧化技术降解甲硝唑：动力学、机制及天然水体基质影响研究
一、 研究团队与发表信息
本研究由来自中国电力建设集团中南勘测设计研究院有限公司的苏荣奎和戴祥荣（并列第一作者）、中南林业科技大学环境科学与工程学院的王汉青、王志祥、李子实、陈永华、欧阳丹霞，以及湖南第一师范学院的罗艺婷（通讯作者）共同完成。研究成果以题为《Metronidazole degradation by UV and UV/H2O2 advanced oxidation processes: kinetics, mechanisms, and effects of natural water matrices》的论文形式，于2022年9月28日发表于国际期刊《International Journal of Environmental Research and Public Health》（IJERPH 2022, 19, 12354）。
二、 学术背景与研究目的
本研究属于环境科学与工程领域，具体聚焦于水处理技术中的高级氧化工艺（Advanced Oxidation Processes, AOPs）。研究的背景源于日益严峻的水环境中抗生素残留问题。抗生素作为一类与人类活动密切相关的有机污染物，通过各种途径进入环境水体，因其化学性质稳定、残留时间长、难生物降解等特点，对生态系统和人类健康构成潜在威胁。传统的污水处理工艺（如活性污泥法）难以有效去除这类新兴污染物。因此，开发高效去除残留抗生素的技术成为亟待解决的难题。
以羟基自由基（•OH）为主要活性物质的高级氧化技术，因其氧化能力强、反应快速、无选择性等优点，在降解难降解有机污染物方面展现出巨大潜力。其中，紫外光活化过氧化氢（UV/H2O2）技术因其操作简便、环境友好、应用前景广阔而备受关注。然而，该技术在应用过程中仍面临污染物降解机制难以确定、复杂水体基质影响难以预测等问题。
本研究选取甲硝唑（Metronidazole, MTZ）作为代表性抗生素模型污染物，因其使用量大且环境相关性高。研究旨在系统性地比较MTZ在单独UV和UV/H2O2体系中的降解行为和机理。具体目标包括：1）测定光化学反应系统的有效光强和有效光程；2）测量MTZ的直接光解参数（摩尔吸光系数和量子产率）；3）对比MTZ在两种体系中的降解动力学；4）通过竞争动力学法和基于稳态假说的动力学模型，测定MTZ与•OH的二级反应速率常数（k_MTZ,•OH）及•OH的稳态浓度（[•OH]_ss），并量化直接光解和•氧化的贡献率；5）研究H2O2投加量和复杂水体基质（如有机物、无机阴离子）对MTZ降解的影响。研究成果对于理解抗生素的氧化降解机制、预测其在复杂环境中的归趋，以及优化基于•OH的高级氧化工艺在实际水体（如河流、湖泊、地下水）中的应用具有重要的理论和实践意义。
三、 详细研究流程与方法
本研究包含一系列紧密衔接的实验与理论分析步骤，流程严谨，方法学基础扎实。
第一步：光化学系统表征。 这是所有后续光化学实验定量化的基础。研究采用了两种液相化学光测定法来精确测量反应系统的关键参数。首先，采用铁草酸钾化学光测定法。该方法基于Fe³⁺在酸性条件和光照下被还原为Fe²⁺的反应（2[Fe(C2O4)3]³⁻ + hν → 2[Fe(C2O4)2]²⁻ + C2O4²⁻ + 2CO2），生成的Fe²⁺与1,10-邻菲啰啉形成红色络合物，在510 nm处测定吸光度，通过标准曲线计算Fe²⁺浓度。根据生成的Fe²⁺量与光照时间的线性关系（拟合方程斜率k0），结合已知的量子产率（φ = 1.25 mol Einstein⁻¹ @ 254 nm），利用公式 I0 = k0/φ 计算得到系统的有效体积平均光强（I0）。其次，采用过氧化氢化学光测定法进行交叉验证。该法基于H2O2在254 nm光下的光解。对于高浓度H2O2（满足2.303ελlc > 2），其光解遵循零级动力学，光解速率公式简化为 -d[H2O2]/dt = φλ I0，由此可直接计算I0。对于低浓度H2O2（满足2.303ελlc < 0.02），其光解遵循一级动力学，速率常数k_obs与I0和有效光程（l）相关：k_obs = 2.303 φλ I0 ελ l。通过测定低浓度H2O2的一级光解速率常数，并结合已知的ελ和φλ，可反算出有效光程l。这两种独立方法测得的结果相互印证，确保了参数的可靠性。
第二步：甲硝唑直接光解特性研究。 在表征好的光化学系统中（I0 = 7.50 × 10⁻⁶ Einstein L⁻¹ s⁻¹， l = 0.935 cm），研究MTZ在单独UV照射下的降解行为。将初始浓度为10 μM的MTZ溶液（pH 3.00，磷酸盐缓冲体系）置于系统中进行光照实验，定时取样并通过超高效液相色谱（UPLC）分析MTZ浓度。结果表明MTZ的降解符合准一级反应动力学。通过测定MTZ溶液在254 nm处的吸光度，根据朗伯-比尔定律（A = ε [MTZ] l）计算其在该波长下的摩尔吸光系数（ε）。进一步，结合测得的直接光解初始速率（r_uv）、有效光强I0、摩尔吸光系数ε和光程l，利用公式 φ_MTZ = r_uv / { I0 [1 - 10^(-ε l [MTZ])] } 计算MTZ在254 nm下的量子产率（φ）。这一步骤量化了MTZ自身对紫外光的吸收和转化能力。
第三步：UV/H2O2体系中MTZ降解动力学与机理探究。 此部分是本研究的核心。首先，在UV体系中加入100 μM H2O2，研究MTZ的降解动力学，发现其降解速率显著提升，为准一级动力学。通过暗对照实验排除了H2O2直接氧化的贡献。为了证实•OH的产生，研究使用了羟基自由基探针分子对苯二甲酸（TPA）。TPA与•OH反应生成具有强荧光的2-羟基对苯二甲酸（HTPA），通过监测荧光强度的增长，直接证明了UV/H2O2体系中•OH的生成。为了定量MTZ与•OH的反应活性，采用了两种方法：1）竞争动力学法：以4-氯苯甲酸（pCBA）作为参比物（其k_pCBA,•OH已知），将MTZ与pCBA置于同一UV/H2O2体系中竞争•OH。通过测定两者降解速率常数的变化，根据公式 k_MTZ,•OH / k_pCBA,•OH = (k_tot,MTZ - k_uv,MTZ) / (k_tot,pCBA - k_uv,pCBA) 计算k_MTZ,•OH。2）基于稳态假说的动力学模型法：这是本研究方法学的亮点。研究者建立了一个包含•OH生成、消耗（与MTZ、H2O2、缓冲盐组分等反应）的详细动力学模型（模型涉及反应及速率常数见表2）。在稳态假设下（•OH的生成速率等于消耗速率），通过联立求解方程，可以同时计算出k_MTZ,•OH和体系中的[•OH]_ss。模型输入参数包括：测得的有效光强I0、光程l、H2O2和MTZ的摩尔吸光系数与量子产率、MTZ的总降解速率常数（k_tot）和直接光解速率常数（k_uv），以及其他背景物质的浓度和已知的速率常数。最后，利用计算出的k_MTZ,•OH和[•OH]_ss，结合公式 r_tot = r_uv + r_•oh = k_uv[MTZ] + k_MTZ,•OH [•OH]_ss [MTZ]，可以量化直接光解和•OH氧化对MTZ降解的贡献率。
第四步：操作条件与水体基质影响研究。 利用建立的稳态动力学模型和实验验证，系统考察了H2O2投加量以及天然水体中常见组分对MTZ降解的影响。1）H2O2投加量：在50-400 μM范围内改变H2O2初始浓度，通过实验测定MTZ的降解速率常数（k_exp,obs），并与模型预测值（k_cal,obs）进行比较。模型还能分别预测不同H2O2浓度下，直接光解（k_cal,uv）和•OH氧化（k_cal,•oh）的贡献率。2）有机质：以腐殖酸（HA）作为天然有机质（NOM）的代表，研究其浓度（0-2.0 mgC L⁻¹）的影响。通过实验测定k_exp,obs，并与模型预测比较。模型通过设定腐殖酸的吸光系数（ε_NOM）和其与•OH的二级反应速率常数（k_NOM,•oh），可以模拟其两种抑制作用：竞争紫外光吸收和淬灭•OH。通过对比忽略其中一种作用时的模型预测值，可以评估两种抑制机制的相对重要性。3）无机阴离子：研究了Cl⁻、SO₄²⁻、NO₃⁻（0-5 mM）在pH 3.00条件下的影响。由于这些阴离子在254 nm处吸光度极低，对直接光解无影响，因此它们主要通过影响•OH及其后续反应来发挥作用。实验测定了添加不同阴离子后的k_exp,obs，并与不含阴离子的对照组进行比较。
数据分析：所有实验均进行三次重复，数据以均值±标准差表示。使用Microsoft Office Excel 2016进行统计分析，使用OriginPro 2018绘制图表。动力学拟合采用线性回归方法。
四、 主要研究结果与逻辑关系
本研究获得了一系列定量化的关键结果，层层递进，逻辑严密。
首先，光化学系统参数测定结果显示，铁草酸钾法和过氧化氢法测得的有效光强高度一致，为7.50 (±0.23~0.25) × 10⁻⁶ Einstein L⁻¹ s⁻¹。过氧化氢法测得的有效光程为0.935 (±0.04) cm。这为后续所有光化学反应的定量计算提供了精确的基准。
其次，MTZ直接光解特性结果表明，在给定光强下，MTZ的直接光解初始速率仅为0.147 μM min⁻¹。究其原因，测得MTZ在254 nm处的摩尔吸光系数为2645.44 M⁻¹ cm⁻¹，处于较低水平；其量子产率为5.9 × 10⁻³ mol Einstein⁻¹，也较低。这两个参数共同决定了MTZ对254 nm紫外光的吸收能力和光转化效率都不高，从而导致其直接光解速率缓慢。这一结果解释了为何需要借助高级氧化工艺来强化其降解。
第三，UV/H2O2体系降解动力学与机理结果显示，加入100 μM H2O2后，MTZ的降解速率大幅提升至1.47 μM min⁻¹，是单独UV时的10倍。TPA荧光实验直接证实了体系中•OH的生成。竞争动力学法测得k_MTZ,•OH为2.79 (±0.32) × 10⁹ M⁻¹ s⁻¹。稳态动力学模型法计算得到的k_MTZ,•OH为2.43 (±0.25) × 10⁹ M⁻¹ s⁻¹，两者高度吻合，验证了模型的可靠性。该模型同时计算出体系中[•OH]_ss为2.36 × 10⁻¹³ M。基于这些参数，贡献率分析得出：在UV/H2O2体系（[H2O2]=100 μM, pH 3.00）中，•OH氧化对MTZ降解的贡献率高达90.1%，而直接光解仅占9.9%。这确凿地证明了在UV/H2O2体系中，•OH是降解MTZ的主导活性物种，解决了“谁是主要作用者”的问题。这一结论也为后续研究H2O2和基质的影响提供了理论基点——它们主要通过影响•OH的生成与消耗来发挥作用。
第四，操作条件与水体基质影响结果如下： - H2O2投加量：实验与模型预测结果吻合良好。当H2O2浓度从50 μM增至400 μM时，MTZ的降解速率常数（k_exp,obs）从1.75 × 10⁻³ s⁻¹增至9.15 × 10⁻³ s⁻¹。模型预测显示，在此范围内，•OH的贡献率始终高于90%。然而，模型进一步预测，当H2O2浓度超过约5 mM时，降解速率会因过量H2O2淬灭•OH（H2O2 + •OH → HO₂• + H₂O）而下降。这为实际应用中确定最佳H2O2投加量提供了理论指导。 - 有机质（腐殖酸）：随着腐殖酸浓度增加，k_exp,obs从3.01 × 10⁻³ s⁻¹逐渐降至2.43 × 10⁻³ s⁻¹，表明其抑制降解。模型分析揭示，腐殖酸的抑制作用主要来自其对•OH的淬灭作用（作为•OH的清除剂），其次才是竞争紫外光吸收。实验值略高于模型预测值，暗示腐殖酸与•OH反应可能产生了一些具有氧化能力的次级自由基，对MTZ降解有微弱贡献。 - 无机阴离子：在pH 3.00条件下，NO₃⁻对MTZ降解有轻微促进作用。这是因为NO₃⁻在紫外光照下可光解产生•OH（尽管量子产率低），而其淬灭•OH的能力很弱。SO₄²⁻和Cl⁻则对MTZ降解影响不大。对于SO₄²⁻，其淬灭•OH的作用可忽略。对于Cl⁻，情况更为复杂：虽然它能与•OH快速反应生成Cl•、ClOH•⁻、Cl₂•⁻等次级氯自由基，消耗•OH，但这些次级自由基本身也具有氧化性，并能通过链式反应重新生成•OH（特别是在较高pH下），从而抵消了部分淬灭效应。实验结果表明，在本研究条件下，Cl⁻的净效应不明显。
五、 研究结论与价值
本研究系统阐明了甲硝唑在UV及UV/H2O2高级氧化体系中的降解行为与机制，并成功预测了操作参数和天然水体基质的影响。主要结论如下：1）MTZ的直接光解因其较低的摩尔吸光系数和量子产率而缓慢；2）在UV/H2O2体系中，•OH是降解MTZ的主导活性物种，贡献率超过90%；3）通过竞争动力学和稳态动力学模型两种方法，可靠地测定了MTZ与•OH的反应速率常数及•OH稳态浓度；4）H2O2投加量存在最佳范围，过低则•OH产量不足，过高则会淬灭•OH；5）天然水体基质中，有机质（以腐殖酸为代表）主要通过淬灭•OH抑制污染物降解，而无机阴离子的影响则因种类而异（NO₃⁻轻微促进，SO₄²⁻和Cl⁻影响不大）。
本研究的科学价值在于：建立了一套结合实验测定与稳态动力学模型的研究方法，能够定量解析AOPs中污染物降解的多种途径贡献，并预测复杂基质的影响。所获得的MTZ降解动力学参数和机制认识，丰富了抗生素类污染物环境转化行为的基础数据库。应用价值在于：为评估和优化UV/H2O2工艺去除自然水体（如河流、湖泊、地下水）中抗生素等微污染物的效率提供了重要的理论依据和预测工具，有助于水处理工艺的选择与设计。
六、 研究亮点
七、 其他有价值内容
本研究详细列出了UV/H2O2体系中涉及的主要反应及其速率常数（表2），包括•OH与H2O2、缓冲盐（磷酸盐）、天然有机质（NOM）、Cl⁻、碳酸盐/碳酸氢盐、NO₃⁻等的反应，这为其他研究者构建类似的动力学模型提供了宝贵的参考。此外，研究还讨论了在pH 3.00条件下未考虑CO₃²⁻/HCO₃⁻影响的原因，体现了实验设计的严谨性。论文最后指出，这种实验与模型相结合的方法可推广应用于测量和预测其他抗生素的降解动力学，为在筛选层面确定优先控制的持久性污染物提供了可行方案。