基于混合效应模型评估节能改造对室内气候条件的影响：方法论与R语言实现

《Atmosphere》：Assessing the Impact of Energy Retrofits on Indoor Climate Conditions Using Mixed Effects Models: Methodology and R Implementation

　　

能源改造干预措施在全球建筑部门追求脱碳目标和提升能效的背景下日益关键。然而，由于室外气候波动和居住者行为等混淆变量的存在，建立改造措施对室内气候条件影响的稳健因果推断仍然面临挑战。本研究提出了一种利用线性混合效应（Linear Mixed Effects, LME）模型分析改造前后室内气候数据的方法论框架，该框架能够在控制环境及行为因素的同时，显式考虑建筑层面的变异性。研究以爱尔兰住宅在节能改造干预前后水蒸气分压（Partial Pressure of Water Vapor, Pw）的案例分析进行方法演示。分析采用标准化系数（Standardized Coefficients）评估预测因子的相对重要性，并通过逐步移除非显著项实现模型简约性。研究提供了完整的R代码以供其他研究人员改编使用。结果表明，LME模型能够提供改造效应的无偏估计，同时避免困扰简单分析方法的聚合偏差（Aggregation Bias）。本研究既作为方法论参考，也为希望严格评估建筑改造在多样室内气候参数下有效性的实践人员提供了实操指南。

本研究发表于《Atmosphere》期刊，旨在解决能源改造对室内气候条件影响评估中的方法论挑战。随着全球建筑部门面临脱碳压力，欧盟及各国政府投入大量资金推进住宅节能改造，但现有研究对改造后室内环境质量和健康影响的实证证据仍然参差不齐且高度依赖情境。传统的前后对比或有限协变量控制的横截面回归方法难以有效分离改造效应与天气、建筑存量及居住者行为等背景变异性，导致估计偏差。为此，研究人员开发了一套基于线性混合效应模型的可转移、可复现分析框架，并以23户爱尔兰住宅在节能改造前后水蒸气分压的监测数据为案例进行演示。

研究背景与问题：建筑能耗占全球温室气体排放的重要份额，欧盟" Renovation Wave"计划到2030年改造3500万栋建筑。改造措施包括围护结构保温增强、门窗升级、气密性提升及化石燃料供暖系统替换等，主要目标是降低供暖能耗。然而，这些措施同时影响室内环境，进而影响居住者健康与福祉。现有证据表明，改造后热舒适度和潮湿状况可能改善，但室内空气质量结果混杂：部分研究发现高密闭性住宅中污染物浓度升高、氡水平上升或夏季过热等问题。因此，改造计划不能仅依据能源绩效评估，亟需建筑尺度的 robust 证据，说明特定改造策略在真实运行条件下如何影响室内湿度、温度和空气质量。这要求方法论能够区分改造效应与天气、建筑存量及居住者行为的背景变异性。

研究人员开展的研究与结论：研究人员构建并验证了一个完整的LME建模框架，通过固定效应控制测得的混淆因子（室外气候、季节、房间类型、居住者数量），同时利用建筑层面的随机截距和随机斜率捕捉建筑间基线条件和对外部驱动因素敏感性的差异。案例分析采用合成数据集（约28,000条小时观测记录，代表原始数据约5%的分层子样本），以水蒸气分压（Pw，单位kPa）为主要结局变量而非更常用的相对湿度（RH），原因在于Pw具有温度独立性，可避免改造后供暖模式改变导致的温度-湿度耦合混淆。

模型选择过程显示，包含改造状态、室外水蒸气分压、室外温度、室内温度、房间类型、季节和居住者数量的模型为最优结构，并包含建筑层面的随机截距以及改造状态和室外水蒸气分压的随机斜率。该模型的边际R²（固定效应解释方差）为0.75，条件R²（固定与随机效应共同解释方差）为0.82，表明模型捕获了室内湿度的大部分变异。

研究结果部分按以下结构展开：

样本特征：23户爱尔兰住宅（1994年或2000年建造，含连排、半独立式等类型），改造前监测期为2015年2–6月，改造后为2015年9月至2017年2月，涵盖供暖季与非供暖季。室内测点包括厨房、客厅和多间卧室，室外气候数据来自都柏林机场气象站。室内平均水蒸气分压改造前为0.62±0.18 kPa，改造后略高。

混合效应模型：通过两阶段模型选择，最终确定模型结构为 RoomPw ~ RetrofitStatus + ExternalPw + ExternalT + RoomT + RoomType + Season + Occupants + (1 + RetrofitStatus + ExternalPw | HomeID)。包含建筑类型、建造年份和收入变量的扩展模型因收敛问题被舍弃；移除季节和居住者的简化模型被似然比检验拒绝（χ²(2)=16.34, p<0.001）。

预测因子相对重要性：标准化系数显示，室外水蒸气分压（β?=0.63，95% CI [0.59, 0.66]）为最强预测因子，其次为室内温度（β?=0.46），室外温度呈中等负效应（β?=?0.18），房间类型（β?=0.41，其他房间高于客厅）、季节（β?=0.04）和居住者数量（β?=0.08，p=0.051）效应较小。改造效应为正且显著（β?=0.35，95% CI [0.25, 0.44]），表明控制协变量后改造后室内水蒸气分压升高。

