CN103698333A - 一种自动定量分析煤岩光片裂隙的系统和方法 - Google Patents

一种自动定量分析煤岩光片裂隙的系统和方法 Download PDF

Info

Publication number
CN103698333A
CN103698333A CN201410029309.3A CN201410029309A CN103698333A CN 103698333 A CN103698333 A CN 103698333A CN 201410029309 A CN201410029309 A CN 201410029309A CN 103698333 A CN103698333 A CN 103698333A
Authority
CN
China
Prior art keywords
image
platform
crack
microscope
coal petrography
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
CN201410029309.3A
Other languages
English (en)
Other versions
CN103698333B (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.)
China University of Geosciences Beijing
Original Assignee
China University of Geosciences Beijing
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 China University of Geosciences Beijing filed Critical China University of Geosciences Beijing
Priority to CN201410029309.3A priority Critical patent/CN103698333B/zh
Publication of CN103698333A publication Critical patent/CN103698333A/zh
Application granted granted Critical
Publication of CN103698333B publication Critical patent/CN103698333B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)

Abstract

本发明公开了一种自动定量分析煤岩光片裂隙的系统和方法,涉及煤能源开发技术领域,旨在提供一种操作简便、分析快捷、测试结果更加精确的自动定量分析煤岩光片裂隙的系统和方法,所述系统包括显微镜自动平台控制系统和显微图像分析系统;所述方法的步骤包括:将显微镜自动平台控制系统和显微图像分析系统链接荧光显微镜;调整荧光显微镜;定标设置;图像的全景扫描与采集;图像参数设置;煤岩光片裂隙统计处理分析。本发明是将显微镜自动平台控制系统和显微图像分析系统链接荧光显微镜,这样测试结果将更精确,不存在裂隙重复统计的情况。该方法操作简便,分析快捷,减少了人工耗时,增加了分析项目,剔除了人为误差,使得测试结果更加丰富精确。

Description

