CN111006753B - 一种相位反馈控制的光纤干涉超低频振动测量装置及方法 - Google Patents

一种相位反馈控制的光纤干涉超低频振动测量装置及方法 Download PDF

Info

Publication number
CN111006753B
CN111006753B CN201911397524.8A CN201911397524A CN111006753B CN 111006753 B CN111006753 B CN 111006753B CN 201911397524 A CN201911397524 A CN 201911397524A CN 111006753 B CN111006753 B CN 111006753B
Authority
CN
China
Prior art keywords
phase
signal
optical fiber
frequency
feedback
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
CN201911397524.8A
Other languages
English (en)
Other versions
CN111006753A (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.)
Zhuhai Renchi Photoelectric Technology Co ltd
Original Assignee
Zhuhai Renchi Photoelectric Technology Co ltd
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 Zhuhai Renchi Photoelectric Technology Co ltd filed Critical Zhuhai Renchi Photoelectric Technology Co ltd
Priority to CN201911397524.8A priority Critical patent/CN111006753B/zh
Publication of CN111006753A publication Critical patent/CN111006753A/zh
Application granted granted Critical
Publication of CN111006753B publication Critical patent/CN111006753B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H9/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
    • G01H9/004Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means using fibre optic sensors

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Instruments For Measurement Of Length By Optical Means (AREA)

Abstract

本公开提供了一种相位反馈控制的光纤干涉超低频振动测量装置及方法,包括信号发生器、激光器、隔离器、光纤耦合器、芯轴式振动探头、反馈控制模块、相位延迟补偿模块和驱动器、相位控制器、法拉第旋镜、光电探测器、相位解调模块、滤波器和处理器;本公开利用反馈控制模块、相位延迟补偿模块、驱动器、相位控制器组成的反馈回路,控制光纤干涉仪参考臂的光纤长度变化,改变了超低频振动测量系统的响应特性,提高了系统对超低频特环境振动的测量灵敏度;在闭环控制系统基础上,带反馈的相位信息经过滤波器实现滤波处理,提高了振动测量精度,改善了工作频带内的平坦度。

Description

一种相位反馈控制的光纤干涉超低频振动测量装置及方法
技术领域
本公开涉及振动测量技术领域,特别涉及一种相位反馈控制的光纤干涉超低频振动测量装置及方法。
背景技术
本部分的陈述仅仅是提供了与本公开相关的背景技术,并不必然构成现有技术。
现代地震仪通过探测地震波,得到地震记录来研究震源、地球内部结构和地震波本身,并实现临震预报。但是,传统的地震波探测系统普遍存在着灵敏度低、动态范围小、漏电、供电困难等问题,限制了地震预报技术的发展。鉴于光纤传感技术耐腐蚀、抗电磁干扰、动态范围大、工作频带宽、灵敏度高、绝缘性好、利于组网及长距离传输等优点,光纤振动仪可以明显解决传统电动式地震仪的应用难题。目前,光纤振动仪主要包括光纤光栅类和光纤干涉类,前者利用光纤光栅测量环境振动,但由于动态范围小,灵敏度低,且无法测量长周期信号,一直未成功应用到自然地震监测中;后者采用非平衡迈克尔逊光纤干涉仪将环境振动转化为干涉仪相位变化。干涉仪输出的信号经过光电转换后,用采集卡转换为数字信号输入计算机,最后利用相位解调装置,得到与被测振动有关的相位信息,可以实现超低频振动信号的测量,代表了新一代地震仪的发展方向。
本公开发明人发现,现有的光纤干涉振动测量装置及方法,虽然相比较于传统电动式甚宽频地震仪具有诸多应用优势,但受激光器相位噪声水平的制约,其加速度最小可测水平在超低频范围内(低于0.1Hz)高于地球低背景噪声模型(NLNM),使系统工作频带下限的截止频率较高,很难测量到超低频的自然地震信号。
发明内容
为了解决现有技术的不足,本公开提供了一种相位反馈控制的光纤干涉超低频振动测量装置及方法,实现了自然地震的无源、高灵敏度和高精度测量。
为了实现上述目的,本公开采用如下技术方案:
本公开第一方面提供了一种相位反馈控制的光纤干涉超低频振动测量装置。
一种相位反馈控制的光纤干涉超低频振动测量装置,包括信号发生器、激光器、隔离器、光纤耦合器、芯轴式振动探头、相位控制器、法拉第旋镜、光电探测器、相位解调模块、滤波器和处理器;
所述信号发生器输出三路信号,分别为调频信号、一倍频信号和二倍频信号,调频信号传输给激光器,所述激光器产生带有波长调制的窄线宽连续光通过隔离器后传输到光纤耦合器的第一端口;
光纤耦合器的第三端口和第四端口分别连接传感臂光纤和参考臂光纤,传感臂光纤分为不敏感部分和敏感部分,不敏感部分光纤作为光信号的传输介质,敏感部分光纤缠绕在芯轴式振动探头上用于感知环境振动;
参考臂光纤分为不敏感部分和敏感部分,不敏感部分光纤仅作为光信号的传输介质,敏感部分光纤缠绕在相位控制器上,实现参考臂光纤轴向长度变化的检测;
光纤耦合器的第三端口和第四端口接收传感臂和参考臂返回的光信号,并在第二端口形成干涉光,输出给光电探测器;
光电探测器将干涉光信号转化为干涉电信号,输出给相位解调模块,相位解调模块利用输入的一倍频信号、二倍频信号和干涉电信号,解调得到能够反映环境振动的相位信号,并通过滤波器滤波后传输给处理器;
处理器根据输入的相位信号以及相位信号相对于环境振动的响应特性,计算得到环境振动的速度、加速度和位移测量数据。
作为可能的一些实现方式,所述芯轴式振动探头缠绕有传感光纤,将环境振动转化为传感光纤的轴向长度变化,所述芯轴式振动探头内置法拉第旋镜,连接在传感光纤的末端,用于使传感臂光信号反射回光纤耦合器。
作为可能的一些实现方式,所述压电陶瓷光纤拉伸器连接有法拉第旋镜,用于使参考臂光信号反射回光纤耦合器。
作为可能的一些实现方式,还包括反馈控制模块、相位延迟补偿模块和驱动器;
光电探测器将干涉光信号转化为干涉电信号并输出给反馈控制模块,反馈控制模块根据环境振动的测量需求输出反馈信号给相位延迟补偿模块;
相位延迟补偿模块用来补偿相位解调模块和反馈控制模块引起的相位延迟,使反馈信号和相位信号同步,输出延迟补偿后的反馈信号给驱动器;
驱动器针对参考臂的相位控制器特性,产生合适的驱动信号给相位控制器,使相位控制器引起的参考臂光纤轴向长度变化与延迟补偿后的反馈信号在工作频带内呈现幅度和相位平坦响应特征。
作为进一步的限定,所述芯轴式振动探头、光纤耦合器、光电探测器、相位解调模块、反馈控制模块、相位延迟补偿模块、驱动器、相位控制器和法拉第旋镜组成以环境振动为输入、以相位信号为输出的闭环控制系统;
所述反馈控制模块、相位延迟补偿模块、驱动器、相位控制器组成上述闭环控制系统的反馈回路。
作为可能的一些实现方式,所述滤波器用于对相位信号进行进一步滤波处理,滤除非工作频带的噪声,改善工作频带内的平坦度;
作为可能的一些实现方式,所述调频信号的角频率为ωc,一倍频信号和二倍频信号相位同步,角频率分别为ωc和2ωc,且与调频信号之间存在固定的相位延迟
Figure BDA0002346699820000041
作为可能的一些实现方式,所述隔离器的输出光送给光纤耦合器的第一端口,用于隔离激光器与后面光纤链路,防止光纤耦合器第一端口的反射光影响激光器的稳定性。
作为可能的一些实现方式,所述相位控制器为压电陶瓷光纤拉伸器,敏感部分光纤放置在压电陶瓷光纤拉伸器内部,所述驱动器为压电陶瓷驱动器。
本公开第二方面提供了一种相位反馈控制的光纤干涉超低频振动测量装置。
一种相位反馈控制的光纤干涉超低频振动测量装置,包括
包括信号发生器、激光器、隔离器、第一光纤耦合器、第二光纤耦合器、盘片式振动探头、相位控制器、光电探测器、相位解调模块、滤波器和处理器;
所述信号发生器输出三路信号,分别为调频信号、一倍频信号和二倍频信号,调频信号传输给激光器,所述激光器产生带有波长调制的窄线宽连续光通过隔离器后传输到第一光纤耦合器的第一端口;
第一光纤耦合器的第三端口和第四端口分别连接传感臂光纤和参考臂光纤,传感臂光纤分为不敏感部分和敏感部分,不敏感部分光纤作为光信号的传输介质,敏感部分光纤缠绕在盘片式振动探头上,感知环境振动,并将光信号传输给第二光纤耦合器;
参考臂光纤分为不敏感部分和敏感部分,不敏感部分光纤仅作为光信号的传输介质,敏感部分光纤缠绕在相位控制器上,实现参考臂光纤轴向长度变化的检测,并将光信号传输给第二光纤耦合器;
第二光纤耦合器的第一端口和第二端口接收传感臂和参考臂返回的光信号,并在第四端口形成干涉光,输出给光电探测器;
光电探测器将干涉光信号转化为干涉电信号,输出给相位解调模块,相位解调模块利用输入的一倍频信号、二倍频信号和干涉电信号,解调得到能够反映环境振动的相位信号,并通过滤波器滤波后传输给处理器;
处理器根据输入的相位信号以及相位信号相对于环境振动的响应特性,计算得到环境振动的速度、加速度和位移测量数据。
作为可能的一些实现方式,还包括反馈控制模块、相位延迟补偿模块和驱动器;
光电探测器将干涉光信号转化为干涉电信号并输出给反馈控制模块,反馈控制模块根据环境振动的测量需求输出反馈信号给相位延迟补偿模块;
相位延迟补偿模块用来补偿相位解调模块和反馈控制模块引起的相位延迟,使反馈信号和相位信号同步,输出延迟补偿后的反馈信号给驱动器;
驱动器针对参考臂的相位控制器特性,产生合适的驱动信号给相位控制器,使相位控制器引起的参考臂光纤轴向长度变化与延迟补偿后的反馈信号在工作频带内呈现幅度和相位平坦响应特征。
作为进一步的限定,所述盘片式振动探头、第一光纤耦合器、第二光纤耦合器、光电探测器、相位解调模块、反馈控制模块、相位延迟补偿模块、驱动器和相位控制器组成以环境振动为输入、以相位信号为输出的闭环控制系统;且反馈控制模块、相位延迟补偿模块、驱动器和相位控制器组成了上述闭环控制系统的反馈回路。
作为可能的一些实现方式,所述滤波器用于对相位信号进行进一步滤波处理,滤除非工作频带的噪声,改善工作频带内的平坦度;
作为可能的一些实现方式,所述调频信号的角频率为ωc,一倍频信号和二倍频信号相位同步,角频率分别为ωc和2ωc,且与调频信号之间存在固定的相位延迟
Figure BDA0002346699820000061
作为可能的一些实现方式,所述隔离器的输出光送给光纤耦合器的第一端口,用于隔离激光器与后面光纤链路,防止光纤耦合器第一端口的反射光影响激光器的稳定性。
作为可能的一些实现方式,所述相位控制器为压电陶瓷光纤拉伸器,敏感部分光纤放置在压电陶瓷光纤拉伸器内部,所述驱动器为压电陶瓷驱动器。
本公开第三方面提供了一种相位反馈控制的光纤干涉超低频振动测量方法。
一种相位反馈控制的光纤干涉超低频振动测量方法,利用本公开第一方面或者第二方面所述的相位反馈控制的光纤干涉超低频振动测量装置,包括以下步骤:
激光器输出带有频率调制的,适合光纤传输的连续光;
连续光经过包括传感臂和参考臂的光纤干涉仪,在光纤耦合器输出端口形成干涉光信号;
干涉光经过相位解调模块,得到能够反映环境振动的相位信号;
相位信息经过反馈回路,使缠绕在参考臂中相位控制器的光纤轴向长度发生变化;
在反馈回路的作用下,组成以环境振动为输入、以相位信号为输出的闭环控制系统;
带反馈的相位信息经过滤波,滤除工作频带外的噪声,改善工作频带内的平坦度;
根据相位信号,以及相位信号相对于环境振动的响应特性,计算得到环境振动的速度、加速度、位移等测量数据。
作为可能的一些实现方式,闭环控制系统输出的相位信息相对于环境振动加速度的传递函数被调整为:
Figure BDA0002346699820000071
其中,
Figure BDA0002346699820000072
表示闭环控制系统输出的相位信号,
Figure BDA0002346699820000073
为相位解调模块输出的反映环境振动的相位信号,
Figure BDA0002346699820000074
表示环境振动的加速度,A1和A2表示直流增益,
Figure BDA0002346699820000075
表示二阶响应模型,Lrb为参考臂光纤长度变化。
本公开第四方面提供了一种光纤干涉地震仪,包括本公开第一方面或者第二方面所述的相位反馈控制的光纤干涉超低频振动测量装置。
与现有技术相比,本公开的有益效果是:
1、本公开利用反馈控制模块、相位延迟补偿模块、驱动器、相位控制器组成的反馈回路,控制光纤干涉仪参考臂的光纤长度变化,改变了超低频振动测量系统的响应特性,提高了系统对超低频特环境振动的测量灵敏度。
2、本公开在闭环控制系统基础上,带反馈的相位信息经过滤波器实现滤波处理,提高了振动测量精度,改善了工作频带内的平坦度。
附图说明
图1为本公开实施例1提供的通用的相位反馈控制的光纤干涉超低频振动测量装置的结构示意图图。
图2为本公开实施例3提供相位反馈控制的光纤干涉超低频振动测量装置的结构示意图。
图3为本公开实施例2或实施例4或实施例6提供的无反馈传递函数结构图。
图4为本公开实施例2或实施例4或实施例6提供的有反馈传递函数结构图。
图5为本公开实施例5提供相位反馈控制的光纤干涉超低频振动测量装置的结构示意图。
具体实施方式
应该指出,以下详细说明都是示例性的,旨在对本公开提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本公开所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本公开的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。
实施例1:
如图1所示,本公开实施例1提供了一种通用的相位反馈控制的光纤干涉超低频振动测量装置具体包括:信号发生器、激光器、隔离器、光纤耦合器、振动探头、相位控制器、法拉第旋镜、光电探测器、相位解调模块、反馈控制模块、相位延迟补偿模块、驱动器、滤波器和处理器。
信号发生器输出3路信号,包括调频信号,1倍频信号和2倍频信号。调频信号的角频率为ωc,输出给激光器;1倍频信号和2倍频信号相位同步,角频率分别为ωc和2ωc,且与调频信号之间存在固定的相位延迟
Figure BDA0002346699820000091
输出到相位解调模块。
激光器应具有窄线宽特性,在信号发生器调频信号的激励下,产生带有波长调制的窄线宽连续光,输出给隔离器。
隔离器的输出光送给光纤耦合器的1端口,起到隔离激光器与后面光纤链路的作用,防止光纤耦合器1端口的反射光影响激光器的稳定性。
光纤耦合器的3、4端口分别连接传感臂光纤和参考臂光纤。传感臂光纤分为不敏感部分和敏感部分。不敏感部分光纤仅作为光信号的传输介质;敏感部分光纤缠绕在振动探头上,感知环境振动。
振动探头缠绕有传感光纤,将环境振动转化为传感光纤的轴向长度变化。振动探头内置法拉第旋镜,连接在传感光纤的末端,使传感臂光信号反射回光纤耦合器。
参考臂光纤分为不敏感部分和敏感部分。不敏感部分光纤仅作为光信号的传输介质;敏感部分光纤缠绕在相位控制器上,实现参考臂光纤轴向长度变化的检测。
相位控制器连接有法拉第旋镜,使参考臂光信号反射回光纤耦合器。
光纤耦合器的3、4端口接收传感臂和参考臂返回的光信号,并在2端口形成干涉光,输出给光电探测器。
光电探测器将干涉光信号转化为干涉电信号,输出给相位解调模块。
相位解调模块利用输入的1倍频信号、2倍频信号和干涉电信号,使用经典的微分交叉相乘(DCM)或者反正切解调方法,解调得到能够反映环境振动的相位信号,并输出给滤波器和反馈控制模块。
反馈控制模块根据环境振动的测量需求,拥有合适的二阶或者更高阶输入输出响应特性,输出反馈信号给相位延迟补偿模块。
相位延迟补偿模块用来补偿相位解调模块和反馈控制模块引起的相位延迟,使反馈信号和相位信号同步,输出延迟补偿后的反馈信号给驱动器。
驱动器针对参考臂的相位控制器特性,产生合适的驱动信号给相位控制器,使相位控制器引起的参考臂光纤轴向长度变化与延迟补偿后的反馈信号在工作频带内呈现幅度和相位平坦响应特征。
其中,振动探头、光纤耦合器、光电探测器、相位解调模块、反馈控制模块、相位延迟补偿模块、驱动器、相位控制器、法拉第旋镜组成以环境振动为输入、以相位信号为输出的闭环控制系统。反馈控制模块、相位延迟补偿模块、驱动器、相位控制器组成了该闭环控制系统的反馈回路,用于改变闭环控制系统的响应特性,提高系统对特定频率范围内环境振动的测量灵敏度。
滤波器对相位信号进行进一步滤波处理,滤除非工作频带的噪声,改善工作频带内的平坦度,输出给处理器。
处理器根据输入的相位信号,以及相位信号相对于环境振动的响应特性,计算得到环境振动的速度、加速度、位移等测量数据,并进行记录。
实施例2:
本公开实施例2提供了一种相位反馈控制的光纤干涉超低频振动测量方法,包括以下步骤:
步骤一:激光器输出带有频率调制的,适合光纤传输的波长为1550nm或1330nm的连续光,该连续光的电场表达式为:
Figure BDA0002346699820000111
Figure BDA0002346699820000112
其中,EL(t)表示激光器输出连续光的电场强度,EC表示电场强度的幅度值,vL(t)带有频率调制的光频,
Figure BDA0002346699820000113
表示随时间随机变化的激光器相位噪声,vC表示光频的直流分量,Kv表示光频调制系数,
Figure BDA0002346699820000114
表示光频调制的初始相位,t表示时间。
步骤二:连续光经过包括传感臂和参考臂的光纤干涉仪,在光纤耦合器输出端口形成干涉光信号:
Figure BDA0002346699820000115
其中,Iinter(t)表示干涉光,Kct和Kcc分别表示光纤耦合器的直通系数和耦合系数,KM表示振动探头中和参考臂的法拉第旋镜的反射系数,n表示光纤折射率,LA和LB分别表示传感臂和参考臂光纤的初始长度,Lra(t)表示环境振动通过振动探头引起的传感臂光纤轴向长度变化,t1表示光波在传感臂光纤来回程的时间延迟,t2表示光波在参考臂光纤来回程的时间延迟。
步骤三:干涉光经过相位解调模块,得到能够反映环境振动的相位信号,此时相位信号相对于环境振动加速度的传递函数为:
Figure BDA0002346699820000121
其中,
Figure BDA0002346699820000122
表示相位信号,
Figure BDA0002346699820000123
表示环境振动的加速度,A1和A2表示直流增益,
Figure BDA0002346699820000124
表示二阶响应模型,其中ζ是阻尼系数,ωn是固有频率,s=jω是拉普拉斯算子,ω表示角频率。A1、ζ、ωn由振动探头的物理参数,以及缠绕在振动探头上的光纤长度决定,均为常量。
Figure BDA0002346699820000125
由光频的直流分量vC、光纤折射率n和光速c决定,也为常量。
需要进一步指出的是,相位信号相对于环境振动加速度的传递函数,可以分解为传感臂光纤长度变化相对于环境振动加速度传递函数(公式4-1),以及相位信号相对于传感臂光纤长度变化传递函数(公式4-2)的叠加。
Figure BDA0002346699820000126
Figure BDA0002346699820000127
步骤四:相位信息经过反馈回路,使缠绕在参考臂中相位控制器的光纤轴向长度发生变化。参考臂光纤长度变化相对于相位信息的传递函数为:
Figure BDA0002346699820000128
其中,m和n分别表示传递函数分母和分子项的最高阶次,αm、αm-1、αm-2、…、α1、α0表示分母项中对应阶次的系数,βn、βn-1、βn-2、β1、β0表示分子项中对应阶次的系数。该传递函数由反馈回路的反馈控制模块、相位延迟补偿模块、驱动器、相位控制器的参数决定。
步骤五:在反馈回路的作用下,振动探头、光纤耦合器、光电探测器、相位解调模块、反馈控制模块、相位延迟补偿模块、驱动器、相位控制器、法拉第旋镜组成以环境振动为输入、以相位信号为输出的闭环控制系统。该闭环控制系统输出的相位信息相对于环境振动加速度的传递函数被调整为:
Figure BDA0002346699820000131
通过选择合适的反馈回路参数,可以提高工作频带内相位信息相对于环境振动加速度的灵敏度。
步骤六:带反馈的相位信息经过滤波器,滤除工作频带外的噪声,改善工作频带内的平坦度。然后根据相位信号,以及相位信号相对于环境振动的响应特性,可以计算得到环境振动的速度、加速度、位移等测量数据。
实施例3:
如图3所示,本公开实施例3提供了一种相位反馈控制的光纤干涉超低频振动测量装置,将实施例1中的相位控制器替换为压电陶瓷光纤拉伸器,敏感部分光纤放置在压电陶瓷光纤拉伸器内部,实现参考臂光纤轴向长度变化的检测;将驱动器替换为压电陶瓷驱动器,其他结构设置与实施例1相同,这里不再赘述。
实施例4:
本公开实施例4提供了一种相位反馈控制的光纤干涉超低频振动测量方法,将实施例2中的相位控制器替换为压电陶瓷光纤拉伸器,敏感部分光纤放置在压电陶瓷光纤拉伸器内部,实现参考臂光纤轴向长度变化的检测;将驱动器替换为压电陶瓷驱动器,其他方法步骤与实施例2相同,这里不再赘述。
实施例5:
如图5所示,本公开实施例5提供了一种相位反馈控制的光纤干涉超低频振动测量装置,具体包括:信号发生器、激光器、隔离器、第一光纤耦合器、第二光纤耦合器、盘片式振动探头、压电陶瓷光纤拉伸器、光电探测器、相位解调模块、反馈控制模块、相位延迟补偿模块、压电陶瓷驱动器、滤波器和处理器。
信号发生器输出三路信号,包括调频信号,一倍频信号和二倍频信号。调频信号的角频率为ωc,输出给激光器;一倍频信号和二倍频信号相位同步,角频率分别为ωc和2ωc,且与调频信号之间存在固定的相位延迟
Figure BDA0002346699820000141
输出到相位解调模块。
激光器应具有窄线宽特性,在信号发生器调频信号的激励下,产生带有波长调制的窄线宽连续光,输出给隔离器。
隔离器的输出光送给第一光纤耦合器的1端口,起到隔离激光器与后面光纤链路的作用,防止第一光纤耦合器1端口的反射光影响激光器的稳定性。
第一光纤耦合器的3、4端口分别连接传感臂光纤和参考臂光纤。传感臂光纤分为不敏感部分和敏感部分。不敏感部分光纤仅作为光信号的传输介质;敏感部分光纤缠绕在盘片式振动探头上,感知环境振动。
盘片式振动探头缠绕有传感光纤,将环境振动转化为传感光纤的轴向长度变化,传感臂光信号传输到第二光纤耦合器的1端口。
参考臂光纤分为不敏感部分和敏感部分。不敏感部分光纤仅作为光信号的传输介质;敏感部分光纤放置在压电陶瓷光纤拉伸器内部,实现参考臂光纤轴向长度变化的检测,参考臂光信号传输到第二光纤耦合器的2端口。
第二光纤耦合器的4端口形成干涉光,输出给光电探测器。
光电探测器将干涉光信号转化为干涉电信号,输出给相位解调模块。
相位解调模块利用输入的一倍频信号、二倍频信号和干涉电信号,使用经典的微分交叉相乘(DCM)或者反正切解调方法,解调得到能够反映环境振动的相位信号,并输出给滤波器和反馈控制模块。
反馈控制模块根据环境振动的测量需求,拥有合适的二阶或者更高阶输入输出响应特性,输出反馈信号给相位延迟补偿模块。
相位延迟补偿模块用来补偿相位解调模块和反馈控制模块引起的相位延迟,使反馈信号和相位信号同步,输出延迟补偿后的反馈信号给压电陶瓷驱动器。
压电陶瓷驱动器针对参考臂的相位控制器特性,产生合适的驱动信号给压电陶瓷光纤拉伸器,使压电陶瓷光纤拉伸器引起的参考臂光纤轴向长度变化与延迟补偿后的反馈信号在工作频带内呈现幅度和相位平坦响应特征。
其中,盘片式振动探头、第一光纤耦合器、第二光纤耦合器、光电探测器、相位解调模块、反馈控制模块、相位延迟补偿模块、压电陶瓷驱动器和压电陶瓷光纤拉伸器组成以环境振动为输入、以相位信号为输出的闭环控制系统。
反馈控制模块、相位延迟补偿模块、压电陶瓷驱动器、压电陶瓷光纤拉伸器组成了该闭环控制系统的反馈回路,用于改变闭环控制系统的响应特性,提高系统对特定频率范围内环境振动的测量灵敏度。
滤波器对相位信号进行进一步滤波处理,滤除非工作频带的噪声,改善工作频带内的平坦度,输出给处理器。
处理器根据输入的相位信号,以及相位信号相对于环境振动的响应特性,计算得到环境振动的速度、加速度、位移等测量数据,并进行记录。
实施例6:
本公开实施例6提供了一种相位反馈控制的光纤干涉超低频振动测量方法,包括以下步骤:
步骤一:激光器输出带有频率调制的,适合光纤传输的波长为1550nm或1330nm的连续光。该连续光的电场表达式为:
Figure BDA0002346699820000161
Figure BDA0002346699820000162
其中,EL(t)表示激光器输出连续光的电场强度,EC表示电场强度的幅度值,vL(t)带有频率调制的光频,
Figure BDA0002346699820000163
表示随时间随机变化的激光器相位噪声,vC表示光频的直流分量,Kv表示光频调制系数,
Figure BDA0002346699820000164
表示光频调制的初始相位,t表示时间。
步骤二:连续光经过包括传感臂和参考臂的光纤干涉仪,在第二光纤耦合器输出端口形成干涉光信号:
Figure BDA0002346699820000165
其中,Iinter(t)表示干涉光,Kct和Kcc分别表示第二光纤耦合器的直通系数和耦合系数,KM表示振动探头中和参考臂的法拉第旋镜的反射系数,n表示光纤折射率,LA和LB分别表示传感臂和参考臂光纤的初始长度,Lra(t)表示环境振动通过振动探头引起的传感臂光纤轴向长度变化,t1表示光波在传感臂光纤传输与在探头和第二光纤耦合器之间传输的时间延迟,t2表示光波在参考臂光纤传输与在探头和第二光纤耦合器之间传输的时间延迟。
步骤三:干涉光经过相位解调模块,得到能够反映环境振动的相位信号,此时相位信号相对于环境振动加速度的传递函数为:
Figure BDA0002346699820000171
其中,
Figure BDA0002346699820000172
表示相位信号,
Figure BDA0002346699820000173
表示环境振动的加速度,A1和A2表示直流增益,
Figure BDA0002346699820000174
表示二阶响应模型,其中ζ是阻尼系数,ωn是固有频率,s=jω是拉普拉斯算子,ω表示角频率。A1、ζ、ωn由振动探头的物理参数,以及缠绕在振动探头上的光纤长度决定,均为常量。
Figure BDA0002346699820000175
由光频的直流分量νc、光纤折射率n和光速c决定,也为常量。
需要进一步指出的是,相位信号相对于环境振动加速度的传递函数,可以分解为传感臂光纤长度变化相对于环境振动加速度传递函数(公式10-1),以及相位信号相对于传感臂光纤长度变化传递函数(公式10-2)的叠加。该传递函数结构图如图3所示。
Figure BDA0002346699820000176
Figure BDA0002346699820000177
步骤四:相位信息经过反馈回路,使缠绕在参考臂中相位控制器的光纤轴向长度发生变化。
参考臂光纤长度变化相对于相位信息的传递函数为:
Figure BDA0002346699820000181
其中,m和n分别表示传递函数分母和分子项的最高阶次,αm、αm-1、αm-2、…、α1、α0表示分母项中对应阶次的系数,βn、βn-1、βn-2、β1、β0表示分子项中对应阶次的系数。该传递函数由反馈回路的反馈控制模块、相位延迟补偿模块、压电陶瓷驱动器、压电陶瓷光纤拉伸器的参数决定。
步骤五:在反馈回路的作用下,芯轴式振动探头、光纤耦合器、光电探测器、相位解调模块、反馈控制模块、相位延迟补偿模块、压电陶瓷驱动器、压电陶瓷光纤拉伸器、法拉第旋镜组成以环境振动为输入、以相位信号为输出的闭环控制系统。该闭环控制系统输出的相位信息相对于环境振动加速度的传递函数被调整为:
Figure BDA0002346699820000182
通过选择合适的反馈回路参数,可以提高工作频带内相位信息相对于环境振动加速度的灵敏度。该传递函数结构图如图4所示。
步骤六:带反馈的相位信息经过滤波器,滤除工作频带外的噪声,改善工作频带内的平坦度。然后根据相位信号,以及相位信号相对于环境振动的响应特性,可以计算得到环境振动的速度、加速度、位移等测量数据。
本公开利用反馈控制模块、相位延迟补偿模块、驱动器、相位控制器组成的反馈回路,控制光纤干涉仪参考臂的光纤长度变化,改变了超低频振动测量系统的响应特性,提高了系统对超低频特环境振动的测量灵敏度。在闭环控制系统基础上,带反馈的相位信息经过滤波器实现滤波处理,提高了振动测量精度,改善了工作频带内的平坦度。
实施例7:
本公开实施例7提供了一种光纤干涉地震仪,包括本公开实施例1或者实施例3或者实施例5所述的相位反馈控制的光纤干涉超低频振动测量装置。
以上所述仅为本公开的优选实施例而已,并不用于限制本公开,对于本领域的技术人员来说,本公开可以有各种更改和变化。凡在本公开的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本公开的保护范围之内。

Claims (16)

1.一种相位反馈控制的光纤干涉超低频振动测量装置,其特征在于,包括信号发生器、激光器、隔离器、光纤耦合器、芯轴式振动探头、相位控制器、光电探测器、相位解调模块、滤波器和处理器;
所述信号发生器输出三路信号,分别为调频信号、一倍频信号和二倍频信号,调频信号传输给激光器,所述激光器产生的连续光通过隔离器后传输到光纤耦合器的第一端口;
光纤耦合器的第三端口和第四端口分别连接传感臂光纤和参考臂光纤,传感臂光纤分为不敏感部分和敏感部分,不敏感部分光纤作为光信号的传输介质,敏感部分光纤缠绕在芯轴式振动探头上用于感知环境振动;参考臂光纤分为不敏感部分和敏感部分,不敏感部分光纤仅作为光信号的传输介质,敏感部分光纤缠绕在相位控制器上,实现参考臂光纤轴向长度变化的检测;
光纤耦合器的第三端口和第四端口接收传感臂和参考臂返回的光信号,并在第二端口形成干涉光,输出给光电探测器;光电探测器将干涉光信号转化为干涉电信号,输出给相位解调模块,相位解调模块利用输入的一倍频信号、二倍频信号和干涉电信号,解调得到能够反映环境振动的相位信号,并通过滤波器滤波后传输给处理器;
处理器根据输入的相位信号以及相位信号相对于环境振动的响应特性,计算得到环境振动的速度、加速度和位移测量数据;
还包括反馈控制模块、相位延迟补偿模块和驱动器;
光电探测器将干涉光信号转化为干涉电信号并输出给反馈控制模块,反馈控制模块根据环境振动的测量需求输出反馈信号给相位延迟补偿模块;
相位延迟补偿模块用来补偿相位解调模块和反馈控制模块引起的相位延迟,使反馈信号和相位信号同步,输出延迟补偿后的反馈信号给驱动器;
驱动器针对参考臂的相位控制器特性,产生合适的驱动信号给相位控制器,使相位控制器引起的参考臂光纤轴向长度变化与延迟补偿后的反馈信号在工作频带内呈现幅度和相位平坦响应特征。
2.如权利要求1所述的相位反馈控制的光纤干涉超低频振动测量装置,其特征在于,所述芯轴式振动探头、光纤耦合器、光电探测器、相位解调模块、反馈控制模块、相位延迟补偿模块、驱动器、相位控制器和法拉第旋镜组成以环境振动为输入、以相位信号为输出的闭环控制系统;
所述反馈控制模块、相位延迟补偿模块、驱动器、相位控制器组成上述闭环控制系统的反馈回路。
3.如权利要求1所述的相位反馈控制的光纤干涉超低频振动测量装置,其特征在于,所述滤波器用于对相位信号进行进一步滤波处理,滤除非工作频带的噪声,改善工作频带内的平坦度。
4.如权利要求1所述的相位反馈控制的光纤干涉超低频振动测量装置,其特征在于,所述调频信号的角频率为ωc,一倍频信号和二倍频信号相位同步,角频率分别为ωc和2ωc,且与调频信号之间存在固定的相位延迟
Figure FDA0003316782070000021
5.如权利要求1所述的相位反馈控制的光纤干涉超低频振动测量装置,其特征在于,所述隔离器的输出光送给光纤耦合器的第一端口,用于隔离激光器与后面光纤链路,防止光纤耦合器第一端口的反射光影响激光器的稳定性。
6.如权利要求1所述的相位反馈控制的光纤干涉超低频振动测量装置,其特征在于,所述相位控制器为压电陶瓷光纤拉伸器,敏感部分光纤放置在压电陶瓷光纤拉伸器内部,所述驱动器为压电陶瓷驱动器。
7.如权利要求1所述的相位反馈控制的光纤干涉超低频振动测量装置,其特征在于,所述芯轴式振动探头缠绕有传感光纤,将环境振动转化为传感光纤的轴向长度变化,所述芯轴式振动探头内置法拉第旋镜,连接在传感光纤的末端,用于使传感臂光信号反射回光纤耦合器。
8.如权利要求6所述的相位反馈控制的光纤干涉超低频振动测量装置,其特征在于,所述压电陶瓷光纤拉伸器连接有法拉第旋镜,用于使参考臂光信号反射回光纤耦合器。
9.一种相位反馈控制的光纤干涉超低频振动测量装置,其特征在于,包括信号发生器、激光器、隔离器、第一光纤耦合器、第二光纤耦合器、盘片式振动探头、相位控制器、光电探测器、相位解调模块、滤波器和处理器;
所述信号发生器输出三路信号,分别为调频信号、一倍频信号和二倍频信号,调频信号传输给激光器,所述激光器产生带有波长调制的窄线宽连续光通过隔离器后传输到第一光纤耦合器的第一端口;
第一光纤耦合器的第三端口和第四端口分别连接传感臂光纤和参考臂光纤,传感臂光纤分为不敏感部分和敏感部分,不敏感部分光纤作为光信号的传输介质,敏感部分光纤缠绕在盘片式振动探头上,感知环境振动,并将光信号传输给第二光纤耦合器;
参考臂光纤分为不敏感部分和敏感部分,不敏感部分光纤仅作为光信号的传输介质,敏感部分光纤缠绕在相位控制器上,实现参考臂光纤轴向长度变化的检测,并将光信号传输给第二光纤耦合器;
第二光纤耦合器的第一端口和第二端口接收传感臂和参考臂返回的光信号,并在第四端口形成干涉光,输出给光电探测器;
光电探测器将干涉光信号转化为干涉电信号,输出给相位解调模块,相位解调模块利用输入的一倍频信号、二倍频信号和干涉电信号,解调得到能够反映环境振动的相位信号,并通过滤波器滤波后传输给处理器;
处理器根据输入的相位信号以及相位信号相对于环境振动的响应特性,计算得到环境振动的速度、加速度和位移测量数据;
还包括反馈控制模块、相位延迟补偿模块和驱动器;
光电探测器将干涉光信号转化为干涉电信号并输出给反馈控制模块,反馈控制模块根据环境振动的测量需求输出反馈信号给相位延迟补偿模块;
相位延迟补偿模块用来补偿相位解调模块和反馈控制模块引起的相位延迟,使反馈信号和相位信号同步,输出延迟补偿后的反馈信号给驱动器;
驱动器针对参考臂的相位控制器特性,产生合适的驱动信号给相位控制器,使相位控制器引起的参考臂光纤轴向长度变化与延迟补偿后的反馈信号在工作频带内呈现幅度和相位平坦响应特征。
10.如权利要求9所述的相位反馈控制的光纤干涉超低频振动测量装置,其特征在于,所述滤波器用于对相位信号进行进一步滤波处理,滤除非工作频带的噪声,改善工作频带内的平坦度。
11.如权利要求9所述的相位反馈控制的光纤干涉超低频振动测量装置,其特征在于,所述调频信号的角频率为ωc,一倍频信号和二倍频信号相位同步,角频率分别为ωc和2ωc,且与调频信号之间存在固定的相位延迟
Figure FDA0003316782070000041
12.如权利要求9所述的相位反馈控制的光纤干涉超低频振动测量装置,其特征在于,所述隔离器的输出光送给光纤耦合器的第一端口,用于隔离激光器与后面光纤链路,防止光纤耦合器第一端口的反射光影响激光器的稳定性。
13.如权利要求9所述的相位反馈控制的光纤干涉超低频振动测量装置,其特征在于,所述相位控制器为压电陶瓷光纤拉伸器,敏感部分光纤放置在压电陶瓷光纤拉伸器内部,所述驱动器为压电陶瓷驱动器。
14.一种相位反馈控制的光纤干涉超低频振动测量方法,其特征在于,利用权利要求1-13任一项所述的相位反馈控制的光纤干涉超低频振动测量装置,包括以下步骤:
激光器输出带有频率调制的,适合光纤传输的连续光;
连续光经过包括传感臂和参考臂的光纤干涉仪,在光纤耦合器输出端口形成干涉光信号;
干涉光经过相位解调模块,得到能够反映环境振动的相位信号;
相位信息经过反馈回路,使缠绕在参考臂中相位控制器的光纤轴向长度发生变化;
在反馈回路的作用下,组成以环境振动为输入、以相位信号为输出的闭环控制系统;
带反馈的相位信息经过滤波,滤除工作频带外的噪声,改善工作频带内的平坦度;
根据相位信号,以及相位信号相对于环境振动的响应特性,计算得到环境振动的速度、加速度、位移测量数据。
15.如权利要求14所述的相位反馈控制的光纤干涉超低频振动测量方法,其特征在于,闭环控制系统输出的相位信息相对于环境振动加速度的传递函数被调整为:
Figure FDA0003316782070000051
其中,
Figure FDA0003316782070000061
表示闭环控制系统输出的相位信号,
Figure FDA0003316782070000062
为相位解调模块输出的反映环境振动的相位信号,
Figure FDA0003316782070000063
表示环境振动的加速度,A1和A2表示直流增益,
Figure FDA0003316782070000064
表示二阶响应模型,Lrb为参考臂光纤长度变化。
16.一种光纤干涉地震仪,其特征在于,包括权利要求1-8任一项所述的相位反馈控制的光纤干涉超低频振动测量装置;
或者,包括权利要求9-13任一项所述的相位反馈控制的光纤干涉超低频振动测量装置。
CN201911397524.8A 2019-12-30 2019-12-30 一种相位反馈控制的光纤干涉超低频振动测量装置及方法 Active CN111006753B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911397524.8A CN111006753B (zh) 2019-12-30 2019-12-30 一种相位反馈控制的光纤干涉超低频振动测量装置及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911397524.8A CN111006753B (zh) 2019-12-30 2019-12-30 一种相位反馈控制的光纤干涉超低频振动测量装置及方法

