CN108362483B - 一种管路系统的监测方法和监测系统 - Google Patents
一种管路系统的监测方法和监测系统 Download PDFInfo
- Publication number
- CN108362483B CN108362483B CN201711304174.7A CN201711304174A CN108362483B CN 108362483 B CN108362483 B CN 108362483B CN 201711304174 A CN201711304174 A CN 201711304174A CN 108362483 B CN108362483 B CN 108362483B
- Authority
- CN
- China
- Prior art keywords
- vibration signal
- pipeline
- pipeline system
- phi
- monitoring
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M13/00—Testing of machine parts
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Pipeline Systems (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明适用于数据处理领域,提供了一种管路系统的监测方法和监测系统,所述监测方法包括:获取预设于待监测管路系统的关键位置的振动加速度传感器传输的振动信号,将振动信号转换为频域振动信号,将所述频域振动信号添加到公式Φ(0,s)=UΦ(z,s)+Q(z,s)中,根据所述公式的结果,判断所述待监测管路系统的健康程度。实施本发明实施例,监测系统通过管路系统的关键点的振动参数即可对整个管路系统进行监测,提升了监测的诊断速度和适应性。
Description
技术领域
本发明属于数据处理领域,尤其涉及一种管路系统的监测方法和监测系统。
背景技术
大型复杂管路系统船舶与海洋工程、石油化工、能源与电力工业等诸多领域中应用广泛。管路系统中的泵、阀门等设备及管路附件出现故障时,常会造成管路系统及与之相连的其他元件和精密仪器的破坏,影响管路系统安全和管路动力系统的正常运行,严重时会造成巨大的经济损失,因此对管路系统运行状态进行监测并实现故障定位一直是国内外研究的热点问题。
目前在管路运输领域应用管路系统监测的较多,监测方法包括负压波检漏法、管线附近地面震动监测法和沿管线敷设光线的信号扰动测量法等,其中相对成熟的是负压波检漏法,该方法是利用流体泄漏时因局部瞬时压力差和速度差而产生的负压波作为监测对象,通过拾取压力波的梯度特征和压力变化率的时间差和信号相关处理确定泄露程度和泄漏位置。这一方法准确度受管路背景压力、介质密度、温度等因素影响较大,并且不适用于含有多个激励源的大型复杂动力管路系统的运行监测。
管路系统健康运行监测的核心问题在于管路系统动力学计算方法。目前常见的管路系统的动力学特性的预测方法有特征线方法,有限元法,特征线-有限元法方法和传递矩阵方法。与其它方法相比,传递矩阵方法是通过波动方程直接描述管路系统的运动。与其他方法相比,传递矩阵方法求解过程简单,易于编程且计算效率高,因而在管路结构振动,管内流体压力脉动及管路系统动力学预测等问题中具有较多应用。
针对传递矩阵方法的计算不稳定、积累数值误差大及数值溢出等“病态问题”,常见的几种处理方法包括传递矩阵的无量纲化,特征值问题的改进算法以及Riccati传递矩阵方法等,此外,通过一个一维数组和两个变量来存储一个数据的方法解决了超精度数据的存储问题该方法在一定程度上能够减少由于矩阵变化带来的计算误差,上述这些的改进措施,在一定程度上均提高了传递矩阵在大型复杂管路系统动力学特性预测精度和预测结果的稳定性,但利用这些方法在计算管路系统时并无优势,难以实现模块化编程及预测。
由于受到管道系统动力学计算方法的制约,目前大型复杂管路系统健康运行监测系统发展缓慢。
发明内容
本发明实施例的目的在于提供一种管路系统的监测方法和监测系统,以解决现有技术判断管路系统健康状态较慢的问题。
本发明实施例是这样实现的,一种管路系统的监测方法,所述监测方法包括:
获取预设于待监测管路系统的关键位置的振动加速度传感器传输的振动信号;
将所述振动信号转换为频域振动信号;
将所述频域振动信号添加到公式Φ(0,s)=UΦ(z,s)+Q(z,s)中,根据所述公式的结果,判断所述待监测管路系统的健康程度,所述Φ(0,s)为管路始端的状态向量,Φ(z,s)为管路任意点的状态向量,U为管段场传递矩阵,Q(z,s)是由外部激励力、重力以及场传递矩阵等共同作用组成的状态向量。
本发明实施例的另一目的在于提供一种管路系统的监测系统,所述监测系统包括:
振动信号获取单元,用于获取预设于待监测管路系统的关键位置的振动加速度传感器传输的振动信号;
转换单元,用于将所述振动信号转换为频域振动信号;
健康程度判断单元,用于将所述频域振动信号添加到公式Φ(0,s)=UΦ(z,s)+Q(z,s)中,根据所述公式的结果,判断所述待监测管路系统的健康程度,所述Φ(0,s)为管路始端的状态向量,Φ(z,s)为管路任意点的状态向量,U为管段场传递矩阵,Q(z,s)是由外部激励力、重力以及场传递矩阵等共同作用组成的状态向量。
本发明实施例,获取预设于待监测管路系统的关键位置的振动加速度传感器传输的振动信号,将振动信号转换为频域振动信号,将所述频域振动信号添加到公式Φ(0,s)=UΦ(z,s)+Q(z,s)中,根据所述公式的结果,判断所述待监测管路系统的健康程度,使得监测系统通过管路系统的关键点的振动参数即可对整个管路系统进行监测,提升了监测的诊断速度和适应性。
附图说明
图1为本发明一示例性实施例示出的一种管路系统的监测方法的流程图;
图2为本发明一示例性实施例示出的一种管路系统的监测系统的结构图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
为了说明本发明所述的技术方案,下面通过具体实施例来进行说明。
如图1所示为本发明一示例性实施例示出的一种管路系统的监测方法的流程图,所述监测方法包括以下步骤:
步骤S101,获取预设于待监测管路系统的关键位置的振动加速度传感器传输的振动信号。
在本发明实施例中,根据待监测管路系统的组成和管路系统的实际布置形式,可以选择出多个关键位置,在这些关键位置布放振动加速度传感器,通过振动加速度传感器传输的振动信号,监测系统可以实时了解待监测管路系统的振动情况,需要指出的是,振动加速度传感器放置的位置、数量等,可根据实际使用的需要自行设置,在此不作赘述。
优选的,关键位置包括但不限于:泵源、阀门、管路系统考察点。
步骤S102,将所述振动信号转换为频域振动信号。
在本发明实施例中,振动加速度传感器传输的振动信号是时域信号,时域信号并不能直接用于计算,因此需要将振动信号转换为频域振动信号。将时域信号转换为频域信号的方法为现有技术,在此不做赘述。
优选的,所述振动信号通过傅里叶变换转换为频域振动信号。
步骤S103,将所述频域振动信号添加到公式Φ(0,s)=UΦ(z,s)+Q(z,s)中,根据所述公式的结果,判断所述待监测管路系统的健康程度,所述Φ(0,s)为管路始端的状态向量,Φ(z,s)为管路任意点的状态向量,U为管段场传递矩阵,Q(z,s)是由外部激励力、重力以及场传递矩阵等共同作用组成的状态向量。
在本发明实施例中,监测系统中预设有一运行参数库,该运行参数库记录了管路系统健康程度的判断标准,将通过公式Φ(0,s)=UΦ(z,s)+Q(z,s)计算的结果与运行参数库中的判断标准进行比较,根据比较的结果即可判断待监测管路系统的健康程度。
通过上述公式,在测量任意点的振动后,即可获取管路始端的状态向量,将该状态向量与运行参数库进行比较,即可判断管路系统的健康状态。
需要指出的是,运行参数库可以在系统运行前植入,也可以在系统投入正常运行时,通过采集管路系统的健康工作状态的数据取得。
所述判断所述待监测管路系统的健康程度,包括:
如果所述Φ(0,s)超过所述运行参数库的健康阈值,则所述待监测管路系统为不健康;以及,
如果所述Φ(0,s)未超过所述运行参数库的健康阈值,则所述待监测管路系统为健康。
本发明实施例,获取预设于待监测管路系统的关键位置的振动加速度传感器传输的振动信号,将振动信号转换为频域振动信号,将所述频域振动信号添加到公式Φ(0,s)=UΦ(z,s)+Q(z,s)中,根据所述公式的结果,判断所述待监测管路系统的健康程度,使得监测系统通过管路系统的关键点的振动参数即可对整个管路系统进行监测,提升了监测的诊断速度和适应性。
作为本发明的一个可选实施例,在所述将所述频域振动信号添加到公式Φ(0,s)=UΦ(z,s)+Q(z,s)中,根据所述公式的结果,判断所述待监测管路系统的健康程度的步骤之后,所述监测方法还包括:
在显示设备上显示所述判断的结果。
在本发明实施例中,监测系统可以将健康判断的结果显示在显示设备上,监测人员即可快速获取健康监测的结果。
作为本发明的另一个可选实施例,在所述将所述频域振动信号添加到公式Φ(0,s)=UΦ(z,s)+Q(z,s)中,根据所述公式的结果,判断所述待监测管路系统的健康程度的步骤之后,所述监测方法还包括:
当所述待监测管路系统为不健康时,进行报警处理。
如图2所示为本发明一示例性实施例示出的一种管路系统的监测系统的结构图,所述监测系统包括:
振动信号获取单元201,用于获取预设于待监测管路系统的关键位置的振动加速度传感器传输的振动信号。
在本发明实施例中,根据待监测管路系统的组成和管路系统的实际布置形式,可以选择出多个关键位置,在这些关键位置布放振动加速度传感器,通过振动加速度传感器传输的振动信号,监测系统可以实时了解待监测管路系统的振动情况,需要指出的是,振动加速度传感器放置的位置、数量等,可根据实际使用的需要自行设置,在此不作赘述。
优选的,关键位置包括但不限于:泵源、阀门、管路系统考察点。
转换单元202,用于将所述振动信号转换为频域振动信号。
在本发明实施例中,振动加速度传感器传输的振动信号是时域信号,时域信号并不能直接用于计算,因此需要将振动信号转换为频域振动信号。将时域信号转换为频域信号的方法为现有技术,在此不做赘述。
优选的,所述振动信号通过傅里叶变换转换为频域信号。
健康程度判断单元203,用于将所述频域振动信号添加到公式Φ(0,s)=UΦ(z,s)+Q(z,s)中,根据所述公式的结果,判断所述待监测管路系统的健康程度,所述Φ(0,s)为管路始端的状态向量,Φ(z,s)为管路任意点的状态向量,U为管段场传递矩阵,Q(z,s)是由外部激励力、重力以及场传递矩阵等共同作用组成的状态向量。
在本发明实施例中,监测系统中预设有一运行参数库,该运行参数库记录了管路系统健康程度的判断标准,将通过公式Φ(0,s)=UΦ(z,s)+Q(z,s)计算的结果与运行参数库中的判断标准进行比较,根据比较的结果即可判断待监测管路系统的健康程度。
通过上述公式,在测量任意点的振动后,即可获取管路始端的状态向量,将该状态向量与运行参数库进行比较,即可判断管路系统的健康状态。
需要指出的是,运行参数库可以在系统运行前植入,也可以在系统投入正常运行时,通过采集管路系统的健康工作状态的数据取得。
所述判断所述待监测管路系统的健康程度,包括:
如果所述Φ(0,s)超过所述运行参数库的健康阈值,则所述待监测管路系统为不健康;以及,
如果所述Φ(0,s)未超过所述运行参数库的健康阈值,则所述待监测管路系统为健康。
本发明实施例,获取预设于待监测管路系统的关键位置的振动加速度传感器传输的振动信号,将振动信号转换为频域振动信号,将所述频域振动信号添加到公式Φ(0,s)=UΦ(z,s)+Q(z,s)中,根据所述公式的结果,判断所述待监测管路系统的健康程度,使得监测系统通过管路系统的关键点的振动参数即可对整个管路系统进行监测,提升了监测的诊断速度和适应性。
作为本发明的一个可选实施例,所述监测系统还包括:
显示单元,用于在显示设备上显示所述判断的结果。
在本发明实施例中,监测系统可以将健康判断的结果显示在显示设备上,监测人员即可快速获取健康监测的结果。
作为本发明的另一个可选实施例,所述监测系统还包括:
报警单元,用于当所述待监测管路系统为不健康时,进行报警处理。
本领域普通技术人员可以理解为上述实施例所包括的各个单元只是按照功能逻辑进行划分的,但并不局限于上述的划分,只要能够实现相应的功能即可;另外,各功能单元的具体名称也只是为了便于相互区分,并不用于限制本发明的保护范围。
本领域普通技术人员还可以理解,实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,所述的程序可以在存储于一计算机可读取存储介质中,所述的存储介质,包括ROM/RAM、磁盘、光盘等。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种管路系统的监测方法,其特征在于,所述监测方法包括:
获取预设于待监测管路系统的关键位置的振动加速度传感器传输的振动信号;
将所述振动信号转换为频域振动信号;
将所述频域振动信号添加到公式Φ(0,s)=UΦ(z,s)+Q(z,s)中,根据所述公式的结果,判断所述待监测管路系统的健康程度,如果所述Φ(0,s)超过运行参数库的健康阈值,则所述待监测管路系统为不健康,如果所述Φ(0,s)未超过所述运行参数库的健康阈值,则所述待监测管路系统为健康,所述Φ(0,s)为管路始端的状态向量,Φ(z,s)为管路任意点的状态向量,U为管段场传递矩阵,Q(z,s)是由外部激励力、重力以及场传递矩阵共同作用组成的状态向量。
2.如权利要求1所述的监测方法,其特征在于,所述关键位置包括:泵源、阀门、管路系统考察点。
3.如权利要求1所述的监测方法,其特征在于,所述振动信号通过傅里叶变换转换为频域振动信号。
4.如权利要求1~3任一项所述的监测方法,其特征在于,在所述将所述频域振动信号添加到公式Φ(0,s)=UΦ(z,s)+Q(z,s)中,根据所述公式的结果,判断所述待监测管路系统的健康程度的步骤之后,所述监测方法还包括:
在显示设备上显示所述判断的结果。
5.一种管路系统的监测系统,其特征在于,所述监测系统包括:
振动信号获取单元,用于获取预设于待监测管路系统的关键位置的振动加速度传感器传输的振动信号;
转换单元,用于将所述振动信号转换为频域振动信号;
健康程度判断单元,用于将所述频域振动信号添加到公式Φ(0,s)=UΦ(z,s)+Q(z,s)中,根据所述公式的结果,判断所述待监测管路系统的健康程度,如果所述Φ(0,s)超过运行参数库的健康阈值,则所述待监测管路系统为不健康,如果所述Φ(0,s)未超过所述运行参数库的健康阈值,则所述待监测管路系统为健康,所述Φ(0,s)为管路始端的状态向量,Φ(z,s)为管路任意点的状态向量,U为管段场传递矩阵,Q(z,s)是由外部激励力、重力以及场传递矩阵共同作用组成的状态向量。
6.如权利要求5所述的监测系统,其特征在于,所述关键位置包括:泵源、阀门、管路系统考察点。
7.如权利要求5所述的监测系统,其特征在于,所述振动信号通过傅里叶变换转换为频域振动信号。
8.如权利要求5~7任一项所述的监测系统,其特征在于,所述监测系统还包括:显示单元,用于在显示设备上显示所述判断的结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711304174.7A CN108362483B (zh) | 2017-12-11 | 2017-12-11 | 一种管路系统的监测方法和监测系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711304174.7A CN108362483B (zh) | 2017-12-11 | 2017-12-11 | 一种管路系统的监测方法和监测系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108362483A CN108362483A (zh) | 2018-08-03 |
CN108362483B true CN108362483B (zh) | 2021-05-18 |
Family
ID=63009886
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711304174.7A Active CN108362483B (zh) | 2017-12-11 | 2017-12-11 | 一种管路系统的监测方法和监测系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108362483B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09159574A (ja) * | 1995-12-04 | 1997-06-20 | Olympus Optical Co Ltd | 光学性能評価装置及び光学性能評価方法 |
JPH10221146A (ja) * | 1997-02-04 | 1998-08-21 | Yokogawa Electric Corp | コリオリ質量流量計 |
CA2694936A1 (en) * | 2007-08-06 | 2009-02-12 | Rosemount Inc. | Process variable transmitter with acceleration sensor |
CN103049670A (zh) * | 2013-01-14 | 2013-04-17 | 哈尔滨工程大学 | 管道激励源识别及其振动响应预测方法 |
CN104634870A (zh) * | 2014-12-24 | 2015-05-20 | 同济大学 | 一种基于振动响应测试的隧道结构损伤识别装置 |
CN104898023A (zh) * | 2015-06-10 | 2015-09-09 | 国网上海市电力公司 | 一种基于传递函数曲线的接地网故障诊断方法 |
CN106950131A (zh) * | 2017-03-16 | 2017-07-14 | 招商局重庆交通科研设计院有限公司 | 锚固结构钢束工作应力及同束不均匀度的无损检测方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5160980B2 (ja) * | 2008-07-08 | 2013-03-13 | オークマ株式会社 | 振動抑制方法及び装置 |
CN102095573A (zh) * | 2009-12-11 | 2011-06-15 | 上海卫星工程研究所 | 星载旋转设备机械部件状态监测与诊断报警方法 |
CN103674578B (zh) * | 2012-09-14 | 2016-07-06 | 中车青岛四方机车车辆股份有限公司 | 高速列车运行动力学性能状态的检测方法 |
CN103245513B (zh) * | 2013-04-10 | 2015-06-10 | 暨南大学 | 一种汽车产品整车装配的动态质量检测方法 |
CN103528844B (zh) * | 2013-10-23 | 2016-02-17 | 青岛理工大学 | 基于经验模态分解的结构损伤预警方法 |
CN104573255B (zh) * | 2015-01-22 | 2017-10-31 | 天津理工大学 | 一种基于改进多体系统传递矩阵的机械臂建模与求解方法 |
CN105241660B (zh) * | 2015-11-09 | 2019-04-23 | 西南交通大学 | 基于健康监测数据的高铁大型桥梁性能测试方法 |
CN105260806B (zh) * | 2015-11-25 | 2018-12-11 | 中国船舶重工集团公司第七一九研究所 | 一种管路系统的流固耦合动力学特性预测方法 |
CN106092565A (zh) * | 2016-06-08 | 2016-11-09 | 成都阜特科技股份有限公司 | 一种振动故障分析方法及其系统 |
CN107014482B (zh) * | 2017-03-29 | 2020-05-01 | 西北工业大学 | 振动状态的在线监测装置和方法 |
CN107192446B (zh) * | 2017-07-31 | 2020-12-29 | 上海电气风电集团有限公司 | 风力发电机组塔筒固有频率监测方法 |
-
2017
- 2017-12-11 CN CN201711304174.7A patent/CN108362483B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH09159574A (ja) * | 1995-12-04 | 1997-06-20 | Olympus Optical Co Ltd | 光学性能評価装置及び光学性能評価方法 |
JPH10221146A (ja) * | 1997-02-04 | 1998-08-21 | Yokogawa Electric Corp | コリオリ質量流量計 |
CA2694936A1 (en) * | 2007-08-06 | 2009-02-12 | Rosemount Inc. | Process variable transmitter with acceleration sensor |
CN103049670A (zh) * | 2013-01-14 | 2013-04-17 | 哈尔滨工程大学 | 管道激励源识别及其振动响应预测方法 |
CN104634870A (zh) * | 2014-12-24 | 2015-05-20 | 同济大学 | 一种基于振动响应测试的隧道结构损伤识别装置 |
CN104898023A (zh) * | 2015-06-10 | 2015-09-09 | 国网上海市电力公司 | 一种基于传递函数曲线的接地网故障诊断方法 |
CN106950131A (zh) * | 2017-03-16 | 2017-07-14 | 招商局重庆交通科研设计院有限公司 | 锚固结构钢束工作应力及同束不均匀度的无损检测方法 |
Non-Patent Citations (2)
Title |
---|
"基于DBN的故障特征提取及诊断方法研究";赵光权等;《仪器仪表学报》;20160915;全文 * |
"基于变分模态分解和模糊C均值聚类的滚动轴承故障诊断";刘长良等;《中国电机工程学报》;20150505;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108362483A (zh) | 2018-08-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112483908A (zh) | 一种基于声波的管道监测系统及监测方法 | |
CN104976518A (zh) | 一种海底管线泄漏监测系统 | |
CN108050396A (zh) | 一种流体管道泄漏源监测定位系统及方法 | |
CN108369118B (zh) | 使用光纤传感器对明渠中的流体流的监测 | |
CN109723977A (zh) | 一种成品油管道的停输泄漏判定方法、设备以及系统 | |
Abdulla et al. | Pipeline leak detection using artificial neural network: Experimental study | |
WO2011070343A2 (en) | Monitoring fluid flow in a conduit | |
BR112021003369A2 (pt) | método para identificar a fonte de transiente de pressão em uma rede de água | |
Islam et al. | Leak detection and location pinpointing in water pipeline systems using a wireless sensor network | |
Oven | Leak Detection in Pipelines by the use of State and Parameter Estimation | |
CN108362483B (zh) | 一种管路系统的监测方法和监测系统 | |
GB2563351A (en) | Piping diagnostic device, piping diagnostic method, discriminant-problem optimizing device, discriminant-problem optimizing method, criteria distribution ge | |
Cai et al. | A data-driven early micro-leakage detection and localization approach of hydraulic systems | |
CN110275084B (zh) | 一种针对脐带缆泄漏故障的特征识别方法 | |
CN103928066B (zh) | 核电厂凝汽器传热管检漏方法及压差检漏系统 | |
Santoso et al. | Pipeline leak detection in two phase flow based on fluctuation pressure difference and artificial neural network (ANN) | |
Mishra et al. | Testbed for real-time monitoring of leak in low pressure gas pipeline | |
Yi-bo et al. | Leakage detection and location for long range oil pipeline using negative pressure wave technique | |
KR20230137021A (ko) | 인공지능을 이용한 배관이상 탐지 방법 및 시스템 | |
Junru et al. | Hybrid modelling for leak detection of long-distance gas transport pipeline | |
Zhao et al. | Detection method for submarine oil pipeline leakage under complex sea conditions by unmanned underwater vehicle | |
CN118090043B (zh) | 一种压力表异常状态定位分析方法、系统、终端及介质 | |
Yang et al. | First Loop Loose Parts Field Diagnosis | |
Isermann et al. | Leak detection of pipelines | |
Grądzki et al. | Analysis of possibility of phase shift model usage for identification of gas leakage from pipeline with corrective elements |
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 |