CN101718563A - 基于3×3光纤耦合器的相移白光干涉测量方法 - Google Patents

基于3×3光纤耦合器的相移白光干涉测量方法 Download PDF

Info

Publication number
CN101718563A
CN101718563A CN200910249630A CN200910249630A CN101718563A CN 101718563 A CN101718563 A CN 101718563A CN 200910249630 A CN200910249630 A CN 200910249630A CN 200910249630 A CN200910249630 A CN 200910249630A CN 101718563 A CN101718563 A CN 101718563A
Authority
CN
China
Prior art keywords
interferometer
road
phase
signal
phase information
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.)
Pending
Application number
CN200910249630A
Other languages
English (en)
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.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN200910249630A priority Critical patent/CN101718563A/zh
Publication of CN101718563A publication Critical patent/CN101718563A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Instruments For Measurement Of Length By Optical Means (AREA)

Abstract

本发明涉及光纤白光干涉测量方法,特别涉及基于3×3光纤耦合器的相移白光干涉测量方法,属于光纤传感技术领域。该方法由一只2×2光纤耦合器和一只3×3光纤耦合器构成光纤M-Z干涉仪,波长扫描光注入到光纤M-Z干涉仪中。输出的三路白光干涉光谱信号幅值相等,在相位上两两相差120°。当波长扫描光从λ1扫描到λ2时,得到三路输出信号I1,I2和I3,通过计算或构造
Figure 200910249630.1_AB_2
两路正交的信号,再利用求导交叉相乘方法解调出δ;计算出δ后再对其进行解包裹运算解调出线性相位信息
Figure 200910249630.1_AB_3
。再利用
Figure 200910249630.1_AB_4
就可以计算出干涉仪的绝对光程差。本发明能够高精度、高速度、绝对测量光纤干涉仪的光程差。

Description