Publications (2)

Publication Number Publication Date
CN111006753A CN111006753A (zh) 2020-04-14
CN111006753B true CN111006753B (zh) 2021-12-03

Family

ID=70119469

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911397524.8A Active CN111006753B (zh) 2019-12-30 2019-12-30 一种相位反馈控制的光纤干涉超低频振动测量装置及方法

Country Status (1)

Country Link
CN (1) CN111006753B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111678583B (zh) * 2020-06-17 2022-02-18 珠海任驰光电科技有限公司 一种光源噪声改善的光纤振动测量装置及方法
CN112097812B (zh) * 2020-09-15 2022-04-19 安徽大学 通过均衡滤波扩展干涉型光纤传感系统工作带宽的方法
CN115224581B (zh) * 2021-04-15 2024-04-05 北京科益虹源光电技术有限公司 激光器的能量补偿方法及装置
CN114384798B (zh) * 2021-12-08 2023-08-08 北京理工大学 减振控制装置及应用其的减振系统
CN117192313B (zh) * 2023-11-08 2024-02-27 国网天津市电力公司电力科学研究院 检测气体绝缘开关设备局部放电的光纤传感系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4486657A (en) * 1982-05-27 1984-12-04 The United States Of America As Represented By The Secretary Of The Navy Phase-lock fiber optic interferometer
CN105092015A (zh) * 2015-05-18 2015-11-25 清华大学 一种非接触式光纤振动传感系统及方法
WO2017041174A1 (en) * 2015-09-08 2017-03-16 Institut National De La Recherche Scientifique System and method for phase-readout and active stabilization of optical interferometers

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6522797B1 (en) * 1998-09-01 2003-02-18 Input/Output, Inc. Seismic optical acoustic recursive sensor system
US6388739B1 (en) * 2001-01-18 2002-05-14 The Boeing Company Self-referencing microdoppler ladar receiver and associated detection method
CN100439859C (zh) * 2007-02-09 2008-12-03 北京交通大学 利用光纤光栅的光纤干涉型在线微位移测量系统
CN100554870C (zh) * 2007-09-25 2009-10-28 北京交通大学 利用光纤光栅及波分复用技术的台阶高度在线测量系统
CN101126629A (zh) * 2007-09-25 2008-02-20 北京交通大学 利用光纤光栅的合成波干涉台阶高度在线测量系统
CN101825432B (zh) * 2010-04-01 2012-04-25 北京交通大学 双波长光纤干涉大量程高分辨率位移测量系统
CN202267808U (zh) * 2011-06-20 2012-06-06 深圳职业技术学院 干涉型光纤传感器数字解调装置
CN102564564A (zh) * 2012-03-09 2012-07-11 大连理工大学 一种基于非平衡迈克尔逊光纤干涉仪的非接触微振动测量系统
CN103208161B (zh) * 2013-03-19 2016-02-03 石家庄供电公司 有源探测型光纤光栅电缆隧道安防监控系统
CN104697609B (zh) * 2015-03-24 2019-02-15 吉林大学 光纤干涉水位传感器
CN105300290B (zh) * 2015-09-28 2017-11-14 北京交通大学 一种基于波数分辨的低相干干涉绝对距离测量系统
CN105300507B (zh) * 2015-10-26 2018-08-28 南京航空航天大学 光纤振动传感器及其m-z传感臂光路结构
CN107179097B (zh) * 2016-03-10 2019-08-09 南京理工大学 基于移相信号调制的光纤干涉型传感调制解调方法及装置
CN106482822B (zh) * 2016-12-07 2023-05-09 吉林大学 一种基于双重同源外差相干检测的相位敏感光时域反射系统
CN107703369A (zh) * 2017-09-29 2018-02-16 北京计算机技术及应用研究所 一种基于磁致伸缩的光纤微弱磁场检测系统
CN108180978A (zh) * 2018-01-30 2018-06-19 广州天长通信技术有限公司 一种结合pgc技术和φ-otdr技术探测光纤振动的方法及装置
CN108663138B (zh) * 2018-05-16 2020-12-01 湖北三江航天万峰科技发展有限公司 一种分布式光纤温度及振动的传感系统及方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4486657A (en) * 1982-05-27 1984-12-04 The United States Of America As Represented By The Secretary Of The Navy Phase-lock fiber optic interferometer
CN105092015A (zh) * 2015-05-18 2015-11-25 清华大学 一种非接触式光纤振动传感系统及方法
WO2017041174A1 (en) * 2015-09-08 2017-03-16 Institut National De La Recherche Scientifique System and method for phase-readout and active stabilization of optical interferometers

