大气中有机过氧自由基的自反应与交叉反应对臭氧形成贡献研究综述

大气中有机过氧自由基的自反应与交叉反应对臭氧形成贡献研究综述大气中有机过氧自由基（RO₂）是对流层光化学循环的关键中间物种，其自反应与交叉反应过程直接调控臭氧（O₃）和二次气溶胶的生成效率，是准确理解区域光化学污染形成机制的核心环节。RO₂主要由挥发性有机物（VOCs）的光氧化产生，在NOₓ浓度较高的城市地区，传统机制认为RO₂主要与NO反应生成NO₂，进而光解产生O₃；但在郊区、森林等低NOₓ环境中，RO₂的自反应（相同RO₂物种之间的反应）和交叉反应（不同RO₂物种之间的反应）占比显著提升，其反应路径与产物分配会改变HOₓ自由基循环效率，最终对O₃生成潜势产生不可忽视的影响。近年来随着外场观测技术和量子化学计算方法的发展，越来越多的研究发现，即使在中等NOₓ浓度的污染区域，RO₂自/交叉反应的贡献也常被低估，这也成为当前模式模拟O₃浓度偏差的重要来源之一。一、有机过氧自由基自反应与交叉反应的基本机制RO₂的自反应通式可表示为2RO₂→2RO+O₂，或2RO₂→ROH+R'CHO+O₂，其反应速率常数与RO₂的结构密切相关。小分子烷基RO₂（如甲基过氧自由基CH₃O₂）的自反应速率常数约为10⁻¹³cm³·molecule⁻¹·s⁻¹量级，而含有羰基、羟基等官能团的取代型RO₂的反应速率常数可提升1-2个数量级。RO₂的交叉反应速率通常介于两种反应物自反应速率的几何平均值附近，但当两种RO₂的反应活性差异较大时，会出现显著的非理想混合效应。比如在异戊二烯氧化体系中，异戊二烯初级过氧自由基（ISOPO₂）与羟基取代的RO₂的交叉反应速率比理论预测值高3倍以上，这类偏差目前尚未被主流化学机制完全参数化。RO₂自/交叉反应的产物分配直接决定自由基的再生效率。当反应以生成烷氧自由基（RO）为主时，RO可进一步通过分解、异构化或与O₂反应生成HO₂，维持HOₓ自由基链循环，从而促进O₃的持续生成；当反应以生成醇和羰基化合物等稳定产物为主时，会导致自由基链终止，降低O₃生成效率。最新的实验研究表明，RO₂的反应路径受温度、压力和相对湿度的共同调控：在298K以下，低温会促进稳定产物的生成，而当温度高于310K时，RO生成通道的占比可提升至60%以上；高湿环境下，部分含羟基的RO₂会与水分子形成络合物，改变反应过渡态的能量，使稳定产物通道的占比提高15%-25%，这一效应在夏季高湿度污染过程中尤为明显。二、不同大气环境下RO₂自/交叉反应的贡献占比在森林等生物源VOCs主导的低NOₓ区域（NO浓度<50pptv），RO₂的自/交叉反应是其主要汇，占RO₂总消耗的比例可达40%-80%。亚马逊雨林的外场观测显示，异戊二烯排放产生的RO₂中，约65%通过自/交叉反应消耗，对应的O₃生成潜势比仅考虑RO₂+NO反应的模拟结果低30%左右，这解释了部分热带雨林地区虽然VOCs排放强度高，但O₃浓度维持在较低水平的现象。我国西双版纳热带森林的观测结果也显示，干季NOₓ浓度较低时，RO₂自/交叉反应对RO₂总消耗的贡献可达55%，使该区域O₃生成效率从2.5ppbv/ppbv（VOCs反应量）降至1.8ppbv/ppbv。在城乡结合部等中等NOₓ区域（NO浓度50-500pptv），RO₂自/交叉反应与RO₂+NO反应的贡献相当，两类反应的竞争关系显著影响O₃的生成敏感性。2022年夏季在北京郊区开展的综合观测发现，上午光照增强阶段，RO₂浓度快速上升至1.5×10⁸molecule·cm⁻³，此时自/交叉反应对RO₂消耗的贡献约为35%，若忽略该过程，模式会高估O₃生成速率约22%。长三角地区的模拟研究也表明，在中等NOₓ污染条件下，RO₂自/交叉反应会使O₃生成对VOCs的敏感性降低，对NOₓ的敏感性增强，这意味着该区域的O₃控制策略需要更精准的NOₓ减排力度，若仅依靠VOCs减排可能达不到预期效果。在城市等NOₓ高值区域（NO浓度>500pptv），传统认为RO₂+NO反应占绝对主导，但近年的观测发现，在早晚高峰VOCs高强度排放时段，RO₂浓度峰值可达3×10⁸molecule·cm⁻³，此时自/交叉反应的贡献仍可达10%-15%。2023年广州市中心城区的观测显示，早高峰时段汽油车排放的烯烃类VOCs快速氧化产生大量短寿命RO₂，其交叉反应速率可达10⁻¹¹cm³·molecule⁻¹·s⁻¹，贡献了约12%的RO₂消耗，对应的O₃生成量比传统模拟结果低8%左右。这类偏差会导致城市O₃峰值浓度的模拟误差，尤其在连续污染过程中，误差会随着反应时间累积，最终使3天平均O₃浓度的模拟偏差超过15%。三、RO₂结构对反应活性及O₃贡献的调控机制不同结构的RO₂的自/交叉反应活性存在数量级差异，是模式模拟不确定性的重要来源。烷基RO₂的反应活性随碳链长度增加而提升，C1-C4烷基RO₂的自反应速率常数从1.2×10⁻¹³cm³·molecule⁻¹·s⁻¹上升至8.5×10⁻¹²cm³·molecule⁻¹·s⁻¹，碳链长度超过5后，反应速率趋于稳定。含不饱和键的烯基RO₂的反应活性是同碳数烷基RO₂的2-3倍，因为双键的共轭效应可稳定反应过渡态，降低反应能垒。芳香族RO₂的自反应速率常数普遍在10⁻¹¹cm³·molecule⁻¹·s⁻¹量级，其交叉反应与小分子RO₂的反应速率更高，是城市地区RO₂反应体系中不可忽视的组分。含氧官能团对RO₂反应活性的调控作用更为显著。α位含羟基的RO₂的自反应速率比无取代的烷基RO₂高1-2个数量级，因为羟基可形成分子内氢键，稳定反应生成的四元环过渡态。β位含羰基的RO₂会发生分子内异构化，生成具有更高反应活性的双官能团RO₂，其交叉反应速率可达10⁻¹⁰cm³·molecule⁻¹·s⁻¹量级。最新的量子化学计算表明，含有过氧羟基（-OOH）的多官能团RO₂的自反应会发生分子内氢转移，直接生成H₂O₂和稳定的二羰基化合物，该路径的占比可达40%以上，会导致自由基的净消耗，显著降低O₃生成效率。四、RO₂自/交叉反应对O₃生成贡献的量化方法目前量化RO₂自/交叉反应对O₃贡献的方法主要包括观测约束法、模式敏感性分析和同位素示踪法三类。观测约束法基于实测的RO₂、NO、HO₂等物种浓度，直接计算各反应的速率，进而估算其对O₃生成的影响。该方法的优势是直接反映实际大气状况，但受RO₂准确观测的限制，目前仅能在少数观测站点应用。近年来发展的化学电离质谱（CIMS）技术可实现多种RO₂物种的同步测量，时间分辨率可达1秒，为量化复杂体系下RO₂各反应路径的贡献提供了可能。模式敏感性分析通过在化学传输模式中开启或关闭RO₂自/交叉反应模块，对比O₃浓度的差异来量化其贡献。该方法可反映区域尺度的贡献特征，但受模式中反应参数和排放清单的不确定性影响较大。目前主流的CB05、MCMv3.3、RADM2等化学机制中，RO₂自/交叉反应的参数化方案存在显著差异，其中MCMv3.3详细考虑了不同结构RO₂的反应速率和产物分配，其模拟结果比CB05高10%-20%。最新的研究通过将实验测定的200多种RO₂的反应参数纳入模式，使区域O₃模拟的平均偏差从22%降至13%，显著提升了模拟性能。同位素示踪法利用¹⁴C或¹⁸O标记的前体物，追踪RO₂反应生成的O₃中的同位素组成，进而区分不同RO₂反应路径的贡献。该方法可在复杂的大气本底中准确分离生物源和人为源RO₂的贡献，尤其适用于城乡交错区域的研究。2024年发表的一项研究利用¹⁸O标记的异戊二烯，量化得出我国华北平原夏季生物源RO₂的自/交叉反应贡献了约18%的O₃生成量，比之前模式估算的结果高7%。五、当前研究的不足与未来方向尽管近年对RO₂自/交叉反应的研究取得了显著进展，但仍存在三个关键科学问题有待解决：一是复杂大气环境下RO₂的反应机制尚未完全明确，尤其是多官能团RO₂的交叉反应路径和产物分配的实验数据十分缺乏，现有参数化方案的不确定性超过50%；二是实际大气中RO₂的组成和浓度观测不足，当前观测技术仅能识别30%左右的RO₂物种，大量高摩尔质量的多官能团RO₂尚未被准确测量，导致其对O₃生成的贡献被低估；三是区域模式中RO₂反应的参数化过于简化，多数机制仅将RO₂分为有限的几类，无法反映不同结构RO₂的反应活性差异，是O₃模拟偏差的重要来源。未来研究需要重点关注三个方向：首先是发展高时间分辨率、高物种分辨度的RO₂观测技术，实现对C1-C10范围内不同官能团RO₂的同步定量测量，明确实际大气中RO₂的组成和时空分布特征；其次是开展接近实际大气条件的模拟实验，测定不同温度、湿度、压力下RO₂自/交叉反应的速率常数和产物分配，建立更精准的反应参数数据库；最后是改进化学传输模式中的RO₂机制，将结构化的反应参数纳入模式，提升O₃模拟的准确性，为区域光化学污染防控提供科学支撑。RO₂自/交叉反应对O₃形成的贡献是一个多尺度、多因素耦合的复杂问题，需要结合外场观测、实验室实验