CN110110278A - The calculation method of difference accumulation algorithm interval line number in fiber-optic vibration detection system - Google Patents
The calculation method of difference accumulation algorithm interval line number in fiber-optic vibration detection system Download PDFInfo
- Publication number
- CN110110278A CN110110278A CN201910396840.7A CN201910396840A CN110110278A CN 110110278 A CN110110278 A CN 110110278A CN 201910396840 A CN201910396840 A CN 201910396840A CN 110110278 A CN110110278 A CN 110110278A
- Authority
- CN
- China
- Prior art keywords
- vibration
- accumulation algorithm
- difference
- difference accumulation
- interval
- 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
Classifications
-
- 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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention discloses the calculation methods of difference accumulation algorithm interval line number in fiber-optic vibration detection system, for determining that two rows subtract each other the interval line number N between data when extracting vibration signal using difference accumulation algorithm, comprising: determine the dominant frequency period of vibration signal；Determine that the detection cycle of degree and light pulse is subtracted each other at the interval in difference accumulation algorithm；Subtract each other degree according to interval to be split the sinusoidal signal of a cycle, obtains segmentation number；Vibration signal is divided equally according to segmentation number, the time required for obtaining vibration difference once；Division arithmetic is carried out with the detection cycle of vibration difference primary required time and light pulse, obtains required light pulse number；Required light pulse number is adjusted, the interval line number N that difference accumulation algorithm can be obtained includes.The present invention is optimized difference accumulation algorithm there are problems that data volume is big and the high cost of the long bring of transmission time, large volume in fiber-optic vibration detection system.
Description
Technical field
The present invention relates to difference accumulation algorithms in optical fiber vibration sensing technical field more particularly to fiber-optic vibration detection system
It is spaced the calculation method of line number.
Background technique
The sensing mode of traditional sensing technology, sensing sensitivity and in terms of there is the case where bottleneck period
Under, optical fiber sensing technology is as a kind of New Sensing Technology, using optical fiber as sensor information, have high sensitivity, propagate with it is roomy,
The advantages such as strong real-time, in the fields such as architectural engineering, petrochemical industry, biologic medical, application prospect is huge, while Fibre Optical Sensor skill
Research of the art in vibration detection field is even more to receive the favor of researcher.
Current fiber-optic vibration assay device structures are complicated, for indexs such as the spatial resolution and the signal-to-noise ratio that guarantee system,
It needs to carry out algorithm process to the information of acquisition, can just extract vibration signal.The extraction algorithm of vibration signal can be divided into time domain
Processing Algorithm and frequency domain algorithm.Common algorithm has difference accumulation algorithm, cumulative mean algorithm etc. in time domain, in frequency domain often
The algorithm seen has wavelet transformation, Short Time Fourier Transform etc..Although these algorithms can extract vibration signal, use
The problem that continuous sampling mode brings that data volume is huge and transmission time is long, causes that system cost is high, volume is heavy.If
Although data volume can be reduced using triggering sampling and shorten transmission time, it will appear loss of data, time irreversibility etc.
Problem.Therefore to the optimization of signal extracting algorithm by biggish concern and research.
Traditional difference accumulation algorithm principle is the Rayleigh beacon signal acquired in fiber-optic vibration detection system
When amplitude data amount reaches data volume set by user (show that a frame waveform need data volume), by the data of acquisition according to
The length of optical fiber carries out branch, and the data of every a line indicate that light pulse is transmitted once in sensor fibre.By each row of data
Difference is carried out with the data every N row and goes signed magnitude arithmetic(al), is finally carried out accumulating operation and is obtained a frame data.By between frame
Variation obtains the specific location in oscillation point.
The selection of N value artificially tests obtained value generally by multiple in traditional difference accumulation algorithm.While N value
Choose extremely important, can direct decision systems detect the specific location in oscillation point.Artificial value obtains the method needs of N value
It takes a significant amount of time and tests, and the not specific foundation of N value choosing method, not can prove that N value at this time is best
Value.The problem that the continuous sampling mode that simultaneity factor uses brings that data volume is huge and transmission time is long, cause system at
The problems such as this is high, volume is heavy.
Summary of the invention
In view of this, the object of the present invention is to provide the meters of difference accumulation algorithm interval line number in fiber-optic vibration detection system
Calculation method is reduced in the case where guaranteeing the performance indicators such as detection system dynamic range, signal-to-noise ratio, dynamic response, frequency range
The data volume of system acquisition, thus the volume for reducing the cost of system, reducing system.
The present invention solves above-mentioned technical problem by following technological means:
The calculation method of difference accumulation algorithm interval line number in fiber-optic vibration detection system, for being calculated using difference is cumulative
Method extract vibration signal when determine that two rows subtract each other the interval line number N between data, by fiber-optic vibration detect in vibration signal regard
For a sinusoidal signal, the sinusoidal signal of a cycle is 360 °, comprising:
Determine the dominant frequency period of vibration signal；
Determine that the detection cycle of degree and light pulse is subtracted each other at the interval in difference accumulation algorithm；
Subtract each other degree according to interval to be split the sinusoidal signal of a cycle, obtains segmentation number；
Vibration signal is divided equally according to segmentation number, the time required for obtaining vibration difference once；
Division arithmetic is carried out with the detection cycle of vibration difference primary required time and light pulse, obtains required light arteries and veins
Rush number；
Required light pulse number is adjusted, the interval line number N of difference accumulation algorithm can be obtained.
Further, it is one in 5 °, 6 °, 7 ° or 8 ° that degree is subtracted each other at the interval.
Further, the time loss referring to vibration signal in acquisition, processing and transmission process, to required light pulse
Number is adjusted, to obtain the interval line number N of difference accumulation algorithm.
Beneficial effects of the present invention: there are data volumes under continuous sampling mode for fiber-optic vibration detection system by the present invention
N value calculating in difference accumulation algorithm is optimized in the problem of big and high cost of the long bring of transmission time, large volume.Phase
Compared under traditional continuous sampling mode by repeatedly debug obtain N value method, the N value that the present invention is obtained by theory analysis
The mode that can be applied to triggering sampling, had both solved the above problem existing for fiber-optic vibration detection system, while realizing
It triggers under sample mode, carries out the function of vibration signal detection.System cost is reduced, equipment volume is reduced, is convenient for optical fiber
The carrying and installation of vibration detecting system.
Detailed description of the invention
Fig. 1 is the calculating side of difference accumulation algorithm interval line number in fiber-optic vibration detection system provided in an embodiment of the present invention
The flow chart of method；
Fig. 2 is that N value chooses schematic illustration during difference provided in an embodiment of the present invention is tired；
Fig. 3 is the variation coordinate diagram of sinusoidal signal provided in an embodiment of the present invention sin value when 1 °~10 ° incremental；
Fig. 4 is the schematic diagram of continuous sampling mode provided in an embodiment of the present invention；
Fig. 5 is the schematic diagram of triggering sample mode provided in an embodiment of the present invention.
Specific embodiment
Below with reference to the drawings and specific embodiments, the present invention is described in detail:
As shown in Figure 1, in fiber-optic vibration detection system difference accumulation algorithm interval line number calculation method, for using
Difference accumulation algorithm determines that two rows subtract each other the interval line number N between data when extracting vibration signal.By Fourier transformationIt is found that any one signal can be combined shape by numerous sine wave
At, and the main energetic of signal concentrates on its dominant frequency sine wave, therefore by fiber-optic vibration detect in vibration signal can regard
For a sinusoidal signal.The sinusoidal signal of a cycle is 360 °, in difference accumulation algorithm by each row of data with every N row data into
Row calculus of differences it is practical be sinusoidal signal subtract each other every X °, as shown in Fig. 2, by fiber-optic vibration detect in vibration signal be considered as one
A sinusoidal signal, period are 360 °.Specific step includes:
Determine the dominant frequency period of vibration signal；
Determine that the detection cycle of degree and light pulse is subtracted each other at the interval in difference accumulation algorithm；
Subtract each other degree according to interval to be split the sinusoidal signal of a cycle, obtains segmentation number；
Vibration signal is divided equally according to segmentation number, the time required for obtaining vibration difference once；
Division arithmetic is carried out with the detection cycle of vibration difference primary required time and light pulse, obtains required light arteries and veins
Rush number；
Required light pulse number is adjusted, the interval line number N of difference accumulation algorithm can be obtained.
Specifically, abscissa expression is opened from 0 ° as shown in figure 3, indicating the variation of sinusoidal signal sin value when 1 °~10 °
Begin according to X ° of incremental ten point, ordinate indicates the corresponding sin value of abscissa, and different curves indicates sin value according to difference
It is that degree is incremented by as a result, dotted line therein indicates reference line, can be used to judge the linearity of different curves.It can be with from Fig. 3
Significantly find out that the linearity of sinusoidal signal sin value transformation when being spaced 5 °, 6 °, 7 °, 8 ° is good and variation is big.Therefore, of the invention
It is one in 5 °, 6 °, 7 ° or 8 ° that degree is subtracted each other at the interval of middle selection.
Carrying out citing below, the present invention will be described:
When vibration source frequency is 15Hz, the selection of difference accumulation algorithm is subtracted each other every 5 °, when the period of light pulse is 50us, N value
30 should be taken.Be explained in detail below: sinusoidal signal is subtracted each other by 5 ° can be divided into 72 parts, and the period of vibration signal is about 67ms, press
Divide equally according to 72 parts, show that vibration signal subtracts the primary time and takes around 0.93ms, and light pulse is primary by sensor fibre
Time is 50us, show that needing the time of about 20 light pulses can be only achieved subtracts each other primary 0.93ms.Simultaneously in view of letter
Time loss number in acquisition, transmission and treatment process, therefore N value takes 30.
The present invention is directed to fiber-optic vibration detection system long band of and transmission time big there are data volume under continuous sampling mode
Come high cost, large volume the problem of, in difference accumulation algorithm N value calculating be optimized.It is continuously adopted compared to traditional
By repeatedly debugging the method for obtaining N value under sample loading mode, the N value that the present invention is obtained by theory analysis can be applied to triggering and adopt
The mode of sample had both solved the above problem existing for fiber-optic vibration detection system, while realizing in the case where triggering sample mode, into
The function of row vibration signal detection.System cost is reduced, equipment volume is reduced, convenient for the carrying of fiber-optic vibration detection system
And installation.
The comparison of two kinds of sample modes is as shown in Figure 4 and Figure 5:
Traditional continuous sampling mode as shown in figure 4, traditional continuous sampling principle is acquired since first sampling point,
In order to guarantee the integrality of data, need continual acquisition extremely to the data volume for having acquired M row fiber lengths.Wherein in optical fiber
The period of light pulse needs to be adjusted to light pulse and is transmitted in sensor fibre once, is just can guarantee in this way by the data of acquisition
Equal timesharing, the initial position of the corresponding sensor fibre of the beginning of every data line.But since uninterrupted sampling makes the letter of data
Breath amount is huge, not only increases the cost of system, and increases equipment volume.Traditional difference N value selection is by multiple
What debugging obtained, the time cost of the user increased, and lack theoretical convincingness.
Modified triggering sample mode is as shown in figure 5, by the period phase of the frequency of trigger signal and light pulse in optical fiber
Together, the initial position of the corresponding sensor fibre of the beginning of every data line collected in this way.The period need to only be more than or equal to simultaneously
Light pulse is transmitted the primary time in sensor fibre, and the data volume for acquiring transmission primaries is a line sensor fibre
The data of length, such data acquisition modes can be converted to triggering sampling by continuous sampling.Sampling is triggered simultaneously not needing to connect
Continuous uninterrupted sampling, greatly reduces the data volume of system, reduces equipment volume.Pass through the difference N value selection side after optimization
Method specifies the selection thinking of N value, and making system, there are reasonability, convenient for the carrying and installation of fiber-optic vibration detection system.
The above examples are only used to illustrate the technical scheme of the present invention and are not limiting, although referring to preferred embodiment to this hair
It is bright to be described in detail, those skilled in the art should understand that, it can modify to technical solution of the present invention
Or equivalent replacement should all cover without departing from the objective and range of technical solution of the present invention in claim of the invention
In range.Technology not described in detail in the present invention, shape, construction portion are well-known technique.
Claims (3)
1. the calculation method of difference accumulation algorithm interval line number in fiber-optic vibration detection system, for using difference accumulation algorithm
Extract vibration signal when determine that two rows subtract each other the interval line number N between data, by fiber-optic vibration detect in vibration signal be considered as
One sinusoidal signal, the sinusoidal signal of a cycle are 360 ° characterized by comprising
Determine the dominant frequency period of vibration signal；
Determine that the detection cycle of degree and light pulse is subtracted each other at the interval in difference accumulation algorithm；
Subtract each other degree according to interval to be split the sinusoidal signal of a cycle, obtains segmentation number；
Vibration signal is divided equally according to segmentation number, the time required for obtaining vibration difference once；
Division arithmetic is carried out with the detection cycle of vibration difference primary required time and light pulse, obtains required light pulse
Number；
Required light pulse number is adjusted, the interval line number N of difference accumulation algorithm can be obtained.
2. the calculation method of difference accumulation algorithm interval line number in fiber-optic vibration detection system according to claim 1,
Be characterized in that: it is one in 5 °, 6 °, 7 ° or 8 ° that degree is subtracted each other at the interval.
3. the calculation method of difference accumulation algorithm interval line number in fiber-optic vibration detection system according to claim 2,
It is characterized in that: referring to time loss of the vibration signal in acquisition, processing and transmission process, required light pulse number being carried out
Adjustment, to obtain the interval line number N of difference accumulation algorithm.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910396840.7A CN110110278B (en) | 2019-05-14 | 2019-05-14 | Method for calculating number of interval lines of differential accumulation algorithm in optical fiber vibration detection system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910396840.7A CN110110278B (en) | 2019-05-14 | 2019-05-14 | Method for calculating number of interval lines of differential accumulation algorithm in optical fiber vibration detection system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110110278A true CN110110278A (en) | 2019-08-09 |
| CN110110278B CN110110278B (en) | 2022-11-11 |
Family
ID=67489949
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910396840.7A Active CN110110278B (en) | 2019-05-14 | 2019-05-14 | Method for calculating number of interval lines of differential accumulation algorithm in optical fiber vibration detection system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110110278B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0882573A (en) * | 1994-09-12 | 1996-03-26 | Shin Etsu Chem Co Ltd | Optical fiber testing method |
| US20050099634A1 (en) * | 2003-11-07 | 2005-05-12 | Marc Dubois | System and method to reduce laser noise for improved interferometric laser ultrasound detection |
| CN105681240A (en) * | 2015-12-25 | 2016-06-15 | 上海物联网有限公司 | High-precision frequency offset estimation method suitable for low signal-to-noise ratio environment |
| US20160275788A1 (en) * | 2015-03-16 | 2016-09-22 | University Of Electronic Science And Technology Of China | Online traffic volume monitoring system and method based on phase-sensitive optical time domain reflectometry |
| US20170122798A1 (en) * | 2014-04-10 | 2017-05-04 | Korea Research Institute Of Standards And Science | Method for measuring vibration displacement using state variation principle |
| CN107015765A (en) * | 2017-03-22 | 2017-08-04 | 北京石油化工学院 | A kind of method of data acquisition and procession in distributed fiberoptic sensor |
| CN108645500A (en) * | 2018-06-11 | 2018-10-12 | 桂林电子科技大学 | Method and device for vehicle vibration sensing sensitivity enhancement |
-
2019
- 2019-05-14 CN CN201910396840.7A patent/CN110110278B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0882573A (en) * | 1994-09-12 | 1996-03-26 | Shin Etsu Chem Co Ltd | Optical fiber testing method |
| US20050099634A1 (en) * | 2003-11-07 | 2005-05-12 | Marc Dubois | System and method to reduce laser noise for improved interferometric laser ultrasound detection |
| US20170122798A1 (en) * | 2014-04-10 | 2017-05-04 | Korea Research Institute Of Standards And Science | Method for measuring vibration displacement using state variation principle |
| US20160275788A1 (en) * | 2015-03-16 | 2016-09-22 | University Of Electronic Science And Technology Of China | Online traffic volume monitoring system and method based on phase-sensitive optical time domain reflectometry |
| CN105681240A (en) * | 2015-12-25 | 2016-06-15 | 上海物联网有限公司 | High-precision frequency offset estimation method suitable for low signal-to-noise ratio environment |
| CN107015765A (en) * | 2017-03-22 | 2017-08-04 | 北京石油化工学院 | A kind of method of data acquisition and procession in distributed fiberoptic sensor |
| CN108645500A (en) * | 2018-06-11 | 2018-10-12 | 桂林电子科技大学 | Method and device for vehicle vibration sensing sensitivity enhancement |
Non-Patent Citations (10)
| Title |
|---|
| HUIJUAN WU 等: "Multi-Point Disturbance Detection and High-Precision Positioning of Polarization-Sensitive Optical Time-Domain Reflectometry", 《JOURNAL OF LIGHTWAVE TECHNOLOGY》 * |
| ZHEN ZHONG 等: "Dynamic Measurement Based on the Linear Characteristic of Phase Change in --OTDR", 《IEEE PHOTONICS TECHNOLOGY LETTERS》 * |
| 尚静等: "光纤分布式扰动传感器信号检测系统研究", 《电子技术应用》 * |
| 尚静等: "小波阈值降噪法在Φ-OTDR扰动传感系统中的应用", 《现代电子技术》 * |
| 李常等: "分布式光纤振动检测系统信号后处理算法研究", 《山东科学》 * |
| 熊显名 等: "φ-OTDR全分布式光纤振动传感系统的车辆识别方法改进", 《激光杂志》 * |
| 罗俊: "提高分布式光纤Ф-OTDR安防系统性能的研究", 《中国优秀博硕士学位论文全文数据库（硕士） 工程科技I辑》 * |
| 贺梦婷 等: "基于IQ解调的相位敏感OTDR的研究", 《光通信技术》 * |
| 韩俊玲: "光学探针的近场光波振幅和相位分布的时域有限差分算法研究", 《中国优秀博硕士学位论文全文数据库（硕士） 基础科学辑》 * |
| 魏源 等: "一种减小Φ-OTDR光纤扰动传感系统中频率漂移误差的新定位算法", 《半导体光电》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110110278B (en) | 2022-11-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5103402A (en) | Method and apparatus for identifying, saving, and analyzing continuous frequency domain data in a spectrum analyzer | |
| Sedlak et al. | New automatic localization technique of acoustic emission signals in thin metal plates | |
| CN101806832B (en) | Measuring method for frequencies of low-frequency signals | |
| CN112985639B (en) | Distributed optical fiber sensing method based on chirp removal and time domain sub-chirp pulse extraction | |
| CN102043091B (en) | Digitized high-precision phase detector | |
| Sun et al. | Target location method for pipeline pre-warning system based on HHT and time difference of arrival | |
| CN106405230A (en) | Frequency measuring method and apparatus | |
| CN105865616A (en) | Modulation spectrum fast zoom method based on FFT | |
| CN104089699A (en) | Substation equipment sound reconstruction algorithm | |
| CN106845334A (en) | A kind of innovative noise extracting method based on mathematical morphology | |
| CN110672996A (en) | Urban cable multidimensional data integrated live detection device | |
| JP5237939B2 (en) | Method for instantaneous determination of signal distortion rate in AC distribution network and related apparatus | |
| CN112362304A (en) | Method for identifying target optical cable in multiple optical cables and corresponding system | |
| CN109009087A (en) | A kind of rapid detection method of R wave of electrocardiosignal | |
| CN110110278A (en) | The calculation method of difference accumulation algorithm interval line number in fiber-optic vibration detection system | |
| CN110703336B (en) | Cable positioning method based on weak magnetic detection technology | |
| CN106970265B (en) | A method of harmonic parameters are estimated using the incomplete S-transformation of Multiple Time Scales | |
| CN101281115B (en) | Correlation analysis recognition method of detection signal in closed spacing corrosion | |
| CN113848256B (en) | Ultrasonic first arrival wave real-time detection method | |
| JP3685991B2 (en) | Signal measuring method and position measuring method having noise removal function | |
| CN106872777B (en) | Harmonic and inter-harmonic separation analysis method | |
| Becker et al. | Detection of repetitive electromagnetic signals | |
| CN109239477A (en) | The detection method and detection system of electricity piece | |
| CN202661619U (en) | Nuclear magnetic resonance signal processing device | |
| CN106501603A (en) | A kind of harmonic analysis method and its device for power system |
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 |