确定平屋面节能模拟与规范中气候特异性有效PIR导热系数及厚度取值的方法论——以平屋面为例

《Journal of Building Engineering》：Methodology for Determining Climate-Specific Effective PIR Thermal Conductivity and Thickness Values for Energy Modeling and Code Prescription: A Case of Flat Roofs

　　

建筑物贡献了加拿大12%的温室气体(Greenhouse Gas, GHG)排放，其中商业及公共机构建筑约占建筑排放的52%，使得屋面成为提升能效的重要目标。准确模拟屋面的传热性能（尤其是采用聚异氰脲酸酯(Polyisocyanurate, PIR)保温层的屋面，其导热系数随温度呈非线性变化）对评估建筑节能策略至关重要。本文建立瞬态传热模型并以现场实验数据验证。对于所研究的供暖主导气候，考虑PIR低温段导热系数升高后，按温度依赖性导热系数计算的比按规范规定的24°C恒定导热系数计算的供热负荷更大。为弥合两种算法差异，提出确定有效平均保温工作温度(Mean Insulation Operating Temperature, MIOTeff)的方法，在该温度下恒定有效导热系数(keff)可等效代表温度依赖型保温性能。演示该方法在五个加拿大城市的应用，最终给出用于能耗模拟与规范制定的气候特异性有效PIR导热系数及厚度取值图表。

论文解读：确定平屋面PIR保温气候特异性有效导热系数与厚度取值的方法及应用

该研究由Fitsum Tariku、Yina Shang与Sudhakar Molleti（Simon Fraser University可持续能源工程学院）完成，发表于《Journal of Building Engineering》。建筑部门占加拿大温室气体(Greenhouse Gas, GHG)排放12%，其中商业及公共机构建筑占建筑排放约52%，平屋面作为低层商业建筑围护结构重要组成部分，其保温(Insulation)性能直接影响建筑能耗。现行建筑能源规范（如加拿大国家建筑能源规范National Energy Code of Canada for Buildings, NECB）以标准平均温度24°C下测得的单一R值（热阻）标定保温材料热工性能，但聚异氰脲酸酯(Polyisocyanurate, PIR)保温的导热系数k具有温度依赖性——呈"U形"非线性关系，最低k出现在约15°C，低于10°C时因发泡剂冷凝导致k升高，高于最优温区k也随温升增大。传统以k₂₄（24°C下导热系数）进行全年能耗计算会高估寒冷地区PIR屋面的实际保温效果，造成设计—实际性能偏差。为此，研究人员通过卡尔加里现场足尺实验结合瞬态RC网络(Resistance-Capacitance Network, RC-Network)模型，提出有效平均保温工作温度(Mean Insulation Operating Temperature, MIOT_eff)概念及确定方法，给出不同加拿大气候区下匹配温度依赖性性能的等效恒定有效导热系数(k_eff)及若沿用k₂₄时应增补的保温厚度，弥补规范与实践差距。

主要关键技术方法：

研究人员在加拿大卡尔加里某商业建筑设置两组平屋面足尺测试段——R1（黑色EPDM(乙烯丙烯二烯单体,Ethylene Propylene Diene Monomer)膜＋63.5 mm PIR）与R2（白色TPO(热塑性聚烯烃,Thermoplastic Polyolefin)膜＋139.7 mm PIR），埋设T型热电偶与热流计全年采集温度与热流；NRC(加拿大国家研究委员会,National Research Council)热流仪测定PIR样品在2°C、12°C、24°C、42°C下的k?T（温度依赖导热系数）并拟合三阶多项式；开发双层可变电阻—单电容的RC?Network瞬态传热模型，以CWEC(加拿大能源计算典型气象年,Canadian Weather Year for Energy Calculation)数据驱动，通过与有限元(Finite Element Model, FEM)及实测数据对比验证(RMSE温度预测0.62°C、热流0.22 W/m2)；以全年变温k_T计算真实能荷Q^v，遍历恒定保温温度T_c求取使∣Q^v－Q^c_L∣最小的MIOT_eff及对应k_eff，并反算若用k₂₄达到规范预期能荷所需增加厚度dL，应用于温哥华、渥太华、卡尔加里、温尼伯、怀特霍斯五城市NECB分区。

