CN102590092A - 用于激光吸收光谱技术的吸收光程延长的装置及方法 - Google Patents
用于激光吸收光谱技术的吸收光程延长的装置及方法 Download PDFInfo
- Publication number
- CN102590092A CN102590092A CN2012100583608A CN201210058360A CN102590092A CN 102590092 A CN102590092 A CN 102590092A CN 2012100583608 A CN2012100583608 A CN 2012100583608A CN 201210058360 A CN201210058360 A CN 201210058360A CN 102590092 A CN102590092 A CN 102590092A
- Authority
- CN
- China
- Prior art keywords
- detector
- absorption
- light path
- porosint
- core
- 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
Links
- 238000010521 absorption reaction Methods 0.000 title claims abstract description 53
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000005516 engineering process Methods 0.000 title claims abstract description 16
- 230000003287 optical effect Effects 0.000 title abstract description 9
- 238000001285 laser absorption spectroscopy Methods 0.000 title abstract 3
- 239000007789 gas Substances 0.000 claims description 45
- 238000000862 absorption spectrum Methods 0.000 claims description 15
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 12
- PNEYBMLMFCGWSK-UHFFFAOYSA-N Alumina Chemical compound [O-2].[O-2].[O-2].[Al+3].[Al+3] PNEYBMLMFCGWSK-UHFFFAOYSA-N 0.000 claims description 6
- 229910052757 nitrogen Inorganic materials 0.000 claims description 6
- GWEVSGVZZGPLCZ-UHFFFAOYSA-N Titan oxide Chemical compound O=[Ti]=O GWEVSGVZZGPLCZ-UHFFFAOYSA-N 0.000 claims description 4
- MCMNRKCIXSYSNV-UHFFFAOYSA-N Zirconium dioxide Chemical compound O=[Zr]=O MCMNRKCIXSYSNV-UHFFFAOYSA-N 0.000 claims description 4
- 239000000463 material Substances 0.000 claims description 4
- 230000003321 amplification Effects 0.000 claims description 3
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 3
- 239000004408 titanium dioxide Substances 0.000 claims description 2
- 238000001514 detection method Methods 0.000 abstract description 6
- 239000011148 porous material Substances 0.000 abstract 1
- 230000001105 regulatory effect Effects 0.000 description 3
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000012925 reference material Substances 0.000 description 1
- 238000001028 reflection method Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
用于激光吸收光谱技术的吸收光程延长的装置及方法,属于气体测量领域,本发明为解决现有的激光吸收光谱技术吸收光程延长方法所导致的气体容器体积大、检测响应速度慢的问题。本发明装置包括气体池、多孔材料芯、探测器、一维平移台和放大器,吸收光程延长方法为:在一维平移台带着探测器移动的方向上建立一维坐标系,该一维坐标系的原点为探测器初始位置,驱动一维平移台带着探测器自初始位置点开始移动了n个位置点,每个位置点的坐标为Xi,测量探测器在第i个位置对应的吸收光程为Leff(Xi),根据n个位置点及其对应吸收光程,通过二次多项式拟合获得吸收光程Leff(Xi)与探测器位置坐标x的关系式f(x)=Leff(x),进而通过调整探测器的位置来达到预期的吸收光程。
Description
技术领域
本发明涉及用于激光吸收光谱技术的吸收光程延长的装置及方法,属于气体测量领域。
背景技术
激光吸收光谱技术已经成为气体物质含量检测的一种常用的技术手段,为了实现痕量气体的检测,一个常用的直接方法是增加光通过被测气体的光程长度,从而产生更强的吸收来达到提高信噪比的目的。目前用于延长吸收光程的方法主要包括多通池多次反射法、高反射腔增强法、积分球漫反射法等。这些方法尽管可以获得几十倍至几万倍的吸收光程放大倍数,但是基于这些方法的气体容器体积大,不利于测量系统的小型化,同时降低了气体的更换速度,最终导致整个系统的检测响应速度慢。
发明内容
本发明目的是为了解决现有的激光吸收光谱技术吸收光程延长方法所导致的气体容器体积大、检测响应速度慢的问题,提供了一种用于激光吸收光谱技术的吸收光程延长的装置及方法。
本发明所述用于激光吸收光谱技术的吸收光程延长的装置,它包括气体池、多孔材料芯、探测器、一维平移台和放大器,气体池内设置有多孔材料芯,激光光束入射至气体池,该激光光束穿过多孔材料芯后入射至探测器的光敏面,探测器由一维平移台带动沿探测器的光敏面所在平面做一维移动,探测器将探测到的光信号转换成电信号,并通过放大器放大后输出。
基于上述装置的吸收光程延长方法为:在一维平移台带着探测器移动的方向上建立一维坐标系,该一维坐标系的原点为探测器初始位置,初始位置点的坐标X0=0,驱动一维平移台带着探测器自初始位置点开始移动了n个位置点,每个位置点的坐标为Xi,i=0,1,2,...,n,n为自然数,且n=7~20,
测量探测器在第i个位置对应的吸收光程为Leff(Xi),根据n个位置点及其对应吸收光程,通过二次多项式拟合获得吸收光程Leff(Xi)与探测器位置坐标x的关系式f(x)=Leff(x),
进而通过调整探测器的位置来达到预期的吸收光程。
在实际应用中通过调整探测器的位置来调整有效吸收光程的大小,实施方法为根据f(x)=Leff(x)获得的多项式,得到其反函数x(Leff),将所需要得到的有效吸收光程Leff=Ln代入此反函数中即可获得探测器所对应的坐标位置x(Ln),通过调节一维平移台来将探测器移到此位置点即可获得预期的有效吸收光程。
本发明的优点:本发明在采用小体积气体池的情况下仍然可以获得较大的吸收光程放大倍数,也就是在保证吸收光谱测量系统小型化的前提下显著提高测量结果的信噪比,从而提高气体检测的灵敏度。此外,可以在不对气体池本身做任何变动的情况下,通过移动探测器的位置来调节吸收光程的放大倍数。
附图说明
图1是本发明所述用于激光吸收光谱技术的吸收光程延长的装置的结构示意图;
图2是多孔材料芯的形状为长方体时,该长方体的激光光束入射面的主视图;
图3是图2的侧视图;
图4是多孔材料芯的形状为圆柱体时,该圆柱体的激光光束入射面的主视图;
图5是图4的侧视图;
图6是激光光束入射至多孔材料芯表面入射角度示意图。
具体实施方式
具体实施方式一:下面结合图1说明本实施方式,本实施方式所述用于激光吸收光谱技术的吸收光程延长的装置,它包括气体池1、多孔材料芯2、探测器3、一维平移台4和放大器5,气体池1内设置有多孔材料芯2,激光光束入射至气体池1,该激光光束穿过多孔材料芯2后入射至探测器3的光敏面,探测器3由一维平移台4带动沿探测器3的光敏面所在平面做一维移动,探测器3将探测到的光信号转换成电信号,并通过放大器5放大后输出。
具体实施方式二:本实施方式对实施方式一作进一步说明,多孔材料芯2的组成物质为氧化铝、氧化锆或氧化钛中的一种或几种,材料孔隙率大于30%,材料孔平均直径小于10μm。
具体实施方式三:下面结合图2和图3说明本实施方式,本实施方式对实施方式一作进一步说明,多孔材料芯2的形状为长方体,该长方体的激光光束入射面的长度为a,宽度为b;垂直于激光光束入射面的长方体的厚度为d,且满足a≥b>3d。
具体实施方式四:下面结合图4和图5说明本实施方式,本实施方式对实施方式一作进一步说明,多孔材料芯2的形状为圆柱体,该圆柱体的激光光束入射圆面的直径为e,垂直于激光光束入射面的圆柱体的厚度为f,且满足e>3f。
具体实施方式五:基于实施方式一所述装置的吸收光程延长的方法为:在一维平移台4带着探测器3移动的方向上建立一维坐标系,该一维坐标系的原点为探测器3初始位置,初始位置点的坐标X0=0,驱动一维平移台4带着探测器3自初始位置点开始移动了n个位置点,每个位置点的坐标为Xi,i=0,1,2,...,n,n为自然数,且n=7~20,
测量探测器3在第i个位置对应的吸收光程为Leff(Xi),根据n个位置点及其对应吸收光程,通过二次多项式拟合获得吸收光程Leff(Xi)与探测器3位置坐标x的关系式f(x)=Leff(x),
进而通过调整探测器3的位置来达到预期的吸收光程。
在实际应用中通过调整探测器的位置来调整有效吸收光程的大小,实施方法为根据f(x)=Leff(x)获得的多项式,得到其反函数x(Leff),将所需要得到的有效吸收光程Leff=Ln代入此反函数中即可获得探测器3所对应的坐标位置x(Ln),通过调节一维平移台4来将探测器3移到此位置点即可获得预期的有效吸收光程。
具体实施方式六:本实施方式对实施方式五作进一步说明,测出任一位置点Xi时的吸收光程Leff(Xi)的过程为:
测量探测器3在第i个位置对应的吸收光程Leff(Xi)的过程为:
步骤一、探测器3固定在该位置点,在气体池1内充入缓冲气氮气,调整激光光束入射至气体池1,放大器5将探测器3采集的光强信号放大后输出,根据该放大后的信号计算获得多孔材料芯2散射出的无吸收光强度I0;
步骤二、在气体池1内充入已知浓度的样品气体,调整激光光束入射至气体池1,放大器5将探测器3采集的光强信号放大后输出,根据该放大后的信号计算获得吸收光强度It;
步骤三、根据步骤一获取的氮气光强度I0和步骤二获取的样品光强度It按公式
计算获得探测器3位于该位置时对应的光程Leff(Xi),
式中,N为样品气体的浓度,σ为样品气体的吸收截面。该截面可由光谱数据库查得。
工作原理:激光光束入射至多孔材料芯2,经过入射表面多孔结构散射后光子会延不同的反向散射开来,一些光子以不同的散射方向进入到多孔材料芯2内部,在内部多孔结构继续的散射作用下,沿着不同的方向前进,最后在多孔材料芯2的不同表面的不同位置处逸出。由于光子在多孔材料芯2内部不是沿着直线传播的,而是在多孔结构的散射作用下曲折式前进,相比于直线式的传播方式所经历的吸收光程可以扩大几倍至几万倍。在气体池内充入气体后,气体会渗入到多孔材料芯2内部并很快达到平衡状态,使得气体池1内部各个位置的气体浓度是一致的。探测器3所接收到的光子是经由不同路径抵达的,因此所携带的吸收信号大小各不相同。探测器3的响应时间不足以分辨单个光子的独立吸收过程,因此最终探测到的光强I是一种平均效果,由公式所得到的光程是一种表征所探测到的光子集体吸收行为的有效光程。
有效吸收光程不是简单的对各个光子吸收光程取平均值,对它的精确理论计算过程是一个复杂的过程,这里是采用标准物定标的方法获得的。先在多孔材料芯2内充入不含气体样品的缓冲气氮气,测得无吸收时的光强I0;再将多孔材料芯2内的气体置换为具有已知浓度N的样品气体,测得此时经过气体吸收后的光强It,则根据公式可以获得此位置探测器3所测量得到的光子吸收对应的有效光程。
当探测器3处于不同位置的时候,接收到的光子所经历的吸收路径是不同的,相应的有效吸收光程也是不同的。理论和实践均表明,在逐步改变探测器3的位置的情况下,所获得的有效光程是按照确定规律变化的。以探测器3的初始位置为原点建立探测器位置的一维坐标系,通过一维平移台4不断调整探测器3至不同的位置Xi,获得相应的有效光程Leff(Xi),通过二次多项式拟合获得有效吸收光程Leff(Xi)与探测器3位置x的关系式f(x)=Leff(x)。在实际应用中可以根据需要通过调整探测器3的位置来调整有效吸收光程的大小。
具体实施方式七:本实施方式对实施方式六作进一步说明,入射至气体池1的激光光束与多孔材料芯2入射面的法线所成的角度θ范围为0~45度。
Claims (7)
1.用于激光吸收光谱技术的吸收光程延长的装置,其特征在于,它包括气体池(1)、多孔材料芯(2)、探测器(3)、一维平移台(4)和放大器(5),气体池(1)内设置有多孔材料芯(2),激光光束入射至气体池(1),该激光光束穿过多孔材料芯(2)后入射至探测器(3)的光敏面,探测器(3)由一维平移台(4)带动沿探测器(3)的光敏面所在平面做一维移动,探测器(3)将探测到的光信号转换成电信号,并通过放大器(5)放大后输出。
2.根据权利要求1所述用于激光吸收光谱技术的吸收光程延长的装置,其特征在于,多孔材料芯(2)的组成物质为氧化铝、氧化锆或氧化钛中的一种或几种,材料孔隙率大于30%,材料孔平均直径小于10μm。
3.根据权利要求1所述用于激光吸收光谱技术的吸收光程延长的装置,其特征在于,多孔材料芯(2)的形状为长方体,该长方体的激光光束入射面的长度为a,宽度为b;垂直于激光光束入射面的长方体的厚度为d,且满足a≥b>3d。
4.根据权利要求1所述用于激光吸收光谱技术的吸收光程延长的装置,其特征在于,多孔材料芯(2)的形状为圆柱体,该圆柱体的激光光束入射圆面的直径为e,垂直于激光光束入射面的圆柱体的厚度为f,且满足e>3f。
5.基于权利要求1所述装置的吸收光程延长的方法,其特征在于,吸收光程延长的方法为:在一维平移台(4)带着探测器(3)移动的方向上建立一维坐标系,该一维坐标系的原点为探测器(3)初始位置,驱动一维平移台(4)带着探测器(3)自初始位置点开始移动了n个位置点,每个位置点的坐标为Xi,i=0,1,2,...,n,其中初始位置点的坐标X0=0,n为自然数,且n=7~20,
测量探测器(3)在第i个位置对应的吸收光程为Leff(Xi),根据n个位置点及其对应吸收光程,通过二次多项式拟合获得吸收光程Leff(Xi)与探测器(3)位置坐标x的关系式f(x)=Leff(x),
进而通过调整探测器(3)的位置来达到预期的吸收光程。
6.根据权利要求5所述用于激光吸收光谱技术的吸收光程延长的方法,其特征在于,测量探测器(3)在第i个位置对应的吸收光程Leff(Xi)的过程为:
步骤一、探测器(3)固定在该位置点,在气体池(1)内充入缓冲气氮气,调整激光光束入射至气体池(1),放大器(5)将探测器(3)采集的光强信号放大后输出,根据该放大后的信号计算获得多孔材料芯(2)散射出的氮气光强度I0;
步骤二、在气体池(1)内充入已知浓度的样品气体,调整激光光束入射至气体池(1),放大器(5)将探测器(3)采集的光强信号放大后输出,根据该放大后的信号计算获得多孔材料芯(2)散射出的样品光强度It;
步骤三、根据步骤一获取的氮气光强度I0和步骤二获取的样品光强度It按公式
计算获得探测器(3)位于该位置时对应的光程Leff(Xi),
式中,N为样品气体的浓度,σ为样品气体的吸收截面。
7.根据权利要求6所述用于激光吸收光谱技术的吸收光程延长的方法,其特征在于,入射至气体池(1)的激光光束与多孔材料芯(2)入射面的法线所成的角度θ范围为0~45度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201210058360 CN102590092B (zh) | 2012-03-07 | 2012-03-07 | 用于激光吸收光谱技术的吸收光程延长的装置及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201210058360 CN102590092B (zh) | 2012-03-07 | 2012-03-07 | 用于激光吸收光谱技术的吸收光程延长的装置及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102590092A true CN102590092A (zh) | 2012-07-18 |
CN102590092B CN102590092B (zh) | 2013-09-25 |
Family
ID=46479034
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201210058360 Expired - Fee Related CN102590092B (zh) | 2012-03-07 | 2012-03-07 | 用于激光吸收光谱技术的吸收光程延长的装置及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102590092B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105572797A (zh) * | 2016-02-15 | 2016-05-11 | 欧阳征标 | 一种太赫兹波脉冲调幅信号与光脉冲调幅信号变换放大器 |
EP3401666A1 (de) | 2017-05-11 | 2018-11-14 | Mettler-Toledo GmbH | Gasmesssystem |
CN109870414A (zh) * | 2019-04-08 | 2019-06-11 | 大连理工大学 | 一种散射增强型气体传感探头 |
CN110940632A (zh) * | 2019-10-31 | 2020-03-31 | 河南农业大学 | 一种基于tdlas的甲烷气体浓度的检测装置及检测方法 |
US11243161B1 (en) | 2020-11-20 | 2022-02-08 | Industrial Technology Research Institute | Gas measurement device and gas measurement method |
CN116519622A (zh) * | 2023-02-03 | 2023-08-01 | 湖北工业大学 | 基于光程可调光谱检测的复杂混合气体检测装置及方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5222389A (en) * | 1990-04-02 | 1993-06-29 | Gaztech International Corporation | Multi-channel gas sample chamber |
JPH06300685A (ja) * | 1993-04-16 | 1994-10-28 | Liquid Gas:Kk | ガス検知器 |
JP2002139425A (ja) * | 2000-11-02 | 2002-05-17 | Anritsu Corp | 校正用光検出器 |
CN1808100A (zh) * | 2005-12-28 | 2006-07-26 | 华东师范大学 | 便携式红外半导体激光吸收式瓦斯气体检测方法及其装置 |
WO2011042439A1 (de) * | 2009-10-06 | 2011-04-14 | Hochschule Regensburg | Miniaturisierte online-spurenanalytik |
CN102216755A (zh) * | 2008-11-14 | 2011-10-12 | 株式会社Ihi | 气体成分浓度测定装置 |
-
2012
- 2012-03-07 CN CN 201210058360 patent/CN102590092B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5222389A (en) * | 1990-04-02 | 1993-06-29 | Gaztech International Corporation | Multi-channel gas sample chamber |
JPH06300685A (ja) * | 1993-04-16 | 1994-10-28 | Liquid Gas:Kk | ガス検知器 |
JP2002139425A (ja) * | 2000-11-02 | 2002-05-17 | Anritsu Corp | 校正用光検出器 |
CN1808100A (zh) * | 2005-12-28 | 2006-07-26 | 华东师范大学 | 便携式红外半导体激光吸收式瓦斯气体检测方法及其装置 |
CN102216755A (zh) * | 2008-11-14 | 2011-10-12 | 株式会社Ihi | 气体成分浓度测定装置 |
WO2011042439A1 (de) * | 2009-10-06 | 2011-04-14 | Hochschule Regensburg | Miniaturisierte online-spurenanalytik |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105572797A (zh) * | 2016-02-15 | 2016-05-11 | 欧阳征标 | 一种太赫兹波脉冲调幅信号与光脉冲调幅信号变换放大器 |
EP3401666A1 (de) | 2017-05-11 | 2018-11-14 | Mettler-Toledo GmbH | Gasmesssystem |
CN110621980A (zh) * | 2017-05-11 | 2019-12-27 | 梅特勒-托莱多有限公司 | 气体测量系统 |
US11327008B2 (en) | 2017-05-11 | 2022-05-10 | Mettler-Toledo Gmbh | Gas measurement system |
CN110621980B (zh) * | 2017-05-11 | 2023-12-15 | 梅特勒-托莱多有限公司 | 气体测量系统 |
CN109870414A (zh) * | 2019-04-08 | 2019-06-11 | 大连理工大学 | 一种散射增强型气体传感探头 |
CN110940632A (zh) * | 2019-10-31 | 2020-03-31 | 河南农业大学 | 一种基于tdlas的甲烷气体浓度的检测装置及检测方法 |
CN110940632B (zh) * | 2019-10-31 | 2022-04-26 | 河南农业大学 | 一种基于tdlas的甲烷气体浓度的检测装置及检测方法 |
US11243161B1 (en) | 2020-11-20 | 2022-02-08 | Industrial Technology Research Institute | Gas measurement device and gas measurement method |
CN116519622A (zh) * | 2023-02-03 | 2023-08-01 | 湖北工业大学 | 基于光程可调光谱检测的复杂混合气体检测装置及方法 |
CN116519622B (zh) * | 2023-02-03 | 2023-10-10 | 湖北工业大学 | 基于光程可调光谱检测的复杂混合气体检测装置及方法 |
US11933724B1 (en) | 2023-02-03 | 2024-03-19 | Hubei University Of Technology | Device of complex gas mixture detection based on optical-path-adjustable spectrum detection and method therefor |
Also Published As
Publication number | Publication date |
---|---|
CN102590092B (zh) | 2013-09-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102590092B (zh) | 用于激光吸收光谱技术的吸收光程延长的装置及方法 | |
CN103499521B (zh) | 纳米颗粒关键几何特征量的测量方法 | |
CN104596645B (zh) | 复杂环境全偏振多角度散射模拟及测试系统 | |
CN103048653B (zh) | 一种微脉冲激光雷达的系统常数定标方法 | |
CN103454203A (zh) | 一种大气颗粒物粒径和化学成分的实时在线测量系统及测量方法 | |
US9804087B2 (en) | Hemispherical scanning optical scatterometer | |
CN103499391A (zh) | 光谱测量系统 | |
CN102565008B (zh) | 一种利用积分球测量材料透射比的方法及装置 | |
CN103363926A (zh) | 测速光幕靶光幕参数的检测装置及方法 | |
CN110243729A (zh) | 粒子计数器 | |
CN102486402A (zh) | 一种测量脉冲激光能量的方法及系统 | |
CN204594848U (zh) | 一种大气颗粒物浓度的监测装置 | |
CN103499393A (zh) | 光谱的测量方法 | |
EP2843394A1 (en) | Turbidity measuring sensor and method | |
CN104237085A (zh) | 动态光散射多角度可调光纤探头检测装置与方法 | |
CN102384783B (zh) | 一种高能激光半积分球阵列衰减器 | |
CN103528528A (zh) | 紧凑型精密激光三角测距仪 | |
CN205449753U (zh) | 一种基于片光源的二维烟雾浓度场测量装置 | |
CN105092426A (zh) | 纳米颗粒90度散射光谱的测量方法 | |
CN104849724B (zh) | 一种气溶胶激光雷达比的测量方法及装置 | |
CN202133468U (zh) | 一种用于测量脉冲激光能量的系统 | |
CN106769731A (zh) | 颗粒物浓度的测量方法及装置 | |
CN116242790A (zh) | 基于非分光红外原理的长短双光路测量系统及方法 | |
US20070146841A1 (en) | Instrument for measuring particle parameters | |
CN201903495U (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130925 |