CN106501310B - 基于采空区空气参数实测的遗煤氧化升温模拟试验方法 - Google Patents

基于采空区空气参数实测的遗煤氧化升温模拟试验方法 Download PDF

Info

Publication number
CN106501310B
CN106501310B CN201610940766.7A CN201610940766A CN106501310B CN 106501310 B CN106501310 B CN 106501310B CN 201610940766 A CN201610940766 A CN 201610940766A CN 106501310 B CN106501310 B CN 106501310B
Authority
CN
China
Prior art keywords
goaf
coal
gas
air
air parameter
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
Application number
CN201610940766.7A
Other languages
English (en)
Other versions
CN106501310A (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 Mining and Technology Beijing CUMTB
Original Assignee
China University of Mining and Technology Beijing CUMTB
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 Mining and Technology Beijing CUMTB filed Critical China University of Mining and Technology Beijing CUMTB
Priority to CN201610940766.7A priority Critical patent/CN106501310B/zh
Publication of CN106501310A publication Critical patent/CN106501310A/zh
Application granted granted Critical
Publication of CN106501310B publication Critical patent/CN106501310B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N25/00Investigating or analyzing materials by the use of thermal means
    • G01N25/20Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity
    • G01N25/22Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity on combustion or catalytic oxidation, e.g. of components of gas mixtures

Abstract

本发明公开一种基于采空区空气参数实测的遗煤氧化升温模拟试验方法,属于煤矿采空区遗煤自燃防治研究领域,其首先通过束管监测方法测量得到采空区真实的空气参数,绘制采空区空气参数随时间变化的曲线,利用计算机编制伺服控制程序,然后进行遗煤氧化升温模拟试验,得到采空区遗煤在工作面推进过程中的温度和指标气体浓度变化曲线,进而可以分析采空区遗煤氧化升温特征及变化规律。本发明可以对采空区遗煤氧化自燃特性进行科学的观测和准确的分析,有利于对于煤矿采空区遗煤自燃灾害防治的深入研究。

Description

基于采空区空气参数实测的遗煤氧化升温模拟试验方法
技术领域
本发明属于煤矿采空区遗煤自燃防治研究领域,具体涉及一种基于采空区空气参数实测的遗煤氧化升温模拟试验方法。
背景技术
煤矿行业生产过程中,采空区遗煤自燃给煤层安全开采造成极大威胁。采空区遗煤氧化升温自燃,导致采空区火灾,高温和大量的有毒有害气体严重损害井下工作人员的生命健康,并损坏工作面设备,造成巨大的经济损失。研究采空区遗煤氧化升温特性,有利于指导采空区遗煤自燃的防治,促进工作面的安全回采。煤矿采空区在工作推采方向上,氧气、温度均有一定的变化规律,由此使得采空区没内遗煤氧化升温存在一定的分布规律,划分出采空区散热带、氧化带、窒息带,但从煤矿开采工程经验来看,但各矿区之间差别较大,即便同一矿区的不同采区,亦有各自的特点。
本研究领域内对煤体氧化升温特征的现有研究方法,大多是于实验室中,人为设定实验参数条件,进行煤的氧化升温特性分析。而对于某一具体开采条件下的工作面采空区遗煤的氧化升温过程,通过传统的束管监测,难以实现精确的分析,且采空区内环境复杂多变,测量误差给测定造成了较大干扰。
本发明针对这一现状,致力于发明一种基于采空区空气参数实测的遗煤氧化升温模拟试验方法,能够基于采空区空气条件,进行采空区煤样氧化升温模拟试验。
发明内容
一种基于采空区空气参数实测的遗煤氧化升温模拟试验方法,其特征在于:主要操作步骤如下:
步骤一、采空区空气参数采集,其包括以下分步骤:
a根据工作面长度,沿工作面布置若干等距节点I1,I2,…,In
b结合煤矿采空区束管监测技术,沿工作面布置一条束管,每个节点处对应埋设一路测线,所述测线包括气样采集管、气流传感器线路。在每个节点位置处,束管设置一个接口,接口通过快速接头与一保护筒相连,所述保护筒内包含气样采集端头、气流传感器,保护筒周身布设有多个气孔,便于气流流通和采集。束管在工作面回风巷道内的末端连接采气站,节点位置处的采空区空气通过气样采集管流至采气站,节点位置处的采空区漏风风流速度由气流传感器测定,并将电信号通过气流传感器线路传至采气站处的电子显示表;
c在工作面推进过程中,每日定时采集采空区空气参数数据,所述采空区空气参数包括采空区内空气的氧气浓度、空气湿度、漏风风流速度,其中,氧气浓度和空气湿度由井下工作人员每日定时从束管采气站中抽取气样,带至地面进行化验获得,漏风风流速度由工作人员通过采气站处的电子显示表读数获得;
步骤二、绘制采空区空气参数变化曲线
根据步骤一所得到的各节点处的空气参数数据和采集时间,绘制不同时刻采空区I1,I2,…,In节点位置处的采空区空气参数对应表,利用origin软件进行数据拟合,得到各节点位置处的空气参数随时间的变化曲线;
步骤三、编制采空区空气参数变化伺服程序
依据步骤二中得到的采空区空气参数随时间的变化曲线,利用计算机编程方法,编制模拟采空区空气参数变化的伺服控制程序,并导入计算机;
步骤四、进行煤样氧化升温试验,其包括以下分步骤:
a准备实验硬件,所述实验硬件为遗煤氧化升温模拟试验台,包括控制计算机、试验箱、加热器、氧气瓶、加湿机、鼓风机、数据监测系统、气相色谱分析仪;
所述计算机装载有伺服控制程序,可通过编制伺服控制程序,实现对试验系统各部分装置的伺服控制,由此实现对试验箱内氧气浓度、温度、湿度、气流速度的调节;
所述试验箱为一立方体封闭箱体,其外壳由绝热材料制作,顶部有圆形箱盖,可旋接于箱体实现密闭,试验箱右侧为进风口,左侧为出风口,进风口与鼓风机相连,出风口设置取气点,与气相色谱分析仪相连;
所述数据监测系统由多个传感器及线路组成,所述传感器包括氧气浓度传感器、湿度传感器、温度传感器;
所述气相色谱分析仪可对出风口处的气样进行收集和分析,实时监测和显示试验箱体内煤自燃指标气体浓度,所述煤自燃指标气体包括乙烯、一氧化碳、乙炔;
b将煤矿工作面现场采集的松散煤样置于试验箱中,启动计算机中的伺服控制程序,调整试验箱内空气参数至初始值,其中,试验箱体内的氧气由氧气瓶提供,湿度由加湿机通过连接在试验箱内的喷头喷射水雾进行调节,初始温度值设定为工作面温度,由计算机控制加热器进行调节;
c由计算机按照采空区内实际条件,连续模拟采空区空气参数变化;
d通过温度传感器实时监测煤样温度,生成数据变化曲线,通过I1,I2,…,In多个节点位置处对应的煤样氧化升温变化曲线;通过气相色谱分析仪对气样的分析,得到采空区遗煤氧化过程中的指标气体浓度变化曲线;
本发明所带来的有益成果在于:
本发明提供的一种基于采空区空气参数实测的遗煤氧化升温模拟试验方法,通过现场测量采空区真实的空气参数,利用计算机编制伺服控制程序,在实验室内实现自动化的升温氧化试验,操作方便快捷。通过氧化升温试验,进行连续的遗煤氧化升温模拟和观测,可深入分析遗煤的氧化升温特性和机理,以及各项遗煤自燃指标气体浓度变化特征,相比采空区的复杂环境,试验箱体内的环境相对封闭,环境干扰因素少,方便对采空区遗煤氧化自燃特性进行科学的观测和准确的分析,有利于对于煤矿采空区遗煤自燃灾害防治的深入研究。
附图说明
下面结合附图,对本发明作进一步清楚、完整的说明:
图1为基于采空区空气参数实测的遗煤氧化升温模拟试验方法流程图;
图2为束管接口处结构示意图;
图3为I1节点位置空气参数随时间的变化曲线;
图4为遗煤氧化升温模拟试验系统结构示意图;
图中:1、气样采集管,2、气流传感器线路,3、快速接头,4、保护筒,5、气样采集端头,6、气流传感器,7、气孔,8、计算机,9、试验箱,10、加热器,11、氧气瓶,12、加湿机,13、鼓风机,14、气相色谱分析仪,15、箱盖,16、进风口,17、出风口,18、氧气浓度传感器,19、湿度传感器,20、温度传感器。
具体实施方式
本发明提供了一种基于采空区空气参数实测的遗煤氧化升温模拟试验方法,下面通过附图和实施例,以某矿S1煤层3302工作面开采为工程背景,实际应用本技术。对本发明的技术方案做进一步的详细描述。
一种基于采空区空气参数实测的遗煤氧化升温模拟试验方法,其特征在于:主要操作步骤如下:
步骤一、采空区空气参数采集,其包括以下分步骤:
a 3302工作面长160m,结合工作面开采和设备条件,将工作面以20m为间距,划分得到9个节点I1,I2,…,I9
b沿3302工作面布置一条束管,每个节点处对应埋设一路测线,所述测线包括气样采集管1、气流传感器线路2。在每个节点位置处,束管设置一个接口,接口通过快速接头3与一保护筒4相连,所述保护筒4内包含气样采集端头5、气流传感器6,保护筒4周身布设有多个气孔7,便于气流流通和采集。束管在工作面回风巷道内的末端连接采气站,节点位置处的采空区空气通过气样采集管1流至采气站,节点位置处的采空区漏风风流速度由气流传感器6测定,并将电信号通过气流传感器线路2传至采气站处的电子显示表;
c在工作面推进过程中,每日定时采集采空区空气参数数据,所述采空区空气参数包括采空区内空气的氧气浓度、空气湿度、漏风速度,其中,氧气浓度和空气湿度由井下工作人员每日定时从束管采气站中抽取气样,带至地面进行化验获得,漏风速度由工作人员通过采气站处的电子显示表读数获得;
步骤二、绘制采空区空气参数变化曲线
根据步骤一所得到的各节点处的空气参数数据和采集时间,绘制不同时刻采空区I1,I2,…,I9位置处的采空区空气参数对应表,利用origin软件进行数据拟合,得到各位置处的空气参数随时间的变化曲线,以I1节点为例进行说明,I1节点位置处的空气参数变化曲线如图3所示;
步骤三、编制采空区空气参数变化伺服程序
依据步骤二中的采空区空气参数随时间的变化曲线,利用计算机C语言编程方法,编制模拟采空区空气参数变化的伺服控制程序,并导入计算机8内;
步骤四、进行煤样氧化升温试验,其包括以下分步骤:
a准备实验硬件,所述实验硬件为遗煤氧化升温模拟试验台,包括控制计算机8、试验箱9、加热器10、氧气瓶11、加湿机12、鼓风机13、数据监测系统、气相色谱分析仪14;
所述计算机8装载有伺服控制程序,可通过编辑伺服控制程序,实现对试验系统各部分装置的伺服控制,由此实现对试验箱9内氧气浓度、温度、湿度、气流速度的调节;
所述试验箱9为一立方体封闭箱体,其外壳由绝热材料制作,顶部有圆形箱盖15,可旋接箱体实现密闭,试验箱9右侧为进风口16,左侧为出风口17,进风口16与鼓风机13相连,出风口17设置取气点与气相色谱分析仪14相连;
所述数据监测系统由多个传感器及线路组成,所述传感器包括氧气浓度传感器18、湿度传感器19、温度传感器20;
所述气相色谱分析仪14可对出风口17处的气样进行收集和分析,实时监测和显示试验箱9内煤自燃指标气体浓度,所述煤自燃指标气体包括乙烯、一氧化碳、乙炔;
b将煤矿工作面现场采集的松散煤样置于试验箱9中,启动计算机中的伺服控制程序,调整试验箱9内空气参数至初始值,其中,试验箱9内的氧气由氧气瓶11提供,湿度由加湿机12通过连接在试验箱内的喷头喷射水雾进行调节,初始温度设定为工作面温度19℃,由计算机8控制加热器10进行设置;
c由计算机按照采空区内实际条件,连续模拟采空区空气参数变化;
d通过温度传感器实时监测煤样温度,生成数据变化曲线,通过I1,I2,…,I9多个节点位置处对应的煤样氧化升温变化曲线;通过气相色谱分析仪14对气样的分析,得到上述各节点位置处遗煤氧化过程中的指标气体浓度变化曲线。

Claims (1)

1.一种基于采空区空气参数实测的遗煤氧化升温模拟试验方法,其特征在于,主要操作步骤如下:
步骤一、采空区空气参数采集,其包括以下分步骤:
a根据工作面长度,沿工作面布置若干等距节点I1,I2,…,In
b结合煤矿采空区束管监测技术,沿工作面布置一条束管,每个节点处对应埋设一路测线,所述测线包括气样采集管、气流传感器线路;在每个节点位置处,束管设置一个接口,接口通过快速接头与一保护筒相连,所述保护筒内包含气样采集端头、气流传感器,节点位置处的采空区漏风风流速度由气流传感器测定,并将电信号通过气流传感器线路传至采气站处的电子显示表;
c在工作面推进过程中,每日定时采集采空区空气参数数据,所述采空区空气参数包括采空区内空气的氧气浓度、空气湿度、漏风风流速度,其中,氧气浓度和空气湿度由井下工作人员每日定时从束管采气站中抽取气样,带至地面进行化验获得,漏风风流速度由工作人员通过采气站处的电子显示表读数获得;
步骤二、绘制采空区空气参数变化曲线
根据步骤一所得到的各节点处的空气参数数据和采集时间,绘制不同时刻采空区I1,I2,…,In节点位置处的采空区空气参数对应表,利用origin软件进行数据拟合,得到各节点位置处的空气参数随时间的变化曲线;
步骤三、编制采空区空气参数变化伺服程序
依据步骤二中的采空区空气参数随时间的变化曲线,利用计算机编程方法,编制模拟采空区空气参数变化的伺服控制程序,并导入计算机;
步骤四、进行煤样氧化升温试验,其包括以下分步骤:
a准备实验硬件,所述实验硬件为遗煤氧化升温模拟试验台,包括控制计算机、试验箱、加热器、氧气瓶、加湿机、鼓风机、数据监测系统、气相色谱分析仪;
所述计算机装载有伺服控制程序,可通过编制伺服控制程序,实现对试验系统各部分装置的伺服控制,由此实现对试验箱内氧气浓度、温度、湿度、气流速度的调节;
所述试验箱为一立方体封闭箱体,其外壳由绝热材料制作,顶部有圆形箱盖,可旋接箱体实现密闭,试验箱右侧为进风口,左侧为出风口,进风口与鼓风机相连,出风口设置取气点,与气相色谱分析仪相连;
所述数据监测系统由多个传感器及线路组成,所述传感器包括氧气浓度传感器、湿度传感器、温度传感器;
所述气相色谱分析仪可对出风口处的气样进行收集和分析,实时监测和显示试验箱体内煤自燃指标气体浓度,所述煤自燃指标气体包括乙烯、一氧化碳、乙炔;
b将煤矿工作面现场采集的松散煤样置于试验箱中,启动计算机中的伺服控制程序,调整试验箱内空气参数至初始值,其中,试验箱体内的氧气由氧气瓶提供,湿度由加湿机通过连接在试验箱内的喷头喷射水雾进行调节,初始温度值设定为工作面温度,由计算机控制加热器进行调节;
c由计算机按照采空区内实际条件,连续模拟采空区空气参数变化;
d通过温度传感器实时监测煤样温度,生成数据变化曲线,通过I1,I2,…,In多个节点位置处对应的煤样氧化升温变化曲线;通过气相色谱分析仪对气样的分析,得到采空区遗煤氧化过程中的指标气体浓度变化曲线。
CN201610940766.7A 2016-10-25 2016-10-25 基于采空区空气参数实测的遗煤氧化升温模拟试验方法 Expired - Fee Related CN106501310B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610940766.7A CN106501310B (zh) 2016-10-25 2016-10-25 基于采空区空气参数实测的遗煤氧化升温模拟试验方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610940766.7A CN106501310B (zh) 2016-10-25 2016-10-25 基于采空区空气参数实测的遗煤氧化升温模拟试验方法

Publications (2)

Publication Number Publication Date
CN106501310A CN106501310A (zh) 2017-03-15
CN106501310B true CN106501310B (zh) 2019-02-12

Family

ID=58321725

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610940766.7A Expired - Fee Related CN106501310B (zh) 2016-10-25 2016-10-25 基于采空区空气参数实测的遗煤氧化升温模拟试验方法

Country Status (1)

Country Link
CN (1) CN106501310B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110687149A (zh) * 2019-09-02 2020-01-14 中国地质大学(北京) 一种地下煤田火区典型产物动态析出释放模拟实验系统
CN111504368B (zh) * 2020-04-15 2022-03-01 河南理工大学 一种高应力条件下采空区煤自燃三带划分及实验方法
CN112083142A (zh) * 2020-09-21 2020-12-15 辽宁工程技术大学 一种模拟采空区火灾指标气体产生的实验装置及方法
CN112798707B (zh) * 2020-12-29 2023-01-13 国家能源集团宁夏煤业有限责任公司 一种矿用监测设备及采空区气体监测方法
CN113295769B (zh) * 2021-05-20 2023-05-30 山东科技大学 煤自燃高温异常区域声电信号传输特性测试的模拟装置及方法
CN114166969A (zh) * 2021-11-30 2022-03-11 西安科技大学 模拟微生物代谢煤中可燃性气体的反应装置及其试验方法
CN114496104B (zh) * 2022-04-02 2022-06-28 中国矿业大学(北京) 一种采煤条件下采空区遗煤自燃程度评估方法及系统

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU516951A1 (ru) * 1974-07-25 1976-06-05 Всесоюзный Научно-Исследовательский И Проектно-Конструкторский Угольный Институт Устройство дл определени температуры воспламенени твердых углеродистых материалов
JPS58140633A (ja) * 1981-07-24 1983-08-20 Central Res Inst Of Electric Power Ind 貯炭パイルの自然発火監視方法
CN1595129A (zh) * 2004-07-16 2005-03-16 中国矿业大学 一种模拟煤炭自燃过程的绝热测试方法及其装置
CN201589761U (zh) * 2009-12-23 2010-09-22 煤炭科学研究总院沈阳研究院 一种基于程序控温的煤自燃模拟装置
CN103399035A (zh) * 2013-07-30 2013-11-20 辽宁工程技术大学 多功能煤样自燃氧化与升温实验装置
CN204374123U (zh) * 2014-12-11 2015-06-03 河南工程学院 一种煤自燃特征试验系统
CN104943995A (zh) * 2015-06-16 2015-09-30 西安科技大学 自加热煤样罐及用该罐进行煤自燃特性测试加热煤的方法
CN105823796A (zh) * 2016-03-11 2016-08-03 华电电力科学研究院 一种煤场储煤自燃倾向性的测定方法
CN105842281A (zh) * 2016-04-28 2016-08-10 煤炭科学技术研究院有限公司 一种煤矿采空区温度场分布规律相似模拟试验装置

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU516951A1 (ru) * 1974-07-25 1976-06-05 Всесоюзный Научно-Исследовательский И Проектно-Конструкторский Угольный Институт Устройство дл определени температуры воспламенени твердых углеродистых материалов
JPS58140633A (ja) * 1981-07-24 1983-08-20 Central Res Inst Of Electric Power Ind 貯炭パイルの自然発火監視方法
CN1595129A (zh) * 2004-07-16 2005-03-16 中国矿业大学 一种模拟煤炭自燃过程的绝热测试方法及其装置
CN201589761U (zh) * 2009-12-23 2010-09-22 煤炭科学研究总院沈阳研究院 一种基于程序控温的煤自燃模拟装置
CN103399035A (zh) * 2013-07-30 2013-11-20 辽宁工程技术大学 多功能煤样自燃氧化与升温实验装置
CN204374123U (zh) * 2014-12-11 2015-06-03 河南工程学院 一种煤自燃特征试验系统
CN104943995A (zh) * 2015-06-16 2015-09-30 西安科技大学 自加热煤样罐及用该罐进行煤自燃特性测试加热煤的方法
CN105823796A (zh) * 2016-03-11 2016-08-03 华电电力科学研究院 一种煤场储煤自燃倾向性的测定方法
CN105842281A (zh) * 2016-04-28 2016-08-10 煤炭科学技术研究院有限公司 一种煤矿采空区温度场分布规律相似模拟试验装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《煤自燃全过程测试和指标气体的研究与应用》;许延辉.;《中国优秀博硕士学位论文全文数据库 (硕士) 工程科技Ⅰ辑》;20060715;正文第5-46页 *

Also Published As

Publication number Publication date
CN106501310A (zh) 2017-03-15

Similar Documents

Publication Publication Date Title
CN106501310B (zh) 基于采空区空气参数实测的遗煤氧化升温模拟试验方法
CN103114827B (zh) 多场耦合煤层气抽采模拟试验方法
CN205844228U (zh) 一种高温矿井巷道通风降温的实验装置
CN102230375A (zh) 煤层瓦斯参数实时监测装置及方法
CN110925008B (zh) 一种煤矿局部通风设备智能调节测试实验平台
CN110905580B (zh) 一种综掘工作面通风优化方法
CN105954459B (zh) 一种检测煤矿采空区瓦斯抽采对煤自燃影响的装置
CN104990687A (zh) 研究综掘面巷道中粉尘运移规律的可视化试验装置及方法
CN102003199A (zh) 煤矿井下热湿环境模拟控制实验系统
CN111650242A (zh) 检测煤自燃产热产气特性对瓦斯浓度场影响的装置及方法
CN207992175U (zh) 二氧化碳地质封存泄露人工控制模拟装置
CN105510205A (zh) 一种氡在煤岩体多孔介质中迁移参数的采集装置
CN103940723B (zh) 实验室快速测定煤矿井下工作面前方煤体渗透率的方法
CN111141866B (zh) 模拟煤层自燃过程及高温异常区域探测的试验装置及方法
CN111120003A (zh) 一种模拟瓦斯与煤自燃共生灾害的试验装置及工作方法
CN104749347A (zh) 研究地铁区间隧道土壤温湿度场演化规律的实验装置
CN201924953U (zh) 煤矿井下热湿环境模拟控制实验系统
CN204346498U (zh) 一种矿井通风机风速、全压监测监控装置
CN203745454U (zh) 煤层瓦斯含量自动测试装置
CN210948653U (zh) 一种快速评价页岩气井产能的装置
CN206111158U (zh) 稠油油藏热采线性物理模拟测控装置
CN111720101B (zh) 页岩气井压裂返排的模拟系统以及模拟方法
CN108825306A (zh) 一种高地温隧道的温度监测装置及监测方法
CN206557188U (zh) 一种矿井巷道烟气流动模拟测试系统
CN110685678A (zh) 一种快速评价页岩气井产能的方法及装置

Legal Events

Date Code Title Description
C06 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: 20190212

Termination date: 20191025