CN116298379A - 测量室外风速与辐射的方法、装置、设备及存储介质 - Google Patents
测量室外风速与辐射的方法、装置、设备及存储介质 Download PDFInfo
- Publication number
- CN116298379A CN116298379A CN202310580503.XA CN202310580503A CN116298379A CN 116298379 A CN116298379 A CN 116298379A CN 202310580503 A CN202310580503 A CN 202310580503A CN 116298379 A CN116298379 A CN 116298379A
- Authority
- CN
- China
- Prior art keywords
- wind speed
- thermometer
- temperature
- ball
- radiation
- 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
Links
- 230000005855 radiation Effects 0.000 title claims abstract description 89
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000003860 storage Methods 0.000 title claims abstract description 11
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 19
- 238000001914 filtration Methods 0.000 claims abstract description 12
- 230000008569 process Effects 0.000 claims description 15
- 238000009529 body temperature measurement Methods 0.000 claims description 14
- 238000005259 measurement Methods 0.000 claims description 11
- 238000004590 computer program Methods 0.000 claims description 10
- 238000002474 experimental method Methods 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 9
- 230000008878 coupling Effects 0.000 claims description 7
- 238000010168 coupling process Methods 0.000 claims description 7
- 238000005859 coupling reaction Methods 0.000 claims description 7
- 230000004907 flux Effects 0.000 claims description 7
- 238000001514 detection method Methods 0.000 claims description 6
- 239000012774 insulation material Substances 0.000 claims description 6
- 238000004381 surface treatment Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 238000010219 correlation analysis Methods 0.000 claims description 4
- 239000003973 paint Substances 0.000 claims description 4
- 239000000523 sample Substances 0.000 claims description 4
- 229910052782 aluminium Inorganic materials 0.000 claims description 3
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims description 3
- 229910052751 metal Inorganic materials 0.000 claims description 3
- 239000002184 metal Substances 0.000 claims description 3
- 238000001228 spectrum Methods 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 abstract description 4
- 238000004364 calculation method Methods 0.000 description 13
- 230000007613 environmental effect Effects 0.000 description 6
- 230000004044 response Effects 0.000 description 6
- 238000012546 transfer Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 239000006185 dispersion Substances 0.000 description 3
- 230000001052 transient effect Effects 0.000 description 3
- 238000009826 distribution Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- VYZAMTAEIAYCRO-UHFFFAOYSA-N Chromium Chemical compound [Cr] VYZAMTAEIAYCRO-UHFFFAOYSA-N 0.000 description 1
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052804 chromium Inorganic materials 0.000 description 1
- 239000011651 chromium Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000005485 electric heating Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 239000011810 insulating material Substances 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000007747 plating Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005507 spraying Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
- G01P5/10—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring thermal variables
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/10—Radiation pyrometry, e.g. infrared or optical thermometry using electric radiation detectors
- G01J5/12—Radiation pyrometry, e.g. infrared or optical thermometry using electric radiation detectors using thermoelectric elements, e.g. thermocouples
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J5/00—Radiation pyrometry, e.g. infrared or optical thermometry
- G01J5/80—Calibration
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Abstract
本发明公开了测量室外风速与辐射的方法、装置、设备及存储介质,属于室外热环境监测领域,包括以下步骤:标定频域阈值‑计算风速和辐射:利用双色球形温度计组中的黑球温度计和亮面球温度计分别采集黑球全日温度时间序列、亮面球全日温度时间序列,同时借助气温测量设备采集室外空气全日温度时间序列‑以频域阈值为滤波阈值,以滤波方法对各温度时间序列进行塔式信号分解,进而计算风速‑根据高阶热平衡理论,带入风速中频分量,计算平均辐射温度。本发明采用上述测量室外风速与辐射的方法、装置、设备及存储介质,仅需测量温度这一种形式的信号,使黑球温度计的应用摆脱对独立风速计的依赖。
Description
技术领域
本发明涉及室外热环境监测技术领域,尤其涉及测量室外风速与辐射的方法、装置、设备及存储介质。
背景技术
平均辐射温度(Mean Radiant Temperature,MRT)是评估室外热舒适与热安全的关键环境参数。其为一个理论假设参数,物理定义是:与实际环境对人体产生的总辐射(包括长波辐射与短波辐射)传热作用相等的假想黑体环境的平均表面温度。上述定义可用以下公式表示。
式中,是斯蒂芬玻尔兹曼常数,/>;/>是平均辐射温度,单位为/>;/>是平均短波辐射通量,单位为/>;L是平均长波辐射通量,单位为/>。、/>、/>分别为被研究对象的短波发射率、长波发射率和环境发射率。
一方面,黑球温度计是测量日间室外环境平均辐射温度的常用工具,该方法用黑色球体模拟人体,基于黑球系统瞬态热平衡假设计算平均辐射温度。瞬态热平衡假设是指认为每一个瞬间球壳外表面的对流换热与辐射传热之和为零,其表达式如下:
式中,是球体外表面强迫对流换热系数,其根据Hey, E在《Small GlobeThermometers. Journal of Physics E: Scientific Instruments》给出的球体对流换热经验公式/>计算可得;/>是黑球空腔气温,单位为/>;/>是室外空气温度,单位为/>;
上述三个参数属于温度与风速两种不同形式的信号,除温度记录装置外还需单独的风速计。平均辐射温度的测量过程在仪器布置与信号采集方面都颇为不便。
另一方面,现有技术条件下通过黑球温度计测算出的室外平均辐射温度离散误差非常明显(黑球温度计方法将剧烈波动的室外风速与相对平稳的辐射强行带入热平衡公式b是造成平均辐射温度计算结果离散误差过大的根本原因)。这是由于室外风速与辐射的主导周期明显不同,室外风速是一种相对高频信号,而辐射是一种相对低频信号。受制造工艺限制黑球温度计的系统热容较大,且室外风速波动非常剧烈,瞬态热平衡假设在室外条件下根本不能成立,真实的黑球表面换热情况应为非稳态过程,其表达式如下:
可知,黑球温度计对室外环境参数的热响应是一个非常复杂的过程,黑球空腔气温()作为风速与辐射的因变量携带丰富的隐含信息。比如通过简单的观测可以发现,环境辐射越强黑球温度信号的局部震荡越剧烈、波动幅度越大。这是因为白天太阳光照射在黑球温度计上使其温度升高,模拟热球风速仪的电加热功能,黑球温度计对风速的波动更加敏感。由此可见/>本身携带了风速的信息,但是其表现被辐射调节。
综上可知,现有利用黑球温度计测量平均辐射温度的技术对黑球温度的时间序列缺乏具体的、分层次的认识,具有以下缺陷:
1、忽略了黑球温度时间序列中表现出的风速信息,尤其是缺乏对其频域特性的分解与利用;
2、需依赖单独的风速计,测量过程复杂;
3、测量结果精度低,离散误差显著。
发明内容
为解决上述问题,本发明提供一种测量室外风速与辐射的方法、装置、设备及存储介质,仅需测量温度这一种形式的信号,使黑球温度计的应用摆脱对独立风速计的依赖,从而可直接用于定点室外风速、辐射测量;且由于其硬件系统简单经济,也适用于大范围、多点位的室外热环境评估与监控。
为实现上述目的,本发明提供了测量室外风速与辐射的方法,包括以下步骤:
S2、计算风速和辐射:
S21、利用双色球形温度计组中的黑球温度计和亮面球温度计分别采集黑球全日温度时间序列、亮面球全日温度时间序列,同时借助气温测量设备采集室外空气全日温度时间序列,采样频率均为0.1Hz;
S23、根据高阶热平衡理论,带入风速中频分量,计算平均辐射温度。
优选的,步骤S1具体包括以下步骤:
S111、将双色球形温度计组移至风洞中,利用风洞模拟室外风速与辐射条件;
S112、记录风速与黑球温度;
S113、对时间序列做不同窗宽的滑动平均后下采样处理,生成频率与黑球温度相同的时间序列;
S114、计算黑球温度与时间序列的Spearman相关系数;
S122、对黑球温度和风速时间序列使用Morlet小波做交叉小波能量谱,以红噪声检验标定95%置信区间所在频域范围为二者的共有周期;
S123、对黑球温度和风速时间序列做小波相干,标定95%置信区间所在频域范围为二者的趋势互锁频域;
S124、将共有周期与趋势互锁频域中的共同部分标记出来,此频域范围为黑球温度和风速时间序列的耦合频域;
优选的,步骤S22具体包括以下步骤:
S223、计算平均风速:
定义:
式中,脚标、/>分别表示参数经塔式分解后以/>为阈的中低频分量和次高频分量;脚标/>表示参数经塔式分解后以/>为阈的低频分量;/>、/>为经风洞实验拟合的中频风速系数;/>、/>为经风洞实验拟合的次高频风速系数;
则有:
优选的,步骤S23具体包括以下步骤:
测量室外风速与辐射的方法的装置,包括双色球形温度计组和气温测量设备,双色球形温度计组包括黑球温度计和亮面球温度计;
黑球温度计和亮面球温度计均包括由上到下依次设置的球壳、隔热材料和支撑杆,球壳内部设置有测温探点,测温探点与测温元件的一端连接,测温元件的另一端依次穿过隔热材料和支撑杆后连接温度记录仪;
其需满足以下条件:
条件一:黑球温度计和亮面球温度计的结构、尺寸相同、系统热容均相等;
条件二:黑球温度计与亮面球温度计的表面发射率满足以下公式:
条件三:黑球温度计的球壳表面和亮面球温度计的球壳表面分别施哑光黑漆、亮面表面处理工艺,且施加表面处理工艺后的黑球温度计与亮面球温度计的表面发射率满足以下公式:
条件四:黑球温度计和亮面球温度计在每个时刻的对流换热系数均相同,且仅与该时刻的空气流速有关;
条件五:黑球温度计的球壳的球心和亮面球温度计的球壳的球心位于同一水平面上,且间距范围为5~10D,D为球壳的直径,二者角系数小于1%;
条件六:双色球形温度计组和气温测量设备的采样频率均为0.1Hz,精度均不低于±0.5℃。
优选的,测温元件为k型热电偶,测温探点位于球壳的球心;
球壳由导热系数大于铝的金属制成;球壳的球壁厚度小于1mm;
当双色球形温度计组仅用于测量室外平均风速时,球壳直径D≤0.02m;
当双色球形温度计组用于测量室外平均风速和室外平均辐射温度时,球壳直径为0.04m。
一种计算机设备,该计算机设备包括:至少一个处理器、存储器以及存储在存储器中并可在至少一个处理器上运行的计算机程序,处理器执行计算机程序时实现上述测量室外风速与辐射的方法。
一种计算机可读存储介质,计算机可读存储介质存储有计算机程序,计算机程序被处理器执行时实现上述测量室外风速与辐射的方法。
本发明具有以下有益效果:
1、使用的硬件装置无加热、转动或其他控制元件,使得其在室外应用的稳定性更强、能耗与成本更低。
2、计算结果中的风速中频分量可作为平均辐射温度计算式的输入参数,使MRT计算结果离散误差大幅降低,精度提高。
3、仅需测量温度这一种形式的信号,使黑球温度计的应用摆脱对独立风速计的依赖。
4、可同时测算室外热环境的三个关键参数(风速、气温、平均辐射温度),有潜力发展为室外热环境监测领域亟需的专业测量工具。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
图1为本发明的测量室外风速与辐射的方法的流程图。
图2为本发明的测量室外风速与辐射的方法的傅里叶周期分布图;
图3为本发明的测量室外风速与辐射的方法的塔式信号分解示意图;
图4为本发明的双色球形温度计组的结构示意图;
图5为本发明的实施例的黑球温度信号分解结果与风速计算结果的比较图;
图6为本发明的实施例的2分钟平均风速计算值与超声风速计测量参考值的比较图;
图7为本发明的实施例的计算误差的概率密度统计结果图。
其中:1、球壳;2、测温探点;3、隔热材料;4、支撑杆;5、测温元件。
具体实施方式
为了使本发明实施例公开的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明实施例进行进一步详细说明。应当理解,此处描述的具体实施例仅仅用以解释本发明实施例,并不用于限定本发明实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本申请保护的范围。所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。
需要说明的是,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或服务器不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
相似的标号和字母在下面的附图中表示类似项,因此,一旦某一项在一个附图中被定义,则在随后的附图中不需要对其进行进一步定义和解释。
在本发明的描述中,需要说明的是,术语“上”、“下”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,或者是该发明产品使用时惯常摆放的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
在本发明的描述中,还需要说明的是,除非另有明确的规定和限定,术语“设置”、“安装”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
图1为本发明的测量室外风速与辐射的方法的流程图,如图1所示,测量室外风速与辐射的方法,包括以下步骤:
图2为本发明的测量室外风速与辐射的方法的傅里叶周期分布图,由图2可知,根据信号分量的傅里叶周期将整个频域空间分成以下四个频段:
优选的,步骤S1中所述的、/>、/>的具体数值由黑球温度计的系统热容决定(黑球温度计的系统热容与球壳的材料、尺寸和工艺有关,系统热容越大热惯性越大,能响应的高频分量越少,响应区间朝频域的低频方向移动),其与时间常数有关。
式中,是平均短波辐射通量;/>是平均长波辐射通量;脚标/>表示各参数经塔式分解后以/>为阈的中低频分量;/>、/>、/>分别为被研究对象的短波发射率、长波发射率和环境发射率;/>为球体高阶热平衡的对流换热系数,其需根据各环境参数以/>为阈的高频分量重新拟合。
经以上分析可知,黑球温度信号中已经包含了风速信息。根据雷诺分解原理,风速是高频分量、次高频分量与中频分量的叠加。黑球温度计受热惯性制约仅能响应其中的中频分量与次高频分量,且具有不同的热响应规律。其中,风速中频分量满足高阶热平衡关系,作为平均辐射温度的输入参数,可克服计算结果离散过大的问题;次高低频分量与中频分量叠加后其本质为时间窗宽的室外风速平均值,而日间平均风速恰恰是对室外人体热舒适与热安全产生直接影响的关键环境参数之一,具有重要应用价值。
优选的,步骤S1具体包括以下步骤:
S111、将双色球形温度计组移至风洞中,利用风洞模拟室外风速与辐射条件;
S112、记录风速与黑球温度;
S113、对时间序列做不同窗宽的滑动平均后下采样处理,生成频率与黑球温度相同的时间序列;
S114、计算黑球温度与时间序列的Spearman相关系数;
S122、对黑球温度和风速时间序列使用Morlet小波做交叉小波能量谱,以红噪声检验标定95%置信区间所在频域范围为二者的共有周期;
S123、对黑球温度和风速时间序列做小波相干,标定95%置信区间所在频域范围为二者的趋势互锁频域;
S124、将共有周期与趋势互锁频域中的共同部分标记出来,此频域范围为黑球温度和风速时间序列的耦合频域;
S2、计算风速和辐射:
S21、利用双色球形温度计组中的黑球温度计和亮面球温度计分别采集黑球全日温度时间序列、亮面球全日温度时间序列,同时借助气温测量设备采集室外空气全日温度时间序列,采样频率均为0.1Hz;
图3为本发明的测量室外风速与辐射的方法的塔式信号分解示意图,如图3所示,S22、以、/>、/>为滤波阈值,将采集的黑球全日温度时间序列、亮面球全日温度时间序列以及室外空气全日温度时间序列,以滤波方法进行塔式信号分解,生成若干个相对低频分量与相对高频分量(此处的相对低频分量为approximation(信号相似),相对高频分量为detail(信号细节),且下一级的approximation和detail的和等于上一级信号),进而计算风速;
优选的,步骤S22具体包括以下步骤:
S223、计算平均风速:
定义:
则有:
S23、根据高阶热平衡理论,带入风速中频分量,计算平均辐射温度。
优选的,步骤S23具体包括以下步骤:
图4为本发明的双色球形温度计组的结构示意图,需要说明的是,图4中左侧视图为亮面球温度计外形图,右侧图为黑球温度计剖视图,如图4所示,测量室外风速与辐射的方法的装置,包括双色球形温度计组和气温测量设备,双色球形温度计组包括黑球温度计和亮面球温度计;
黑球温度计和亮面球温度计均包括由上到下依次设置的球壳、隔热材料和支撑杆,球壳1内部设置有测温探点2,测温探点2与测温元件5的一端连接,测温元件5的另一端依次穿过隔热材料3和支撑杆4后连接温度记录仪,温度记录仪为至少包含3个输入通道的数据记录仪;
其需满足以下条件:
条件一:黑球温度计和亮面球温度计的结构、尺寸相同、系统热容均相等;
条件二:黑球温度计与亮面球温度计的表面发射率满足以下公式:
条件三:黑球温度计的球壳1表面和亮面球温度计的球壳1表面分别施哑光黑漆、亮面表面处理工艺,且施加表面处理工艺后的黑球温度计与亮面球温度计的表面发射率满足以下公式:
条件四:黑球温度计和亮面球温度计在每个时刻的对流换热系数均相同,且仅与该时刻的空气流速有关;
条件五:黑球温度计的球壳1的球心和亮面球温度计的球壳1的球心位于同一水平面上,且间距范围为5~10D,D为球壳的直径,二者角系数小于1%;
条件六:双色球形温度计组和气温测量设备的采样频率均为0.1Hz,精度均不低于±0.5℃。
优选的,测温元件5为k型热电偶,测温探点2位于球壳的球心;
球壳1由导热系数大于铝的金属制成;球壳1的球壁厚度小于1mm;
当双色球形温度计组仅用于测量室外平均风速时,球壳1直径D≤0.02m;
当双色球形温度计组用于测量室外平均风速和室外平均辐射温度时,球壳1直径为0.04m。
优选的,还包括与黑球温度计或者亮面球温度计的尺寸、结构、系统热容相同的白球温度计,用于根据计算出的风速中频分量作为输入量,进而测算室外环境的长波辐射通量L和短波辐射通量S。计算公式记载于在先申请202310159425.6所述的测量室外辐射的方法、配套硬件装置、设备及存储介质中,故在此不做赘述。
一种计算机设备,该计算机设备包括:至少一个处理器、存储器以及存储在存储器中并可在至少一个处理器上运行的计算机程序,处理器执行计算机程序时实现上述测量室外风速与辐射的方法。
一种计算机可读存储介质,计算机可读存储介质存储有计算机程序,计算机程序被处理器执行时实现上述测量室外风速与辐射的方法。
实施例
本实施例使用的双色球形温度计组物理结构为:铜制,外直径为2cm,壁厚0.7mm;系统热容,为方便计算此热容为平均到球体的单位外表面积的热容;时间常数/>约为85s;测温元件为k型热电偶。黑球表面喷涂哑光黑漆,分光光度计测得长波发射率/>、短波发射率/>;亮面球温度计采用表面镀铬工艺,使其呈镜面效果,分光光度计测得长波发射率/>、短波发射率/>,即/>。
上述双色球形温度计组硬件成本约为40元人民币。
在含辐射、风速双变量的风洞内完成以下频域阈值标定工作:
根据相关性分析标定:对风速采用滑动平均下采样处理后与黑球温度TB计算Spearman相关系数,当滑动平均窗宽度为114s时存在最小Spearman相关系数(-0.87),故该黑球温度计所能响应的最短傅里叶周期/>。114秒与2分钟相近,此双色球形温度计组可用于测量室外2分钟平均风速。
本例采用离散小波分解进行滤波,使用小波Symlet16,第6层对应的傅里叶周期为15.75min(),第9层对应的傅里叶周期为126min(/>)。次高频分量(/>)取前6层信号细节之和,中低频分量(/>)取第6层信号近似,低频分量(/>)取第9层信号近似。
经风洞实验拟合该黑球温度计风速计算式:
风速中频分量:
风速次高频分量:
以8月某日某地的室外2分钟平均风速,9:00-17:00的风速计算值为例,以超声风速计测量结果作为参考值,计算结果如图5所示(需要说明的是,图5中左列为黑球温度信号分解的中低频和次高频部分,右列为据此计算出的风速信号中频和次高频部分),同时结合图6和图7,可知本实施例所述方法在室外平均风速测量方面的精度满足室外热环境研究领域的观测应用要求。
因此,本发明采用上述测量室外风速与辐射的方法、装置、设备及存储介质,具有灵活性强、拓展度高、结构简单等优点,同时其通过测量温度(气温和球温),解析黑球温度信号提取风速,计算室外平均风速,使黑球温度计应用摆脱对独立风速计的依赖,并提高平均辐射温度测算精度。
最后应说明的是:以上实施例仅用以说明本发明的技术方案而非对其进行限制,尽管参照较佳实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对本发明的技术方案进行修改或者等同替换,而这些修改或者等同替换亦不能使修改后的技术方案脱离本发明技术方案的精神和范围。
Claims (10)
1.测量室外风速与辐射的方法,其特征在于:包括以下步骤:
S2、计算风速和辐射:
S21、利用双色球形温度计组中的黑球温度计和亮面球温度计分别采集黑球全日温度时间序列、亮面球全日温度时间序列,同时借助气温测量设备采集室外空气全日温度时间序列,采样频率均为0.1Hz;
S23、根据高阶热平衡理论,带入风速中频分量,计算平均辐射温度。
3.根据权利要求1所述的测量室外风速与辐射的方法,其特征在于:步骤S1具体包括以下步骤:
S111、将双色球形温度计组移至风洞中,利用风洞模拟室外风速与辐射条件;
S112、记录风速与黑球温度;
S113、对时间序列做不同窗宽的滑动平均后下采样处理,生成频率与黑球温度相同的时间序列;
S114、计算黑球温度与时间序列的Spearman相关系数;
S122、对黑球温度和风速时间序列使用Morlet小波做交叉小波能量谱,以红噪声检验标定95%置信区间所在频域范围为二者的共有周期;
S123、对黑球温度和风速时间序列做小波相干,标定95%置信区间所在频域范围为二者的趋势互锁频域;
S124、将共有周期与趋势互锁频域中的共同部分标记出来,此频域范围为黑球温度和风速时间序列的耦合频域;
4.根据权利要求1所述的测量室外风速与辐射的方法,其特征在于:步骤S22具体包括以下步骤:
S223、计算平均风速:
定义:
式中,脚标、/>分别表示参数经塔式分解后以/>为阈的中低频分量和次高频分量;脚标/>表示参数经塔式分解后以/>为阈的低频分量;/>、/>为经风洞实验拟合的中频风速系数;/>、/>为经风洞实验拟合的次高频风速系数;
则有:
6.一种如上述权利要求1-5任一项所述的测量室外风速与辐射的方法的装置,包括双色球形温度计组和气温测量设备,其特征在于:双色球形温度计组包括黑球温度计和亮面球温度计;
黑球温度计和亮面球温度计均包括由上到下依次设置的球壳、隔热材料和支撑杆,球壳内部设置有测温探点,测温探点与测温元件的一端连接,测温元件的另一端依次穿过隔热材料和支撑杆后连接温度记录仪;
其需满足以下条件:
条件一:黑球温度计和亮面球温度计的结构、尺寸相同、系统热容均相等;
条件二:黑球温度计与亮面球温度计的表面发射率满足以下公式:
条件三:黑球温度计的球壳表面和亮面球温度计的球壳表面分别施哑光黑漆、亮面表面处理工艺,且施加表面处理工艺后的黑球温度计与亮面球温度计的表面发射率满足以下公式:
条件四:黑球温度计和亮面球温度计在每个时刻的对流换热系数均相同,且仅与该时刻的空气流速有关;
条件五:黑球温度计的球壳的球心和亮面球温度计的球壳的球心位于同一水平面上,且间距范围为5~10D,D为球壳的直径,二者角系数小于1%;
条件六:双色球形温度计组和气温测量设备的采样频率均为0.1Hz,精度均不低于±0.5℃。
7.根据权利要求6所述的装置,其特征在于:测温探点位于球壳的球心;
球壳由导热系数大于铝的金属制成;球壳的球壁厚度小于1mm;
当双色球形温度计组仅用于测量室外平均风速时,球壳直径D≤0.02m;
当双色球形温度计组用于测量室外平均风速和室外平均辐射温度时,球壳直径为0.04m。
9.一种计算机设备,该计算机设备包括:至少一个处理器、存储器以及存储在存储器中并可在至少一个处理器上运行的计算机程序,其特征在于:处理器执行计算机程序时实现上述权利要求1-5任意一项的测量室外风速与辐射的方法。
10.一种计算机可读存储介质,计算机可读存储介质存储有计算机程序,其特征在于:计算机程序被处理器执行时实现上述权利要求1-5任意一项的测量室外风速与辐射的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310580503.XA CN116298379B (zh) | 2023-05-23 | 2023-05-23 | 测量室外风速与辐射的方法、装置、设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310580503.XA CN116298379B (zh) | 2023-05-23 | 2023-05-23 | 测量室外风速与辐射的方法、装置、设备及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116298379A true CN116298379A (zh) | 2023-06-23 |
CN116298379B CN116298379B (zh) | 2023-08-04 |
Family
ID=86815324
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310580503.XA Active CN116298379B (zh) | 2023-05-23 | 2023-05-23 | 测量室外风速与辐射的方法、装置、设备及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116298379B (zh) |
Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4592661A (en) * | 1983-11-18 | 1986-06-03 | Dobbie Instruments (Australia) Pty. Ltd. | Heat stress monitor |
JPS62288437A (ja) * | 1986-06-05 | 1987-12-15 | Mitsubishi Heavy Ind Ltd | 空気調和機の制御装置 |
JP2003232676A (ja) * | 2002-02-06 | 2003-08-22 | Takeo Saito | 3次元cpc型放射計 |
JP2004317015A (ja) * | 2003-04-16 | 2004-11-11 | Yamatake Corp | グローブ温度検出器及びこれを用いたグローブ温度検出方法並びにこれを利用するコンピュータ |
JP2010164550A (ja) * | 2008-12-19 | 2010-07-29 | Panasonic Electric Works Co Ltd | 赤外線放射素子及び当該赤外線放射素子を備えた赤外線式ガス検知器及び当該赤外線放射素子の製造方法 |
JP2013228374A (ja) * | 2012-03-27 | 2013-11-07 | Institute Of National Colleges Of Technology Japan | 日射状況予測システム,日射状況予測装置及び日射状況予測方法 |
CN104048763A (zh) * | 2013-03-13 | 2014-09-17 | 凯尔西-海耶斯公司 | 叶形开口辐射传感器 |
CN105371897A (zh) * | 2015-12-04 | 2016-03-02 | 深圳市建筑科学研究院股份有限公司 | 室外热舒适度监测系统及其监测方法 |
DE102015216780A1 (de) * | 2014-09-29 | 2016-03-31 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung zur Bestimmung eines Streustrahlungsbeitrags für eine Streustrahlungskorrektur eines Röntgenbildes |
KR20190109368A (ko) * | 2019-09-18 | 2019-09-25 | 창원대학교 산학협력단 | 미기상 정보수집용 종합환경계측기 |
JP2020056598A (ja) * | 2018-09-28 | 2020-04-09 | 三機工業株式会社 | 熱中症指標測定装置 |
CN211148956U (zh) * | 2019-10-30 | 2020-07-31 | 长春气象仪器研究所有限责任公司 | 基于黑球补偿的探空仪温度测量系统 |
JP2021096084A (ja) * | 2019-12-13 | 2021-06-24 | 学校法人東京理科大学 | 測定機器 |
CN114600057A (zh) * | 2019-07-05 | 2022-06-07 | 法国圣-戈班玻璃公司 | 用以生成车辆的热舒适度图的无线系统和方法 |
CN114811896A (zh) * | 2022-05-11 | 2022-07-29 | 中国电建集团城市规划设计研究院有限公司 | 一种工作区域空调系统的控制方法和装置 |
US20220365245A1 (en) * | 2019-10-01 | 2022-11-17 | Kyoto Electronics Manufacturing Co., Ltd. | Output method for electronic wbgt meter, and electronic wbgt meter |
CN115839769A (zh) * | 2023-02-24 | 2023-03-24 | 烟台大学 | 测量室外辐射的方法、配套硬件装置、设备及存储介质 |
-
2023
- 2023-05-23 CN CN202310580503.XA patent/CN116298379B/zh active Active
Patent Citations (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4592661A (en) * | 1983-11-18 | 1986-06-03 | Dobbie Instruments (Australia) Pty. Ltd. | Heat stress monitor |
JPS62288437A (ja) * | 1986-06-05 | 1987-12-15 | Mitsubishi Heavy Ind Ltd | 空気調和機の制御装置 |
JP2003232676A (ja) * | 2002-02-06 | 2003-08-22 | Takeo Saito | 3次元cpc型放射計 |
JP2004317015A (ja) * | 2003-04-16 | 2004-11-11 | Yamatake Corp | グローブ温度検出器及びこれを用いたグローブ温度検出方法並びにこれを利用するコンピュータ |
JP2010164550A (ja) * | 2008-12-19 | 2010-07-29 | Panasonic Electric Works Co Ltd | 赤外線放射素子及び当該赤外線放射素子を備えた赤外線式ガス検知器及び当該赤外線放射素子の製造方法 |
JP2013228374A (ja) * | 2012-03-27 | 2013-11-07 | Institute Of National Colleges Of Technology Japan | 日射状況予測システム,日射状況予測装置及び日射状況予測方法 |
CN104048763A (zh) * | 2013-03-13 | 2014-09-17 | 凯尔西-海耶斯公司 | 叶形开口辐射传感器 |
DE102015216780A1 (de) * | 2014-09-29 | 2016-03-31 | Siemens Aktiengesellschaft | Verfahren und Vorrichtung zur Bestimmung eines Streustrahlungsbeitrags für eine Streustrahlungskorrektur eines Röntgenbildes |
CN105371897A (zh) * | 2015-12-04 | 2016-03-02 | 深圳市建筑科学研究院股份有限公司 | 室外热舒适度监测系统及其监测方法 |
JP2020056598A (ja) * | 2018-09-28 | 2020-04-09 | 三機工業株式会社 | 熱中症指標測定装置 |
CN114600057A (zh) * | 2019-07-05 | 2022-06-07 | 法国圣-戈班玻璃公司 | 用以生成车辆的热舒适度图的无线系统和方法 |
KR20190109368A (ko) * | 2019-09-18 | 2019-09-25 | 창원대학교 산학협력단 | 미기상 정보수집용 종합환경계측기 |
US20220365245A1 (en) * | 2019-10-01 | 2022-11-17 | Kyoto Electronics Manufacturing Co., Ltd. | Output method for electronic wbgt meter, and electronic wbgt meter |
CN211148956U (zh) * | 2019-10-30 | 2020-07-31 | 长春气象仪器研究所有限责任公司 | 基于黑球补偿的探空仪温度测量系统 |
JP2021096084A (ja) * | 2019-12-13 | 2021-06-24 | 学校法人東京理科大学 | 測定機器 |
CN114811896A (zh) * | 2022-05-11 | 2022-07-29 | 中国电建集团城市规划设计研究院有限公司 | 一种工作区域空调系统的控制方法和装置 |
CN115839769A (zh) * | 2023-02-24 | 2023-03-24 | 烟台大学 | 测量室外辐射的方法、配套硬件装置、设备及存储介质 |
Non-Patent Citations (4)
Title |
---|
TAIRUI ZHANG: "Determination of the proof strength and flow properties of materials from spherical indentation tests: An analytical approach based on the expanding cavity model", JOURNAL OF STRAIN ANALYSIS * |
刘国丹: "基于太阳辐射影响的室内热环境特性", 暖通空调 * |
唐培傑: "夏热冬暖地区室外平均辐射温度估算修正", 哈尔滨工业大学学报 * |
张磊: "室外热环境评价指标湿球黑球温度简化计算方法", 重庆建筑大学学报 * |
Also Published As
Publication number | Publication date |
---|---|
CN116298379B (zh) | 2023-08-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN207662532U (zh) | 高精度温度检测装置 | |
CN111413364A (zh) | 一种建筑墙体内混凝土蓄热系数原位无损检测方法及系统 | |
CN115452180B (zh) | 一种高焓气流恢复温度测量方法及测量装置 | |
CN103983365B (zh) | 多测头瞬态辐射热流计及热辐射热流密度的测定方法 | |
Holdaway et al. | Inclusion of linearized moist physics in NASA’s Goddard Earth Observing System data assimilation tools | |
CN112595420A (zh) | 一种红外体温筛选仪及校正方法 | |
CN104792423A (zh) | 一种可校准测温参数的测温方法和装置 | |
CN116298379B (zh) | 测量室外风速与辐射的方法、装置、设备及存储介质 | |
CN115839769B (zh) | 测量室外辐射的方法、配套硬件装置、设备及存储介质 | |
Pal et al. | An acoustic hotspot tracking algorithm for highly centralized gas temperature distribution | |
CN111307481B (zh) | 一种基于非线性人工神经网络的动态热壁热流逆估计方法 | |
CN102621180B (zh) | 一种节能门窗性能测试方法 | |
CN203677065U (zh) | 一种红外体温监测自动校准系统 | |
Pan et al. | Intelligent compensation method of infrared temperature measurement for multiple interference factors | |
CN100439900C (zh) | 一种材料黑度系数的测试方法 | |
Liu et al. | Analysis of digital filters used in time-series small heat flux measurement | |
Yang et al. | Development of an air temperature measuring instrument used in the field of meteorological detection | |
Loomans et al. | Simulation and measurement of the stationary and transient characteristics of the hot sphere anemometer | |
Yang et al. | Multiphysical field analysis of a temperature sensor for meteorological observation | |
CN213932844U (zh) | 一种红外体温筛选仪 | |
Cárdenas-García et al. | Blackbody for metrological control of ear thermometers | |
CN204501718U (zh) | 一种可缩短测温时间的测温装置和智能奶瓶套 | |
OLESEN et al. | Measurements of the physical parameters of the thermal environment | |
Kaiho et al. | Neural network estimation of eardrum temperature using multiple sensors integrated on a wristwatch-sized device | |
CN204479185U (zh) | 一种测温装置 |
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 |