CN107063317A - 一种多芯光纤Bragg光栅曲率传感器的解调方法 - Google Patents

一种多芯光纤Bragg光栅曲率传感器的解调方法 Download PDF

Info

Publication number
CN107063317A
CN107063317A CN201710289611.6A CN201710289611A CN107063317A CN 107063317 A CN107063317 A CN 107063317A CN 201710289611 A CN201710289611 A CN 201710289611A CN 107063317 A CN107063317 A CN 107063317A
Authority
CN
China
Prior art keywords
optical fiber
fiber bragg
mrow
bragg raster
centre wavelength
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
CN201710289611.6A
Other languages
English (en)
Other versions
CN107063317B (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.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong University
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 Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN201710289611.6A priority Critical patent/CN107063317B/zh
Publication of CN107063317A publication Critical patent/CN107063317A/zh
Application granted granted Critical
Publication of CN107063317B publication Critical patent/CN107063317B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D5/00Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
    • G01D5/26Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
    • G01D5/32Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
    • G01D5/34Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
    • G01D5/353Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre
    • G01D5/3537Optical fibre sensor using a particular arrangement of the optical fibre itself
    • G01D5/3538Optical fibre sensor using a particular arrangement of the optical fibre itself using a particular type of fiber, e.g. fibre with several cores, PANDA fiber, fiber with an elliptic core or the like

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Optical Transform (AREA)

Abstract

本发明涉及一种多芯光纤Bragg光栅曲率传感器的解调方法,采用余弦型宽带光源将光纤Bragg光栅中心波长的变化转换为其反射功率的变化,利用推导出的光纤Bragg光栅反射功率与中心波长的关系表达式,通过测量光纤Bragg光栅的反射功率计算出对应的中心波长,实现对光纤Bragg光栅中心波长的解调。多芯光纤Bragg光栅曲率传感器包含两个光纤Bragg光栅,分别位于两个相对于多芯光纤横截面中心对称的纤芯中,利用所提出的FBG中心波长解调方法获得两个光纤Bragg光栅中心波长的变化量,进而得到对应的曲率值。该曲率传感器解调技术具有灵敏度高、实现方式简单和成本低等优点。

Description

一种多芯光纤Bragg光栅曲率传感器的解调方法
技术领域
本发明属于曲率传感器解调技术领域,具体涉及一种多芯光纤Bragg光栅曲率传感器的解调方法。
背景技术
曲率是描述物体形状的一个重要参数,通过对曲率的测量可以了解物体形状的变化趋势。曲率传感器在结构体健康监测、表面形状测量等方面具有广泛的应用前景。
基于光纤光栅的曲率传感器具有体积小、重量轻、抗电磁干扰能力强和准分布式测量等优点,特别适合大型结构体、恶劣环境下的曲率测量(M.J.Gander,“Bendmeasurement using Bragg gratings in multi-core fiber,”Electronics Letters 36(2),120~121,2000.)。近年来,为进一步改善光纤光栅型曲率传感器的性能,基于多芯光纤光栅的曲率传感器引起了科研人员的广泛关注。如图1所示,现有多芯光纤光栅曲率传感器的结构包括宽带光源、光耦合器1、第一环形器2、第二环行器5、多芯光纤4及光谱仪,多芯光纤4的两个中心对称的纤芯中对应刻有第一光纤Bragg光栅和第二光纤Bragg光栅;宽带光源与光耦合器1输入端相连;第一环形器2的a端与光耦合器1的第一输出端相连,b端通过多芯光纤4连接器与刻有第一Bragg光栅的纤芯相连,c端与光谱仪相连;第二环形器的a端与光耦合器1的第二输出端相连,b端通过多芯光纤4连接器与刻有第二Bragg光栅的纤芯相连,c端与光谱仪相连。该多芯光纤光栅曲率传感器在工作时,宽带光源首先由3dB光耦合器1分为两路信号,两路信号光分别经第一环形器2和第二环形器的a端进入第一环形器2和第二环形器,然后分别经第一环形器2和第二环形器的b端进入多芯光纤4连接器,分别被多芯光纤4中的两个Bragg光栅反射,反射信号再次通过多芯光纤4连接器后由第一环形器2和第二环形器的c端分别输出,通过光谱分析仪测量光纤光栅的反射波长变化量,进而得到光纤光栅处的曲率值,即传感器的性能很大程度上是由光谱分析仪的性能(即频谱分辨率)所决定的,而光谱分析仪昂贵的价格直接导致解调成本增高,并且采用光谱分析仪的解调方式复杂。
发明内容
本发明的目的在于提供一种解调方式简单、解调成本低的解调方法。
为达到上述要求,本发明采取的技术方案是:提供一种多芯光纤Bragg光栅曲率传感器的解调方法,采用余弦型宽带光源将光纤Bragg光栅中心波长的变化转换为其反射功率的变化,通过测量光纤Bragg光栅的反射功率值实现对其中心波长的解调,具体包括以下步骤:
S1、所述余弦型宽带光源的输出光谱注入到包含第一光纤Bragg光栅和第二光纤Bragg光栅的多芯光纤中,利用测得的第一光纤Bragg光栅和第二光纤Bragg光栅的反射功率,由于光纤Bragg光栅的反射功率值随其中心波长的不同呈现余弦变化,可以分别得到弯曲时第一光纤Bragg光栅和第二光纤Bragg光栅的中心波长。
S2、将弯曲时的第一光纤Bragg光栅中心波长和无弯曲时的第一光纤Bragg光栅中心波长相减得到第一光纤Bragg光栅的中心波长漂移量;将弯曲时的第二光纤Bragg光栅中心波长和无弯曲时的第二光纤Bragg光栅中心波长相减得到第二光纤Bragg光栅的中心波长漂移量;
S3、根据所述第一光纤Bragg光栅的中心波长漂移量计算得到第一纤芯轴向应力,根据所述第二光纤Bragg光栅的中心波长漂移量计算得到第二纤芯轴向应力;
S4、根据所述第一纤芯轴向应力、第二纤芯轴向应力及两纤芯间距计算得到曲率值。
与现有技术相比,本发明具有以下优点:
(1)采用余弦型宽带光源,将光纤Bragg光栅中心波长的变化转换为其反射功率的变化,利用推导出的光纤Bragg光栅反射功率与中心波长的关系表达式,通过测量光纤Bragg光栅的反射功率计算出对应的中心波长,通过中心波长漂移量计算获得曲率值,实现对曲率传感器的解调,无需使用昂贵的光谱分析仪,使得解调方式大大简化,有效降低解调成本;
(2)本解调方式扩展性强,通过对多对光纤Bragg光栅解调,可以实现基于多芯光纤Bragg光栅的二维曲率传感和形状传感。
附图说明
此处所说明的附图用来提供对本申请的进一步理解,构成本申请的一部分,在这些附图中使用相同的参考标号来表示相同或相似的部分,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:
图1为本发明传感器的结构示意图;
图2为本发明解调方法的流程图。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚,以下结合附图及具体实施例,对本申请作进一步地详细说明。为简单起见,以下描述中省略了本领域技术人员公知的某些技术特征。
一种多芯光纤Bragg光栅曲率传感器的解调方法,采用余弦型宽带光源替换图1中的宽带光源,采用光功率计、光电探测器等能检测到反射功率的仪器替换光谱仪,实现对光纤Bragg光栅中心波长的解调,如图2所示,包括以下步骤:
S1、余弦型宽带光源的输出光谱注入到包含第一光纤Bragg光栅和第二光纤Bragg光栅的多芯光纤中,根据余弦型宽带光源输出光谱S(λ)、第一光纤Bragg光栅的反射谱F1(λ)及光功率计测得的第一反射功率值P1计算得出弯曲时的第一光纤Bragg光栅中心波长λ1,根据余弦型宽带光源输出光谱S(λ)、第二光纤Bragg光栅的反射谱F2(λ)及光功率计测得的第二反射功率值P2计算得出弯曲时的第二光纤Bragg光栅中心波长λ2
计算第一光纤Bragg光栅中心波长λ1和第二光纤Bragg光栅中心波长λ2的公式为:
其中,P为光功率计测得的反射功率值,S0为已知的频谱峰值功率,λI为已知的宽带光源的中心波长,ΔλI为已知的宽带光源的频谱周期;R为已知的峰值反射率,Δλ为已知的光纤Bragg光栅的反射谱的半极大全宽度,λ0为光纤Bragg光栅的反射谱中心波长,即为待计算值。
化简上述公式得到根据上述已知值,可以计算出弯曲时光纤Bragg光栅的反射谱中心波长。
S2、将弯曲时的第一光纤Bragg光栅中心波长λ1和无弯曲时的第一光纤Bragg光栅中心波长λB相减得到第一光纤Bragg光栅的中心波长漂移量ΔλB1;将弯曲时的第二光纤Bragg光栅中心波长λ2和无弯曲时的第二光纤Bragg光栅中心波长λB相减得到第二光纤Bragg光栅的中心波长漂移量ΔλB2
S3、根据第一光纤Bragg光栅的中心波长漂移量ΔλB1计算得到第一纤芯轴向应力εFBG1,根据第二光纤Bragg光栅的中心波长漂移量ΔλB2计算得到第二纤芯轴向应力εFBG2
第一纤芯轴向应力εFBG1和第二纤芯轴向应力εFBG2的计算公式为:
其中,ε为弯曲所引入的纤芯轴向应力,ΔλB为光纤Bragg光栅中心波长漂移量,λB为无弯曲时光纤Bragg光栅的中心波长,Pε为弹光系数。
S4、根据第一纤芯轴向应力、第二纤芯轴向应力及两个刻有Bragg光栅的纤芯的距离计算得到曲率值;
曲率值C的计算公式为:
其中,d为两个刻有Bragg光栅的纤芯距离。
该余弦型宽带光源的消光比大于10dB,大消光比可以提高传感器的灵敏度。
以上实施例仅表示本发明的几种实施方式,其描述较为具体和详细,但并不能理解为对本发明范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明保护范围。因此本发明的保护范围应该以权利要求为准。

Claims (5)

1.一种多芯光纤Bragg光栅曲率传感器的解调方法,其特征在于,采用余弦型宽带光源进行解调,包括以下步骤:
S1、所述余弦型宽带光源的输出光谱注入到包含第一光纤Bragg光栅和第二光纤Bragg光栅的多芯光纤中,利用测得的第一光纤Bragg光栅和第二光纤Bragg光栅的反射功率,分别得到弯曲时第一光纤Bragg光栅和第二光纤Bragg光栅的中心波长。
S2、将弯曲时的第一光纤Bragg光栅中心波长和无弯曲时的第一光纤Bragg光栅中心波长相减得到第一光纤Bragg光栅的中心波长漂移量;将弯曲时的第二光纤Bragg光栅中心波长和无弯曲时的第二光纤Bragg光栅中心波长相减得到第二光纤Bragg光栅的中心波长漂移量;
S3、根据所述第一光纤Bragg光栅的中心波长漂移量计算得到第一纤芯轴向应力,根据所述第二光纤Bragg光栅的中心波长漂移量计算得到第二纤芯轴向应力;
S4、根据所述第一纤芯轴向应力、第二纤芯轴向应力及两纤芯间距计算得到曲率值。
2.根据权利要求1所述的多芯光纤Bragg光栅曲率传感器的解调方法,其特征在于,所述步骤S1中计算光纤Bragg光栅中心波长的公式为:
<mrow> <mi>P</mi> <mo>=</mo> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mi>&amp;infin;</mi> </msubsup> <mi>S</mi> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>)</mo> </mrow> <mi>F</mi> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>&amp;lambda;</mi> </mrow>
其中,P为光功率计测得的反射功率值,S(λ)为余弦型宽带光源的输出光谱,S0为频谱峰值功率,λI为宽带光源的中心波长,ΔλI为宽带光源的频谱周期;F(λ)为光纤Bragg光栅的反射谱,R为峰值反射率,Δλ为反射谱的半极大全宽度,λ0为光纤Bragg光栅的反射谱中心波长,即为待计算值。
3.根据权利要求1所述的多芯光纤Bragg光栅曲率传感器的解调方法,其特征在于,所述步骤S3中计算第一纤芯轴向应力和第二纤芯轴向应力的公式为:
<mrow> <mi>&amp;epsiv;</mi> <mo>=</mo> <mfrac> <mrow> <msub> <mi>&amp;Delta;&amp;lambda;</mi> <mi>B</mi> </msub> </mrow> <mrow> <msub> <mi>&amp;lambda;</mi> <mi>B</mi> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>P</mi> <mi>&amp;epsiv;</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>;</mo> </mrow>
其中,ε为弯曲所引入的纤芯轴向应力,ΔλB为光纤Bragg光栅中心波长漂移量,λB为无弯曲时光纤Bragg光栅的中心波长,Pε为弹光系数。
4.根据权利要求1所述的多芯光纤Bragg光栅曲率传感器的解调方法,其特征在于,所述步骤S4中计算曲率值的公式为:
<mrow> <mi>C</mi> <mo>=</mo> <mfrac> <mrow> <msub> <mi>&amp;epsiv;</mi> <mrow> <mi>F</mi> <mi>B</mi> <mi>G</mi> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>&amp;epsiv;</mi> <mrow> <mi>F</mi> <mi>B</mi> <mi>G</mi> <mn>2</mn> </mrow> </msub> </mrow> <mi>d</mi> </mfrac> </mrow>
其中,C为曲率值,εFBG1为第一纤芯轴向应力,εFBG2为第二纤芯轴向应力,d为两个刻有Bragg光栅的纤芯距离。
5.根据权利要求1所述的多芯光纤Bragg光栅曲率传感器的解调方法,其特征在于,所述余弦型宽带光源的消光比大于10dB。
CN201710289611.6A 2017-04-27 2017-04-27 一种多芯光纤Bragg光栅曲率传感器的解调方法 Active CN107063317B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710289611.6A CN107063317B (zh) 2017-04-27 2017-04-27 一种多芯光纤Bragg光栅曲率传感器的解调方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710289611.6A CN107063317B (zh) 2017-04-27 2017-04-27 一种多芯光纤Bragg光栅曲率传感器的解调方法

Publications (2)

Publication Number Publication Date
CN107063317A true CN107063317A (zh) 2017-08-18
CN107063317B CN107063317B (zh) 2019-12-31

Family

ID=59605203

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710289611.6A Active CN107063317B (zh) 2017-04-27 2017-04-27 一种多芯光纤Bragg光栅曲率传感器的解调方法

Country Status (1)

Country Link
CN (1) CN107063317B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108225207A (zh) * 2017-12-27 2018-06-29 武汉华星光电半导体显示技术有限公司 柔性显示装置及柔性显示面板的弯曲状态量的获取方法
CN109000687A (zh) * 2018-07-19 2018-12-14 西南交通大学 一种基于偏芯相移光纤光栅的曲率解调装置及其方法
CN113446962A (zh) * 2021-07-13 2021-09-28 西南交通大学 一种基于强耦合多芯光纤的温度不敏感曲率传感器、曲率测量装置和方法
CN113984097A (zh) * 2021-12-27 2022-01-28 之江实验室 用于多芯光纤三维形状传感的片上解调系统及承载设备

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1737676A (zh) * 2005-08-26 2006-02-22 天津大学 光纤Bragg光栅传感解调装置及解调方法
CN101840018A (zh) * 2010-04-09 2010-09-22 哈尔滨工程大学 多芯长周期光纤光栅制造方法及光纤旋转定位装置
CN202692951U (zh) * 2012-07-26 2013-01-23 中铁二院工程集团有限责任公司 用于柔性防护网受力变形监测的光纤光栅弯曲传感器
CN105928469A (zh) * 2016-06-07 2016-09-07 华中科技大学 一种高灵敏可判别弯曲方向的无温度交叉灵敏的光纤曲率传感器
CN205655805U (zh) * 2016-05-04 2016-10-19 中国计量大学 基于包层模式干涉结合布拉格光纤光栅曲率传感器
CN106546187A (zh) * 2016-10-19 2017-03-29 暨南大学 一种长周期光纤光栅曲率传感器及其检测方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1737676A (zh) * 2005-08-26 2006-02-22 天津大学 光纤Bragg光栅传感解调装置及解调方法
CN101840018A (zh) * 2010-04-09 2010-09-22 哈尔滨工程大学 多芯长周期光纤光栅制造方法及光纤旋转定位装置
CN202692951U (zh) * 2012-07-26 2013-01-23 中铁二院工程集团有限责任公司 用于柔性防护网受力变形监测的光纤光栅弯曲传感器
CN205655805U (zh) * 2016-05-04 2016-10-19 中国计量大学 基于包层模式干涉结合布拉格光纤光栅曲率传感器
CN105928469A (zh) * 2016-06-07 2016-09-07 华中科技大学 一种高灵敏可判别弯曲方向的无温度交叉灵敏的光纤曲率传感器
CN106546187A (zh) * 2016-10-19 2017-03-29 暨南大学 一种长周期光纤光栅曲率传感器及其检测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
M.J.GANDER ETC: ""Bend measurement using Bragg gratings in multicore fibre"", 《ELECTRONICS LETTERS》 *
YONG-SEN YU ECT: ""bend sensor using an embedded etched fiber bragg grating"", 《MICROWAVE AND OPPTICAL TECHNOLOGY LETTERS》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108225207A (zh) * 2017-12-27 2018-06-29 武汉华星光电半导体显示技术有限公司 柔性显示装置及柔性显示面板的弯曲状态量的获取方法
WO2019127808A1 (zh) * 2017-12-27 2019-07-04 武汉华星光电半导体显示技术有限公司 柔性显示装置及柔性显示面板的弯曲状态量的获取方法
CN109000687A (zh) * 2018-07-19 2018-12-14 西南交通大学 一种基于偏芯相移光纤光栅的曲率解调装置及其方法
CN109000687B (zh) * 2018-07-19 2021-01-01 西南交通大学 一种基于偏芯相移光纤光栅的曲率解调装置及其方法
CN113446962A (zh) * 2021-07-13 2021-09-28 西南交通大学 一种基于强耦合多芯光纤的温度不敏感曲率传感器、曲率测量装置和方法
CN113984097A (zh) * 2021-12-27 2022-01-28 之江实验室 用于多芯光纤三维形状传感的片上解调系统及承载设备
CN113984097B (zh) * 2021-12-27 2022-03-15 之江实验室 用于多芯光纤三维形状传感的片上解调系统及承载设备

Also Published As

Publication number Publication date
CN107063317B (zh) 2019-12-31

Similar Documents

Publication Publication Date Title
CN107063317A (zh) 一种多芯光纤Bragg光栅曲率传感器的解调方法
CN103604446B (zh) 一种基于单探测器的多通道光纤光栅绝对波长解调系统的解调方法
CN107024236A (zh) F‑p/fbg光纤传感器解调系统
CN102589617B (zh) 一种基于啁啾光纤光栅的全光纤型多参量监测系统
CN105716755B (zh) 一种基于Loyt-Sagnac干涉仪的灵敏度增强型传感器
CN104215368A (zh) 一种f-p腔光纤压力传感装置及其解调方法
CN102169027A (zh) 准分布式光纤温度、应力传感器及检测仪
CN103528609A (zh) 复合干涉型的多参量光纤传感器
CN105973282A (zh) 光纤f-p传感器腔长小波相位提取解调方法
CN104297208A (zh) 基于光子晶体光纤的干涉型光纤传感器
CN107860405A (zh) 一种基于游标效应的光谱解调方法及其解调装置
CN108387251B (zh) 一种光纤光栅解调仪器及方法
CN105241482B (zh) 一种有源光纤光栅传感器波长解调系统及方法
CN103940363A (zh) 基于小波互相关技术的高精度光纤应变低频传感解调方法
CN104390594B (zh) 光纤微结构位移传感器
CN105758434A (zh) 基于线阵InGaAs扫描FBG反射谱的传感解调方法
CN204963883U (zh) 一种利用啁啾光纤光栅实现激光拍频的测量装置
CN105180977A (zh) 一种单光纤迈克尔逊干涉传感器及传感系统
CN109556756B (zh) 基于多波长光纤激光器游标效应的温度传感器
CN111811554A (zh) 基于光腔衰荡大范围高精度光纤光栅传感方法及装置
CN203083927U (zh) 基于单模-细芯-多模-单模结构的光纤折射率传感器
CN104848879B (zh) 基于线性工作匹配光栅法的光纤布拉格光栅传感器信号解调方法
CN103453940A (zh) 基于多模结构的光纤传感器
CN105953919B (zh) 一种全光纤傅里叶光谱分析仪
CN108204827A (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