CN107917768A - 一种基于低频声波的空气温度测量装置及方法 - Google Patents

一种基于低频声波的空气温度测量装置及方法 Download PDF

Info

Publication number
CN107917768A
CN107917768A CN201710971935.8A CN201710971935A CN107917768A CN 107917768 A CN107917768 A CN 107917768A CN 201710971935 A CN201710971935 A CN 201710971935A CN 107917768 A CN107917768 A CN 107917768A
Authority
CN
China
Prior art keywords
air themperature
module
sound
sound wave
velocity
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
CN201710971935.8A
Other languages
English (en)
Other versions
CN107917768B (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.)
North China Electric Power University
Original Assignee
North China Electric Power University
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 North China Electric Power University filed Critical North China Electric Power University
Priority to CN201710971935.8A priority Critical patent/CN107917768B/zh
Publication of CN107917768A publication Critical patent/CN107917768A/zh
Application granted granted Critical
Publication of CN107917768B publication Critical patent/CN107917768B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K11/00—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
    • G01K11/22—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using measurement of acoustic effects
    • G01K11/24—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using measurement of acoustic effects of the velocity of propagation of sound

Landscapes

  • Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
  • Measuring Temperature Or Quantity Of Heat (AREA)

Abstract

本发明公开了属于空气温度测量技术范围的一种基于低频声波的空气温度测量装置及方法。所述温度测量装置具有由扬声器、一对传声器、湿度计和固定装置组成的传感器模块；传声器接收扬声器发出的低频声波，经过数据采集模块、互相关分析模块计算出声波传播时间；然后经校准模块、声速计算模块准确计算出声速，最后空气温度计算模块根据空气中声速和温度的关系结合声速计算模块求得的声速、湿度计测量的湿度和气压计提供的大气压值准确计算出空气温度，本发明相比于现有的接触式温度测量，可以非侵入式的测量环境温度，所提供的空气温度测量装置结构简单、成本低廉、易于安装和维护，适用于空气温度的快速、准确、非接触式测量，具有广阔的应用前景。

Description

