CN113419341A - 一种光纤干涉仪光场的复数化重构方法及装置 - Google Patents

一种光纤干涉仪光场的复数化重构方法及装置 Download PDF

Info

Publication number
CN113419341A
CN113419341A CN202110688836.5A CN202110688836A CN113419341A CN 113419341 A CN113419341 A CN 113419341A CN 202110688836 A CN202110688836 A CN 202110688836A CN 113419341 A CN113419341 A CN 113419341A
Authority
CN
China
Prior art keywords
signal
interferometer
optical fiber
optical
pulse
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
CN202110688836.5A
Other languages
English (en)
Other versions
CN113419341B (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.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
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 National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN202110688836.5A priority Critical patent/CN113419341B/zh
Publication of CN113419341A publication Critical patent/CN113419341A/zh
Application granted granted Critical
Publication of CN113419341B publication Critical patent/CN113419341B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G02—OPTICS
    • G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/0012—Optical design, e.g. procedures, algorithms, optimisation routines

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optics & Photonics (AREA)
  • Instruments For Measurement Of Length By Optical Means (AREA)

Abstract

本发明涉及光学领域，涉及到光纤干涉仪信号解调、光场变换、信号串扰、偏振抑制等技术内容，具体为一种光纤干涉仪光场的复数化重构方法及装置，所述方法包括：S1：判定问询脉冲信号形式；S2：根据问询脉冲信号形式确认干涉信号采样方式；S3:干涉光场复数化重构。本发明可将光电探测器实际采集的光纤光干涉仪干涉信号准确重构出复数表达形式，这种干涉信号复数表达形式可更加便捷高效地应用于信号基本运算过程，是偏振合成、串扰抑制等技术实现必须的信号处理环节。

Description

一种光纤干涉仪光场的复数化重构方法及装置
技术领域
本发明涉及光学领域，涉及到光纤干涉仪信号解调、光场变换、信号串扰、偏振抑制等技术内容，具体为一种光纤干涉仪光场的复数化重构方法及装置。
背景技术
光纤干涉仪在光纤陀螺、光纤水声探测、地震波检测等领域有重要的应用。由于采用光学相干检测方式，具有灵敏度高、动态范围大等突出优势，结合光纤传感***自身体积小、重量轻、抗电磁干扰等固有优势优点，在民用和军用领域都得到了广泛应用。但是，由于光纤干涉仪灵敏度高，导致很多环境扰动、光路串扰等都会影响相位，通过各种信号处理手段去除扰动因素影响是基础技术之一，且信号处理技术一般较为复杂。
在很多干涉信号处理过程中，往往需要将多个干涉结果进行乘除法运算达到相位信息相消或相长的目的。文献“Reduction of crosstalk in inline sensor arraysusing inverse scattering”(Proc.SPIE,2008,Vol.704,paper 70044Z)和专利“Methodand apparatus for suppression of crosstalk and noise in time-divisionmultiplexed interferometer sensor systems(专利号:US7466422B2)”报道了一种反转剥层算法，这种剥层算法需要将两个干涉信号进行减法运算来消除不同干涉通道之间的串扰。文献“Low-Crosstalk and Polarization-Independent Inline InterferometricFiber Sensor Array Based on Fiber Bragg Gratings.”(J.Lightw.Technol.34(2016)4232-4239)报道了一种偏振合成算法，这种算法需要将四个干涉信号进行乘法和减法运算来消除不同干涉通道之间的偏振影响因素。在实际中，光纤干涉仪的干涉结果在采用光电探测器转换后获取的为实数表达结果，可写为(公式1)
Figure BDA0003125700870000011
I(t)表示干涉仪输出的光强度，A表示直流量，B表示交流量幅度，Re表示取复数的实部运算，
Figure BDA0003125700870000012
表示干涉仪的信号项，
Figure BDA0003125700870000013
为干涉仪的噪声。当涉及到干涉光场的加减乘除运算时，直接采用(公式1)计算时会出现多余频率项。对于两个干涉信号I1(t)和I2(t)
Figure BDA0003125700870000014
Figure BDA0003125700870000015
期望获得干涉信号的和频或差频项，即获得
Figure BDA0003125700870000016
或
Figure BDA0003125700870000017
但在实际运算中，直接将I1(t)和I2(t)进行加减乘除运算，无论哪种都会产生多余的相位项，也多与相位项与期望相位项无法分离。例如，进行混频相乘则获得：
Figure BDA0003125700870000021
很容易可以推出，进行加法、减法和除法，都会是多个频率项的耦合结果。一种简单的处理方法是，若能够获取干涉光场的复数表达形式，即
Figure BDA0003125700870000022
Figure BDA0003125700870000023
将两个干涉光场的复数表达形式直接进行乘法和除法运算，即可获得
Figure BDA0003125700870000024
和
Figure BDA0003125700870000025
再取实部运算即可得到
Figure BDA0003125700870000026
和
Figure BDA0003125700870000027
可见干涉光场的复数表达形式在很多应用中是简化算法、提高计算效率的关键，在光纤干涉仪的复杂信号处理中往往起到化繁为简甚至成为确保算法正确性的基础。
由于干涉光场的复杂性，信号往往加载在余弦函数的相位项中，直接通过数学变换由余弦项获得正弦项会带来频谱成分的移相，对干涉光场的复数化重构需要根据干涉光场的信号特征来进行特殊设计。
发明内容
本发明的目的是提供一种光纤干涉仪光场的复数化重构方法及装置，在实际中，光纤干涉仪的干涉结果在采用光电探测器转换后获取的为实数表达结果，本发明通过将根据实际检测的干涉光场实数表达形式重构为复数表达形式，为光纤干涉仪传感***的串扰抑制、偏振抑制等提供可快速处理的信号基础。
为了实现上述技术目的，本发明采用的技术方案为：
一种光纤干涉仪光场的复数化重构方法，包括以下步骤：
S1：判定问询脉冲信号形式：若是问询脉冲信号为一对含有固定频差fm的脉冲，则判定为外差调制信号；若是问询脉冲为三个相邻相位差为120度的脉冲对，则判定为3×3调制信号；
S2：根据问询脉冲信号形式确认干涉信号采样方式：当问询脉冲信号为外差调制信号时，采集一组干涉结果，定义为IH(t)；当问询脉冲信号为3×3调制信号时，采集三组干涉结果，分别记为I1(t)、I2(t)和I3(t)；
S3:干涉光场复数化重构：
S3.1针对S2采集的一组干涉结果IH(t)，采取以下步骤：
S3.1.1将IH(t)分别与cos(2πfmt)、sin(2πfmt)混频，获得两个表达式IC和IS，即：
Figure BDA0003125700870000028
S3.1.2将IC和IS分别进行低通滤波，保留小于fm的频谱成分，获得表达式IC_low和IS_low：
Figure BDA0003125700870000031
Figure BDA0003125700870000032
其中，B表示交流量幅度，
Figure BDA0003125700870000033
表示干涉仪的信号项，
Figure BDA0003125700870000034
为干涉仪的噪声；
S3.1.3生成相位项复光场表达形式Icom：
Figure BDA0003125700870000035
S3.2针对S2采集的三组干涉结果I1(t)、I2(t)和I3(t)采取以下步骤：
S3.2.1将I1(t)、I2(t)和I3(t)相加后取平均，获取直流项：
A＝{I1(t)+I2(t)+I3(t)}/3
S3.2.2将I1(t)、I2(t)和I3(t)分别减去直流项，获取三个交流项，即：
Figure BDA0003125700870000036
Figure BDA0003125700870000037
Figure BDA0003125700870000038
S3.2.3将
Figure BDA0003125700870000039
减去
Figure BDA00031257008700000310
获得正弦项：
Figure BDA00031257008700000311
S3.2.4合成含有干涉仪相位项的干涉结果复数表达形式Icom：
Figure BDA00031257008700000312
优选的，对于含有外差调制信号的问询脉冲信号，通过光学移频方式产生。
优选的，对于含有外差调制信号的问询脉冲信号，通过采用三角波调制光源频率产生。
优选的，对于含有3×3调制信号的问询脉冲信号，通过三组相差分别为
Figure BDA00031257008700000313
0和
Figure BDA00031257008700000314
的脉冲对产生。
优选的，对于含有3×3调制信号的问询脉冲信号，通过3×3耦合器产生。
本发明还提供一种基于上述方法的光纤干涉仪光场的复数化重构装置,包括调制信号加载模块1和信号处理模块2；
所述调制信号加载模块1用于产生问询干涉仪的脉冲序列，并将脉冲序列注入到光纤干涉仪中，包括问询脉冲模块101和光纤干涉仪102，所述问询脉冲模块101用于产生问询干涉仪的脉冲序列，并将脉冲序列注入到光纤干涉仪102中，所述光纤干涉仪102用于产生干涉光场，包括耦合器102A、法拉第旋镜102B和102C，所述问询脉冲模块101通过光纤连接至光纤干涉仪102的输入端口102A1；
所述信号处理模块2用于将光纤干涉仪102产生的干涉光场进行光电转换后进行复数化重构处理，包括光电转换模块201、采集模块202和PC机203；所述光电转换模块201用于光电转换，所述采集模块202用于对光电转换后的电信号进行采样，所述PC机203用于将采样信号存储、显示和复数化重构；所述光电转换模块201通过电缆连接到采集模块202，所述采集模块202通过信号线连接到所述PC机203。
所述制信号加载模块1和所述信号处理模块2之间采用光纤连接，连接方式为所述光纤干涉仪102的输出端口102A2连接到所述光电转换模块201。
本发明具有以下有益效果：本发明可将光电探测器实际采集的光纤光干涉仪干涉信号准确重构出复数表达形式，这种干涉信号复数表达形式可更加便捷高效地应用于信号基本运算过程，是偏振合成、串扰抑制等技术实现必须的信号处理环节。
附图说明
图1为本发明中光纤干涉仪光场的复数化重构方法流程图；
图2为本发明中光纤干涉仪光场的复数化重构装置的组成结构图。
具体实施方式
下面结合附图对本发明作进一步说明，如附图1所示，本发明所述中光纤干涉仪光场的复数化重构方法针对两种不同调制手段的光纤干涉仪光场进行不同的处理包含如下步骤：
S1：判定问询脉冲信号形式
若是问询脉冲信号为一对含有固定频差fm的脉冲，则判定为外差调制信号；若是问询脉冲为三个相邻相位差为120度的脉冲，则判定为3×3调制信号；
S2：根据问询脉冲信号形式确认干涉信号采样方式：当问询脉冲信号为外差调制信号时，采集一组干涉结果，定义为IH(t)，可表达为
Figure BDA0003125700870000041
IH(t)表示经过外差调制后干涉仪输出的光强度，A表示直流量，B表示交流量幅度，Re表示取复数的实部运算，2πfmt为外差调制项，其中fm为调制频率，
Figure BDA0003125700870000042
表示干涉仪的信号项，
Figure BDA0003125700870000043
为干涉仪的噪声。
当问询脉冲信号为3×3调制信号时，采集三组干涉结果，分别记为I1(t)、I2(t)和I3(t)，可表达为
Figure BDA0003125700870000044
Figure BDA0003125700870000045
Figure BDA0003125700870000046
I1(t)、I2(t)和I3(t)表示经过3×3调制干涉仪输出的光强度，A表示直流量，B表示交流量幅度，Re表示取复数的实部运算，
Figure BDA0003125700870000047
表示干涉仪的信号项，
Figure BDA0003125700870000048
为干涉仪的噪声。
Figure BDA0003125700870000049
是干涉仪本身的相位，I1(t)、I2(t)和I3(t)之间由于有3×3调制分别叠加了
Figure BDA00031257008700000410
0和
Figure BDA00031257008700000411
相位。
S3:干涉光场复数化重构：
S3.1针对S2采集的一组干涉结果IH(t)，采取以下步骤：
S3.1.1由于外差调制项2πfmt为人为引入项，将IH(t)分别与cos(2πfmt)、sin(2πfmt)混频，获得两个表达式IC和IS，其中：
Figure BDA0003125700870000051
S3.2.2将IC和IS分别进行低通滤波，保留小于fm的频谱成分，获得：
Figure BDA0003125700870000052
Figure BDA0003125700870000053
S3.3.3合成干涉结果的复数表达形式
Figure BDA0003125700870000054
Figure BDA0003125700870000055
S3.2针对S2采集的三组干涉结果I1(t)、I2(t)和I3(t)采取以下步骤：
S3.2.1将I1(t)、I2(t)和I3(t)相加，获取直流项
Figure BDA0003125700870000056
因此有
A＝{I1(t)+I2(t)+I3(t)}/3
S3.2.2将I1(t)、I2(t)和I3(t)分别减去直流项，获取三个交流项，即：
Figure BDA0003125700870000057
Figure BDA0003125700870000058
Figure BDA0003125700870000059
S3.2.3将
Figure BDA00031257008700000510
减去
Figure BDA00031257008700000511
获得正弦项：
Figure BDA0003125700870000061
S3.2.4合成干涉结果复数表达形式Icom
Figure BDA0003125700870000062
可以看出，采用本发明所述的干涉光复光场合成方法，可以获得干涉信号光场的复数表达形式，也即干涉信号的复光场，为偏振控制和抑制多通道串扰等提供信号基础。

Claims (6)

1.一种光纤干涉仪光场的复数化重构方法，其特征在于，包括以下步骤：
S1：判定问询脉冲信号形式：若是问询脉冲信号为一对含有固定频差fm的脉冲，则判定为外差调制信号；若是问询脉冲为三个相邻相位差为120度的脉冲对，则判定为3×3调制信号；
S2：根据问询脉冲信号形式确认干涉信号采样方式：当问询脉冲信号为外差调制信号时，采集一组干涉结果，定义为IH(t)；当问询脉冲信号为3×3调制信号时，采集三组干涉结果，分别记为I1(t)、I2(t)和I3(t)；
S3:干涉光场复数化重构：
S3.1针对S2采集的一组干涉结果IH(t)，采取以下步骤：
S3.1.1将IH(t)分别与cos(2πfmt)、sin(2πfmt)混频，获得两个表达式IC和IS，即：
Figure FDA00031257008600000113
S3.1.2将IC和IS分别进行低通滤波，保留小于fm的频谱成分，获得表达式IC_low和IS_low：
Figure FDA0003125700860000011
Figure FDA0003125700860000012
其中，B表示交流量幅度，
Figure FDA0003125700860000013
表示干涉仪的信号项，
Figure FDA0003125700860000014
为干涉仪的噪声；
S3.1.3生成相位项复光场表达形式Icom：
Figure FDA0003125700860000015
S3.2针对S2采集的三组干涉结果I1(t)、I2(t)和I3(t)采取以下步骤：
S3.2.1将I1(t)、I2(t)和I3(t)相加后取平均，获取直流项：
A＝{I1(t)+I2(t)+I3(t)}/3
S3.2.2将I1(t)、I2(t)和I3(t)分别减去直流项，获取三个交流项，即：
Figure FDA0003125700860000016
Figure FDA0003125700860000017
Figure FDA0003125700860000018
S3.2.3将
Figure FDA0003125700860000019
减去
Figure FDA00031257008600000110
获得正弦项：
Figure FDA00031257008600000111
S3.2.4合成含有干涉仪相位项的干涉结果复数表达形式Icom：
Figure FDA00031257008600000112
2.一种根据权利要求1所述光纤干涉仪光场的复数化重构方法，其特征在于：对于含有外差调制信号的问询脉冲信号，通过光学移频方式产生。
3.一种根据权利要求2所述光纤干涉仪光场的复数化重构方法，其特征在于：对于含有外差调制信号的问询脉冲信号，通过两个光学移频器产生。
4.一种根据权利要求1所述光纤干涉仪光场的复数化重构方法，其特征在于：对于含有3×3调制信号的问询脉冲信号，通过三组相差分别为-
Figure FDA0003125700860000021
0和
Figure FDA0003125700860000022
的脉冲对产生。
5.一种根据权利要求4所述光纤干涉仪光场的复数化重构方法，其特征在于：对于含有3×3调制信号的问询脉冲信号，通过3×3耦合器产生。
6.一种基于权利要求1至5任一项所述方法的光纤干涉仪光场的复数化重构装置，其特征在于：所述装置包括调制信号加载模块(1)和信号处理模块(2)；
所述调制信号加载模块(1)用于产生问询干涉仪的脉冲序列，并将脉冲序列注入到光纤干涉仪中，包括问询脉冲模块(101)和光纤干涉仪(102)，所述问询脉冲模块(101)用于产生问询干涉仪的脉冲序列，并将脉冲序列注入到光纤干涉仪(102)中，所述光纤干涉仪(102)用于产生干涉光场，包括耦合器(102A)、法拉第旋镜(102B)和(102C)，所述问询脉冲模块(101)通过光纤连接至光纤干涉仪(102)的输入端口(102A1)；
所述信号处理模块(2)用于将光纤干涉仪(102)产生的干涉光场进行光电转换后进行复数化重构处理，包括光电转换模块(201)、采集模块(202)和PC机(203)；所述光电转换模块(201)用于光电转换，所述采集模块(202)用于对光电转换后的电信号进行采样，所述PC机(203)用于将采样信号存储、显示和复数化重构；所述光电转换模块(201)通过电缆连接到采集模块(202)，所述采集模块(202)通过信号线连接到所述PC机(203)；
所述制信号加载模块(1)和所述信号处理模块(2)之间采用光纤连接，连接方式为所述光纤干涉仪(102)的输出端口(102A2)连接到所述光电转换模块(201)。
CN202110688836.5A 2021-06-22 2021-06-22 一种光纤干涉仪光场的复数化重构方法及装置 Active CN113419341B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110688836.5A CN113419341B (zh) 2021-06-22 2021-06-22 一种光纤干涉仪光场的复数化重构方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110688836.5A CN113419341B (zh) 2021-06-22 2021-06-22 一种光纤干涉仪光场的复数化重构方法及装置

Publications (2)

Publication Number Publication Date
CN113419341A true CN113419341A (zh) 2021-09-21
CN113419341B CN113419341B (zh) 2022-11-18

Family

ID=77789693

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110688836.5A Active CN113419341B (zh) 2021-06-22 2021-06-22 一种光纤干涉仪光场的复数化重构方法及装置

Country Status (1)

Country Link
CN (1) CN113419341B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1752836A (zh) * 2005-09-27 2006-03-29 北京理工大学 一种光纤干涉信号解调方法
CN101718563A (zh) * 2009-12-09 2010-06-02 北京理工大学 基于3×3光纤耦合器的相移白光干涉测量方法
CN102353393A (zh) * 2011-07-01 2012-02-15 浙江大学 基于π/2相位调制的干涉型光传感器的正交解调装置
US20140176937A1 (en) * 2011-08-18 2014-06-26 Tiegen Liu Distributed disturbance sensing device and the related demodulation method based on polarization sensitive optical frequency domain reflectometry
CN104932112A (zh) * 2015-06-26 2015-09-23 中国科学院上海光学精密机械研究所 一种实时光场重构结构和重构方法
CN106643813A (zh) * 2016-10-14 2017-05-10 中国人民解放军国防科学技术大学 一种含有相位载波调制信号的干涉光复光场合成方法
CN106679790A (zh) * 2016-12-05 2017-05-17 华南理工大学 一种提高分布式光纤振动传感灵敏度的互相关解调方法
CN108020314A (zh) * 2016-11-01 2018-05-11 北京大学 光纤水听器阵列***和加速度传感器阵列***及测量方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1752836A (zh) * 2005-09-27 2006-03-29 北京理工大学 一种光纤干涉信号解调方法
CN101718563A (zh) * 2009-12-09 2010-06-02 北京理工大学 基于3×3光纤耦合器的相移白光干涉测量方法
CN102353393A (zh) * 2011-07-01 2012-02-15 浙江大学 基于π/2相位调制的干涉型光传感器的正交解调装置
US20140176937A1 (en) * 2011-08-18 2014-06-26 Tiegen Liu Distributed disturbance sensing device and the related demodulation method based on polarization sensitive optical frequency domain reflectometry
CN104932112A (zh) * 2015-06-26 2015-09-23 中国科学院上海光学精密机械研究所 一种实时光场重构结构和重构方法
CN106643813A (zh) * 2016-10-14 2017-05-10 中国人民解放军国防科学技术大学 一种含有相位载波调制信号的干涉光复光场合成方法
CN108020314A (zh) * 2016-11-01 2018-05-11 北京大学 光纤水听器阵列***和加速度传感器阵列***及测量方法
CN106679790A (zh) * 2016-12-05 2017-05-17 华南理工大学 一种提高分布式光纤振动传感灵敏度的互相关解调方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
代欣学: "干涉型光纤水听器信号解调技术研究", 《中国优秀博硕士学位论文全文数据库(硕士)基础科学辑》 *

Also Published As

Publication number Publication date
CN113419341B (zh) 2022-11-18

Similar Documents

Publication Publication Date Title
Chen et al. High-fidelity distributed fiber-optic acoustic sensor with fading noise suppressed and sub-meter spatial resolution
CN103954310B (zh) 一种干涉型光纤传感器的大动态信号解调装置及解调方法
Kemao Windowed Fourier transform for fringe pattern analysis
US5313266A (en) Demodulators for optical fiber interferometers with [3×3] outputs
CN104677398B (zh) 一种基于频漂动态补偿的φ‑otdr增敏方法及装置
CN105067017A (zh) 一种改进的生成载波相位pgc解调方法
Hou et al. Large dynamic range and high sensitivity PGC demodulation technique for tri-component fiber optic seismometer
CN101686087B (zh) 基于3×3耦合器两路输出的波长解调***及方法
Jiang et al. High-fidelity acoustic signal enhancement for phase-OTDR using supervised learning
US4881813A (en) Passive stabilization of a fiber optic nonlinear interferometric sensor
Chen et al. Denoising method of the Φ-OTDR system based on EMD-PCC
CN109631962A (zh) 一种基于pgc方案的多频载波偏振调制解调***及方法
CN113419341B (zh) 一种光纤干涉仪光场的复数化重构方法及装置
US6363034B1 (en) Methods and apparatus for phase angle demodulation
Li et al. Denoising method based upon ECEMD and selective addition for phi-OTDR
CN106643813A (zh) 一种含有相位载波调制信号的干涉光复光场合成方法
Duan et al. Simple lock-in detection technique utilizing multiple harmonics for digital PGC demodulators
Zhang et al. Denoising using 3× 3 coupler demodulation
Zhang et al. Investigation on the maximum signal handling capability of fiber optic interferometric sensor based on the digital heterodyne demodulation scheme
CN111579047A (zh) 一种光纤矢量水听器的信号解调方法
Ren et al. A Robust Phase Unwrapping Algorithm With High Measurement Range Based on Unscented Kalman Filter in Interferometric Optical Fiber Sensing
Liu et al. Novel Method for Improving Temperature Resolution of Fiber Optic Sensor Based on Variational Mode Decomposition
CN114323092B (zh) 一种计算与消除内调制pgc信号检测中伴生调幅的方法
Yang et al. Influence of carrier signal frequency error on the demodulation phase of the fiber optic hydrophone
Wu et al. Phase unwrapping error identification and suppression method in φ-OTDR systems based on PELT-VMD-ARIMA

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