一种自动定量分析煤岩光片裂隙的系统和方法
技术领域
本发明涉及能源开发技术领域,特别是一种自动定量分析煤岩光片裂隙的系统和方法。
背景技术
煤中的裂隙常被统称为割理,它是煤层气渗流的主要通道,裂隙的发育程度直接控制煤储层的渗透率,关系到煤层气井产能的大小,因此对煤岩裂隙的定量表征是煤储层物性研究的一个重要组成部分。前人对煤岩裂隙测量统计的主要方法是:首先将煤岩样品抛光制作成规格为30mm×30mm煤岩光片,然后在显微镜下将该煤岩光片划分成10mm×10mm的9个微区,分别将每个微区内出现的裂隙按以下4种类型进行人工统计加权而得。
在长期显微镜观察裂隙的实践工作中,我们发现在反射光显微镜下不能清楚地区别煤岩光片中的裂隙或划道,容易将因制片而产生的划道当成裂隙统计,导致错误结果。同时,这种现有的煤岩裂隙测量统计方法分析工作繁琐,人工耗时量大,镜下观察时需要人工移动样品以切换到下一个视域。而且裂隙统计结果信息量小,分析项目有限,难以满足煤储层物性参数精细定量表征的需要。同时,仪器操作及人为主观判断对测试结果有很大影响,可能导致煤岩裂隙统计结果的人为误差。因此,这种对煤岩裂隙的表征方法存在着很大的局限性和不精确性,亟需改进和替代。
发明内容
为了克服上述缺陷,本发明要解决的技术问题是:提供一种操作简便、分析快捷、减少人为误差、测试结果更加精确的自动定量分析煤岩光片裂隙的系统和方法。
本发明为解决上述技术问题采用的技术方案是:一种自动定量分析煤岩光片裂隙的系统,包括显微镜自动平台控制系统和/或显微图像分析系统;
所述显微镜自动平台控制系统包括电源控制器、显微镜、CCD、计算机和自动载物台控制器,所述电源控制器控制显微镜和计算机,所述显微镜自动载物台控制器控制自动扫描台,通过CCD保存数据并将数据传输给计算机,计算机软件对数据进行处理。
所述显微镜自动平台控制系统包括以下几个区域:工具栏区、图像显示区、平台控制及坐标点区;
其中,所述工具栏区包括图像采集的前期参数设置和定标设置,以及活动图像的打开、冻结、保存和全屏显示;图像显示区显示动态图像及当前保存的静态图像;
所述平台控制包括控制平台的移动方式、移动速度以及移动步距,所述坐标点区能够显示平台的当前状态及当前图像所对应的坐标点;
利用该系统实现样品表面的全景扫描的过程如下:获取起始和终止位置,确定平台的移动速度,开始扫描,扫描的走向为从左到右,即X起始到X终止,从上到下,即Y起始到Y终止,每一步都采集视域下的一个图像,然后将每个小的图像进行拼接,最终得到一个完整的图像。
所述显微图像分析系统包括显微镜下的图像显示、图像处理和图像分析;
所述图像显示包括图像的打开、保存和关闭,以及图像的放大、缩小与定倍显示;
所述图像处理包括亮度、对比度、色度、饱和度、反像显示、图像旋转与阴影校正;
所述图像分析包括图像的灰度直方图的制作、二值化分析、微粒的清除、线的连接和切割以及线条和微粒的手动添加和删除;
所述图像分析系统包括自动扫描台、CCD,所述自动扫描台和CCD获取图像之后,进行图像的处理和分析测量,获得结果及报告及进行加载标尺标注后,进行保存。
通过该系统,能够实现图像的直观显示及分析处理、数据的整理归类和对比分析,最终得到煤岩光片裂隙条数、长度、宽度、面积和周长,达到煤岩光片裂隙自动定量表征的目的。
所述图像显示区,在这个区域上移动鼠标上的滚轮能够调节平台移动的速度,用鼠标进行拖动操作,在平台移动过程中,能单击右键停止平台的移动或暂停连续移动。
一种自动定量分析煤岩光片裂隙的方法,具体步骤如下:
(1)将显微镜自动平台控制系统和显微图像分析系统链接荧光显微镜;
(2)调整荧光显微镜:将煤岩光片粘到载玻片上,压平后,放到载物台上,选择放大倍数,调试荧光显微镜对焦,使荧光显微镜下的图片清晰可辨;
(3)定标设置:打开荧光显微镜自动平台控制系统,进行图像采集的定标设置;
(4)图像的全景扫描与采集:选择图像起点——右上角,使该点置零,通过平台控制移动载物台,选择终点——左下角,然后显示平台框架,在平台框架内进行图像扫描、拼接,最终得到一个完整的图像;
(5)图像参数设置:打开显微图像分析系统,在系统中打开拼接好的完整图像,进行图像的参数设置,调节色彩、对比度、饱和度,使得图像清晰,便于后期图像处理;
(6)煤岩光片裂隙统计处理分析:首先进行图像色彩分割,选择灰度,进行二值化,将煤岩基质和裂隙分割开来;然后清除微粒,去除非裂隙部分——孔隙带来的误差,同时手动删除、或者连接裂隙以降低自动统计的误差;手动测量裂隙,图像上裂隙会显示出不同的颜色,便于区分;最后对光片上的所有裂隙进行汇总统计;
所述煤岩光片的规格为30mm×30mm。
所述煤岩光片放到载物台上时,煤岩光片的一边要与载物台的边缘平行。
本发明的有益效果是:与现有技术相比,本发明是将显微镜自动平台控制系统和显微图像分析系统链接荧光显微镜,这是因为煤岩光片上长短各异的裂隙完全被胶体充填,而胶体具有荧光显示,则在荧光下可以清晰地辨认各类型裂隙。使用显微镜自动平台控制系统控制荧光显微镜载物台的移动方向和距离,完成显微镜下图像的视域自动转换以及图像的全景扫描及采集;同时利用显微图像分析系统将显微镜下的图像显示在计算机显示器上,然后进行处理和分析,最终实现煤岩光片裂隙的自动定量表征,且可分析的裂隙参数增多,统计结果中不仅包括裂隙的条数、长度和宽度,还增加了裂隙周长,面积以及每条裂隙面积所占的比例等参数。针对有些裂隙贯穿几个微区,存在重复统计的可能性,发明首先利用计算机图像全景扫描技术,得到一个完整的煤岩光片图像之后,再进行图像处理。这样测试结果将更加精确,不存在裂隙重复统计的情况。该方法操作简便,分析快捷,减少了人工耗时,增加了分析项目,剔除了人为误差,使得测试结果更加丰富精确。
附图说明
以下结合附图来具体说明本发明所述的一种自动定量分析煤岩光片裂隙的系统。
图1是显微镜自动平台控制系统的结构示意图;
图2是显微图像分析系统的结构示意图。
具体实施方式
本发明采用的原理如下:本发明是将显微镜自动平台控制系统和/或显微图像分析系统链接荧光显微镜,使用显微镜自动平台控制系统控制荧光显微镜载物台的移动方向和距离,完成显微镜下图像的视域自动转换以及图像的全景扫描及采集;同时利用显微图像分析系统将显微镜下的图像显示在计算机显示器上,然后进行处理和分析,最终实现煤岩光片裂隙的自动定量表征。
显微镜自动平台控制系统包括以下几个区域:工具栏区、图像显示区、平台控制及坐标点区;其中,所述工具栏区包括图像采集的前期参数设置和定标设置,以及活动图像的打开、冻结、保存和全屏显示;图像显示区显示动态图像及当前保存的静态图像,在这个区域上移动鼠标上的滚轮能够调节平台移动的速度,用鼠标进行拖动操作,在平台移动过程中,能单击右键停止平台的移动或暂停连续移动;所述平台控制包括控制平台的移动方式、移动速度以及移动步距,所述坐标点区能够显示平台的当前状态及当前图像所对应的坐标点;使用该系统,可以实现样品表面的全景扫描,只要获取起始和终止位置,确定平台的移动速度,然后开始扫描,扫描的走向为从左到右,即X起始到X终止,从上到下,即Y起始到Y终止,每一步都采集视域下的一个图像,然后将每个小的图像进行拼接,最终得到一个完整的图像。
所述显微镜自动平台控制系统包括电源控制器、显微镜、CCD、计算机和自动载物台控制器,所述电源控制器控制显微镜和计算机,所述显微镜自动载物台控制器控制自动扫描台,通过CCD保存数据并将数据传输给计算机,计算机软件对数据进行处理。
显微图像分析系统包括显微镜下的图像显示、图像处理和图像分析;所述图像显示包括图像的打开、保存和关闭,以及图像的放大、缩小与定倍显示;所述图像处理包括亮度、对比度、色度、饱和度、反像显示、图像旋转与阴影校正;所述图像分析包括图像的灰度直方图的制作、二值化分析、微粒的清除、线的连接和切割以及线条和微粒的手动添加和删除;通过该系统,能够实现图像的直观显示及分析处理、数据的整理归类和对比分析,最终得到煤岩光片裂隙条数、长度、宽度、面积和周长,达到煤岩光片裂隙自动定量表征的目的。
所述图像分析系统包括自动扫描台、CCD,所述自动扫描台和CCD获取图像之后,进行图像的处理和分析测量,获得结果及报告及进行加载标尺标注后,进行保存。
利用上述系统进行自动定量分析煤岩光片裂隙的方法,具体步骤如下:
(1)调整荧光显微镜:将规格为30mm×30mm的煤岩光片粘到载玻片上,压平后,放到载物台上,煤岩光片的一边要与载物台的边缘平行。选择放大倍数,调试荧光显微镜对焦,使荧光显微镜下的图片清晰可辨。
(2)定标设置:打开荧光显微镜自动平台控制系统,进行图像采集的定标设置。由于计算机上的图像是通过显微镜光电转换后得到的,通常情况下与显微镜下观察到的图像倍数偏大,这就需要对系统进行校准,以便设置好的定标系数能与对应的显微镜物镜匹配使用。
(3)图像的全景扫描与采集:选择图像起点——右上角,使该点置零,通过平台控制移动载物台,选择终点——左下角,然后显示平台框架,在平台框架内进行图像扫描、拼接,最终得到一个完整的30mm×30mm的图像。
(4)图像参数设置:打开显微图像分析系统,在系统中打开拼接好的完整图像,进行图像的参数设置,调节色彩、对比度、饱和度,使得图像清晰,便于后期图像处理。
(5)煤岩光片裂隙统计处理分析:首先进行图像色彩分割,选择灰度,进行二值化,将煤岩基质和裂隙分割开来;然后清除微粒,去除非裂隙部分——孔隙带来的误差,同时能够手动删除、或者连接一些裂隙以降低自动统计的误差;还能够手动测量典型裂隙,图像上裂隙会显示出不同的颜色,便于区分;最后对光片上的所有裂隙进行汇总统计。
应当理解的是,本发明的上述具体实施方式仅仅用于示例性说明或解释本发明的原理,而不构成对发明的限制。因此,在不偏离本发明的精神和范围的情况下所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。此外,本发明所附权利要求旨在涵盖落入所附权利要求范围和边界、或者这种范围和边界的等同形式内的全部变化和修改例。

