CN105758328B - 一种纳应变量级的动态光纤应变传感装置 - Google Patents

一种纳应变量级的动态光纤应变传感装置 Download PDF

Info

Publication number
CN105758328B
CN105758328B CN201610326091.7A CN201610326091A CN105758328B CN 105758328 B CN105758328 B CN 105758328B CN 201610326091 A CN201610326091 A CN 201610326091A CN 105758328 B CN105758328 B CN 105758328B
Authority
CN
China
Prior art keywords
mrow
msub
optical fiber
laser
strain
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
CN201610326091.7A
Other languages
English (en)
Other versions
CN105758328A (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.)
Zhonglian Jinguan Information Technology Beijing Co ltd
Original Assignee
Anhui Normal 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 Anhui Normal University filed Critical Anhui Normal University
Priority to CN201610326091.7A priority Critical patent/CN105758328B/zh
Publication of CN105758328A publication Critical patent/CN105758328A/zh
Application granted granted Critical
Publication of CN105758328B publication Critical patent/CN105758328B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/16Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge

Abstract

本发明揭示了一种纳应变量级的动态光纤应变传感装置,激光器连接第一耦合器,所述第一耦合器的两个分光端分别连接参考光纤和传感光纤的一端,所述参考光纤上安装有偏振控制器,所述参考光纤和传感光纤的另一端连接第二耦合器,所述第二耦合器连接数据采集模块,所述数据采集模块输出所采集的信号至数据处理模块,同步触发模块的同步触发输出口分别连接激光器和数据采集模块。

Description