一种基于低频声波的空气温度测量装置及方法
技术领域
本发明属于环境温度测量技术领域，特别涉及一种基于低频声波的空气温度测量装置及方法。
背景技术
环境温度是工业过程中的一个重要测量参数，很多工业过程需要对空气温度进行严格监控，以保障生产顺利进行，保证产品质量。例如，仓库内的空气温度测量可以帮助保证产品储存质量、发现潜在的热源；温室内的空气温度信息可以帮助控制温室内环境稳定，促进植物生长，减少能源消耗；建筑室内的空气温度信息可以帮助设计建筑内气流组织、空调和采暖，有效节约能源，提高生活质量。
目前空气温度的测量主要通过接触式测量仪表和非接触式仪表实现。接触式测温仪表如热电偶、热电阻等只能提供点测温结果并且需要与被测介质接触。非接触式测温仪表如红外测温仪等受物体表面发射率、测量距离、烟尘和水汽等外界因素的影响，测量误差较大。声波在空气中传播时，声波的速度主要由空气温度决定，通过测量空气中声波的速度可推断空气温度。本发明提出了一种基于低频声波的空气温度测量装置及方法，通过声源发出低频声波，然后对两个传声器接收到的声波信号进行分析处理得到声速信息，最后根据空气中声速和温度、湿度、大气压的物理关系获得空气温度信息。此方法是一种非侵入式的测温技术，可以快速、准确地测量空气温度。
发明内容
本发明的目的是提供一种基于低频声波的空气温度测量装置及方法，其特征在于：所述基于低频声波的空气温度测量装置的传感器模块1由扬声器2、上游传声器3、下游传声器4、湿度计5和固定装置6组成；固定装置6为倒п字形，扬声器2固定在п字形左边，上游传声器3、下游传声器4和湿度计5依次间隔固定在固定装置6的底板上，声源信号生成模块7、功率放大器8和扬声器2串联；上游传声器3、下游传声器4分别连接数据采集模块9，数据采集模块9、互相关分析模块10、声速计算模块12和空气温度计算模块14串联；空气温度计算模块14再分别连接湿度计5与气压计13；校准模块11分别连接互相关分析模块10和声速计算模块12；空气温度计算模块14输出温度信号。
所述扬声器2的中心和上游传声器3、下游传声器4的感应元件处于同一高度，扬声器2、上游传声器3、下游传声器4和湿度计5位于同一直线上。
所述基于低频声波的空气温度测量装置的空气温度测量方法，其特征在于，包括如下步骤：
1)首先传感器模块根据上游传声器3、下游传声器4间的距离对扬声器2的发射功率进行设定，以保证下游传声器4接收的声波没有完全衰减，并且具有高的信噪比；
2)声源信号生成模块7，生成0.1秒的200-1500Hz线性扫频信号，经过功率放大器8驱动扬声器2以发出低频声波；
3)传声器接收扬声器发出的低频声波，经过数据采集模块、互相关分析模块计算出声波传播时间；
4)将传感器模块置于恒温室中，测得不同温度下的声波传播时间，由校准模块计算出等效声波路径长度和系统时延；
5)经校准模块、声速计算模块准确计算出声速；
6)最后空气温度计算模块根据空气中声速和温度的关系，结合声速计算模块求得的声速、湿度计测量的湿度和气压计提供的大气压值准确计算出空气温度。
所述数据采集模块对两路声信号进行高速同步采集，采样频率大于1Msps，保证温度引起的微小的声波传播时间变化能够检测出。
所述互相关分析模块只使用传声器采集的声信号的直达声部分，避免声波反射引起的声波传播时间计算误差。
所述步骤5)中声速计算模块准确计算出声速需要校准模块提供两个传声器间的等效声波路径长度和系统时延，如果传声器位置发生变化或者测量地点发生变化，系统需要重新校准。
所述空气温度计算模块根据测得的声速计算温度时，需要输入湿度计同步测量的湿度信息以及气压计提供的实时大气压值。
声源信号生成模块7，生成0.1秒的200-1500Hz线性扫频信号，经过功率放大器8驱动扬声器2以发出低频声波；上游传声器3和下游传声器4采集空气中的声波并由数据采集模块9同步采集转化为数字信号，使用互相关分析模块10对两路传声器接收的声信号的直达声部分进行互相关运算，测得声波在上游传声器3和下游传声器4间的传播时间t；将传感器模块1置于恒温室中并测得不同空气温度下的声波传播时间，经过校准模块11得到实际声波路径长度L和系统时延τ；声速计算模块12根据互相关分析模块10计算的声波传播时间t、校准模块11提供的实际声波路径长度L和系统时间τ可准确得到空气中声速c；空气温度计算模块14根据声速c、湿度计5同步测量的空气湿度信息以及气压计提供的大气压值来计算出空气温度T
所述基于低频声波的空气温度测量方法是根据空气中声速c和空气温度T的物理关系：
其中，γ为比热比，Z为压缩性因子，R为通用气体常数，Ma为干空气摩尔质量，Mw为水蒸汽摩尔质量，xw为水蒸汽摩尔分数，RH为相对湿度，f为水蒸汽强化因子，psv为饱和蒸汽压，p为大气压；当空气温度T升高时，声速c增加，两个传声器间的声波传播时间降低，反之声波传播时间升高；实际声波路径长度L和系统时间τ可准确得到空气中声速c；
空气温度计算模块14根据式(3)声速c、湿度计5同步测量的空气湿度信息以及气压计提供的大气压值，以式(1)(2)来计算出空气温度T。
本发明的有益效果是：本发明通过扬声器发出低频扫频声波，对两个传声器接收的声波进行互相关测得声波传播时间，通过校准获得实际声波路径长度和系统时间进而提高声速测量精度。将测量的声速、湿度计同步测量的空气湿度以及大气压信息提供给空气温度计算模块，即可计算出空气温度T。
所提供的空气温度测量装置结构简单、成本低廉、易于安装和维护，适用于空气温度的快速、准确、非接触式测量，具有广阔的应用前景。
附图说明
图1为基于低频声波的空气温度测量装置原理结构示意图。
图中：1.传感器模块；2.扬声器；3.上游传声器；4.下游传声器；5.湿度计；6.固定装置；7.声源信号生成模块；8.功率放大器；9.数据采集模块；10.互相关分析模块；11.校准模块；12.声速计算模块；13.气压计；14.空气温度计算模块。
具体实施方式
本发明提供一种基于低频声波的空气温度测量装置及方法。下面结合附图和实施例对本发明予以说明。
图1所示为基于低频声波的空气温度测量装置原理结构示意图。图中所示的基于低频声波的空气温度测量装置，其传感器模块1由扬声器2、上游传声器3、下游传声器4、湿度计5和固定装置6组成；固定装置6为倒п字形，扬声器2固定在п字形左边，上游传声器3、下游传声器4和湿度计5依次间隔固定在固定装置6的底板上，声源信号生成模块7、功率放大器8和扬声器2串联；上游传声器3、下游传声器4分别连接数据采集模块9，数据采集模块9、互相关分析模块10、声速计算模块12和空气温度计算模块14串联；空气温度计算模块14再分别连接湿度计5与气压计13；校准模块11分别连接互相关分析模块10和声速计算模块12；空气温度计算模块14输出温度信号。其中扬声器2的中心和上游传声器3、下游传声器4的感应元件处于同一高度，扬声器2、上游传声器3、下游传声器4和湿度计5位于同一直线上。
传感器模块1放置在测温环境中，声源信号生成模块7生成0.1秒的200-1500Hz线性扫频信号，经过功率放大器8驱动扬声器2以发出低频声波；其扬声器2的发射功率根据上游传声器3、下游传声器4间的距离进行设定，以保证下游传声器4接收的声波没有完全衰减并且具有较高的信噪比。
上游传声器3和下游传声器4采集空气中的声波并由数据采集模块9同步采集转化为数字信号，使用互相关分析模块10对两路传声器接收的声信号的直达声部分进行互相关运算，测得声波在上游传声器3和下游传声器4间的传播时间t；将传感器模块1置于恒温室中并测得不同空气温度下的声波传播时间，经过校准模块11得到实际声波路径长度L和系统时延τ；声速计算模块12根据互相关分析模块10计算的声波传播时间t、校准模块11提供的实际声波路径长度L和系统时间τ可准确得到空气中声速c；空气温度计算模块14根据声速c、湿度计5同步测量的空气湿度信息以及气压计提供的大气压值来计算出空气温度T。
数据采集模块9的采样频率大于1Msps。声速计算模块12需要校准模块11提供上游传声器3、下游传声器4间的等效声波路径长度L和系统时间τ，如果上游传声器3、下游传声器4位置发生变化或者测量地点发生变化，系统需要重新校准。
本发明基于低频声波的空气温度测量方法是根据空气中声速c和空气温度T的物理关系：
其中，γ为比热比，Z为压缩性因子，R为通用气体常数，Ma为干空气摩尔质量，Mw为水蒸汽摩尔质量，xw为水蒸汽摩尔分数，RH为相对湿度，f为水蒸汽强化因子，psv为饱和蒸汽压，p为大气压。
当空气温度T升高时，声速c增加，上游传声器3、下游传声器4间的声波传播时间t降低，反之声波传播时间t升高。本发明通过扬声器2发出低频扫频声波，对上游传声器3、下游传声器4接收的声波进行互相关测得声波传播时间t，通过校准获得实际声波路径长度L和系统时间τ，声速c可计算为：
利用声速c、湿度计5同步测量的空气湿度信息RH以及气压计提供的大气压值p可根据式(1)、(2)计算出空气温度T。