研究结果：

2. Experimental Study（实验研究）

R1与R2全年MIOT（平均保温工作温度）分别介于-4.6~51.3°C与-14.1~34.1°C，高于24°C的时间仅占17.5%（R1）和5.5%（R2）。PIR的k随MIOT呈U形变化，最低k位于MIOT≈15°C。黑膜(R1)夏季有向内得热，白膜(R2)全年向外散热；厚保温使R2热流波动更小。实测证实实际MIOT大幅偏离标准24°C，且k₂₄全年仅约4%时段出现，直接用k₂₄不能反映真实传热。

3.1. Numerical Model and Validation（数值模型与验证）

将PIR层拆为上、下两半分别赋予温度依赖热阻R_t、R_b及中心节点热容C建立RC?Network常微分方程(Ordinary Differential Equation, ODE)，逐时求解MIOT与热流。模型与FEM及R2全年实测吻合良好（温度RMSE 0.62°C vs FEM 0.58°C；热流0.22 vs 0.20 W/m2），证明简化RC模型可可靠用于后续气候区分析。

3.2. Methodology to Determine the Effective PIR Thermal Conductivity Values（确定PIR有效导热系数取值方法）

以变温k_T算全年累计能荷Q^v，对一系列假定恒定保温温度T_c算恒定k下的Q^c_L，取ΔQ＝Q^v－Q^c_L＝0对应的T_c定义为MIOT_eff，该温度下k(MIOT_eff)即为等效恒定有效导热系数k_eff。因PIR的k?T曲线U形对称，同一k_eff可能对应低温侧与高温侧两个MIOT_eff（如卡尔加里226 mm厚时MIOT_eff为3°C或30°C，k_eff＝0.0287 W/(m·K)）。绘制不同保温厚度L下MIOT_eff–L关系曲线及年能荷–L曲线，可从图读取：①用k_eff时应取哪个MIOT查厂商标定k；②若坚持用厂商标定k₂₄，需将规范厚度增至使实际能荷等于k₂₄设计能荷的新厚度L＋dL。

3.3. Determination of Effective PIR Thermal Conductivity Values for Different Insulation Thicknesses and Climates Zones（不同保温厚度与气候区的PIR有效导热系数确定）

对NECB分区城市取规范隐含PIR厚度（温哥华174 mm、渥太华198 mm、卡尔加里226 mm、温尼伯226 mm、怀特霍斯249 mm），变温模型算出实际年总能荷较k₂₄设计值偏高0.1%（温哥华）至8.9%（怀特霍斯）。为等价变温性能，MIOT_eff分别为：温哥华6°C（或24°C）、渥太华3°C（或30°C）、卡尔加里3°C（或30°C）、温尼伯0°C（或35°C）、怀特霍斯-1°C（或38°C），对应k_eff略高于k₂₄。若仍用k₂₄，需增厚度dL分别为0 mm、5 mm、6 mm、13 mm、20 mm。给出152–279 mm各厚度下五城市MIOT_eff与k_eff对照表及能荷偏差图。

讨论与结论总结（译自Conclusion）：

在寒冷供暖主导气候下，PIR低温段k升高致使按温度依赖导热系数算得的年供热负荷高于按24°C恒定k₂₄计算结果，偏差最高达8.9%（怀特霍斯），表明冷区须采用温度修正导热值。通过最小化变温与恒温算法能荷差确定的MIOT_eff及k_eff可使恒导热模型等价反映实际PIR屋面热工表现；若沿用厂商标定k₂₄，需按气候区适当增加保温厚度（温哥华0 mm，渥太华约5 mm，卡尔加里约6 mm，温尼伯约13 mm，怀特霍斯约20 mm）。文末提供气候特异性MIOT_eff–保温厚度图及k_eff–厚度表，供规范制定者与设计师确定匹配实际性能的PIR保温取值或补偿厚度。建议后续结合厂商提供的宽温区（尤低温段）k?T试验数据做同类分析。