CN110702236A - 一种用于光纤波长标准器的高精度光谱定标方法 - Google Patents
一种用于光纤波长标准器的高精度光谱定标方法 Download PDFInfo
- Publication number
- CN110702236A CN110702236A CN201911084805.8A CN201911084805A CN110702236A CN 110702236 A CN110702236 A CN 110702236A CN 201911084805 A CN201911084805 A CN 201911084805A CN 110702236 A CN110702236 A CN 110702236A
- Authority
- CN
- China
- Prior art keywords
- wavelength
- light beam
- optical fiber
- optical
- optical power
- 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
- 239000013307 optical fiber Substances 0.000 title claims abstract description 68
- 238000000034 method Methods 0.000 title claims abstract description 30
- 238000001228 spectrum Methods 0.000 title abstract description 7
- 238000010521 absorption reaction Methods 0.000 claims abstract description 67
- 230000003287 optical effect Effects 0.000 claims abstract description 60
- 230000003595 spectral effect Effects 0.000 claims abstract description 16
- BJQHLKABXJIVAM-UHFFFAOYSA-N bis(2-ethylhexyl) phthalate Chemical compound CCCCC(CC)COC(=O)C1=CC=CC=C1C(=O)OCC(CC)CCCC BJQHLKABXJIVAM-UHFFFAOYSA-N 0.000 claims description 9
- 238000012544 monitoring process Methods 0.000 claims description 7
- 239000000835 fiber Substances 0.000 claims description 5
- 239000000203 mixture Substances 0.000 claims 1
- 239000007789 gas Substances 0.000 description 20
- 238000005259 measurement Methods 0.000 description 12
- 238000000862 absorption spectrum Methods 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 238000001514 detection method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000005622 photoelectricity Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010183 spectrum analysis Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J9/00—Measuring optical phase difference; Determining degree of coherence; Measuring optical wavelength
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M11/00—Testing of optical apparatus; Testing structures by optical methods not otherwise provided for
- G01M11/02—Testing optical properties
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
本发明提出了一种用于光纤波长标准器的高精度光谱定标方法，包括以下步骤：S1、利用可调谐激光器发射激光；S2、利用光纤分束器对激光进行分束处理，得到第一光束和第二光束；S3、将第一光束通过待定标的光纤波长标准器，利用光功率计分别测量第一光束和第二光束的光功率，利用光波长计测量第二光束的波长；S4、基于第一光束和第二光束的光功率、第二光束的波长获得吸收峰扫描曲线，计算吸收峰波长；S5、将计算得到的吸收峰波长赋予待定标的光纤波长标准器，实现光纤波长标准器的光谱定标。本发明方法可以为光纤波长标准器提供高精度定标，为高分辨率、高精度光纤波长提供了量值溯源的可靠途径与保障。
Description
技术领域
本发明涉及一种高精度的光纤波长标准器定标方法，属于光纤波长测量技术领域。
背景技术
光纤波长参数是决定密集波分复用系统性能、光纤传感探测与解调精度、相干探测与扫描精度、各种光源波长与光谱分析等先进技术的重要参数，广泛应用于大容量光纤传输系统、光纤传感系统、信息通信与光电检测、医学生物等诸多领域。光纤波长标准器是光纤波长测量领域常用的标准装置，通过光纤波长标准器可以准确测量光纤的波长，得到高分辨率、高精度的光纤波长参数。
目前市面上通用的光纤波长标准器为光纤连接的气体吸收室，气体吸收室内会填充不同的气体，比如HCN气体吸收池，其吸收光谱范围覆盖1525nm～1565nm波段。在使用过程中，气体吸收室受到多种因素影响，其测量结果可能出现偏差，此时必须为气体吸收室赋予准确的波长值，实现气体吸收室的光谱定标，修正气体吸收室的偏差。
发明内容
针对现有技术中光纤波长标准器需要进行校准的问题，本发明提出了一种用于光纤波长标准器的高精度光谱定标方法，将波长不同的激光通过待定标的光纤波长标准器，实时监测输出光的功率和波长，绘制吸收峰扫描曲线，用吸收峰波长作为光纤波长标准器的标准波长，实现光纤波长标准器的高精度定标。
为解决上述技术问题，本发明采用了如下技术手段：
一种用于光纤波长标准器的高精度光谱定标方法，包括以下步骤：
S1、利用可调谐激光器发射不同波长的激光；
S2、利用光纤分束器对激光进行分束处理，得到第一光束和第二光束；
S3、将第一光束通过待定标的光纤波长标准器，并利用光功率计分别测量第一光束和第二光束的光功率，利用光波长计测量第二光束的波长；
S4、基于第一光束和第二光束的光功率、第二光束的波长获得吸收峰扫描曲线，计算吸收峰波长；
S5、将计算得到的吸收峰波长赋予待定标的光纤波长标准器，实现光纤波长标准器的光谱定标。
进一步的，所述步骤S2中的光纤分束器采用90/10分束器，所述第一光束具有90％的激光，所述第二光束具有10％的激光。
进一步的，所述步骤S3的具体操作如下：
S31、将所述第一光束通过待定标的光纤波长标准器，利用光功率计实时测量输出光的光功率；
S32、利用光纤分束器对所述第二光束进行分束处理得到第三光束和第四光束；
S33、利用光功率计实时监测第三光束的光功率，并利用光波长计实时测量第四光束的波长。
进一步的，所述步骤S32中的光纤分束器采用50/50分束器。
进一步的，步骤S4的具体操作如下：
S41、计算第一光束的光功率与第三光束的光功率的比值，进行光功率修正；
S42、利用光功率的比值与第四光束的波长绘制吸收峰扫描曲线；
S43、基于线性拟合算法获得吸收峰扫描曲线的拟合结果，计算吸收峰扫描曲线的中心波长，即吸收峰波长。
进一步的，所述吸收峰扫描曲线包括洛伦兹线型、高斯线型和高斯洛伦兹混合线型。
采用以上技术手段后可以获得以下优势：
本发明提出了一种用于光纤波长标准器的高精度光谱定标方法，根据待定标的光纤波长标准器的吸收光谱范围，利用可调谐激光器发射不同波长的激光，进行吸收峰扫描，将90％的激光通过光纤波长标准器，可以保证有足够多的光功率进入标准器被吸收，从而增加吸收峰的测量准确性。正常情况下，第一光束的光功率存在漂移影响，无法直接用于绘制吸收谱线，本发明利用第三光束的光功率去除光功率漂移的影响，修正吸收谱线，并根据吸收峰扫描曲线的线型进行数据拟合，减小数据拟合误差，可以计算出高精度的吸收峰波长，确保定标的准确性。本发明方法可以通过自动控制程序实现光源波长调谐、光功率与波长监测、吸收峰数据处理的流程，大大减小吸收峰波长测量时间，使得背景温度变化、光源与波长计功率/波长的变化对结果的影响降低到可以忽略的水平。本发明方法为高分辨率、高精度光纤波长提供了量值溯源的可靠途径与保障，可以应用于国家计量部门、科研院所、校准检测企业、光通信设备及器件生产企业，促进我国新一代光纤通信、信息与传感、光电、医学生物等相关技术领域的发展。
附图说明
图1为本发明一种用于光纤波长标准器的高精度光谱定标方法的步骤流程图。
图2为本发明一种用于光纤波长标准器的高精度光谱定标方法的原理图。
图3为光纤波长标准器的结构示意图。
图4为本发明实施例中吸收峰扫描曲线的拟合结果示意图。
图中，1是光纤，2是气体吸收管。
具体实施方式
下面结合附图对本发明的技术方案作进一步说明：
一种用于光纤波长标准器的高精度光谱定标方法，如图1、2所示，具体包括以下步骤：
S1、利用可调谐激光器发射不同波长的激光。本发明方法是对气体吸收管类型的光纤波长标准器的吸收峰波长进行高精度测量的方法，如图3所示，气体吸收管类型的光纤波长标准器通常是两根光纤1连接气体吸收管2的两端，待测光从一根光纤进入气体吸收管，待测光被吸收管内的气体吸收衰减后输出，根据吸收情况计算待测光的波长。这类标准器内的气体不同，其吸收光谱范围(即吸收峰的波长范围)也会不同，比如HCN气体，其吸收光谱范围覆盖1525nm～1565nm波段。根据待定标的光纤波长标准器的吸收光谱范围，利用程序控制可调谐激光器发射吸收光谱范围内的激光。
S2、利用光纤分束器对激光进行分束处理，得到第一光束和第二光束，其中，光纤分束器采用90/10分束器，第一光束具有90％的激光，第二光束具有10％的激光。由于吸收峰的谱宽在pm量级，同时要求扫描步长尽量小，根据不同分子、不同气压，步长一般在0.1pm～10pm，因此在吸收峰的扫描过程中需要保证光源功率稳定性与波长稳定性足够好，具有90％的第一光束可以保证有足够多的光功率进入标准器被吸收，确保后续测量和计算的准确性。
S3、将第一光束通过待定标的光纤波长标准器，并利用光功率计分别测量第一光束和第二光束的光功率，利用光波长计测量第二光束的波长；具体操作如下：
S31、将第一光束通过待定标的光纤波长标准器，利用光功率计实时测量光纤波长标准器输出光的光功率。
S32、利用光纤分束器对第二光束进行分束处理，得到第三光束和第四光束，此处的光纤分束器采用50/50分束器，第三光束和第四光束分别获得50％的第二光束的激光。
S33、利用光功率计实时监测第三光束的光功率，并利用光波长计实时测量第四光束的波长。由于第三光束与第四光束是1:1关系，两者的光功率应该相同，对第三光束和第四光束的实时测量可以保证在每个吸收峰的扫描过程中，实时监测光波长与光功率的变化情况，方便后续的计算。
S4、基于第一光束和第二光束的光功率、第二光束的波长获得吸收峰扫描曲线，计算吸收峰波长；具体操作如下：
S41、用步骤S31中测量到的第一光束的光功率除以步骤S33中测量到的第三光束的光功率，计算第一光束的光功率与第三光束的光功率的比值，消除第一光束中光功率漂移的影响，进行光功率修正。
S42、利用光功率的比值与第四光束的波长绘制吸收峰扫描曲线，吸收峰扫描曲线的横坐标是波长，纵坐标是光功率比值。
S43、基于线性拟合算法获得吸收峰扫描曲线的拟合结果，计算吸收峰扫描曲线的中心波长，即吸收峰波长。由于不同气压的气体吸收峰的光谱展宽线型不同，因此在绘制完吸收峰扫描曲线后，需要对扫描曲线进行不同线型的拟合结果的比较，确定最好的拟合结果，从而减小了数据拟合误差。根据研究发现，吸收峰扫描曲线包括洛伦兹线型、高斯线型和高斯洛伦兹混合线型。在本发明实施例中，通过对不同吸收线型的数据拟合算法的比较，确定了洛伦兹线型的拟合结果，如图4所示，根据拟合结果计算中心波长，即图中最低点的波长值，xc＝1554.58965，该中心波长就是本发明实施例的吸收峰波长。
S5、将计算得到的吸收峰波长赋予待定标的光纤波长标准器，实现光纤波长标准器的光谱定标。
本发明方法通过二级分束监测操作，保证了光源功率的稳定性与波长的稳定性，第一次分束确保光纤波长标准器吸收峰测量的准确性，第二次分束用于对监测结果进行修正和处理，获得稳定重复的吸收峰扫描曲线，直接保障了吸收峰波长测量结果的准确。确保了测量结果的准确性之后，本发明还通过不同吸收线型的拟合算法进行拟合处理，计算吸收峰波长，通过对不同拟合结果的对比，减小数据拟合误差，确保计算出高精度的吸收峰波长。本发明方法可以通过自动控制程序实现光源波长调谐、光功率与波长监测、吸收峰数据处理的流程，大大减小吸收峰波长测量时间，使得背景温度变化、光源与波长计功率/波长的变化对结果的影响降低到可以忽略的水平。本发明方法为高分辨率、高精度光纤波长提供了量值溯源的可靠途径与保障，可以应用于国家计量部门、科研院所、校准检测企业、光通信设备及器件生产企业，促进我国新一代光纤通信、信息与传感、光电、医学生物等相关技术领域的发展。
上面结合附图对本发明的实施方式作了详细地说明，但是本发明并不局限于上述实施方式，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下做出各种变化。
Claims (6)
1.一种用于光纤波长标准器的高精度光谱定标方法，其特征在于，包括以下步骤：
S1、利用可调谐激光器发射不同波长的激光；
S2、利用光纤分束器对激光进行分束处理，得到第一光束和第二光束；
S3、将第一光束通过待定标的光纤波长标准器，并利用光功率计分别测量第一光束和第二光束的光功率，利用光波长计测量第二光束的波长；
S4、基于第一光束和第二光束的光功率、第二光束的波长获得吸收峰扫描曲线，计算吸收峰波长；
S5、将计算得到的吸收峰波长赋予待定标的光纤波长标准器，实现光纤波长标准器的光谱定标。
2.根据权利要求1所述的一种用于光纤波长标准器的高精度光谱定标方法，其特征在于，所述步骤S2中的光纤分束器采用90/10分束器，所述第一光束具有90％的激光，所述第二光束具有10％的激光。
3.根据权利要求1所述的一种用于光纤波长标准器的高精度光谱定标方法，其特征在于，所述步骤S3的具体操作如下：
S31、将所述第一光束通过待定标的光纤波长标准器，利用光功率计实时测量输出光的光功率；
S32、利用光纤分束器对所述第二光束进行分束处理得到第三光束和第四光束；
S33、利用光功率计实时监测第三光束的光功率，并利用光波长计实时测量第四光束的波长。
4.根据权利要求3所述的一种用于光纤波长标准器的高精度光谱定标方法，其特征在于，所述步骤S32中的光纤分束器采用50/50分束器。
5.根据权利要求4所述的一种用于光纤波长标准器的高精度光谱定标方法，其特征在于，步骤S4的具体操作如下：
S41、计算第一光束的光功率与第三光束的光功率的比值，进行光功率修正；
S42、利用光功率的比值与第四光束的波长绘制吸收峰扫描曲线；
S43、基于线性拟合算法获得吸收峰扫描曲线的拟合结果，计算吸收峰扫描曲线的中心波长，即吸收峰波长。
6.根据权利要求5所述的一种用于光纤波长标准器的高精度光谱定标方法，其特征在于，所述吸收峰扫描曲线包括洛伦兹线型、高斯线型和高斯洛伦兹混合线型。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911084805.8A CN110702236B (zh) | 2019-11-08 | 2019-11-08 | 一种用于光纤波长标准器的高精度光谱定标方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911084805.8A CN110702236B (zh) | 2019-11-08 | 2019-11-08 | 一种用于光纤波长标准器的高精度光谱定标方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110702236A true CN110702236A (zh) | 2020-01-17 |
| CN110702236B CN110702236B (zh) | 2021-01-26 |
Family
ID=69204666
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911084805.8A Expired - Fee Related CN110702236B (zh) | 2019-11-08 | 2019-11-08 | 一种用于光纤波长标准器的高精度光谱定标方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110702236B (zh) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1146011A (zh) * | 1995-05-11 | 1997-03-26 | 仓敷纺绩株式会社 | 光谱测定法、光学测量法和装置 |
| JP2000304695A (ja) * | 1999-04-23 | 2000-11-02 | Mitsubishi Heavy Ind Ltd | 排ガス中の三酸化硫黄濃度測定装置 |
| US20020121592A1 (en) * | 2001-03-05 | 2002-09-05 | Blazo Stephen F. | Extended range frequency calibration device |
| CN101782513A (zh) * | 2009-01-20 | 2010-07-21 | 中国计量科学研究院 | 直接测量气体吸收谱形的平衡检测电路及检测方法 |
| CN102629731A (zh) * | 2012-02-14 | 2012-08-08 | 浙江嘉莱光子技术有限公司 | 一种激光器波长和功率同时稳定的控制方法及其控制装置 |
| CN102735644A (zh) * | 2012-07-06 | 2012-10-17 | 北京大方科技有限责任公司 | 一种原位式激光气体分析仪的在线标定方法 |
| CN102967580A (zh) * | 2012-11-09 | 2013-03-13 | 山东微感光电子有限公司 | 一种基于vcsel的低功耗气体检测方法及装置 |
| CN104535528A (zh) * | 2014-11-26 | 2015-04-22 | 东南大学 | Bp神经网络实时提取tdlas气体吸收光谱吸光度的方法 |
| CN105021310A (zh) * | 2015-07-28 | 2015-11-04 | 上海交通大学 | 高精度光纤光栅温度传感系统 |
| CN105890779A (zh) * | 2016-06-07 | 2016-08-24 | 哈尔滨工业大学 | 一种用于波长扫描系统的实时波长标定装置及标定方法 |
| CN107356915A (zh) * | 2017-09-11 | 2017-11-17 | 南京信息工程大学 | 一种中红外差分吸收激光雷达的标定方法和标定系统 |
-
2019
- 2019-11-08 CN CN201911084805.8A patent/CN110702236B/zh not_active Expired - Fee Related
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1146011A (zh) * | 1995-05-11 | 1997-03-26 | 仓敷纺绩株式会社 | 光谱测定法、光学测量法和装置 |
| JP2000304695A (ja) * | 1999-04-23 | 2000-11-02 | Mitsubishi Heavy Ind Ltd | 排ガス中の三酸化硫黄濃度測定装置 |
| US20020121592A1 (en) * | 2001-03-05 | 2002-09-05 | Blazo Stephen F. | Extended range frequency calibration device |
| CN101782513A (zh) * | 2009-01-20 | 2010-07-21 | 中国计量科学研究院 | 直接测量气体吸收谱形的平衡检测电路及检测方法 |
| CN102629731A (zh) * | 2012-02-14 | 2012-08-08 | 浙江嘉莱光子技术有限公司 | 一种激光器波长和功率同时稳定的控制方法及其控制装置 |
| CN102735644A (zh) * | 2012-07-06 | 2012-10-17 | 北京大方科技有限责任公司 | 一种原位式激光气体分析仪的在线标定方法 |
| CN102967580A (zh) * | 2012-11-09 | 2013-03-13 | 山东微感光电子有限公司 | 一种基于vcsel的低功耗气体检测方法及装置 |
| CN104535528A (zh) * | 2014-11-26 | 2015-04-22 | 东南大学 | Bp神经网络实时提取tdlas气体吸收光谱吸光度的方法 |
| CN105021310A (zh) * | 2015-07-28 | 2015-11-04 | 上海交通大学 | 高精度光纤光栅温度传感系统 |
| CN105890779A (zh) * | 2016-06-07 | 2016-08-24 | 哈尔滨工业大学 | 一种用于波长扫描系统的实时波长标定装置及标定方法 |
| CN107356915A (zh) * | 2017-09-11 | 2017-11-17 | 南京信息工程大学 | 一种中红外差分吸收激光雷达的标定方法和标定系统 |
Non-Patent Citations (2)
| Title |
|---|
| R. CHAUX, B. LAVOREL: "elative line intensity measurement in absorption spectra using a tunable diode laser at 1.6µm: application to the determination of 13CO2/12CO2 isotope ratio", 《APPLIED PHYSICS B》 * |
| 沈超 等: "光学气体吸收池住吸收光谱技术中的发展与应用", 《红外》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110702236B (zh) | 2021-01-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5983779B2 (ja) | ガス吸収分光装置及びガス吸収分光方法 | |
| CN101923162B (zh) | 拉曼激光雷达标定装置及其标定方法 | |
| US11209356B2 (en) | Gas absorption spectrometer and gas absorption spectroscopy method | |
| CN110068548B (zh) | 用于离轴积分腔系统中激光器的波长锁定装置及其锁定方法 | |
| CN101290289A (zh) | 紫外差分烟气浓度测量系统校准方法及实施装置 | |
| CN106500911B (zh) | 一种基于气体吸收谱线压力展宽效应的压力计校准方法 | |
| CN109342022B (zh) | 一种可调谐激光器波长动态标定装置和方法 | |
| CN102252982A (zh) | 一种用于激光气体分析仪的波长漂移补偿方法 | |
| WO2021036167A1 (zh) | 利用回音壁模式激光光源测fp透过率曲线的装置和方法 | |
| CN101813519B (zh) | 光谱仪杂散光校正方法 | |
| CN109579820B (zh) | 一种提高光纤陀螺仪标度因数性能的方法 | |
| CN112903628B (zh) | 一种负压状态下痕量气体检测方法 | |
| CN110702236B (zh) | 一种用于光纤波长标准器的高精度光谱定标方法 | |
| CN103996971A (zh) | Co2探测激光雷达中脉冲差频激光器稳频装置及方法 | |
| CN101915748B (zh) | 一种激光气体浓度检测的拟合动态寻峰方法 | |
| CN108548787B (zh) | 一种光谱仪精确测量硫化氢气体浓度的方法 | |
| CN108760653B (zh) | 一种光谱仪精确测量二氧化硫气体浓度的方法 | |
| CN115882325A (zh) | 一种中心波长锁定方法、装置及钠导星激光器 | |
| CN109115720B (zh) | 一种多程气体吸收池有效程长与气体浓度同时测量方法 | |
| CN114878515A (zh) | 一种可调谐二极管激光吸收信号时间-频率转换方法 | |
| CN102538967A (zh) | 一种高精度光谱仪多波长校准系统 | |
| CN106533551B (zh) | 一种基于光谱分辨率校正提高光信噪比测量精度的方法 | |
| CN107091730B (zh) | 一种估算光电倍增管的绝对光响应率的装置 | |
| CN217687505U (zh) | 一种激光光束的稳定分光装置及激光功率计标定的装置 | |
| CN110926621A (zh) | 针对傅立叶变换型太赫兹源波长测量仪的校准装置及方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | 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 |
Granted publication date: 20210126 Termination date: 20211108 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |