CN102221535A - 三真空管气体折射率测量仪 - Google Patents

三真空管气体折射率测量仪 Download PDF

Info

Publication number
CN102221535A
CN102221535A CN 201110067880 CN201110067880A CN102221535A CN 102221535 A CN102221535 A CN 102221535A CN 201110067880 CN201110067880 CN 201110067880 CN 201110067880 A CN201110067880 A CN 201110067880A CN 102221535 A CN102221535 A CN 102221535A
Authority
CN
China
Prior art keywords
light beam
vacuum tube
gas
light
lambda
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
CN 201110067880
Other languages
English (en)
Other versions
CN102221535B (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.)
Tsinghua University
Original Assignee
Tsinghua University
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 Tsinghua University filed Critical Tsinghua University
Priority to CN201110067880.0A priority Critical patent/CN102221535B/zh
Publication of CN102221535A publication Critical patent/CN102221535A/zh
Application granted granted Critical
Publication of CN102221535B publication Critical patent/CN102221535B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

三真空管气体折射率测量仪属气体折射率精密测量技术领域,其特征在于以正交偏振双频激光器为光源,依次设置在该光源发射端光路上的平板分光镜、圆形玻璃板及贴在其中心的1/4波片、气体池及放置在其中的三支不同长度的真空管、一维位移台、角锥反射棱镜、偏振分光棱镜和平面反射镜,分别设置在上述分光棱镜和平面反射镜光路上的两个光电探测器,用于对上述两个光电探测器进行信号放大、采集、相位比较的信号处理单元,及用于系统控制与数据处理的计算机控制系统。利用该装置可对多种气体的折射率进行精确测量。

Description

