CN106643987B - 一种多声道超声波气体流量计故障诊断及补偿方法 - Google Patents
一种多声道超声波气体流量计故障诊断及补偿方法 Download PDFInfo
- Publication number
- CN106643987B CN106643987B CN201610988681.6A CN201610988681A CN106643987B CN 106643987 B CN106643987 B CN 106643987B CN 201610988681 A CN201610988681 A CN 201610988681A CN 106643987 B CN106643987 B CN 106643987B
- Authority
- CN
- China
- Prior art keywords
- sound channel
- sound
- velocity
- measurement
- formula
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F25/00—Testing or calibration of apparatus for measuring volume, volume flow or liquid level or for metering by volume
- G01F25/10—Testing or calibration of apparatus for measuring volume, volume flow or liquid level or for metering by volume of flowmeters
- G01F25/15—Testing or calibration of apparatus for measuring volume, volume flow or liquid level or for metering by volume of flowmeters specially adapted for gas meters
Abstract
本发明公开了一种多声道超声波气体流量计故障诊断及补偿方法。本发明基于四声道超声波气体流量测量装置，包括四声道超声波传感器管段、传感器前端处理模块、数据采集控制模块和计算机。本发明首先通过在线声速测量确定异常声道的位置，然后舍弃异常声道的测量数据，最后通过调整其它正常工作声道对应的权重系数，利用正常工作声道的测量数据计算管道内气体的流速和流量，以弥补声道出现故障后流量计量的偏差，使得多声道超声波气体流量计在声道故障发生之后到故障排除之前仍能保证较高的计量精度。本发明为多声道超声波气体流量计的在线故障诊断及补偿提供了一种有效的解决思路，所提出的补偿方法可以明显降低因声道故障带来的流量计量的偏差。
Description
技术领域
本发明涉及多声道超声波气体流量测量领域，尤其涉及一种多声道超声波气体流量计故障诊断及补偿方法。
背景技术
在天然气、煤气等气体输送过程中，稳定可靠地计量是保证贸易交接公平、公正的重要条件。在气体计量领域，超声波气体流量计凭借测量精度高、量程比宽、无压力损耗、可双向测量等优点，在欧美等发达国家天然气贸易交接计量中获得广泛应用。国内超声波气体流量计的研究和应用起步较晚，在西气东输工程中，天然气贸易交接和管道流量计量采用了高精度的超声波气体流量计，这在国内大口径输气管道计量系统中尚属首例。超声波气体流量计根据声道数量的不同可以分为单声道超声波气体流量计和多声道超声波气体流量计。多声道超声波气体流量计与单声道相比，能更好地适应动态流场分布，具有较高的测量精度，广泛应用于大口径天然气等气体贸易计量。
在气体流量计量过程中，特别是在气体贸易交接计量等工业应用过程中，流量计必须保证持续稳定地运行，以保证气体流量计量的准确性，维护贸易交接的公平。在多声道超声波气体流量计运行的过程中，换能器、电路或其它原因故障均有可能导致流量计的某一声道测量出现异常。由于管道内气体流量是通过对各声道流速测量值加权积分计算得到的，因此当某一声道流量测量出现异常情况时，流量计量将出现较大偏差，这将造成巨大的经济损失。当出现声道故障问题后，一方面由于气体流量计量的不间断性不能马上切断管路进行故障的诊断和排除，另一方面流量计的生产厂家维修和更换流量计部件需要一定的时间。为了保证流量计在声道出现故障到问题解决这一时间段内仍能提供较高的计量精度，必须采取一定的补偿方法降低声道发生故障后流量计量的偏差。
本发明针对多声道超声波气体流量计在某一声道出现故障后会导致较大的流量计量偏差这一问题，提出了相应的故障诊断及补偿方法。该发明首先通过在线测量声速确定异常声道的位置，然后通过调整正常工作声道对应的权重系数，利用正常工作声道的测量数据进行加权计算流量，最后能够实现气体流量的正常测量。本发明可以弥补声道故障带来的流量计量偏差，能够使流量计在无法及时更换和维修的情况下仍能保证较高的计量精度，这为多声道超声波气体流量计的稳定可靠连续运行提供了一种解决思路。
发明内容
本发明的目的是克服现有技术的不足，提供一种多声道超声波气体流量计故障诊断及补偿方法。
本发明的技术方案如下：
多声道超声波气体流量计故障诊断及补偿方法包括以下步骤：
1)四声道超声波气体流量计在管道内不同高度的水平截面上布置着四对超声波换能器，四个声道沿管道中心所在水平面对称布置，即水平面上下各布置两个声道，每个声道上的两个换能器的连线与管道轴线呈一定角度。通过每对超声波换能器交替发射和接收超声波信号测得各个声道上超声波信号在顺流中的渡越时间td和在逆流中的渡越时间tu，然后根据时差法计算得到对应声道上的平均流速v。通过依次计算四个声道上的流速测量值v，获取圆形管道内不同流层上的流速信息。
2)通过判断单一声道上最大测量声速差是否在允许的范围内监测各个声道的工作状态。各个声道上的测量声速为其中，i是声道编号(1-4)，Cfi是第i声道上的测量声速，Li是第i声道长度，tdi是第i声道上超声波的顺流渡越时间，tui是第i声道上超声波的逆流渡越时间。在声道长度一定的情况下，各声道的测量声速只与该声道上超声波信号在顺、逆流中的渡越时间有关。在多声道超声波气体流量计正常运行的过程中，单一声道上测量得到的测量声速不会有明显的改变。当某一声道出现故障后，该声道上测量得到的顺、逆流渡越时间将偏离正常值，这将导致该声道上的测量声速将偏离正常范围。当单一声道测量声速差大于0.5m/s时，可判定该声道出现故障。
3)通过判断测量声速和理论计算声速的偏差是否在允许的范围内来监测各个声道的工作状态。工况条件下的理论计算声速为其中，T是实际测量中气体的热力学温度。当某一声道出现故障后，该声道上测量得到的顺、逆流渡越时间将偏离正常值，这将导致测量声速与理论计算声速之间出现较大偏差。通过将理论计算声速和同一条件下多声道超声波气体流量计的各声道测量声速进行比较，当某一声道测量声速和理论计算声速的相对偏差大于0.5％时，可判定该声道出现故障。
4)预先计算四声道超声波气体流量计的任意三个声道对应的最优权重系数，以调整任意单个声道出现故障后其他正常工作声道对应的权重系数。权重系数调整具体方法如下：在声道A出现故障不能正常测量的情况下，利用其它三个正常测量声道B、C和D的测量数据计算管道内气体的平均流速。通过建立式(1)所示的最小二乘函数，当fA(w)取到最小值时，计算得到的wB、wC和wD是最优的权重系数。此时，利用三个正常工作声道的测量数据加权计算得到的流速和参考流速的误差平方和最小。同理，在声道B、C、D分别出现故障不能正常测量的情况下，其它正常工作声道对应的加权权重系数可以分别通过建立如式(2)-式(4)所示的最小二乘函数计算得到。式(1)-式(4)中，vAi、vBi、vCi、vDi分别是声道A、B、C、D第i次测量的流速，wA、wB、wC、wD分别是声道A、B、C、D对应的加权权重系数，Vi为第i次测量的参考流速，n是测量的总次数。式(1)中wB、wC、wD满足wB+wC+wD＝1，式(2)中wA、wC、wD满足wA+wC+wD＝1，式(3)中wA、wB、wD满足wA+wB+wD＝1，式(4)中wA、wB、wC满足wA+wB+wC＝1。
5)当监测到某一声道出现故障不能正常测量时，舍弃该异常声道的测量数据，利用其余正常工作声道的测量数据和对应的调整后的权重系数计算管道内的气体流速，弥补声道出现故障后流速测量的偏差。
优选的，所述的四声道超声波气体流量测量传感器的安装位置由Gauss-Jacobi积分方法确定，即水平面上下各布置两个声道，每个声道上的两个换能器的连线与管道轴线呈一定角度。
本发明具有的有益效果：
1)探索了多声道超声波气体流量计故障诊断及补偿方法，在流量计无法及时更换和维修的情况下仍能保证较高的计量精度，为多声道超声波气体流量计的稳定可靠连续运行提供了一种解决思路。
2)本发明对多声道超声波气体流量计故障诊断及补偿方法进行了实验验证，能取得较好的效果。
附图说明
图1是四声道超声波气体流量计的声道分布左视图。
图2是四声道超声波气体流量计的声道分布俯视图。
图3是单个声道上时差法测量原理图。
具体实施方式
下面结合实施例和说明书附图对本发明作进一步说明。
如图1所示，四声道超声波气体流量计在管道内不同高度的水平截面上布置着四对超声波换能器，四个声道沿管道中心所在水平面对称布置，即水平面上下各布置两个声道。如图2所示，每个声道上的两个换能器的连线与管道轴线呈一定角度。如图3所示，通过每对超声波换能器交替发射和接收超声波信号测得各个声道上超声波信号在顺流中的渡越时间td和在逆流中的渡越时间tu，然后根据时差法计算得到对应声道上的平均流速v。通过依次计算四个声道上的流速测量值v，获取圆形管道内不同流层上的流速信息。
本发明对多声道超声波气体流量计故障诊断及补偿方法进行了研究，具体操作步骤如下：
1)由四声道超声波气体流量计依次获取四个声道上超声波信号在顺流中的渡越时间td和在逆流中的渡越时间tu，利用式(5)分别计算各个声道上的流速测量值v，其中L为声道长度，α为声道与管道轴线之间的夹角。
2)通过判断单一声道上最大测量声速差是否在允许的范围内监测各个声道的工作状态。各个声道上的测量声速由式(6)计算得到，其中，i是声道编号(1-4)，Cfi是第i声道上的测量声速，Li是第i声道长度，tdi是第i声道上超声波的顺流渡越时间，tui是第i声道上超声波的逆流渡越时间。在声道长度一定的情况下，各声道的测量声速只与该声道上超声波信号在顺、逆流中的渡越时间有关。在多声道超声波气体流量计正常运行的过程中，单一声道上测量得到的测量声速不会有明显的改变。当某一声道出现故障后，该声道上测量得到的顺、逆流渡越时间将偏离正常值，这将导致该声道上的测量声速偏离正常范围。当单一声道测量声速差大于0.5m/s时，可判定该声道出现故障。
3)通过判断测量声速和理论计算声速的偏差是否在允许的范围内监测各个声道的工作状态。工况条件下的理论计算声速由式(7)计算得到，其中，T是实际测量中气体的热力学温度。当某一声道出现故障后，该声道上测量得到的顺、逆流渡越时间将偏离正常值，这将导致测量声速与理论计算声速之间出现较大偏差。通过将理论计算声速和同一条件下多声道超声波气体流量计的各声道测量声速进行比较，当某一声道测量声速和理论计算声速的相对偏差大于0.5％时，可判定该声道出现故障。
4)预先计算四声道超声波气体流量计的任意三个声道对应的最优权重系数，以调整任意单个声道出现故障后其他正常工作声道对应的权重系数。权重系数调整具体方法如下：在声道A出现故障不能正常测量的情况下，利用其它三个正常测量声道B、C和D的测量数据计算管道内气体的平均流速。通过建立式(8)所示的最小二乘函数，当fA(w)取到最小值时，计算得到的wB、wC和wD是最优的权重系数。此时，利用三个正常工作声道的测量数据加权计算得到的流速和参考流速的误差平方和最小。同理，在声道B、C、D分别出现故障不能正常测量的情况下，其它正常工作声道对应的加权权重系数可以分别通过建立如式(9)-式(11)所示的最小二乘函数计算得到。式(8)-式(11)中，vAi、vBi、vCi、vDi分别是声道A、B、C、D第i次的流速测量值，wA、wB、wC、wD分别是声道A、B、C、D对应的加权权重系数，Vi为第i次测量的参考流速，n是测量的总次数。式(8)中wB、wC、wD满足wB+wC+wD＝1，式(9)中wA、wC、wD满足wA+wC+wD＝1，式(10)中wA、wB、wD满足wA+wB+wD＝1，式(11)中wA、wB、wC满足wA+wB+wC＝1。
5)当监测到某一声道出现故障不能正常测量时，舍弃该异常声道的测量数据，利用其余正常工作声道的测量数据和对应的调整后的权重系数计算管道内的气体流速，弥补声道出现故障后流速测量的偏差。当声道A、B、C、D分别发生故障时，管道内平均流速计算公式分别如式(12)-(15)所示，其中vA、vB、vC、vD分别是声道A、B、C、D上的流速测量值，wA、wB、wC、wD分别是声道A、B、C、D上流速的加权权重系数。
VA＝wBvB+wCvC+wDvD (12)
VB＝wAvA+wCvC+wDvD (13)
VC＝wAvA+wBvB+wDvD (14)
VD＝wAvA+wBvB+wCvC (15)
已利用四声道超声波气体流量计在300mm的管径中进行了气体实流实验，在0m/s～8m/s流速范围内选取了9个流速点，每个流速点重复测量300次。实验结果显示，利用本发明中所提到的方法，在某一声道出现声道故障问题后，补偿后的流速测量平均值和四个声道正常工作下流速测量平均值的相对误差小于1.2％。这表明本发明提出的声道故障诊断及补偿方法是有效的。
Claims (5)
1.一种多声道超声波气体流量计故障诊断及补偿方法，其特征在于包括以下步骤：
1)四声道超声波气体流量计在管道内不同高度的水平截面上布置着四对超声波换能器，通过每对超声波换能器交替发射和接收超声波信号测得各个声道上超声波信号在顺、逆流中的渡越时间，然后根据时差法计算得到对应声道上的平均流速；
2)通过判断单一声道上最大测量声速差是否在允许的范围内监测各个声道的工作状态，判断声道是否出现故障；
3)通过将理论计算声速和同一条件下各声道的测量声速进行比较，根据各个声道上测量声速和理论计算声速之间的偏差是否在允许的范围内，判断声道是否出现故障，当故障发生时定位出现故障的声道；
4)预先计算四声道超声波气体流量计的任意三个声道对应的最优权重系数；
5)当某一声道出现故障时，利用其余正常工作声道的测量数据和对应的三声道最优权重系数计算管道内的气体流速，弥补声道出现故障后流速测量的偏差，管道内的气体流速计算公式如式(1)-式(4)所示：
VA＝wBvB+wCvC+wDvD (1)
VB＝wAvA+wCvC+wDvD (2)
VC＝wAvA+wBvB+wDvD (3)
VD＝wAvA+wBvB+wCvC (4)
其中，VA、VB、VC、VD分别为声道A、B、C、D发生故障时补偿后的管道内的气体流速，vA、vB、vC、vD分别是声道A、B、C、D上的流速测量值，wA、wB、wC、wD分别是声道A、B、C、D上流速的加权权重系数。
2.根据权利要求1所述的多声道超声波气体流量计故障诊断及补偿方法，其特征在于所述的四对超声波换能器的安装位置由Gauss-Jacobi积分方法确定。
3.根据权利要求1所述的多声道超声波气体流量计故障诊断及补偿方法，其特征在于所述的步骤2)具体为：
在声道长度一定的情况下，各声道的测量声速只与该声道上超声波信号在顺、逆流中的渡越时间有关，在多声道超声波气体流量计正常运行的过程中，单一声道上测得的测量声速不会有明显的改变，当某一声道出现故障后，该声道上测量得到的顺、逆流渡越时间将偏离正常值，这将导致该声道上的测量声速偏离正常范围，当单一声道测量声速差超过允许的范围时，可判定该声道测量出现故障。
4.根据权利要求1所述的多声道超声波气体流量计故障诊断及补偿方法，其特征在于所述的步骤3)具体为：
当某一声道出现故障后，该声道上测量得到的顺、逆流渡越时间将偏离正常值，这将导致测量声速与理论计算声速之间出现较大偏差，通过将理论计算声速和同一条件下多声道超声波气体流量计的各声道测量声速进行比较，根据各个声道上测量声速和理论计算声速之间的偏差是否在允许的范围内可以准确定位故障声道的位置；
其中，理论计算声速为T是实际测量中气体的热力学温度；测量声速为其中，i是声道编号1-4，Cfi是第i声道上的测量声速，Li是第i声道长度，tdi是第i声道上超声波的顺流渡越时间，tui是第i声道上超声波的逆流渡越时间。
5.根据权利要求1所述的多声道超声波气体流量计故障诊断及补偿方法，其特征在于所述的步骤4)具体为：
在声道A出现故障不能正常测量的情况下，利用其他三个正常测量声道B、C和D的测量数据计算管道内气体的平均流速，通过建立式(5)所示的最小二乘函数，当fA(w)取到最小值时，计算得到的wB、wC和wD是最优的权重系数，此时，利用三个正常工作声道的测量数据计算得到的流速和参考流速的误差平方和最小；同理，在声道B、C、D分别出现故障不能正常测量的情况下，其它正常工作声道对应的最优的权重系数分别可以通过建立如式(6)-式(8)所示的最小二乘函数计算得到，在式(5)-式(8)中，vAi、vBi、vCi、vDi分别是声道A、B、C、D第i次测量的流速，wA、wB、wC、wD分别是声道A、B、C、D对应的加权权重系数，Vi为第i次测量的参考流速，n是测量的总次数，式(5)中wB、wC、wD满足wB+wC+wD＝1，式(6)中wA、wC、wD满足wA+wC+wD＝1，式(7)中wA、wB、wD满足wA+wB+wD＝1，式(8)中wA、wB、wC满足wA+wB+wC＝1，
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610988681.6A CN106643987B (zh) | 2016-11-10 | 2016-11-10 | 一种多声道超声波气体流量计故障诊断及补偿方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610988681.6A CN106643987B (zh) | 2016-11-10 | 2016-11-10 | 一种多声道超声波气体流量计故障诊断及补偿方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106643987A CN106643987A (zh) | 2017-05-10 |
| CN106643987B true CN106643987B (zh) | 2019-02-22 |
Family
ID=58805672
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610988681.6A Expired - Fee Related CN106643987B (zh) | 2016-11-10 | 2016-11-10 | 一种多声道超声波气体流量计故障诊断及补偿方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106643987B (zh) |
Families Citing this family (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107367306A (zh) * | 2017-08-31 | 2017-11-21 | 合肥工业大学 | 一种基于Kalman滤波的气体超声波流量计融合方法和系统 |
| CN107607158A (zh) * | 2017-09-18 | 2018-01-19 | 李言钦 | 声学法测量管内轴向流速分布及流量的系统 |
| CN108120481B (zh) * | 2017-11-10 | 2023-10-31 | 天津新科成套仪表有限公司 | 一种超声流量计量方法与计量处理装置 |
| CN107907172A (zh) * | 2017-11-10 | 2018-04-13 | 陈兵 | 一种超声流量计量监控方法与系统 |
| CN111896080B (zh) * | 2020-07-31 | 2022-08-12 | 中国重汽集团济南动力有限公司 | 一种尿素质量传感器液位测量偏差诊断的计算方法及装置 |
| CN112033581B (zh) * | 2020-08-25 | 2022-07-22 | 博锐尚格科技股份有限公司 | 一种冷热量表在线故障诊断方法及装置 |
| CN113702662B (zh) * | 2021-08-30 | 2024-01-12 | 河北大学 | 一种多声道同步发射流速测量装置及方法 |
| CN114136387B (zh) * | 2021-11-25 | 2022-12-20 | 北京化工大学 | 基于支持向量机模型svm算法的多声道超声流量计误差补偿方法 |
| CN114689104B (zh) * | 2022-03-23 | 2024-01-12 | 天津水泥工业设计研究院有限公司 | 一种大型烟道用超声波传感设备自校准系统及方法 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7942068B2 (en) * | 2009-03-11 | 2011-05-17 | Ge Infrastructure Sensing, Inc. | Method and system for multi-path ultrasonic flow rate measurement |
| EP2386835B1 (de) * | 2010-05-12 | 2015-11-25 | SICK Engineering GmbH | Ultraschallmessung der Strömungsgeschwindigkeit eines Fluids in einer Rohrleitung |
| EP2428777B1 (en) * | 2010-09-09 | 2014-04-30 | SICK Engineering GmbH | A method of operating an ultrasonic gas flow meter and an ultrasonic gas flow meter |
| CN103074873B (zh) * | 2013-02-07 | 2014-12-24 | 天津大学 | 一种非理想流场下的多声道超声流量计声道布置方法 |
| CN103196504B (zh) * | 2013-03-21 | 2015-01-14 | 浙江大学 | 一种多声道超声波流量测量方法及装置 |
| CN105091970B (zh) * | 2015-08-07 | 2018-06-01 | 北京控制工程研究所 | 超声波流量计动态补偿方法 |
| CN106092229B (zh) * | 2016-06-16 | 2018-12-04 | 浙江大学 | 针对超声波气体流量计的有效信号盲提取方法和装置 |
-
2016
- 2016-11-10 CN CN201610988681.6A patent/CN106643987B/zh not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN106643987A (zh) | 2017-05-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106643987B (zh) | 一种多声道超声波气体流量计故障诊断及补偿方法 | |
| RU2594951C2 (ru) | Измеритель скорости потока, работающий по принципу дифференцированного давления, с резервными датчиками давления, позволяющими обнаружить отказ датчиков и снижение производительности | |
| CN103196504B (zh) | 一种多声道超声波流量测量方法及装置 | |
| CN106979808B (zh) | 一种超声与靶式流量计组合式湿天然气流量测量方法 | |
| CN108506740B (zh) | 一种基于流速计的液体管道泄漏区域确定方法及系统 | |
| CN101246035A (zh) | 便携式天然气流量计量核查系统 | |
| CN113551744A (zh) | 一种超声流量计性能在线监测方法及系统 | |
| US10775214B2 (en) | Insertion type ultrasonic flow meter, flow measuring system and method | |
| CN102967349A (zh) | 一种气体超声流量计的网络化在线检验方法 | |
| CN208296916U (zh) | 一种基于超声波时差法的多声道流量测量装置 | |
| CN106323392A (zh) | 抗干扰的超声流量传感器的测量管段 | |
| CN102207398A (zh) | 用于燃油终端结算的超声波流量测量的装置和方法 | |
| AU2015286219A1 (en) | Measurement of flow through pipelines | |
| CN202158878U (zh) | 用于燃油终端结算的超声波流量测量的装置 | |
| CN111059477B (zh) | 一种基于双层架构反向寻源化工管道泄漏检测与定位方法 | |
| CN102288265B (zh) | 双通道法超声流量计时间差检测装置 | |
| CN209558048U (zh) | 一种基于分布式次声波传感器的管道泄漏定位系统 | |
| CN102095889B (zh) | 三通道超声时差流速测量方法 | |
| US11698281B2 (en) | Real-time measurement of two-phase mass flow rate and enthalpy using pressure differential devices | |
| Zheng et al. | Research on flow diagnosis of multipath ultrasonic flowmeter | |
| CN101650212B (zh) | 一种流速仪法流量测量数据准确度评定方法 | |
| CN114563065A (zh) | 气体超声波流量计校准方法 | |
| JPH0961283A (ja) | 配管漏洩監視装置 | |
| CN107621293A (zh) | 地下用高精度气体超声流量测量装置及测量方法 | |
| CN201859001U (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 |
Granted publication date: 20190222 Termination date: 20201110 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |