CN113822569B - 一种空调系统负荷计算方法、装置、设备及介质 - Google Patents
一种空调系统负荷计算方法、装置、设备及介质 Download PDFInfo
- Publication number
- CN113822569B CN113822569B CN202111097093.0A CN202111097093A CN113822569B CN 113822569 B CN113822569 B CN 113822569B CN 202111097093 A CN202111097093 A CN 202111097093A CN 113822569 B CN113822569 B CN 113822569B
- Authority
- CN
- China
- Prior art keywords
- building
- load
- heat
- heat accumulator
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004378 air conditioning Methods 0.000 title claims abstract description 67
- 238000004364 calculation method Methods 0.000 title claims abstract description 44
- 238000012546 transfer Methods 0.000 claims abstract description 56
- 238000000034 method Methods 0.000 claims abstract description 26
- 230000005855 radiation Effects 0.000 claims description 33
- 238000005338 heat storage Methods 0.000 claims description 13
- 230000002068 genetic effect Effects 0.000 claims description 12
- 238000004422 calculation algorithm Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 5
- 230000002093 peripheral effect Effects 0.000 claims description 2
- 239000000463 material Substances 0.000 description 18
- 238000005265 energy consumption Methods 0.000 description 10
- 238000005457 optimization Methods 0.000 description 10
- 230000017525 heat dissipation Effects 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 230000006870 function Effects 0.000 description 7
- 238000005286 illumination Methods 0.000 description 5
- 238000004134 energy conservation Methods 0.000 description 4
- 230000000704 physical effect Effects 0.000 description 4
- 238000009423 ventilation Methods 0.000 description 4
- 238000004590 computer program Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 230000003204 osmotic effect Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000010438 heat treatment Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000035772 mutation Effects 0.000 description 2
- 230000001681 protective effect Effects 0.000 description 2
- 235000002918 Fraxinus excelsior Nutrition 0.000 description 1
- 244000181980 Fraxinus excelsior Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000005034 decoration Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Air Conditioning Control Device (AREA)
Abstract
本发明公开了一种空调系统负荷计算方法、装置、设备及介质，其中，所述方法包括：获取建筑的蓄热体参数和内部热源参数；根据所述内部热源参数计算所述建筑的内部热源负荷；根据所述蓄热体参数以及预先设置的建筑蓄热体传热模型，计算所述建筑的蓄热体负荷；基于所述内部热源负荷和所述蓄热体负荷，确定所述建筑的空调系统总负荷。上述方法所提供的建筑蓄热体传热模型通用性强，对于不同的建筑空调系统，无需重新构建模型，仅需获取目标建筑的相关参数便可快速计算目标建筑的空调系统负荷，计算效率更高。
Description
技术领域
本发明涉及暖通空调技术领域，特别是涉及一种空调系统负荷计算方法、装置、设备及介质。
背景技术
空调系统的负荷计算对于其设计、控制优化和故障诊断都是至关重要的。随着建筑信息化模型(BIM)的广泛使用，空调系统设计所需的信息被高度数字化，已有研究和软件可以实现利用BIM进行负荷计算的功能。然而，采用BIM对空调系统进行逐时负荷计算非常繁琐且耗时，在这个过程中，无论是仿真模型EnergyPlus还是eQuest，都需要建立一个完整的建筑模型；另一方面，对于不同的建筑还需对应建立新的模型。
因此，提出一种具有通用性且计算效率高的空调系统负荷估算方法是十分必要的。
发明内容
针对上述技术问题，本发明提供一种空调系统负荷计算方法、装置、设备及介质，能够更高效地计算空调系统的总负荷。
第一方面，本发明提供一种空调系统负荷计算方法，包括：
获取建筑的蓄热体参数和内部热源参数；
根据所述内部热源参数计算所述建筑的内部热源负荷；
根据所述蓄热体参数以及预先设置的建筑蓄热体传热模型，计算所述建筑的蓄热体负荷；
基于所述内部热源负荷和所述蓄热体负荷，确定所述建筑的空调系统总负荷。
可选的，所述预先设置的建筑蓄热体传热模型基于RC热网络模型构建；其中，所述建筑蓄热体传热模型的参数通过遗传算法进行优化。
可选的，所述根据所述内部热源负荷和所述蓄热体负荷确定所述建筑的空调系统总负荷，具体包括：
计算所述内部热源负荷和所述蓄热体负荷的负荷总和；
根据所述负荷总和及所述建筑的内部热容得到所述建筑的空调系统总负荷。
可选的，所述建筑的蓄热体包括外围护结构和内部蓄热体。
可选的，所述内部热源参数包括：内部热源的逐时分项计量数据和实地采集数据。
第二方面，本发明提供一种空调系统负荷计算装置，包括：
获取模块，用于获取建筑的蓄热体参数和内部热源参数；
第一计算模块，用于根据所述内部热源参数计算所述建筑的内部热源负荷；根据所述蓄热体参数以及预先设置的建筑蓄热体传热模型，计算所述建筑的蓄热体负荷；
第二计算模块，用于基于所述内部热源负荷和所述蓄热体负荷，确定所述建筑的空调系统总负荷。
可选的，所述第二计算模块具体用于：
计算所述内部热源负荷和所述蓄热体负荷的负荷总和；
根据所述负荷总和及所述建筑的内部热容得到所述建筑的空调系统总负荷。
可选的，所述建筑蓄热体传热模型基于RC热网络模型构建；其中，所述建筑蓄热体传热模型的参数通过遗传算法进行优化。
第三方面，本发明提供一种数据处理设备，包括处理器，所述处理器和存储器耦合，所述存储器存储有程序，所述程序由所述处理器执行，使得所述数据处理设备执行第一方面所述的空调系统负荷计算方法。
第四方面，本发明还提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时，实现如上述第一方面所述的空调系统负荷计算方法。
相比现有技术，本发明的有益效果在于：
本发明提供的一种通用性强的空调系统负荷计算方法，不同于现有技术需利用EnergyPlus，eQuest、或TRNSYS等软件建立一个完整的二维或三维模型后，再进行空调负荷模拟，而仅需获取目标建筑的相关参数，并将所述参数输入至预设的建筑蓄热体传热模型中，便可快速计算得到目标建筑的空调系统的负荷，计算效率更高。
相应地，本发明还提供一种空调系统负荷计算装置、设备及存储介质。
附图说明
为了更清楚地说明本发明的技术方案，下面将对实施方式中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施方式，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的空调系统负荷计算方法流程图；
图2是本发明实施例提供的建筑空调系统负荷计算模型示意图；
图3是本发明实施例提供的建筑内部蓄热体传热模型示意图；
图4是本发明实施例提供的空调系统负荷计算装置结构框图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
空调系统是建筑中的主要能耗来源，因此，空调节能是建筑节能研究的一个重要方向。在空调节能的研究中，空调负荷计算对于空调节能而言意义重大，可用于确定空调系统的设备容量、系统参数及控制方案等，对此，第一方面，本发明一个实施例提供了一种空调系统负荷计算方法，如图1所示，所述方法具体包括下述步骤。
S1：获取建筑的蓄热体参数和内部热源参数。
建筑中的蓄热体无论在冬季还是夏季对于调节室内热环境和减少建筑能耗都有积极的意义。在冬季，建筑的蓄热构件可在白天储存太阳能，并在夜间室外温度非常低时将热量再慢慢释放到室内环境中，在夏季，结合夜间的通风降温，蓄热体可有效降低室内白天空气峰值温度或减小白天建筑的峰值负荷；建筑的内部热源则包括照明、电力设备、人员、新风和渗透风等；本实施例利用建筑蓄热体及内部热源的相关参数进行空调系统负荷的计算。
具体地，建筑的蓄热体包括外围护结构和内部蓄热体，蓄热体参数则对应细分为外围护结构相关参数以及内部蓄热体相关参数，其中，外围护结构相关参数包括外围护结构的面积、温度和材料物理性质参数，内部蓄热体相关参数则包括内部蓄热体的面积和材料物理性质参数。
可以理解的是，在本实施例中，建筑的外围护结构具体包括屋面、外墙和外窗，建筑的内部蓄热体包括楼板、内墙、家具、地毯和各种装饰等。
具体地，内部热源参数包括内部热源的逐时分项计量数据和实地采集数据，其中，内部热源的逐时分项计量数据具体为实时采集的电耗数据，实地采集数据则包括人员的数量和渗透风量等。
在一个实施例中，所获取的建筑的蓄热体参数和内部热源参数具体如下表1所示。
表1
其中，气象参数通过建立小型气象站获取，用于计算建筑外空气的焓值；建筑围护结构的外壁面温度、室外综合温度和室内温度及相对湿度等数据均为实际测量获得；建筑围护结构和内部结构的材料以及对应的物理性质参数通过调研和搜集资料获得。
具体地，建筑外围护结构中，外墙以及屋面材料及其物理性质参数如下表2所示。
表2
S2：根据所述内部热源参数计算所述建筑的内部热源负荷；根据所述蓄热体参数以及预先设置的建筑蓄热体传热模型，计算所述建筑的蓄热体负荷。
通过建筑内部热源的相关参数可计算建筑内部热源负荷，具体包括照明散热、电力设备散热、人员和新风负荷等。
在本实施例中，所述预先设置的建筑蓄热体传热模型基于RC热网络模型构建；其中，所述建筑蓄热体传热模型的参数通过遗传算法进行优化。通过输入蓄热体参数至建筑蓄热体传热模型中，可计算得到建筑的蓄热体负荷。
RC(热阻热容)热网络模型是一种常见的灰箱模型，其原理是利用简化的物理模型来模拟建筑能源系统能量流动，具有参数少、直观物理意义、高准确度的优点，可有效提高空调系统负荷的计算效率。
S3：基于所述内部热源负荷和所述蓄热体负荷，确定所述建筑的空调系统总负荷。
具体地，通过上述步骤得到内部热源负荷和蓄热体负荷后，计算二者的负荷总和，再根据所述负荷总和及所述建筑的内部热容得到建筑的空调系统总负荷。
在一个实施例中，空调系统总负荷可表示为下式：
其中，C为热容，单位为J/m2·K；T为温度；Q为负荷，单位为W/m2，Qest为计算而得的空调系统负荷，Qrf为屋面负荷，Qei为外墙负荷，Qwin为外窗负荷，Qim为建筑内部蓄热体产生的负荷，Qconv为照明灯具、电力设备和人员的对流换热负荷，Qfr为新风和渗透风产生的空调负荷，Qla为人员产生的潜热负荷；式中下标rf、ei、win、im和in分别代表屋面、外墙、外窗、建筑内部蓄热体和建筑内部。
在本实施例中，屋面负荷Qrf、外墙负荷Qei和外窗负荷Qwin均通过3R2C模型求解得到，建筑内部蓄热体产生的负荷Qim则通过改进后的2R2C模型求解获得，下面将对其求解过程进行详细介绍。
参照图2，可以理解的是，建筑外围护结构传热模型具体包括屋面和外墙的传热模型，以及外窗传热模型，其中，屋面和外墙具有储能作用，故采用简化的3R2C模型，且不同方向的外墙需分别计算；而对于外窗传热模型，该模型忽略其储能，故将其看作纯电阻。
具体地，建筑外围护结构传热模型的传热方程表示为下述公式：
其中，R为热阻，单位为m2·K/W。关于数字下标，R1，R3，R5分别代表外部、中部和内部蓄热体热阻；C2，C4则分别代表外部蓄热体热阻R1与中部蓄热体热阻R3间、中部蓄热体热阻R3与内部蓄热体热阻R5间的热容；由于板壁内外表面的空气边界层可视为纯热阻，故热容为零。
对于屋面和外墙的3R2C模型中的参数3个R和2个C，可通过遗传算法优化计算得到，具体地，如果已知屋顶和外墙每一层详细的材料和物性参数，则使用理论模型对其进行优化，使3R2C模型与理论模型的幅频特性与相频特性尽量相近；如果无法得到屋顶和外墙详细的材料及其物性参数，则采用实测数据对3R2C模型进行优化，即：通过实测的热流数据来推测3R2C模型中的参数。
以下将通过一个实施例描述利用遗传算法优化建筑外围护结构传热模型参数的过程。
具体地，建筑外围护结构的传热属于一维传热过程，因此，其传热模型可以表示为如下所示的关于变量s的拉普拉斯变化方程。
其中，T表示温度，q表示热流，M(s)表示整个壁面所有材料层的总系数矩阵，A(s)～D(s)表示矩阵中的元素，即传热方程组的系数。
可以理解的是，若将所有壁面材料层包括内外空气层分开表示，则可表示为下述方程。
若材料热容很小，则上式可简化为：其中，Ri为该材料的热阻。
因此，建筑与空气的内外表面的换热热阻分别为Rin和Rout，所以有：
因为建筑外围护结构传热模型方程有唯一解，即关于变量的拉普拉斯变化方程有唯一解，则有：
令：
则生成方程GX(jω)、GY(jω)和GZ(jω)，分别用于表示理论的建筑外围护结构传热模型中，建筑外围护结构的材料外壁面、材料内部和材料内壁面的频率特性，该频率特性可以表示为幅值特性和相位特性。
在本实施例中，用简化的3R2C模型表示建筑外墙和屋面的传热特性；具体地，简化模型包括外部、中部和内部三个部分，对简化模型进行拉普拉斯变换计算：
其中：
A'＝1+(C4R5+C2R3+C2R5)s+C4C2R5R3s2
B'＝-(R5+R3+R1)-(C4R3R1+C2R3R1+C2R5R1+C4R5R3)s-C4C2R5R3R1s2
C'＝-(C4+C2)s+C4C2R3s2
D'＝1+(C4R1+C4R3+C2R1)s+C4C2R3R1s2
所述3R2C模型也可以表达为下述形式：
其中：
G'X(s)＝A'(s)/B'(s)
G'Y(s)＝1/B'(s)
G'Z(s)＝D'(s)/B'(s)
令生成方程GX'(jω)、GY'(jω)和GZ'(jω)，这三个方程分别是简化的3R2C传热模型中，建筑外围护结构的外部、中部和内部的频率特性。
将GX'(jω)、GY'(jω)和GZ'(jω)与理论模型的GX(jω)、GY(jω)和GZ(jω)进行匹配，得到建筑外围护结构传热模型中的各参数值。
在本实施例中，建筑外围护结构传热模型优化的目的是使得简化的3R2C模型的传热幅频特性及相频特性与实际结构理论模型的传热幅频特性及相频特性尽量接近，因此，可推得目标函数为：
其中，PL表示相位，N表示频域，W为关于外部、中部和内部材料层传热的幅频特性与相频特性的权重比例；在本实施例中，W的三个权重比例均取值为1，以得到最好的优化效果。
具体地，优化上述目标函数的约束条件如下所示：
其中，R和C表示所有材料层总的热阻和热容，频域N的取值范围为在本实施例中，N等于10(n1-n2)+1，其中n1和n2分别取值为8和3。
进一步地，对建筑外围护结构的简化3R2C模型进行优化计算时，所需输入的参数包括：外墙和屋面的各组成层材料的热阻、比热容、厚度、密度，以及外部壁面温度和内部壁面温度。
为了求解最小函数值的参数，本实施例定义适应度函数为f3R2C(R1,R5,C4)＝J3R2C(R1,R5,C4)。
在本实施例中，使用遗传算法进行优化的基本运算过程如下：
1)选取参数随机值：给参数随机赋值，并且判断是否在约束条件范围内。
2)初始化：设置进化代数计数器i＝0，设置最大进化代数Max，随机生成若干个个体作为初始种群。
3)个体评价：计算初始种群中各个个体的适应度。
4)选择运算：将选择算子作用于群体；选择目的是将优化的个体直接遗传到下一代，或通过配对交叉产生新的个体再遗传到下一代。选择操作建立在群体中个体的适应度评估基础上。
5)交叉运算：将交叉算子作用于群体，交叉算子是遗传算法的核心。
6)变异运算：将变异算子作用于群体，即是对群体中的个体的某些基因的基因值作变动。
7)群体经过选择、交叉、变异运算之后得到下一代群体。
8)终止条件判断：若i＝Max，则以进化过程中所得到的具有最大适应度个体作为最优解输出，终止计算。
参照图2和图3，建筑的内部蓄热体传热模型具体划分为三种，分别为模型1、模型2和模型3。
对于模型1，可理解为对建筑内部蓄热体进行分类，对蓄热体接受的辐射情况不进行分类。在模型1中，根据材料的热惰性，建筑内部蓄热体被分为了三种类型：轻型蓄热体、中型蓄热体和重型蓄热体；由于模型1没有对蓄热体接受到的辐射得热进行分类，故认为所得的太阳辐射、人员、照明和设备的辐射是均匀分布在不同种类的蓄热体上，因此，模型1存在以下关系：
Rim,11≠Rim,12≠Rim,13
Cim,11≠Cim,12≠Cim,13
此外，还有Qr,11＝Qr,12＝Qr,13，即三类蓄热体辐射得热分别等于室内总辐射得热的三分之一。
具体地，模型1的传热方程表示为下述公式：
其中，第一位数字下标1表示模型1，第二位数字分别代表三类不同类型的蓄热体。
对于模型2，可理解为对蓄热体接收的辐射情况进行分类，对建筑内部蓄热体不进行分类。在模型2中，根据内部蓄热体所接收到的辐射得热类型，可将建筑内部蓄热体分为三类：仅接受太阳辐射的蓄热体、接受其他辐射(人员、照明灯具和电力设备等产生的辐射)的蓄热体以及不接受辐射的蓄热体；由于模型2未考虑内部蓄热体本身具有的热惰性差别，故存在下述关系：
Rim,21＝Rim,22＝Rim,23
Cim,21＝Cim,22＝Cim,23
此外，还有Qr,21≠Qr,22，其中，Qr,21等于太阳直射辐射，Qr,22等于太阳直射辐射加上室内辐射。
具体地，模型2的传热方程表示为下述公式：
其中，第一位数字下标2代表模型2，第二位数字表示不同类型的蓄热体。
对于模型3，可理解为对建筑内部蓄热体和蓄热体所接收的辐射情况都进行分类。在模型3中，根据材料的热惰性将建筑内部蓄热体分为两类：轻型蓄热体和重型蓄热体，不同于模型1，模型3不仅考虑到蓄热体可能出现温度分布不均匀现象(如蓄热体表面不同部位温度不一致，或蓄热体表面温度与蓄热体内部温度不一致等)，还考虑到不同类型的蓄热体所接收到的辐射类型不同，故设置为串并联混合模型。
具体地，模型3认为太阳辐射只作用在小部分重型蓄热体上(如靠近窗户的地面)，而其他辐射(人员、照明灯具和电力设备等产生的辐射)则作用在其余重型蓄热体和轻型蓄热体上，故存在下述关系：
Rim,31＝Rim,33≠Rim,35，Rim,32＝Rim,34≠Rim,36
Cim,31＝Cim,33≠Cim,35，Cim,32＝Cim,34≠Cim,36
此外，还有Qr,31≠Qr,32＝Qr,33，其中，Qr,31等于太阳直射辐射，Qr,32和Qr,33分别等于太阳直射辐射与室内辐射的总和的一半。
具体地，模型3的传热方程表示为下述公式：
/>
其中，第一位数字下标3代表模型3；第二位数字代表不同类型的蓄热体，具体地，R和C的下标1、3、5分别代表蓄热体表面，下标2、4、6则分别代表蓄热体内部。
上述建筑内部蓄热体传热模型中，模型1、模型2和模型3的参数R和C均通过遗传算法优化计算得到，从而使模型得到的负荷与实测得到的实际负荷尽可能相近，在另一实施例中，可通过部分实测数据来估算建筑内蓄热体传热模型的相关参数。
具体地，建筑内部蓄热体的参数优化目的是使建筑空调系统计算的负荷值与实际负荷值之间的均方根误差尽可能小，所以优化目标函数设置为:
其中，Qact表示建筑空调系统实际的逐时负荷值，Qest表示建筑空调系统计算的逐时负荷值，Rim,1～Rim,m和Cim,1～Cim,n分别表示建筑内部蓄热体的热阻和热容，m和n的取值可根据选取的模型确定。
需要说明的是，建筑内部蓄热体传热模型可采用实测的冷量数据对模型参数进行预估，并采用遗传算法进行优化，具体优化过程与上述建筑外围护传热模型参数的优化过程思路一致，本发明不再赘述。
在实际应用中，选择建筑内部蓄热体的三类模型中效果最好(即平均相对误差最小)的一类进行建筑空调系统负荷的计算。
经分析，实测数据的优化结果对辐射得热的分布敏感，对蓄热体的分类不敏感，推荐使用模型2；而模拟数据的优化结果则对辐射得热的分布不敏感，对蓄热体的分类非常敏感，推荐使用模型1。
考虑更加复杂的模型并没有表现出更好的估算效果，因而简单的3R3C并联模型便可以对建筑内部蓄热体进行很好的描述。
对于内部热源负荷，本实施例转为计算建筑内部热源得热，具体包括两个部分，一是照明及电力设备散热计算，二是人员、新风和渗透风负荷计算。
其中，对于照明及电力设备散热的计算，由于建筑内所有电力设备最终都将电耗转化为内部负荷，故可以将逐时分项计量数据作为内部的热量计算，则建筑不同日类型下的逐时分项能耗数据计算模型为：
其中，QL为照明负荷；EL为照明分项能耗或动力分项能耗；a为能耗平均值；2πωn为以天为周期下的谐波频率；ε为误差；h表示小时；δn和ηn均为系数。
需要说明的是，电力设备负荷QE的计算与照明负荷QL相同，即：
具体地，根据回归模型对不同日类型下的逐时分项能耗数据进行计算；对于同一日类型，则选取能耗最大的值作为基准值，其他时刻的能耗值与基准值的比值为对应时刻的同时使用系数，利用所述同时使用系数计算各时刻的能耗值。
需要说明的是，照明灯具和电力设备的对流换热得热以及辐射得热的比例可根据样品手册确定，若没有样品手册作为参考，则可参考ASHRAE Handbook或者相关研究数据。
具体地，本实施例采用的照明灯具的对流换热比例为20％，辐射比例为80％；电力设备的对流换热比例为80％，辐射比例为20％，依据该比例可计算出相应的对流换热得热。
示例性的，有QLC＝QL·γ1，QEC＝QE·γ1；其中，QLC和QEC分别表示照明灯具和电力设备的对流换热得热，γ1表示照明灯具或电力设备的对流换热得热所占比例值。
在人员负荷的计算中，具体地，对于有人员出入的建筑，可选取人员数量的高峰值作为基准值，其他时刻的人员数量与基准值的比值作为该时刻的人员同时在室系数，利用所述人员同时在室系数计算各时刻的人体负荷；对于没有人员进出的建筑，则可以进行短期现场调研或者参考相关规范。
需要说明的是，人员散热功率的选取可参考相关标准，例如，显热量、潜热量以及对流换热得热和辐射得热的比例可参考ASHRAE Handbook，取值分别为40％、20％和40％；又如，可根据Braun的研究，将人员散热的对流换热比例设置为70％，辐射比例设置为30％。
示例性的，有QPC＝QP·γ2,其中，QP表示人体散热量，n表示建筑内人员数量，/>表示群集系数，q表示不同室温和劳动性质时的成年男子散热量，QPC表示人体的对流换热部分负荷，γ2表示人体对流换热得热比例。
综上，有Qconv＝QLC+QEC+QPC。
在新风负荷和渗透风负荷的计算中，具体地，新风量可根据实际空调系统的新风设置情况进行设置；渗透风量的设置可根据建筑类型和分区在0.1～0.5次/小时范围内选取。
具体地，有：
Qi＝0.28mVρ(tw-tn)
Qf＝G(hw-hn)
Qfr＝Qi+Qf
其中，m为渗透风换气次数，V为房间体积，ρ为室外干球温度下的空气密度，tw为室外干球温度；tn为室内干球温度，Qi则为渗透风负荷；G为新风量，Qf新风负荷；Qfr则为新风和渗透风产生的空调负荷；hw为室外空气焓；hn为室内空气焓。
本发明上述实施例在现有的RC热网络模型基础上，建立了建筑外围护结构传热模型，并对建筑内部蓄热体的简化模型进行了三种不同类型的改进，同时还使用内部热源的实时分项计量数据对建筑内部的照明、电力设备等内部得热进行描述，相对比现有技术能够有效提高计算效率，另一方面，所构建的建筑蓄热体传热模型具备通用性，对于不同的建筑空调系统，仅需输入对应参数便可得到负荷计算结果。
第二方面，参照图4，本发明另一实施例提供一种空调系统负荷计算装置，包括获取模块101、第一计算模块102和第二计算模块103。
获取模块101用于获取建筑的蓄热体参数和内部热源参数。
第一计算模块102用于根据所述内部热源参数计算所述建筑的内部热源负荷；根据所述蓄热体参数以及预先设置的建筑蓄热体传热模型，计算所述建筑的蓄热体负荷。
第二计算模块103用于基于所述内部热源负荷和所述蓄热体负荷，确定所述建筑的空调系统总负荷。
上述装置内的各模块之间信息交互、执行过程等内容，由于与本发明第一方面的方法实施例基于同一构思，具体内容可参见本发明方法实施例中的叙述，此处不再赘述。
以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的模块可以是或者也可以不是物理上分开的，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方法的目的。
第三方面，本发明提供一种数据处理设备，包括处理器，所述处理器和存储器耦合，所述存储器存储有程序，所述程序由所述处理器执行，使得所述数据处理设备执行第一方面所述的空调系统负荷计算方法。
第四方面，本发明还提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时，实现如上述第一方面所述的空调系统负荷计算方法。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于计算机可监听存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory，ROM)或随机存储记忆体(Random AccessMemory，RAM)等。
以上所述是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也视为本发明的保护范围。
Claims (10)
1.一种空调系统负荷计算方法，其特征在于，包括：
获取建筑的蓄热体参数和内部热源参数；
根据所述内部热源参数计算所述建筑的内部热源负荷；
根据所述蓄热体参数以及预先设置的建筑蓄热体传热模型，计算所述建筑的蓄热体负荷；
基于所述内部热源负荷和所述蓄热体负荷，确定所述建筑的空调系统总负荷；其中，空调系统总负荷的计算包括对建筑内部蓄热体产生的负荷的计算；建筑内部蓄热体产生的负荷通过建筑的内部蓄热体传热模型的三种模型中效果最好的一种确定；
其中，所述建筑的内部蓄热体传热模型具体划分为三种，分别为模型1、模型2和模型3；模型1表征对建筑内部蓄热体进行分类且对蓄热体接受的辐射情况不进行分类的情况；模型2表征对建筑内部蓄热体不进行分类且对蓄热体接收的辐射情况进行分类的情况；模型3表征对建筑内部蓄热体和蓄热体所接收的辐射情况都进行分类的情况。
2.根据权利要求1所述的空调系统负荷计算方法，其特征在于，所述预先设置的建筑蓄热体传热模型是基于RC热网络模型构建；
其中，所述建筑蓄热体传热模型的参数通过遗传算法进行优化。
3.根据权利要求1所述的空调系统负荷计算方法，其特征在于，所述根据所述内部热源负荷和所述蓄热体负荷确定所述建筑的空调系统总负荷，具体包括：
计算所述内部热源负荷和所述蓄热体负荷的负荷总和；
根据所述负荷总和及所述建筑的内部热容得到所述建筑的空调系统总负荷。
4.根据权利要求1～3任一所述的空调系统负荷计算方法，其特征在于，所述建筑的蓄热体包括外围护结构和内部蓄热体。
5.根据权利要求1～3任一所述的空调系统负荷计算方法，其特征在于，所述内部热源参数包括：
内部热源的逐时分项计量数据和实地采集数据。
6.一种空调系统负荷计算装置，其特征在于，包括：
获取模块，用于获取建筑的蓄热体参数和内部热源参数；
第一计算模块，用于根据所述内部热源参数计算所述建筑的内部热源负荷；根据所述蓄热体参数以及预先设置的建筑蓄热体传热模型，计算所述建筑的蓄热体负荷；
第二计算模块，用于基于所述内部热源负荷和所述蓄热体负荷，确定所述建筑的空调系统总负荷；其中，空调系统总负荷的计算包括对建筑内部蓄热体产生的负荷的计算；建筑内部蓄热体产生的负荷通过建筑的内部蓄热体传热模型的三种模型中效果最好的一种确定；其中，所述建筑的内部蓄热体传热模型具体划分为三种，分别为模型1、模型2和模型3；模型1表征对建筑内部蓄热体进行分类且对蓄热体接受的辐射情况不进行分类的情况；模型2表征对建筑内部蓄热体不进行分类且对蓄热体接收的辐射情况进行分类的情况；模型3表征对建筑内部蓄热体和蓄热体所接收的辐射情况都进行分类的情况。
7.根据权利要求6所述的空调系统负荷计算装置，其特征在于，所述第二计算模块具体用于：
计算所述内部热源负荷和所述蓄热体负荷的负荷总和；
根据所述负荷总和及所述建筑的内部热容得到所述建筑的空调系统总负荷。
8.根据权利要求6所述的空调系统负荷计算装置，其特征在于，所述建筑蓄热体传热模型基于RC热网络模型构建；其中，
所述建筑蓄热体传热模型的参数通过遗传算法进行优化。
9.一种数据处理设备，其特征在于，包括：
处理器，所述处理器和存储器耦合，所述存储器存储有程序，所述程序由所述处理器执行，使得所述数据处理设备执行如权利要求1～5中任一项所述的空调系统负荷计算方法。
10.一种计算机存储介质，其特征在于，所述计算机存储介质存储有计算机指令，所述计算机指令用于执行上述权利要求1～5中任一项所述的空调系统负荷计算方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111097093.0A CN113822569B (zh) | 2021-09-17 | 2021-09-17 | 一种空调系统负荷计算方法、装置、设备及介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111097093.0A CN113822569B (zh) | 2021-09-17 | 2021-09-17 | 一种空调系统负荷计算方法、装置、设备及介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113822569A CN113822569A (zh) | 2021-12-21 |
| CN113822569B true CN113822569B (zh) | 2024-05-24 |
Family
ID=78914815
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111097093.0A Active CN113822569B (zh) | 2021-09-17 | 2021-09-17 | 一种空调系统负荷计算方法、装置、设备及介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113822569B (zh) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2011163710A (ja) * | 2010-02-12 | 2011-08-25 | Ei Engineering:Kk | 原単位算出システム及びこれを実行させるためのプログラム並びこのプログラムを記録した記録媒体 |
| CN105605733A (zh) * | 2015-12-25 | 2016-05-25 | 深圳达实智能股份有限公司 | 一种空调冷机的电网响应方法及装置 |
| CN105627506A (zh) * | 2015-12-25 | 2016-06-01 | 深圳达实智能股份有限公司 | 一种空调冷机的建筑冷负荷预测方法及装置 |
| US10339232B1 (en) * | 2015-02-25 | 2019-07-02 | Clean Power Research, L.L.C. | Computer-implemented system and method for modeling building heating energy consumption |
| CN110097217A (zh) * | 2019-04-16 | 2019-08-06 | 天津大学 | 一种基于等效rc模型的建筑动态室温预测方法 |
| CN110489908A (zh) * | 2019-08-27 | 2019-11-22 | 上海电气分布式能源科技有限公司 | 建筑物冷热负荷评估方法及系统、设备和介质 |
| CN112054505A (zh) * | 2020-07-30 | 2020-12-08 | 上海理工大学 | 一种基于建筑基础信息的建筑电力需求弹性快速量化方法 |
| CN113363974A (zh) * | 2021-06-16 | 2021-09-07 | 广东电网有限责任公司 | 基于累积电量低频采样的居民负荷构成解析方法及装置 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2020107851A1 (zh) * | 2018-11-29 | 2020-06-04 | 天津大学 | 基于既有大型公共建筑空调系统的低成本调适方法及调适系统 |
-
2021
- 2021-09-17 CN CN202111097093.0A patent/CN113822569B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2011163710A (ja) * | 2010-02-12 | 2011-08-25 | Ei Engineering:Kk | 原単位算出システム及びこれを実行させるためのプログラム並びこのプログラムを記録した記録媒体 |
| US10339232B1 (en) * | 2015-02-25 | 2019-07-02 | Clean Power Research, L.L.C. | Computer-implemented system and method for modeling building heating energy consumption |
| CN105605733A (zh) * | 2015-12-25 | 2016-05-25 | 深圳达实智能股份有限公司 | 一种空调冷机的电网响应方法及装置 |
| CN105627506A (zh) * | 2015-12-25 | 2016-06-01 | 深圳达实智能股份有限公司 | 一种空调冷机的建筑冷负荷预测方法及装置 |
| CN110097217A (zh) * | 2019-04-16 | 2019-08-06 | 天津大学 | 一种基于等效rc模型的建筑动态室温预测方法 |
| CN110489908A (zh) * | 2019-08-27 | 2019-11-22 | 上海电气分布式能源科技有限公司 | 建筑物冷热负荷评估方法及系统、设备和介质 |
| CN112054505A (zh) * | 2020-07-30 | 2020-12-08 | 上海理工大学 | 一种基于建筑基础信息的建筑电力需求弹性快速量化方法 |
| CN113363974A (zh) * | 2021-06-16 | 2021-09-07 | 广东电网有限责任公司 | 基于累积电量低频采样的居民负荷构成解析方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113822569A (zh) | 2021-12-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wang et al. | Parameter estimation of internal thermal mass of building dynamic models using genetic algorithm | |
| Homod | Review on the HVAC system modeling types and the shortcomings of their application | |
| Balaras | The role of thermal mass on the cooling load of buildings. An overview of computational methods | |
| Wang et al. | Simplified building model for transient thermal performance estimation using GA-based parameter identification | |
| Li et al. | Fast bidirectional building performance optimization at the early design stage | |
| Kalogirou | Applications of artificial neural-networks for energy systems | |
| Homod et al. | Energy savings by smart utilization of mechanical and natural ventilation for hybrid residential building model in passive climate | |
| Magnier et al. | Multiobjective optimization of building design using TRNSYS simulations, genetic algorithm, and Artificial Neural Network | |
| CN108006887A (zh) | 既有建筑空调用能模型构建方法及空调能耗预测方法 | |
| CN102238763B (zh) | 基站能耗计算系统及方法 | |
| Carlucci et al. | An optimization procedure based on thermal discomfort minimization to support the design of comfortable net zero energy buildings | |
| CN110097217B (zh) | 一种基于等效rc模型的建筑动态室温预测方法 | |
| CN111310257A (zh) | 一种bim环境下的区域建筑能耗预测方法 | |
| Arendt et al. | Influence of input data on airflow network accuracy in residential buildings with natural wind-and stack-driven ventilation | |
| Shen et al. | Exploring potential for residential energy saving in New York using developed lightweight prototypical building models based on survey data in the past decades | |
| Xu et al. | Implicit method for solving building heat transfer model and its application in energy-saving materials | |
| CN114037165A (zh) | 一种考虑多源不确定性的建筑冷热负荷预测方法 | |
| CN113822569B (zh) | 一种空调系统负荷计算方法、装置、设备及介质 | |
| Sreshthaputra et al. | Improving building design and operation of a Thai Buddhist temple | |
| CN110737937A (zh) | 基于反射涂层和垂直绿化的既有建筑外墙优化应用方法 | |
| He et al. | Simplified building thermal model used for optimal control of radiant cooling system | |
| Pan et al. | Multi-objective optimization for building performance design considering thermal comfort and energy consumption | |
| Kalogeras et al. | Verification and validation of a simulation model for energy use in buildings | |
| Alam et al. | Energy simulation to estimate building energy consumption using EnergyPlus | |
| CN113341792A (zh) | 一种基于bim的建筑数据采集方法及系统 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant |