一种用于供热系统质调节的加权外温的确定方法与流程

1.本发明涉及一种用于供热系统质调节的加权外温的确定方法，属于能源利用的工程领域。


背景技术：

2.建筑热负荷随室外温度动态变化。对此，为保证热源端热量供给与建筑用户端热量需求的相互平衡，避免对用户端“过量供热”而导致的不必要能源消耗，需对供热系统的热量供给采取调节措施，以实现按需供热。运行调节是确保供热系统安全、稳定和经济运行的关键环节。质调节是最常用的供热系统调节方式，是指保持供热系统的供水流量不变，而通过改变供水温度来保证供热需求的调节方式。
3.气候补偿技术是在我国应用最广泛的供热质调节方法，20世纪90年代就已经开始应用。气候补偿器是调节热源出力与用户端热量需求相互平衡的设备，该设备的基本工作原理是根据室外温度调节热源出力，从而将系统二次管网出口的供水温度控制在合理的范围内，以满足用户端的热量需求，实现热量的供需平衡。然而，实际应用中气候补偿器并不能恰好满足用户端的热量需求，过量供热现象依然存在。究其原因，主要在于气候补偿器是基于稳态传热的原理，对室外温度变化所产生的影响进行线性补偿。然而，由于建筑围护结构热惯性的作用，室外温度对室温的影响表现出衰减和延迟特性，也即室温随外温并非呈现出线性变化。缺乏对外温影响衰减和延迟特性的合理反映，会导致对室内温度的估计存在偏差，相应无法对外温变化做到合理的补偿。
4.因此，对供热系统供水温度的质调节，不应只对如外温当前时刻的变化进行补偿，还需考虑围护结构热惯性作用下外温影响的衰减和延迟特性，同时补偿当前时刻及之前历史时刻外温对当前时刻室温的影响，由此避免对室内温度不合理的估计，进而保证热源端热量供给与用户端热量需求的相一致，即实现热量的“供需平衡”。


技术实现要素：