Claims (9)

1.一种基于低频声波的空气温度测量装置，其特征在于：所述基于低频声波的空气温度测量装置的传感器模块(1)由扬声器(2)、上游传声器(3)、下游传声器(4)、湿度计(5)和固定装置(6)组成；固定装置(6)为倒п字形，扬声器(2)固定在п字形左边，上游传声器(3)、下游传声器(4)和湿度计(5)依次间隔固定在固定装置(6)的底板上，声源信号生成模块(7)、功率放大器(8)和扬声器(2)串联；上游传声器(3)、下游传声器(4)分别连接数据采集模块(9)，数据采集模块(9)、互相关分析模块(10)、声速计算模块(12)和空气温度计算模块(14)串联；空气温度计算模块(14)再分别连接湿度计(5)与气压计(13)；校准模块(11)分别连接互相关分析模块(10)和声速计算模块(12)；空气温度计算模块(14)输出温度信号。
2.根据权利要求1所述基于低频声波的空气温度测量装置，其特征在于：所述扬声器(2)的中心和上游传声器(3)、下游传声器(4)的感应元件处于同一高度，扬声器(2)、上游传声器(3)、下游传声器(4)和湿度计(5)位于同一直线上。
3.一种基于低频声波的空气温度测量装置的空气温度测量方法，其特征在于，包括如下步骤：
1)首先传感器模块根据上游传声器(3)、下游传声器(4)间的距离对扬声器(2)的发射功率进行设定，以保证下游传声器(4)接收的声波没有完全衰减，并且具有高的信噪比；
2)声源信号生成模块(7，生成0.1秒的200-1500Hz线性扫频信号，经过功率放大器(8驱动扬声器(2以发出低频声波；
3)传声器接收扬声器发出的低频声波，经过数据采集模块、互相关分析模块计算出声波传播时间；
4)将传感器模块置于恒温室中，测得不同温度下的声波传播时间，由校准模块计算出等效声波路径长度和系统时延；
5)经校准模块、声速计算模块准确计算出声速；
6)最后空气温度计算模块根据空气中声速和温度的关系，结合声速计算模块求得的声速、湿度计测量的湿度和气压计提供的大气压值准确计算出空气温度。
4.根据权利要求3所述基于低频声波的空气温度测量装置的空气温度测量方法，其特征在于，所述数据采集模块对两路声信号进行高速同步采集，采样频率大于1Msps，保证温度引起的微小的声波传播时间变化能够检测出。
5.根据权利要求3所述基于低频声波的空气温度测量装置的空气温度测量方法，其特征在于，所述互相关分析模块只使用传声器采集的声信号的直达声部分，避免声波反射引起的声波传播时间计算误差。
6.根据权利要求3所述基于低频声波的空气温度测量装置的空气温度测量方法，其特征在于，所述步骤5)中声速计算模块准确计算出声速需要校准模块提供两个传声器间的等效声波路径长度和系统时延，如果传声器位置发生变化或者测量地点发生变化，系统需要重新校准。
7.根据权利要求3所述基于低频声波的空气温度测量装置的空气温度测量方法，其特征在于，所述空气温度计算模块根据测得的声速计算温度时，需要输入湿度计同步测量的湿度信息以及气压计提供的实时大气压值。
8.根据权利要求3所述基于低频声波的空气温度测量装置的空气温度测量方法，其特征在于，声源信号生成模块(7)，生成0.1秒的200-1500Hz线性扫频信号，经过功率放大器(8)驱动扬声器(2)以发出低频声波；上游传声器(3)和下游传声器(4)采集空气中的声波并由数据采集模块(9)同步采集转化为数字信号，使用互相关分析模块(10)对两路传声器接收的声信号的直达声部分进行互相关运算，测得声波在上游传声器(3)和下游传声器(4)间的传播时间t；将传感器模块(1)置于恒温室中并测得不同空气温度下的声波传播时间，经过校准模块(11)得到实际声波路径长度L和系统时延τ；声速计算模块(12)根据互相关分析模块(10)计算的声波传播时间t、校准模块(11)提供的实际声波路径长度L和系统时间τ可准确得到空气中声速c；空气温度计算模块(14)根据声速c、湿度计(5)同步测量的空气湿度信息以及气压计提供的大气压值来计算出空气温度T。
9.根据权利要求3所述基于低频声波的空气温度测量装置的空气温度测量方法，其特征在于，所述基于低频声波的空气温度测量方法是根据空气中声速c和空气温度T的物理关系：
其中，γ为比热比，Z为压缩性因子，R为通用气体常数，Ma为干空气摩尔质量，Mw为水蒸汽摩尔质量，xw为水蒸汽摩尔分数，RH为相对湿度，f为水蒸汽强化因子，psv为饱和蒸汽压，p为大气压；当空气温度T升高时，声速c增加，两个传声器间的声波传播时间降低，反之声波传播时间升高；实际声波路径长度L和系统时间τ可准确得到空气中声速c；
空气温度计算模块14根据式(3)声速c、湿度计5同步测量的空气湿度信息以及气压计提供的大气压值，以式(1)(2)来计算出空气温度T。
CN201710971935.8A 2017-10-18 2017-10-18 一种基于低频声波的空气温度测量装置及方法 Expired - Fee Related CN107917768B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710971935.8A CN107917768B (zh) 2017-10-18 2017-10-18 一种基于低频声波的空气温度测量装置及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710971935.8A CN107917768B (zh) 2017-10-18 2017-10-18 一种基于低频声波的空气温度测量装置及方法

Publications (2)

Publication Number Publication Date
CN107917768A true CN107917768A (zh) 2018-04-17
CN107917768B CN107917768B (zh) 2020-02-21

Family

ID=61895695

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710971935.8A Expired - Fee Related CN107917768B (zh) 2017-10-18 2017-10-18 一种基于低频声波的空气温度测量装置及方法

Country Status (1)

Country Link
CN (1) CN107917768B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109443515A (zh) * 2018-09-13 2019-03-08 中国船舶重工集团公司第七〇五研究所 一种小尺寸水听器灵敏度空气中测试系统及方法
CN110231084A (zh) * 2019-06-25 2019-09-13 中国科学院声学研究所 一种用于低温低气压的声速测量装置及方法
CN111289035A (zh) * 2018-12-06 2020-06-16 罗斯蒙特航天公司 声学空气数据传感器和系统
CN112595431A (zh) * 2020-10-22 2021-04-02 华帝股份有限公司 一种用于烹饪设备的温度检测方法
CN114279591A (zh) * 2021-12-22 2022-04-05 重庆长安汽车股份有限公司 一种基于声传播的气温测量装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002054999A (ja) * 2000-08-11 2002-02-20 Takenaka Komuten Co Ltd 音波発生検出装置、環境状態測定装置及び方法
CN101403643A (zh) * 2008-11-14 2009-04-08 华北电力大学 一种利用声学测量炉膛烟气温度的装置
CN201307039Y (zh) * 2008-11-18 2009-09-09 华北电力大学 声学测量炉膛烟气温度的装置
CN102818652A (zh) * 2012-08-15 2012-12-12 华北电力大学 一种基于声学测温的铝包内铝水温度在线监测系统及方法
CN104132749A (zh) * 2014-08-25 2014-11-05 哈尔滨理大晟源科技开发有限公司 声波测定炉膛温度装置及测定方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002054999A (ja) * 2000-08-11 2002-02-20 Takenaka Komuten Co Ltd 音波発生検出装置、環境状態測定装置及び方法
CN101403643A (zh) * 2008-11-14 2009-04-08 华北电力大学 一种利用声学测量炉膛烟气温度的装置
CN201307039Y (zh) * 2008-11-18 2009-09-09 华北电力大学 声学测量炉膛烟气温度的装置
CN102818652A (zh) * 2012-08-15 2012-12-12 华北电力大学 一种基于声学测温的铝包内铝水温度在线监测系统及方法
CN104132749A (zh) * 2014-08-25 2014-11-05 哈尔滨理大晟源科技开发有限公司 声波测定炉膛温度装置及测定方法

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109443515A (zh) * 2018-09-13 2019-03-08 中国船舶重工集团公司第七〇五研究所 一种小尺寸水听器灵敏度空气中测试系统及方法
CN109443515B (zh) * 2018-09-13 2021-05-14 中国船舶重工集团公司第七一五研究所 一种小尺寸水听器灵敏度空气中测试系统及方法
CN111289035A (zh) * 2018-12-06 2020-06-16 罗斯蒙特航天公司 声学空气数据传感器和系统
CN111289035B (zh) * 2018-12-06 2023-10-31 罗斯蒙特航天公司 声学空气数据传感器和系统
CN110231084A (zh) * 2019-06-25 2019-09-13 中国科学院声学研究所 一种用于低温低气压的声速测量装置及方法
CN110231084B (zh) * 2019-06-25 2024-05-28 中国科学院声学研究所 一种用于低温低气压的声速测量装置及方法
CN112595431A (zh) * 2020-10-22 2021-04-02 华帝股份有限公司 一种用于烹饪设备的温度检测方法
CN114279591A (zh) * 2021-12-22 2022-04-05 重庆长安汽车股份有限公司 一种基于声传播的气温测量装置