Claims (7)

1.一种自动定量分析煤岩光片裂隙的系统,包括显微镜自动平台控制系统和/或显微图像分析系统; 
所述显微镜自动平台控制系统包括电源控制器、显微镜、CCD、计算机和自动载物台控制器,所述电源控制器控制显微镜和计算机,所述显微镜自动载物台控制器控制自动扫描台,通过CCD保存数据并将数据传输给计算机,计算机软件对数据进行处理。 
2.根据权利要求1所述的自动定量分析煤岩光片裂隙的系统,所述显微镜自动平台控制系统包括以下区域:工具栏区、图像显示区、平台控制及坐标点区;所述工具栏区包括图像采集的前期参数设置和定标设置,以及活动图像的打开、冻结、保存和全屏显示;图像显示区显示动态图像及当前保存的静态图像;所述平台控制包括控制平台的移动方式、移动速度以及移动步距,所述坐标点区能够显示平台的当前状态及当前图像所对应的坐标点。 
3.根据权利要求1所述的自动定量分析煤岩光片裂隙的系统,所述显微图像分析系统包括自动扫描台、CCD,所述自动扫描台和CCD获取图像之后,进行图像的处理和分析测量,获得结果及报告及进行加载标尺标注后,进行保存。 
4.根据权利要求1所述的自动定量分析煤岩光片裂隙的系统,所述显微图像分析系统包括显微镜下的图像显示、图像处理和图像分析;所述图像显示包括图像的打开、保存和关闭,以及图像的放大、缩小与定倍显示;所述图像处理包括亮度、对比度、色度、饱和度、反像显示、图像旋转与阴影校正;所述图像分析包括图像的灰度直方图的制作、二值化分析、微粒的清除、线的连接和切割以及线条和微粒的手动添加和删除。 
5.根据权利要求1所述的自动定量分析煤岩光片裂隙的系统,所述图像显示区,在这个区域上移动鼠标上的滚轮调节平台移动的速度,用鼠标进行拖动操作,在平台移动过程中,单击右键停止平台的移动或暂停连续移动。 
6.一种自动定量分析煤岩光片裂隙的方法,具体步骤如下: 
(1)将显微镜自动平台控制系统和显微图像分析系统链接荧光显微镜; 
(2)调整荧光显微镜:将煤岩光片粘到载玻片上,压平后,放到载物台上,选择放大倍数,调试荧光显微镜对焦; 
(3)定标设置:打开荧光显微镜自动平台控制系统,进行图像采集的定标设置; 
(4)图像的全景扫描与采集:选择图像起点——右上角,使该点置零,通过平台控制移动载物台,选择终点——左下角,然后显示平台框架,在平台框架内进行图像扫描、拼接,得到一个完整的图像; 
(5)图像参数设置:打开显微图像分析系统,在系统中打开拼接好的完整图像,进行图像的参数设置,调节色彩、对比度、饱和度,使得图像清晰,便于后期图像处理; 
(6)煤岩光片裂隙统计处理分析:首先进行图像色彩分割,选择灰度,进行二值化,将煤岩基质和裂隙分割开来;然后清除微粒,去除非裂隙部分——孔隙带来的误差,同时手动删除、或者连接裂隙以降低自动统计的误差;手动测量裂隙,图像上裂隙会显示出不同的颜色,便于区分;最后对光片上的所有裂隙进行汇总统计。 
7.根据权利要求3所述的自动定量分析煤岩光片裂隙的方法,所述煤岩光片放到载物台上时,煤岩光片的一边要与载物台的边缘平行。 
CN201410029309.3A 2014-01-22 2014-01-22 一种自动定量分析煤岩光片裂隙的系统和方法 Expired - Fee Related CN103698333B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410029309.3A CN103698333B (zh) 2014-01-22 2014-01-22 一种自动定量分析煤岩光片裂隙的系统和方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410029309.3A CN103698333B (zh) 2014-01-22 2014-01-22 一种自动定量分析煤岩光片裂隙的系统和方法

Publications (2)

Publication Number Publication Date
CN103698333A true CN103698333A (zh) 2014-04-02
CN103698333B CN103698333B (zh) 2015-11-18

Family

ID=50359936

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410029309.3A Expired - Fee Related CN103698333B (zh) 2014-01-22 2014-01-22 一种自动定量分析煤岩光片裂隙的系统和方法

Country Status (1)

Country Link
CN (1) CN103698333B (zh)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104198493A (zh) * 2014-09-18 2014-12-10 中国石油大学(北京) 岩石的各向异性测量系统及其使用方法
CN104815805A (zh) * 2015-03-19 2015-08-05 南昌大学 一种led芯片的自动筛选系统及筛选方法
WO2015176615A1 (zh) * 2014-05-22 2015-11-26 马立强 一种煤岩裂隙发育红外辐射监测试验方法
CN105629453A (zh) * 2016-03-10 2016-06-01 德诺杰亿(北京)生物科技有限公司 荧光显微镜
CN106204463A (zh) * 2015-05-05 2016-12-07 中国石油化工股份有限公司 一种岩心扫描图像的处理方法
CN108303515A (zh) * 2018-01-30 2018-07-20 河南理工大学 一种煤矿井下煤储层裂隙系统观测统计方法
CN117571733A (zh) * 2023-11-21 2024-02-20 核工业北京地质研究院 一种机械破岩损伤特征测试装置及其测试方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4959549A (en) * 1986-11-28 1990-09-25 Commonwealth Scientific And Industrial Research Organisation Determination of properties of coal
CN1635360A (zh) * 2003-12-26 2005-07-06 宝钢集团上海梅山有限公司 一种煤岩显微组分的测定系统及测定方法
CN103163134A (zh) * 2013-02-16 2013-06-19 肖建新 煤储层裂隙观测亚微观分析技术
CN103454150A (zh) * 2013-09-16 2013-12-18 安徽理工大学 一种全视角三维可调的多用途细观可视装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4959549A (en) * 1986-11-28 1990-09-25 Commonwealth Scientific And Industrial Research Organisation Determination of properties of coal
CN1635360A (zh) * 2003-12-26 2005-07-06 宝钢集团上海梅山有限公司 一种煤岩显微组分的测定系统及测定方法
CN103163134A (zh) * 2013-02-16 2013-06-19 肖建新 煤储层裂隙观测亚微观分析技术
CN103454150A (zh) * 2013-09-16 2013-12-18 安徽理工大学 一种全视角三维可调的多用途细观可视装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
姚艳斌等: "沁水盆地煤储层微裂隙发育的煤岩学控制机理", 《中国矿业大学学报》, vol. 39, no. 1, 31 January 2010 (2010-01-31), pages 6 - 13 *
张慧等: "煤中显微裂隙的成因类型及其研究意义", 《岩石矿物学杂质》, vol. 21, no. 3, 30 September 2002 (2002-09-30), pages 278 - 284 *
王越等: "型煤冷态强度影响因素的光学显微分析", 《煤炭转化》, vol. 36, no. 2, 30 April 2013 (2013-04-30), pages 79 - 82 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015176615A1 (zh) * 2014-05-22 2015-11-26 马立强 一种煤岩裂隙发育红外辐射监测试验方法
CN104198493A (zh) * 2014-09-18 2014-12-10 中国石油大学(北京) 岩石的各向异性测量系统及其使用方法
CN104198493B (zh) * 2014-09-18 2017-02-15 中国石油大学(北京) 岩石的各向异性测量系统及其使用方法
CN104815805A (zh) * 2015-03-19 2015-08-05 南昌大学 一种led芯片的自动筛选系统及筛选方法
CN106204463A (zh) * 2015-05-05 2016-12-07 中国石油化工股份有限公司 一种岩心扫描图像的处理方法
CN105629453A (zh) * 2016-03-10 2016-06-01 德诺杰亿(北京)生物科技有限公司 荧光显微镜
CN108303515A (zh) * 2018-01-30 2018-07-20 河南理工大学 一种煤矿井下煤储层裂隙系统观测统计方法
CN117571733A (zh) * 2023-11-21 2024-02-20 核工业北京地质研究院 一种机械破岩损伤特征测试装置及其测试方法

