CN107941754B - 一种气体折射率的测量方法 - Google Patents

一种气体折射率的测量方法 Download PDF

Info

Publication number
CN107941754B
CN107941754B CN201711121495.3A CN201711121495A CN107941754B CN 107941754 B CN107941754 B CN 107941754B CN 201711121495 A CN201711121495 A CN 201711121495A CN 107941754 B CN107941754 B CN 107941754B
Authority
CN
China
Prior art keywords
gas
optical fiber
refractive index
optical
expressed
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.)
Active
Application number
CN201711121495.3A
Other languages
English (en)
Other versions
CN107941754A (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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN201711121495.3A priority Critical patent/CN107941754B/zh
Publication of CN107941754A publication Critical patent/CN107941754A/zh
Application granted granted Critical
Publication of CN107941754B publication Critical patent/CN107941754B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/41Refractivity; Phase-affecting properties, e.g. optical path length
    • G01N21/45Refractivity; Phase-affecting properties, e.g. optical path length using interferometric methods; using Schlieren methods

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开了一种基于气体折射率测量装置的测量方法,该装置包括宽谱光源、电光调制器和计算器,宽谱光源输出端连接第一光纤耦合器,第一光纤耦合器的一个输出端与气体测量池连接,气体测量池的输出光和第一光纤耦合器的输出光经过第二光纤耦合器合束,第一光纤耦合器、气体测量池和第二光纤耦合器构成一个马赫曾德干涉仪,马赫曾德干涉仪的输出端连接电光调制器,电光调制器输出的调制信号经过色散光纤后入射到高速光电探测器上,高速光电探测器将光信号装换成微波信号并通过低噪放放大,低噪放输出端连接微波功分器;本发明可实现高精度折射率的测量,同时也适用于测量透明液体和固体的折射率。

Description

一种气体折射率的测量方法
技术领域
本发明涉及一种测量方法,具体为一种气体折射率的测量方法。
背景技术
折射率是光在真空(因为在空气中与在真空中的传播速度差不多,所以一般用在空气的传播速度)中的速度与光在该材料中的速度之比率。材料的折射率越高,使入射光发生折射的能力也就越强。折射率越高,镜片越薄,即镜片中心厚度相同,相同度数同种材料,折射率高的比折射率低的镜片边缘更薄。折射率与介质的电磁性质密切相关,折射率是表征物质光学物质特性的基本物理量,该参数是决定物质合成、制造和在各领域中应用的重要条件。
目前气体折射率的测量方法有光干涉法、表面等离子提共振法、光纤传感方法和临界角法,这些方法往往存在测量分辨率不高的缺陷,不能满足精确测量的要求,因此本发明提出一种新的基于微波光子技术的微波气体折射率测量方法,已解决上述问题。
发明内容
本发明要解决的技术问题是克服现有气体折射率测量装置测量分辨率不高的缺陷,提供一种气体折射率的测量装置,从而解决上述问题。
为了解决上述技术问题,本发明提供了如下的技术方案:
本发明为一种气体折射率的测量装置,包括宽谱光源、电光调制器和计算器,宽谱光源输出端连接第一光纤耦合器,第一光纤耦合器的一个输出端与气体测量池连接,气体测量池的输出光和第一光纤耦合器的输出光经过第二光纤耦合器合束,第一光纤耦合器、气体测量池和第二光纤耦合器构成一个马赫曾德干涉仪,马赫曾德干涉仪的输出端连接电光调制器,电光调制器输出的调制信号经过色散光纤后入射到高速光电探测器上,高速光电探测器将光信号装换成微波信号并通过低噪放放大,低噪放输出端连接微波功分器,微波功分器将一部分微波信号注入到电光调制器中,同时将另一部分微波信号输入频谱仪,频谱仪末端连接计算机。
作为本发明的一种优选技术方案,电光调制器、色散光纤、高速光电探测器、低噪放和微波功分器组成光电振荡器环路,且光电振荡器环路输入端与马赫曾德干涉仪的输出端相连接,能够将马赫曾德干涉仪输出端产生的正弦梳状谱注入光电振荡器环路中,并通过光电振荡器环路产生微波信号。
作为本发明的一种优选技术方案,气体测量池包括左自聚焦透镜、右自聚焦透镜和进气阀,两个自聚焦透镜之间的距离为d,且各器件都密封安装在气体测量池的铝外壳中,通过在气体测量池中注入待测气体导致马赫曾德干涉仪光程差的改变,从而改变光电振荡器输出的微波信号的中心频率,根据微波信号中心频率的变化量来得到待测气体的折射率。
作为本发明的一种优选技术方案,宽谱光源可采用高斯型或矩形光源作为发射光源,使得光源发射装置的选择性更高。
本发明所达到的有益效果是:本发明为一种气体折射率的测量装置,通过第一光纤耦合器,气体测量池和第二光纤耦合器构成马赫曾德干涉仪,使得宽谱光源经该干涉仪后,当干涉仪两臂光程差在光源相干范围内时,在干涉仪的输出端将产生干涉条纹,该干涉条纹在频域上为一正弦梳状谱;通过电光调制器,色散光纤,高速光电探测器,低噪放和微波功分器构成光电振荡器环路,将光信号装换成微波信号并测量输出微波信号的中心频率;通过设置的气体测量池,能够改变马赫曾德干涉仪光程差,从而改变光电振荡器输出的微波信号的中心频率,根据微波信号中心频率的变化量来得到待测气体的折射率;本发明提出的测试原理和方法可实现高精度折射率的测量,同时也适用于测量透明液体和固体的折射率。
附图说明
附图用来提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例一起用于解释本发明,并不构成对本发明的限制。
在附图中:
图1是本发明系统框架结构示意图;
图2是本发明气体测量池结构示意图。
图中标号:101:宽谱光源;102:光纤起偏器;103:光纤耦合器;104:气体测量池;105:电光调制器;106:光纤耦合器;107:色散光纤;108:高速光电探测器;109:低噪放;201:微波功分器;202:频谱仪;203:计算机。
具体实施方式
以下结合附图对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
在本发明的描述中,需要说明的是,术语“竖直”、“上”、“下”、“水平”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
在本发明的描述中,还需要说明的是,除非另有明确的规定和限定,术语“设置”、“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。
实施例:宽谱光源101(可采用高斯型或矩形宽谱光源)经光纤起偏器102后进入一光纤耦合器103。该光纤耦合器的一个输出端和气体测量池104连接。光纤耦合器103的另一路光进入一电光调制器105,该调制器的输出光和气体测量池104的输出光经过光纤耦合器106合束。光纤耦合器103,气体测量池104,电光调制器105和光纤耦合器106构成一个马赫曾德干涉仪。该干涉仪的两臂的两路光通过光纤耦合器合路后进入色散补偿光纤107,两路光经色散补偿光纤后将产生延时,经过延时的光信号通过高速光电探测器108实现光电转换并通过低噪放109放大,放大后的微波信号经过微波功分器201后分成两路,一路注入到电光调制器105中,从而将电光调制器105,光纤耦合器106,色散补偿光纤107,高速光电探测器108,低噪放109和功分器201组成一个光电振荡器环路,使得环路有微波信号产生,并且该微波信号的输出频率与马赫曾德干涉仪的两臂光程差有关。由光电振荡器产生的微波信号通过电光调制器105被调制到光域上,该光载微波信号经过色散补偿光纤107后入射到高速光电探测器108上,该探测器将光信号装换成微波信号,该微波信号通过低噪放109放大后经过一个微波功分器201功分后,一部分微波信号注入到电光调制器105中,一部分信号通过频谱仪202来测量光电振荡器输出微波信号的中心频率并通过计算机203来记录该微波信号中心频率的变化。
气体测量池结构如图2所示:
气体测量池通过进气阀303将气体注入进测量池中,301和302分别为通光的自聚焦透镜,形成光路的对射,两个自聚焦透镜之间的距离为d,所有器件都密封安装在气体测量池的铝外壳304中。
具体测量原理如下:
本测量方法的原理是通过在气体测量池中注入待测气体导致马赫曾德干涉仪光程差的改变,从而改变光电振荡器输出的微波信号的中心频率,根据微波信号中心频率的变化量来得到待测气体的折射率。宽谱光源经过马赫曾德干涉仪后,将发生干涉,干涉条纹的输出在频域上可表示为:
宽谱光源的电场可表示为:
Figure BDA0001467435280000041
式中ω为光源频率。则光源的光功率谱密度可表示为:
T(ω)=|E(ω)|2 (2)
光源经干涉以后,光路1中的气体测量池在注入待测气体前后,由于光路中气体折射率的变化,导致各光谱成份产生了一定的延时,在频域上可表示为:
E1(ω)=A1E(ω)ejωΔτ (3)
上式中A1为光路1的幅度衰减系数,Δτ为待测气体相对与测量池中无气体时折射率的变化而引入的延迟量,可表示为:
Figure BDA0001467435280000042
式中c为光速,n为待测气体的折射率,d为气体测量池中两个自聚焦透镜之间的距离。干涉仪的光路2被射频信号调制,被调制后的光在频域上可表示为:
Figure BDA0001467435280000043
上式中A2为光路2的幅度衰减系数,ξ为射频信号的角频率,
Figure BDA0001467435280000044
为调制的载波和边带之间的相位差。两路光经过干涉仪的第二个耦合器后,输出为:
Figure BDA0001467435280000045
干涉仪的输出经过一段色散光纤延时后,光载波发生了一个时延,该时延线的电场传递函数可表示为:
H(ω)=|H(ω)|e-jφ(ω) (7)
φ(ω)为色散光纤延迟引入的相位,根据泰勒级数展开,该相位可表示为:
Figure BDA0001467435280000046
式中,τ(ω0)为中心频率为ω0时的群时延,β为光纤的色散,其单位为ps2/km,β可表示为:
Figure BDA0001467435280000047
式中D(ps/km/nm)为光纤的色散系数,λ0为光源波长。
光电振荡器输出的响应函数可表示为:
HRF(ξ)=∫T(ω)[H*(ω)H(ω+ξ)+H(ω)H*(ω-ξ)]dω (10)
由式(6)——(9)可得响应函数为:
Figure BDA0001467435280000048
上式中H(ω)为理想的光电振荡器的响应函数,可表示为:H(ξ)=∫T(ω)exp[-jξβL(ω-ω0)]dω(11)
由此可知,光电振荡器输出的微波信号的中心频率可表示为:
Figure BDA0001467435280000051
根据式(4)可得,待测气体的折射率为:
Figure BDA0001467435280000052
由上式可知,根据光电振荡器输出的射频信号的频率,光源的中心波长,气体测量池中两自聚焦透镜之间的距离,色散光纤的色散值和长度就可以得到待测气体的折射率。本发明提出的测量系统需要先确定光电振荡器中的色散光纤的长度和色散值,光源的中心频率,根据上述参数计算两自聚焦透镜之间的距离使得光电振荡器输出的微波信号的频率在常用频谱仪测量频段范围内(不失为一般性,一般频谱仪的频率带宽为几十KHz——26.5GHz)。本系统的分辨率可通过设置式(13)中的各项参数来改变。由式(12)可知光电振荡器输出的微波信号的中心频率3dB带宽可达80MHz左右,不失为一般性通过调节干涉仪两臂光程差使得系统输出的微波信号的频率分辨率为100MHz,色散光纤为1km,色散系数为-150ps/km/nm时,光源中心波长为1550nm,气体测量池中两自聚焦透镜之间的距离为10cm时,测试系统的折射率测量分辨率可达0.00024,从而可实现高精度折射率的测量。本发明提出的测试原理和方法也同样适用于测量透明液体和固体的折射率。
本发明提出的气体折射率测量系统的工作流程如下所示:
上电后,调制器驱动板通过程序自动控制强度型光调制器工作在线性工作点。调制器工作点确定后,先不将待测气体充入到气体测量池中,此时记录光电振荡器输出的微波信号的中心频率为f1。将待测气体充入到气体测量池中,再次记录光电振荡器输出的微波信号的中心频率为f2。根据式(13)就可以得到待测光学元器件的厚度为
Figure BDA0001467435280000053
最后应说明的是:以上仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (1)

1.一种气体折射率的测量方法,该测量装置包括宽谱光源、电光调制器和计算器,宽谱光源输出端连接第一光纤耦合器,第一光纤耦合器的一个输出端与气体测量池连接,第一光纤耦合器的另一路进入一电光调制器,所述电光调制器的输出端和所述气体测量池的输出端经过经过第二光纤耦合器合束,第一光纤耦合器、气体测量池、电光调制器和第二光纤耦合器构成一个马赫曾德干涉仪,马赫曾德干涉仪的输出端连接色散补偿光纤,两路光经色散补偿光纤后产生延时,经过延时的光信号入射到高速光电探测器上,高速光电探测器将光信号装换成微波信号并通过低噪放放大,低噪放输出端连接微波功分器,微波功分器将一部分微波信号注入到电光调制器中,同时将另一部分微波信号输入频谱仪,频谱仪末端连接计算机;
其特征在于:通过电光调制器,第二光纤耦合器,色散补偿光纤,高速光电探测器,低噪放和微波功分器构成光电振荡器环路,将光信号装换成微波信号并测量输出微波信号的中心频率;通过设置的气体测量池,能够改变马赫曾德干涉仪光程差,从而改变光电振荡器输出的微波信号的中心频率,根据微波信号中心频率的变化量来得到待测气体的折射率;
通过在气体测量池中注入待测气体导致马赫曾德干涉仪光程差的改变,从而改变光电振荡器输出的微波信号的中心频率,根据微波信号中心频率的变化量来得到待测气体的折射率,气体折射率测量系统的工作流程如下:上电后,调制器驱动板通过程序自动控制强度光调制器工作在线性工作点,调制器工作点确定后,先不将待检测气体充入到气体检测池中,此时记录光电振荡器输出的微波信号的中心频率为f1,将待检测气体充入到气体检测池中,再次记录光电振荡器输出的微波信号的中心频率为f2;
宽谱光源经过马赫曾德干涉仪后,将发生干涉,干涉条纹的输出在频域上可表示为:
宽谱光源的电场可表示为:
Figure FDF0000010904080000011
式中ω为光源频率,则光源的光功率谱密度可表示为:
T(ω)=|E(ω)|2 (2)
光源经干涉以后,光路1中的气体测量池在注入待测气体前后,由于光路中气体折射率的变化,导致各光谱成份产生了一定的延时,在频域上可表示为:
E1(ω)=A1E(ω)ejωΔτ (3)
上式中A1为光路1的幅度衰减系数,Δτ为待测气体相对与测量池中无气体时折射率的变化而引入的延迟量,可表示为:
Figure FDF0000010904080000012
式中c为光速,n为待测气体的折射率,d为气体测量池中两个自聚焦透镜之间的距离,干涉仪的光路2被射频信号调制,被调制后的光在频域上可表示为:
Figure FDF0000010904080000013
上式中A2为光路2的幅度衰减系数,ξ为射频信号的角频率,
Figure FDF0000010904080000021
为调制的载波和边带之间的相位差,两路光经过干涉仪的第二个耦合器后,输出为:
Figure FDF0000010904080000022
干涉仪的输出经过一段色散光纤延时后,光载波发生了一个时延,该时延线的电场传递函数可表示为:
H(ω)=|H(ω)|e-jφ(ω) (7)
φ(ω)为色散光纤延迟引入的相位,根据泰勒级数展开,该相位可表示为:
Figure FDF0000010904080000023
式中,τ(ω0)为中心频率为ω0时的群时延,β为光纤的色散,其单位为ps2/km,β可表示为:
Figure FDF0000010904080000024
式中D(ps/km/nm)为光纤的色散系数,λ0为光源波长;
光电振荡器输出的响应函数可表示为:
HRF(ξ)=∫T(ω)[H*(ω)H(ω+ξ)+H(ω)H*(ω-ξ)]dω (10)
由式(6)——(9)可得响应函数为:
Figure FDF0000010904080000025
上式中H(ω)为理想的光电振荡器的响应函数,可表示为
H(ξ)=∫T(ω)exp[-jξβL(ω-ω0)]dω (11)
由此可知,光电振荡器输出的微波信号的中心频率可表示为:
Figure FDF0000010904080000026
根据式(4)可得,待测气体的折射率为:
Figure FDF0000010904080000027
由上式可知,根据光电振荡器输出的射频信号的频率,光源的中心波长,气体测量池中两自聚焦透镜之间的距离,色散光纤的色散值和长度就能够得到待测气体的折射率。
CN201711121495.3A 2017-11-14 2017-11-14 一种气体折射率的测量方法 Active CN107941754B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711121495.3A CN107941754B (zh) 2017-11-14 2017-11-14 一种气体折射率的测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711121495.3A CN107941754B (zh) 2017-11-14 2017-11-14 一种气体折射率的测量方法

Publications (2)

Publication Number Publication Date
CN107941754A CN107941754A (zh) 2018-04-20
CN107941754B true CN107941754B (zh) 2021-06-01

Family

ID=61934002

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711121495.3A Active CN107941754B (zh) 2017-11-14 2017-11-14 一种气体折射率的测量方法

Country Status (1)

Country Link
CN (1) CN107941754B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107941752B (zh) * 2017-11-14 2020-10-13 奥菲(泰州)光电传感技术有限公司 一种气体折射率的测量装置
CN108844454A (zh) * 2018-06-25 2018-11-20 中国科学院上海光学精密机械研究所 一种干涉仪相移装置
CN108955940B (zh) * 2018-07-19 2021-01-29 南京航空航天大学 一种光纤光栅温度传感解调方法
CN112857609B (zh) * 2020-12-30 2023-09-22 上海第二工业大学 一种面向光纤端面接触液体折射率测量的微波光子滤波器系统及方法
CN112834070B (zh) * 2020-12-30 2023-09-22 上海第二工业大学 一种利用微波光子滤波器测量光纤端面接触气体温度的方法
CN114354542B (zh) * 2021-10-10 2024-06-04 上海第二工业大学 一种利用马赫-曾德尔干涉仪结构的微波光子滤波器测量液体折射率的装置及方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101576488A (zh) * 2009-06-05 2009-11-11 西南石油大学 光电混合集成硫化氢气体浓度传感器装置及其测试方法
CN101936879A (zh) * 2010-07-28 2011-01-05 山东大学 一种基于马赫曾德干涉仪的光声光谱气体检测系统
CN103022857A (zh) * 2012-12-10 2013-04-03 贵州大学 高稳定单模微波光电振荡器
CN103278150A (zh) * 2013-05-10 2013-09-04 浙江大学 一种检测角速度的光载微波陀螺方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201837582U (zh) * 2010-10-19 2011-05-18 天津天狮生物发展有限公司 一体化光学臭氧产量检测装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101576488A (zh) * 2009-06-05 2009-11-11 西南石油大学 光电混合集成硫化氢气体浓度传感器装置及其测试方法
CN101936879A (zh) * 2010-07-28 2011-01-05 山东大学 一种基于马赫曾德干涉仪的光声光谱气体检测系统
CN103022857A (zh) * 2012-12-10 2013-04-03 贵州大学 高稳定单模微波光电振荡器
CN103278150A (zh) * 2013-05-10 2013-09-04 浙江大学 一种检测角速度的光载微波陀螺方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《下册》;上海交通大学物理教研室编;《大学物理教程》;20110731;第182页 *
《宽带可调谐光电振荡器研究》;刘辰钧;《中国优秀硕士学位论文全文数据库》;20140101;第26-29页 *

Also Published As

Publication number Publication date
CN107941754A (zh) 2018-04-20

Similar Documents

Publication Publication Date Title
CN107941754B (zh) 一种气体折射率的测量方法
CN107872274B (zh) 一种光纤色散系数的测量方法
CN102607720B (zh) 一种测量光程的方法和系统
CN103091072B (zh) 基于光单边带调制的光器件测量方法、测量装置
CN107894327B (zh) 一种光纤色散系数的测量装置
CN104316186A (zh) 一种基于光学频率梳的光谱测量方法
CN108827601A (zh) 一种光纤干涉仪臂长差的测量装置
CN105091776A (zh) 基于单边带扫频调制的光纤激光静态应变拍频解调系统
Zhang et al. A PGC-DCDM demodulation scheme insensitive to phase modulation depth and carrier phase delay in an EOM-based SPM interferometer
CN107917669B (zh) 一种光纤位移传感器解调方法
CN103414513B (zh) 一种具有高动态范围的脉冲光动态消光比测量装置及方法
CN104950311A (zh) 基于oeo的自校准的大量程、高精度绝对距离测量系统
CN107941752B (zh) 一种气体折射率的测量装置
CN108955939B (zh) 一种光纤光栅温度传感解调系统
CN108844717A (zh) 一种光纤干涉仪臂长差的测量方法
CN107976300B (zh) 一种保偏光纤拍长的测量方法
CN107631694B (zh) 一种光学元器件厚度的测量方法
CN107796422B (zh) 一种光纤位移传感器解调装置
CN107907058B (zh) 一种光学元器件厚度的测量装置
CN107806981B (zh) 一种保偏光纤拍长的测量装置
CN101592526A (zh) 一种光平均波长的测量方法及装置
CN105353210A (zh) 一种高灵敏大带宽光子微波频率测量装置及方法
CN107917877B (zh) 一种光纤氢气传感器解调方法
CN108955940B (zh) 一种光纤光栅温度传感解调方法
CN107402118B (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