CN104739389B - 曲面应变片、脉搏信号提取装置及方法和脉象诊疗*** - Google Patents

曲面应变片、脉搏信号提取装置及方法和脉象诊疗*** Download PDF

Info

Publication number
CN104739389B
CN104739389B CN201510177740.7A CN201510177740A CN104739389B CN 104739389 B CN104739389 B CN 104739389B CN 201510177740 A CN201510177740 A CN 201510177740A CN 104739389 B CN104739389 B CN 104739389B
Authority
CN
China
Prior art keywords
curved surface
pulse
fiber grating
surface strain
pulse signal
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
CN201510177740.7A
Other languages
English (en)
Other versions
CN104739389A (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.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN201510177740.7A priority Critical patent/CN104739389B/zh
Publication of CN104739389A publication Critical patent/CN104739389A/zh
Application granted granted Critical
Publication of CN104739389B publication Critical patent/CN104739389B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0002—Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network
    • A61B5/0015—Remote monitoring of patients using telemetry, e.g. transmission of vital signals via a communication network characterised by features of the telemetry system
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/021—Measuring pressure in heart or blood vessels
    • A61B5/02108—Measuring pressure in heart or blood vessels from analysis of pulse wave characteristics
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/7253—Details of waveform analysis characterised by using transforms
    • A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235—Details of waveform analysis
    • A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • Public Health (AREA)
  • Animal Behavior & Ethology (AREA)
  • Surgery (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Artificial Intelligence (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Signal Processing (AREA)
  • Psychiatry (AREA)
  • Mathematical Physics (AREA)
  • Cardiology (AREA)
  • Evolutionary Computation (AREA)
  • Fuzzy Systems (AREA)
  • Vascular Medicine (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)

Abstract

本发明属于光学生物脉搏信息监测及远程医疗领域，具体涉及光纤光栅曲面应变片、脉搏信号提取装置及方法和脉象诊疗***。用于脉搏信号提取的光纤光栅曲面应变片，包括光纤光栅和曲面应变片；本发明还提供了一种脉搏信号提取装置，由腕带、光纤探头组成，所述光纤探头包括上述的光纤光栅曲面应变片、外壳和光纤跳线；所述光纤光栅曲面应变片置于外壳内，所述光纤跳线连接光纤光栅曲面应变片的光纤光栅，并延伸至出所述外壳；所述光纤探头埋入所述腕带内。本发明利用光学检测灵敏度高、噪声低的特点，基于光纤微震动传感器精确提取人体脉搏信息，充分减少人为因素影响，并利用光纤传感可通过远程传输采集数据的特点实现“悬丝诊脉”。

Description

曲面应变片、脉搏信号提取装置及方法和脉象诊疗***
技术领域
本发明属于光学生物脉搏信息监测及远程医疗领域，具体涉及光纤光栅曲面应变片、脉搏信号提取装置及方法和脉象诊疗***。
背景技术
中医诊脉是一门古老的技术，采用仪器辅助诊脉是一项技术难题，这是由于人体脉象的复杂性、易受干扰性及个体差异性等特点，即使对于经验丰富的老中医而言，脉象仍是只可意会难以言传，一般检测手段很难采集到丰富而详尽的脉象信息，并且非常容易受到人体其它行为的干扰，因此无法给出客观全面的诊脉结果。
针对中医脉象信息进行提取，现有技术主要是采用配备有单探头换能器的脉象仪。人们利用这类仪器可以初步识别浮、沉、虚、实、滑、涩等实际中脉象，但是无法精确测量出多种脉象混叠的数据，因此只能针对某些病患提出试探性的病症诊断。这样的医疗仪器是无法适应中西医结合需求的。光纤传感器是目前能够采用的精度最高的传感器。无论是压电晶体式、光电探测式、压敏电阻式还是超声波式，都无法达到光纤传感器的检测精度。针对漫射光层析成像***的复杂算法和蒙特卡洛模拟的不确定性，光纤传感器的可行性更高，可以测量出的脉象数据预计也能远远超过现有水平。***的体积和重量决定应用场景。目前国内外在相关领域，会使用心电图、脑电图等方式进行脉搏信号的检测、体积较大、价格昂贵，同时在人体上夹持的部位过多。例如心脏手术，如果需要监测脉搏信号，采用心电图显然是不现实的。
在光纤传感器方面，目前医学上应用的主要是传光型光纤传感器。但在目前脉诊客观化的应用上，光纤主要用来传递信号而非感应获取信号，现有的光纤传感装置实质上仍为光电式装置。目前也有提出采用双光纤光栅来消除环境干扰测量脉搏信息的方法，但只是停留在信号的初步提取没有进一步分析及数据处理，停留在能测出脉搏波形的可行性上。该方法是在一根光纤上刻了两个光栅，一个光栅用来测量脉搏信号，另一个用来测量环境干扰。
发明内容
本发明针对的技术问题，为了提高脉搏信号的测量精度以及使传感器小型化，我们选择了光纤光栅。针对现有的光纤光栅传感器没有测量脉搏这种低频、微弱的信号，提出了用光纤光栅曲面应变片来提高光纤光栅传感器与脉搏位置的贴合度与灵敏度，同时一定程度上减弱了干扰。通过光纤光栅传感器测量的精确信号，使得相近脉象的区分更加明显。所采用的光纤为普通单模光纤，与现有光缆网络组网方便，提高了远程医疗的可靠性。实现了脉诊的客观化、精确化、便捷化、可靠化。除了确定脉象外，为方便医者的诊断，同时提供体温、心率的参数，为临床医疗以及远程医疗提供可靠依据。
具体技术方案如下：
一种用于脉搏信号提取的光纤光栅曲面应变片，包括光纤光栅和曲面应变片；所述曲面应变片由正面应变片和背面应变片两片对称的应变片粘合而成，正面应变片镂空成带有双十字型槽，背面应变片镂空成一字型槽，背面应变片镂空位置与正面应变片镂空部分对应设置；每片应变片均由两端为矩形形状面板和中部结构为双T形状面板的一体化结构组成，且沿纵向轴对称，两端沿纵向轴对称向中间弯曲；所述正面应变片纵向轴线上的槽用于放置光纤光栅；光纤光栅整体处于曲面应变片内；所述光纤光栅上的光栅部分置于槽内，剩余部分固定于曲面应变片上；所述双T形状面板上的沿纵向轴的镂空槽是相通的。
进一步地，所述光纤光栅剩余部分的固定方式为无影胶封装的方式或焊料点焊的方式。
本发明还提供了一种脉搏信号提取装置，由腕带、光纤探头组成，所述光纤探头包括上述的光纤光栅曲面应变片、外壳和光纤跳线；所述光纤光栅曲面应变片置于外壳内，所述光纤跳线连接光纤光栅曲面应变片的光纤光栅，并延伸至出所述外壳；所述光纤探头埋入所述腕带内。
进一步地，所述外壳和腕带采用硅胶或橡胶材料制成。
本发明的一种脉象诊疗***，包括家庭终端，医疗终端和传输光缆；所述家庭终端包括计算机和所述的脉搏信号提取装置；所述医疗终端包括计算机、光纤光栅解调***、数据处理分析***和脉象数据库，家庭终端采集信号，传输光缆传递信号至医疗终端进行解调、处理和分析后，输出脉象信息。
本发明提供的一种脉搏波信号处理方法，包括如下步骤：
(1)用上述的脉搏信号提取装置测量脉搏波信息；
(2)对脉搏波信息的反射中心波长变化做谱图分析，获得特征参量，建
立脉象的概率神经网络(PNN)模型；
(3)将脉象的概率神经网络(PNN)模型与数据库中的脉象数据比对，得到脉象结果。
采用本发明获得的技术效果：本发明利用光学检测灵敏度高、噪声低的特点，基于光纤微震动传感器精确提取人体脉搏信息，充分减少人为因素影响，并利用光纤传感可通过远程传输采集数据的特点实现“悬丝诊脉”。在此基础上通过大量数据采集，对中医脉象建立数据模型，建立起准确而完善的光纤脉象分析诊疗***。
说明书附图
图1是光纤光栅曲面应变片结构正面示意图；
图2是光纤光栅曲面应变片结构背面示意图；
图3是光纤光栅曲面应变片结构侧面示意图；
图4是脉搏信号提取装置结构示意图；
图5是光纤光栅曲面应变片结构参数设置示意图；
图6是本发明中的脉象诊疗***框图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。
如图1、图2、图3所示，一种光纤光栅曲面应变片1，包括光纤光栅和曲面应变片；所述曲面应变片由正面应变片和背面应变片两片对称的应变片粘合而成，正面应变片镂空成带有双十字型槽4，背面应变片镂空成一字型槽5，背面应变片镂空位置与正面应变片镂空部分对应设置；每片应变片均由两端为矩形形状面板2和中部结构为双T形状面板3的一体化结构组成，且沿纵向轴对称，两端沿纵向轴对称向中间弯曲；所述正面应变片纵向轴线上的槽用于放置光纤光栅；光纤光栅整体处于曲面应变片内；所述光纤光栅上的光栅部分置于槽内，剩余部分固定于曲面应变片上；所述双T形状面板上的沿纵向轴的镂空槽是相通的。双T形状面板中间部分部跨槽连接厚度为槽厚度的三分之一，起固定光栅作用。
如图4所示，一种脉搏信号提取装置，由腕带6、光纤探头7组成，所述光纤探头7包括光纤光栅曲面应变片1、外壳和光纤跳线；所述光纤光栅曲面应变片置于外壳内，所述光纤跳线连接光纤光栅曲面应变片的光纤光栅，并延伸至出所述外壳；所述光纤探头埋入所述腕带内。
光纤光栅曲面应变片弯曲一定程度圆弧，本实施例弯曲半径为为11.1cm的圆弧。可以调整为如图5中所示的a、b、c、d参数(本实施例的参数值是a＝7mm，b＝3mm，c＝1.5mm，d＝8mm)，提高其灵敏度。光纤光栅曲面应变片制作材料选择是金属或硅胶材料。
光纤探头中运用的光纤光栅曲面应变片，主要作用是：
脉搏波属于低频，微弱信号，裸露的光纤光栅在测量脉搏时，灵敏度不能达到要求，故采用曲面应变片进行增敏处理。在一定程度上增加了信噪比，减小了温度，肌肉抖动的误差。将测得的信号进行脉搏震动时间分布函数的转化，以此实时的监测被检测活体脉搏的震动周期，振幅，并通过其中的拓扑变化关系，为后期的脉搏——中医信号分析做了硬件上的支持。值得注意的是，整个过程中需要比较精准合适的信号放大处理，以便从复杂的背景中提取出有用的脉搏信息。
如图4所示，所述外壳和腕带采用硅胶材料制成。经过腕带包裹的光纤探头可以受到很好的保护，并可以固定在被检测者的手腕上，并且通过弹性放缩可以使光纤探头很好的贴近脉搏部位，利于提取脉搏信号。
如图6所示，本发明还提供了一种脉象诊疗***，包括家庭终端8，传输光缆9和医疗终端10；所述家庭终端包括所述的脉搏信号提取装置；所述医疗终端包括计算机、光纤光栅解调***、数据处理与分析***和脉象数据库；家庭终端采集到信号，传输光缆传递信号至医疗终端进行解调、处理和分析后，输出脉象信息。
光纤光栅解调***：对光纤光栅中心反射的微小偏移进行精确测量解调。
数据处理与分析***：接收到的信号通过数据的初步处理，形成信号波形，分析产生脉搏信号。
本发明提供的一种脉搏波信号处理方法，包括如下步骤：
(1)用本发明中的脉搏信号提取装置测量脉搏波信息；
(2)对脉搏波信息的反射中心波长变化做谱图分析，获得特征参量，建立脉象的概率神经网络(PNN)模型；
(3)将脉象的概率神经网络(PNN)模型与数据库中的脉象数据比对，得到脉象结果。
通过建立脉象的概率神经网络(PNN)模型最终确定并输出脉象。概率神经网络(PNN)模型的识别分类方法是概率统计学中被广泛接受的一种决策方法。该信号分析识别方法，采用三种谱图分析，分别是功率谱、倒谱、传递函数，从中确定八个特征参量。这八个特征参量是由测量光波中心波长变化通过谱图分析确定。
特征参量如下：
1.fi：经过傅里叶变换后的基频等信息。它反映了心率。
2.N：功率谱谐波个数。反映脉搏节律。
3.SER：谱能比。反映脉搏的能量随频率的分布。
4.Cx0：倒谱零分量。反映脉搏的强度大小。
5.M：传递函数的共振峰个数。反映脉搏***的共振特征。
6.△F：共振峰平均间距。
7.R.C.：倒谐波的特征。
8.Cx1/Cx0:倒谱第一谐波幅值与倒谱零分量之比。反映脉象的流利程度。
据此建立脉象识别的概率神经网络(PNN)模型，将上述8个特征参量依次简记为xi(i＝1～8)，将处理后的数据与数据库中的脉象数据比对，最终对脉象进行确定。
以上仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理前提下的若干改进和润饰，应视为本发明的保护范围。

Claims (6)

1.一种用于脉搏信号提取的光纤光栅曲面应变片，其特征在于：包括光纤光栅和曲面应变片；所述曲面应变片由正面应变片和背面应变片两片对称的应变片粘合而成，正面应变片镂空成带有双十字型槽，背面应变片镂空成一字型槽，背面应变片镂空位置与正面应变片镂空部分对应设置；每片应变片均由两端为矩形形状面板和中部结构为双T形状面板的一体化结构组成，且沿纵向轴对称，两端沿纵向轴对称向中间弯曲；所述正面应变片纵向轴线上的槽用于放置光纤光栅；光纤光栅整体处于曲面应变片内；所述光纤光栅上的光栅部分置于槽内，剩余部分固定于曲面应变片上；所述双T形状面板上的沿纵向轴的镂空槽是相通的。
2.如权利要求1所述的用于脉搏信号提取的光纤光栅曲面应变片，其特征在于，所述光纤光栅剩余部分的固定方式为无影胶封装的方式或焊料点焊的方式。
3.一种脉搏信号提取装置，其特征在于：由腕带、光纤探头组成，所述光纤探头包括外壳、光纤跳线和如权利要求1所述的光纤光栅曲面应变片；所述光纤光栅曲面应变片置于外壳内，所述光纤跳线连接光纤光栅曲面应变片的光纤光栅，并延伸至出所述外壳；所述光纤探头埋入所述腕带内。
4.如权利要求3所述的一种脉搏信号提取装置，所述外壳和腕带采用硅胶或橡胶材料制成。
5.一种脉象诊疗***，其特征在于包括家庭终端，医疗终端和传输光缆；所述家庭终端包括计算机和如权利要求3至4任一项所述的脉搏信号提取装置；所述医疗终端包括计算机、光纤光栅解调***、数据处理分析***和脉象数据库，家庭终端采集信号，传输光缆传递信号至医疗终端进行解调、处理和分析后，输出脉象信息。
6.一种脉搏波信号处理方法，其特征在于，包括如下步骤：
(1)用权利要求3所述的脉搏信号提取装置测量脉搏波信息；
(2)对脉搏波信息的反射中心波长变化做谱图分析，获得特征参量，建立脉象的概率神经网络模型；
(3)将脉象的概率神经网络模型与数据库中的脉象数据比对，得到脉象结果。
CN201510177740.7A 2015-04-15 2015-04-15 曲面应变片、脉搏信号提取装置及方法和脉象诊疗*** Expired - Fee Related CN104739389B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510177740.7A CN104739389B (zh) 2015-04-15 2015-04-15 曲面应变片、脉搏信号提取装置及方法和脉象诊疗***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510177740.7A CN104739389B (zh) 2015-04-15 2015-04-15 曲面应变片、脉搏信号提取装置及方法和脉象诊疗***

Publications (2)

Publication Number Publication Date
CN104739389A CN104739389A (zh) 2015-07-01
CN104739389B true CN104739389B (zh) 2017-04-26

Family

ID=53580056

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510177740.7A Expired - Fee Related CN104739389B (zh) 2015-04-15 2015-04-15 曲面应变片、脉搏信号提取装置及方法和脉象诊疗***

Country Status (1)

Country Link
CN (1) CN104739389B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106595510A (zh) * 2017-01-06 2017-04-26 成都翰光科技有限公司 一种应变片及光纤光栅应变传感器
CN108670209A (zh) * 2018-03-29 2018-10-19 中国科学院微电子研究所 一种自动识别中医脉象的方法和***
CN108922610A (zh) * 2018-08-08 2018-11-30 湖南川马医疗科技有限公司 一种基于大数据的智能定位脉诊分析仪及分析方法
CN110575148A (zh) * 2019-09-29 2019-12-17 北京一药良心信息科技有限公司 一种中医远程诊脉方法、装置和服务器
CN110742618A (zh) * 2019-10-29 2020-02-04 南通大学 一种吞咽智能检测***
CN114983363B (zh) * 2022-05-07 2024-06-28 华中科技大学 一种可穿戴锥形光纤血压监测装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102772200A (zh) * 2012-08-06 2012-11-14 山东省科学院激光研究所 脉诊信号传感器探头及光纤脉诊信息采集***
CN104367324A (zh) * 2013-08-16 2015-02-25 上海微创电生理医疗科技有限公司 压力传感装置及其制造方法、模具以及医疗导管

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8182433B2 (en) * 2005-03-04 2012-05-22 Endosense Sa Medical apparatus system having optical fiber load sensing capability
KR100908997B1 (ko) * 2007-09-12 2009-07-22 한국 한의학 연구원 광섬유격자 센서를 이용한 맥진 시스템

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102772200A (zh) * 2012-08-06 2012-11-14 山东省科学院激光研究所 脉诊信号传感器探头及光纤脉诊信息采集***
CN104367324A (zh) * 2013-08-16 2015-02-25 上海微创电生理医疗科技有限公司 压力传感装置及其制造方法、模具以及医疗导管

Also Published As

Publication number Publication date
CN104739389A (zh) 2015-07-01

Similar Documents

Publication Publication Date Title
CN104739389B (zh) 曲面应变片、脉搏信号提取装置及方法和脉象诊疗***
CN101330868B (zh) 用于监测用户心率和/或心率变化的监测设备、包括这种监测设备的手表
CN102772200B (zh) 脉诊信号传感器探头及光纤脉诊信息采集***
CN107205671A (zh) 至少部分基于脉搏波形的自动诊断
CN203953636U (zh) 一种血压检测装置、智能腕带、智能手表和通信***
CN104027108A (zh) 一种新型光学心电、脉搏综合检测装置
CN100518631C (zh) 三维动态脉象检测***
WO2006050385A2 (en) Cardiac motion characterization by strain measurement
CN103648373A (zh) 个人健康数据采集
CN105380618A (zh) 生物体信息解析***及生物体信息处理***
KR20140060737A (ko) 생체 신호 전송기, 생체 신호 수신기, 및 생체 신호 송수신 방법
CN106264504A (zh) 基于手指动脉波的无创血压检测***及方法
CN103932684B (zh) 一种新型光学脉诊仪***
CN106955090A (zh) 诊脉装置、气体检测装置、模拟脉搏装置和中医诊疗***
US10912482B2 (en) Method for determining at least one type and/or condition of cells and system
CN106361299B (zh) 一种三重薄膜心血管检测传感器
CN102355879A (zh) 用于向用户指示刺激信号的***、方法和计算机程序产品
CN217932165U (zh) 一种柔性微纳光纤耦合器及微应变传感应用***
DE102008042115A1 (de) Messvorrichtung sowie Messverfahren
CN110507295B (zh) 光纤感应组件及生命体征监测装置
CN104523260A (zh) 一种微循环血流速度测量装置及方法
CN207654141U (zh) 诊脉装置、气体检测装置、模拟脉搏装置和中医诊疗***
CN108158563A (zh) 一种基于低频可闻声波的中医脉诊仪
CN102711624B (zh) 使用身体上或身体中通信技术来获得关于用户的身体成分的信息的***和方法
CN112120679A (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
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: 20170426

Termination date: 20180415