基于3×3光纤耦合器的相移白光干涉测量方法
技术领域
本发明涉及光纤白光干涉测量方法,特别涉及基于3×3光纤耦合器的相移白光干涉测量方法,属于光纤传感技术领域。
背景技术
干涉型光纤传感器,如法布里-珀罗(Fabry-Perot,法珀)、马赫-曾德尔(Mach-Zehnder)、迈克尔逊(Michelson)、萨格纳克(Sagnac)等,具有很高的测量灵敏度,得到广泛的应用。待测量,如应变、温度、位移或者压力等,作用在干涉仪上,引起干涉仪光程差的变化,从干涉仪的输出干涉信号中解调出光程差,从而实现对被测物理量的测量。
干涉测量法是测量很多参数(如温度,压力,应变等)非常重要的一种方法,经常使用干涉仪的一个臂作为参考光(未经过含有所需测量参数的环境),另一个臂作为信号光(经过含有所需测量参数的环境),解调出两臂之间的相位差,就可以获得需要测量的参数。因此干涉信号的解调技术是光纤干涉测量术的核心技术。干涉解调技术主要分为相对测量技术和绝对测量技术。相对测量技术只能测量信号的变化量,不能测量静态或缓变信号,而绝对测量可以测量出干涉仪的光程差,因而可以测量静态和缓变信号。
绝对测量一般使用光纤白光干涉测量法(WLI)。传统的白光干涉测量法使用扫描干涉仪来补偿光程,但这样的***稳定性差,分辨率低。光谱域白光干涉测量法利用宽带光源或波长扫描光源,探测传感干涉仪的输出白光光谱,通过光谱中相位相差2π的两个波长来绝对测量干涉仪的光程差。若使用宽带光注入干涉仪,则用光谱分析仪来探测干涉仪的白光光谱;若使用波长扫描光注入干涉仪,则用一只光电二极管就可以获得干涉仪的白光光谱。目前已报道了一些光谱域白光干涉测量法,比如傅里叶变换白光干涉测量法(Yijiang,可得一个正弦信号
Figure G2009102496301D0000021
2倍的第2路信号减去1和3路信号之和后除以系数3可得到一个余弦信号
Figure G2009102496301D0000022
与上面所得到的g1信号相位相差90°;将g1和g2这两路正交信号相除后取反余弦可得到干涉仪的相位信息 δ = arctg ( - 3 I 1 - I 3 2 I 2 - I 1 - I 3 ) .
4)步骤3中提取出相位信息δ是通过反正切算法得到的,相位范围在
Figure G2009102496301D0000024
之间,不连续。为了获取连续相位信息,需要对相位进行解包裹运算。方法是:先给解调出的相位信息δ都乘以2,将其范围扩展到【-π,π】之间,再进行相位解包裹运算,然后再除以2,恢复到原始的线性相位信息
Figure G2009102496301D0000025
根据线性相位信息
Figure G2009102496301D0000026
与波长之间的微分关系
Figure G2009102496301D0000027
干涉仪的光程差D就可通过下式直接得到:
Figure G2009102496301D0000028
如上所述的基于3×3光纤耦合器的相移白光干涉测量方法,还可通过下述方法实现:
方法二,具体实现步骤如下:
a)与方法一中的步骤1)和2)相同;
b)将步骤a中得到的三路干涉信号中的1和3两路进行如下处理可得一个正弦信号
Figure G2009102496301D0000029
2倍的第2路信号减去1和3路信号之和后除以系数3可得到一个余弦信号
Figure G2009102496301D00000210
与上面所得到的g1信号相位相差90°;分别对g1和g2求导后可得g1=b(λ)δ′cosδ,g2′=-b(λ)δ′(t)sinδ;将g1、g2、g1、“Fourier Transform White_Light Interferometry for theMeasurement of Fiber_Optic Extrinsic Fabry_PerotInterferometric Sensors,”Opt Lett.,Vol.33,NO.16),软件解调法(Yijiang,“Wavelength_scanning white-light interferometrywith a 3*3coupler based interferometer,”OptLett.,Vol.20,NO.2),这些方法测量精度高,动态范围大,但是算法比较复杂,测量一次时间约为130us。
发明内容
本发明的目的是为了更好的解决快速、高效解调干涉仪绝对光程差的问题,提出一种简单高速的方法来测量光纤干涉仪绝对光程差,具有精度高、速度快、能绝对测量等优点。
本发明的基于3×3光纤耦合器的相移白光干涉测量方法,是通过下述技术方案实现的:
方法一,具体实现步骤如下:
1)应用一只2×2耦合器和一只3×3耦合器构成光纤M-Z干涉仪;将波长扫描光入射到干涉仪,使干涉仪的输出端得到三路白光干涉光谱。
2)若3×3耦合器对称,干涉仪的三路输出端互成120°。定义3×3耦合器输出的三路干涉信号为:
Figure G2009102496301D0000031
其中,定义K=1时为1路信号,K=2时为2路信号,K=3时为3路信号,a(λ)为光源光谱轮廓引入的背景光(即直流信号),b(λ)是干涉条纹的对比度,线性相位信息
Figure G2009102496301D0000032
(式中m取值为整数,-π≤δ≤π),D是干涉仪两臂之间的光程差。
3)将步骤2中得到的三路干涉信号中的1和3两路做如下处理g2′交叉相乘后再相减得到x=g2g1′-g1g2′=b2(λ)(sin2δ+cos2δ)δ′=b2(λ)δ′;将g1和g2分别平方后相加得y=b2(λ);x除以y后得到
Figure G2009102496301D0000041
对z积分后即可解调出相位信息δ;
c)对步骤b中得到的相位信息δ进行相位解包裹运算,恢复原始线性相位信息
Figure G2009102496301D0000042
线性相位信息
Figure G2009102496301D0000043
是光源波长从λ1扫描到λ2所引起的干涉仪的相位变化,根据相位与波长之间的微分关系
Figure G2009102496301D0000044
干涉仪的光程差D就可通过(2)式直接得到:
通过上述两种方法实现了绝对测量光程差的目的。
有益效果:
本发明方法所提出的基于3×3耦合器的相移白光干涉测量术具有灵敏度高,精度高,动态范围大,算法简单,速度快等特点,满足高速绝对测量的要求。
附图说明
图1为本发明方法应用于基于3×3耦合器的白光扫描干涉测量的一个具体实施方案;
图2为三路白光干涉光谱信号中的一路;
图3为用来校正波长的标准具的光谱;
图4为方法一解调出来随波长变化的线性相位信息
Figure G2009102496301D0000046
图5为构造出来的I1-I3和2I2-I1-I3两路正交信号的李萨茹图;
图6为方法二解调出来的线性相位信息
Figure G2009102496301D0000047
图7为本发明的两种方法解调出的相位与傅立叶变换白光干涉测量法的比较;
图中,1-宽带光源、2-光隔离器、3-光纤可调谐法珀滤波器、4-2×2耦合器、5-弹性柱、6-3×3耦合器、7-光电探测器、8-信号采集与处理***、9-锯齿波发生器、10-法珀标准具、11-光纤光栅、12-光电探测器。
具体实施方案
下面结合附图和实施例对本发明做进一步说明。
本发明的测量方法通过实验进行了验证。图1是本发明方法的一个实施原理图,结合Mach-Zehnder干涉仪的两臂之间光程差的测量对本发明方法进行说明。
宽带光源1发出的光经光隔离器2输入到光纤可调谐法珀滤波器3,光纤可调谐法珀滤波器3在锯齿波发生器9的驱动下输出波长扫描窄带光,波长扫描窄带光被2×2耦合器4分成两路,其中一路光通过法珀标准具10、光纤光栅11,光电探测器12,用以探测法珀标准具10和光纤光栅11的透射谱;另外一路光注入干涉仪;干涉仪中,一路光作为参考光,另一路作为信号光(包含有待测量信息的光),两路光在3×3耦合器6处发生干涉,三路输出的干涉信号经过光电探测器7转化为电信号,信号采集与处理***8用来采集探测器7和探测器12的输出信号,并进行处理,同时控制锯齿波发生器9输出驱动信号。
图1实例中,宽带光源1为ASE光源,光谱范围覆盖从1525nm到1565nm,光纤可调谐法珀滤波器3的自由光谱区(FSR)宽度为65nm,精细度为200,所以带宽为0.325nm;法珀标准具10的自由光谱区宽度为0.8nm(100GHz),精细度为14,法珀标准具10的波长热稳定性为从0℃到70℃的变化小于0.7GHz;一个带宽为0.7nm的中心波长为1535.027nm的光纤光栅11被串联到法珀标准具10后面,用于抹掉法珀标准具10输出梳状光谱中的一个峰值波长,作为波长标记,而这个被抹掉的峰值波长就是光纤光栅的11中心波长。法珀标准具10用于将采集到的干涉信号从时间域转化为波长域,进而取代笨重昂贵的光谱分析仪。
采集的三路白光干涉光谱的幅值相等,相位相互相差120°。其中一路的白光光谱如图2所示,实际上就是在ASE光源的轮廓上叠加了一个白光干涉条纹,所以在处理信号的时候,先对干涉信号做归一化处理,除掉光源轮廓的影响,减小光源波动带给测量的误差。波长校正时,选取了法珀标准具10上1525.649nm和1564.676nm这两个透射峰值点之间的一段信号对干涉信号来进行波长校正,法珀标准具10的输出信号如图3所示,同时确定了波长扫描的起止波长为1525.649nm和1564.676nm。
然后分别利用上述的两种方法进行相位解调。利用方法一解调出来的相位信息如图4所示。利用方法二,构造出来的两路正交信号的李萨茹图如图5所示,从图中可看出获得的两路信号是正交的,解调出来的相位信息如图6所示。
与江毅等提出的“傅里叶变换白光干涉术”(专利申请号:200710177837.3)进行比较,三种方案解调出来的相位随扫描波长变化如图7所示,使用公式(2),可以计算出来的干涉仪的绝对光程差为2996.19um,运行一次时间约为30us。可见本申请提出的方法获得的结果与其它方法测量得到的结果完全相同,但是运行速度却远远比其它方法提高了四五倍,更易于实现实时高速测量的要求。