三真空管气体折射率测量仪
技术领域
本发明属于气体折射率测量技术领域,特别适用于在光学精密测长、光学设计等技术领域中监测环境气体折射率。
背景技术
气体折射率是光学精密干涉测量、光学设计等领域中的重要参数。工作在空气环境中的仪器在设计以及数据处理时均需考虑空气折射率的影响。通常情况下空气折射率约为1.00027,且变化范围在10-6,因此一般将空气折射率近似为1。但在高精度的长度测量时,这一项近似所引入的误差将不容忽视,为实现纳米甚至亚纳米的测量不确定度,需要将空气折射率值精确测量到万分之一、十万分之一,甚至是更高精度。其次,某些光学仪器需要在保护气体(例如氮气等惰性气体)的环境下工作,为了得到满足设计要求的成像效果,需要对保护气体的折射率进行精确测量。此外,在大气环境监测、温室气体效应监测中,也需要对大气折射率进行精确测量。
目前常见的气体折射率测量方法可分为两类:间接测量法和直接测量法。其中,间接测量法是在由Edlen于1966年推导得出的经验公式基础上发展而来,它通过测量环境的温度、湿度、压力和CO2含量计算得到空气折射率。该方法的最大问题在于它只适用于空气折射率的测量,对于其他气体则完全不适用。直接测量法是根据折射率的定义对其进行直接测量。根据实现方法不同,又可细分为表面等离子波测量法、瑞利干涉仪法、抽/充气干涉测量法、梯形真空腔法、多波长干涉测量法和双真空室干涉测量法等。这些测量方法各有特点,同时也各自存在不足:要么测量准确度有限,要么测量时间过长,要么仪器加工困难,要么测量范围受限等。例如,清华大学殷纯永等人发明的双真空室双频测相空气折射率干涉仪(CN1095542C),该仪器使用两根真空管,这就大大限制了仪器的测量范围,使得测量范围有限,只能测量3×10-5范围内的空气折射率变化(例如,1.00024~1.00030)不适用于环境条件变化较为剧烈的场合,且仅适用于空气折射率的测量,不能测量除空气以外的其他气体。
发明内容
本发明提出的三真空管气体折射率测量仪,有效弥补了上述方法的不足,它的主要特点为:气体折射率的测量范围大、可测多种气体、测量速度快,测量准确度高等。
本发明提出了一种基于三真空管的气体折射率测量仪,该测量仪测量范围更大,可达4×10-4,可满足不同场合下气体折射率的测量;另外,本发明不但可用于测量空气折射率,更可用于测量其他气体(如氮气,二氧化碳等)的折射率。另外,本发明还具有测量速度快、测量准确度高的优点。
本发明的目的在于提出一种新型的气体折射率测量方法和装置,有效解决现有方法的不足,该发明可对环境气体的折射率进行精确、实时监控,从而有望在光学精密测长、光学设计和大气监测中获得应用。
根据本发明的一个方面,提供了一种气体折射率测量仪,其特征在于包括:
平板分光镜,用于把来自一个正交偏振双频激光光源的光分成正交偏振的第一光束和第二光束,
气体池,用于容纳待测气体,
并行设置在所述气体池中的第一真空管、第二真空管、第三真空管,
角锥棱镜,用于反射所述第一光束和第二光束,其中所述第二光束比所述第一光速更靠近所述角锥棱镜的轴线。
附图说明
图1是根据本发明的一个实施例的三真空管气体折射率测量仪的系统图。
图2是根据本发明的一个实施例的圆形玻璃板(31)和1/4波片的左视图。
图3是根据本发明的一个实施例的等效合成波长链构造过程示意图。
图4是根据本发明的一个实施例的利用等效合成波长链进行折射率测量的原理图。
图5是根据本发明的一个实施例的气体折射率测量的流程图。
具体实施方式
本发明的原理如图1所示。图1中箭头和圆环分别表示由激光器(1)发出的一束偏振方向平行和垂直于纸面的两个正交偏振光,它们的频率分别为f1、f2,经平板分光镜(2)分为两束,分别为第一光束(101)和第二光束(102),其中,第一光束(101)透过圆形玻璃板(31)和气体池(4)的观察窗(41)后沿第一真空管(5)的外部传播并经观察窗(42)出射,,第二光束(102)透过粘贴在所述圆形玻璃板(31)中心位置处的1/4波片(32)(圆形玻璃板(31)和1/4波片的侧视图如图2所示)和气体池(4)的第一观察窗(41)后沿第一真空管(5)的内部传播,并经第二观察窗(42)透射。第一光束和第二光束经角锥棱镜(9)反射后,再次透过第二观察窗(42)并通过第一真空管(5)的外部和内部,经第一观察窗(41)透射后,在平板分光镜(2)处合光成第三光束(103),产生干涉。通过第一真空管(5)内部的第二光束(102)两次通过圆形玻璃板(3)上的1/4波片(32),偏振方向将旋转90°,通过第一真空管(5)外部的第一光束(101)不经过1/4波片,因而偏振方向不发生改变。因此在平板分光镜(2)处通过第一真空管(5)内外的不同频率的第一光束(101)和第二光束(102)将产生干涉并形成第三光束(103)。利用偏振分光棱镜(10)和平板反射镜(11)将第三光束(103)的两个偏振态分开,即得到偏振态垂直于纸面的第四光束(104)和偏振态平行于纸面的第五光束(105),分别用光电探测器(12)、(13)对第四光束(104)和第五光束(105)进行探测,再经信号处理单元(14)对光电探测器(12)和光电探测器(13)的探测信号进行相位比较后,通过计算即可得到气体池(4)中待测气体的折射率。
探测器(12)、(13)探测到的两路干涉信号可表示为
Figure BDA0000051288760000031
式中I01、I02为干涉信号强度,
Figure BDA0000051288760000032
为第一光束(101)和第二光束(102)的初始相位,
Figure BDA0000051288760000033
为第一真空管(5)内外光束的光程差引入的干涉信号相位差。探测器(12)探测到的是f1激光通过第一真空管(5)外、f2激光通过第一真空管(5)内形成的干涉信号,探测器(13)探测到的是f2激光通过第一真空管(5)外、f1激光通过第一真空管(5)内形成的干涉信号。探测器(12)和探测器(13)探测到的信号对应着同一个光程差,因此干涉信号表达式中含有的相位
Figure BDA0000051288760000034
大小相等但符号相反。假设气体池(4)中待测气体的折射率为n,则相位差与气体折射率之间的关系可表示为
式(2)中L为第一真空管(1)的长度,λ为激光器(1)的工作波长。
如果测得探测器(12)和(13)的相位差且将它表示成干涉信号整数级次N和小数ε的和,则(2)式可变为:
( n - 1 ) = ( N + ϵ ) · λ 4 L - - - ( 3 )
定义
Figure BDA0000051288760000038
则(3)式变为:
(n-1)=(N+ε)·λs    (4)
由(4)式可知,对于固定的激光光源和固定长度的真空管,λs是常量,通过测量它的干涉级次的整数和小数部分,即可得到待测气体的折射率。我们称λs为“等效合成波长”。
等效合成波长链的构建过程如图3所示。第二真空管(6)和第三真空管(7)仅在长度上与第一真空管(5)有差别,设第一真空管(5)、第二真空管(6)和第三真空管(7)的长度分别为L1、L2和L3,且L1>L2>L3。三支真空管平行放置在一维平移台(8)上,通过控制一维平移台(8)的运动可以依次将三支真空管送入第一光束(101)和第二光束(102)的光路中,需要注意的是,每次有且仅有一支真空管在光路中。利用一维位移台(8)依次将第二真空管(6)、第三真空管(7)送入光路中,重复第一真空管(5)在光路中的测量过程,则可得到三个与(4)式类似的等式,等效合成波长分别为λs1、λs2和λs3,我们称之为零级等效合成波长,并用λs1 0、λs2 0和λs3 0重新表示它们(上脚标表示级次,下脚标表示次序),如式(5)所示:
λ S 1 0 = λ / 4 L 1 λ S 2 0 = λ / 4 L 2 λ S 3 0 = λ / 4 L 3 - - - ( 5 )
将式(5)中三个等式两两相减,可以得到两个类似于(4)式的等式,此时的等效合成波长称为一级等效合成波长,记为λs1 1、λs2 1,如式(6)所示:
λ S 1 1 = λ 4 ( L 1 - L 2 ) λ S 2 1 = λ 4 ( L 2 - L 3 ) - - - ( 6 )
利用同样的方法(将式(6)中两等式相减)还可得到二级等效合成波长λs1 2,如式(7)所示:
λ S 1 2 = λ 4 ( L 1 - 2 L 2 + L 3 ) - - - ( 7 )
通过上述处理即可构建一个如图3所示的等效合成波长链。
实际进行气体折射率测量时,与上述合成波长链的构造过程刚好相反,其过程如图4所示(图中int表示对括号内的函数做取整运算)。记三支真空管分别在光路中时测得的相位差分为Ns1 0s1 0、Ns2 0s2 0和Ns3 0s3 0,则根据定义,一级和二级等效合成波长的小数部分分别为:
ϵ s 1 1 = ϵ s 1 0 - ϵ s 2 0 ϵ s 2 1 = ϵ s 2 0 - ϵ s 3 0 ϵ S 1 2 = ϵ s 1 0 - 2 ϵ s 2 0 + ϵ s 3 0 - - - ( 8 )
实际中直接相位测量通常只能得到小数部分ε,不能得到整数N。然而,通过上述等效合成波长链的构造过程可知,高级次的等效合成波长要远远大于低级次的等效合成波长。根据(4)式可知,当真空管长度选择恰当时,可构造出合适的高级次等效合成波长(例如λs1 2),使得Ns1 2=0,从而利用直接测得的小数部分即可初步确定折射率n。利用由λs1 2得到的折射率作为初值,结合级间过渡条件(由高级次向低级次传递时保证N的单值性),可得到更精确的折射率值,当传递到最低级等效合成波长时,即可得到最高准确度的n。
激光光源可采用真空波长为632.991372nm、频差约为2.2MHz的正交偏振He-Ne激光器,三支真空管的设计长度分别为L1=165mm、L2=158mm和L3=151.5mm,实际长度可通过万能工具显微镜测量,长度测量不确定度约为2μm。使用的光电探测器是带宽3MHz的PIN管,相位计的测相准确度优于1°。
测量流程图如图5所示,详细测量过程描述如下:
A)根据待测气体的测量范围,选择三支长度合适的真空管,放置并固定在一维位移台(8)上;
B)在气体池(4)中充入待测气体;
C)利用一维位移台(8)将所述第一真空管(5)送入由第一光束(101)和第二光束(102)组成的光路,利用信号处理单元(14)测量得到相位差的小数部分εs1 0;重复上述过程,可得到当第二真空管(6)和第三真空管(7)在光路中时的相位差的小数部分εs2 0和εs3 0
D)利用步骤(3)测得的εs1 0、εs2 0和εs3 0,根据公式(8)可计算得到等效合成波长对应的相位差εs1 1、εs2 1和εS1 2
E)利用等效合成波长理论,根据(5)、(6)、(7)式计算得到等效合成波长,利用图(4)即可依次计算得到2级等效合成波长测得的折射率1级等效合成波长测得的折射率
Figure BDA0000051288760000052
0级等效合成波长测得的折射率
Figure BDA0000051288760000053
其中
Figure BDA0000051288760000054
就是待测气体折射率的最终测量结果。
本发明的特点包括:
(1)测量范围大;通过合理选择真空管长度,气体折射率的测量范围可达3×10-4
(2)测量对象范围广;通过在气体池中充入不同气体,可对多种气体的折射率进行测量,例如,空气,氮气,二氧化碳等。
(3)测量准确度高;合理选择最长真空管的长度,可使测量仪的相对测量不确定度优于10-8
(4)测量速度快;单次测量时间可控制在1分钟之内,适用于在线测量。

