CN103257114A - 基于光纤Bragg光栅传感网络的井下瓦斯检测方法 - Google Patents
基于光纤Bragg光栅传感网络的井下瓦斯检测方法 Download PDFInfo
- Publication number
- CN103257114A CN103257114A CN2013101676888A CN201310167688A CN103257114A CN 103257114 A CN103257114 A CN 103257114A CN 2013101676888 A CN2013101676888 A CN 2013101676888A CN 201310167688 A CN201310167688 A CN 201310167688A CN 103257114 A CN103257114 A CN 103257114A
- Authority
- CN
- China
- Prior art keywords
- bragg grating
- optical fiber
- microprocessor
- fiber bragg
- light
- 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
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种基于光纤Bragg光栅传感网络的井下瓦斯检测方法,包括有光纤检测模块、数据采集中心、监控中心,光纤检测模块包括有光源、光纤Bragg光栅传感器,数据采集中心包括有依次连接的光电探测器、信号放大器、AD转换器,监控中心包括有依次连接的微处理器、计算机,光纤Bragg光栅传感器连接光电探测器,AD转换器连接微处理器;通过监测点的光信号被光纤Bragg光栅传感器的Bragg光栅滤波分成两束后被光电探测器分别接收,进行放大滤波后的数字信号量采用去极值数字滤波法进行数据处理转换后由微处理器分析处理并计算出瓦斯浓度。本发明很大程度上提升了系统的准确性和稳定性,确保监控系统的连续性、可靠性及有效性,最终实现“监控预防”的目的。
Description
技术领域
本发明涉及到井下瓦斯监控系统的传感技术领域,具体为一种基于光纤Bragg光栅传感网络的井下瓦斯检测方法。
背景技术
信息技术数据的采集主要依靠传感器和探测器进行,因此新型的传感技术在数据采集中尤其是在恶劣的井下环境中采集而来的数据在井下瓦斯监控系统中就显的尤为突出,作为煤矿安全监控系统中一项重要的指标,传感器网络自井下采集来的数据直接关系到监测系统的监控效率,所以传感器的稳定性、可靠性、以及传感器采集来的数据在瓦斯监控整个系统中就起到至关重要的作用。目前国内开发的的传感器虽然基本满足了井下瓦斯抽采系统的需求,但是面对恶劣的井下工作环境,现有的传感器无论是在敏感元件的质量上还是在传感器测量数据的准确性上都有待提高,长期工作稳定性差,调校频繁、复杂以及工作寿命短已经成为现有传感器无法回避的问题,所以本文提出一种光纤Bragg传感器代替传统的瓦斯传感器以应对现有系统对传感器使用寿命、稳定性、调节周期,精确度等提出的更高要求的现状。
为了保证煤矿安全监控系统的有效性和稳定性,研制开发一种基于光纤Bragg光栅传感网络的新型井下瓦斯监控系统很有必要,这样不仅可以保证井下数据的精确性,监测系统的稳定连续性,煤矿安全监控系统的正常运行,而且可以最大程度上避免因事故而带来的损失。
发明内容
本发明要解决的技术问题是针对现有矿井安全监控系统中传感系统的不足,提供一种采用差分吸收检测原理的基于光纤Bragg光栅传感网络的井下瓦斯检测方法,以解决现有技术检测方法精度低,易受干扰,实时性差等的问题。
为了达到上述目的,本发明所采用的技术方案为:
基于光纤Bragg光栅传感网络的井下瓦斯检测方法,其特征在于:包括有瓦斯监测系统,所述的瓦斯监测系统包括有光纤检测模块、数据采集中心、监控中心,所述的光纤检测模块包括有光源、光纤Bragg光栅传感器,数据采集中心包括有依次连接的光电探测器、信号放大器、AD转换器,监控中心包括有依次连接的微处理器、计算机,所述的光纤Bragg光栅传感器的信号输出端连接光电探测器,AD转换器的信号输出端连接微处理器,微处理器还连接有数据存储器、声光报警器,计算机还连接有LCD显示模块;通过监测点的光信号被光纤Bragg光栅传感器的Bragg光栅滤波分成两束后被光电探测器分别接收,进行放大滤波后的数字信号量采用去极值数字滤波法进行数据处理转换后由微处理器分析处理并计算出瓦斯浓度,具体包括以下步骤:
(1)根据待测区域不同的地域情况,采用多个光纤Bragg光栅传感器分别检测不同待检测区域的待测参数,待测参数以光信号为载体返回控制室;
(2)所述监测点的光信号返回控制室被光纤Bragg光栅传感器的Bragg光栅滤波分成两束被光电探测器接收,电信号经过信号放大器处理后再经AD转换器转换为数字信号送入微处理器进行处理;
(3)在步骤(2)中所述微处理器内部操作系统中创建多个系统任务,用来保证系统各个模块功能的连续进行,所述的微处理器包含1个16通道差分输入多路选择器和1个10位SAR ADC并对接收来自A/D采样转换后的数字量采用去极值数字滤波法做了数据处理后,采用差分吸收法检测计算分析出结果,所述待检测区域的瓦斯浓度通过如下公式计算出来:
有比尔-朗伯定律:
I(λ)=K(λ)I0(λ)exp[-α(λ)LC+β(λ)+γ(λ)L+δ(λ)]
在此公式中,I(λ)为透射光的光强;I0(λ)为入射光的光强;α(λ)为待测气体的吸收系数;K(λ)为整个系统的光效率;β(λ)为瑞利散射系数;γ(λ)为米氏散射系数;C为待测气体浓度;L为待测气体的整个吸收光程;δ(λ)为气体密度造成波动的吸收系数;
波长相差较小分别为λ1、λ2且在吸收系数上有差别,对应待测气体的参考和工作波长,得到瓦斯气体浓度C:
所述的基于光纤Bragg光栅传感网络的井下瓦斯检测方法,其特征在于:所述的光源采用InGaAsP/InP超辐射的发光二极管所提供的红外光波。
所述的基于光纤Bragg光栅传感网络的井下瓦斯检测方法,其特征在于:所述的光纤Bragg光栅传感器通过选定的用户分配各自唯一的地址码,同时进行扩频处理,数据接收时可以利用相关运算恢复初始指定信号并对井下多接触面多点进行分布式的定点测量。
所述的基于光纤Bragg光栅传感网络的井下瓦斯检测方法,其特征在于:所述的微处理器采用C8051F330。
本发明的有益效果在于:
本发明可以在井下复杂多变的环境下仍然可以在最大程度上保证监测系统的稳定性和监测数据的准确性,在面对各种干扰的情况下,基于光纤Bragg光栅传感网络瓦斯检测方法同时把编码技术和传感技术相结合并融合于一体,基本实现了在待测地域监测点的多点分布式测量,本系统同时还对环境中的例如温度、暗电流的造成的噪声等变化均有较强的抗干扰能力,很大程度上提升了系统的准确性和稳定性,确保监控系统的的连续性、可靠性及有效性,最终实现“监控预防”的目的。
附图说明
图1为本发明的系统图。
图2为本发明微处理器C8051F330建立的主程序流程图。
具体实施方式
如图1所示。基于光纤Bragg光栅传感网络的井下瓦斯检测方法,包括有瓦斯监测系统,瓦斯监测系统包括有光纤检测模块1、数据采集中心2、监控中心3,光纤检测模块1包括有光源4、光纤Bragg光栅传感器5,数据采集中心2包括有依次连接的光电探测器6、信号放大器7、AD转换器8,监控中心3包括有依次连接的微处理器9、计算机10,光纤Bragg光栅传感器5的信号输出端连接光电探测器6,AD转换器8的信号输出端连接微处理器9,微处理器9还连接有数据存储器11、声光报警器12,计算机10还连接有LCD显示模块13;通过监测点的光信号被光纤Bragg光栅传感器5的Bragg光栅滤波分成两束后被光电探测器分别接收,进行放大滤波后的数字信号量采用去极值数字滤波法进行数据处理转换后由微处理器分析处理并计算出瓦斯浓度,具体包括以下步骤:
(1)根据待测区域不同的地域情况,采用多个光纤Bragg光栅传感器5分别检测不同待检测区域的待测参数,待测参数以光信号为载体返回控制室;
(2)所述监测点的光信号返回控制室被光纤Bragg光栅传感器5的Bragg光栅滤波分成两束被光电探测器6接收,电信号经过信号放大器7处理后再经AD转换器8转换为数字信号送入微处理器9进行处理;
(3)在步骤(2)中所述微处理器9内部操作系统中创建多个系统任务,用来保证系统各个模块功能的连续进行,所述的微处理器9包含1个16通道差分输入多路选择器和1个10位SAR ADC并对接收来自A/D采样转换后的数字量采用去极值数字滤波法做了数据处理后,采用差分吸收法检测计算分析出结果,所述待检测区域的瓦斯浓度通过如下公式计算出来:
有比尔-朗伯定律:
I(λ)=K(λ)I0(λ)exp[-α(λ)LC+β(λ)+γ(λ)L+δ(λ)]
在此公式中,I(λ)为透射光的光强;I0(λ)为入射光的光强;α(λ)为待测气体的吸收系数;K(λ)为整个系统的光效率;β(λ)为瑞利散射系数;γ(λ)为米氏散射系数;C为待测气体浓度;L为待测气体的整个吸收光程;δ(λ)为气体密度造成波动的吸收系数;
波长相差较小分别为λ1、λ2且在吸收系数上有差别,对应待测气体的参考和工作波长,得到瓦斯气体浓度C:
本发明采用美国Cygnal公司的C8051F330单片机作为系统的控制核心,一方面根据不同的待测区域不同的地域情况,采用多个光纤Bragg光栅传感器分别用来检测待检测区域的待测参数,来自InGaAsP/InP超辐射的发光二极管的红外光波经过待测地域后被接收转换为电信号,对电信号进行处理后给微处理器。另一方面微处理器先采用去极值数字滤波方法对采样转换来的数据进行处理,同时微处理器通过用户分配唯一的地址码和对的Bragg传感器进行扩频处理,以保证系统资源的合理利用以及增多传感器个数以加大区域监测范围,还可以确保当数据接收中断时可以利用相关运算恢复初始信号,监控主机采用以太网完成和数据库服务器数据之间的通讯;远程客户终端的WEB浏览服务由本地数据库的服务器提供,以便随时随地了解当前系统监控信息。
1. 差分吸收检测原理
当一束光通过某种介质时,就会有光的一部分会被这种介质吸收,少量的光波通过介质会被介质散射,其余的光波则沿着原来的方向继续的传播,所以有比尔-朗伯定律可知,介质的吸收长度、透射光的光强、光的吸收系数与物质浓度之间满足:
I=I0exp[-αLC+β+γL+δ]
式中,I,I0分别为透射和入射光强;L为待测气体的吸收光程、;C为待测气体的浓度;α为某单位波长下的浓度、单位长度介质的吸收系数;β为瑞利散射系数;δ为气体密度波动造成的吸收系数;γ为米氏散射系数。用光来检测气体时,应该还要包括比例系数K,为此公式可改写为
I=KI0exp[-αLC+β+γL+δ]
但是仅利用上式确定待测气体的浓度C时还是不足的,因为δ的反应的为平均变化数,随着时间的变化而不停的发生变化,所以它是随机量。所以用两个吸收系数有较大差距但是波长相近的两束波长为λ1、λ2的单色光在极短的时间通过将要测量的气体,有:
I(λ1)=K(λ1)I0(λ1)exp[-α(λ1)LC+β(λ1)+γ(λ1)L+δ(λ1)]
I(λ2)=K(λ2)I0(λ2)exp[-α(λ2)LC+β(λ2)+γ(λ2)L+δ(λ2)]
气体浓度关系表达式为:
由于λ1、λ2数值相差较小,同时单色光也几乎是同一时间接近然后通过待测 气体的,则:
β(λ1)≈β(λ2),γ(λ1)≈γ(λ2),δ(λ1)≈δ(λ2),
所以
对光学系统内部进行调节使的:
K(λ1)I0(λ1)=K(λ2)I0(λ2)
可得
有光学上的相近条件可知:
[α(λ1)-α(λ2)]L≤1
在实际监控过程中,波长为λ1、λ2分别对应着不同吸收谱线,波长为λ2的光作为参考波长不能被待测气体吸收,检测过程中为空值,因此I(λ1)<I(λ2),I(λ2)/I(λ1)>1和I(λ1)/I(λ2)<1。对ln[I(λ2)/I(λ1)]进行泰勒展开:
所以气体的浓度为
在波长λ1、λ2、α1、α2已知情况下就可以通过I(λ2)-I(λ1)和I(λ2)关系式来测量待测气体浓度。
1)去极值的数字滤波方法
系统软件主要用来控制数据的分析与处理以及数据的保存,分析计算出的结果的显示以及需要保存的数据经由RS-232接口传递给上位机等功能更加强大的模块进行。启动脉冲来临前,首先进行初始化处理,对A/D、定时存储器等完成前期的初始化工作然后等待中断请求的到来,为了提高检测精度和消除随机噪声干扰,采用去极值的数字滤波方法对A/D采样转换来的数据进行处理。将事先标定好的线性系数存储在系统的ROM存储器中,当需要做数据处理时,根据当前分析所得数据获得相关的系数,求取当前瓦斯气体的浓度数据送显示器显示,当待测企图浓度超过设定的限定值时,处理器发出声光报警。软件部分采用模块化设计,程序的主流程图如图2所示
为了达到提高待测气体浓度数据精度的目的,对经过A/D采样转换后的数字量采用去极值数字滤波法进行数据处理。此方法首先需要建立一个具有一定容量的窗口平台,窗口内包含的数据量为n。将这n个数据放到一个数组当中送到微处理器进行循环分析处理得到n个结果,将n个结果进行比较获得该组检测数据当中待测值最小值与最大值,之后再计算出所有数据结果的累加和,在这个累加和的数值当中减去循环得到的最大值和最小值,这样的累加和的数值就组成变成了由原来的n个数据变成了n-2个数据,把得到n-2个数据进行累加再除以n-2得到的均值系统就当做一个监测数据。在一组数据的累加和中去除去最小值和最大值可以有效的消除一组监测数据中可能出现的变化较大的随机信号变量对最终测量结果的影响,从而达到了消除脉冲干扰,提高精确率的良好效果。根据误差理论的原则,相同的物理量在等同等精度测量初始条件下条件下,当测量次数越大时,其结果所得到的均值就越逼近真实值,当测量次数>10次之后,测量精度的提升将会放缓,这样就可以在确保检测精度和满足实时检测需要的这两个标准的要求下,可以采用选择连续采样12次的检测数据利用去极值数字滤波法处理后得到的数据做为当前瓦斯浓度检测的最终结果。
Claims (4)
1.基于光纤Bragg光栅传感网络的井下瓦斯检测方法,其特征在于:包括有瓦斯监测系统,所述的瓦斯监测系统包括有光纤检测模块、数据采集中心、监控中心,所述的光纤检测模块包括有光源、光纤Bragg光栅传感器,数据采集中心包括有依次连接的光电探测器、信号放大器、AD转换器,监控中心包括有依次连接的微处理器、计算机,所述的光纤Bragg光栅传感器的信号输出端连接光电探测器,AD转换器的信号输出端连接微处理器,微处理器还连接有数据存储器、声光报警器,计算机还连接有LCD显示模块;通过监测点的光信号被光纤Bragg光栅传感器的Bragg光栅滤波分成两束后被光电探测器分别接收,进行放大滤波后的数字信号量采用去极值数字滤波法进行数据处理转换后由微处理器分析处理并计算出瓦斯浓度,具体包括以下步骤:
(1)根据待测区域不同的地域情况,采用多个光纤Bragg光栅传感器分别检测不同待检测区域的待测参数,待测参数以光信号为载体返回控制室;
(2)所述监测点的光信号返回控制室被光纤Bragg光栅传感器的Bragg光栅滤波分成两束被光电探测器接收,电信号经过信号放大器处理后再经AD转换器转换为数字信号送入微处理器进行处理;
(3)在步骤(2)中所述微处理器内部操作系统中创建多个系统任务,用来保证系统各个模块功能的连续进行,所述的微处理器包含1个16通道差分输入多路选择器和1个10位SAR ADC并对接收来自A/D采样转换后的数字量采用去极值数字滤波法做了数据处理后,采用差分吸收法检测计算分析出结果,所述待检测区域的瓦斯浓度通过如下公式计算出来:
有比尔-朗伯定律:
I(λ)=K(λ)I0(λ)exp[-α(λ)LC+β(λ)+γ(λ)+δ(λ)]
在此公式中,I(λ)为透射光的光强;I0(λ)为入射光的光强;α(λ)为待测气体的吸收系数;K(λ)为整个系统的光效率;β(λ)为瑞利散射系数;γ(λ)为米氏散射系数;C为待测气体浓度;L为待测气体的整个吸收光程;δ(λ)为气体密度造成波动的吸收系数;
波长相差较小分别为λ1、λ2且在吸收系数上有差别,对应待测气体的参考和工作波长,得到瓦斯气体浓度C:
2.根据权利要求1所述的基于光纤Bragg光栅传感网络的井下瓦斯检测方法,其特征在于:所述的光源采用InGaAsP/InP超辐射的发光二极管所提供的红外光波。
3.根据权利要求1所述的基于光纤Bragg光栅传感网络的井下瓦斯检测方法,其特征在于:所述的光纤Bragg光栅传感器通过选定的用户分配各自唯一的地址码,同时进行扩频处理,数据接收时可以利用相关运算恢复初始指定信号并对井下多接触面多点进行分布式的定点测量。
4.根据权利要求1所述的基于光纤Bragg光栅传感网络的井下瓦斯检测方法,其特征在于:所述的微处理器采用C8051F330。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310167688.8A CN103257114B (zh) | 2013-05-08 | 2013-05-08 | 基于光纤Bragg光栅传感网络的井下瓦斯检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310167688.8A CN103257114B (zh) | 2013-05-08 | 2013-05-08 | 基于光纤Bragg光栅传感网络的井下瓦斯检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103257114A true CN103257114A (zh) | 2013-08-21 |
CN103257114B CN103257114B (zh) | 2015-04-22 |
Family
ID=48961161
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310167688.8A Expired - Fee Related CN103257114B (zh) | 2013-05-08 | 2013-05-08 | 基于光纤Bragg光栅传感网络的井下瓦斯检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103257114B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103926209A (zh) * | 2013-12-26 | 2014-07-16 | 安徽理工大学 | 一种基于doas双采样技术的煤层瓦斯联测系统 |
CN111380670A (zh) * | 2020-02-28 | 2020-07-07 | 浙江工业大学 | 一种体状及微棱镜镜片的棱镜度智能检测系统 |
CN111397852A (zh) * | 2020-02-28 | 2020-07-10 | 浙江工业大学 | 一种体状及微棱镜镜片的棱镜度智能检测方法 |
CN114200875A (zh) * | 2020-07-02 | 2022-03-18 | 广东电网有限责任公司江门供电局 | 一种超辐射双光源外破隐患巡查系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101030822A (zh) * | 2006-02-27 | 2007-09-05 | 华为技术有限公司 | 光网络单元和光线路终端及无源光纤接入网及其传输方法 |
CN102207458A (zh) * | 2011-04-08 | 2011-10-05 | 东北大学 | 基于光纤光栅的近红外吸收浓度传感器信号解调方法 |
-
2013
- 2013-05-08 CN CN201310167688.8A patent/CN103257114B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101030822A (zh) * | 2006-02-27 | 2007-09-05 | 华为技术有限公司 | 光网络单元和光线路终端及无源光纤接入网及其传输方法 |
CN102207458A (zh) * | 2011-04-08 | 2011-10-05 | 东北大学 | 基于光纤光栅的近红外吸收浓度传感器信号解调方法 |
Non-Patent Citations (2)
Title |
---|
周孟然等: "差分光学吸收光谱法精确测量煤矿瓦斯浓度", 《煤炭科学技术》, vol. 32, no. 3, 31 March 2004 (2004-03-31) * |
陈伟: "基于光纤布拉格光栅传感技术的瓦斯监测系统设计", 《煤矿安全》, vol. 42, no. 7, 31 July 2011 (2011-07-31) * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103926209A (zh) * | 2013-12-26 | 2014-07-16 | 安徽理工大学 | 一种基于doas双采样技术的煤层瓦斯联测系统 |
CN111380670A (zh) * | 2020-02-28 | 2020-07-07 | 浙江工业大学 | 一种体状及微棱镜镜片的棱镜度智能检测系统 |
CN111397852A (zh) * | 2020-02-28 | 2020-07-10 | 浙江工业大学 | 一种体状及微棱镜镜片的棱镜度智能检测方法 |
CN114200875A (zh) * | 2020-07-02 | 2022-03-18 | 广东电网有限责任公司江门供电局 | 一种超辐射双光源外破隐患巡查系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103257114B (zh) | 2015-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN205374298U (zh) | 基于tdlas的痕量气体浓度检测装置 | |
CN102052930B (zh) | 光纤光栅分布式应变传感器及其应变监测方法 | |
CN102292621B (zh) | 分布式光纤感测中的改进 | |
CN104343466B (zh) | 一种全光纤煤矿安全监测系统 | |
CN103439232A (zh) | 一种减光式烟尘颗粒浓度测量方法及其装置 | |
CN103196465B (zh) | 一种相敏光时域反射仪传感信号噪声分离及信号提取方法 | |
CN101625247B (zh) | 基于dsp的大量程高速光纤光栅传感解调装置与解调方法 | |
CN102914321A (zh) | 一种极弱光纤光栅传感系统及其查询方法 | |
CN102798411A (zh) | 基于布里渊散射的分布式光纤传感测量系统及测量方法 | |
CN103257114A (zh) | 基于光纤Bragg光栅传感网络的井下瓦斯检测方法 | |
CN103674916B (zh) | 新型荧光信号解调装置以及荧光信号解调方法 | |
CN203275291U (zh) | 基于光纤传感技术的煤矿瓦斯气体安全监测系统 | |
CN105136909B (zh) | 一种基于阵列波导光栅的多通道声发射传感解调系统 | |
CN103439233A (zh) | 烟道粉尘浓度检测系统 | |
CA2661092A1 (en) | Optical device for measuring a physical parameter in a hydrogen contaminated sensing zone | |
CN101619807A (zh) | 机载天然气管道泄漏监测系统及监测方法 | |
CN102928003A (zh) | 一种具有实时参考的光纤光栅解调系统 | |
CN105841723A (zh) | 基于小波奇异性检测的fbg传感解调寻峰方法 | |
CN202453113U (zh) | 一种基于fbg的测量海水温度和盐度的传感器 | |
CN102270265B (zh) | 一种混合料水分在线检测方法 | |
CN101949825A (zh) | 光开放环境下的叶片水分近红外无损检测装置及方法 | |
CN204302180U (zh) | 一种痕量气体传感器光路系统及气室 | |
CN111982189B (zh) | 基于光子探测技术同时实现温度与应变的高精度传感方法 | |
CN103558177B (zh) | 生肉组织水分近红外检测装置及检测方法 | |
CN102279096A (zh) | 激光器寿命在线测试方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150422 Termination date: 20160508 |
|
CF01 | Termination of patent right due to non-payment of annual fee |