5.本发明的目的在于：考虑建筑围护结构热惯性作用下外温影响的衰减和延迟特性，构建房间室温与当前时刻及之前历史时刻外温相互关系的数学模型，再应用建筑热环境模拟软件dest计算获得房间全年的逐时室温，进一步与逐时外温构成统计样本数据组，进而应用多元线性回归的统计分析方法计算确定当前及之前历史时刻外温对应的加权系数，由此构建一套加权外温的分析确定方法。该方法以外温的加权形式反映出当前及之前历史时刻外温对当前室内温度环境的共同影响，体现出建筑围护结构热惯性作用下外温影响的衰减和延迟特性。为达到发明目的，本发明采用一种用于供热系统质调节的加权外温的确定方法，所述方法包含以下步骤：
6.(1)房间室温的模拟计算
7.针对某建筑研究对象，考虑只有外温作用的情况，应用建筑热环境模拟软件dest进行外温作用下房间室温全年的逐时模拟计算，并与逐时外温构成统计样本数据组，
8.(2)室温随外温的变化分析
9.根据《建筑环境系统模拟分析方法
‑
dest》中公式(2
‑
30)和(2
‑
35)，得到房间逐时室温与逐时外温的相互关系如下：
10.δt(τ)＝∑
i
δt
i
(τ)
ꢀꢀꢀꢀꢀ
(1.1)
[0011][0012]
式中，i为房间热特性和热反应系数的序号；λ
i
为房间热特性系数；是外温作用于房间的热反应系数；u
t
为外温，℃；δt为房间室温，℃；δt
i
为对应于各个特征系数的房间室温的分解分量，℃；τ为时间，h；δτ为时间间隔，在此为1h；
[0013][0014]
公式(1.2)表示出当前时刻室温分量与前一时刻室温分量、当前时刻和前一时刻外温的相互关系，以此类推，可得到前一时刻室温分量与之前两个时刻室温分量、前一时刻和之前两个时刻外温的相互关系，之前两个时刻室温分量与之前三个时刻室温分量、之前两个时刻和之前三个时刻外温的相互关系，
……
，以及之前n
‑
1时刻的室温分量与之前n时刻室温分量、之前n
‑
1时刻和之前n时刻外温的相互关系，具体如公式(2.1)～(2.n
‑
1)所示，
[0015][0016]
将公式(2.1)～(2.n
‑
1)中各个时刻的室温分量依次顺序代入公式(1.2)，并进行合并同类项，得公式(3)，
[0017][0018]
房间热特性系数λ
i
为小于0的常数，这使得n数量足够大时，将趋于0，当小于0.01时，认为去除公式(3)右侧的项仍可满足计算精度要求，这表示历史时刻距离当前时刻足够长时，可以忽略外温的历史状态对当前时刻室温的影响，进而公式(3)简化为房间室温的当前分量δt
i
(τ)与各个时刻的外温u
t
(τ)，u
t
(τ
‑
δτ)，u
t
(τ
‑
2δτ)，
……
u
t
(τ
‑
(n
‑
1)δτ)的线性关系，
[0019]
再依据公式(1.1)，将与各个序号i对应的房间室温分量δt
i
求和，最终得到当前时刻房间室温δt(τ)与各个时刻外温u
t
(τ)，u
t
(τ
‑
δτ)，u
t
(τ
‑
2δτ)，
……
u
t
(τ
‑
(n
‑
1)δτ)的线性关系，如公式(4)所示，
[0020]
δt(τ)＝c0u
t
(τ) c1u
t
(τ
‑
δτ) c2u
t
(τ
‑
2δτ)
……
c
nn
‑‑1u
t
(τ
‑
(n
‑
1)δτ)
ꢀꢀꢀꢀ
(4)
[0021]
式中，c0，c1，c2，
……
c
n
‑1定义为各个时刻外温对室温的贡献率，
[0022]
采用建筑热环境模拟软件dest，计算获得全年逐时室温及相对应的外温，进而根据公式(4)构造出室温与多个外温相互关系的多元线性方程组，对此，采用多元线性回归的统计分析方法，计算获得公式(4)中各个时刻外温对室温的贡献率c0，c1，c2，
……
c
n
‑1。
[0023]
(3)加权外温的提出及构造
[0024]
公式(4)表明，当前时刻室温是当前及之前历史时刻外温共同影响的结果，这体现出建筑围护结构热惯性作用下，外温对室内温度环境影响的衰减和延迟特性，进一步，假设一个稳定的室外空气温度t
w
，则应用公式(4)可得该稳定外温t
w
作用下的房间室温为：δr(τ)＝c0t
w
(τ) c1t
w
(τ
‑
δτ) c2t
w
(τ
‑
2δτ)
……
c
n
‑1t
w
(τ
‑
(n
‑
1)δτ)
[0025]
＝ct
w
(τ)
ꢀꢀꢀ
(5)
[0026]
式中，
[0027]
若这个稳定的室外温度对室温的影响与实际动态变化的室外温度对室温的影响相同，可将这个稳定的室外温度t
w
定义为实际动态变化外温的加权外温，进一步结合公式(4)和公式(5)，转换得到加权外温t
w
的计算式如下：
[0028][0029][0030]
式中，d0，d1，d2，
……
d
n
‑1为与各个时刻外温相对应的加权系数，
[0031]
可以看出，加权外温体现出外温对室内温度影响的衰减和延迟特性，进而真实反映出当前及之前历史时刻外温对当前时刻室温的共同影响，由此可在供热系统质调节的分析过程中将繁琐的计算过程加以简化，并同时具有明确的物理意义。
[0032]
(4)加权外温的分析确定
[0033]
采用建筑热环境模拟软件dest计算得到外温作用下的房间逐时室温，进而与逐时外温构成统计样本数据组，在此基础上，采用多元线性回归的统计分析方法，以外温对室温贡献率小于5％的前n
‑
1时刻为起始点，将该时刻之前历史时刻的外温去除，进而将当前时刻室温与当前时刻及其之前n
‑
1时刻的共n个外温数据建立线性模型，进一步重新进行外温对室温贡献率的分配计算，直至去除所有贡献率不明显的各个之前历史时刻的外温，最后，应用逐时室温和外温所构成的统计样本数据组，采用多元线性回归的统计分析方法，计算得到各个时刻外温对室温的贡献率c
i
，进而按照公式(7)计算确定与各个时刻外温相对应的加权系数d
i
，由此构造获得外温的加权形式。
[0034]
本发明所述的基于建筑热环境模拟计算和多元线性回归统计分析的加权外温确定方法的有益效果主要体现如下：
[0035]
该方法以加权外温的形式，体现出建筑围护结构热惯性作用下外温对房间室温影响的衰减和延迟特性，真实反映出当前时刻及之前历史时刻外温对当前时刻室温的共同影
响和作用，除室外温度外，太阳辐射得热、室内人员、灯光和设备等热扰同时影响室内温度环境以及维持室内温度恒定所需要的建筑热负荷，对这些扰量可采用与外温相一致的分析方法构造建立相对应的加权扰量，基于此可在供热系统质调节的过程中，合理调节供水温度，以对外温等扰量的影响做出相匹配的补偿，同时本专利所述的方法可操作性强，应用流程和步骤清晰，保证了本发明的可实施性。
附图说明
[0036]
图1为本发明方法的流程示意图。
[0037]
图2为房间室温和外温随时间的逐时变化关系图。
具体实施方式
[0038]
以北京地区某建筑的项层北向房间为研究对象，已知该建筑的平立面布局和建筑围护结构构造，房间室内外无通风换气，室内无人员、灯光和设备等热扰。对此研究对象，结合附图1，说明本发明的具体应用思路和应用步骤如下：
[0039]
1.房间室温的模拟计算
[0040]
在建筑热环境模拟软件dest中进一步将房间太阳辐射得热设为0，进而计算获得项层北向房间的逐时室温，并与逐时外温构成统计样本数据组。对于本案例，由于不存在通风换气、室内人员、灯光和设备等热扰以及太阳辐射得热对室内温度环境的影响和作用，外温成为决定房间室温的唯一扰量。房间室温和外温随时间的逐时变化关系如图2所示。
[0041]
2.室温随外温变化的多元线性回归
[0042]
取n初始值为10，建立当前时刻室温与当前时刻至之前9个时刻外温之间的多元线性回归模型，采用ibm spss statistics 25软件进行拟合计算，得到当前时刻以及之前9个时刻的室外温度对当前时刻室温的贡献率，拟合结果如式(8)：
[0043]
δt(τ)＝0.19u
t
(τ) 0.24u
t
(τ
‑
δτ) 0.15u
t
(τ
‑
2δτ) 0.13u
t
(τ
‑
3δτ)
[0044]
0.09u
t
(τ
‑
4δτ) 0.07u
t
(τ
‑
5δτ) 0.05u
t
(τ
‑
6δτ)
ꢀꢀꢀꢀ
(8)
[0045]
0.04u
t
(τ
‑
7δτ) 0.03u
t
(τ
‑
8δτ) 0.02u
t
(τ
‑
9δτ)
[0046]
3.加权系数计算及加权外温的构造
[0047]
根据公式(8)的计算结果，在室温随外温变化的线性模型中，保留贡献率大于等于5％的当前以及之前6个时刻的室外温度，并重新拟合计算这7个时刻的室外温度对当前时刻室温的贡献率，分别为0.16、0.20、0.18、0.14、0.11、0.08和0.12。在此基础上，根据公式(7)计算得出当前时刻以及之前6个时刻外温的加权系数值分别为0.17、0.20、0.18、0.14、0.11、0.08和0.12。由此，构造得到加权外温的表达式，如公式(9)。t
w
(τ)＝0.17u
t
(τ) 0.20u
t
(τ
‑
δτ) 0.18u
t
(τ
‑
2δτ) 0.14u
t
(τ
‑
3δτ)
[0048]
0.11u
t
(τ
‑
4δτ) 0.08u
t
(τ
‑
5δτ) 0.12u
t
(τ
‑
6δτ)
ꢀꢀꢀꢀꢀꢀꢀ
(9)。