Claims (7)

1.气体折射率测量仪,其特征在于包括:
平板分光镜(2),用于把来自一个正交偏振双频激光光源(1)的光分成正交偏振的第一光束(101)和第二光束(102),
气体池(4),用于容纳待测气体,
并行设置在所述气体池(4)中的第一真空管(5)、第二真空管(6)、第三真空管(7),
角锥棱镜(9),用于反射所述第一光束(101)和第二光束(102),其中所述第二光束(102)比所述第一光速(101)更靠近所述角锥棱镜(9)的轴线。
2.根据权利要求1的气体折射率测量仪,其特征在于进一步包括:
平移台(8),所述第一真空管(5)、第二真空管(6)、第三真空管(7)被并行设置在所述平移台(8)上,从而能够通过所述平移台(8)的移动而把所述第一真空管(5)、第二真空管(6)、第三真空管(7)中选定的一个置于所述第一光束(101)和第二光束(102)的光路中,使得第一光束(101)透过气体池(4)的第一观察窗(41)和第二观察窗(42)后沿第一真空管(5)的外部传播,同时使得第二光束(102)透过气体池(4)的第一观察窗(41)和第二观察窗(42)后沿第一真空管(5)的内部传播。
3.根据权利要求2的气体折射率测量仪,其特征在于进一步包括:
设置在第二光束(102)的光路中的1/4波片(32),以及设置在第一光束(101)和第二光束(102)光路中的圆形玻璃板(31),所述1/4波片(32)通过胶合的方式紧贴在圆形玻璃板(31)的中心位置。
4.根据权利要求3的气体折射率测量仪,其特征在于进一步包括:
偏振分光棱镜(10);
将从所述气体池(4)出射的所述第一光束(101)和第二光束(102)相干形成的出射光第三光束(103)分成两部分,第一部分光(104)沿原光路出射,第二部分光(105)沿与入射光路垂直的方向反射。
5.根据权利要求4的气体折射率测量仪,其特征在于进一步包括:
平面反射镜(11),用于反射第一部分光(104)。
6.根据权利要求5的气体折射率测量仪,其特征在于进一步包括:
设置在所述第一部分光(104)的光路上的第一光电探测器(12),
设置在所述第二部分光(105)的光路上的第二光电探测器(13),
信号处理单元(14),用于处理第一光电探测器(12)和第二光电探测器(13)所探测的信号。
7.气体折射率测量方法,其特征在于该方法基于如权利要求1-6所述的气体折射率测量仪,并包括:
A)根据待测气体的测量范围,选择三支长度合适的真空管,放置并固定在一维位移台(8)上;
B)在气体池(4)中充入待测气体;
C)利用一维位移台(8)将所述第一真空管(5)送入由第一光束(101)和第二光束(102)组成的光路,利用信号处理单元(14)测量得到相位差的小数部分εs1 0;重复上述过程,可得到当第二真空管(6)和第三真空管(7)在光路中时的相位差的小数部分εs2 0和εs3 0
D)利用步骤C)测得的εs1 0、εs2 0和εs3 0,根据公式
ϵ s 1 1 = ϵ s 1 0 - ϵ s 2 0 ϵ s 2 1 = ϵ s 2 0 - ϵ s 3 0 ϵ S 1 2 = ϵ s 1 0 - 2 ϵ s 2 0 + ϵ s 3 0
计算得到等效合成波长对应的相位差εs1 1、εs2 1和εS1 2
E)根据
λ S 1 0 = λ / 4 L 1 λ S 2 0 = λ / 4 L 2 λ S 3 0 = λ / 4 L 3 - - - ( 5 )
λ S 1 1 = λ 4 ( L 1 - L 2 ) λ S 2 1 = λ 4 ( L 2 - L 3 ) - - - ( 6 )
λ S 1 2 = λ 4 ( L 1 - 2 L 2 + L 3 ) - - - ( 7 )
计算得到等效合成波长,利用关系
Figure FDA0000051288750000031
依次计算得到2级等效合成波长测得的折射率
Figure FDA0000051288750000032
1级等效合成波长测得的折射率0级等效合成波长测得的折射率
Figure FDA0000051288750000034
其中
Figure FDA0000051288750000035
就是待测气体折射率的最终测量结果。
CN201110067880.0A 2011-03-21 2011-03-21 三真空管气体折射率测量仪 Expired - Fee Related CN102221535B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201110067880.0A CN102221535B (zh) 2011-03-21 2011-03-21 三真空管气体折射率测量仪

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201110067880.0A CN102221535B (zh) 2011-03-21 2011-03-21 三真空管气体折射率测量仪