模型规范与诊断：方差膨胀因子均<1.5，不存在多重共线性。Q-Q图显示残差近似正态，残差-拟合值图显示轻微异方差性。Durbin-Watson统计量为1.45，提示残留自相关有限。加入AR(1)相关结构的敏感性模型AIC改善（?44,763 vs. ?44,426），但为保持可解释性保留更简单的LME设定。剔除5%–95%百分位极值的敏感性分析显示，改造效应从0.090 kPa略降至0.079 kPa（降幅12.1%），但方向一致，核心推論不变。

随机效应与建筑异质性：组内相关系数（ICC）为0.29，表明29%的总方差来自建筑间系统性差异。随机截距方差（σ²=0.02）超过残差方差（σ²=0.01）。改造状态和室外水蒸气分压的随机斜率方差接近零，提示建筑对改造和室外条件的响应敏感性相似。然而，随机截距与室外水蒸气分压随机斜率的强负相关（ρ=?0.79）表明，基线湿度较高的建筑对外部 moisture 条件的耦合较弱，且该结构在改造前后保持稳定。

房间类型与季节效应：非供暖季室内湿度高于供暖季，其他房间高于客厅，改造后两种模式均整体抬升。经Tukey校正的成对比较确认，其他房间与客厅的湿度差在改造前后均为0.107 kPa（p<0.001），改造效应在两个房间类型均为0.090 kPa（p<0.001），表明改造作为统一的 moisture 基线抬升作用于所有房间。

留存时间验证：70/30分层随机抽样验证显示，测试集R²=0.817，RMSE=0.11 kPa，与全数据集条件R²=0.824高度一致，预测值与观测值紧密沿1:1线分布，无系统偏差，表明模型未过拟合且固定效应结构对新观测具有良好泛化能力。

讨论部分总结：研究强调了LME框架在控制混淆变量同时适当处理分层数据的 rigorous 性。改造后0.090 kPa的水蒸气分压升高（典型温带室内条件下约对应相对湿度4–5个百分点增幅）与既往文献中密闭性提升未伴随充分机械通风时室内 moisture 增加的发现一致，也与爱尔兰特定情境下气密性改善可能损害室内空气质量的证据相符。该幅度在已有研究中具有实际意义，因持续相对湿度超过60%与霉菌萌发及尘螨增殖条件相关，进而关联哮喘加重和过敏性鼻炎等呼吸健康结局。在爱尔兰凉爽湿润的大西洋气候中，即使适度的湿度增加也可能使住宅从生物风险阈值以下升至以上，凸显将湿度监测和充分通风规范纳入改造计划的必要性。

研究同时指出框架的跨领域适用性：任何涉及重复测量嵌套于实验单元、明确的前后干预结构、以及多混淆变量的研究设计均可采用该方法。框架优势包括：显式分区改造效应与背景条件变化，避免简单前后对比的偏差；正确量化重复测量依赖结构，防止标准误差低估和I类错误膨胀；通过随机效应量化建筑异质性，为针对性干预提供信息；保留完整数据信息，避免聚合导致的信息损失和生态学谬误；处理非平衡纵向数据，适应传感器缺失和退出等现场研究现实；以及通过开源代码提升研究透明度和可复现性。

研究局限性涉及：无法完全排除未测量混淆（如改造意识驱动的行为改变、传感器漂移或重新安置）；小时数据残留短时自相关，AR(1)结构敏感性分析虽改善拟合但未改变核心结论，更复杂时序模型可在关注时间动态时考虑；案例局限于爱尔兰住宅情境，结果向其他气候区、建筑存量或改造策略（如含热回收机械通风的深度改造）的推广需谨慎；数据集对建筑层面随机斜率差异的检测力有限。

主要技术方法：研究采用混合效应模型框架，固定效应控制室外气候（水蒸气分压、温度）、季节、房间类型和居住者数量，随机效应设定为建筑层面的随机截距以及改造状态和室外水蒸气分压的随机斜率；模型选择结合理论驱动的固定效应筛选、方差膨胀因子检验、似然比检验和AIC/BIC信息准则；采用标准化系数（sjPlot包）比较不同尺度预测因子的相对重要性；运用最大似然估计处理缺失数据，进行包含与排除极值（5th–95th percentile）的敏感性分析及AR(1)相关结构的替代模型比较；留存时间验证采用70/30分层随机抽样评估时序泛化能力。

研究结论翻译：能源改造是提升建筑部门能效和实现脱碳目标的关键干预措施。然而，严谨评估其影响需要能够处理混淆、聚类和异质性的适当统计方法。线性混合效应模型提供了这样一个框架，能够在量化建筑层面改造响应变异性的同时，实现改造效应的无偏估计。本文呈现了实施该方法的完整方法论指导和可复现R代码。建筑间改造响应的 substantial 异质性凸显了建筑特异性分析的价值。鼓励研究人员在其自身改造评估研究中采用和改编这一方法。所提供的带有详细注释的R代码大幅降低了实施的技术门槛。希望通过降低 rigorous 改造评估的技术门槛，加速建筑气候科学进展，并支持基于证据的改造政策与实践。