Disclosure of Invention
The invention aims to provide a mobile interference vibration source filtering method based on an image operator, which is used for rapidly filtering a mobile interference vibration source through an established vibration source moving speed estimation model.
In order to achieve the purpose, the technical scheme of the invention is as follows:
a mobile interference vibration source filtering method based on an image operator comprises the steps of receiving an alarm vibration source signal, wherein the alarm vibration source signal is a vibration source signal with the amplitude larger than 0, generating a two-dimensional array by the received vibration source signal, the ordinate of the two-dimensional array is a time axis, the abscissa of the two-dimensional array is a spatial position axis, and obtaining the position of the center of the alarm vibration source signal on the time axis and the position of the center of the alarm vibration source signal on the spatial position axis, wherein the method further comprises the following steps:
the first step is as follows: establishing a two-dimensional angle operator based on a time axis and a spatial position axis, and providing an angle distribution template within a semicircular range with the radius of r from the center of the operator by the operator;
the second step is that: acquiring angle distribution data between the vibration source signal and the vibration source signal center by using an angle distribution template;
the third step: establishing a vibration source moving speed estimation model according to the angle distribution data;
the fourth step: extracting the moving speed of a vibration source signal through a vibration source moving speed estimation model, and determining the type of the vibration source according to the moving speed of the vibration source: the vibration source with the moving speed greater than the fixed speed threshold is identified as a harmless interference vibration source, and the vibration source with the moving speed less than the fixed speed threshold is identified as an alarm vibration source;
wherein:
the first step is as follows: the two-dimensional angle operator is established, and the operator provides an angle distribution template within a semicircular range with the radius of r from the center of the operator, wherein the angle distribution template comprises the following steps:
firstly, the time axis and the spatial position axis are de-unitized through a de-unitization formula,
(formula for removing units)
Wherein:
t represents system time in units of s;
tnowrepresents the current running time of the system in units of s;
turepresents the de-unitization factor, take tu=1s;
x represents the coordinate of the time axis after the unitization, and is dimensionless;
z represents the monitored source distance in m;
zmaxrepresents the maximum distance monitored, in m;
zurepresenting the spatial de-unitization factor, take zu=10m;
y represents the space axis coordinate after the unitization, and is dimensionless;
secondly, establishing a two-dimensional array to form a two-dimensional angle operator omega
Wherein the content of the first and second substances,
θi,j(i is 1,2, …, r, j is 1,2, …,2r) is the corresponding angle value in the ith row and j column in the two-dimensional angle operator, and the value range is [ 0-180%],
Wherein the content of the first and second substances,
r is the effective radius of the two-dimensional angle operator;
round means rounding;
nan denotes a non-number;
the second step is that: the angle distribution data between the vibration source signal and the vibration source signal center is obtained by the angle distribution template:
firstly, extracting a two-dimensional time and space vibration source signal matrix A near an alarm point
Wherein the content of the first and second substances,
a represents the characteristic value of the vibration source signal, and the value range is [0, + ∞ ];
secondly, calculating an angle distribution matrix B:
wherein:
represents the direct product of two matrices;
sign () is a sign function;
thirdly, counting the number of each angle in the angle distribution matrix B to obtain an angle number vector
Wherein the content of the first and second substances,
representing angles in the angle distribution matrix B as
Obtaining the angle number variation curve along with the angle;
the third step: the establishment of the vibration source moving speed estimation model according to the angle distribution data comprises the following steps: according to the angle number vector N, establishing a vibration source moving speed estimation model,
(model for estimating the moving speed of vibration source)
Wherein the content of the first and second substances,
representing an angle vector
In the middle, the angle value with the largest number of angles satisfies
v represents the moving speed of the vibration source, m/s.
The scheme is further as follows: the fixed speed threshold is 1 m/s.
Compared with the prior art, the invention has the following advantages:
1. the operation efficiency is high, real-time data processing of the system is avoided, and only vibration source type identification is needed after alarm generation, so that the time complexity of the system is reduced, and the operation efficiency is improved;
2. the identification effect is obvious, the vehicle interference of parallel optical fiber movement can be accurately identified, the system alarm type identification rate is improved, the harmless interference vibration source alarm quantity is reduced, and the system effective alarm quantity is improved.
The invention is described in detail below with reference to the figures and examples.
Detailed Description
A mobile interference vibration source filtering method based on image operators is a method for judging a mobile interference vibration source by simulating a moving speed, and comprises the steps of receiving an alarm vibration source signal, judging whether the alarm vibration source signal is a vibration source mobile interference source signal or not, if so, not giving an alarm, and judging the vibration source mobile interference source signal by the following steps: the method comprises the following steps of obtaining optical fiber vibration signal data, preprocessing the vibration signal data, determining vibration characteristic data at corresponding positions of optical fibers, and establishing a two-dimensional array based on time and space arrangement, namely: generating a two-dimensional array from the received vibration source signal, where the ordinate of the two-dimensional array is a time axis, and the abscissa of the two-dimensional array is a spatial position axis, and obtaining the position of the center of the alarm vibration source signal on the time axis and the position of the center of the alarm vibration source signal on the spatial position axis, where fig. 1 is a schematic diagram of vibration signal data, where the method further includes:
the first step is as follows: establishing a two-dimensional angle operator based on a time axis and a spatial position axis, and providing an angle distribution template within a semicircular range with the radius of r from the center of the operator by the operator;
the second step is that: acquiring angle distribution data between the vibration source signal and the vibration source signal center by using an angle distribution template;
the third step: establishing a vibration source moving speed estimation model according to the angle distribution data;
the fourth step: extracting the moving speed of a vibration source signal through a vibration source moving speed estimation model, and determining the type of the vibration source according to the moving speed of the vibration source: the vibration source with the moving speed greater than the fixed speed threshold is identified as a harmless interference vibration source, and the vibration source with the moving speed less than the fixed speed threshold is identified as an alarm vibration source;
wherein:
the first step is as follows: the two-dimensional angle operator is established, and the operator provides an angle distribution template within a semicircular range with the radius of r from the center of the operator, wherein the angle distribution template comprises the following steps:
firstly, the time axis and the spatial position axis are de-unitized through a de-unitization formula,
(formula for Deunitization, formula 1)
Wherein:
t represents system time in units of s;
tnowrepresents the current running time of the system in units of s;
turepresents the de-unitization factor, take tu=1s;
x represents the coordinate of the time axis after the unitization, and is dimensionless;
z represents the monitored source distance in m;
zmaxrepresents the maximum distance monitored, in m;
zurepresenting the spatial de-unitization factor, take zu=10m;
y represents the space axis coordinate after the unitization, and is dimensionless;
secondly, establishing a two-dimensional array to form a two-dimensional angle operator omega
Wherein the content of the first and second substances,
θi,j(i is 1,2, …, r, j is 1,2, …,2r) is the corresponding angle value in the ith row and j column in the two-dimensional angle operator, and the value range is [ 0-180%],
Wherein the content of the first and second substances,
r is the effective radius of the two-dimensional angle operator, and simultaneously represents the row number of the operator matrix in the formula 2 to explain the dimension of the matrix;
round means rounding;
nan denotes a non-number;
the second step is that: the angle distribution data between the vibration source signal and the vibration source signal center is obtained by the angle distribution template:
firstly, extracting a two-dimensional time and space vibration source signal matrix A near an alarm point
Wherein the content of the first and second substances,
a represents the characteristic value of the vibration source signal, and the value range is [0, + ∞ ];
secondly, calculating an angle distribution matrix B:
wherein:
represents the direct product of two matrices;
sign () is a sign function;
thirdly, counting the number of each angle in the angle distribution matrix B to obtain an angle number vector
Wherein the content of the first and second substances,
representing angles in the angle distribution matrix B as
Obtaining the angle number variation curve along with the angle;
the third step: the establishment of the vibration source moving speed estimation model according to the angle distribution data comprises the following steps: according to the angle number vector N, establishing a vibration source moving speed estimation model,
(model of vibration source moving speed, equation 7)
Wherein the content of the first and second substances,
representing an angle vector
In the middle, the angle value with the largest number of angles satisfies
v represents the moving speed of the vibration source, m/s.
Wherein: the fixed speed threshold value ranges from 1m/s to 2m/s, and is 1m/s in the embodiment. The above method process can be further explained by the flow chart of fig. 2, and the process is as follows:
s201: carrying out signal detection and data acquisition on the optical fiber vibration signal;
s202: preprocessing the acquired data to obtain a vibration signal characteristic signal shown in figure 1;
s203: judging whether an alarm condition is met according to the existing method, and if the alarm condition is not met, continuing to perform signal detection operation;
s204: adding 1 to the system running time;
s205: establishing a two-dimensional angle operator according to the formula (2) and the formula (3), wherein a two-dimensional array stored in the r x 2r dimension is in omega:
s206: if the alarm condition is met, the alarm point position information obtained according to the existing method is shown as a five-pointed star in figure 1, and the alarm point position information is calculated according to the received alarm point coordinate (x) through an angle statistical model between the vibration signal and the alarm point0,y0) Counting vibration signal distribution angles near the alarm points;
firstly, a two-dimensional time and space vibration signal matrix A near an alarm point is extracted by a formula (4):
then, an angle distribution matrix B is calculated from equation (5):
finally, counting the number of each angle in the angle distribution matrix B to obtain an angle number vector
Wherein the content of the first and second substances,
representing angles in the angle distribution matrix B as
The angle number-dependent angle change curve shown in fig. 6 is obtained.
S207: detecting the moving speed v of the vibration source according to the formula (7) through a vibration source moving speed estimation model:
s208: and (3) judging whether the vibration source is a mobile interference vibration source or not by combining a formula (8) according to the detected vibration source speed:
wherein f represents a vibration source identification mapping model function; v. of0Indicating a speed threshold for determining whether the mobile interfering source is present.
S209: if the vibration source speed meets the mobile vibration source, the alarm caused by the mobile interference vibration source is judged, and the alarm is automatically filtered.
S210: if the vibration source speed does not meet the characteristics of the mobile vibration source, the system judges that the vibration source is a suspicious dangerous vibration source and carries out system alarm.
The two-dimensional angle operator established in the embodiment is shown in fig. 3, and the unitized radius adopted in the embodiment is 200, which can be adjusted as required; the alarm points and the distribution conditions of the vibration signals nearby the alarm points are shown in fig. 4, wherein the position of a pentagram is the position of the alarm point given by the existing method, and the signals in the annular wire frame are partial signals participating in angle distribution statistics; fig. 5 shows the distribution angle situation between the vibration signal and the alarm point obtained by the angle statistical model, and the curve of the number of angles obtained by statistics along with the change of the angle is shown in fig. 6; fig. 7 is a graph showing a relationship between a vibration source moving speed curve and a vibration signal obtained after the velocity estimation model. Therefore, the moving interference vibration source identification algorithm based on the image operator can well detect the moving speed of the vibration source, and has strong vibration source moving speed detection capability.
The method avoids real-time data processing of the system, and carries out vibration source type identification after alarm generation by a method of simulating the moving speed, thereby reducing the time complexity of the system and improving the operation efficiency.