摘要：采用控温箱—热流计法对墙体进行检测，此过程中必会受到外界噪声的影响，同时在实际测试的过程中很难有长时间的测试数据来计算出墙体传热系数，于是采用了扩展卡尔曼滤波和粒子滤波这两种方法来预测估计，比较两种方法所得的预测值与真实测的值，仿真结果表明,粒子滤波精度和计算效率优于扩展卡尔曼滤波，同时通过计算的估计值来计算墙体传热系数。
关键字：建筑围护结构 控温箱-热流计法 粒子滤波
1引言
当今社会建筑能耗占总能耗的很大部分，因此建筑节能是现在势在必行的一项艰巨任务。建筑设计水平中的建筑热工节能设计的最重要设计参数是传热系数。建筑构件的传热系数,可在实验室标准条件下对其试件进行测试。建筑物围护结构建造完成后, 其传热系数需要现场检测才能确定。现场检测围护结构的传热系数是对新型建筑保温隔热技术及系统产品的热工性能进行检验、鉴别的必要手段, 通过检测节能建筑的实际热性能,来判定新型建筑保温隔热系统的产品、技术是否符合节能设计要求,鉴定新系统的产品、技术的优缺点等。
2建筑围护结构墙体传热系数计算方法
2.1建筑稳态传热理论
当室内外温度不等时，热量总是从温度高的一侧传向较低的一侧。在建筑热工学中，为了简化计算，墙体、屋面、地面是同一个问题—平壁稳定传热，当室内温度高于室外温度时墙体的传热由墙体内表面吸热、墙体导热、墙体外表面散热三个过程组成。连续传热过程属于稳态传热[1]。
建筑围护结构的墙体传热系数的求解就是通过控温箱—热流计法人工达到一个稳态的传热过程后得出的结果。控制热箱的加热功率，对箱内加热，墙体两面出现温差，热量开始从内墙表面传向外墙表面，连续加热箱体，最后内外墙之间的传热达到一个稳态。
2.2建筑围护结构墙体系数获取方法
控温箱-热流计检测法是一套自动控温装置，控温设备由双层框构成，具有致冷和加热功能，根据季节进行双向切换使用，夏季高温时期用致冷运行方式，春秋季用加热方式运行。根据热流计法计算公式(2.1)、(2.2)计算被测部位的热阻和传热系数[2]。
其中，K为传热系数，为热流计读数，为热流计侧头系数，为被测物的热阻，为被测物冷端表面温度，为被测物热端表面温度，为内表面换热阻，为表面换热阻。
本文采用无线式墙体传热系数现场检测装置，将热流计和温度传感器置于内墙表面及热箱内，将温度传感器置于外墙表面，在规定的测试时间内，采集数据。经过对数据进行运算与分析得出被测墙体及屋面的传热系数。建筑物围护结构传热系数及采暖供热量检测方法中规定的被测墙面的热阻的计算公式为：
3 数据处理的两种方法
3.1 扩展卡尔曼算法（EKF）
卡尔曼滤波[3]是用反馈控制的方法估计过程状态，估计某一时刻的状态，然后以（含噪声的）测量变量的方式获得反馈，利用前一时刻地估计值和现时刻的观测值来更新对状态变量的估计，求现时刻的估计值。该算法为两个部分，时间更新方程和测量更新方程。时间更新方程负责向前推算当前状态变量和误差协方差估计的值，为下一个时间状态构造先验估计。测量更新方程负责反馈，它将先验估计和新的测量变量结合以构造改进的后验估计。如果被估计的过程和观测变量与过程的关系是非线性的，提出了扩展卡尔曼滤波。将期望和方差线性化的卡尔曼滤波器称作扩展卡尔曼滤波器。
面对非线性关系时，通过求过程和量测方程的偏导来计算当前估计。将卡尔曼滤波中的公式换一种方式表示
随机变量 和代表过程激励噪声和观测噪声。差分方程式（3.1）中的非线性函数f 将上一时刻的状态映射到当前时刻k 的状态。非线性函数h 反映了状态变量 和观测变量 的关系。扩展卡尔曼也是有预测和更新两个过程。时间更新方程将状态和协方差估计从时刻向前推算到时刻。
EKF是一种常用的非线性滤波算法，在这种算法中，非线性因子的存在对滤波稳定性和状态估计精度都有影响，对于维数比较大的非线性系统，估计的过程、量测噪声协方差容易出现异常现象，只有当系统的动态模型和量测模型都接近线性时，EKF才有可能接近真实值。
3.2 粒子滤波算法(PF)
3.2.1 粒子滤波算法原理
粒子滤波是一种基于序贯蒙特卡洛模拟的近似贝叶斯滤波方法。粒子滤波是指通过寻找一组在状态空间中传播的随机样本对概率密度函数进行近似，以样本均值代替积分运算, 从而获得状态最小方差估计的过程,这些样本即称为“粒子”。
动态系统的状态模型表示为：
图4-2 两种算法热流值状态真实值、误差值的比较图
从实验数据看出热箱加热一个小时后，热箱系统已经达到一个相对稳定的系统。图4-1图分别得到粒子滤波和扩展卡尔曼滤波处理后的状态值和误差值，从误差图中可以明显看出粒子滤波的误差小。图4-2看出粒子滤波要比扩展卡尔曼滤波的效果好。同时实验中得到的均方误差粒子滤波比扩展卡尔曼小。从而得出粒子滤波得出的状态真实值更接近于真实状态。根据粒子滤波得到的箱内外温度差和热流值的状态真实值以及控温箱—热流计法原理可以计算出墙体传热系数。首先根据式(2.3)计算墙体的热阻，再根据式(2.2)计算出传热系数值。
参考文献
[1]Jones, G.F.,Jones, R.W. Steady-state heat transfer in an insulated, reinforced concrete wall: Theory, numerical simulations, and experiments, Energy and Buildings,?v 29, 293-305,?Mar 1999
[2]Yu Guoqing, Gao Gangfeng, Lin Xinfeng, Identification of building envelope overall coefficient of heat transfer using recursive least squares algorithm, Advanced Materials Research,?v 250-253,?p 3239-3244,?2011
[3]G. Welch and G. Bishop. An introduction to the Kalman Filter., Technical Report TR 95-041, University of North Carolina Department of Computer Science, 1995.
[4]Fasheng Wang?, Yuejin Lin, Improving particle filter with a new sampling strategy, Proceedings of 2009 4th International Conference on Computer Science and Education, ICCSE 2009, p 408-412, 2009
论文作者:王小琴1,杨勇2,孙光倩3
论文发表刊物:《建筑学研究前沿》2018年第3期
论文发表时间:2018/6/25
标签：系数论文; 墙体论文; 粒子论文; 热流论文; 状态论文; 建筑论文; 过程论文; 《建筑学研究前沿》2018年第3期论文;