一种纳应变量级的动态光纤应变传感装置
技术领域
本发明涉及应变传感领域,尤其涉及光纤应变传感领域。
背景技术
传统的应变传感器是以应变-电量为基础的,而电信号及其装置易受电磁干扰、湿度等环境影响,在某些情况下甚至不能正常工作。光纤应变传感则几乎不受电磁干扰,环境适应性极强,而且灵敏度高,甚至可以达到nε级别,在超精密监测中具有无可比拟的优势。因此,光纤应变传感系统是一个非常有应用前景和实际意义的课题。目前的纳应变量级动态光纤应变传感装置结构复杂,制作成本高,难以满足生产应用需求。
发明内容
本发明所要解决的技术问题是实现一种结构简单、测量精度高的动态光纤应变传感装置。
为了实现上述目的,本发明采用的技术方案为:一种纳应变量级的动态光纤应变传感装置,激光器连接第一耦合器,所述第一耦合器的两个分光端分别连接参考光纤和传感光纤的一端,所述参考光纤上安装有偏振控制器,所述参考光纤和传感光纤的另一端连接第二耦合器,所述第二耦合器连接数据采集模块,所述数据采集模块输出所采集的信号至数据处理模块,同步触发模块的同步触发输出口分别连接激光器和数据采集模块。
所述激光器为频率可调的稳频窄线宽连续激光器。
所述参考光纤和传感光纤长度不相等。
基于纳应变量级的动态光纤应变传感装置的传感方法:同步触发模块驱动可调谐激光器改变激光频率,同时驱动数据采集模块在每一个激光频率下采集第二耦合器输出的信号,数据处理模块通过对不同时刻的各个频率下采集到的信号做互相关运算,找出使互相关值最大的频率差,解调获得传感光纤探测到的应变。
步骤1、同步触发模块驱动激光器发出的激光产生f0的频移,同时驱动数据采集模块采集M次每次采集N个点,记为data(N,M),并将数据送入数据处理模块经平均运算后得到一个点,记为P1(f0);
步骤2、重复步骤1,使激光依次产生2f0,3f0,…,kf0的频移,得到P1(2f0), P1(3f0),…,P1(kf0);
步骤3、若需探测外界干扰作用在传感光纤上产生的应变,则重复上述步骤 1-3,得到P2(m),其中m=f0,2f0,3f0,…,kf0
步骤4、令
其中,n=tf0,(t+1)f0,…,kf0;t=1,2,3,…,k-q;q是常数,计算P1' 与P2',P3'与P4'的互相关系数,记为:
T1(t)=corrcoef(P1',P2')
T2(t)=corrcoef(P3',P4')
步骤6、找出T1(t)与T2(t)最大时对应的t,若有多个这样的t,则取最小的记为:
t1=min(find(max(T1)==T1))
t2=min(find(max(T2)==T2))
步骤7、获得最优频移:
其中,ξ是传感光纤应变光学校正系数,n是传感光纤纤芯折射率,c是真空中光速;
步骤8、探测到的光纤应变为:
其中,ΔLi是第i次探测到的光纤形变长度,L是探测应变时,发生形变光纤的长度,λ是激光波长,Lh为参考光纤和传感光纤长度之差。
本发明纳应变量级的动态光纤应变传感装置结构简单,几乎不受电磁干扰,环境适应性强,具有nε级别的灵敏度,而且可以动态的实时监测应变变化,可用于大型工程现场以及重大政治、经济、军事基地的周界安防。
附图说明
下面对本发明说明书中每幅附图表达的内容及图中的标记作简要说明:
图1为纳应变量级的动态光纤应变传感装置结构示意图;
上述图中的标记均为:1、激光器;2、第一耦合器;3、偏振控制器;4、参考光纤;5、传感光纤;6、第二耦合器;7、数据采集模块;8、数据处理模块;9、同步触发模块。
具体实施方式
如图1所示,可调谐激光器1连接第一耦合器2,第一耦合器2的两个分光端分别连接参考光纤4和传感光纤5,其中参考光纤4安装有偏振控制器3,第二耦合器6将参考光纤4和传感光纤5合为一束后连接数据采集模块7,数据采集模块7连接数据处理模块8,同步触发模块9的同步触发输出口分别连接可调谐激光器1和数据采集模块7。可调谐激光器1为频率可调的稳频窄线宽连续激光器1。参考光纤4和传感光纤5长度不相等,二者长度之差记为Lh,Lh依据实际应用而定(如Lh=2米)。
上述纳应变量级的动态光纤应变传感装置结构简单,几乎不受电磁干扰,环境适应性强,具有nε级别的灵敏度,而且可以动态的实时监测应变变化。
基于上述纳应变量级的动态光纤应变传感装置,传感方法为:同步触发模块9驱动可调谐激光器1改变激光频率,同时驱动数据采集模块7在每一个激光频率下采集第二耦合器6输出的信号,数据处理模块8通过对不同时刻的各个频率下采集到的信号做互相关运算,找出使互相关值最大的频率差,进而解调出传感光纤5探测到的应变。
具体来说:按本发明搭建传感装置后,调整偏振控制器3,使第二耦合器6 输出的信号最大且稳定,同步触发模块9驱动可调谐激光器1发出的激光产生 f0(如1MHz)的频移,同步触发模块9同时驱动数据采集模块7采集M次(如 500次),每次采集N个(如500个)点,记为data(N,M),将数据送入数据处理模块8经平均运算后得到一个点,记为P1(f0)
重复上述步骤,使激光依次产生2f0,3f0,…,kf0的频移,得到P1(2f0),P1(3f0),…,P1(kf0)。
若需探测外界干扰作用在传感光纤5上产生的应变,则重复上述步骤,得到 P2(m),其中m=f0,2f0,3f0,…,kf0,令
其中,n=tf0,(t+1)f0,…,kf0;t=1,2,3,…,k-q;q是常数(如q=3),计算P1'与P2',P3'与P4'的互相关系数,记为
T1(t)=corrcoef(P1',P2')
T2(t)=corrcoef(P3',P4')
找出T1(t)与T2(t)最大时对应的t,若有多个这样的t,则取最小的记为
t1=min(find(max(T1)==T1))
t2=min(find(max(T2)==T2))
则最优频移即为
其中,ξ是传感光纤5应变光学校正系数,n是传感光纤5纤芯折射率,c是真空中光速,则探测到的光纤应变为
其中,ΔLi是第i次探测到的光纤形变长度,L是探测应变时,发生形变光纤的长度,λ是激光波长。
上面结合附图对本发明进行了示例性描述,显然本发明具体实现并不受上述方式的限制,只要采用了本发明的方法构思和技术方案进行的各种非实质性的改进,或未经改进将本发明的构思和技术方案直接应用于其它场合的,均在本发明的保护范围之内。

Claims (1)