Also Published As

Publication number Publication date
CN103698333B (zh) 2015-11-18

Similar Documents

Publication Publication Date Title
CN103698333B (zh) 一种自动定量分析煤岩光片裂隙的系统和方法
CN110383038B (zh) 用于对空气样本进行自动分析的系统和方法
US11915407B2 (en) Automated system and method for clarity measurements and clarity grading
CN103207187B (zh) 外观检查装置及外观检查方法
CN106312062A (zh) 一种检验铺粉质量的方法及增材制造设备
JP2024069188A (ja) 空気サンプルの自動分析をするためのシステムおよび方法
CN103499303B (zh) 一种羊毛细度自动测量方法
CN107643295A (zh) 一种基于机器视觉的布匹缺陷在线检测的方法和系统
CN107004123A (zh) 迭代式缺陷滤除工艺
CN115908432B (zh) 一种物料输出品质检测系统及预测方法
Dan et al. Application of deep learning-based image recognition technology to asphalt–aggregate mixtures: Methodology
CN102854195B (zh) 彩色滤光片上缺陷坐标的检测方法
CN101556251A (zh) 一种基于数字信号处理器的ctp制版质量检测方法
CN103839271A (zh) 一种基于最佳匹配的图像纹理合成方法
CN104614383A (zh) 一种羽毛缺陷检测系统及方法
TW201512649A (zh) 偵測晶片影像瑕疵方法及其系統與電腦程式產品
US9972079B2 (en) Wafer appearance inspection apparatus
JP6049052B2 (ja) ウエハ外観検査装置及びウエハ外観検査装置における感度しきい値設定方法
CN108765391A (zh) 一种基于深度学习的平板玻璃异物影像分析方法
CN105928952B (zh) Aoi控制系统及其控制方法
CN105572133A (zh) 瑕疵检测方法及其装置
CN204389406U (zh) 蓝宝石长晶缺陷及表面缺陷光学检测系统
CN111986802A (zh) 一种肺腺癌病理分化等级的辅助判定系统及方法
Moaveni Advanced image analysis and techniques for degradation characterization of aggregates
CN104165696A (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: 20151118

CF01 Termination of patent right due to non-payment of annual fee