Claims (2)

1.)一种基于3×3光纤耦合器的相移白光干涉测量方法,其特征在于具体实现步骤如下:
1)应用一只2×2耦合器和一只3×3耦合器构成光纤M-Z干涉仪;将波长扫描光入射到干涉仪,使干涉仪的输出端得到三路白光干涉光谱;
2)若3×3耦合器对称,干涉仪的三路输出端互成120°,定义3×3耦合器输出的三路干涉信号为:
k=1,2,3
其中,定义K=1时为1路信号,K=2时为2路信号,K=3时为3路信号,a(λ)为光源光谱轮廓引入的背景光即直流信号,b(λ)是干涉条纹的对比度,线性相位信息
Figure F2009102496301C0000012
式中m取值为整数,-π≤δ≤π,D是干涉仪两臂之间的光程差;
3)将步骤2)中得到的三路干涉信号中的1和3两路做如下处理可得一个正弦信号
Figure F2009102496301C0000013
2倍的第2路信号减去1和3路信号之和后除以系数3可得到一个余弦信号与上面所得到的g1信号相位相差90°;将g1和g2这两路正交信号相除后取反余弦可得到干涉仪的相位信息
Figure F2009102496301C0000015
4)步骤3)中提取出相位信息δ是通过反正切算法得到的,相位范围在
Figure F2009102496301C0000016
之间,不连续;为了获取连续相位信息,需要对相位进行解包裹运算,方法是:先给解调出的相位信息δ都乘以2,将其范围扩展到【-π,π】之间,再进行相位解包裹运算,然后再除以2,恢复到原始的线性相位信息
Figure F2009102496301C0000021
根据线性相位信息
Figure F2009102496301C0000022
与波长之间的微分关系
Figure F2009102496301C0000023
干涉仪的光程差D就可通过下式直接得到:
2.根据权利要求1所述的一种基于3×3光纤耦合器的相移白光干涉测量方法,其特征在于还可通过下述方法实现:
a)前两个步骤同权利要求1的1)和2);
b)将步骤a)中得到的三路干涉信号中的1和3两路进行如下处理可得一个正弦信号
Figure F2009102496301C0000025
2倍的第2路信号减去1和3路信号之和后除以系数3可得到一个余弦信号
Figure F2009102496301C0000026
与上面所得到的g1信号相位相差90°;分别对g1和g2求导后可得g1′=b(λ)δ′cosδ,g2′=-b(λ)δ′(t)sinδ;将g1、g2、g1′、g2′交叉相乘后再相减得到x=g2g1′-g1g2′=b2(λ)(sin2δ+cos2δ)δ′=b2(λ)δ′;将g1和g2分别平方后相加得y=b2(λ);x除以y后得到
Figure F2009102496301C0000027
对z积分后即可解调出相位信息δ;
c)对步骤b)中得到的相位信息δ进行相位解包裹运算,恢复原始线性相位信息
Figure F2009102496301C0000028
根据相位与波长之间的微分关系
Figure F2009102496301C0000029
干涉仪的光程差D为
CN200910249630A 2009-12-09 2009-12-09 基于3×3光纤耦合器的相移白光干涉测量方法 Pending CN101718563A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN200910249630A CN101718563A (zh) 2009-12-09 2009-12-09 基于3×3光纤耦合器的相移白光干涉测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200910249630A CN101718563A (zh) 2009-12-09 2009-12-09 基于3×3光纤耦合器的相移白光干涉测量方法

