CN207611485U - An Embedded Road Vehicle Type Recognition System - Google Patents
An Embedded Road Vehicle Type Recognition System Download PDFInfo
- Publication number
- CN207611485U CN207611485U CN201721385512.XU CN201721385512U CN207611485U CN 207611485 U CN207611485 U CN 207611485U CN 201721385512 U CN201721385512 U CN 201721385512U CN 207611485 U CN207611485 U CN 207611485U
- Authority
- CN
- China
- Prior art keywords
- optical fiber
- vibration
- vibrating sensing
- vehicle type
- road vehicle
- 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
Links
- 239000013307 optical fiber Substances 0.000 abstract description 115
- 238000000034 method Methods 0.000 abstract description 34
- 239000000835 fiber Substances 0.000 abstract description 19
- 238000005516 engineering process Methods 0.000 abstract description 5
- 230000003287 optical effect Effects 0.000 abstract description 5
- 230000008447 perception Effects 0.000 abstract description 2
- 230000001953 sensory effect Effects 0.000 abstract 1
- 238000004458 analytical method Methods 0.000 description 23
- 238000012549 training Methods 0.000 description 16
- 230000007704 transition Effects 0.000 description 15
- 238000005259 measurement Methods 0.000 description 11
- 238000001845 vibrational spectrum Methods 0.000 description 9
- 230000004044 response Effects 0.000 description 8
- 239000004567 concrete Substances 0.000 description 7
- 238000012706 support-vector machine Methods 0.000 description 7
- 230000006870 function Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 5
- 239000004568 cement Substances 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000000354 decomposition reaction Methods 0.000 description 4
- 230000005284 excitation Effects 0.000 description 4
- 238000000605 extraction Methods 0.000 description 4
- 239000002184 metal Substances 0.000 description 4
- 238000004804 winding Methods 0.000 description 4
- 229910000831 Steel Inorganic materials 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000012850 discrimination method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000007774 longterm Effects 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 239000010959 steel Substances 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 239000011384 asphalt concrete Substances 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 230000006698 induction Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000006855 networking Effects 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 238000011056 performance test Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
Description
技术领域technical field
本实用新型涉及道路工程领域,特别是涉及一种埋入式道路车辆类型识别系统。The utility model relates to the field of road engineering, in particular to an embedded road vehicle type identification system.
背景技术Background technique
道路行驶车辆的自动检测和识别是智能交通系统中重要组成部分。在交通规划阶段,对调查交通流进行车型识别可为交通统计提供更可靠的依据,也有利于制定更科学合理的交通规划;在道路交通监测与控制领域,收费站、停车场等方面也对车型识别有着大量的应用需求;此外车型识别还可为交通事件的处理、车辆的跟踪提供证据与帮助。The automatic detection and identification of road vehicles is an important part of the intelligent transportation system. In the stage of traffic planning, the identification of vehicle types for surveyed traffic flow can provide a more reliable basis for traffic statistics, and is also conducive to the formulation of more scientific and reasonable traffic planning; in the field of road traffic monitoring and control, toll stations, parking lots, etc. Vehicle type identification has a large number of application requirements; in addition, vehicle type identification can also provide evidence and help for the handling of traffic incidents and vehicle tracking.
为构建一个高效、便捷的车辆识别系统,必须保证较高的识别精度、耐久性以及安装的便捷性。目前针对自动化车型识别技术,主要集中在三类方法:基于图像的识别,基于声音特征的识别和基于埋入式传感器的识别。图像识别采用普通摄像头或红外摄像头对车辆轮廓进行采集和识别,能够较为准确地判定车辆类型,但在夜间和低可视度的条件下效果不佳;声音特征的识别采用麦克风等设备进行数据采集和识别,但其受环境噪声影响较大,无法适用于多车道、大流量状况下的车型识别;基于埋入式传感器的方法是目前主流的车型识别方法,包括地感线圈、压电传感器、地磁传感器以及加速度传感器,此类传感器中,地感线圈较为廉价,但是识别精度有限,另几种传感器则存在组网不便、寿命有限以及价格昂贵等缺陷,可靠性不佳。In order to build an efficient and convenient vehicle identification system, it is necessary to ensure high identification accuracy, durability and ease of installation. At present, for automatic vehicle recognition technology, there are mainly three types of methods: image-based recognition, sound feature-based recognition and embedded sensor-based recognition. Image recognition uses ordinary cameras or infrared cameras to collect and recognize vehicle outlines, which can accurately determine the type of vehicle, but the effect is not good at night and in low-visibility conditions; sound feature recognition uses microphones and other equipment for data collection and recognition, but it is greatly affected by environmental noise and cannot be applied to vehicle identification under multi-lane and high-flow conditions; the method based on embedded sensors is the current mainstream vehicle identification method, including ground induction coils, piezoelectric sensors, Geomagnetic sensors and acceleration sensors. Among these sensors, the ground induction coil is relatively cheap, but the recognition accuracy is limited. The other types of sensors have defects such as inconvenient networking, limited life, and high price, and their reliability is not good.
因此,为便于实时交通管理与控制,需要一种可靠性更高、组网更便捷的道路车辆类型识别系统和方法。Therefore, in order to facilitate real-time traffic management and control, a road vehicle type identification system and method with higher reliability and more convenient networking are needed.
实用新型内容Utility model content
鉴于以上所述现有技术的缺点,本实用新型的目的在于提供一种埋入式道路车辆类型识别系统和方法,用于解决现有技术中的问题。In view of the above-mentioned shortcomings of the prior art, the purpose of this utility model is to provide an embedded road vehicle type identification system and method for solving the problems in the prior art.
为实现上述目的及其他相关目的,本实用新型第一方面提供一种埋入式道路车辆类型识别系统,包括路面结构本体,所述路面结构本体中设有振动传感光纤,所述振动传感光纤包括一个或多个振动传感光纤传感段,各振动传感光纤传感段之间通过振动传感光纤过渡段连接,还包括振动光纤分析装置,所述振动光纤分析装置与振动传感光纤通过光纤引出线相连接。In order to achieve the above purpose and other related purposes, the first aspect of the utility model provides an embedded road vehicle type identification system, including a road structure body, a vibration sensing optical fiber is arranged in the road structure body, and the vibration sensor The optical fiber includes one or more vibration sensing optical fiber sensing sections, each vibration sensing optical fiber sensing section is connected by a vibration sensing optical fiber transition section, and also includes a vibration optical fiber analysis device, the vibration optical fiber analysis device is connected with the vibration sensing The optical fibers are connected by optical fiber pigtails.
在本实用新型一些实施方式中,所述振动传感光纤为单模光纤。In some embodiments of the present utility model, the vibration sensing optical fiber is a single-mode optical fiber.
在本实用新型一些实施方式中,所述振动传感光纤为金属铠装光纤。In some embodiments of the present utility model, the vibration sensing optical fiber is a metal armored optical fiber.
在本实用新型一些实施方式中,所述振动传感光纤的直径为2~5mm。In some embodiments of the present utility model, the diameter of the vibration sensing optical fiber is 2-5 mm.
在本实用新型一些实施方式中,所述振动传感光纤的长期允许拉伸力≥600N,短期允许拉伸力≥1500N。In some embodiments of the present utility model, the long-term allowable tensile force of the vibration sensing optical fiber is ≥600N, and the short-term allowable tensile force is ≥1500N.
在本实用新型一些实施方式中,所述振动传感光纤的衰减≤0.2db/Km。In some embodiments of the present invention, the attenuation of the vibration sensing optical fiber is ≤0.2db/Km.
在本实用新型一些实施方式中,所述振动传感光纤包括多个振动传感光纤传感段,各振动传感光纤传感段之间的间距≥0.2m且≤0.5m。In some embodiments of the present utility model, the vibration sensing optical fiber includes a plurality of vibration sensing optical fiber sensing sections, and the distance between each vibration sensing optical fiber sensing section is ≥0.2m and ≤0.5m.
在本实用新型一些实施方式中,所述振动传感光纤传感段为螺旋形,螺旋形的长轴竖直分布,振动传感光纤传感段长轴方向的高度为5~20mm,振动传感光纤传感段的直径为250-350mm,每个振动传感光纤传感段中光纤的长度≥4m且≤6m。In some embodiments of the present invention, the sensing section of the vibration sensing optical fiber is spiral, and the long axis of the spiral is distributed vertically. The diameter of the sensing optical fiber sensing section is 250-350mm, and the length of the optical fiber in each vibration sensing optical fiber sensing section is ≥4m and ≤6m.
在本实用新型一些实施方式中,所述路面结构本体中设有多排振动传感光纤传感段。In some embodiments of the present utility model, multiple rows of vibration sensing optical fiber sensing sections are arranged in the road surface structure body.
在本实用新型一些实施方式中,所述振动传感光纤传感段在车道宽度方向均匀分布,覆盖率为2~3个/米。In some embodiments of the present utility model, the vibration sensing optical fiber sensing segments are evenly distributed in the width direction of the lane, and the coverage rate is 2-3 per meter.
在本实用新型一些实施方式中,所述路面结构本体中按车道的长度方向布设有2~3排振动传感光纤传感段,各排振动传感光纤传感段之间的间距为3~10m。In some embodiments of the present utility model, 2 to 3 rows of vibration-sensing optical fiber sensing sections are arranged in the road surface structure body according to the length direction of the lane, and the distance between each row of vibration-sensing optical fiber sensing sections is 3-3. 10m.
在本实用新型一些实施方式中,所述路面结构本体为水泥混凝土铺面、沥青混凝土铺面或复合型铺面结构中的一种或多种的组合;In some embodiments of the present utility model, the pavement structure body is a combination of one or more of cement concrete pavement, asphalt concrete pavement or composite pavement structure;
在本实用新型一些实施方式中,所述路面结构本体中设有多排振动传感光纤传感段,至少部分的振动传感光纤传感段依次串联。In some embodiments of the present utility model, multiple rows of vibration-sensing optical fiber sensing sections are arranged in the road surface structure body, and at least part of the vibration-sensing optical fiber sensing sections are connected in series in sequence.
本实用新型第二方面提供一种道路车辆类型识别方法,使用所述的埋入式道路车辆类型识别系统埋入式道路车辆类型识别系统,包括如下步骤:The second aspect of the present invention provides a road vehicle type identification method, using the embedded road vehicle type identification system The embedded road vehicle type identification system includes the following steps:
1)采用所述的埋入式车型识别系统,通过振动光纤分析装置采集车辆经过时的路面结构振动数据;1) Using the embedded vehicle type identification system, the vibration data of the road surface structure when the vehicle passes through the vibration optical fiber analysis device is collected;
2)对各振动数据组进行特征提取;2) Carry out feature extraction to each vibration data group;
3)依据上述车辆的特征信息,对行驶车辆进行分类。3) According to the characteristic information of the above-mentioned vehicles, classify the driving vehicles.
在本实用新型一些实施方式中,所述步骤2)中,所述特征选自行驶车辆的移动速率、轴型、轮型、轴距及振动频谱分布特征中的一种或多种的组合。In some implementations of the present utility model, in the step 2), the features are selected from one or more combinations of moving speed, axle type, wheel type, wheelbase and vibration spectrum distribution characteristics of the traveling vehicle.
在本实用新型一些实施方式中,所述步骤3)中,用于对行驶车辆进行分类的特征信息为轴型和/或振动频谱分布特性。In some embodiments of the present utility model, in the step 3), the feature information used to classify the driving vehicle is the shaft type and/or the vibration spectrum distribution characteristics.
在本实用新型一些实施方式中,所述步骤2)中,对各振动数据组进行特征提取的方法具体包括如下步骤:In some embodiments of the present invention, in the step 2), the method for feature extraction of each vibration data group specifically includes the following steps:
a)将同一横断面内测点所观测的原始振动数据编入同一振动数据组,并进行经验模态分解(EMD),获得多个经EMD处理后的振动分量;a) Compile the original vibration data observed at the measuring points in the same cross-section into the same vibration data group, and perform empirical mode decomposition (EMD) to obtain multiple vibration components after EMD processing;
b)叠加特定阶数的振动分量,获得叠加后的振动曲线fr(t)(x轴为时间,y轴为振动强度):b) Superimpose vibration components of a specific order to obtain the superimposed vibration curve fr(t) (x-axis is time, y-axis is vibration intensity):
其中IMFi为EMD处理后第i阶的振动分量,n和m为特定阶数的最低阶数和最高阶数。Among them, IMF i is the vibration component of the i-th order after EMD treatment, and n and m are the lowest and highest orders of specific orders.
c)计算单位时间(τ)内振动数据的平方和作为该段时间内的短时能量,依此得到振动短时能量的时程曲线E(T)(x轴为时间,y轴为振动短时能量):c) Calculate the sum of the squares of the vibration data per unit time (τ) as the short-term energy within this period, and thus obtain the time-history curve E(T) of the short-term vibration energy (the x-axis is time, and the y-axis is the short-term vibration energy time energy):
其中单位时间τ即为短时能量分析时间帧长度,初始值可以设置为0.5s;The unit time τ is the length of the short-term energy analysis time frame, and the initial value can be set to 0.5s;
d)设定车辆振动能量判定阈值,用于判定车辆是否经过的短时能量时程曲线中满足该阈值要求的部分视为因行驶车辆激励产生的振动时段[T1,T2]。随后从叠加的振动曲线fr(t)截取出该时段内的振动曲线fr’(t);d) Set the vehicle vibration energy judgment threshold, and the part of the short-term energy time-history curve used to determine whether the vehicle passes through that meets the threshold requirements is regarded as the vibration period [T 1 , T 2 ] generated by the excitation of the driving vehicle. Then intercept the vibration curve fr'(t) in this period from the superimposed vibration curve fr(t);
e)根据车辆在相邻测量断面激励起振的时间差,求得该车辆移动速率为:e) According to the time difference between the excitation and vibration of the vehicle at the adjacent measurement section, the moving speed of the vehicle is obtained as:
其中,L为相邻断面的间距;Δt为时间差;Among them, L is the distance between adjacent sections; Δt is the time difference;
f)依据行驶车辆的最小轴距(即下面公示中的S),计算确定用于轴型特征识别的短时能量分析时间帧长度为:f) According to the minimum wheelbase of the driving vehicle (that is, S in the publicity below), calculate and determine the short-term energy analysis time frame length for axle type feature recognition as:
其中,v表示车辆移动速率;S为联轴轴距,通常取1m;Among them, v represents the moving speed of the vehicle; S is the coupling wheelbase, usually 1m;
g)对截取的振动曲线fr’(t)进行归一化处理,得到归一化后的振动曲线fn(t):g) Normalize the intercepted vibration curve fr'(t) to obtain the normalized vibration curve fn(t):
依据用于轴型特征识别的短时能量分析时间帧长度,计算归一化后的振动曲线fn(t)的短时能量时程曲线E’(T)。设定轴型识别阈值并确定满足阈值要求的曲线波峰数量,依据波峰数量判断行驶车辆轴型,再依据波峰间的距离确定轴距;According to the short-term energy analysis time frame length used for axial feature recognition, the short-term energy time-history curve E'(T) of the normalized vibration curve fn(t) is calculated. Set the axle type recognition threshold and determine the number of curve peaks that meet the threshold requirements, judge the axle type of the driving vehicle based on the number of peaks, and then determine the wheelbase based on the distance between the peaks;
h)参照步骤a)-g),计算车辆前后轴经过同一测量断面时各振动传感光纤传感段内的截取后的短时能量分布曲线,在单个测量断面内,将振动传感光纤传感段依据车辆行驶方向从左往右进行排序,并依此定义振动光纤传感段的序数,分别计算不同轴经过时的最大短时能量值,记为E(i,j),其中i为车轴的序数,j为振动光纤传感段的序数;h) With reference to steps a)-g), calculate the intercepted short-term energy distribution curves in each vibration sensing optical fiber sensing section when the front and rear axles of the vehicle pass through the same measurement section. The sensing segments are sorted from left to right according to the driving direction of the vehicle, and the ordinal numbers of the vibrating optical fiber sensing segments are defined accordingly, and the maximum short-term energy values when passing by different axes are calculated respectively, which are denoted as E(i,j), where i is the ordinal number of the axle, and j is the ordinal number of the vibrating optical fiber sensing section;
i)以前轴为基准,计算每个振动光纤传感段内车辆后轴与前轴最大短时能量值的比值,记为Er(i),其中i=1,2,3……n,n为单个测量断面上传感段的数量,计算Er(i)的标准差以衡量前后轴在同一断面内引起的振动差异,设定轮型判定阈值,标准差超过该阈值,则说明该轴型为双轮组,反之则为单轮组;i) Based on the front axle, calculate the ratio of the maximum short-term energy value of the vehicle rear axle to the front axle in each vibration optical fiber sensing segment, which is recorded as Er(i), where i=1,2,3...n,n For the number of sensing segments on a single measurement section, calculate the standard deviation of Er(i) to measure the vibration difference caused by the front and rear axles in the same section, and set the wheel type judgment threshold. If the standard deviation exceeds this threshold, it means that the axle type is Double wheel group, and vice versa for single wheel group;
j)采用时频分析手段获得截取的振动曲线fr’(t)的频谱分布,获得信号的幅频特性在时间上的分布特性S(f,t)(x轴为时间,y轴为频率,z轴为振动幅值),对幅频分布数据在时间方向(x轴)上进行叠加,叠加结果作为车辆振动响应的幅频曲线S(f)(x轴为频率,y轴为振动幅值);j) Obtain the frequency spectrum distribution of the intercepted vibration curve fr'(t) by means of time-frequency analysis, and obtain the distribution characteristic S(f,t) of the amplitude-frequency characteristic of the signal in time (the x-axis is time, the y-axis is frequency, The z-axis is the vibration amplitude), and the amplitude-frequency distribution data is superimposed on the time direction (x-axis), and the superposition result is used as the amplitude-frequency curve S(f) of the vehicle vibration response (the x-axis is the frequency, and the y-axis is the vibration amplitude );
k)截取对车辆振动响应敏感的频段为特征频段,计算该频段内的加权频率fw以表征车辆振动响应的频谱分布特征,计算方法为:k) Intercept the frequency band sensitive to the vehicle vibration response as the characteristic frequency band, and calculate the weighted frequency f w in this frequency band to characterize the spectrum distribution characteristics of the vehicle vibration response. The calculation method is:
其中,fh和fl分别为特征频段的上下截止频率。Among them, f h and f l are the upper and lower cut-off frequencies of the characteristic frequency band respectively.
在本实用新型一些实施方式中,所述道路车辆类型识别方法中,还对行驶车辆进行车型识别先验参数训练;In some embodiments of the present utility model, in the road vehicle type identification method, prior parameter training for vehicle type identification is also performed on the driving vehicle;
在本实用新型一些实施方式中,所述道路车辆类型识别方法中,还对行驶车辆进行车型识别概率判别。In some embodiments of the present utility model, in the road vehicle type identification method, the vehicle type identification probability is also determined for the driving vehicle.
在本实用新型一些实施方式中,对行驶车辆进行车型识别先验参数训练过程中,将获得到的行驶车辆移动速率、轴距、振动频谱分布特征过支持向量机进行参数训练。In some embodiments of the present invention, during the prior parameter training of the vehicle type identification for the driving vehicle, the acquired moving speed, wheelbase, and vibration spectrum distribution characteristics of the driving vehicle are used for parameter training through the support vector machine.
在本实用新型一些实施方式中,对行驶车辆进行车型识别概率判别的过程中,依据计算获得的行驶车辆轴的数量、轴距将行驶车辆分为“中小型客车”、“大型货车”、“大客车或中小型货车”三类,提取行驶车辆的移动速率、轴距、振动频谱特征,按照联合概率判别方法进行判别,将“大客车或中小型货车”进一步区分为“双轴大客车”和“双轴货车”两类。In some embodiments of the present invention, in the process of identifying the probability of vehicle type identification of the driving vehicle, the driving vehicle is divided into "small and medium-sized passenger cars", "large trucks", " According to the three categories of large buses or small and medium-sized trucks, the moving speed, wheelbase, and vibration spectrum characteristics of the driving vehicles are extracted, and the joint probability discrimination method is used to distinguish "large buses or small and medium-sized trucks" into "double-axle buses" and "twin-axle trucks".
附图说明Description of drawings
图1显示为本实用新型埋入式振动传感光纤主视图。Fig. 1 shows the front view of the embedded vibration sensing optical fiber of the utility model.
图2显示为本实用新型埋入式振动传感光纤俯视图。Fig. 2 shows a top view of the embedded vibration sensing optical fiber of the utility model.
图3显示为本实用新型埋入式道路车辆类型识别系统结构示意图。Fig. 3 is a schematic diagram showing the structure of the embedded road vehicle type identification system of the present invention.
图4显示为本实用新型埋入式车型识别方法流程图。Fig. 4 shows the flow chart of the utility model embedded vehicle identification method.
图5显示为本实用新型埋入式车型识别方法行驶车辆移动速率提取流程图。Fig. 5 shows a flow chart for extracting the moving speed of a driving vehicle in the embedded vehicle type identification method of the present invention.
图6显示为本实用新型埋入式车型识别方法行驶车辆轴型及轴距提取流程图。Fig. 6 shows a flowchart for extracting the axle type and wheelbase of the driving vehicle in the embedded vehicle type identification method of the present invention.
图7显示为本实用新型埋入式车型识别方法轮型提取流程图。Fig. 7 shows the flow chart of extracting the wheel type of the embedded vehicle type identification method of the present invention.
图8显示为本实用新型埋入式车型识别振动频谱特性提取流程图。Fig. 8 shows a flow chart of extracting vibration spectrum characteristics of the embedded vehicle type identification of the utility model.
图9显示为实施例中EMD分解示意图。Fig. 9 is a schematic diagram showing the decomposition of EMD in the embodiment.
图10显示为实施例1中一辆双轴货车及四轴货车的短时能量分布曲线图。Fig. 10 shows the short-term energy distribution curves of a two-axle truck and a four-axle truck in Example 1.
图11显示为实施例2中埋入式振动测量装置的排布示意图。FIG. 11 is a schematic diagram showing the layout of the embedded vibration measurement device in Embodiment 2.
元件标号说明Component designation description
1 振动传感光纤1 Vibration sensing fiber
11 振动传感光纤传感段11 Vibration sensing fiber optic sensing section
12 振动传感光纤过渡段12 Vibration sensing optical fiber transition section
2 振动光纤分析装置2 Vibration fiber optic analysis device
3 光纤引出线3 Optical fiber lead-out wires
具体实施方式Detailed ways
以下由特定的具体实施例说明本实用新型的实施方式,熟悉此技术的人士可由本说明书所揭露的内容轻易地了解本实用新型的其他优点及功效。The implementation of the present utility model is illustrated by specific specific examples below, and those skilled in the art can easily understand other advantages and effects of the present utility model from the content disclosed in this specification.
请参阅图1至图9。须知,本说明书所附图式所绘示的结构、比例、大小等,均仅用以配合说明书所揭示的内容,以供熟悉此技术的人士了解与阅读,并非用以限定本实用新型可实施的限定条件,故不具技术上的实质意义,任何结构的修饰、比例关系的改变或大小的调整,在不影响本实用新型所能产生的功效及所能达成的目的下,均应仍落在本实用新型所揭示的技术内容得能涵盖的范围内。同时,本说明书中所引用的如“上”、“下”、“左”、“右”、“中间”及“一”等的用语,亦仅为便于叙述的明了,而非用以限定本实用新型可实施的范围,其相对关系的改变或调整,在无实质变更技术内容下,当亦视为本实用新型可实施的范畴。See Figures 1 through 9. It should be noted that the structures, proportions, sizes, etc. shown in the drawings attached to this specification are only used to match the content disclosed in the specification, for those who are familiar with this technology to understand and read, and are not used to limit the implementation of the utility model Therefore, it has no technical substantive meaning. Any modification of structure, change of proportional relationship or adjustment of size shall still fall within the scope of The technical content disclosed by the utility model must be within the scope covered. At the same time, terms such as "upper", "lower", "left", "right", "middle" and "one" quoted in this specification are only for the convenience of description and are not used to limit this specification. The practicable range of the utility model, and the change or adjustment of its relative relationship, without any substantial change in the technical content, shall also be regarded as the practicable scope of the utility model.
如图1-图3所示,本实用新型提供一种埋入式道路车辆类型识别系统,包括路面结构本体,所述路面结构本体中设有振动传感光纤1,所述振动传感光纤包括一个或多个振动传感光纤传感段11,各振动传感光纤传感段11之间通过振动传感光纤过渡段12连接,还包括振动光纤分析装置2,所述振动光纤分析装置2与振动传感光纤1通过光纤引出线3相连接。As shown in Figures 1-3, the utility model provides an embedded road vehicle type identification system, which includes a road surface structure body, and a vibration sensing optical fiber 1 is arranged in the road surface structure body, and the vibration sensing optical fiber includes One or more vibration sensing optical fiber sensing sections 11, each vibration sensing optical fiber sensing section 11 is connected by a vibration sensing optical fiber transition section 12, and also includes a vibration optical fiber analysis device 2, and the vibration optical fiber analysis device 2 and Vibration sensing optical fibers 1 are connected through optical fiber lead-out lines 3 .
本实用新型所提供的埋入式道路车辆类型识别系统中,所述路面结构本体1可以是水泥混凝土铺面、沥青混凝土铺面或复合型铺面结构等中的一种或多种的组合,将振动传感光纤1铺设于路面结构本体1中,从而可以采集通过路面的车辆所造成的结构振动信息。In the embedded road vehicle type identification system provided by the utility model, the road surface structure body 1 can be one or more combinations of cement concrete pavement, asphalt concrete pavement or composite pavement structure, etc. The sensing optical fiber 1 is laid in the road structure body 1 so as to collect structural vibration information caused by vehicles passing the road.
本实用新型所提供的埋入式道路车辆类型识别系统中,本领域技术人员可选择合适的光纤种类用于采集车辆经过时的路面结构振动数据,例如,所述振动传感光纤1可以为单模光纤(Single Mode Fiber),所述单模光纤通常是指只能传送一种模式的光的光纤,所述振动传感光纤1可以为金属铠装光纤,所述振动传感光纤1的直径(含外套)可以为2~5mm,所述振动传感光纤1的长期允许拉伸力通常需要≥600N,短期允许拉伸力通常需要≥1500N(光纤短期和长期拉伸性能测试按照YD/T 769-2003《核心网用光缆-中心管式通信用室外光缆》进行),从而具有一定的抗拉强度,以保证光纤被拉伸时不会受损或断裂,所述振动传感光纤1的衰减≤0.2db/Km,以保证DOVS的顺利运行。In the embedded road vehicle type identification system provided by the utility model, those skilled in the art can select a suitable type of optical fiber for collecting the vibration data of the road surface structure when the vehicle passes by, for example, the vibration sensing optical fiber 1 can be a single Mode fiber (Single Mode Fiber), the single mode fiber usually refers to the optical fiber that can only transmit the light of one mode, the vibration sensing fiber 1 can be a metal armored fiber, the diameter of the vibration sensing fiber 1 (including the jacket) can be 2 to 5mm, the long-term allowable tensile force of the vibration sensing optical fiber 1 usually needs to be ≥ 600N, and the short-term allowable tensile force usually needs to be ≥ 1500N (the short-term and long-term tensile performance tests of optical fibers are based on YD/T 769-2003 "Optical Cables for Core Networks - Outdoor Optical Cables for Central Tube Communications"), so that it has a certain tensile strength to ensure that the optical fiber will not be damaged or broken when it is stretched. The vibration sensing optical fiber 1 Attenuation ≤0.2db/Km to ensure the smooth operation of DOVS.
本实用新型所提供的埋入式道路车辆类型识别系统中,所述振动传感光纤1通常包括多个振动传感光纤传感段11,各振动传感光纤传感段11之间的间距通常≥0.2m且≤0.5m,所述振动传感光纤传感11段通常为螺旋形(helix),螺旋形的长轴可以竖直分布(通常为相对于路面而言),振动传感光纤传感段11长轴方向的高度可以为5~30mm,振动传感光纤传感段11的直径可以为250-350mm,每个振动传感光纤传感段11中光纤的长度通常≥4m且≤6m。In the embedded road vehicle type identification system provided by the utility model, the vibration sensing optical fiber 1 usually includes a plurality of vibration sensing optical fiber sensing sections 11, and the distance between each vibration sensing optical fiber sensing section 11 is usually ≥0.2m and ≤0.5m, the 11 sections of the vibration sensing fiber optic sensor are usually helix, the long axis of the helix can be distributed vertically (usually relative to the road surface), the vibration sensing fiber optic sensor The height of the long-axis direction of the sensing section 11 can be 5-30mm, the diameter of the vibration-sensing optical fiber sensing section 11 can be 250-350mm, and the length of the optical fiber in each vibration-sensing optical fiber sensing section 11 is usually ≥ 4m and ≤ 6m .
本实用新型所提供的埋入式道路车辆类型识别系统中,所述路面结构本体中设有多排振动传感光纤传感段11,所述振动传感光纤传感段11可以在车道宽度方向均匀分布,覆盖率为2~3个/米,测量时,所述路面结构本体中按车道的长度方向可以布设有多排振动传感光纤传感段11,更具体的,可以布设有2~3排振动传感光纤传感段11,各排振动传感光纤传感段11之间的间距可以为3~10m,至少部分或全部的振动传感光纤传感段11依次串联。本领域技术人员可根据需要调整振动传感光纤过渡段12的参数,振动传感光纤过渡段12的参数整体上可以与振动传感光纤传感段11基本相同,各车型识别装置之间通常需保持一定长度的振动传感光纤过渡段12,从而保证振动感知具有足够的空间分辨率,例如,在本实用新型的一实施方式中,各振动传感光纤传感段11之间的振动传感光纤过渡段12可以≥0.2米。In the embedded road vehicle type identification system provided by the utility model, multiple rows of vibration-sensing optical fiber sensing sections 11 are arranged in the road surface structure body, and the vibration-sensing optical fiber sensing sections 11 can Evenly distributed, the coverage rate is 2 to 3 per meter. When measuring, multiple rows of vibration sensing optical fiber sensing sections 11 can be arranged in the road surface structure body according to the length direction of the lane. More specifically, 2 to 3 There are three rows of vibration-sensing optical fiber sensing sections 11, the distance between each row of vibration-sensing optical fiber sensing sections 11 can be 3-10m, and at least some or all of the vibration-sensing optical fiber sensing sections 11 are connected in series in sequence. Those skilled in the art can adjust the parameters of the vibration-sensing optical fiber transition section 12 as required, and the parameters of the vibration-sensing optical fiber transition section 12 can be basically the same as the vibration-sensing optical fiber sensing section 11 on the whole. Keep a certain length of the vibration sensing optical fiber transition section 12, so as to ensure that the vibration sensing has sufficient spatial resolution, for example, in an embodiment of the present utility model, the vibration sensing between each vibration sensing optical fiber sensing section 11 The optical fiber transition section 12 may be ≥0.2 meters.
本实用新型所提供的埋入式道路车辆类型识别系统中,所述振动光纤分析装置2可选用本领域内各种分布式光纤振动解调仪,其可以检测所连入光纤沿线的振动信号,例如可以是上海拜安传感技术有限公司生产的FT630-02光纤振动传感分析仪等。In the embedded road vehicle type identification system provided by the utility model, the vibration optical fiber analysis device 2 can be selected from various distributed optical fiber vibration demodulators in the field, which can detect the vibration signal along the connected optical fiber, For example, it can be the FT630-02 fiber optic vibration sensor analyzer produced by Shanghai Baian Sensing Technology Co., Ltd.
本实用新型还提供一种道路车辆类型识别方法,使用如上所述的埋入式道路车辆类型识别系统,测量时,可以将所述道路车辆类型识别系统布设于道路结构中,当车辆经过时,通过振动光纤解调设备采集光纤的振动数据,保存作为原始的振动数据,如图4所示,具体可以包括如下步骤:The utility model also provides a road vehicle type identification method, using the above-mentioned embedded road vehicle type identification system, when measuring, the road vehicle type identification system can be arranged in the road structure, when the vehicle passes by, Collect the vibration data of the optical fiber through the vibration fiber demodulation device, and save it as the original vibration data, as shown in Figure 4, which may specifically include the following steps:
1)采用如上所述的埋入式车型识别系统,通过振动光纤分析装置采集车辆经过时的路面结构振动数据,获得原始振动数据,具体为通过振动光缆解调车辆经过光缆时,设备采集光缆的振动数据,并保存作为原始的振动数据;1) Using the embedded vehicle identification system as described above, the vibration data of the road surface structure when the vehicle passes through the vibration optical fiber analysis device is used to obtain the original vibration data. Vibration data, and save as the original vibration data;
2)对各振动数据组进行特征提取,所述特征可以包括行驶车辆的移动速率、轴型、轮型、轴距及振动频谱分布中的一种或多种;2) feature extraction is carried out to each vibration data group, and described feature can comprise one or more in the speed of movement of driving vehicle, axle type, wheel type, wheelbase and vibration spectrum distribution;
3)依据上述车辆的特征信息,用于对行驶车辆进行分类的特征信息为轴型和/或振动频谱分布特性,在本实用新型一实施方式中,可以对行驶车辆进行分类,将行驶车辆分为“中小型客车”“大型货车”“大客车或中小型货车”三类,判断方法可以是例如:首先对车辆轴型进行判断,若振动时段内曲线波峰为2个,即该车辆有2根车轴,类型为“中小型客车”或“大客车或中小型货车”;若波峰有3个及以上,则类型为“大型货车”,其次,依据车辆前后轴轴距确定轴距判定阈值,若轴距小于该阈值,则其为“中小型客车”,反之则为“大客车或中小型货车”;3) According to the characteristic information of the above-mentioned vehicles, the characteristic information used to classify the driving vehicles is the axis type and/or the vibration spectrum distribution characteristics. In one embodiment of the present invention, the driving vehicles can be classified, and the driving vehicles There are three types of "small and medium passenger cars", "large trucks" and "big buses or small and medium trucks". The judgment method can be, for example: first judge the axle type of the vehicle. The root axle, the type is "small and medium-sized passenger car" or "big bus or small and medium-sized truck"; if there are 3 or more peaks, the type is "large truck". If the wheelbase is smaller than the threshold, it is a "small and medium-sized passenger car", otherwise it is "large bus or small and medium-sized truck";
4)对行驶车辆进行车型识别先验参数训练,更具体可以是对“大客车或中小型货车”的行驶车辆进行车型识别先验参数训练,对行驶车辆进行车型识别先验参数训练过程中,可以将获得到的行驶车辆移动速率、轴数、轴距、振动频谱分布特征过支持向量机进行参数训练;4) Carry out vehicle type recognition prior parameter training on driving vehicles, more specifically, it may be to carry out vehicle type recognition prior parameter training on "buses or small and medium-sized trucks" driving vehicles, during the process of performing vehicle type recognition prior parameter training on driving vehicles, The obtained moving speed, number of axles, wheelbase, and vibration spectrum distribution characteristics of the driving vehicle can be used for parameter training through the support vector machine;
5)对行驶车辆进行车型识别概率判别,更具体可以是进行“大客车或中小型货车”行驶车辆的车型识别概率判别,对行驶车辆进行车型识别概率判别的过程中,可以提取行驶车辆的移动速率、轴距、振动频谱特征,按照联合概率判别方法进行判别。5) Carry out vehicle type identification probability discrimination for driving vehicles, more specifically, it can be used to determine the vehicle type identification probability of "buses or small and medium-sized trucks". During the process of vehicle type identification probability discrimination, the movement The speed, wheelbase, and vibration spectrum features are judged according to the joint probability discrimination method.
本实用新型所提供的的道路车辆类型识别方法中,所述步骤2)中,对各振动数据组进行特征提取的方法具体包括如下步骤:In the road vehicle type identification method provided by the utility model, in the step 2), the method for feature extraction of each vibration data group specifically includes the following steps:
a)将同一横断面内测点所观测的原始振动数据编入同一振动数据组,并进行经验模态分解(EMD),获得多个经EMD处理后的振动分量;a) Compile the original vibration data observed at the measuring points in the same cross-section into the same vibration data group, and perform empirical mode decomposition (EMD) to obtain multiple vibration components after EMD processing;
b)叠加特定阶数(依据实测的振动数据调整阶数(频率区间)的区间范围和数量,例如可以是2~9阶或3~10阶)的振动分量,获得叠加后的振动曲线fr(t)(x轴为时间,y轴为振动强度):b) Superimpose the vibration components of a specific order (according to the measured vibration data to adjust the interval range and quantity of the order (frequency interval), for example, it can be 2-9 or 3-10) vibration components, and obtain the superimposed vibration curve fr( t) (x-axis is time, y-axis is vibration intensity):
其中IMFi为EMD处理后第i阶的振动分量,n和m为特定阶数的最低阶数和最高阶数。Among them, IMF i is the vibration component of the i-th order after EMD treatment, and n and m are the lowest and highest orders of specific orders.
c)计算单位时间(τ)内振动数据的平方和作为该段时间内的短时能量,依此得到振动短时能量的时程曲线E(T)(x轴为时间,y轴为振动短时能量):c) Calculate the sum of the squares of the vibration data per unit time (τ) as the short-term energy within this period, and thus obtain the time-history curve E(T) of the short-term vibration energy (the x-axis is time, and the y-axis is the short-term vibration energy time energy):
其中单位时间τ即为短时能量分析时间帧长度,初始值设置可以为0.5s。The unit time τ is the length of the short-term energy analysis time frame, and the initial value can be set to 0.5s.
d)设定车辆振动能量判定阈值,用于判定车辆是否经过的短时能量时程曲线中满足该阈值要求的部分视为因行驶车辆激励产生的振动时段[T1,T2],随后从叠加的振动曲线fr(t)截取出该时段内的振动曲线fr’(t);d) Set the vehicle vibration energy judgment threshold, which is used to judge whether the part of the short-term energy time-history curve that the vehicle passes through meets the threshold requirement as the vibration period [T 1 , T 2 ] generated by the excitation of the driving vehicle, and then from The superimposed vibration curve fr(t) intercepts the vibration curve fr'(t) within this period;
e)根据车辆在相邻测量断面激励起振的时间差,求得该车辆移动速率为:e) According to the time difference between the excitation and vibration of the vehicle at the adjacent measurement section, the moving speed of the vehicle is obtained as:
其中,L为相邻断面的间距;Δt为时间差。Among them, L is the distance between adjacent sections; Δt is the time difference.
f)依据行驶车辆的最小轴距(即下面公示中的S),计算确定用于轴型特征识别的短时能量分析时间帧长度为:f) According to the minimum wheelbase of the driving vehicle (that is, S in the publicity below), calculate and determine the short-term energy analysis time frame length for axle type feature recognition as:
其中,v表示车辆移动速率;S为联轴轴距,通常取1m;Among them, v represents the moving speed of the vehicle; S is the coupling wheelbase, usually 1m;
g)对截取的振动曲线fr’(t)进行归一化处理,得到归一化后的振动曲线fn(t):g) Normalize the intercepted vibration curve fr'(t) to obtain the normalized vibration curve fn(t):
依据用于轴型特征识别的短时能量分析时间帧长度,计算归一化后的振动曲线fn(t)的短时能量时程曲线E’(T)。设定轴型识别阈值并确定满足阈值要求的曲线波峰数量,依据波峰数量判断行驶车辆轴型,再依据波峰间的距离确定轴距;According to the short-term energy analysis time frame length used for axial feature recognition, the short-term energy time-history curve E'(T) of the normalized vibration curve fn(t) is calculated. Set the axle type recognition threshold and determine the number of curve peaks that meet the threshold requirements, judge the axle type of the driving vehicle based on the number of peaks, and then determine the wheelbase based on the distance between the peaks;
h)依据上述步骤(步骤a-g),计算车辆前后轴经过同一测量断面时各振动传感光纤传感段内的截取后的短时能量分布曲线,在单个测量断面内,将振动传感光纤传感段依据车辆行驶方向从左往右进行排序,并依此定义振动光纤传感段的序数,分别计算不同轴经过时的最大短时能量值,记为E(i,j),其中i为车轴的序数,j为振动光纤传感段的序数;h) According to the above steps (steps a-g), calculate the intercepted short-term energy distribution curves in each vibration sensing optical fiber sensing section when the front and rear axles of the vehicle pass through the same measurement section, and in a single measurement section, the vibration sensing optical fiber sensor The sensing segments are sorted from left to right according to the driving direction of the vehicle, and the ordinal numbers of the vibrating optical fiber sensing segments are defined accordingly, and the maximum short-term energy values when passing by different axes are calculated respectively, which are denoted as E(i,j), where i is the ordinal number of the axle, and j is the ordinal number of the vibrating optical fiber sensing section;
i)以前轴为基准,计算每个振动光纤传感段内车辆后轴与前轴最大短时能量值的比值(参照上述E(i,j)),记为Er(i),其中i=1,2,3……n,n为单个测量断面上传感段的数量,计算Er(i)的标准差以衡量前后轴在同一断面内引起的振动差异,设定轮型判定阈值,标准差超过该阈值,则说明该轴型为双轮组,反之则为单轮组。i) Based on the front axle, calculate the ratio of the maximum short-term energy value of the vehicle rear axle to the front axle in each vibration optical fiber sensing section (refer to the above E(i,j)), denoted as Er(i), where i= 1,2,3...n, n is the number of sensing segments on a single measurement section, calculate the standard deviation of Er(i) to measure the vibration difference caused by the front and rear axles in the same section, set the wheel type judgment threshold, standard deviation If the threshold is exceeded, it means that the shaft type is a double-wheel set, otherwise it is a single-wheel set.
j)采用时频分析手段获得截取的振动曲线fr’(t)的频谱分布,时频分析手段可以为小波分解、短时傅里叶变换或S变换,最终获得信号的幅频特性在时间上的分布特性S(f,t)(x轴为时间,y轴为频率,z轴为振动幅值),对幅频分布数据在时间方向(x轴)上进行叠加,叠加结果作为车辆振动响应的幅频曲线S(f)(x轴为频率,y轴为振动幅值)。j) Obtain the spectral distribution of the intercepted vibration curve fr'(t) by means of time-frequency analysis, which can be wavelet decomposition, short-time Fourier transform or S-transform, and finally obtain the amplitude-frequency characteristics of the signal in time The distribution characteristics of S(f,t) (the x-axis is time, the y-axis is frequency, and the z-axis is vibration amplitude), the amplitude-frequency distribution data is superimposed on the time direction (x-axis), and the superposition result is used as the vehicle vibration response The amplitude-frequency curve S(f) (x-axis is frequency, y-axis is vibration amplitude).
k)截取对车辆振动响应敏感的频段为特征频段,计算该频段内的加权频率fw以表征车辆振动响应的频谱分布特征,计算方法为:k) Intercept the frequency band sensitive to the vehicle vibration response as the characteristic frequency band, and calculate the weighted frequency f w in this frequency band to characterize the spectrum distribution characteristics of the vehicle vibration response. The calculation method is:
其中,fh和fl分别为特征频段的上下截止频率。Among them, f h and f l are the upper and lower cut-off frequencies of the characteristic frequency band respectively.
针对现有车型识别手段的缺点,本实用新型通过DOVS(分布式光纤振动感知系统)针对车型识别技术的需要,提供了一种新型的、准确的、可大范围识别的基于振动感知的道路车辆类型识别方法、装置及其系统。Aiming at the shortcomings of the existing vehicle identification methods, the utility model provides a new, accurate, and widely identifiable road vehicle based on vibration perception through DOVS (Distributed Optical Vibration Sensing System) to meet the needs of vehicle identification technology. Type identification method, device and system thereof.
实施例1Example 1
浦东外环线洲海路段辅道为水泥混凝土路面,该路段交通组成复杂,包括小客车、大客车以及多种类型的货车。The auxiliary road of the Zhouhai section of the Pudong Outer Ring Road is a cement concrete pavement, and the traffic composition of this section is complex, including passenger cars, buses and various types of trucks.
使用所述的基于DOVS的埋入式车型识别系统对行驶车辆类型进行探测和识别,所述系统采用装配式施工方式进行安装。首先依据图1-图2所示缠绕方式构建埋入式车型识别装置,实施例中所使用的识别装置信息如下:采用一根400m长的振动传感光纤构建识别装置,振动光纤传感段采用单模光纤,护套采用金属铠装,总直径(含护套)为3mm。光纤传感段采用圆环形缠绕,缠绕直径为300mm,每次缠绕4圈,约3.8m长。The DOVS-based embedded vehicle type identification system is used to detect and identify the type of driving vehicle, and the system is installed in a prefabricated construction method. Firstly, the embedded vehicle identification device is constructed according to the winding method shown in Fig. 1-Fig. 2. The information of the identification device used in the embodiment is as follows: a 400m long vibration sensing optical fiber is used to construct the identification device, and the vibration optical fiber sensing section adopts Single-mode optical fiber, the sheath is metal armored, the total diameter (including the sheath) is 3mm. The optical fiber sensing section is wound in a circular shape, with a winding diameter of 300mm, 4 turns each time, and a length of about 3.8m.
然后通过图3所示布设方式搭建埋入式车型识别系统,实施例中振动光纤过渡段为直径3mm的单模光纤,每个过渡段长度相同,均为0.2m,所述多个振动测量装置通过过渡段光纤串联后,经扎带固定于装配式混凝土铺面板中的钢筋网上。整块板共布设24个振动测量装置,其中沿板长度方向共布设3排,每排8个,每一排间距约2.2m。光纤线圈绑扎于钢筋网之后,随混凝土浇筑于装配式铺面板内,光纤位置距离装配板顶面约7cm。该系统在埋设于装配板之后,采用吊装的方式安装于既有路面中。系统随装配板搭建完成后,通过光纤引出线导出并连接于振动分析设备,光纤引出线类型、参数与传感段及过渡段相同,长度为25m。振动分析设备采用上海拜安传感技术有限公司生产的FT630-02光纤振动传感分析仪。Then build an embedded vehicle identification system through the layout shown in Figure 3. In the embodiment, the transition section of the vibration optical fiber is a single-mode optical fiber with a diameter of 3mm, and the length of each transition section is the same, which is 0.2m. The multiple vibration measurement devices After the optical fiber is connected in series through the transition section, it is fixed on the steel mesh in the prefabricated concrete paving panel by cable ties. A total of 24 vibration measuring devices are arranged on the whole board, among which there are 3 rows arranged along the length direction of the board, each row has 8 devices, and the distance between each row is about 2.2m. After the optical fiber coil is bound to the steel mesh, it is poured with concrete into the prefabricated paving panel, and the position of the optical fiber is about 7cm away from the top surface of the prefabricated panel. After the system is buried in the assembly plate, it is installed in the existing road surface by means of hoisting. After the system is built with the assembly board, it is exported and connected to the vibration analysis equipment through the optical fiber lead-out line. The type and parameters of the optical fiber lead-out line are the same as those of the sensing section and the transition section, and the length is 25m. The vibration analysis equipment adopts the FT630-02 fiber optic vibration sensor analyzer produced by Shanghai Bayan Sensing Technology Co., Ltd.
依据埋设的车型识别系统,采集了超过5小时的路面振动信号,其中包括217辆行驶车辆产生的振动信号。获取原始振动信号之后,采用埋入式车型识别方法(如图4~8所示)对不同车辆类型进行识别,整个流程采用MATLAB软件进行实现。[l1]According to the buried vehicle identification system, the vibration signals of the road surface were collected for more than 5 hours, including the vibration signals generated by 217 moving vehicles. After obtaining the original vibration signal, the embedded vehicle type identification method (as shown in Figure 4-8) is used to identify different vehicle types, and the whole process is realized by MATLAB software. [l1]
捕获原始振动信号之后,依据图4所示步骤进行车型识别。在步骤2中,依据该路段捕获的振动信号,原始振动信号在EMD分解之后,叠加2~9阶分量获得相应时程曲线(图9),即:After the original vibration signal is captured, the vehicle model is identified according to the steps shown in Figure 4. In step 2, according to the captured vibration signal of this road section, after the original vibration signal is decomposed by EMD, the 2nd to 9th order components are superimposed to obtain the corresponding time history curve (Fig. 9), namely:
其中,IMFi为EMD处理后第i阶的振动分量。Among them, IMF i is the i-th order vibration component after EMD treatment.
获取叠加后的振动时程曲线后,计算单位时间(τ)内振动数据的平方和作为该段时间内的短时能量,依此得到振动短时能量的时程曲线E(T)(x轴为时间,y轴为振动短时能量):After obtaining the superimposed vibration time-history curve, calculate the sum of the squares of the vibration data per unit time (τ) as the short-term energy within this period, and thus obtain the time-history curve E(T) of the vibration short-term energy (x-axis is time, and the y-axis is vibration short-term energy):
其中单位时间τ即为短时能量分析时间帧长度,取值为0.5s。The unit time τ is the length of the short-term energy analysis time frame, and the value is 0.5s.
依据经验设定车辆振动判定阈值为无车状态下短时能量均值的二倍,并依此判定短时能量超出车辆振动判定阈值的时段为车辆行驶经过的时段,从而截取车辆经过时段的振动曲线fr’(t)。Based on experience, the vehicle vibration judgment threshold is set to be twice the average value of the short-term energy in the no-car state, and the time period when the short-term energy exceeds the vehicle vibration judgment threshold is determined as the time period when the vehicle is driving, so as to intercept the vibration curve of the vehicle passing time period fr'(t).
以车辆前轴经过相邻两个测量断面时,短时能量曲线峰值对应的时刻的差值作为车辆经过两个测量断面的时间差,依此计算车辆的行驶车速,前后两个测量断面分别选择第一排和第三排光纤的埋设断面。When the front axle of the vehicle passes through two adjacent measurement sections, the difference of the moment corresponding to the peak value of the short-term energy curve is used as the time difference of the vehicle passing through the two measurement sections, and the driving speed of the vehicle is calculated accordingly. The buried sections of the first and third rows of optical fibers.
依据行驶车速计算每辆车的短时能量分析时间帧长度,该实施例中测量车辆的车速主要分布于8~25m/s,对应的短时能量分析时间帧长度分布于0.01s~0.06s。The length of the short-term energy analysis time frame of each vehicle is calculated according to the driving speed. In this embodiment, the speed of the measured vehicle is mainly distributed in the range of 8-25m/s, and the corresponding short-term energy analysis time frame length is distributed in the range of 0.01s-0.06s.
采用计算得到的时间帧长度计算对应时程曲线的短时能量,得到振动数据在时域上的短时能量分布曲线。采用MATLAB内的findpeaks函数对波峰进行提取,设置轴型识别阈值为3(归一化处理后的阈值),并依此确定行驶车辆的轴数。依据轴数可区分出“大型货车”这一类车型。图10所示为一辆双轴货车及四轴货车的短时能量分布曲线。The short-term energy corresponding to the time history curve is calculated by using the calculated time frame length, and the short-term energy distribution curve of the vibration data in the time domain is obtained. Use the findpeaks function in MATLAB to extract the peaks, set the axle type recognition threshold to 3 (threshold after normalization processing), and determine the number of axles of the driving vehicle accordingly. According to the number of axles, a category of "large trucks" can be distinguished. Figure 10 shows the short-term energy distribution curves of a two-axle truck and a four-axle truck.
计算单辆车相邻波峰的时间差,乘以行驶车速可计算行驶车辆的轴距。设定轴距判定阈值为3.4m,即轴距小于3.4m的行驶车辆可判定为“中小型客车”,反之则为“大客车或中小型货车”。Calculate the time difference between adjacent peaks of a single vehicle and multiply it by the speed of the vehicle to calculate the wheelbase of the vehicle. Set the wheelbase judgment threshold to 3.4m, that is, a vehicle with a wheelbase less than 3.4m can be judged as a "small and medium-sized passenger car", and vice versa as a "large bus or small and medium-sized truck".
大客车或中小型货车的分类采用车型识别先验参数训练和车型识别概率判别方法进行。车型识别先验参数训练采用支持向量机方法进行,提取的训练参数为频谱分布特征,行驶车速和前后轴轴距。首先对截取的振动曲线采用短时傅里叶变换:The classification of buses or small and medium-sized trucks is carried out by the prior parameter training of vehicle type identification and the probability discrimination method of vehicle type identification. The prior parameter training of car model recognition is carried out by support vector machine method, and the extracted training parameters are spectral distribution features, driving speed and front and rear axle distance. First, short-time Fourier transform is used for the intercepted vibration curve:
其中γ(t)为窗函数,本实施例中选用的为矩形窗,窗宽度为0.256s。S(t,f)为幅值与时间和频率的函数(x轴为时间,y轴为频率,z轴为振动幅值)。再对S(t,f)在时间方向(x轴)上进行叠加,叠加结果作为车辆振动响应的幅频曲线S(f)(x轴为频率,y轴为振动幅值)。Wherein γ(t) is a window function, a rectangular window is selected in this embodiment, and the window width is 0.256s. S(t,f) is a function of amplitude versus time and frequency (time on the x-axis, frequency on the y-axis, and vibration amplitude on the z-axis). Then S(t, f) is superimposed on the time direction (x-axis), and the superposition result is used as the amplitude-frequency curve S(f) of the vehicle vibration response (x-axis is frequency, y-axis is vibration amplitude).
以0~40Hz作为本实施例中行驶车辆的特征频段,计算该频段内的加权频率fw以表征车辆振动响应的频谱分布特征,计算方法为:Taking 0-40Hz as the characteristic frequency band of the driving vehicle in this embodiment, the weighted frequency fw in this frequency band is calculated to characterize the spectrum distribution characteristics of the vehicle vibration response. The calculation method is:
获取得到每一辆行驶车辆的频谱分布特征后,选择60个行驶车辆的加权频率、车速、轴距进行支持向量机训练,构造支持向量机网络。支持向量机训练过程中的核函数采用径向基函数,其公式如下:After obtaining the spectrum distribution characteristics of each driving vehicle, select the weighted frequency, vehicle speed, and wheelbase of 60 driving vehicles for support vector machine training, and construct a support vector machine network. The kernel function in the training process of the support vector machine adopts the radial basis function, and its formula is as follows:
其中,x’为支持向量的值,x为需要分类的样本值,σ为核函数宽度,在该向量机网络中取值为0.5。在该实施例中,支持向量机训练和验证均采用MATLAB软件内LIBSVM工具箱实现。进行随机抽选60份行驶车辆数据作为训练集,其余数据作为验证集,采用训练集训练生成的向量机网络进行分类验证。Among them, x' is the value of the support vector, x is the sample value to be classified, σ is the width of the kernel function, and the value is 0.5 in the vector machine network. In this embodiment, the training and verification of the support vector machine are realized by using the LIBSVM toolbox in the MATLAB software. Randomly select 60 pieces of driving vehicle data as the training set, and the rest of the data as the verification set, and use the vector machine network generated by training the training set for classification verification.
依据上述步骤对车辆进行分类,分类结果如表1。可以看出,本方面所提供的埋入式车型识别装置、系统及方法能有效地识别不同车型,避免了传统识别方式成本高、可靠性不佳的弊端,实现了水泥混凝土路面上车型的实时识别。According to the above steps, the vehicles are classified, and the classification results are shown in Table 1. It can be seen that the embedded vehicle type identification device, system and method provided by this aspect can effectively identify different vehicle types, avoid the disadvantages of high cost and poor reliability of traditional identification methods, and realize real-time identification of vehicle types on cement concrete roads. identify.
表1车辆类型识别结果Table 1 Vehicle type recognition results
实施例2Example 2
上海市浦东新区华京路-台北东路T形交叉口,铺面结构为混凝土结构,交通组成包含小客车、公交车、小型货车等。At the T-shaped intersection of Huajing Road-Taipei East Road, Pudong New District, Shanghai, the pavement structure is a concrete structure, and the traffic composition includes passenger cars, buses, and minivans.
采用所述的基于DOVS的埋入式车型识别系统对交叉口内的车辆类型进行探测和识别,所述系统采用装配式施工方式进行安装。首先依据图1所示缠绕方式构建埋入式车型识别装置,使用的识别装置信息如下:采用一根1200m长的振动传感光纤构建识别装置,振动光纤传感段采用单模光纤,护套采用金属铠装,总直径(含护套)为3mm。光纤传感段采用圆环形缠绕,缠绕直径为300mm,每次缠绕4圈,约4m长。The DOVS-based embedded vehicle type identification system is used to detect and identify the types of vehicles in the intersection, and the system is installed in a prefabricated construction method. First, build the embedded vehicle identification device according to the winding method shown in Figure 1. The information of the identification device used is as follows: a 1200m long vibration sensing optical fiber is used to construct the identification device, the vibration optical fiber sensing section adopts single-mode optical fiber, and the sheath adopts Metal armored, with an overall diameter (including sheath) of 3mm. The optical fiber sensing section is wound in a circular shape, with a winding diameter of 300mm, 4 turns each time, and a length of about 4m.
然后通过图3所示布设方式搭建埋入式车型识别系统,每块铺面板中振动光纤过渡段为直径3mm的单模光纤,每个过渡段长度相同,均为0.2m,所述多个识别装置通过过渡段光纤进行串联后,经扎带固定于钢筋网上。有7块板中布设21个识别装置,沿板长度方向布设3排,每排7个,即3×7布设;有2块板布设数量为7个(1×7);有1块板布设数量为14个(2×7);有1块板布设数量为36个(6×6)。每块板中均通过光纤引出线导出振动光纤并置于板边预留的光纤盒内。光纤引出线类型、参数与传感段及过渡段相同,长度为25m。现场所有预制板装配完毕后,将相邻铺面板间光纤盒内振动光纤进行熔接,并通过光纤引出线与振动分析设备相连。形成整体13块板的埋入式车型识别系统(如图11所示)。振动分析设备采用上海拜安传感技术有限公司生产的FT630-02光纤振动传感分析仪。Then build an embedded vehicle identification system through the layout shown in Figure 3. The transition section of the vibrating optical fiber in each paving panel is a single-mode optical fiber with a diameter of 3mm. The length of each transition section is the same, both being 0.2m. After the device is connected in series through the optical fiber in the transition section, it is fixed on the steel mesh by cable ties. There are 21 identification devices arranged in 7 boards, and 3 rows are arranged along the length direction of the board, with 7 devices in each row, that is, 3×7 layout; 2 boards are arranged with 7 pieces (1×7); there is 1 board arranged The quantity is 14 (2×7); there is 1 board and the quantity is 36 (6×6). The vibrating optical fiber is led out from each board through the optical fiber lead-out line and placed in the fiber box reserved on the edge of the board. The type and parameters of the optical fiber lead-out line are the same as those of the sensing section and the transition section, and the length is 25m. After all the prefabricated panels are assembled on site, the vibrating optical fibers in the optical fiber boxes between adjacent paving panels are fused and connected to the vibration analysis equipment through the optical fiber lead-out lines. Form the embedded vehicle type identification system (as shown in Figure 11) of 13 boards as a whole. The vibration analysis equipment adopts the FT630-02 fiber optic vibration sensor analyzer produced by Shanghai Bayan Sensing Technology Co., Ltd.
最后采用埋入式车型识别方法(如图4~8所示)对40辆车辆的类型进行测量和识别,由于数据量较少,仅依据轴型参数,对中小型客车、大客车或中小型货车、大型货车三类车型进行识别,识别过程中数据处理方法和流程与实施例1相同。三种类型车辆的总体识别正确率为95%。Finally, the embedded vehicle type identification method (as shown in Figure 4-8) is used to measure and identify the types of 40 vehicles. Due to the small amount of data, only based on the axle type parameters, small and medium-sized buses, buses or small and medium-sized Trucks and large trucks are identified, and the data processing method and flow in the identification process are the same as those in Embodiment 1. The overall recognition accuracy of the three types of vehicles is 95%.
综上所述,本实用新型有效克服了现有技术中的种种缺点而具高度产业利用价值。To sum up, the utility model effectively overcomes various shortcomings in the prior art and has high industrial application value.
上述实施例仅例示性说明本实用新型的原理及其功效,而非用于限制本实用新型。任何熟悉此技术的人士皆可在不违背本实用新型的精神及范畴下,对上述实施例进行修饰或改变。因此,举凡所属技术领域中具有通常知识者在未脱离本实用新型所揭示的精神与技术思想下所完成的一切等效修饰或改变,仍应由本实用新型的权利要求所涵盖。The above-mentioned embodiments only illustrate the principles and effects of the present utility model, but are not intended to limit the present utility model. Anyone familiar with this technology can modify or change the above-mentioned embodiments without departing from the spirit and scope of the present invention. Therefore, all equivalent modifications or changes made by those with ordinary knowledge in the technical field without departing from the spirit and technical ideas disclosed in the utility model should still be covered by the claims of the utility model.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201721385512.XU CN207611485U (en) | 2017-10-25 | 2017-10-25 | An Embedded Road Vehicle Type Recognition System |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201721385512.XU CN207611485U (en) | 2017-10-25 | 2017-10-25 | An Embedded Road Vehicle Type Recognition System |
Publications (1)
Publication Number | Publication Date |
---|---|
CN207611485U true CN207611485U (en) | 2018-07-13 |
Family
ID=62793974
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201721385512.XU Active CN207611485U (en) | 2017-10-25 | 2017-10-25 | An Embedded Road Vehicle Type Recognition System |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN207611485U (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107730895A (en) * | 2017-10-25 | 2018-02-23 | 同济大学 | A kind of flush type road vehicle type identifying system and method |
-
2017
- 2017-10-25 CN CN201721385512.XU patent/CN207611485U/en active Active
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107730895A (en) * | 2017-10-25 | 2018-02-23 | 同济大学 | A kind of flush type road vehicle type identifying system and method |
CN107730895B (en) * | 2017-10-25 | 2024-01-30 | 同济大学 | Buried road vehicle type identification system and method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107730895B (en) | Buried road vehicle type identification system and method | |
CN102622783B (en) | Multilane free-flow electronic toll collection method based on accurate position match | |
CN103164962B (en) | Mountain road sharp turn section real-time vehicle speed early warning method | |
CN107256635A (en) | Vehicle identification method based on distributing optical fiber sensing in a kind of intelligent transportation | |
CN109523780A (en) | A kind of mountainous area highway snowfall environment early warning system and method | |
CN110060484A (en) | A kind of real-time early warning system and method violating the regulations of the highway passenger transportation based on block chain | |
CN106600981B (en) | Road section dual-way vehicle speed estimation method based on distributed sensing information | |
CN105448106A (en) | Vehicle detection device based on geomagnetic sensor array | |
CN104881897B (en) | ETC tracks vehicle detecting system and detection method | |
CN105206051B (en) | A kind of city road network Rationality Assessment method and device | |
CN105654720B (en) | Loop detector layout method based on urban road congestion identification | |
CN103198707B (en) | A vehicle control method based on the characteristics of dangerous traffic flow in sunny conditions | |
CN107240291A (en) | A kind of parking lot intelligent berth inducible system | |
CN114624726B (en) | Axle identification system and axle identification method | |
CN108109395A (en) | A kind of vehicle-state detecting system based on fiber grating | |
CN106767855B (en) | Method and system for updating electronic map data | |
CN106023326A (en) | Expressway non-stop toll collection method based on face recognition | |
CN101131796A (en) | Road traffic parameter checking device and method thereof | |
CN207611485U (en) | An Embedded Road Vehicle Type Recognition System | |
Shi et al. | Evaluation of road service performance based on human perception of vibration while driving vehicle | |
CN215298466U (en) | Site-unrestricted expressway entrance overload control and dissuasion system | |
CN103793985A (en) | Automatic lifting type portal frame control method based on geomagnetism detection | |
CN106601006B (en) | An urban intelligent traffic guidance device based on RFID and WeChat platform | |
CN202013661U (en) | Traffic driving guidance system based on floating car GPS data | |
CN116740973A (en) | Long tunnel group full-chain cooperative dynamic vehicle speed control method and system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
GR01 | Patent grant | ||
GR01 | Patent grant |