Publications (2)

Publication Number Publication Date
CN102221535A true CN102221535A (zh) 2011-10-19
CN102221535B CN102221535B (zh) 2014-07-09

Family

ID=44778146

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201110067880.0A Expired - Fee Related CN102221535B (zh) 2011-03-21 2011-03-21 三真空管气体折射率测量仪

Country Status (1)

Country Link
CN (1) CN102221535B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103076304A (zh) * 2013-01-05 2013-05-01 浙江理工大学 调制式激光干涉空气折射率测量方法及装置
CN103196868A (zh) * 2013-04-10 2013-07-10 上海理工大学 光刻胶折射率的测定方法
CN108254339A (zh) * 2017-12-22 2018-07-06 北京航天计量测试技术研究所 一种大量程气体折射率实时测量系统
CN108318419A (zh) * 2017-12-20 2018-07-24 北京航天计量测试技术研究所 一种用于空气折射率实时测量的真空腔

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN86107252A (zh) * 1986-10-21 1988-05-04 北京机械工业管理学院分部 空气折射率测量装置
DE3911471A1 (de) * 1989-04-08 1990-10-11 Kerner Anna Heterodynes verfahren zur refraktometrischen messung und stabilisierung der wellenlaenge
CN1103161A (zh) * 1993-11-26 1995-05-31 清华大学 高精度激光腔变位移/折射率测量方法及其装置
CN1174327A (zh) * 1997-07-23 1998-02-25 中国计量科学研究院 双瑞利空气折射率干涉仪
CN1260484A (zh) * 2000-02-25 2000-07-19 清华大学 双真空室双频测相空气折射率干涉仪

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN86107252A (zh) * 1986-10-21 1988-05-04 北京机械工业管理学院分部 空气折射率测量装置
DE3911471A1 (de) * 1989-04-08 1990-10-11 Kerner Anna Heterodynes verfahren zur refraktometrischen messung und stabilisierung der wellenlaenge
CN1103161A (zh) * 1993-11-26 1995-05-31 清华大学 高精度激光腔变位移/折射率测量方法及其装置
CN1174327A (zh) * 1997-07-23 1998-02-25 中国计量科学研究院 双瑞利空气折射率干涉仪
CN1260484A (zh) * 2000-02-25 2000-07-19 清华大学 双真空室双频测相空气折射率干涉仪

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《自然科学进展》 20040331 张书练等 正交线偏振激光器原理与应用(II)--物理现象研究 全文 1-7 第14卷, 第3期 *
《计量学报》 20040731 陈强华 等 基于等效合成波方法的双真空管空气折射率测量仪 第203-206页,附图1 1-7 第25卷, 第3期 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103076304A (zh) * 2013-01-05 2013-05-01 浙江理工大学 调制式激光干涉空气折射率测量方法及装置
CN103076304B (zh) * 2013-01-05 2015-01-14 浙江理工大学 调制式激光干涉空气折射率测量方法及装置
CN103196868A (zh) * 2013-04-10 2013-07-10 上海理工大学 光刻胶折射率的测定方法
CN108318419A (zh) * 2017-12-20 2018-07-24 北京航天计量测试技术研究所 一种用于空气折射率实时测量的真空腔
CN108254339A (zh) * 2017-12-22 2018-07-06 北京航天计量测试技术研究所 一种大量程气体折射率实时测量系统