Publications (1)

Publication Number Publication Date
CN101718563A true CN101718563A (zh) 2010-06-02

Family

ID=42433170

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200910249630A Pending CN101718563A (zh) 2009-12-09 2009-12-09 基于3×3光纤耦合器的相移白光干涉测量方法

Country Status (1)

Country Link
CN (1) CN101718563A (zh)

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102889900A (zh) * 2012-09-11 2013-01-23 天津大学 基于相位斜率确定单色频率干涉级次的低相干干涉解调方法
CN103033202A (zh) * 2012-05-18 2013-04-10 天津大学 一种相移式高速低相干干涉解调装置及方法
CN104457583A (zh) * 2014-11-24 2015-03-25 上海光亮光电科技有限公司 一种光纤长度测量装置及方法
CN106197741A (zh) * 2016-07-14 2016-12-07 盐城工学院 基于微纳长周期光纤光栅传感器的温度检测装置及方法
CN107765090A (zh) * 2017-10-16 2018-03-06 上海传输线研究所(中国电子科技集团公司第二十三研究所) 一种3×3耦合器相位次序标定装置及方法
CN109084883A (zh) * 2018-08-09 2018-12-25 合肥工业大学 基于相位-botdr光纤分布式布里渊振动传感测量方法
CN109781240A (zh) * 2019-01-30 2019-05-21 武汉理工大学 基于双波长回归分析的光纤振动传感器相位解调装置与方法
CN110243763A (zh) * 2019-07-02 2019-09-17 东北大学 一种非接触光声成像装置及方法
CN110412038A (zh) * 2019-07-17 2019-11-05 天津大学 一种基于单光纤光栅和神经网络的结构损伤位置识别***
CN111121661A (zh) * 2019-12-17 2020-05-08 南京理工大学 用于光滑表面形貌测量的窄带非单色光n+1幅相移测试算法
CN112834072A (zh) * 2021-02-08 2021-05-25 广东海洋大学 一种检测条纹对比度变化的迈克尔逊干涉光纤温度传感器
CN113419341A (zh) * 2021-06-22 2021-09-21 中国人民解放军国防科技大学 一种光纤干涉仪光场的复数化重构方法及装置
CN113654766A (zh) * 2021-08-13 2021-11-16 北京航空航天大学 单模光纤耦合器相移测量装置、方法、存储介质
CN114440947A (zh) * 2021-12-22 2022-05-06 北京自动化控制设备研究所 基于3×3耦合器的干涉信号相位解调方法及***
CN114739432A (zh) * 2022-04-07 2022-07-12 合肥工业大学 基于几何移相器的双波长正交相位偏置锁定干涉测量方法
CN116907627A (zh) * 2023-09-13 2023-10-20 之江实验室 基于光程差辅助的大动态范围分布式相位传感方法和装置
CN114739432B (zh) * 2022-04-07 2024-06-28 合肥工业大学 基于几何移相器的双波长正交相位偏置锁定干涉测量方法

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103033202A (zh) * 2012-05-18 2013-04-10 天津大学 一种相移式高速低相干干涉解调装置及方法
CN103033202B (zh) * 2012-05-18 2015-06-03 天津大学 一种相移式高速低相干干涉解调装置及方法
CN102889900A (zh) * 2012-09-11 2013-01-23 天津大学 基于相位斜率确定单色频率干涉级次的低相干干涉解调方法
CN102889900B (zh) * 2012-09-11 2015-10-28 天津大学 相位斜率确定单色频率干涉级次的低相干干涉解调方法
CN104457583A (zh) * 2014-11-24 2015-03-25 上海光亮光电科技有限公司 一种光纤长度测量装置及方法
CN106197741A (zh) * 2016-07-14 2016-12-07 盐城工学院 基于微纳长周期光纤光栅传感器的温度检测装置及方法
CN106197741B (zh) * 2016-07-14 2018-08-28 盐城工学院 基于微纳长周期光纤光栅传感器的温度检测装置及方法
CN107765090A (zh) * 2017-10-16 2018-03-06 上海传输线研究所(中国电子科技集团公司第二十三研究所) 一种3×3耦合器相位次序标定装置及方法
CN109084883A (zh) * 2018-08-09 2018-12-25 合肥工业大学 基于相位-botdr光纤分布式布里渊振动传感测量方法
CN109781240A (zh) * 2019-01-30 2019-05-21 武汉理工大学 基于双波长回归分析的光纤振动传感器相位解调装置与方法
CN110243763A (zh) * 2019-07-02 2019-09-17 东北大学 一种非接触光声成像装置及方法
CN110412038A (zh) * 2019-07-17 2019-11-05 天津大学 一种基于单光纤光栅和神经网络的结构损伤位置识别***
CN111121661A (zh) * 2019-12-17 2020-05-08 南京理工大学 用于光滑表面形貌测量的窄带非单色光n+1幅相移测试算法
CN112834072A (zh) * 2021-02-08 2021-05-25 广东海洋大学 一种检测条纹对比度变化的迈克尔逊干涉光纤温度传感器
CN112834072B (zh) * 2021-02-08 2021-09-24 广东海洋大学 一种检测条纹对比度变化的迈克尔逊干涉光纤温度传感器
CN113419341A (zh) * 2021-06-22 2021-09-21 中国人民解放军国防科技大学 一种光纤干涉仪光场的复数化重构方法及装置
CN113654766A (zh) * 2021-08-13 2021-11-16 北京航空航天大学 单模光纤耦合器相移测量装置、方法、存储介质
CN113654766B (zh) * 2021-08-13 2022-04-08 北京航空航天大学 单模光纤耦合器相移测量装置、方法、存储介质
CN114440947A (zh) * 2021-12-22 2022-05-06 北京自动化控制设备研究所 基于3×3耦合器的干涉信号相位解调方法及***
CN114739432A (zh) * 2022-04-07 2022-07-12 合肥工业大学 基于几何移相器的双波长正交相位偏置锁定干涉测量方法
CN114739432B (zh) * 2022-04-07 2024-06-28 合肥工业大学 基于几何移相器的双波长正交相位偏置锁定干涉测量方法
CN116907627A (zh) * 2023-09-13 2023-10-20 之江实验室 基于光程差辅助的大动态范围分布式相位传感方法和装置
CN116907627B (zh) * 2023-09-13 2023-12-19 之江实验室 基于光程差辅助的大动态范围分布式相位传感方法和装置