Also Published As

Publication number Publication date
CN111006753A (zh) 2020-04-14

Similar Documents

Publication Publication Date Title
CN111006753B (zh) 一种相位反馈控制的光纤干涉超低频振动测量装置及方法
Li et al. Ultra-high sensitive quasi-distributed acoustic sensor based on coherent OTDR and cylindrical transducer
Wang et al. Distributed OTDR-interferometric sensing network with identical ultra-weak fiber Bragg gratings
CN106680535B (zh) 基于光纤布喇格光栅反射谱特性实现激光拍频的差动型光学加速度计
Yang et al. Wideband fiber-optic Fabry-Perot acoustic sensing scheme using high-speed absolute cavity length demodulation
Pan et al. Distributed optical fiber vibration sensor based on Sagnac interference in conjunction with OTDR
Gangopadhyay et al. Vibration: history and measurement with an extrinsic Fabry–Perot sensor with solid-state laser interferometry
CN109724685A (zh) 基于Fizeau干涉的光纤光栅水声传感阵列解调方法及系统
CN104792402A (zh) 一种基于光纤萨格纳克干涉仪的声波传感测量装置
CN107389978A (zh) 一种弱反射布拉格光栅加速度计及其传感方法
CN108344880B (zh) 一种不等臂长迈克尔逊光纤加速度计及其传感方法
CN111308125B (zh) 一种基于光纤Sagnac干涉仪的加加速度探测方法及加加速度计
CN105092015A (zh) 一种非接触式光纤振动传感系统及方法
Duo et al. Self-suppression of common-mode noises of the different fiber optic interferometric accelerometers
Ouyang et al. Nano-g accelerometer with differential Fabry–Pérot interferometer for low-frequency noise suppression
Liu et al. Pico-strain resolution fiber-optic sensor with white-light interferometry
Li et al. Fiber optic jerk sensor
CN103076082A (zh) 基于单模-多模-单模光纤模间干涉振动、应力传感装置
Hong et al. Location of a wideband perturbation using a fiber Fox–Smith interferometer
Huang et al. Ultra-high-resolution phase demodulation based miniature fiber-optic accelerometer at low and medium frequencies
CN104048617B (zh) 对偏振态变化不敏感的高精度光纤光栅传感解调方法
Shang et al. Discussion on the sensitivity of optical cables based on distributed acoustic sensing
Liang et al. Novel fiber Bragg grating sensing method based on the sidelobe modulation for ultrasound detection
Sui et al. Broadband acoustic vibration sensor based on cladding-mode resonance of double-cladding fiber
CN108444464B (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