Also Published As

Publication number Publication date
CN107917768B (zh) 2020-02-21

Similar Documents

Publication Publication Date Title
CN107917768A (zh) 一种基于低频声波的空气温度测量装置及方法
Kaimal et al. A continuous wave sonic anemometer-thermometer
CN105044676A (zh) 一种基于能量的声源定位方法
CN104048808A (zh) 一种动态熵探针
CN114001804B (zh) 一种基于时间差法超声波计量装置的校准方法及系统
CN102062636A (zh) 便携式现场激光能量测量装置
CN102269769A (zh) 超声波三维测风方法和三维超声波风速风向仪
CN102539009A (zh) 电站锅炉水冷壁向火侧壁温在线监测系统及监测方法
CN103926018B (zh) 一种液体温度测量装置及方法
CN202886398U (zh) 一种利用热电偶的热球式风速计
CN103336052B (zh) 一种室内相对湿度在线监测系统及湿度计算方法
CN201583657U (zh) 带温度和湿度补偿的超声波测距装置
CN107228630A (zh) 一种双金属片形变测试装置
CN102980941B (zh) 一种利用声速变化的超声气体相对湿度检测方法及装置
CN109253781A (zh) 节流式压差流量传感器的校准方法及校准系统
CN102539111B (zh) 一种基于dsPIC单片机的温度脉动仪
CN201772928U (zh) 一种数字化红外冠层温度传感器
CN103345003B (zh) 探空仪防辐射装置以及利用探空仪获取中空温度的方法
CN202770855U (zh) 一种基于声学的锅炉炉膛烟气流速监测系统
Huang et al. Temperature measurement system based on ultrasonic phase-shift method
TW202305407A (zh) 大氣亂流偵測方法及大氣亂流偵測裝置
CN203688579U (zh) 风速计
CN203164156U (zh) 一种利用声速变化的超声气体相对湿度检测装置
CN202453111U (zh) 电站锅炉水冷壁向火侧壁温在线监测系统
CN104965103A (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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20200221

Termination date: 20201018