1.一种纳应变量级的动态光纤应变传感装置,其特征在于:激光器连接第一耦合器,所述第一耦合器的两个分光端分别连接参考光纤和传感光纤的一端,所述参考光纤上安装有偏振控制器,所述参考光纤和传感光纤的另一端连接第二耦合器,所述第二耦合器连接数据采集模块,所述数据采集模块输出所采集的信号至数据处理模块,同步触发模块的同步触发输出口分别连接激光器和数据采集模块;
所述激光器为频率可调的稳频窄线宽连续激光器;
所述参考光纤和传感光纤长度不相等;
同步触发模块驱动可调谐激光器改变激光频率,同时驱动数据采集模块在每一个激光频率下采集第二耦合器输出的信号,数据处理模块通过对不同时刻的各个频率下采集到的信号做互相关运算,找出使互相关值最大的频率差,解调获得传感光纤探测到的应变;
所述的纳应变量级的动态光纤应变传感装置的传感方法:
步骤1、同步触发模块驱动激光器发出的激光产生f0的频移,同时驱动数据采集模块采集M次每次采集N个点,记为data(N,M),并将数据送入数据处理模块经平均运算后得到一个点,记为P1(f0);
<mrow> <msub> <mi>P</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>f</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>N</mi> </mfrac> <mo>&amp;CenterDot;</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mrow> <mo>(</mo> <mfrac> <mn>1</mn> <mi>M</mi> </mfrac> <mo>&amp;CenterDot;</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>M</mi> </munderover> <mi>d</mi> <mi>a</mi> <mi>t</mi> <mi>a</mi> <mo>(</mo> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> </mrow>
步骤2、重复步骤1,使激光依次产生2f0,3f0,…,kf0的频移,得到P1(2f0),P1(3f0),…,P1(kf0);
步骤3、若需探测外界干扰作用在传感光纤上产生的应变,则重复上述步骤1-3,得到P2(m),其中m=f0,2f0,3f0,…,kf0
步骤4、令
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>P</mi> <mn>1</mn> <mo>&amp;prime;</mo> </msubsup> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mo>(</mo> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> <msub> <mi>f</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>P</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>P</mi> <mn>2</mn> <mo>&amp;prime;</mo> </msubsup> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mo>(</mo> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> <msub> <mi>f</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <msub> <mi>P</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mo>(</mo> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> <msub> <mi>f</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced>
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>P</mi> <mn>3</mn> <mo>&amp;prime;</mo> </msubsup> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mo>(</mo> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> <msub> <mi>f</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>P</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mo>(</mo> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> <msub> <mi>f</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>P</mi> <mn>4</mn> <mo>&amp;prime;</mo> </msubsup> <mrow> <mo>(</mo> <mi>n</mi> <mo>-</mo> <mo>(</mo> <mrow> <mi>t</mi> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> <msub> <mi>f</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>P</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced>
其中,n=tf0,(t+1)f0,…,kf0;t=1,2,3,…,k-q;q是常数,计算P1'与P2',P3'与P4'的互相关系数,记为:
T1(t)=corrcoef(P1',P2')
T2(t)=corrcoef(P3',P4')
步骤5、找出T1(t)与T2(t)最大时对应的t,若有多个这样的t,则取最小的记为:
t1=min(find(max(T1)==T1))
t2=min(find(max(T2)==T2))
步骤6、获得最优频移:
其中,ξ是传感光纤应变光学校正系数,n是传感光纤纤芯折射率,c是真空中光速;
步骤7、探测到的光纤应变为:
<mrow> <mfrac> <mrow> <msub> <mi>&amp;Delta;L</mi> <mi>i</mi> </msub> </mrow> <mi>L</mi> </mfrac> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>&amp;Delta;L</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>+</mo> <mfrac> <mrow> <mi>&amp;lambda;</mi> <mo>&amp;CenterDot;</mo> <msub> <mi>t</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> </msub> <mo>&amp;CenterDot;</mo> <msub> <mi>f</mi> <mn>0</mn> </msub> </mrow> <mrow> <mi>c</mi> <mo>&amp;CenterDot;</mo> <mi>&amp;xi;</mi> </mrow> </mfrac> <mo>&amp;CenterDot;</mo> <msub> <mi>L</mi> <mi>h</mi> </msub> <mo>)</mo> </mrow> <mo>/</mo> <mi>L</mi> </mrow>
其中,ΔLi是第i次探测到的光纤形变长度,L是探测应变时,发生形变光纤的长度,λ是激光波长,Lh为参考光纤和传感光纤长度之差。
CN201610326091.7A 2016-05-17 2016-05-17 一种纳应变量级的动态光纤应变传感装置 Active CN105758328B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610326091.7A CN105758328B (zh) 2016-05-17 2016-05-17 一种纳应变量级的动态光纤应变传感装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610326091.7A CN105758328B (zh) 2016-05-17 2016-05-17 一种纳应变量级的动态光纤应变传感装置

Publications (2)

Publication Number Publication Date
CN105758328A CN105758328A (zh) 2016-07-13
CN105758328B true CN105758328B (zh) 2018-04-06

Family

ID=56324169

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610326091.7A Active CN105758328B (zh) 2016-05-17 2016-05-17 一种纳应变量级的动态光纤应变传感装置

Country Status (1)

Country Link
CN (1) CN105758328B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103175558A (zh) * 2013-01-17 2013-06-26 广东电网公司电力调度控制中心 分布式光纤传感系统的参数测量装置
JP5413931B2 (ja) * 2013-02-20 2014-02-12 株式会社フジクラ 光ファイバ位置特定のための光学マーキング部を備えた光ファイバセンサおよび光ファイバセンサの計測方法と光ファイバセンサ装置
CN104316217A (zh) * 2014-10-23 2015-01-28 国网河南省电力公司周口供电公司 一种基于布里渊散射自测温光纤复合导线在线监测系统
CN104568219A (zh) * 2015-01-15 2015-04-29 厦门大学 一种基于单通带微波光子滤波器的温度测量装置及方法
CN104613889A (zh) * 2015-02-03 2015-05-13 中国计量学院 一种基于光纤环形激光器的弯曲传感测量系统
CN104819770A (zh) * 2015-05-14 2015-08-05 中国人民解放军国防科学技术大学 基于外差检测和相位解调的相位光时域反射装置及方法
CN205808360U (zh) * 2016-05-17 2016-12-14 安徽师范大学 一种纳应变量级的动态光纤应变传感装置

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103175558A (zh) * 2013-01-17 2013-06-26 广东电网公司电力调度控制中心 分布式光纤传感系统的参数测量装置
JP5413931B2 (ja) * 2013-02-20 2014-02-12 株式会社フジクラ 光ファイバ位置特定のための光学マーキング部を備えた光ファイバセンサおよび光ファイバセンサの計測方法と光ファイバセンサ装置
CN104316217A (zh) * 2014-10-23 2015-01-28 国网河南省电力公司周口供电公司 一种基于布里渊散射自测温光纤复合导线在线监测系统
CN104568219A (zh) * 2015-01-15 2015-04-29 厦门大学 一种基于单通带微波光子滤波器的温度测量装置及方法
CN104613889A (zh) * 2015-02-03 2015-05-13 中国计量学院 一种基于光纤环形激光器的弯曲传感测量系统
CN104819770A (zh) * 2015-05-14 2015-08-05 中国人民解放军国防科学技术大学 基于外差检测和相位解调的相位光时域反射装置及方法
CN205808360U (zh) * 2016-05-17 2016-12-14 安徽师范大学 一种纳应变量级的动态光纤应变传感装置

Also Published As

Publication number Publication date
CN105758328A (zh) 2016-07-13

Similar Documents

Publication Publication Date Title
CN105606193B (zh) 基于双可调谐光源的光纤法珀声振动传感装置及解调方法
CN105488935A (zh) 一种基于非对称双Mach-Zehnder干涉的分布式光纤扰动定位系统及其定位方法
CN107063431A (zh) 一种基于双光路的光纤振动传感系统及方法
CN202547766U (zh) 光纤布拉格光栅振动传感测量系统
CN108692749A (zh) 基于偏振干涉的马赫-曾德干涉仪光程差测量装置及方法
CN109655008A (zh) 一种光纤弯曲传感器及包括该光纤弯曲传感器的智能手套
CN206974448U (zh) 双端探测的联合拉曼和布里渊散射的分布式光纤传感装置
CN103017887A (zh) 一种光纤振动传感系统及其检测方法
CN103308148B (zh) 直线型多通道光纤声振动传感装置及方法
CN105758328B (zh) 一种纳应变量级的动态光纤应变传感装置
CN105783763B (zh) 一种高精度动态光纤应变传感装置及其传感方法
CN108627233A (zh) 一种光纤振动探测系统
CN205785095U (zh) 一种高精度动态光纤应变传感装置
CN106679935A (zh) 一种单模光纤双折射测量的装置及方法
CN201903351U (zh) 一种光纤光栅动态变化解调装置
CN107314823A (zh) 干涉型光纤传感器的相位测量方法及装置
CN205808360U (zh) 一种纳应变量级的动态光纤应变传感装置
CN101324447A (zh) 基于ccd和长周期光纤光栅的布拉格光栅传感解调系统
CN103292831B (zh) 光纤光栅反射滤波结构的单端布里渊光时域分析传感器
CN106404243B (zh) 一种基于光纤光栅偏振信息检测的高频动态信息解调系统以及方法
CN104848879B (zh) 基于线性工作匹配光栅法的光纤布拉格光栅传感器信号解调方法
CN108759982B (zh) 一种基于光混沌的精度可调的光纤液位测量装置及方法
CN207317965U (zh) 基于布里渊散射并融合混沌光信号的分布式传感系统
CN109633235A (zh) 一种多参量传感的光纤电流互感器
CN206177328U (zh) 光纤微腔传感器解调装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20220524

Address after: 3-106-29, 1st floor, building 6, qingshangyuan, Haidian District, Beijing 100085

Patentee after: Zhonglian Jinguan information technology (Beijing) Co.,Ltd.

Address before: 241000 Wuhu Road, Yijiang District, Anhui,

Patentee before: ANHUI NORMAL University

TR01 Transfer of patent right