CN113591410A - 针对博物馆展陈空间气流组织可视化的监测数据同化方法 - Google Patents

针对博物馆展陈空间气流组织可视化的监测数据同化方法 Download PDF

Info

Publication number
CN113591410A
CN113591410A CN202110902814.4A CN202110902814A CN113591410A CN 113591410 A CN113591410 A CN 113591410A CN 202110902814 A CN202110902814 A CN 202110902814A CN 113591410 A CN113591410 A CN 113591410A
Authority
CN
China
Prior art keywords
value
condition
parameter
condition parameter
air state
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.)
Granted
Application number
CN202110902814.4A
Other languages
English (en)
Other versions
CN113591410B (zh
Inventor
钱玮昕
唐铭
刘京
陈庆庆
刘洋
丁荣良
杜雷
张然
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NATIONAL MUSEUM OF CHINA
Harbin Institute of Technology
Original Assignee
NATIONAL MUSEUM OF CHINA
Harbin Institute of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by NATIONAL MUSEUM OF CHINA, Harbin Institute of Technology filed Critical NATIONAL MUSEUM OF CHINA
Priority to CN202110902814.4A priority Critical patent/CN113591410B/zh
Publication of CN113591410A publication Critical patent/CN113591410A/zh
Application granted granted Critical
Publication of CN113591410B publication Critical patent/CN113591410B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D21/00—Measuring or testing not otherwise provided for
    • G01D21/02—Measuring two or more variables by means not covered by a single other subclass
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Air Conditioning Control Device (AREA)

Abstract

针对博物馆展陈空间气流组织可视化的监测数据同化方法，涉及一种博物馆内的气流组织的监测数据同化方法。为了解决现有的博物馆内空气状态参数的监测方法使得流场可视化的结果与实际观测值产生较大误差的问题。本发明通过样本的正向模拟计算出边界条件与模拟结果之间的相关关系；然后基于迭代集合卡尔曼滤波，结合有限监测点的观测数据，计算各误差在模拟值与测量值的差值中所占的权重，反向推导出逼近实际工况的边界条件参数，并将优化后的边界条件参数重新代入CFD模型进行计算，循环迭代直至模拟结果与观测数据的差距达到收敛条件，获得能够反映博物馆展陈空间内较为真实的空气状态参数分布结果。主要用于博物馆展陈空间气流组织的监测。

Description

针对博物馆展陈空间气流组织可视化的监测数据同化方法
技术领域
本发明涉及一种博物馆内的气流组织的监测数据同化方法。
背景技术
温湿度等空气状态参数的监测和控制对于博物馆藏品保护非常重要。自2008年起，国内博物馆陆续建立了展陈环境温湿度监测系统，通过放置在展厅、展柜等不同位置的温湿度传感器，对文物保存环境中的温湿度数据进行连续监测和记录。但存在以下问题：展陈空间内能够放置的温湿度传感器数量始终是有限的，无法准确掌握热湿流场的空间全局分布情况，从而无法获知不同藏品所在位置的温湿度及变化情况。目前，可采用计算流体动力学(CFD)模拟来解决上述问题。利用室内环境CFD模拟，即在热湿负荷及空调设计参数已知的基础上实现流场可视化的数值模拟技术。然而在实际应用中，室外气象参数、照明和设备开启率以及观展人数的动态变化等因素导致展厅内的热湿负荷存在波动，同时，不同工况下的空调运行参数也有所变化，从而流场可视化的结果与实际观测产生较大差异。
发明内容
本发明是为了解决现有的博物馆内空气状态参数的监测方法使得流场可视化的结果与实际观测值产生较大误差的问题；进而提供一种针对博物馆展陈空间气流组织可视化的监测数据同化方法。
针对博物馆展陈空间气流组织可视化的监测数据同化方法，包括以下步骤：
S1、针对展陈空间内影响热湿环境的因素，确定模型边界条件中具有不确定性的条件参数变量，统计具有不确定性的条件参数变量的数量，并输入到数据同化算法程序中，同时向数据同化算法程序中输入各条件参数变量对应的条件参数值的区间范围；
根据每个条件参数变量在其变化区间内的一般分布规律，对各条件参数变量分别进行采样，并构建所有条件参数变量的初始集合E；
S2、向数据同化算法程序中输入展陈空间内无线监测系统所采集的观测数据，构建测量值向量；
S3、通过CFD正演模型分别对初始集合中的各元素进行求解，获得观测节点位置对应的空气状态参数模拟值矩阵；
S4、基于各条件参数变量的取样样本值与各空气状态参数模拟值之间的协方差矩阵、各空气状态参数模拟值之间的协方差矩阵以及测量误差方差矩阵，计算卡尔曼增益矩阵；
S5、通过集合卡尔曼滤波法修正集合元素中的条件参数变量的参数值并更新初始集合为E'；
S6、计算空气状态参数模拟值与空气状态参数测量值之间的均方根误差，设定收敛标准，判定计算收敛与否；当所有空气状态参数的RMSE均达到收敛标准时，则判定计算收敛，执行步骤S7；否则，令E＝E'，返回步骤S3，循环执行步骤S3～S6，直至迭代计算收敛；
S7、输出监测数据同化后博物馆展陈空间内气流组织的可视化结果。
进一步地，S1中，初始集合E的计算过程如下：
设模型边界条件中存在m个具有不确定性的条件参数变量xi，i＝1,2,3，…，m；设每个条件参数变量的取样数量为n，即初始集合中元素的数量为n，初始集合中的第j个元素为Zj＝[x1,j，x2,j，x3,j，…，xm,j]，j＝1,2,3,…,n；构成初始集合E＝[Z1，Z2，Z3，···，Zn]，对应的矩阵形式为：
Figure BDA0003200349740000021
式中，xi,j——第i个条件参数变量的第j个取样样本值。
进一步地，S2中，构建测量值向量的计算过程如下：
设数据同化算法程序中所构建的测量值向量为Y，假设在一次测量中，在c个测量点测量d种空气状态参数，则观测数据共包含b＝c×d个空气状态参数测量值ya，a＝1,2,3，…，b；则测量值向量Y为：
Y＝[y1 y2 … yb]T (2)。
进一步地，S3中，构建模拟值矩阵的具体计算过程如下：
将模型中的边界条件设定为集合E中每个元素的样本值，一组边界条件为一种工况，将集合E中的元素Zj分别代入CFD正演模型中，通过CFD模拟，对每个工况分别进行求解计算，获得各个工况对应的展陈空间流场分布模拟结果，构建模拟值矩阵P如下：
Figure BDA0003200349740000023
式中，pa,j——第j个集合元素对应的第a个空气状态参数模拟值。
进一步地，S4中，卡尔曼增益矩阵的具体计算过程如下：
K＝Cxp(Cpp+R)-1 (4)
Figure BDA0003200349740000031
Figure BDA0003200349740000032
Figure BDA0003200349740000033
Figure BDA0003200349740000034
Figure BDA0003200349740000035
Figure BDA0003200349740000036
式中，Cxp——集合中的每个条件参数变量的取样样本值与各空气状态参数模拟值的协方差矩阵；Cxip——集合中的第i个条件参数变量的取样样本值与各空气状态参数模拟值的协方差矩阵；Cpp——各空气状态参数模拟值之间的协方差矩阵；R——测量误差方差矩阵，测量误差方差矩阵是由无线监测系统所采集的各测量值自身的背景噪声所组成的方差矩阵，其中σa——第a个空气状态参数测量值标准差，由测量仪器误差限ηa计算得，σa＝ηa/3；
Figure BDA0003200349740000041
——输入的第i个条件参数变量n个样本的平均值；
Figure BDA0003200349740000042
——第a个空气状态参数n个样本模拟值的平均值。
进一步地，S5中，修正后的初始集合E’、修正后的初始集合元素Z'j和修正后的条件参数变量样本x'i,j的计算过程如下：
E'＝E+K(YH-P) (11)
式中，H——将空气状态参数测量值投影至空气状态参数模拟值的投影矩阵，为1×n的元素均为1的矩阵，从而使空气状态参数测量值矩阵转化为b×n的矩阵；
则集合中的条件参数变量样本xi,j修正并更新如下：
Figure BDA0003200349740000043
Figure BDA0003200349740000044
进一步地，S6中，空气状态参数模拟值与空气状态参数测量值之间的均方根误差和收敛准则计算过程如下：
第a个空气状态参数模拟值与空气状态参数测量值的均方根误差计算如下：
Figure BDA0003200349740000045
设置收敛准则为：
RMSEa≤|Da|2 (15)
式中，Da——第a个空气状态参数的测量值的最大误差。
进一步地，S1中，判定模型边界条件中具有不确定性的条件参数变量的方法为：针对展陈空间内影响热湿环境的各个因素，分析其在展陈空间使用过程中或特定的研究时间段内的数值大小是否在一定区间范围内波动变化，若该因素波动变化，则判定其为具有不确定性的条件参数变量，否则判定其为常量。
进一步地，S1中，所述的条件参数包括展陈空间内设备散热量、人体散热量、围护结构传热量、热源温度及空调运行参数。
进一步地，S1中，所述的条件参数值的取值范围根据历史数据或设计值给出。
有益效果：
对于建筑空间流场的模拟，传统方法是根据室内热湿负荷设计值(或估计值)及空调系统设计参数设置CFD计算模型边界条件，当计算结果与实际测量不符时，根据测量值与模拟值的差异，以试错的方式人工调节模型边界条件等参数。不同于上述方法，本发明将CFD模拟技术与集合卡尔曼滤波法相耦合，基于IEnKF原理，通过将监测系统的测量数据与CFD模拟结果相融合，并且不同于已有技术中直接通过卡尔曼滤波计算修正模拟结果，为了保证流场最终计算结果满足计算流体动力学控制方程的约束，该技术利用集合卡尔曼滤波法反向修正CFD计算模型的边界条件参数，并将CFD正向模拟与集合卡尔曼滤波反向修正边界条件的过程进行循环迭代，最终获得优化后接近实际工况的边界条件参数对应的流场模拟结果。该技术有效提高了模拟计算的精度，优化了博物馆展陈空间气流组织可视化的输出结果，同时，相比于传统方法，大大减少了人工试错的时间成本，具有灵活性强，实用性高等优点。针对对文物保存环境有极高要求的博物馆展陈空间，该技术通过有限数量监测点的数据，掌握温湿度空间场分布规律，更加准确地反映空间内的空气状态参数，实现对展陈空间内不同位置温湿度情况的监测，从而更好地指导展品的位置摆放和空调系统的运行调节，为博物馆藏品预防性保护和风险防控提供技术支撑。
通过实验验证：送风风速的均方根误差减小到0.02m/s，为初始误差的10％，送风温度的均方根误差减小到0.1℃，为原来的2.5％，热源温度均方根误差减小到0.4℃，为原来的8％，说明采用本发明的方法大大减小了流场可视化的结果与实际观测值的误差。
附图说明
图1为基于监测数据同化的展陈空间气流组织可视化技术流程图；
图2为实施方式中所建立的理想空间结构立面图；
图3为边界条件参数数据同化效果图；其中，图3(a)为送风口1风速数据同化结果；图3(b)为送风口2风速数据同化结果；图3(c)为送风口温度数据同化结果；图3(d)为热源温度数据同化结果；
图4为优化前后气流组织可视化结果对比图；其中，图4(a)为优化前样本平均空间立面温度云图；图4(b)为优化后样本平均空间立面温度云图；图4(c)为优化前样本平均空间立面风速云图；图4(d)为优化后样本平均空间立面风速云图。
具体实施方式
如图1所示，本实施方式所述的针对博物馆展陈空间气流组织可视化的监测数据同化方法，包括以下步骤：
S1、针对展陈空间内影响热湿环境的各个因素(冷、热、湿源等)，分析其在展陈空间使用过程中(或特定的研究时间段内)大小是否在一定区间范围内波动变化，该因素若波动变化，则判定其为具有不确定性的条件参数变量，否则判定其为常量，从而确定模型边界条件中所有具有不确定性的条件参数变量，条件参数变量是影响室内环境的边界条件的变量，条件参数包括展陈空间内设备散热量、人体散热量、围护结构传热量、热源温度及空调运行参数(如送风风速、送风风温)等；
统计具有不确定性的条件参数变量的数量，并输入到数据同化算法程序中，根据先验信息向数据同化算法程序中输入各条件参数变量对应的条件参数值的区间范围(各条件参数值的取值范围根据历史数据或设计值给出)；
根据每个条件参数变量在其变化区间内的一般分布规律(通常为高斯分布，对于空调运行参数则是根据历史运行数据总结得到分布规律)，对各条件参数变量分别进行采样，并构建所有条件参数变量的初始集合，初始集合中的每个元素是由每个条件参数变量的单次取样样本共同组成；
设模型边界条件中存在m个具有不确定性的条件参数变量xi，i＝1,2,3，…，m；设每个条件参数变量的取样数量为n(也可以说是取样次数)，即初始集合中元素的数量为n，初始集合中的第j个元素为Zj＝[x1,j，x2,j，x3,j，…，xm,j]，j＝1,2,3,…,n；构成初始集合E＝[Z1，Z2，Z3，···，Zn]，对应的矩阵形式：
Figure BDA0003200349740000061
式中，xi,j——第i个条件参数变量的第j个取样样本值。
S2、向数据同化算法程序中输入展陈空间内无线监测系统所采集的观测数据，构建测量值向量；设数据同化算法程序中所构建的测量值向量为Y，假设在一次测量中，在c个测量点测量d种空气状态参数(空气流速、温度、湿度等)，则观测数据共包含b＝c×d个空气状态参数测量值ya，a＝1,2,3，…，b；则测量值向量Y为：
Y＝[y1 y2 … yb]T (2)。
S3、通过计算流体动力学CFD正演模型，分别对初始集合中的各元素进行求解，获得观测节点位置对应的空气状态参数模拟值矩阵：
将模型中的边界条件设定为集合E中每个元素的样本值，集合中的一个元素对应一组边界条件，一组边界条件为一种工况，也就是说从而在已知的变化范围内产生n种工况，将集合E中的元素Zj分别代入CFD正演模型中，通过CFD模拟，对每个工况分别进行求解计算，获得各个工况对应的展陈空间流场分布模拟结果；每个元素代入CFD计算模型求解之后都会得到空气流速、温度和湿度等结果，从而输出展陈空间内布置监测传感器的位置所对应的空气状态参数模拟值，构建模拟值矩阵P；其中，模拟输出的空气状态参数的种类与监测系统所采集的空气状态参数种类一致，构建模拟值矩阵P如下：
Figure BDA0003200349740000071
式中，pa,j——第j个集合元素对应的第a个空气状态参数模拟值。
S4、基于各条件参数变量的取样样本值与各空气状态参数模拟值之间的协方差矩阵、各空气状态参数模拟值之间的协方差矩阵以及测量误差方差矩阵，计算卡尔曼增益矩阵K；
K＝Cxp(Cpp+R)-1 (4)
Figure BDA0003200349740000072
Figure BDA0003200349740000073
Figure BDA0003200349740000081
Figure BDA0003200349740000082
Figure BDA0003200349740000083
Figure BDA0003200349740000084
式中，Cxp——集合中的每个条件参数变量的取样样本值与各空气状态参数模拟值的协方差矩阵；Cxip——集合中的第i个条件参数变量的取样样本值与各空气状态参数模拟值的协方差矩阵；Cpp——各空气状态参数模拟值之间的协方差矩阵；R——测量误差方差矩阵，测量误差方差矩阵是由无线监测系统所采集的各测量值自身的背景噪声所组成的方差矩阵，其中σa——第a个空气状态参数测量值标准差，由测量仪器误差限ηa计算得，σa＝ηa/3；
Figure BDA0003200349740000085
——输入的第i个条件参数变量n个样本的平均值；
Figure BDA0003200349740000086
——第a个空气状态参数n个样本模拟值的平均值。
S5、通过集合卡尔曼滤波修正集合元素中的条件参数变量的参数值并更新初始集合为E'；
E'＝E+K(YH-P) (11)
式中，H——将空气状态参数测量值投影至空气状态参数模拟值的投影矩阵，为1×n的元素均为1的矩阵，从而使空气状态参数测量值矩阵转化为b×n的矩阵；
则集合中的条件参数变量样本xi,j修正并更新如下：
Figure BDA0003200349740000091
Figure BDA0003200349740000092
S6、计算空气状态参数模拟值与空气状态参数测量值之间的均方根误差(RootMean Squared Error，RMSE)，设定收敛标准，判定计算收敛与否；
第a个空气状态参数模拟值与空气状态参数测量值的均方根误差计算如下：
Figure BDA0003200349740000093
设置收敛准则为：
RMSEa≤|Da|2 (15)
式中，Da——第a个空气状态参数的测量值的最大误差。
当所有空气状态参数的RMSE均达到收敛标准时，则判定计算收敛，执行步骤S7；否则，令E＝E'，返回步骤3，循环执行步骤3～6，直至迭代计算收敛。
S7、输出监测数据同化后博物馆展陈空间内气流组织的可视化结果。
对于建筑空间流场的模拟，传统方法是根据室内热湿负荷设计值(或估计值)及空调系统设计参数设置CFD计算模型边界条件，当计算结果与实际测量不符时，根据测量值与模拟值的差异，以试错的方式人工调节模型边界条件等参数。不同于上述方法，本发明将CFD模拟技术与集合卡尔曼滤波法相耦合，基于IEnKF原理，通过将监测系统的测量数据与CFD模拟结果相融合，并且不同于已有技术中直接通过卡尔曼滤波计算修正模拟结果，为了保证流场最终计算结果满足计算流体动力学控制方程的约束，该技术利用集合卡尔曼滤波法反向修正CFD计算模型的边界条件参数，并将CFD正向模拟与集合卡尔曼滤波反向修正边界条件的过程进行循环迭代，最终获得优化后接近实际工况的边界条件参数对应的流场模拟结果。该技术有效提高了模拟计算的精度，优化了博物馆展陈空间气流组织可视化的输出结果，同时，相比于传统方法，大大减少了人工试错的时间成本，具有灵活性强，实用性高等优点。针对对文物保存环境有极高要求的博物馆展陈空间，该技术通过有限数量监测点的数据，掌握温湿度空间场分布规律，更加准确地反映空间内的空气状态参数，实现对展陈空间内不同位置温湿度情况的监测，从而更好地指导展品的位置摆放和空调系统的运行调节，为博物馆藏品预防性保护和风险防控提供技术支撑。
实施例1：
如图2所示，根据典型空调房间进行简化，建立一个二维立面的理想空间，该空间立面尺寸为宽×高＝5×3m，其顶面为层间楼板，两侧为建筑内墙，底部为地面且设有热源，该房间采用上送下回的气流组织形式；其中，顶板上设有两个送风口，两送风口风速、风温协同变化(即风量成定比，风温一致)，两侧壁面下方分别设有一个回风口。此外，图2中设置了3个室内温度和风速的观测点，分别为观测点1、观测点2和观测点3。
经分析，针对该空间气流组织模拟设置的边界条件中共包含3个独立的条件参数变量，分别为送风风速、送风风温及热源温度，即m＝3，条件参数变量为xi，i＝1,2,3。已知该空间的设计条件和空调系统的历史运行数据，各边界条件参数值的分布规律如表1所示：
表1边界条件参数变量设计值分布规律
Figure BDA0003200349740000101
根据表1中边界条件的分布规律对各条件参数变量分别进行随机采样并构建初始集合E，取样本数量为30，则初始集合E由30个元素构成，即n＝30，那么初始集合E＝[Z1，Z2，Z3，···，Z30]，改写为矩阵形式为：
Figure BDA0003200349740000102
当前实际工况下该空间的空气状态参数测量数据如表2所示，在三个测点位置(观测点1、观测点2和观测点3)分别测量了两种空气状态参数(温度和风速)，共6个空气状态参数测量值，则b＝6，获得测量值向量Y：
Y＝[y1 y2 … y6]T (2)
表2空气状态参数测量数据
Figure BDA0003200349740000104
Figure BDA0003200349740000111
将上述30种边界条件工况分别代入CFD正演模型中进行计算求解，获得30种流场模拟结果，输出每种工况下对应3个测点位置的温度、风速模拟值，得到模拟值矩阵P如下：
Figure BDA0003200349740000112
计算条件参数变量的取样值与各空气状态参数模拟值之间的协方差矩阵Cxp、各空气状态参数模拟值之间的协方差矩阵Cpp以及空气状态参数测量值误差方差矩阵R，计算过程如下：
Figure BDA0003200349740000113
Figure BDA0003200349740000114
Figure BDA0003200349740000115
Figure BDA0003200349740000121
Figure BDA0003200349740000122
Figure BDA0003200349740000123
基于以上协方差矩阵的计算结果计算卡尔曼增益矩阵K，
K＝Cxp(Cpp+R)-1 (10)
依据卡尔曼滤波方程修正集合元素中的条件参数变量的参数值并更新集合为E’，
E'＝E+K(YH-P) (11)
则集合中的条件参数变量的样本值x'i,j修正并更新如下：
Figure BDA0003200349740000124
Figure BDA0003200349740000125
计算空气状态参数模拟值与测量值之间的方根误差RMSE，判断计算是否收敛。若未收敛，令E＝E'，返回步骤3，循环执行步骤3～6，直至迭代计算收敛，最终输出该空间优化的气流组织可视化结果。
各边界条件参数在数据同化过程中的优化效果如图3(a)-图3(d)所示，可见集合中各条件参数变量的30个样本在迭代计算过程中不断向真实值逼近，最终达到收敛。送风风速的均方根误差减小到0.02m/s，为初始误差的10％；送风温度的均方根误差减小到0.1℃，为原来的2.5％；热源温度均方根误差减小到0.4℃，为原来的8％；此外，对集合中30个样本条件参数变量及其对应的模拟结果分别取平均，优化前后的数据对比如表3所示，优化前后气流组织可视化结果对比如图4(a)-图4(d)所示，可见该技术有效提高了空间气流组织可视化结果的精度。
表3优化前后数据对比
Figure BDA0003200349740000131

Claims (10)

1.针对博物馆展陈空间气流组织可视化的监测数据同化方法，其特征在于：包括以下步骤：
S1、针对展陈空间内影响热湿环境的因素，确定模型边界条件中具有不确定性的条件参数变量，统计具有不确定性的条件参数变量的数量，并输入到数据同化算法程序中，同时向数据同化算法程序中输入各条件参数变量对应的条件参数值的区间范围；
根据每个条件参数变量在其变化区间内的一般分布规律，对各条件参数变量分别进行采样，并构建所有条件参数变量的初始集合E；
S2、向数据同化算法程序中输入展陈空间内无线监测系统所采集的观测数据，构建测量值向量；
S3、通过CFD正演模型分别对初始集合中的各元素进行求解，获得观测节点位置对应的空气状态参数模拟值矩阵；
S4、基于各条件参数变量的取样样本值与各空气状态参数模拟值之间的协方差矩阵、各空气状态参数模拟值之间的协方差矩阵以及测量误差方差矩阵，计算卡尔曼增益矩阵；
S5、通过集合卡尔曼滤波法修正集合元素中的条件参数变量的参数值并更新初始集合为E'；
S6、计算空气状态参数模拟值与空气状态参数测量值之间的均方根误差，设定收敛标准，判定计算收敛与否；当所有空气状态参数的RMSE均达到收敛标准时，则判定计算收敛，执行步骤S7；否则，令E＝E'，返回步骤S3，循环执行步骤S3～S6，直至迭代计算收敛；
S7、输出监测数据同化后博物馆展陈空间内气流组织的可视化结果。
2.根据权利要求1所述的针对博物馆展陈空间气流组织可视化的监测数据同化方法，其特征在于：S1中，初始集合E的计算过程如下：
设模型边界条件中存在m个具有不确定性的条件参数变量xi，i＝1,2,3，…，m；设每个条件参数变量的取样数量为n，即初始集合中元素的数量为n，初始集合中的第j个元素为Zj＝[x1,j，x2,j，x3,j，…，xm,j]，j＝1,2,3,…,n；构成初始集合E＝[Z1，Z2，Z3，···，Zn]，对应的矩阵形式为：
Figure FDA0003200349730000011
式中，xi,j——第i个条件参数变量的第j个取样样本值。
3.根据权利要求2所述的针对博物馆展陈空间气流组织可视化的监测数据同化方法，其特征在于：S2中，构建测量值向量的计算过程如下：
设数据同化算法程序中所构建的测量值向量为Y，假设在一次测量中，在c个测量点测量d种空气状态参数，则观测数据共包含b＝c×d个空气状态参数测量值ya，a＝1,2,3，…，b；则测量值向量Y为：
Y＝[y1 y2…yb]T (2)。
4.根据权利要求3所述的针对博物馆展陈空间气流组织可视化的监测数据同化方法，其特征在于：S3中，构建模拟值矩阵的具体计算过程如下：
将模型中的边界条件设定为集合E中每个元素的样本值，一组边界条件为一种工况，将集合E中的元素Zj分别代入CFD正演模型中，通过CFD模拟，对每个工况分别进行求解计算，获得各个工况对应的展陈空间流场分布模拟结果，构建模拟值矩阵P如下：
Figure FDA0003200349730000021
式中，pa,j——第j个集合元素对应的第a个空气状态参数模拟值。
5.根据权利要求4所述的针对博物馆展陈空间气流组织可视化的监测数据同化方法，其特征在于：S4中，卡尔曼增益矩阵的具体计算过程如下：
K＝Cxp(Cpp+R)-1 (4)
Figure FDA0003200349730000022
Figure FDA0003200349730000023
Figure FDA0003200349730000031
Figure FDA0003200349730000032
Figure FDA0003200349730000033
Figure FDA0003200349730000034
式中，Cxp——集合中的每个条件参数变量的取样样本值与各空气状态参数模拟值的协方差矩阵；Cxip——集合中的第i个条件参数变量的取样样本值与各空气状态参数模拟值的协方差矩阵；Cpp——各空气状态参数模拟值之间的协方差矩阵；R——测量误差方差矩阵，测量误差方差矩阵是由无线监测系统所采集的各测量值自身的背景噪声所组成的方差矩阵，其中σa——第a个空气状态参数测量值标准差，由测量仪器误差限ηa计算得，σa＝ηa/3；
Figure FDA0003200349730000035
——输入的第i个条件参数变量n个样本的平均值；
Figure FDA0003200349730000036
——第a个空气状态参数n个样本模拟值的平均值。
6.根据权利要求5所述的针对博物馆展陈空间气流组织可视化的监测数据同化方法，其特征在于：S5中，修正后的初始集合E’、修正后的初始集合元素Z'j和修正后的条件参数变量样本x'i,j的计算过程如下：
E'＝E+K(YH-P) (11)
式中，H——将空气状态参数测量值投影至空气状态参数模拟值的投影矩阵，为1×n的元素均为1的矩阵，从而使空气状态参数测量值矩阵转化为b×n的矩阵；
则集合中的条件参数变量样本xi,j修正并更新如下：
Figure FDA0003200349730000041
Figure FDA0003200349730000042
7.根据权利要求6所述的针对博物馆展陈空间气流组织可视化的监测数据同化方法，其特征在于：S6中，空气状态参数模拟值与空气状态参数测量值之间的均方根误差和收敛准则计算过程如下：
第a个空气状态参数模拟值与空气状态参数测量值的均方根误差计算如下：
Figure FDA0003200349730000043
设置收敛准则为：
RMSEa≤|Da|2 (15)
式中，Da——第a个空气状态参数的测量值的最大误差。
8.根据权利要求1-7任一项权利要求所述的针对博物馆展陈空间气流组织可视化的监测数据同化方法，其特征在于：S1中，判定模型边界条件中具有不确定性的条件参数变量的方法为：针对展陈空间内影响热湿环境的各个因素，分析其在展陈空间使用过程中或特定的研究时间段内的数值大小是否在一定区间范围内波动变化，若该因素波动变化，则判定其为具有不确定性的条件参数变量，否则判定其为常量。
9.根据权利要求8所述的针对博物馆展陈空间气流组织可视化的监测数据同化方法，其特征在于：S1中，所述的条件参数包括展陈空间内设备散热量、人体散热量、围护结构传热量、热源温度及空调运行参数。
10.根据权利要求9所述的针对博物馆展陈空间气流组织可视化的监测数据同化方法，其特征在于：S1中，所述的条件参数值的取值范围根据历史数据或设计值给出。
CN202110902814.4A 2021-08-06 2021-08-06 针对博物馆展陈空间气流组织可视化的监测数据同化方法 Active CN113591410B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110902814.4A CN113591410B (zh) 2021-08-06 2021-08-06 针对博物馆展陈空间气流组织可视化的监测数据同化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110902814.4A CN113591410B (zh) 2021-08-06 2021-08-06 针对博物馆展陈空间气流组织可视化的监测数据同化方法

Publications (2)

Publication Number Publication Date
CN113591410A true CN113591410A (zh) 2021-11-02
CN113591410B CN113591410B (zh) 2024-03-19

Family

ID=78255877

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110902814.4A Active CN113591410B (zh) 2021-08-06 2021-08-06 针对博物馆展陈空间气流组织可视化的监测数据同化方法

Country Status (1)

Country Link
CN (1) CN113591410B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115239769A (zh) * 2022-06-02 2022-10-25 中国国家博物馆 对象的分布位置确定方法、装置、电子设备和介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106646356A (zh) * 2016-11-23 2017-05-10 西安电子科技大学 一种基于卡尔曼滤波定位的非线性系统状态估计方法
CN107545100A (zh) * 2017-07-31 2018-01-05 江苏大学 一种venlo型温室环境的高分辨率降阶建模方法
CN110489835A (zh) * 2019-08-05 2019-11-22 安徽建筑大学 基于Ansys软件的自然通风与墙体蓄热耦合模拟方法
CN113158595A (zh) * 2021-04-23 2021-07-23 重庆中国三峡博物馆 一种博物馆文物展柜内部三维瞬态温度场模型的构建方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106646356A (zh) * 2016-11-23 2017-05-10 西安电子科技大学 一种基于卡尔曼滤波定位的非线性系统状态估计方法
CN107545100A (zh) * 2017-07-31 2018-01-05 江苏大学 一种venlo型温室环境的高分辨率降阶建模方法
CN110489835A (zh) * 2019-08-05 2019-11-22 安徽建筑大学 基于Ansys软件的自然通风与墙体蓄热耦合模拟方法
CN113158595A (zh) * 2021-04-23 2021-07-23 重庆中国三峡博物馆 一种博物馆文物展柜内部三维瞬态温度场模型的构建方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115239769A (zh) * 2022-06-02 2022-10-25 中国国家博物馆 对象的分布位置确定方法、装置、电子设备和介质

Also Published As

Publication number Publication date
CN113591410B (zh) 2024-03-19

Similar Documents

Publication Publication Date Title
Sangogboye et al. Performance comparison of occupancy count estimation and prediction with common versus dedicated sensors for building model predictive control
Jihad et al. Forecasting the heating and cooling load of residential buildings by using a learning algorithm “gradient descent”, Morocco
CN103049612B (zh) 一种基于模型降阶技术的建筑室内环境优化方法
JP6289749B2 (ja) 室内環境モデル作成装置
JP5674193B2 (ja) 空調環境モニタリングシステム
US20170300599A1 (en) System and method for calibrating multi-level building energy simulation
EP2941598A1 (en) On-line optimization scheme for hvac demand response
JP2017532688A (ja) 空調システム、並びに空調システムの動作を制御するためのシステム及び方法
WO2007122677A1 (ja) 熱対流場・物質拡散場の設計支援方法、設計支援システム、および設計支援用プログラム
Alam et al. Uncertainties in neural network model based on carbon dioxide concentration for occupancy estimation
JP6581490B2 (ja) 空調パラメータ生成装置、空調運用評価装置、空調パラメータ生成方法およびプログラム
JP6433588B2 (ja) 空調システム並びに空調システムの動作を制御するためのシステム及び方法
WO2015011446A1 (en) Control device and method for buildings
Chen et al. Fan-independent air balancing method based on computation model of air duct system
Wu et al. Multi-stage regression linear parametric models of room temperature in office buildings
CN112183826B (zh) 基于深度级联生成对抗网络的建筑能耗预测方法及相关产品
CN113591410A (zh) 针对博物馆展陈空间气流组织可视化的监测数据同化方法
Zekar et al. Development and assessment of simplified building representations under the context of an urban energy model: Application to arid climate environment
Arendt et al. Influence of input data on airflow network accuracy in residential buildings with natural wind-and stack-driven ventilation
Vidhyashankar et al. Modelling spatial variations in thermal comfort in indoor open-plan spaces using a whole-building simulation tool
Boukhris et al. Modeling coupled heat transfer and air flow in a partitioned building with a zonal model: Application to the winter thermal comfort
WO2021205512A1 (ja) 推定方法、シミュレーション方法、推定装置、及び推定プログラム
Abdo-Allah et al. Modeling, analysis, and state feedback control design of a multizone HVAC system
JP7484243B2 (ja) 栽培施設の設計支援システム、栽培施設空調システム、栽培施設の設計支援方法、および、栽培施設の設計支援プログラム
Andrade-Cabrera et al. A study on the trade-off between energy forecasting accuracy and computational complexity in lumped parameter building energy models

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
GR01 Patent grant