Also Published As

Publication number Publication date
CN102221535B (zh) 2014-07-09

Similar Documents

Publication Publication Date Title
US11397149B2 (en) Laser radar system apparatus for multi-wavelength measurement of atmospheric carbon dioxide concentration and vertical aerosol profile
US8665452B2 (en) Method and apparatus for measuring the refractive index of air based on the laser synthetic wavelength interferometry
CN102183486B (zh) 一种基于光频梳的气体折射率测量仪及其测量方法
CN104655025B (zh) 激光干涉波长杠杆式绝对距离测量方法与装置
CN104749580B (zh) 绝对探测大气温度的转动拉曼激光雷达系统及探测方法
CN103575701A (zh) 基于频域oct的透明材料折射率及厚度测量方法和装置
CN103954589A (zh) 一种光学材料折射率的精密测量装置及方法
CN102221535B (zh) 三真空管气体折射率测量仪
CN102735646A (zh) 透明介质折射率的测量装置及测量方法
CN102323237B (zh) 一种空气折射率快速高精度绝对测量装置及其测量方法
CN107631993A (zh) 一种基于水汽残差谱分析的红外光谱识别方法
CN104807781B (zh) 一种基于色散干涉法的空气折射率测量装置及测量方法
CN107064067A (zh) 一种双干涉仪的空气折射率廓线的测量系统及方法
CN106500911A (zh) 一种基于气体吸收谱线压力展宽效应的压力计校准方法
CN102636457B (zh) 一种微量液体折射率的测量系统及测量方法
US10190977B2 (en) Method of measurement of multilayer structures
Lin et al. Discovery of New Lines in the R 9 Multiplet of the 2 v 3 Band of CH 12 4
CN102288104A (zh) 六轴4细分干涉仪
CN102200581B (zh) 使用单色仪分光的高精度的水汽拉曼系统及定标方法
CN106537086B (zh) 用渐变孔径光学积分法测量单模光纤模场直径的方法
CN104949938A (zh) 一种基于游标效应的马赫曾德调制型谐振腔传感器
CN205003080U (zh) 基于谱域干涉仪的折射率和厚度同步测量系统
CN101435880A (zh) 高层大气风场静态宽场消色差消温差多方向实时同时探测法
US10761021B2 (en) Apparatus and method for measurement of multilayer structures
CN106226268A (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: 20140709

Termination date: 20150321

EXPY Termination of patent right or utility model