CN113984126A - 基于不同掺杂双芯弱反射fbg阵列的温度应变监测***和方法 - Google Patents
基于不同掺杂双芯弱反射fbg阵列的温度应变监测***和方法 Download PDFInfo
- Publication number
- CN113984126A CN113984126A CN202111297771.8A CN202111297771A CN113984126A CN 113984126 A CN113984126 A CN 113984126A CN 202111297771 A CN202111297771 A CN 202111297771A CN 113984126 A CN113984126 A CN 113984126A
- Authority
- CN
- China
- Prior art keywords
- optical fiber
- core
- laser
- light source
- weak
- 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
- 238000012544 monitoring process Methods 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 38
- 239000013307 optical fiber Substances 0.000 claims abstract description 109
- 239000000835 fiber Substances 0.000 claims abstract description 88
- 230000003287 optical effect Effects 0.000 claims abstract description 41
- 238000005259 measurement Methods 0.000 claims abstract description 30
- 238000012545 processing Methods 0.000 claims abstract description 29
- 238000005516 engineering process Methods 0.000 claims abstract description 23
- 238000002168 optical frequency-domain reflectometry Methods 0.000 claims abstract description 14
- 238000003491 array Methods 0.000 claims abstract description 5
- 230000008569 process Effects 0.000 claims description 24
- 238000007781 pre-processing Methods 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 10
- 230000000694 effects Effects 0.000 claims description 10
- 230000035945 sensitivity Effects 0.000 claims description 7
- 230000006870 function Effects 0.000 claims description 6
- 238000001228 spectrum Methods 0.000 claims description 6
- 238000013528 artificial neural network Methods 0.000 claims description 5
- 230000035559 beat frequency Effects 0.000 claims description 5
- 230000005540 biological transmission Effects 0.000 claims description 5
- 238000003860 storage Methods 0.000 claims description 5
- 238000010408 sweeping Methods 0.000 claims description 5
- 238000002310 reflectometry Methods 0.000 claims description 3
- 230000008901 benefit Effects 0.000 abstract description 5
- 238000011065 in-situ storage Methods 0.000 abstract description 3
- 238000012360 testing method Methods 0.000 abstract description 2
- 238000001069 Raman spectroscopy Methods 0.000 description 7
- 238000000253 optical time-domain reflectometry Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 238000005452 bending Methods 0.000 description 2
- 238000010801 machine learning Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000012549 training Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000010259 detection of temperature stimulus Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000002019 doping agent Substances 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000010287 polarization Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING 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
- G01D21/00—Measuring or testing not otherwise provided for
- G01D21/02—Measuring two or more variables by means not covered by a single other subclass
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Optical Transform (AREA)
Abstract
本发明公开了基于不同掺杂双芯弱反射FBG阵列的温度应变监测***和方法，包括宽带脉冲光源、多个光纤耦合器、反射镜、光开关、双芯光纤扇出器、双路光电探测器、光纤延时线、信号处理单元、计算机、双芯弱FBG阵列以及光纤；将双纤芯时分复用和OFDR技术相结合，利用不同掺杂纤芯在时域上交替实现对同一测试点原位监测，具有不同温度/应变的敏感系数，通过同一波长的脉冲光源实现信号解耦，既保持了弱光栅可以分布式大容量监测的优势，又可以实现基于宽脉冲激光信号的OFDR监测技术，空间分辨率与动态范围两个指标不再相互制约，借助宽脉冲激光的高能量优点，可以实现长距离、高空间分辨率和大动态范围的测量。
Description
技术领域
本发明涉及光纤传感和分布式多参量监测技术领域，更具体地说，涉及一种基于不同掺杂双芯弱反射FBG阵列的温度/应变监测***。
背景技术
随着科学技术的发展和物联网应用要求的提高，光纤传感网络正在向大容量和多参量测量方向发展，基于瑞利散射、布里渊散射、拉曼散射的分布式光纤传感网络为高电压、强磁场干扰、大电流、复杂几何空间、易燃、易爆等恶劣环境的空间上连续分布各点的温度、应变等物理参量的测量提供了可行的新手段。分布式光纤传感技术是伴随着光时域反射技术(OTDR)的产生而发展起来的，如利用光时域反射技术测量瑞利背向散射信号的强度和偏振态来监测温度/应变；利用光时域反射技术测量拉曼背向散射信号的强度来监测温度；利用光时域反射技术测量布里渊散射信号的强度和频移来监测温度/应变。
受限于技术原理，基于布里渊散和基于拉曼散射的分布式光纤传感***响应速度与空间分辨率都较低，复杂昂贵的***也限制了两类分布式测量技术的工程化应用。基于瑞利散射的分布式光纤传感***响应速度快、灵敏度高，但***的信噪比较低，测量精度与空间分辨率低，难以实现对温度和应变的定量检测等。
目前，可实现同时监测温度和应变参数的光纤传感***主要有以下几种：
1.K·贝达德、R·S戴尔、李洁等提出利用单根光纤实现温度和应变的同时监测(K·贝达德，等.用于同时测量温度和应变的光纤,CN109959987A[P].2019.)，该光纤的纤芯中包含多个掺杂剂浓度分布，该发明便于针对光纤使用更简单的制造和组装工艺，但该发明仅仅讨论在同一位置实现两个参量的监测，没有对不同位置的分布式监测进行讨论。
2.利用并行排列的松套和紧套两根光纤，根据布里渊频移计算温度和应变。该方法需要在工程应用前设计传感光纤结构，不适合己经敷设传感光纤的场合。
3.利用布里渊频移和功率对温度和应变都敏感的特性，根据二元一次方程组解算温度和应变。该方法中布里渊散射功率的低信噪比会导致温度和应变测量精度的大幅下降，而且己敷设传感光纤的功率初值难以确定，这些都制约了该方法的工程应用。
4.利用布里渊频移和拉曼散射信号解算温度和应变。该方法可有效区分温度和应变，但使用单模光纤测量温度的拉曼设备价格昂贵，两台设备更导致成本上升，而且两台设备测量的数据存在融合与配准问题。
相比现有基于散射光信号进行分布式探测的光纤传感器，光纤光栅具有更高的信噪比，传感功能也更加丰富。现有弱光纤光栅阵列虽然在动态振动信号测量方面可实现分布式测量，但是，对温度与应变等静态量信号还只是一种准分布式测量。并且由于光纤对温度与应变的交叉敏感，通过另外加一根只感温光纤的方式很难保证该光纤不受到应力的扰动，温度补偿光栅与测量光栅的位置偏差等都会造成测量精度难以保证，在工程应用中存在困难。
发明内容
为了解决上述问题，针对现有技术的缺陷和不足，本发明提供一种基于不同掺杂双芯弱光栅阵列与瑞利散射结合的多参量分布式传感监测***，利用单根光纤实现空间高密度的温度和应变测量。
本方案的技术解决方案是：
一种基于不同掺杂双芯弱反射FBG阵列的温度应变监测***，包括宽带脉冲光源(1)、第二光纤耦合器(2)、第三光纤耦合器(3)、第四光纤耦合器(4)、第九光纤耦合器(9)、反射镜(5)、光开关(6)、双芯光纤扇出器(8)、双路光电探测器(10)、光纤延时线(7)、信号处理单元(11)、计算机(12)、双芯弱FBG阵列(13)以及一定长度的光纤；
所述宽带脉冲光源(1)通过分光比为1:9的第二光纤耦合器(2)将光源分成两路，其中光功率为9/10的激光进入第一支路，功率为1/10比例的激光进入第二支路；所述第四光纤耦合器(4)、第九光纤耦合器(9)和光纤延时线(7)相连接构成马赫曾德尔干涉仪；所述第三光纤耦合器(3)与反射镜(5)和光开关(6)相连接构成迈克尔逊干涉仪；所述双路光电探测器(10)、信号处理单元(11)和计算机(12)相连接构成传感解调及监测部分；所述光开关(6)通过双芯光纤扇出器(8)与双芯弱FBG阵列(13)连接形成通路；所述计算机通过电缆分别与光开关(6)、宽带脉冲光源(1)连通；
宽带脉冲光源由计算机(12)设置参数后，主路光源由第二光纤耦合器(2)、第三光纤耦合器(3)输出后通过时分方式控制光开关(6)，在时间t11和t21期间分别对不同掺杂的第一纤芯、第二纤芯注入宽带脉冲激光，经双芯弱FBG阵列(13)反射后的激光通过第三光纤耦合器(3)与经反射镜(5)反射后的参考光源发生干涉，干涉后的激光输出给双路光电探测器(10)；第二光纤耦合器(2)输出的辅路光源通过第四光纤耦合器(4)、光纤延时线(7)和第九光纤耦合器(9)形成干涉，干涉后的激光输出给双路光电探测器(10)；双路光电探测器(10)输出的电信号送给信号处理单元(11)，信号处理单元(11)的完成两路电信号的采集、存储以及预处理，信号解调的由计算机中的软件完成，计算机软件功能包括：将双纤芯时分复用和OFDR技术相结合，对时间t12和t22期间的测量数据进行融合，在完成各弱FBG寻址的基础上，基于瑞利散射效应完成各FBG处的温度和应变的测量。
所述的监测***，监测***中的宽带脉冲光源受计算机控制，脉冲激光的重复周期、宽度和波长范围与***不匹配时，可以进行相应微调。
所述的监测***，双芯弱FBG阵列(13)包括第一纤芯和第二纤芯，采用不同折射率系数的单模双芯光纤刻制而成，采用相同的掩模版刻写，各FBG的反射光谱具有相同的3dB带宽、中心波长和超低反射率。
所述的监测***，信号解调由计算机中的软件完成，计算机软件具体包括以下步骤：
步骤1：在t11时间内对第一纤芯注入脉冲激光；
步骤2：在t12时间内对反射回来的传感光信号进行光电转换，并完成第一纤芯的信号解调的预处理，包含第一纤芯上各弱FBG的粗略寻址，基于OFDR技术的弱光栅寻址过程如下：
注入的宽脉冲扫频激光的光场表达式为：
其中v0是宽脉冲扫频激光的初始光频，ξ是光源的扫频速率，A0为光强的振幅，φ(t)是光源自带的相位噪声；
上述宽脉冲扫频激光在其中一纤芯中传输时，经过某一弱光栅反射后，得到以下反射光：
其中：L为某一弱光栅与光源的距离，A(L)表示经过距离2L传输后的反射光的振幅，c为光速，ni(i为0或1)为两不同掺杂光纤的有效折射率；
光电探测器输出的电信号表达式为：
其中k是与光电探测器灵敏度有关系的常数，fb是注入激光与反射激光的差频，等于
式(3)中的第二项是差频部分2πfbt，随时间呈线性变化，属于有用项，对上式进行快速傅里叶变换就可以得到拍频信号的频率大小，然后再结合已知的扫频速率、折射率和光速就可以计算出反射点的位置L；
步骤3：在t12时间内对第二纤芯注入脉冲激光；
步骤4：在t22时间内对反射回来的传感光信号进行光电转换，并完成第二纤芯的信号解调的预处理，包含第二纤芯上各弱FBG的粗略寻址，其计算过程和步骤2中的计算方法相同；
步骤5：在t3时间内，对第一纤芯和第二纤芯的数据进行融合，实现各弱FBG的精确寻址，同时利用瑞利散射效应完成各FBG处的温度和应变，在信号处理过程使用神经网路和人工智能现代数字信号处理技术实现快速提取与温度和应变相关的特征参量，以提高测量的实时性和精确性；
步骤6：根据步骤5微调脉冲激光的参数，以更好满足监测***的需求；
步骤7：重复以上步骤1-步骤5的过程。
基于任一所述的温度应变监测***的温度应变监测方法，主路光源由第二光纤耦合器(2)、第三光纤耦合器(3)输出后通过时分方式控制光开关(6)，在时间t11和t21期间分别对不同掺杂的第一纤芯、第二纤芯注入宽带脉冲激光，经双芯弱FBG阵列(13)反射后的激光通过第三光纤耦合器(3)与经反射镜(5)反射后的参考光源发生干涉，干涉后的激光输出给双路光电探测器(10)；第二光纤耦合器(2)输出的辅路光源通过第四光纤耦合器(4)、光纤延时线(7)和第九光纤耦合器(9)形成干涉，干涉后的激光输出给双路光电探测器(10)；双路光电探测器(10)输出的电信号送给信号处理单元(11)，信号处理单元(11)的完成两路电信号的采集、存储以及预处理，信号解调的由计算机中的软件完成，计算机软件功能包括：将双纤芯时分复用和OFDR技术相结合，对时间t12和t22期间的测量数据进行融合，在完成各弱FBG寻址的基础上，基于瑞利散射效应完成各FBG处的温度和应变的测量。
所述的温度应变监测方法，信号解调由计算机中的软件完成，计算机软件具体包括以下步骤：
步骤1：在t11时间内对第一纤芯注入脉冲激光；
步骤2：在t12时间内对反射回来的传感光信号进行光电转换，并完成第一纤芯的信号解调的预处理，包含第一纤芯上各弱FBG的粗略寻址，基于OFDR技术的弱光栅寻址过程如下：
注入的宽脉冲扫频激光的光场表达式为：
其中v0是宽脉冲扫频激光的初始光频，ξ是光源的扫频速率，A0为光强的振幅，φ(t)是光源自带的相位噪声；
上述宽脉冲扫频激光在其中一纤芯中传输时，经过某一弱光栅反射后，得到以下反射光：
其中：L为某一弱光栅与光源的距离，A(L)表示经过距离2L传输后的反射光的振幅，c为光速，ni(i为0或1)为两不同掺杂光纤的有效折射率；
光电探测器输出的电信号表达式为：
其中k是与光电探测器灵敏度有关系的常数，fb是注入激光与反射激光的差频，等于
式(3)中的第二项是差频部分2πfbt，随时间呈线性变化，属于有用项，对上式进行快速傅里叶变换就可以得到拍频信号的频率大小，然后再结合已知的扫频速率、折射率和光速就可以计算出反射点的位置L；
步骤3：在t12时间内对第二纤芯注入脉冲激光；
步骤4：在t22时间内对反射回来的传感光信号进行光电转换，并完成第二纤芯的信号解调的预处理，包含第二纤芯上各弱FBG的粗略寻址，其计算过程和步骤2中的计算方法相同；
步骤5：在t3时间内，对第一纤芯和第二纤芯的数据进行融合，实现各弱FBG的精确寻址，同时利用瑞利散射效应完成各FBG处的温度和应变，在信号处理过程使用神经网路和人工智能现代数字信号处理技术实现快速提取与温度和应变相关的特征参量，以提高测量的实时性和精确性；
步骤6：根据步骤5微调脉冲激光的参数，以更好满足监测***的需求；
步骤7：重复以上步骤1-步骤5的过程。
由于采用了上述的技术方案，本发明具有如下优越性：
1、本发明将双纤芯时分复用和OFDR技术相结合，利用不同掺杂纤芯在时域上交替实现对同一测试点原位监测，具有不同温度/应变的敏感系数，通过同一波长的脉冲光源实现信号解耦，既保持了弱光栅可以分布式大容量监测的优势，又可以实现基于宽脉冲激光信号的OFDR监测技术，空间分辨率与动态范围两个指标不再相互制约，借助宽脉冲激光的高能量优点，可以实现长距离、高空间分辨率和大动态范围的测量。
2、本发明相对于拉曼和布里渊分布式传感***，本***结构简单，响应速度快，空间分辨率高，项目的研究成果将可直接应用于油气管线、高速铁路、周界安防等领域的监测，具有很高的经济与社会价值，应用前景十分广泛。
附图说明
图1为***组成原理框图；
图2为不同掺杂双芯光纤光栅；
图3为宽带脉冲激光注入和信号处理的逻辑时序图；
图4为利用FBG光谱样本训练GPR模型；
图5为利用GPR模型预测温度和应变；
1、宽带脉冲光源、2、第二光纤耦合器、3、第三光纤耦合器、4、第四光纤耦合器、9、第九光纤耦合器、5、反射镜、6、光开关、7、延时线、8、双芯光纤扇出器、10、双路光电探测器、11、信号处理单元、12、计算机、13、双芯弱FBG阵列；
具体实施方式
以下结合具体实施例，对本发明进行详细说明。
如图1所示，一种基于不同掺杂双芯弱反射FBG阵列的温度应变监测***，包括宽带脉冲光源1、多个光纤耦合器(第二光纤耦合器2、第三光纤耦合器3、第四光纤耦合器4、第九光纤耦合器9)、反射镜(5)、光开关(6)、双芯光纤扇出器(8)、双路光电探测器(10)、光纤延时线(7)、信号处理单元(11)、计算机(12)、双芯弱FBG阵列(13)以及一定长度的光纤；
所述宽带脉冲光源1通过分光比为1:9的第二光纤耦合器2将光源分成两路，其中光功率为9/10的激光进入图1中的上支路(第一支路)，功率为1/10比例的激光进入图1中的下支路(第二支路)；所述第四光纤耦合器4、第九光纤耦合器9和光纤延时线7相连接构成马赫曾德尔干涉仪；所述第三光纤耦合器(3)与反射镜5和光开关6相连接构成迈克尔逊干涉仪；所述双路光电探测器10、信号处理单元11和计算机12相连接构成传感解调及监测部分；所述光开关6通过双芯光纤扇出器8与双芯弱FBG阵列13连接形成通路；所述计算机通过电缆分别与光开关6、宽带脉冲光源1连通。
监测***中的宽带脉冲光源受计算机控制，脉冲激光的重复周期、宽度和波长范围与***不匹配时，可以进行相应微调；双芯弱FBG阵列13包括第一纤芯和第二纤芯，采用不同折射率系数的单模双芯光纤刻制而成，如图2所示，采用相同的掩模版刻写，各FBG的反射光谱具有相同的3dB带宽、中心波长和超低反射率。
监测***的工作原理是：宽带脉冲光源由计算机12设置参数后，主路光源由第二光纤耦合器2、第三光纤耦合器3输出后通过时分方式控制光开关6，如图3所示，在时间t11和t21期间分别对不同掺杂的两个纤芯(第一纤芯、第二纤芯)注入宽带脉冲激光，经双芯弱FBG阵列13反射后的激光通过第三光纤耦合器(3)与经反射镜5反射后的参考光源发生干涉，干涉后的激光输出给双路光电探测器10；第二光纤耦合器(2)输出的辅路光源通过第四光纤耦合器(4)、光纤延时线7和第九光纤耦合器(9)形成干涉，干涉后的激光输出给双路光电探测器10；双路光电探测器10输出的电信号送给信号处理单元11，信号处理单元11的完成两路电信号的采集、存储以及预处理，信号解调的由计算机中的软件完成，计算机软件功能主要包括：将双纤芯时分复用和OFDR技术相结合，对时间t12和t22期间的测量数据进行融合，在完成各弱FBG寻址的基础上，基于瑞利散射效应完成各FBG处的温度和应变的测量。具体包括以下步骤：
步骤1：在t11时间内对第一纤芯注入脉冲激光；
步骤2：在t12时间内对反射回来的传感光信号进行光电转换，并完成第一纤芯的信号解调的预处理，包含第一纤芯上各弱FBG的粗略寻址，基于OFDR技术的弱光栅寻址过程如下：
注入的宽脉冲扫频激光的光场表达式为：
其中v0是宽脉冲扫频激光的初始光频，ξ是光源的扫频速率，A0为光强的振幅，φ(t)是光源自带的相位噪声。
上述宽脉冲扫频激光在其中一纤芯中传输时，经过某一弱光栅反射后，得到以下反射光：
其中：L为某一弱光栅与光源的距离，A(L)表示经过距离2L传输后的反射光的振幅，c为光速，ni(i为0或1)为两不同掺杂光纤的有效折射率。
光电探测器输出的电信号表达式为：
其中k是与光电探测器灵敏度有关系的常数，fb是注入激光与反射激光的差频，等于
式(3)中的第二项是差频部分2πfbt，随时间呈线性变化，属于有用项，对上式进行快速傅里叶变换就可以得到拍频信号的频率大小，然后再结合已知的扫频速率、折射率和光速就可以计算出反射点的位置L。
步骤3：在t12时间内对第二纤芯注入脉冲激光；
步骤4：在t22时间内对反射回来的传感光信号进行光电转换，并完成第二纤芯的信号解调的预处理，包含第二纤芯上各弱FBG的粗略寻址，其计算过程和步骤2中的计算方法相同；
步骤5：在t3时间内，对第一纤芯和第二纤芯的数据进行融合，实现各弱FBG的精确寻址，同时利用瑞利散射效应完成各FBG处的温度和应变，在信号处理过程使用神经网路和人工智能现代数字信号处理技术实现快速提取与温度和应变相关的特征参量，以提高测量的实时性和精确性；
步骤6：根据步骤5微调脉冲激光的参数，以更好满足监测***的需求；
步骤7：重复以上步骤1-步骤5的过程。
其中基于机器学习的弱FBG光谱智能识别算法如下：
第1步为训练过程：按照图4的方法，利用实验测试的双芯弱FBG反射光谱的数据进行训练，通过机器学习得到高斯过程回归(GPR)模型；
第2步为预测过程：参考图5，将监测到的弱FBG光谱输入到已经训练完成的GPR模型中，其输出结果为温度或应变。
其中对第一纤芯和第二纤芯的数据进行融合，具体工作分为两步：
步骤1：考虑到第一纤芯1和第二纤芯2的数据是时分测量的，即不是严格意义上的同时测量，需以时间t为基准，对时间滞后的第二纤芯2数据的时间轴添加偏移量Δt，将对准后的得到的光谱信息作为图4中的输入数据；
步骤2：在实际工程布线中，双芯光纤处于任意弯曲状态，弯曲方向具有随意性，在对弱FBG寻址时，为了减少误差，将两个纤芯的位置测量结果的平均值
作为最终弱FBG的寻址位置。
应当理解的是，对本领域普通技术人员来说，可以根据上述说明加以改进或变换，而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (8)
1.一种基于不同掺杂双芯弱反射FBG阵列的温度应变监测***，其特征在于，包括宽带脉冲光源(1)、第二光纤耦合器(2)、第三光纤耦合器(3)、第四光纤耦合器(4)、第九光纤耦合器(9)、反射镜(5)、光开关(6)、双芯光纤扇出器(8)、双路光电探测器(10)、光纤延时线(7)、信号处理单元(11)、计算机(12)、双芯弱FBG阵列(13)以及一定长度的光纤；
所述宽带脉冲光源(1)通过分光比为1:9的第二光纤耦合器(2)将光源分成两路，其中光功率为9/10的激光进入第一支路，功率为1/10比例的激光进入第二支路；所述第四光纤耦合器(4)、第九光纤耦合器(9)和光纤延时线(7)相连接构成马赫曾德尔干涉仪；所述第三光纤耦合器(3)与反射镜(5)和光开关(6)相连接构成迈克尔逊干涉仪；所述双路光电探测器(10)、信号处理单元(11)和计算机(12)相连接构成传感解调及监测部分；所述光开关(6)通过双芯光纤扇出器(8)与双芯弱FBG阵列(13)连接形成通路；所述计算机通过电缆分别与光开关(6)、宽带脉冲光源(1)连通。
2.根据权利要求1所述的监测***，其特征在于，主路光源由第二光纤耦合器(2)、第三光纤耦合器(3)输出后通过时分方式控制光开关(6)，在时间t11和t21期间分别对不同掺杂的第一纤芯、第二纤芯注入宽带脉冲激光，经双芯弱FBG阵列(13)反射后的激光通过第三光纤耦合器(3)与经反射镜(5)反射后的参考光源发生干涉，干涉后的激光输出给双路光电探测器(10)；第二光纤耦合器(2)输出的辅路光源通过第四光纤耦合器(4)、光纤延时线(7)和第九光纤耦合器(9)形成干涉，干涉后的激光输出给双路光电探测器(10)；双路光电探测器(10)输出的电信号送给信号处理单元(11)，信号处理单元(11)的完成两路电信号的采集、存储以及预处理，信号解调的由计算机中的软件完成。
3.根据权利要求1所述的监测***，其特征在于，计算机软件功能包括：将双纤芯时分复用和OFDR技术相结合，对时间t12和t22期间的测量数据进行融合，在完成各弱FBG寻址的基础上，基于瑞利散射效应完成各FBG处的温度和应变的测量。
4.根据权利要求1所述的监测***，其特征在于，监测***中的宽带脉冲光源受计算机控制，脉冲激光的重复周期、宽度和波长范围与***不匹配时，可以进行相应微调。
5.根据权利要求1所述的监测***，其特征在于，双芯弱FBG阵列(13)包括第一纤芯和第二纤芯，采用不同折射率系数的单模双芯光纤刻制而成，采用相同的掩模版刻写，各FBG的反射光谱具有相同的3dB带宽、中心波长和超低反射率。
6.根据权利要求1所述的监测***，其特征在于，信号解调由计算机中的软件完成，计算机软件具体包括以下步骤：
步骤1：在t11时间内对第一纤芯注入脉冲激光；
步骤2：在t12时间内对反射回来的传感光信号进行光电转换，并完成第一纤芯的信号解调的预处理，包含第一纤芯上各弱FBG的粗略寻址，基于OFDR技术的弱光栅寻址过程如下：
注入的宽脉冲扫频激光的光场表达式为：
其中v0是宽脉冲扫频激光的初始光频，ξ是光源的扫频速率，A0为光强的振幅，φ(t)是光源自带的相位噪声；
上述宽脉冲扫频激光在其中一纤芯中传输时，经过某一弱光栅反射后，得到以下反射光：
其中：L为某一弱光栅与光源的距离，A(L)表示经过距离2L传输后的反射光的振幅，c为光速，ni(i为0或1)为两不同掺杂光纤的有效折射率；
光电探测器输出的电信号表达式为：
其中k是与光电探测器灵敏度有关系的常数，fb是注入激光与反射激光的差频，等于
式(3)中的第二项是差频部分2πfbt，随时间呈线性变化，属于有用项，对上式进行快速傅里叶变换就可以得到拍频信号的频率大小，然后再结合已知的扫频速率、折射率和光速就可以计算出反射点的位置L；
步骤3：在t12时间内对第二纤芯注入脉冲激光；
步骤4：在t22时间内对反射回来的传感光信号进行光电转换，并完成第二纤芯的信号解调的预处理，包含第二纤芯上各弱FBG的粗略寻址，其计算过程和步骤2中的计算方法相同；
步骤5：在t3时间内，对第一纤芯和第二纤芯的数据进行融合，实现各弱FBG的精确寻址，同时利用瑞利散射效应完成各FBG处的温度和应变，在信号处理过程使用神经网路和人工智能现代数字信号处理技术实现快速提取与温度和应变相关的特征参量，以提高测量的实时性和精确性；
步骤6：根据步骤5微调脉冲激光的参数，以更好满足监测***的需求；
步骤7：重复以上步骤1-步骤5的过程。
7.基于权利要求1-6任一所述的温度应变监测***的温度应变监测方法，其特征在于，主路光源由第二光纤耦合器(2)、第三光纤耦合器(3)输出后通过时分方式控制光开关(6)，在时间t11和t21期间分别对不同掺杂的第一纤芯、第二纤芯注入宽带脉冲激光，经双芯弱FBG阵列(13)反射后的激光通过第三光纤耦合器(3)与经反射镜(5)反射后的参考光源发生干涉，干涉后的激光输出给双路光电探测器(10)；第二光纤耦合器(2)输出的辅路光源通过第四光纤耦合器(4)、光纤延时线(7)和第九光纤耦合器(9)形成干涉，干涉后的激光输出给双路光电探测器(10)；双路光电探测器(10)输出的电信号送给信号处理单元(11)，信号处理单元(11)的完成两路电信号的采集、存储以及预处理，信号解调的由计算机中的软件完成，计算机软件功能包括：将双纤芯时分复用和OFDR技术相结合，对时间t12和t22期间的测量数据进行融合，在完成各弱FBG寻址的基础上，基于瑞利散射效应完成各FBG处的温度和应变的测量。
8.根据权利要求5所述的温度应变监测方法，其特征在于，信号解调由计算机中的软件完成，计算机软件具体包括以下步骤：
步骤1：在t11时间内对第一纤芯注入脉冲激光；
步骤2：在t12时间内对反射回来的传感光信号进行光电转换，并完成第一纤芯的信号解调的预处理，包含第一纤芯上各弱FBG的粗略寻址，基于OFDR技术的弱光栅寻址过程如下：
注入的宽脉冲扫频激光的光场表达式为：
其中v0是宽脉冲扫频激光的初始光频，ξ是光源的扫频速率，A0为光强的振幅，φ(t)是光源自带的相位噪声；
上述宽脉冲扫频激光在其中一纤芯中传输时，经过某一弱光栅反射后，得到以下反射光：
其中：L为某一弱光栅与光源的距离，A(L)表示经过距离2L传输后的反射光的振幅，c为光速，ni(i为0或1)为两不同掺杂光纤的有效折射率；
光电探测器输出的电信号表达式为：
其中k是与光电探测器灵敏度有关系的常数，fb是注入激光与反射激光的差频，等于
式(3)中的第二项是差频部分2πfbt，随时间呈线性变化，属于有用项，对上式进行快速傅里叶变换就可以得到拍频信号的频率大小，然后再结合已知的扫频速率、折射率和光速就可以计算出反射点的位置L；
步骤3：在t12时间内对第二纤芯注入脉冲激光；
步骤4：在t22时间内对反射回来的传感光信号进行光电转换，并完成第二纤芯的信号解调的预处理，包含第二纤芯上各弱FBG的粗略寻址，其计算过程和步骤2中的计算方法相同；
步骤5：在t3时间内，对第一纤芯和第二纤芯的数据进行融合，实现各弱FBG的精确寻址，同时利用瑞利散射效应完成各FBG处的温度和应变，在信号处理过程使用神经网路和人工智能现代数字信号处理技术实现快速提取与温度和应变相关的特征参量，以提高测量的实时性和精确性；
步骤6：根据步骤5微调脉冲激光的参数，以更好满足监测***的需求；
步骤7：重复以上步骤1-步骤5的过程。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111297771.8A CN113984126B (zh) | 2021-11-04 | 2021-11-04 | 基于不同掺杂双芯弱反射fbg阵列的温度应变监测***和方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111297771.8A CN113984126B (zh) | 2021-11-04 | 2021-11-04 | 基于不同掺杂双芯弱反射fbg阵列的温度应变监测***和方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113984126A true CN113984126A (zh) | 2022-01-28 |
| CN113984126B CN113984126B (zh) | 2024-05-14 |
Family
ID=79746327
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111297771.8A Active CN113984126B (zh) | 2021-11-04 | 2021-11-04 | 基于不同掺杂双芯弱反射fbg阵列的温度应变监测***和方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113984126B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116989832A (zh) * | 2023-09-26 | 2023-11-03 | 北京大学 | 一种可变多空间传感尺度地震勘探方法及*** |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101158591A (zh) * | 2007-10-15 | 2008-04-09 | 北京航空航天大学 | 适用于光纤分布式温度和应力传感装置的探测方法 |
| JP2010054366A (ja) * | 2008-08-28 | 2010-03-11 | Fujikura Ltd | 光ファイバ位置特定のための光学マーキング部を備えた光ファイバセンサおよび光ファイバセンサの計測方法と光ファイバセンサ装置 |
| CA2813869A1 (en) * | 2010-10-14 | 2012-04-19 | Fiber Sensys, Inc. | Interferometer systems |
| CN102980681A (zh) * | 2012-11-16 | 2013-03-20 | 暨南大学 | 一种基于布里渊散射的分布式应变和温度光纤传感器 |
| CN104180831A (zh) * | 2013-05-24 | 2014-12-03 | 无锡万润光子技术有限公司 | 基于双芯光纤灵敏度增强型光时域反射分布式Michelson干涉仪 |
| CN105698871A (zh) * | 2016-03-29 | 2016-06-22 | 天津大学 | 基于光频域反射的分布式应变温度同时测量装置及方法 |
| CN105698831A (zh) * | 2016-01-26 | 2016-06-22 | 武汉理工大学 | 双芯光纤光栅阵列传感网络及分布式传感信息获取方法 |
| CN106840221A (zh) * | 2017-01-06 | 2017-06-13 | 武汉理工大学 | 基于色散马赫曾德尔干涉的光纤光栅解调装置及方法 |
| CN107036733A (zh) * | 2017-04-21 | 2017-08-11 | 武汉理工大学 | 基于暗脉冲光源的双芯弱光栅阵列的多参量分布式测量***及测量方法 |
| CN108007603A (zh) * | 2017-12-18 | 2018-05-08 | 武汉理工大学 | 一种基于非对称双芯光纤的多参量分布测量*** |
| US20190063963A1 (en) * | 2015-11-06 | 2019-02-28 | Febus Optics | Optoelectronic distributed measuring device based on brillouin scattering |
| CN109959403A (zh) * | 2019-03-29 | 2019-07-02 | 武汉理工大学 | 一种多参量大容量传感*** |
-
2021
- 2021-11-04 CN CN202111297771.8A patent/CN113984126B/zh active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101158591A (zh) * | 2007-10-15 | 2008-04-09 | 北京航空航天大学 | 适用于光纤分布式温度和应力传感装置的探测方法 |
| JP2010054366A (ja) * | 2008-08-28 | 2010-03-11 | Fujikura Ltd | 光ファイバ位置特定のための光学マーキング部を備えた光ファイバセンサおよび光ファイバセンサの計測方法と光ファイバセンサ装置 |
| CA2813869A1 (en) * | 2010-10-14 | 2012-04-19 | Fiber Sensys, Inc. | Interferometer systems |
| CN102980681A (zh) * | 2012-11-16 | 2013-03-20 | 暨南大学 | 一种基于布里渊散射的分布式应变和温度光纤传感器 |
| CN104180831A (zh) * | 2013-05-24 | 2014-12-03 | 无锡万润光子技术有限公司 | 基于双芯光纤灵敏度增强型光时域反射分布式Michelson干涉仪 |
| US20190063963A1 (en) * | 2015-11-06 | 2019-02-28 | Febus Optics | Optoelectronic distributed measuring device based on brillouin scattering |
| CN105698831A (zh) * | 2016-01-26 | 2016-06-22 | 武汉理工大学 | 双芯光纤光栅阵列传感网络及分布式传感信息获取方法 |
| CN105698871A (zh) * | 2016-03-29 | 2016-06-22 | 天津大学 | 基于光频域反射的分布式应变温度同时测量装置及方法 |
| CN106840221A (zh) * | 2017-01-06 | 2017-06-13 | 武汉理工大学 | 基于色散马赫曾德尔干涉的光纤光栅解调装置及方法 |
| CN107036733A (zh) * | 2017-04-21 | 2017-08-11 | 武汉理工大学 | 基于暗脉冲光源的双芯弱光栅阵列的多参量分布式测量***及测量方法 |
| CN108007603A (zh) * | 2017-12-18 | 2018-05-08 | 武汉理工大学 | 一种基于非对称双芯光纤的多参量分布测量*** |
| CN109959403A (zh) * | 2019-03-29 | 2019-07-02 | 武汉理工大学 | 一种多参量大容量传感*** |
Non-Patent Citations (3)
| Title |
|---|
| CHUNJIAO JIN; XINGLIN TONG; QIAO XIONG: "An spectral domain OCT system using All-fiber probe", 2019 18TH INTERNATIONAL CONFERENCE ON OPTICAL COMMUNICATIONS AND NETWORKS, 19 December 2019 (2019-12-19), pages 1 - 3 * |
| MARIANO MAHISSI, XINGLIN TONG, CUI ZHANG, ET AL.: "Study on the vibration performances for a high temperature fiber F-P accelerometer", OPTICAL FIBER TECHNOLOGY, vol. 62, pages 1 - 10 * |
| 桂鑫, 李政颖, 王洪海等: "基于大规模光栅阵列光纤的分布式传感技术及应用综述", 应用科学学报, vol. 39, no. 05, pages 747 - 776 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116989832A (zh) * | 2023-09-26 | 2023-11-03 | 北京大学 | 一种可变多空间传感尺度地震勘探方法及*** |
| CN116989832B (zh) * | 2023-09-26 | 2023-12-26 | 北京大学 | 一种可变多空间传感尺度地震勘探方法及*** |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113984126B (zh) | 2024-05-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Sun et al. | Distributed fiber-optic vibration sensor using a ring Mach-Zehnder interferometer | |
| Zhu et al. | Improved Φ-OTDR sensing system for high-precision dynamic strain measurement based on ultra-weak fiber Bragg grating array | |
| CN101764646B (zh) | 利用波长编码的光时域反射测试装置对光纤链路进行的测量方法 | |
| CN102052930B (zh) | 光纤光栅分布式应变传感器及其应变监测方法 | |
| CN102279095B (zh) | 一种减小双折射色散对保偏光纤偏振耦合测量影响的装置 | |
| CN109959403B (zh) | 一种多参量大容量传感*** | |
| CN105973450A (zh) | 光纤Fizeau干涉阵列分布式振动传感***及方法 | |
| CN111442789A (zh) | 一种基于模式复用提高传感***空间分辨率与测量精度的方法 | |
| Tong et al. | High-speed Mach-Zehnder-OTDR distributed optical fiber vibration sensor using medium-coherence laser | |
| CN101634571A (zh) | 光纤脉栅分布传感装置 | |
| CN102809387B (zh) | 一种botdr信号解调方法 | |
| CN113984126B (zh) | 基于不同掺杂双芯弱反射fbg阵列的温度应变监测***和方法 | |
| Wang et al. | Dual-wavelength Michelson interferometer employing time delay estimation for distributed disturbance location | |
| CN113804412B (zh) | 一种基于环行光路结构的光纤器件微量色度色散测量装置 | |
| Xiao et al. | Frequency response enhancement of Φ-OTDR using interval-sweeping pulse equivalent sampling based on compressed sensing | |
| CN105021330A (zh) | 碳纤维增强型智能钢绞线、预应力监测装置及方法 | |
| Xu et al. | Bending-loss-resistant distributed temperature and strain discriminative Brillouin sensor based on 98 mol% germania-doped few-mode fiber | |
| Shaheen et al. | Phase-sensitive optical time domain reflectometry based on geometric phase measurement | |
| CN109141487A (zh) | 一种分布式光纤传感器 | |
| Vasiljević-Toskić et al. | A simple fiber optic sensing system based on dual Sagnac interferometer for disturbance location detection | |
| Magalhães et al. | Curvature sensor based on a long-period grating in a fiber ring resonator interrogated by an OTDR | |
| CN111537010A (zh) | 基于otdr的f-p干涉型传感头多点测量方法及装置 | |
| Tong et al. | Improved distributed optical fiber vibration sensor based on Mach-Zehnder-OTDR | |
| Wang et al. | High sensitivity distributed static strain sensing based on all grating optical fiber in optical frequency domain reflectometry | |
| Sun et al. | Distributed fiber-optic sensor with a ring Mach-Zehnder interferometer |
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 | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20231205 Address after: No. 1 Fuyin Road, Shidao Management District, Rongcheng City, Weihai City, Shandong Province, 264309 Applicant after: Weihai Research Institute of Wuhan University of Technology Address before: 276826 No. 99, Tianjin Road, development zone, Rizhao City, Shandong Province Applicant before: Rizhao Wuhan University of science and Technology Institute of biomedicine and new materials |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |