CN103267543B - 显微粒子图像微流量测量仪及方法 - Google Patents
显微粒子图像微流量测量仪及方法 Download PDFInfo
- Publication number
- CN103267543B CN103267543B CN201310134340.9A CN201310134340A CN103267543B CN 103267543 B CN103267543 B CN 103267543B CN 201310134340 A CN201310134340 A CN 201310134340A CN 103267543 B CN103267543 B CN 103267543B
- Authority
- CN
- China
- Prior art keywords
- flow
- optical glass
- velocity
- dimensional
- glass rotating
- 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
- 238000005259 measurement Methods 0.000 title claims abstract description 36
- 239000002245 particle Substances 0.000 title claims abstract description 29
- 238000000034 method Methods 0.000 title claims abstract description 25
- 239000005304 optical glass Substances 0.000 claims abstract description 37
- 239000012530 fluid Substances 0.000 claims abstract description 17
- 238000010586 diagram Methods 0.000 claims abstract description 9
- 238000002347 injection Methods 0.000 claims abstract description 8
- 239000007924 injection Substances 0.000 claims abstract description 8
- 238000010276 construction Methods 0.000 claims abstract description 7
- 238000011084 recovery Methods 0.000 claims abstract description 5
- 238000012800 visualization Methods 0.000 claims abstract description 5
- 238000004458 analytical method Methods 0.000 claims description 5
- 230000007613 environmental effect Effects 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 9
- 238000000691 measurement method Methods 0.000 description 6
- 230000003287 optical effect Effects 0.000 description 4
- 230000001360 synchronised effect Effects 0.000 description 3
- 238000006073 displacement reaction Methods 0.000 description 2
- 239000011521 glass Substances 0.000 description 2
- 238000002032 lab-on-a-chip Methods 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011549 displacement method Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000005305 interferometry Methods 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Landscapes
- Measuring Volume Flow (AREA)
Abstract
本发明公开了一种显微粒子图像微流量测量仪,包括12位CCD相机(1),荧光滤色镜(2)、三棱镜(3)、激光器(4)、透镜组(5)、激光滤色镜(6)、显微物镜(7)、微流体注射泵(8)、光学玻璃旋转盘(9)、微流体芯片(10)、回收容器(11)、同步器(12)、计算机(13);本发明所采用的微流量测量技术属于非接触式方法。传感元件无须嵌入流体内部,避免了传统的接触式方法对微流体干涉造成的测量误差,不同于其它非接触式的间接方法,可直接测量流体流量参数,提高了测量的准确性;利用计算机图形学理论和方法,绘制三维空间速度矢量图,完成三维流场可视化构建,实现三维流场构建,达到精确测量微流量参数的目的。
Description
技术领域
本发明涉及微流量测量技术,特别是涉及一种显微粒子图像微流量测量仪及方法。
背景技术
微流体的主要特征是流速低、流量微小、多物理现象并存且相互耦合,而当微流体的体积流量落入微纳升/皮升尺度(1pL/s~1μL/s)范围内,存在诸多不为人知的物理现象。微流体微纳升/皮升尺度流量测量在国际上是一个难题,成熟技术少且多沿用接触式测量方法,难以实现流量的精确测量。“芯片实验室(Lab-on-a-chip)”技术,是微纳升/皮升尺度的微流体在生物、化学和医学领域中的重要应用,由于流速极低,流量的精确控制就成为“芯片实验室”实现样品的制备、反应以及合成控制的必要条件,因而对微流体流量测量提出了很高的精度要求。
目前广泛采用的微流量检测方法主要分为接触式和非接触式。
1)传统的测量方式以接触式测量为主,但是当流动尺度缩小至几十微米量级,传统流量测量技术就很难发挥应有作用,而对微纳升/皮升尺度流量测量则更是无能为力了。虽然近40年以来发展起来的基于MEMS技术的微流量传感器研究取得了很大进展,但由于绝大多数传感器必须嵌入流场内部,对流动本身产生了干扰和破坏,应用于微纳升/皮升等极限尺度流量测量时准确性降低,测量结果可信度不高。
2)最常用的非接触式测量方法为称重法,是根据一段时间累积的流体重量计算获得流量。但对于微纳升/皮升尺度流量测量,这显然是相当耗时的测量方法,而且在用于具有分支结构的微通道流量测量时出现了困难。随着光学及光学传感器测量方法的发展,非接触测量方法中又出现了红外光学干涉法、光纤位移法等微纳升/皮升尺度流量测量技术。但由于测量是基于不同流量下被测对象光学特性或光纤位移等参数变化的间接方法,存在诸多不确定的影响因素。
发明内容
本发明的目的在于针对现有技术的不足提供一种显微粒子图像微流量测量仪及方法。
本发明创造的目的是通过下述的技术方案实现的:
一种显微粒子图像微流量测量仪,包括12位CCD相机(1),荧光滤色镜(2)、三棱镜(3)、激光器(4)、透镜组(5)、激光滤色镜(6)、显微物镜(7)、微流体注射泵(8)、光学玻璃旋转盘(9)、微流体芯片(10)、回收容器(11)、同步器(12)、计算机(13);激光器(4)发出脉冲光束,依次经过透镜组(5)、激光滤色镜(6)、三棱镜(3),经三棱镜(3)全反射至显微物镜(7),并经显微物镜(7)、光学玻璃旋转盘(9)照射至被测的微流体芯片(10),微流体注射泵(8)以恒定流量向微流体芯片(10)微通道注射流体,微流体芯片(10)中的被照射荧光粒子产生波长为610nm的荧光,荧光依次经过光学玻璃旋转盘(9)、三棱镜(3)、荧光滤色镜(2)最终被12位的CCD相机(1)采集并形成了粒子图像,采集的粒子图像由计算机(13)进行处理,获取全流场不同流体层的粒子图像,通过速度分析算法计算微流体三维速度分布,然后进行求解;所述光学玻璃旋转盘(9)由直流电机(14)带动进行旋转,光学玻璃旋转盘(9)打磨成楔形,其圆周上的厚度分布为逐渐变化的,一端最薄另一端最厚。
所述的显微粒子图像微流量测量仪进行测量的方法,扫描开始时,直流电机(14)带动光学玻璃旋转盘(9)转动,微流体采集区域所对应的光学玻璃旋转盘(9)厚度不断增加,12位CCD相机(1)同步采集相应厚度情况下的微流体层的数据:当光学玻璃旋转盘(9)旋转半周,从光学玻璃旋转盘(9)最薄处旋转到最厚处时,完成一次扫描,采集的结果送至计算机构建三维流场;
为了构建三维流场,需要计算流体层二维平面速度;将通过扫描周期性地将采集了一定数量粒子图像对,通过系综相关算法求解得到对应的平均的二维平面速度;而对于存在测量平面法线方向上的速度分量的复杂微流动,采用测量得到的二维平面速度通过连续积分方程获得该方向的速度分布,积分常数通过上下壁面滑移边界条件进行确定;为了提高测量的可靠性和精度,对已获得的速度分布计算沿流向的所有截面的流量,得到一组流量测量值,将其平均值作为最终所需要的流量值,见式(1):
其中,u为流向速度分量;i表示沿流动方向的截面编号,j、k分别表示展向及景深方向的离散点编号;ΔAj,k为微元面积,计算过程嵌入速度场算法中;对于流速极低的皮升尺度流量测量,严格控制影响流动稳定性的外界环境因素,并利用多次测量结果的系综平均消除不稳定涨落现象;进而借助现有的计算机图形学理论和方法,绘制三维空间速度矢量图,完成三维流场可视化构建。
本发明与现有技术相比,具有的有益的效果是:
1、本发明所采用的微流量测量技术属于非接触式方法。传感元件无须嵌入流体内部,避免了传统的接触式方法对微流体干涉造成的测量误差,不同于其它非接触式的间接方法,可直接测量流体流量参数,提高了测量的准确性;
2、利用计算机图形学理论和方法,绘制三维空间速度矢量图,完成三维流场可视化构建,实现三维流场构建,达到精确测量微流量参数的目的。
附图说明
图1是显微粒子图像微流量测量仪示意图;
图2是旋转盘扫描机构示意图;
图3是光学玻璃旋转盘14的立体结构示意图;
图4是旋转盘扫描机构原理示意图;
1、12位CCD相机,2、荧光滤色镜,3、三棱镜,4、激光器,5、透镜组,6、激光滤色镜,7、显微物镜,8、微流体注射泵,9、光学玻璃旋转盘,10、微流体芯片,儿、回收容器,12、同步器,13、计算机,14、直流电机。
具体实施方式
以下结合具体实施例,对本发明进行详细说明。
如图1所示:显微粒子图像微流量测量仪包括12位CCD相机1,荧光滤色镜2,三棱镜3,激光器4,透镜组5,激光滤色镜6,显微物镜7,微流体注射泵8,光学玻璃旋转盘9,微流体芯片10,回收容器11,同步器12,计算机13;计算机13分别连接同步器12的一个端子及12位CCD相机1的一个端子,同步器12的其他三个端子分别连接光学玻璃旋转盘9,激光器4和12位CCD相机1的另一个端子,激光器4发出脉冲光束,依次经过透镜组5、激光滤色镜6、三棱镜3,经三棱镜3全反射至显微物镜7,并经显微物镜7、光学玻璃旋转盘9照射至被测的微流体芯片10,微流体注射泵8以恒定流量向微流体芯片10微通道注射流体,微流体芯片10中的被照射荧光粒子产生波长为610nm的荧光,荧光依次经过光学玻璃旋转盘9、三棱镜3、荧光滤色镜2最终被12位的CCD相机1采集并形成了粒子图像,采集的粒子图像由计算机13进行处理,获取全流场不同流体层的粒子图像,通过速度分析算法计算微流体三维速度分布,然后进行求解。
如图2、图3所示:光学玻璃旋转盘9由直流电机14带动进行旋转,其旋转可由同步器12对直流电机14进行同步控制。其中光学玻璃旋转盘9打磨成楔形,其圆周上的厚度分布为逐渐变化的,一端最薄另一端最厚。直流电机14采用同步“锁相”技术,驱动光学玻璃旋转盘9以设定角速度旋转,对流体跨层扫描,激光脉冲及CCD图像采集之间进行同步控制。
具体工作过程如下:
在进行微流量测量过程中,计算机13发出控制信号到同步器12,同步器12控制激光器4、12位CCD相机1和直流电机15进行同步工作。激光器4发出的激光经过透镜组5、激光滤色镜6、三棱镜3,经过三棱镜3反射至显微物镜7,透过光学玻璃旋转盘9的照射到待测量的微流体芯片10中的荧光粒子表面发生反射,荧光粒子发射出的波长为610nm,荧光依次经过光学玻璃旋转盘9、三棱镜3、荧光滤色镜2被12位CCD相机1采集成像得到粒子图像。
在扫描过程中,光学玻璃旋转盘9的不同厚度位置,测量的焦平面不同,从而12位CCD相机1采集的流体层也不同。其过程如图4所示,根据“菲涅尔”折射变焦原理,当两块玻璃的厚度相差△D的时候,所产生的焦平面相差△h,因此,当光学玻璃旋转盘9转动时,随着玻璃盘厚度的变化,测量的焦平面也在不断变化,直到采集到所有的流体层。我们设定光学玻璃旋转盘9最薄处为零位。扫描开始时,电机带动光学玻璃旋转盘9转动,微流体采集区域所对应的光学玻璃盘厚度不断增加,12位CCD相机1同步采集相应厚度情况下的微流体层的数据。当光学玻璃旋转盘9旋转半周,从光学玻璃旋转盘9最薄处旋转到最厚处时,完成一次扫描,采集的结果送至计算机构建三维流场。
为了构建三维流场,需要计算流体层二维平面速度。在每个扫描周期内拍得一对图像,这对图像拍摄的时间间隔是一定的,这对图像就记录了测量二维平面中粒子的位移情况。仅仅一对图像通过速度分析算法所得的结果往往会带来很大的误差,因此,需要将通过扫描周期性地将采集了一定数量粒子图像对,通过系综相关算法(系综相关算法是常用的速度分析算法)求解得到相应的平均的二维平面速度。而对于存在测量平面法线方向上的速度分量的复杂微流动,采用测量得到的二维平面速度通过连续积分方程获得该方向的速度分布,积分常数通过上下壁面滑移边界条件进行确定。为了提高测量的可靠性和精度,对已获得的速度分布计算沿流向的所有截面的流量,得到一组流量测量值,将其平均值作为最终所需要的流量值,见式(1)。
其中,u为流向速度分量;i表示沿流动方向的截面编号,j、k分别表示展向及景深方向的离散点编号:ΔAj,k为微元面积,计算过程嵌入速度场算法中。对于流速极低的皮升尺度流量测量,严格控制影响流动稳定性的外界环境因素,并利用多次测量结果的系综平均消除不稳定涨落现象。进而借助现有的计算机图形学理论和方法,绘制三维空间速度矢量图,完成三维流场可视化构建。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (2)
1.一种显微粒子图像微流量测量仪,其特征在于,包括12位CCD相机(1),荧光滤色镜(2)、三棱镜(3)、激光器(4)、透镜组(5)、激光滤色镜(6)、显微物镜(7)、微流体注射泵(8)、光学玻璃旋转盘(9)、微流体芯片(10)、回收容器(11)、同步器(12)、计算机(13);激光器(4)发出脉冲光束,依次经过透镜组(5)、激光滤色镜(6)、三棱镜(3),经三棱镜(3)全反射至显微物镜(7),并经显微物镜(7)、光学玻璃旋转盘(9)照射至被测的微流体芯片(10),微流体注射泵(8)以恒定流量向微流体芯片(10)微通道注射流体,微流体芯片(10)中的被照射荧光粒子产生波长为610nm的荧光,荧光依次经过光学玻璃旋转盘(9)、三棱镜(3)、荧光滤色镜(2)最终被12位的CCD相机(1)采集并形成了粒子图像,采集的粒子图像由计算机(13)进行处理,获取全流场不同流体层的粒子图像,通过速度分析算法计算微流体三维速度分布,然后进行求解;所述光学玻璃旋转盘(9)由直流电机(14)带动进行旋转,光学玻璃旋转盘(9)打磨成楔形,其圆周上的厚度分布为逐渐变化的,一端最薄另一端最厚。
2.根据权利要求1所述的显微粒子图像微流量测量仪进行测量的方法,其特征在于,扫描开始时,直流电机(14)带动光学玻璃旋转盘(9)转动,微流体采集区域所对应的光学玻璃旋转盘(9)厚度不断增加,12位CCD相机(1)同步采集相应厚度情况下的微流体层的数据;当光学玻璃旋转盘(9)旋转半周,从光学玻璃旋转盘(9)最薄处旋转到最厚处时,完成一次扫描,采集的结果送至计算机构建三维流场;
为了构建三维流场,需要计算流体层二维平面速度;将通过扫描周期性地将采集了一定数量粒子图像对,通过系综相关算法求解得到对应的平均的二维平面速度;而对于存在测量平面法线方向上的速度分量的复杂微流动,采用测量得到的二维平面速度通过连续积分方程获得该方向的速度分布,积分常数通过上下壁面滑移边界条件进行确定;为了提高测量的可靠性和精度,对已获得的速度分布计算沿流向的所有截面的流量,得到一组流量测量值,将其平均值作为最终所需要的流量值,见式(1):
其中,u为流向速度分量;i表示沿流动方向的截面编号,j、k分别表示展向及景深方向的离散点编号;ΔAj,k为微元面积,计算过程嵌入速度场算法中;对于流速极低的皮升尺度流量测量,严格控制影响流动稳定性的外界环境因素,并利用多次测量结果的系综平均消除不稳定涨落现象;进而借助现有的计算机图形学理论和方法,绘制三维空间速度矢量图,完成三维流场可视化构建。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310134340.9A CN103267543B (zh) | 2013-04-12 | 2013-04-12 | 显微粒子图像微流量测量仪及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310134340.9A CN103267543B (zh) | 2013-04-12 | 2013-04-12 | 显微粒子图像微流量测量仪及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103267543A CN103267543A (zh) | 2013-08-28 |
CN103267543B true CN103267543B (zh) | 2016-02-17 |
Family
ID=49011184
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310134340.9A Expired - Fee Related CN103267543B (zh) | 2013-04-12 | 2013-04-12 | 显微粒子图像微流量测量仪及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103267543B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103675333A (zh) * | 2013-12-08 | 2014-03-26 | 中国科学院过程工程研究所 | 一种实时测量微流体速度场的装置及方法 |
CN104950133B (zh) * | 2015-07-06 | 2017-12-15 | 浙江大学 | 微流体流速传感芯片、检测系统及检测方法 |
CN105891178B (zh) * | 2016-05-08 | 2018-06-12 | 重庆科技学院 | 一种集成化全内反射微流控芯片检测一体机使用方法 |
CN105784644B (zh) * | 2016-05-08 | 2018-07-31 | 重庆科技学院 | 一种恒温全内反射微流控芯片检测一体机使用方法 |
CN105784663B (zh) * | 2016-05-08 | 2018-06-12 | 重庆科技学院 | 一种恒温全内反射微流控芯片检测一体机 |
CN107015021B (zh) * | 2017-04-17 | 2019-05-07 | 青岛理工大学 | 一种摩擦副流体薄膜流速的三维观测方法 |
CN106896239B (zh) * | 2017-04-17 | 2019-08-20 | 青岛理工大学 | 一种定制滑块的摩擦副流体薄膜轴向流速成像测量装置 |
CN112898590B (zh) * | 2021-02-11 | 2022-06-21 | 鄂州篆光信息技术有限公司 | 一种金属有机框架微结构激光制造装置和方法 |
CN113238075B (zh) * | 2021-04-22 | 2023-02-14 | 哈尔滨工程大学 | 一种基于光纤光镊技术的流速计 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6653651B1 (en) * | 1998-12-09 | 2003-11-25 | Carl D. Meinhart | Micron resolution particle image velocimeter |
JP4811632B2 (ja) * | 2004-08-17 | 2011-11-09 | 横河電機株式会社 | 3次元共焦点顕微鏡システム |
CN102749475A (zh) * | 2012-07-26 | 2012-10-24 | 中国航空工业集团公司沈阳发动机设计研究所 | 一种单ccd相机三维粒子图像测速方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101052494B1 (ko) * | 2009-06-30 | 2011-07-29 | 한양대학교 산학협력단 | 이머젼 오일 기술을 적용한 마이크로 입자영상속도계 시스템 및 미세유동칩을 이용한 전자 인쇄용 잉크의 유동특성 측정방법 |
-
2013
- 2013-04-12 CN CN201310134340.9A patent/CN103267543B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6653651B1 (en) * | 1998-12-09 | 2003-11-25 | Carl D. Meinhart | Micron resolution particle image velocimeter |
JP4811632B2 (ja) * | 2004-08-17 | 2011-11-09 | 横河電機株式会社 | 3次元共焦点顕微鏡システム |
CN102749475A (zh) * | 2012-07-26 | 2012-10-24 | 中国航空工业集团公司沈阳发动机设计研究所 | 一种单ccd相机三维粒子图像测速方法 |
Non-Patent Citations (2)
Title |
---|
Development of a high-speed scanning micro PIV system using a rotating disc;K P Angele, Y Suzuki, J Miwa and N Kasagi;《MEASUREMENT SCIENCE AND TECHNOLOGY》;20060605;第17卷(第7期);第1639页倒数第1段及图1、2 * |
基于三维速度场构建的微流量测量方法研究;韩巍,王昊利;《仪器仪表学报》;20120229;第33卷(第2期);第415页左栏第1段、右栏倒数第1、2段,第416页左栏第1段及图2、3、4 * |
Also Published As
Publication number | Publication date |
---|---|
CN103267543A (zh) | 2013-08-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103267543B (zh) | 显微粒子图像微流量测量仪及方法 | |
CN103645341B (zh) | 全流场3d可视化测速方法 | |
CN101216422B (zh) | 液体浓度检测装置及检测方法 | |
CN104990500B (zh) | 一种检测金刚石对顶砧中物质体积和折射率的装置及方法 | |
CN108444897A (zh) | 基于光片照明与鞘流技术的免标记微流控细胞仪及方法 | |
CN207907955U (zh) | 一种基于线阵ccd的微流量计 | |
CN115932321B (zh) | 基于粒子图像测速的微观溶蚀可视化装置及方法 | |
CN101710131B (zh) | 离焦数字三维微流场荧光测试仪 | |
CN101660992A (zh) | 一种快速检测藻类细胞沉降速度的方法 | |
CN105136744A (zh) | 一种基于微流控芯片粒子捕获式的单粒子散射测量装置 | |
CN102636457B (zh) | 一种微量液体折射率的测量系统及测量方法 | |
CN107607047B (zh) | 一种在双光束光阱中校对图像传感器轴线失准的方法 | |
CN110160917B (zh) | 接触熔化过程中表面张力及反冲力的间接测量装置及方法 | |
CN105301282A (zh) | 一种流体三维速度场测量系统 | |
CN201780274U (zh) | 一种光学表面亚表层损伤测量装置 | |
CN114280327B (zh) | 基于光纤光镊的高灵敏加速度测量方法及传感器 | |
CN105548135A (zh) | 一种表面增强拉曼微流控芯片及包含该芯片的检测系统 | |
KR20120084450A (ko) | 하이브리드 마이크로/나노 piv 시스템을 이용한 마이크로 및 나노스케일의 3차원 유동 측정 방법 | |
CN105180820B (zh) | 一种透明毛细管内表面面型测试系统及测试方法 | |
Gerbino et al. | Local void fraction and fluid velocity measurements in a capillary channel with a single optical probe | |
CN201152850Y (zh) | 液体浓度检测装置 | |
CN108776087A (zh) | 一种用于测量气泡参数的相位临界角散射测量方法及装置 | |
CN112782059A (zh) | 一种无透镜成像系统中的准球细胞高度方向位移测量方法 | |
Vedyashkina et al. | Optical-electronic complex for investigation of the processes of heat and mass transfer by laser contactless method | |
CN102052278B (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: 20160217 Termination date: 20160412 |