Similar Documents

Publication Publication Date Title
CN101718563A (zh) 基于3×3光纤耦合器的相移白光干涉测量方法
US11181400B2 (en) Fiber Bragg Grating demodulation device capable of supressing fluctuations at variable ambient temperature and demodulation method thereof
CN100541148C (zh) 一种光纤傅立叶变换白光干涉测量法
CN107024236B (zh) F-p/fbg光纤传感器解调***
CN101458101A (zh) 一种基于3×3耦合器的波长扫描白光干涉测量方法
CN108827601A (zh) 一种光纤干涉仪臂长差的测量装置
CN103900797B (zh) 带有光程扫描位置和速度校正的光学相干域偏振测量装置
Jiang et al. Recent developments in fiber optic spectral white-light interferometry
CN113108710B (zh) 基于椭圆拟合的光学低频应变检测***与检测方法
CN108120378A (zh) 基于飞秒光频梳的正弦相位调制干涉绝对测距装置与方法
CN103900798A (zh) 一种带有光程扫描在线校正的光学相干域偏振测量装置
CN108332785A (zh) 一种大规模光纤光栅传感器的测量装置和方法
US11619783B2 (en) Sine-cosine optical frequency detection devices for photonics integrated circuits and applications in LiDAR and other distributed optical sensing
CN108844717A (zh) 一种光纤干涉仪臂长差的测量方法
CN108955939B (zh) 一种光纤光栅温度传感解调***
CN209783610U (zh) 基于色散补偿光栅对的光电振荡器双参量传感装置
Wang et al. Wide measurement range distributed strain sensing with phase-accumulation optical frequency domain reflectometry
CN105091939B (zh) 一种用于光纤Sagnac干涉仪传感器的高分辨率绝对相位解调方法
CN101592526A (zh) 一种光平均波长的测量方法及装置
Jiang et al. Phase-shifted white-light interferometry for the absolute measurement of fiber optic Mach–Zehnder interferometers
CN101169339B (zh) 一种光纤傅立叶变换白光干涉相对测量法
CN105865500A (zh) 一种光纤激光干涉型传感器的移相解调法
US20220247498A1 (en) Multi-point self-calibration for broadband optical sensor interrogator
CN111238553B (zh) 一种光纤光栅传感器波长解调方法及装置
Chen et al. Modified frequency-shifted interferometer: encoding wavelength into phase

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Open date: 20100602