CN111147133A - 一种基于φ-OTDR的车流量实时监测系统及方法 - Google Patents
一种基于φ-OTDR的车流量实时监测系统及方法 Download PDFInfo
- Publication number
- CN111147133A CN111147133A CN201911345757.3A CN201911345757A CN111147133A CN 111147133 A CN111147133 A CN 111147133A CN 201911345757 A CN201911345757 A CN 201911345757A CN 111147133 A CN111147133 A CN 111147133A
- Authority
- CN
- China
- Prior art keywords
- signal
- vibration
- optical
- light
- optical fiber
- 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 82
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000000253 optical time-domain reflectometry Methods 0.000 title 1
- 230000003287 optical effect Effects 0.000 claims abstract description 100
- 239000013307 optical fiber Substances 0.000 claims abstract description 89
- 230000008569 process Effects 0.000 claims abstract description 15
- 238000004891 communication Methods 0.000 claims abstract description 6
- 230000035559 beat frequency Effects 0.000 claims description 22
- 230000001427 coherent effect Effects 0.000 claims description 16
- 238000010586 diagram Methods 0.000 claims description 16
- 239000000835 fiber Substances 0.000 claims description 14
- 238000001514 detection method Methods 0.000 claims description 10
- 239000013598 vector Substances 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 9
- 238000001914 filtration Methods 0.000 claims description 8
- 238000009826 distribution Methods 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 4
- 230000008878 coupling Effects 0.000 claims description 4
- 238000010168 coupling process Methods 0.000 claims description 4
- 238000005859 coupling reaction Methods 0.000 claims description 4
- 230000010287 polarization Effects 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 abstract description 4
- 230000000694 effects Effects 0.000 description 4
- 238000002310 reflectometry Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000009827 uniform distribution Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/07—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems
- H04B10/071—Arrangements for monitoring or testing transmission systems; Arrangements for fault measurement of transmission systems using a reflected signal, e.g. using optical time domain reflectometers [OTDR]
-
- 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
- G01D5/00—Mechanical 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/26—Mechanical 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/32—Mechanical 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/34—Mechanical 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/353—Mechanical 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/35338—Mechanical 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 using other arrangements than interferometer arrangements
- G01D5/35354—Sensor working in reflection
- G01D5/35358—Sensor working in reflection using backscattering to detect the measured quantity
- G01D5/35361—Sensor working in reflection using backscattering to detect the measured quantity using elastic backscattering to detect the measured quantity, e.g. using Rayleigh backscattering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
- G01H9/004—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means using fibre optic sensors
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/90—Determination of colour characteristics
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/0104—Measuring and analyzing of parameters relative to traffic conditions
- G08G1/0108—Measuring and analyzing of parameters relative to traffic conditions based on the source of data
- G08G1/0116—Measuring and analyzing of parameters relative to traffic conditions based on the source of data from roadside infrastructure, e.g. beacons
-
- G—PHYSICS
- G08—SIGNALLING
- G08G—TRAFFIC CONTROL SYSTEMS
- G08G1/00—Traffic control systems for road vehicles
- G08G1/01—Detecting movement of traffic to be counted or controlled
- G08G1/04—Detecting movement of traffic to be counted or controlled using optical or ultrasonic detectors
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/60—Receivers
- H04B10/61—Coherent receivers
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Traffic Control Systems (AREA)
Abstract
Description
技术领域
背景技术
伴随我国经济不断发展,道路交通建设不断完善,对公路上车流量进行实时监测逐渐成为交通领域的研究热点。一方面,城市智能交通已得到社会各界的广泛关注,这样一个稳定可靠能全天24小时对公路路况车流量进行实时监测的系统对于保障交通安全畅通和智能控制提供了基础;另一方面,一个方便的车流量实时监测系统也为广大司乘人员出行提供了参考便利。
发明内容
本发明解决其技术问题所采用的技术方案是:
传感光缆,平行铺设在监测路段的道路两侧,对道路车辆行驶过程中引起的振动信号进行实时监测;
基于的光纤振动监测主机,与传感光缆相连,用于获取光纤振动信号,将光纤振动信号转换为光纤内的瑞利后向散射信号;通过相干检测方法,将后向瑞利散射信号与本地光源信号进行拍频,并提取拍频信号的交流部分信号,对交流部分信号进行采样,得到数字信号,在数字域中实现幅度和相位的解调,解调得到光信号强度;
监控端,用于根据光信号强度建立信号瀑布图,通过信号瀑布图实时监控公路上车辆行驶引起的振动信号。
当传感光缆的光纤线路上发生外界扰动时,光纤相应位置的折射率发生变化,导致该处光相位发生变化,引起后向散射光光强发生变化,进而得到从光纤中散射回来的瑞利后向散射信号。
激光器作为光源提供窄线宽的连续光,由第一耦合器分路为本地光与调制光;
声光调制器将调制光路连续光斩波调制为连续脉冲光;
掺铒光纤放大器将调制光路脉冲光进行放大;
光衰减器用于调节光功率;
调制光路中脉冲光与环形器第一端口连接,第二端口与传感光缆连接,第三端口即瑞利散射信号光;
第二耦合器两输入端为本地光与瑞利散射信号光,两者相干后经过平衡探测器输出为电信号传输至数字采集卡;
数字采集卡对信号进行数据采集,采集数据送入工控板进一步分析。
进一步地,本发明的传感光缆为普通通信单模光缆。
进一步地,本发明的监控端展示车流量形式为光缆振动信号瀑布图,通过对振动信号瀑布图进行分析,得到车辆数量、车辆行驶速度的信息;监控端还包括远程移动终端。
S1、通过平行铺设在监测路段道路两侧的传感光缆实时监测道路车辆行驶过程引起的光纤振动信号;
S2、通过基于的光纤振动监测主机获取光纤振动信号,将光纤振动信号转换为光纤内的瑞利后向散射信号:当传感光缆的光纤线路上发生外界扰动时,光纤相应位置的折射率发生变化,导致该处光相位发生变化,引起后向散射光光强发生变化,进而得到从光纤中散射回来的瑞利后向散射信号;
S3、通过相干检测方法,将后向瑞利散射信号与本地光源信号进行拍频,通过平衡探测器得到拍频信号的交流部分信号,并对交流部分信号进行采样,得到数字信号,在数字域中实现幅度和相位的解调,解调得到光信号强度;
S4、将解调得到的光信号强度发送给监控端,建立信号瀑布图,信号瀑布图为光信号强度在空间、时间域上的实时分布情况;信号瀑布图中横坐标为光纤距离,纵坐标为时间,信号瀑布图中单点颜色反映光信号强度,设定阈值,对超过阈值信号视为振动信号,光信号强度越大与无振动点背景信号颜色对比度越大;通过信号瀑布图实现对车流量、车重量的实时显示。
进一步地,本发明的步骤S2中得到的瑞利后向散射信号具体为:
光纤内光信号后向散射过程为一系列反射镜的反射过程,这些反射镜是在特定长度ΔL的光纤内随机分布的散射体后向散射的矢量和;
在特定长度ΔL的光纤内有M个随机分布的独立的瑞利散射单元,其偏振均相同,M个散射单元的场矢量为:
其中,rk和φk分别为脉冲跨度内第K段光纤的M个后向散射的幅度和相位的矢量和,ai和Ωi分别是ΔL长度光纤内第i个后向散射的幅度和相位;rk、φk、ai和Ωi为随机变量;
振动点散射回来的瑞利后向散射信号包含非振动光场与振动光场信号,后向散射幅值Es表示为:
ES(t)=E0expi[2πft+Φ0(t)]+E1expi[2πft+Φ1(t)]
为了方便解调过程提取特征分量,在光器件上引入频移Δf:
ES(t)=E0expi[2π(f+Δf)t+Φ0(t)]+E1expi[2π(f+Δf)t+Φ1(t)]
其中,E0为平静时光场强度,E1为振动点光场强度,Φ0、Φ1分别表示平静与振动条件光场相位,即第一项为平静状态未调制光信号光场,第二项为振动点受调制光信号光场,此处将模型简化振动点未受调制瑞利散射单元与受调制瑞利散射单元数量相等;
t时刻的本地光光路光场为:EL0(t)=ELexpi2πft
其中,EL为本地光光场强度。
进一步地,本发明的步骤S3中的拍频信号具体为:
瑞利后向散射信号与本地光源信号进行拍频后为:
其中,E0为平静时光场强度,E1为振动点光场强度,EL为本地光光场强度,EL0为t时刻的本地光光路光场,Φ0、Φ1分别表示平静与振动条件光场相位。
进一步地,本发明的步骤S3中在数字域中实现幅度和相位的解调,解调得到光信号强度的具体方法为:
拍频信号与本地光相干接收后表示为:
其中第一项IL=E0 2+El 2+EL 2表示直流分量,直流分量通过交流耦合滤除;
I∝ELO*[E0*cos(2πΔft+Φ0(t))+E1*cos(2πΔft+Φ1(t))]
对信号平方处理得到:
ELO 2[E0 2+E1 2+2E0E1cos[2πΔft+Φ0(t)]cos[Φ1(t)*cos[2πΔft+Φ1(t)]
低通滤波后:
ELO 2*[E0 2+E1 2+2E0E1cos(Φ0(t)-Φ1(t))]
=ELO(t)2*[E0(t)2+E1(t)2+2E0E1cosΔΦ(t)]
其中,ΔΦ(t)为振动导致光信号相位变化;
得到光电流强度在一定范围内与光场幅值线性对应。
进一步地,本发明的步骤S5中通过信号瀑布图实现对车流量、车重量的实时显示的具体方法为:
当有车辆沿与传感光缆平行铺设的马路行驶时,车辆行驶轨迹在瀑布图上展示,形成一条斜线,斜线斜率倒数即为车辆行驶速度;当有多辆车沿线行驶时,瀑布图上信号反映为多条斜线;
瀑布图单条斜线颜色深度表示单点振动强度,斜线粗细表示振动信号影响范围,通过斜线颜色深度和斜线粗细反映经过车辆的重量大小。
本发明产生的有益效果是:本发明的基于的车流量实时监测系统及方法,将分布式光纤传感技术引入车流量监控领域,可与现有车流量监控体系相结合,可以更有效精确实现对公路沿线车辆行驶情况的24小时实时监测;本发明的信号处理方法中引入对瑞利后向散射信号的处理,将其与本地光源信号进行拍频,对拍频信号中的交流部分提取数字信号,进而解调得到光信号强度,本发明的解调方法解调效果好,能够有效提高计算的准确度;并通过信号瀑布图的形式对车流量等信息进行直观展示,能够提高监测的便利性,尤其针对晚上光线暗、大雾和雨天等视线不明情况导致摄像头等其他监测系统受限制情况对马路上双向车辆进行实时监测,为广大司乘人员出行提供便利,为道路交通管理提供基础。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图2为本发明中测震主机结构图;
图3为车流量信号图;
图4为光纤内光信号后向散射的反射镜反射图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
传感光缆,平行铺设在监测路段的道路两侧,对道路车辆行驶过程中引起的振动信号进行实时监测;
基于的光纤振动监测主机,与传感光缆相连,用于获取光纤振动信号,将光纤振动信号转换为光纤内的瑞利后向散射信号;通过相干检测方法,将后向瑞利散射信号与本地光源信号进行拍频,并提取拍频信号的交流部分信号,对交流部分信号进行采样,得到数字信号,在数字域中实现幅度和相位的解调,解调得到光信号强度;获得瑞利后向散射信号的方式为:
当传感光缆的光纤线路上发生外界扰动时,光纤相应位置的折射率发生变化,导致该处光相位发生变化,引起后向散射光光强发生变化,进而得到从光纤中散射回来的瑞利后向散射信号;
监控端,用于根据光信号强度建立信号瀑布图,通过信号瀑布图实时监控公路上车辆行驶引起的振动信号。
图2为分布式测震主机原理图,包括窄线宽激光器、耦合器、声光调制器、掺铒光纤放大器、光衰减器、环形器、平衡探测器、数字采集卡、工控板等主要部件。
各部件功能分别如下:
激光器作为光源提供窄线宽的连续光,由耦合器1分路为本地光与调制光;
AOM声光调制器将调制光路连续光斩波调制为连续脉冲光;
EDFA掺铒光纤放大器将调制光路脉冲光进行放大;
VOA光衰减器可调节光功率;
调制光路中脉冲光与环形器1端口连接,2端口与传感光纤连接,3端口即瑞利散射信号光;
耦合器2两输入端为本地光与瑞利散射信号光,两者相干后经过平衡探测器输出为电信号传输至采集卡;
采集卡对信号进行数据采集,采集数据送入工控板进一步分析。
图3为现场车辆行驶引起平行公路光缆振动的信号瀑布图,可根据信号瀑布图分析车辆数量、车辆行驶速度等信息。
S1、通过平行铺设在监测路段道路两侧的传感光缆实时监测道路车辆行驶过程引起的光纤振动信号;
S2、通过基于的光纤振动监测主机获取光纤振动信号,将光纤振动信号转换为光纤内的瑞利后向散射信号:当传感光缆的光纤线路上发生外界扰动时,由于弹光效应,光纤相应位置的折射率将发生变化,导致该处光相位发生变化,由于干涉作用,相位的变化将引起后向散射光光强发生变化,进而得到从光纤中散射回来的瑞利后向散射信号;
S3、通过相干检测方法,将后向瑞利散射信号与本地光源信号进行拍频,通过平衡探测器得到拍频信号的交流部分信号,并对交流部分信号进行采样,得到数字信号,在数字域中实现幅度和相位的解调,解调得到光信号强度;
S4、将解调得到的光信号强度发送给监控端,建立信号瀑布图,信号瀑布图为光信号强度在空间、时间域上的实时分布情况;信号瀑布图中横坐标为光纤距离,纵坐标为时间,信号瀑布图中单点颜色反映光信号强度,设定阈值,对超过阈值信号视为振动信号,光信号强度越大与无振动点背景信号颜色对比度越大;通过信号瀑布图实现对车流量、车重量的实时显示。
该方法中进行信号采集处理、解调和瀑布图展示的具体方法为:
1.信号采集处理过程
(1)光纤振动信号转换为光纤内瑞利后向散射信号
如图4所示,光纤内光信号后向散射可理解为一系列反射镜反射,这些反射镜可以看作是在一个特定长度光纤内随机分布的散射体后向散射的矢量和。这些散射体的后向散射光的电场具有随机的相位和幅度。因此,这些散射的矢量和将会随机落在一个复平面内。
为了得到这些反射镜的反射率和相位的统计值,我们假设在ΔL内有M个随机分布的独立的瑞利散射元。而且这些散射的偏振都是相同的,M个散射的场矢量和可以表示为:
其中,rk和φk分别为脉冲跨度内第K段光纤的M个后向散射的幅度和相位的矢量和,ai和Ωi分别是ΔL长度光纤内第i个后向散射的幅度和相位。rk、φk、ai和Ωi都假定为随机变量。为了使分析简单,假设ai和Ωi都是独立并对于所有的i分布都是一样的,而且Ωi是均匀分布在区间(-π,π)上的。那么,反射镜的反射率r和相位φ的概率密度函数如下:
其中,σ2近似等于0.5×10-7。随机变量r和φ分别满足瑞利分布和均匀分布。当脉冲宽度为NΔL,从距离mΔL处的后向散射光的干涉场可以表示为:
其中,Pk、rk和φk分别表示第k个反射镜的偏振、反射率和相位,它们都是随机变量。θk则表示入侵扰动引入到第k个反射镜中的相位变化。
一维脉冲响应模型:假设在t=0时刻,我们将一个相干光脉冲入射到光纤中,其脉冲宽度为W,频率为v,在光纤的入射端处后向瑞利散射光波e(t)可以表示为:
其中,ai和τi分别表示第i个散射波的振幅和往返时间,N表示散射点的总数,α表示光纤的衰减常数,c表示真空中的光速,nf表示光纤的折射率。
其中:
第i个散射点往返时间和离光纤入射段之间距离zi可以表示成τi=2nfzi/c。因此后向散射光功率可以表示成如下:
p(t)=|e(t)|2=p1(t)+p2(t)
其中,φij=2πv(τi-τj),表示第i个和第j个散射点的相位差。上式中p1(t)代表N个散射点各自散射功率之和,该项不随着光纤温度或应变的变化,也不随着激光频率的改变。而p2(t)表示了不同散射点之间干涉的结果,使得瑞利散射回波呈现出干涉衰落波形。在p2(t)表达式中有一项cos(φij),其中相邻散射点之间的相位差φij正比于激光器的频率v,折射率nf,散射点间距sij=zi-zj,其关系是可以表示为φij=4πvnfsij/c。因此干涉项p2(t)是v,nf和sij的函数,而nf和sij取决于光纤温度和应变的变化。
当光纤线路上发生外界扰动时,由于弹光效应,光纤相应位置的折射率将发生变化,这将导致该处光相位发生变化,由于干涉作用,相位的变化将引起后向散射光光强发生变化。
振动点散射回来的瑞利后向散射信号包含非振动光场与振动光场信号,后向散射幅值Es表示为:
ES(t)=E0expi[2πft+Φ0(t)]+E1expi[2πft+Φ1(t)]
从光纤中散射回来的瑞利信号(传感光信号)表示为:
为了方便解调过程提取特征分量,在光器件上引入频移Δf:
ES(t)=E0expi[2π(f+Δf)t+Φ0(t)]+E1expi[2π(f+Δf)t+Φ1(t)]
其中,E0为平静时光场强度,E1为振动点光场强度,Φ0、Φ1分别表示平静与振动条件光场相位,即第一项为平静状态未调制光信号光场,第二项为振动点受调制光信号光场,此处将模型简化振动点未受调制瑞利散射单元与受调制瑞利散射单元数量相等;
t时刻的本地光光路光场为:EL0(t)=ELexpi2πft
其中,EL为本地光光场强度。
(2)相干检测方法对原始光信号解调处理
通过相干检测方法,将后向瑞利散射光与本地光进行拍频,拍频之后信号经过一个双平衡探测器接收,得到拍频信号的交流部分。接着利用数据采集卡对得到的交流信号进行采样,获得数字信号就可以在数字域中实现幅度和相位的解调。具体过程如下:
瑞利散射信号为:
与本地光相干接收后可以表示为:
针对解调光信号第一项,通过平衡探测器BPD进行交流耦合滤除,第二项则通过带通滤除,得到交流拍频信号(振动信号)。
其中第一项IL=E0 2+El 2+EL 2表示直流分量,直流分量通过交流耦合滤除;
I∝ELO*[E0*cos(2πΔft+Φ0(t))+E1*cos(2πΔft+Φ1(t))]
对信号平方处理得到:
ELO 2[E0 2+E1 2+2E0E1cos[2πΔft+Φ0(t)]cos[Φ1(t)*cos[2πΔft+Φ1(t)]
低通滤波后:
ELO 2*[E0 2+E1 2+2E0E1cos(Φ0(t)-Φ1(t))]
=ELO(t)2*[E0(t)2+E1(t)2+2E0E1cosΔΦ(t)]
其中,ΔΦ(t)为振动导致光信号相位变化;
得到光电流强度在一定范围内与光场幅值线性对应。
在直接探测系统中,第3项光谱积分为零,信号出自第2项。散射系数很小,约为10-7。而就大的多,约为10-3~10-4量级。因此,采用相干探测方法具有较强信号功率。拍频之后信号经过一个双平衡探测器接收,得到拍频信号的交流部分表示为
上式中相位因子φ(t)=2βz+δ。光信号强度在一定范围内与振动强度线性相关。
2.信号瀑布图
信号瀑布图实际为解调光信号强度(振动信号强度)在空间、时间域上的实时分布情况。
信号瀑布图横坐标为光纤距离(空间域),将整段光纤分区为多个小段采集数据,实现高定位精度(考虑到采样数据量限制,一般为50km,10m空间分辨率)。
信号瀑布图纵坐标为时间,瀑布图下降速度由采样帧数决定。
信号瀑布图中单点颜色反映信号强度,设定阈值,对超过阈值信号视为振动信号,信号强度越大与无振动点背景信号颜色对比度越大。
当有车辆沿与传感光缆平行铺设的马路行驶,车辆行驶轨迹可在瀑布图上展示,形成一条斜线,斜线斜率倒数即为车辆行驶速度;当有多辆车沿线行驶,瀑布图上信号反映为多条斜线;瀑布图单条斜线颜色深度(单点振动强度)以及斜线粗细(振动信号影响范围)反映经过车辆的重量大小。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (10)
传感光缆,平行铺设在监测路段的道路两侧,对道路车辆行驶过程中引起的振动信号进行实时监测;
基于的光纤振动监测主机,与传感光缆相连,用于获取光纤振动信号,将光纤振动信号转换为光纤内的瑞利后向散射信号;通过相干检测方法,将后向瑞利散射信号与本地光源信号进行拍频,并提取拍频信号的交流部分信号,对交流部分信号进行采样,得到数字信号,在数字域中实现幅度和相位的解调,解调得到光信号强度;
监控端,用于根据光信号强度建立信号瀑布图,通过信号瀑布图实时监控公路上车辆行驶引起的振动信号。
3.根据权利要求1所述的基于的车流量实时监测系统,其特征在于,基于的光纤振动监测主机包括:窄线宽激光器、第一耦合器、第二耦合器、声光调制器、掺铒光纤放大器、光衰减器、环形器、平衡探测器、数字采集卡、工控板;其中:
激光器作为光源提供窄线宽的连续光,由第一耦合器分路为本地光与调制光;
声光调制器将调制光路连续光斩波调制为连续脉冲光;
掺铒光纤放大器将调制光路脉冲光进行放大;
光衰减器用于调节光功率;
调制光路中脉冲光与环形器第一端口连接,第二端口与传感光缆连接,第三端口即瑞利散射信号光;
第二耦合器两输入端为本地光与瑞利散射信号光,两者相干后经过平衡探测器输出为电信号传输至数字采集卡;
数字采集卡对信号进行数据采集,采集数据送入工控板进一步分析。
S1、通过平行铺设在监测路段道路两侧的传感光缆实时监测道路车辆行驶过程引起的光纤振动信号;
S2、通过基于的光纤振动监测主机获取光纤振动信号,将光纤振动信号转换为光纤内的瑞利后向散射信号:当传感光缆的光纤线路上发生外界扰动时,光纤相应位置的折射率发生变化,导致该处光相位发生变化,引起后向散射光光强发生变化,进而得到从光纤中散射回来的瑞利后向散射信号;
S3、通过相干检测方法,将后向瑞利散射信号与本地光源信号进行拍频,通过平衡探测器得到拍频信号的交流部分信号,并对交流部分信号进行采样,得到数字信号,在数字域中实现幅度和相位的解调,解调得到光信号强度;
S4、将解调得到的光信号强度发送给监控端,建立信号瀑布图,信号瀑布图为光信号强度在空间、时间域上的实时分布情况;信号瀑布图中横坐标为光纤距离,纵坐标为时间,信号瀑布图中单点颜色反映光信号强度,设定阈值,对超过阈值信号视为振动信号,光信号强度越大与无振动点背景信号颜色对比度越大;通过信号瀑布图实现对车流量、车重量的实时显示。
光纤内光信号后向散射过程为一系列反射镜的反射过程,这些反射镜是在特定长度ΔL的光纤内随机分布的散射体后向散射的矢量和;
在特定长度ΔL的光纤内有M个随机分布的独立的瑞利散射单元,其偏振均相同,M个散射单元的场矢量为:
其中,rk和φk分别为脉冲跨度内第K段光纤的M个后向散射的幅度和相位的矢量和,ai和Ωi分别是ΔL长度光纤内第i个后向散射的幅度和相位;rk、φk、ai和Ωi为随机变量;
振动点散射回来的瑞利后向散射信号包含非振动光场与振动光场信号,后向散射幅值Es表示为:
ES(t)=E0expi[2πft+Φ0(t)]+E1expi[2πft+Φ1(t)]
为了方便解调过程提取特征分量,在光器件上引入频移Δf:
Es(t)=E0expi[2π(f+Δf)t+Φ0(t)]+E1expi[2π(f+Δf)t+Φ1(t)]
其中,E0为平静时光场强度,E1为振动点光场强度,Φ0、Φ1分别表示平静与振动条件光场相位,即第一项为平静状态未调制光信号光场,第二项为振动点受调制光信号光场,此处将模型简化振动点未受调制瑞利散射单元与受调制瑞利散射单元数量相等;
t时刻的本地光光路光场为:EL0(t)=ELexpi2πft
其中,EL为本地光光场强度。
拍频信号与本地光相干接收后表示为:
其中第一项IL=E0 2+El 2+EL 2表示直流分量,直流分量通过交流耦合滤除;
I∝ELO*[E0*cos(2πΔft+Φ0(t))+E1*cos(2πΔft+Φ1(t))]
对信号平方处理得到:
ELO 2[E0 2+E1 2+2E0E1cos[2πΔft+Φ0(t)]cos[Φ1(t)*cos[2πΔft+Φ1(t)]
低通滤波后:
ELO 2*[E0 2+E1 2+2E0E1cos(Φ0(t)-Φ1(t))]
=ELO(t)2*[E0(t)2+E1(t)2+2E0E1cosΔΦ(t)]
其中,ΔΦ(t)为振动导致光信号相位变化;
得到光电流强度在一定范围内与光场幅值线性对应。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911345757.3A CN111147133B (zh) | 2019-12-24 | 2019-12-24 | 一种基于φ-OTDR的车流量实时监测系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911345757.3A CN111147133B (zh) | 2019-12-24 | 2019-12-24 | 一种基于φ-OTDR的车流量实时监测系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111147133A true CN111147133A (zh) | 2020-05-12 |
CN111147133B CN111147133B (zh) | 2021-09-14 |
Family
ID=70519645
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911345757.3A Active CN111147133B (zh) | 2019-12-24 | 2019-12-24 | 一种基于φ-OTDR的车流量实时监测系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111147133B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115063984A (zh) * | 2022-08-05 | 2022-09-16 | 武汉理工大学 | 基于光纤光栅阵列传感网的公路车辆识别和定位方法 |
WO2022221596A1 (en) * | 2021-04-14 | 2022-10-20 | Nec Laboratories America, Inc. | Mapping using optical fiber sensing |
CN115410403A (zh) * | 2022-04-19 | 2022-11-29 | 北京见合八方科技发展有限公司 | 基于无源感知的公路车辆定位追踪方法、装置及可读介质 |
WO2023004182A1 (en) * | 2021-07-23 | 2023-01-26 | Nec Laboratories America, Inc. | Location determination of deployed fiber cables using distributed fiber optic sensing |
WO2024051501A1 (zh) * | 2022-09-05 | 2024-03-14 | 国家石油天然气管网集团有限公司 | 基于光纤的多目标定位方法、处理器、装置及存储介质 |
Citations (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005338023A (ja) * | 2004-05-31 | 2005-12-08 | Taisei Kiso Sekkei Kk | 変状監視システム |
CN101122555A (zh) * | 2007-09-12 | 2008-02-13 | 上海理工大学 | 基于后向光子相关光谱的高浓度超细颗粒测量装置及方法 |
CN102349834A (zh) * | 2011-06-20 | 2012-02-15 | 深圳职业技术学院 | 人体血糖浓度无创检测方法及其系统 |
CN203377882U (zh) * | 2013-06-26 | 2014-01-01 | 南车青岛四方机车车辆股份有限公司 | 轨道车辆光纤网络诊断装置及用该装置的轨道车辆 |
CN103634046A (zh) * | 2013-06-26 | 2014-03-12 | 南车青岛四方机车车辆股份有限公司 | 轨道车辆光纤网络诊断装置和方法及用该装置的轨道车辆 |
CN103926588A (zh) * | 2014-04-28 | 2014-07-16 | 无锡成电光纤传感科技有限公司 | 一种基于φ-otdr的轨道车辆定位和测速系统 |
CN104700624A (zh) * | 2015-03-16 | 2015-06-10 | 电子科技大学 | 基于相敏光时域反射仪的车流量在线监测系统及监测方法 |
CN205038802U (zh) * | 2015-10-12 | 2016-02-17 | 重庆国文科技有限公司 | 高速公路突发事件监测系统 |
CN105865613A (zh) * | 2016-06-03 | 2016-08-17 | 哈尔滨工业大学深圳研究生院 | 海洋立体监测中的水下光学检测与成像传感方法及系统 |
CN106597423A (zh) * | 2016-12-07 | 2017-04-26 | 南京富岛信息工程有限公司 | 一种基于相位敏感光时域反射计的列车行驶轨迹检测方法 |
CN106768281A (zh) * | 2017-04-11 | 2017-05-31 | 光子瑞利科技(北京)有限公司 | 相位敏感型φ‑OTDR的分布式光纤听音器 |
CN107101658A (zh) * | 2017-05-24 | 2017-08-29 | 上海大学 | 相位敏感光时域反射分布式光纤传感系统快速定位方法 |
CN107256635A (zh) * | 2017-07-14 | 2017-10-17 | 浙江大学 | 一种智能交通中基于分布式光纤传感的车辆识别方法 |
CN107505041A (zh) * | 2017-07-31 | 2017-12-22 | 南京发艾博光电科技有限公司 | 一种基于相位敏感光时域反射计的相位解调装置和方法 |
CN107591002A (zh) * | 2017-09-21 | 2018-01-16 | 电子科技大学 | 一种基于分布式光纤的高速公路交通参数实时估计方法 |
CN207337666U (zh) * | 2017-08-29 | 2018-05-08 | 招商局重庆交通科研设计院有限公司 | 交通事故定位系统 |
CN108287017A (zh) * | 2017-12-30 | 2018-07-17 | 武汉理工光科股份有限公司 | 基于φ-OTDR分布式光纤传感的高速公路护栏碰撞事故报警系统 |
CN207894512U (zh) * | 2017-12-30 | 2018-09-21 | 武汉理工光科股份有限公司 | 基于φ-OTDR的振动传感装置 |
CN108645500A (zh) * | 2018-06-11 | 2018-10-12 | 桂林电子科技大学 | 一种车辆振动传感增敏的方法及装置 |
CN108875684A (zh) * | 2018-06-29 | 2018-11-23 | 电子科技大学 | 基于光纤感测时频图处理的列车运行状态参数估计方法 |
CN109039440A (zh) * | 2018-09-14 | 2018-12-18 | 长飞光纤光缆股份有限公司 | 一种高效otdr测试系统及方法 |
CN109405956A (zh) * | 2018-12-11 | 2019-03-01 | 珠海任驰光电科技有限公司 | 一种分布式光纤振动传感系统解调方法及装置 |
CN109752720A (zh) * | 2018-10-09 | 2019-05-14 | 武汉理工光科股份有限公司 | 基于分布式振动传感器振动信息定位方法及系统 |
CN110031082A (zh) * | 2019-05-10 | 2019-07-19 | 山东省科学院激光研究所 | 分布式光纤振动监测系统事件定位方法、装置及应用 |
CN110466568A (zh) * | 2019-08-26 | 2019-11-19 | 南京派光智慧感知信息技术有限公司 | 一种用于车站的列车运行监控系统 |
-
2019
- 2019-12-24 CN CN201911345757.3A patent/CN111147133B/zh active Active
Patent Citations (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005338023A (ja) * | 2004-05-31 | 2005-12-08 | Taisei Kiso Sekkei Kk | 変状監視システム |
CN101122555A (zh) * | 2007-09-12 | 2008-02-13 | 上海理工大学 | 基于后向光子相关光谱的高浓度超细颗粒测量装置及方法 |
CN102349834A (zh) * | 2011-06-20 | 2012-02-15 | 深圳职业技术学院 | 人体血糖浓度无创检测方法及其系统 |
CN203377882U (zh) * | 2013-06-26 | 2014-01-01 | 南车青岛四方机车车辆股份有限公司 | 轨道车辆光纤网络诊断装置及用该装置的轨道车辆 |
CN103634046A (zh) * | 2013-06-26 | 2014-03-12 | 南车青岛四方机车车辆股份有限公司 | 轨道车辆光纤网络诊断装置和方法及用该装置的轨道车辆 |
CN103926588A (zh) * | 2014-04-28 | 2014-07-16 | 无锡成电光纤传感科技有限公司 | 一种基于φ-otdr的轨道车辆定位和测速系统 |
CN104700624A (zh) * | 2015-03-16 | 2015-06-10 | 电子科技大学 | 基于相敏光时域反射仪的车流量在线监测系统及监测方法 |
CN205038802U (zh) * | 2015-10-12 | 2016-02-17 | 重庆国文科技有限公司 | 高速公路突发事件监测系统 |
CN105865613A (zh) * | 2016-06-03 | 2016-08-17 | 哈尔滨工业大学深圳研究生院 | 海洋立体监测中的水下光学检测与成像传感方法及系统 |
CN106597423A (zh) * | 2016-12-07 | 2017-04-26 | 南京富岛信息工程有限公司 | 一种基于相位敏感光时域反射计的列车行驶轨迹检测方法 |
CN106768281A (zh) * | 2017-04-11 | 2017-05-31 | 光子瑞利科技(北京)有限公司 | 相位敏感型φ‑OTDR的分布式光纤听音器 |
CN107101658A (zh) * | 2017-05-24 | 2017-08-29 | 上海大学 | 相位敏感光时域反射分布式光纤传感系统快速定位方法 |
CN107256635A (zh) * | 2017-07-14 | 2017-10-17 | 浙江大学 | 一种智能交通中基于分布式光纤传感的车辆识别方法 |
CN107505041A (zh) * | 2017-07-31 | 2017-12-22 | 南京发艾博光电科技有限公司 | 一种基于相位敏感光时域反射计的相位解调装置和方法 |
CN207337666U (zh) * | 2017-08-29 | 2018-05-08 | 招商局重庆交通科研设计院有限公司 | 交通事故定位系统 |
CN107591002A (zh) * | 2017-09-21 | 2018-01-16 | 电子科技大学 | 一种基于分布式光纤的高速公路交通参数实时估计方法 |
CN108287017A (zh) * | 2017-12-30 | 2018-07-17 | 武汉理工光科股份有限公司 | 基于φ-OTDR分布式光纤传感的高速公路护栏碰撞事故报警系统 |
CN207894512U (zh) * | 2017-12-30 | 2018-09-21 | 武汉理工光科股份有限公司 | 基于φ-OTDR的振动传感装置 |
CN108645500A (zh) * | 2018-06-11 | 2018-10-12 | 桂林电子科技大学 | 一种车辆振动传感增敏的方法及装置 |
CN108875684A (zh) * | 2018-06-29 | 2018-11-23 | 电子科技大学 | 基于光纤感测时频图处理的列车运行状态参数估计方法 |
CN109039440A (zh) * | 2018-09-14 | 2018-12-18 | 长飞光纤光缆股份有限公司 | 一种高效otdr测试系统及方法 |
CN109752720A (zh) * | 2018-10-09 | 2019-05-14 | 武汉理工光科股份有限公司 | 基于分布式振动传感器振动信息定位方法及系统 |
CN109405956A (zh) * | 2018-12-11 | 2019-03-01 | 珠海任驰光电科技有限公司 | 一种分布式光纤振动传感系统解调方法及装置 |
CN110031082A (zh) * | 2019-05-10 | 2019-07-19 | 山东省科学院激光研究所 | 分布式光纤振动监测系统事件定位方法、装置及应用 |
CN110466568A (zh) * | 2019-08-26 | 2019-11-19 | 南京派光智慧感知信息技术有限公司 | 一种用于车站的列车运行监控系统 |
Non-Patent Citations (2)
Title |
---|
SHENG LIANG; XIAOTING ZHAO; RUIJUN LIU; XIAODONG ZHANG; LIANG WA: "Fiber-Optic Auditory Nerve of Ground in the Suburb: For Traffic Flow Monitoring", 《IEEE ACCESS》 * |
冯亚非 郑欢 桑卫兵 葛辉良: "基于瑞利散射的分布式正弦扰动信号检测", 《声学与电子工程》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2022221596A1 (en) * | 2021-04-14 | 2022-10-20 | Nec Laboratories America, Inc. | Mapping using optical fiber sensing |
WO2023004182A1 (en) * | 2021-07-23 | 2023-01-26 | Nec Laboratories America, Inc. | Location determination of deployed fiber cables using distributed fiber optic sensing |
CN115410403A (zh) * | 2022-04-19 | 2022-11-29 | 北京见合八方科技发展有限公司 | 基于无源感知的公路车辆定位追踪方法、装置及可读介质 |
CN115410403B (zh) * | 2022-04-19 | 2023-11-10 | 北京见合八方科技发展有限公司 | 基于无源感知的公路车辆定位追踪方法、装置及可读介质 |
CN115063984A (zh) * | 2022-08-05 | 2022-09-16 | 武汉理工大学 | 基于光纤光栅阵列传感网的公路车辆识别和定位方法 |
CN115063984B (zh) * | 2022-08-05 | 2022-11-15 | 武汉理工大学 | 基于光纤光栅阵列传感网的公路车辆识别和定位方法 |
WO2024051501A1 (zh) * | 2022-09-05 | 2024-03-14 | 国家石油天然气管网集团有限公司 | 基于光纤的多目标定位方法、处理器、装置及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN111147133B (zh) | 2021-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111147133B (zh) | 一种基于φ-OTDR的车流量实时监测系统及方法 | |
CN104819770B (zh) | 基于外差检测和相位解调的相位光时域反射装置及方法 | |
CN105806465B (zh) | 一种基于固定反射点的新型φ-otdr探测装置及其探测方法 | |
CN105509868B (zh) | 相位敏感光时域反射光纤分布式传感系统相位计算方法 | |
EP3098651B1 (en) | Light pulse compression reflectometer | |
CN107917738A (zh) | 一种同时测量温度、应变和振动的分布式光纤传感系统 | |
CN102645268A (zh) | 光频分复用相位敏感光时域反射计 | |
CN106019313B (zh) | 基于偏振双边缘的单像素探测测风激光雷达 | |
Ghafoori-Shiraz et al. | Fault location in optical fibers using optical frequency domain reflectometry | |
CN106940444A (zh) | 基于微波差分增益的相干多普勒测风激光雷达 | |
CN111854921A (zh) | 一种分布式光纤减速带振动预警系统及方法 | |
CN107505041A (zh) | 一种基于相位敏感光时域反射计的相位解调装置和方法 | |
CN113447110B (zh) | 一种分布式光纤振动传感系统及其相位载波解调方法 | |
CN105866795A (zh) | 基于全光纤扫描f-p干涉仪的大动态风场探测激光雷达 | |
CN207036249U (zh) | 一种高灵敏度的分布式光纤振动传感系统 | |
CN105067103A (zh) | 基于光频域反射计的振动检测装置及其方法 | |
CN108387909A (zh) | 基于激光雷达网的区域环境监测系统 | |
US11735038B2 (en) | Non-blind area real-time monitoring and alarming system for accident on freeway | |
CN109120337A (zh) | 一种少模时域反射仪 | |
CN103438982A (zh) | 一种基于布里渊分布式光纤传感的晃动监测系统 | |
CN113189450B (zh) | 基于微气象监测和svm算法的线路覆冰监测系统 | |
CN205120239U (zh) | 基于光频域反射计的振动检测装置 | |
CN109541636A (zh) | 一种无盲区高距离分辨率激光雷达测风系统及方法 | |
CN107545684B (zh) | 一种分布式光纤振动传感系统中的多点扰动事件报警方法 | |
CN209605977U (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 |