欧洲十三种室内微环境中颗粒物致癌风险评估：便携式空气净化器干预前后评估的局限性（TwinAIR）

《Toxics》：Particulate-Matter Carcinogenic Risk in Thirteen European Indoor Micro-Environments: The Limits of Pre–Post Evaluation of Portable Air Cleaning (TwinAIR)

　　

室内颗粒物被列为第一类致癌物，且未确定暴露阈值。目前广泛部署于欧洲学校、医院和办公室的便携式空气净化器，其效果几乎总是通过比较同一房间安装前后的平均浓度变化来评估。本研究应用概率增量终生癌症风险（ILCR）框架，对五个国家十三个有人居住的微环境中的粒径分级颗粒物进行了评估，这些数据来自TwinAIR试点网络，以五分钟分辨率连续监测。随后，研究人员利用所依据的测量数据检验了上述评估设计。结果显示风险普遍较低：所有场所和粒径分级的估计值均低于10-6的最低关注阈值（de minimis threshold），低至一个到三个数量级不等。这一结论在基于实测浓度分布重新计算后依然成立（四十二项估计值低于阈值11.2至353倍），并且在使用儿童暴露参数替代后也保持不变。然而，评估设计本身未能经受住检验。在配备空气净化器的两个房间中，实测PM2.5分别上升了7.8%和21.0%，而解释方差仅1%-6%的热湿回归模型——即概率估计所依据的浓度模型——却预测分别下降6.7%和2.3%。将相同的阶段边界应用于三个未安装设备的微环境中的四个传感器，产生的变化范围为-74.3%至+37.2%，其在幅度上超过了归因于干预措施的对比效应，并在方向上与之相反。因此，本研究的创新性体现在两个方面：一是首次在同一尺度上对十三个真实有人居住的室内空间进行了统一的概率风险评估；二是利用同一网络的数据证明，在没有同期对照的情况下，室内干预前后的设计无法识别空气净化的效果。

论文解读文章

一、研究背景与问题提出

环境颗粒物（PM）已被国际癌症研究机构（IARC）列为第一类人类致癌物，其暴露-反应函数在低至约2 μg/m3的水平下仍呈线性且无确定阈值，意味着不存在无风险的暴露水平。然而，绝大多数证据来源于室外浓度研究，而欧洲人绝大部分时间生活在室内。室内颗粒物浓度受室外渗透和室内源共同控制，与室外监测仅有弱相关性，且不同建筑类型、季节和通风模式下差异可达一个数量级以上。世界卫生组织（WHO）修订后的空气质量指南已将年均PM_2.5指导值降至5 μg/m3，许多有人居住的欧洲室内空间在工作日的大部分时间超过此限值。政策应对措施主要是加强通风，以及自2020年以来大规模部署便携式高效空气净化器。尽管现有证据表明过滤可降低室内颗粒物浓度，但仍存在三个关键空白：第一，干预研究几乎只报告浓度变化而非风险变化，而风险取决于暴露时长、频率和吸入率等参数，并非简单的比例缩放；第二，干预研究仅报告集中趋势，对于一种线性无阈值的致癌物而言，暴露分布的上尾承载了不成比例的累积剂量，仅报告均值可能遗漏重要信息；第三，现有的室内风险评估多为单一地点、单一国家和单一建筑类型，缺乏跨场地可比性。为此，TwinAIR项目设计了统一的监测网络和评估方案，旨在解决上述问题。

二、研究方法概述

研究人员采用了一种结合概率模拟与多变量统计的分析框架。首先，在每个微环境中部署统一的多传感器节点，以五分钟分辨率连续记录温度、相对湿度、CO₂、总挥发性有机物（tVOC）以及PM₁、PM_2.5、PM₄和PM₁₀的质量浓度。暴露建模方面，研究人员并未直接将实测PM浓度作为固定点估计，而是将其建模为温度和相对湿度的线性函数（普通最小二乘回归），以再现真实空间中环境条件与颗粒物水平之间的协方差。随后，依据美国环保署（USEPA）吸入风险评估框架计算慢性每日摄入量（CDI）和增量终生癌症风险（ILCR），其中斜率因子（SF）按粒径分级设定：PM_2.5为0.005 (mg/kg/day)^-1，PM₁₀为0.003，PM₁为0.002。所有风险表达式均通过蒙特卡罗模拟（Monte Carlo simulation）进行概率化评估，每次迭代生成10,000次联合实现。此外，还采用了秩相关敏感性分析和多变量线性回归模型来识别驱动输出变异的关键输入变量。为验证模型可靠性，研究人员直接利用实测的五分钟浓度分布重复了整个暴露评估，并对未安装空气净化器的微环境进行了阴性对照分析。

三、研究结果

基线致癌风险：十三个微环境的基线（PRE）模型化增量终生癌症风险跨越约两个数量级，从瑞典TP4005实验室PM₁的6.63×10^-10到希腊小学教职员室PM_2.5的4.63×10^-8。所有值均低于10^-6的最低关注阈值，低至一到三个数量级，且没有任何场所接近10^-4的可接受上限。粒径排序一致：在所有同时解析了三个粒径分级的微环境中，平均ILCR遵循PM_2.5 > PM₁₀ > PM₁的顺序。

阶段边界前后的模型化平均风险：六个希腊微环境获得了配对的前后（PRE/POST）模型化估计。在安装了空气净化器的两个医院房间中，模型化平均风险在阶段边界后降低了7.3%-14.8%。然而，未安装设备的希腊小学教室同样出现了6.6%-9.6%的降低，而教职员室的降幅最大达45.0%-49.9%。由于这两个房间均未接受干预，这些差异不能归因于空气净化，而与希腊站点供暖季结束的时间重合。

模型化风险的离散度：离散度的变化比均值更为显著。在医院医生办公室，模型化风险的标准差降低了29.9%-37.3%，而均值差异仅为8.0%-14.8%。然而，这种变窄现象也在未安装空气净化器的房间中出现，例如小学教室PM₁₀和PM_2.5的标准差降低了约41%。

PM₁的例外情况：PM₁的行为与较粗粒径分数不同。在小学教室中，PM₁的标准差在阶段边界后反而升高了22.9%，变异系数从0.347升至0.456。由于该模式同时出现在有和无空气净化器的房间，它不可能是过滤的特性，更可能与亚微米颗粒物的产生机制及其对热湿变化的敏感性有关。

气态污染物：总挥发性有机物（tVOC）的变化方向在不同场所间不一致，且与颗粒物结果无关。在Triklala小学教室，平均tVOC下降了近90%；而在同一学校的教职员室，平均值反而上升，最大值增长了五倍以上。CO₂数据仅来自无空气净化器的房间，Triklala教职员室最高超过1400 ppm，都柏林教室最高超过5600 ppm，表明这些地点的通风供应不足。

模型化风险的决定因素：在前干预时期，模型化致癌风险由房间的热湿状态驱动，而非实测PM浓度。在医院医生办公室的回归模型中，没有一个实测颗粒物通道达到显著性，而温度和相对湿度的多个传感器通道却呈现显著关联。值得注意的是，同一节点内四个温度通道对同一结局变量的回归系数方向相反，表明传感元件在室内的空间和热微观位置实质上决定了其检测到的关联。

实测浓度与阴性对照分析：由于ILCR与浓度成线性关系，研究人员可直接利用经验分布重复评估。在安装空气净化器的两个房间中，实测PM_2.5并未下降反而上升：医生办公室增加7.8%，急诊室增加21.0%。模型化估计之所以报告7%-15%的降低，是因为使用了仅解释1%-6%方差的弱预测回归模型替代了实测浓度。将相同的阶段边界应用于未安装任何设备的微环境，实测PM_2.5的变化范围从-74.3%到+37.2%，其幅度超过了干预房间观察到的变化。综合来看，这三项分析表明，本研究采用的干预前后对比无法识别便携式空气净化的效果。

四、讨论与结论

讨论要点：研究人员的三个主要发现中，最持久的成果是方法论层面的。第一，虽然所有场所的颗粒物致癌风险绝对值均较低，但干预前后设计在此设置中无法识别空气净化效果。第二，离散度比均值更重要：对于线性无阈值致癌物，暴露分布的上尾承载着剂量，仅报告均值会丢失关键信息。第三，环境背景而非浓度本身驱动模型化风险，且传感器来源问题值得重视——同一节点内不同温度通道的回归系数方向相反，表明数字孪生风险引擎应传播传感器层级的不确定性而非通过平均来解决。此外，tVOC与PM相关风险之间的近乎零相关证实，颗粒物干预不应期望管理化学暴露。CO₂最高值超过1400甚至5600 ppm表明多个房间通风不足，过滤只能去除颗粒物而不能引入新风，是通风的补充而非替代。

研究结论：在五个国家的十三个有人居住的微环境中，室内颗粒物的模型化增量终生癌症风险均低于10^-6的最低关注阈值，这一结论对暴露度量选择具有稳健性：根据实测的五分钟浓度分布重新计算后，每项估计仍低于阈值11至353倍。在负荷最高的教室中，采用六年暴露期的儿童场景下，根据实测浓度计算的裕度缩小至1.4倍。研究人员无法证明便携式空气净化对该风险存在效果。安装净化器的两个房间中实测浓度在研究窗口内反而上升；概率模型报告的下降低是由于用弱预测的热湿回归替代了实测浓度；将相同阶段边界应用于未接受任何干预的房间产生了幅度相等或更大、方向相反的变化。本研究的持久成果是方法论层面的：在没有同期对照、且暴露由低R²环境回归估计的情况下，室内干预前后设计无法识别空气净化的效果。对此类干预的评估应报告实测分布、同期对照以及暴露的离散度。研究提出三项建议：室内空气干预研究应采用分布感知终点——变异系数和高百分位数——与均值并列作为共同主要终点；便携式过滤应被视为通风的补充而非替代；空气净化效果的评估应使用同期对照房间，并记录洁净空气输送率、房间体积和围护结构完整性。