CN113419341A - 一种光纤干涉仪光场的复数化重构方法及装置 - Google Patents
一种光纤干涉仪光场的复数化重构方法及装置 Download PDFInfo
- 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
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)
I(t)表示干涉仪输出的光强度,A表示直流量,B表示交流量幅度,Re表示取复数的实部运算,表示干涉仪的信号项,为干涉仪的噪声。当涉及到干涉光场的加减乘除运算时,直接采用(公式1)计算时会出现多余频率项。对于两个干涉信号I1(t)和I2(t)
很容易可以推出,进行加法、减法和除法,都会是多个频率项的耦合结果。一种简单的处理方法是,若能够获取干涉光场的复数表达形式,即
将两个干涉光场的复数表达形式直接进行乘法和除法运算,即可获得和再取实部运算即可得到和可见干涉光场的复数表达形式在很多应用中是简化算法、提高计算效率的关键,在光纤干涉仪的复杂信号处理中往往起到化繁为简甚至成为确保算法正确性的基础。
由于干涉光场的复杂性,信号往往加载在余弦函数的相位项中,直接通过数学变换由余弦项获得正弦项会带来频谱成分的移相,对干涉光场的复数化重构需要根据干涉光场的信号特征来进行特殊设计。
发明内容
本发明的目的是提供一种光纤干涉仪光场的复数化重构方法及装置,在实际中,光纤干涉仪的干涉结果在采用光电探测器转换后获取的为实数表达结果,本发明通过将根据实际检测的干涉光场实数表达形式重构为复数表达形式,为光纤干涉仪传感***的串扰抑制、偏振抑制等提供可快速处理的信号基础。
为了实现上述技术目的,本发明采用的技术方案为:
一种光纤干涉仪光场的复数化重构方法,包括以下步骤:
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,即:
S3.1.2将IC和IS分别进行低通滤波,保留小于fm的频谱成分,获得表达式IC_low和IS_low:
S3.1.3生成相位项复光场表达形式Icom:
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)分别减去直流项,获取三个交流项,即:
S3.2.4合成含有干涉仪相位项的干涉结果复数表达形式Icom:
优选的,对于含有外差调制信号的问询脉冲信号,通过光学移频方式产生。
优选的,对于含有外差调制信号的问询脉冲信号,通过采用三角波调制光源频率产生。
优选的,对于含有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),可表达为
当问询脉冲信号为3×3调制信号时,采集三组干涉结果,分别记为I1(t)、I2(t)和I3(t),可表达为
I1(t)、I2(t)和I3(t)表示经过3×3调制干涉仪输出的光强度,A表示直流量,B表示交流量幅度,Re表示取复数的实部运算,表示干涉仪的信号项,为干涉仪的噪声。是干涉仪本身的相位,I1(t)、I2(t)和I3(t)之间由于有3×3调制分别叠加了0和相位。
S3:干涉光场复数化重构:
S3.1针对S2采集的一组干涉结果IH(t),采取以下步骤:
S3.1.1由于外差调制项2πfmt为人为引入项,将IH(t)分别与cos(2πfmt)、sin(2πfmt)混频,获得两个表达式IC和IS,其中:
S3.2.2将IC和IS分别进行低通滤波,保留小于fm的频谱成分,获得:
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)分别减去直流项,获取三个交流项,即:
S3.2.4合成干涉结果复数表达形式Icom
可以看出,采用本发明所述的干涉光复光场合成方法,可以获得干涉信号光场的复数表达形式,也即干涉信号的复光场,为偏振控制和抑制多通道串扰等提供信号基础。
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,即:
S3.1.2将IC和IS分别进行低通滤波,保留小于fm的频谱成分,获得表达式IC_low和IS_low:
S3.1.3生成相位项复光场表达形式Icom:
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)分别减去直流项,获取三个交流项,即:
S3.2.4合成含有干涉仪相位项的干涉结果复数表达形式Icom:
2.一种根据权利要求1所述光纤干涉仪光场的复数化重构方法,其特征在于:对于含有外差调制信号的问询脉冲信号,通过光学移频方式产生。
3.一种根据权利要求2所述光纤干涉仪光场的复数化重构方法,其特征在于:对于含有外差调制信号的问询脉冲信号,通过两个光学移频器产生。
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)。
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)
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 | 北京大学 | 光纤水听器阵列***和加速度传感器阵列***及测量方法 |
-
2021
- 2021-06-22 CN CN202110688836.5A patent/CN113419341B/zh active Active
Patent Citations (8)
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)
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 |