CN104864894A - Brillouin scattering spectrum peak searching method for BOTDR (Brillouin optical time domain reflectometer) - Google Patents
Brillouin scattering spectrum peak searching method for BOTDR (Brillouin optical time domain reflectometer) Download PDFInfo
- Publication number
- CN104864894A CN104864894A CN201510178877.4A CN201510178877A CN104864894A CN 104864894 A CN104864894 A CN 104864894A CN 201510178877 A CN201510178877 A CN 201510178877A CN 104864894 A CN104864894 A CN 104864894A
- Authority
- CN
- China
- Prior art keywords
- goes
- brillouin spectrum
- data
- brillouin
- botdr
- 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
Landscapes
- Optical Communication System (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
The invention discloses a Brillouin scattering spectrum peak searching method for a BOTDR (Brillouin optical time domain reflectometer). According to the method disclosed by the invention, noise grade judgment is carried out on a Brillouin spectrum according to a generation principle of Brillouin back scattering signals, and signal enhancement methods such as filtering, cross-correlation and the like are carried out on Brillouin spectrum signals according to different noise grades, thereby reducing influences imposed on data by noises under a low signal-to-noise ratio condition, and improving the computational efficiency under a high signal-to-noise ratio condition; primary peak searching is carried out by using filtered and enhanced Brillouin spectrum data, the peak searching precision is improved, the peak searching accuracy is increased, and two times of screening are carried out by using original data by being combined with a relation between the pulse width PW and a half-peak bandwidth parameter BW, thereby further reducing the peak searching error, and improving the peak searching precision. The method disclosed by the invention improves the accuracy of a BOTDR product in testing a multi-peak optical fiber and the space resolution of a single-peak optical fiber at a large strain condition under the premise of not changing hardware, and further expands an application field of the BOTDR product.
Description
Technical field
The present invention relates to a kind of BOTDR Brillouin spectrum Peak Search Method.
Background technology
Brillouin light domain reflectometer (BOTDR) relies on the Brillouin shift distribution of backward Brillouin scattering light in measuring optical fiber to calculate the Strain Distribution of optical fiber, pulsed light is incident from optical fiber one end with certain frequency, produce Brillouin scattering after phonon in incident pulsed light and optical fiber interacts, backward Brillouin scattering light turns back to incidence end along the former road of optical fiber.Due to Brillouin scattering optical frequency shift in optical fiber and there is linear relationship between optical fiber axial strain and temperature, the Brillouin frequency shifts therefore measuring optical fiber distributes the Strain Distribution that can calculate in optical fiber.The characteristics such as Brillouin light domain reflectometer has low energy dependence, high environmental resistance, electromagnetism interference, anticorrosive, waterproof, moisture resistant, thermal adaptability are strong, and owing to can single-endedly measure, difficulty of construction is lower and extensively concerned.Brillouin light domain reflectometer can be used for the field such as health monitoring of Geotechnical Engineering health monitoring, geo-hazard early-warning monitoring, cable and pipeline, is that engineering field is for one of the strongest product replacing conventional port sensor.
Because optical fiber material, technique all there are differences, when using BOTDR product to carry out strain detecting, there is the Brillouin spectrum of multiple optical fiber all can there is more than one crest, each crest includes respective strain/temperature information, if crest position each in optical fiber can be searched out accurately, just strain/temperature information that each crest contains can be extracted, just can mutually correct, will effectively promote the strain testing precision of BOTDR product; For some special optical fibers, in Brillouin spectrum, each crest also exists the response coefficient of respective strain/temperature, the position of each crest is accurately found out, and simultaneous analysis is carried out to its strain/temperature information, strain/temperature test function simultaneously of BOTDR product can be realized; Simultaneously, Brillouin spectrum is only existed to the optical fiber of a crest, when the spatial resolution of BOTDR product cannot meet the demands, if the strained increase of tested optical fiber, also tested optical fiber can be caused to test and to occur the test peak that another intensity is lower, now find the position at the Brillouin spectrum peak that intensity is lower accurately, optical fiber institute can be calculated when spatial resolution deficiency strained, effectively can improve the spatial resolution of BOTDR product.
The peak-seeking mode of existing BOTDR product is as Brillouin spectrum crest location using point the highest for intensity in Brillouin spectrum, unimodal searching can only be realized, the crest that in Brillouin spectrum, other intensity is lower cannot be found out, analyzed, and causes the strain information in other crest cannot be used for the lifting of spatial resolution or the lifting of strain precision; And when Brillouin spectrum signal is lower, even with noise intensity close to time, it is Brillouin spectrum peak that traditional unimodal search scheme is easy to noise identification, and the strain data calculated with this also exists great error, have impact on the performance boost of BOTDR product, propagation and employment.
Summary of the invention
For the above-mentioned technical matters existed in prior art, the present invention proposes a kind of BOTDR Brillouin spectrum Peak Search Method, it adopts following technical scheme:
A kind of BOTDR Brillouin spectrum Peak Search Method, comprises the steps:
Step 101, reading Brillouin spectrum data points N and Brillouin spectrum data D [0] ~ D [N-1] frequently;
Step 102, the maximal value DMAX calculated in Brillouin spectrum data D [0] ~ D [N-1];
Step 103, judge Brillouin spectrum data noise grade SL;
Step 104, points being less than noise threshold NT all in Brillouin spectrum data D are set to 0, generate interim Brillouin spectrum data DL [0] ~ DL [N-1];
If step 105 noise grade SL≤2, perform step 106, otherwise perform step 107;
Step 106, data enhancing is carried out to Brillouin spectrum data D [0] ~ D [N-1], the Brillouin spectrum that is enhanced data DZ [0] ~ DZ [N-1];
Step 107, copy Brillouin spectrum data D [0] ~ D [N-1], the Brillouin spectrum that is enhanced data DZ [0] ~ DZ [N-1];
Step 108, foundation strengthen Brillouin spectrum crest flag data BSZ [0] ~ BSZ [N-1], interim Brillouin spectrum crest flag data BSL [0] ~ BSL [N-1], Brillouin spectrum crest flag data BS [0] ~ BS [N-1], and are all initialized as 0;
Step 109, the preliminary wave crest point found in enhancing Brillouin spectrum data DZ;
Step 110, enhancing Brillouin spectrum data DZ all points being labeled as crest to be screened for the first time;
Step 111, reading pulse width parameter PW;
Step 112, according to pulse width PW, calculate half-peak bandwidth reference parameter BW;
Step 113, reading frequency spacing parameter FG;
Step 114, postsearch screening is carried out to the point that all values in interim Brillouin spectrum crest flag data BSL is 1;
Step 115, the value of BSL [I] be corresponding in turn to be assigned to BS [I], 0 < I < N-1, export Brillouin spectrum crest flag data BS, BS data intermediate value be 1 point be crest.
Further, in above-mentioned steps 103, Brillouin spectrum data noise grade SL decision method is as follows:
If step 1031 DMAX≤50dB, goes to step 1032, otherwise goes to step 1033;
Step 1032, noise grade SL are 1, noise threshold NT is 40dB, goes to step 1038;
If step 1,033 50 < is DMAX≤60dB, goes to step 1034, otherwise go to step 1035;
Step 1034, noise grade SL are 2, noise threshold NT is 45dB, goes to step 1038;
If step 1,035 60 < is DMAX≤70dB, goes to step 1036, otherwise go to step 1037;
Step 1036, noise grade SL are 3, noise threshold NT is 50dB, goes to step 1038;
Step 1037, noise grade SL are 4, noise threshold NT is 55dB;
Step 1038, Brillouin spectrum data noise grade SL judge to terminate.
Further, in above-mentioned steps 109, first secondary wave crest decision method is as follows:
Step 1091, initialization I=1,0 < I < N-1;
If step 1092 DZ [I] is the maximal value in DZ [I-2] ~ DZ [I+2], goes to step 1093, otherwise go to step 1094;
Step 1093, be 1 by BSZ [I] assignment;
Step 1094, by the value of I increase by 1;
If step 1095 I >=N-1, goes to step 1096, otherwise, go to step 1092;
Step 1096, first secondary wave crest judge to terminate.
Further, in above-mentioned steps 110, first screening decision method is as follows:
Step 1101, initialization I=1,0 < I < N-1;
If step 1102 BSZ [I] is 1, goes to step 1103, otherwise go to step 1104;
If step 1103 DL [I] is maximal value in DL [I-2] ~ DL [I+2], goes to step 1105, otherwise go to step 1104;
Step 1104, be 0 by BSL [I] assignment, go to step 1106;
Step 1105, be 1 by BSL [I] assignment;
Step 1106, by the value of I increase by 1;
If step 1107 I >=N-1, goes to step 1108, otherwise, go to step 1103;
Step 1108, first screening judge to terminate.
Further, in above-mentioned steps 114, postsearch screening decision method is as follows:
Step 1141, initialization I=1,0 < I < N-1;
If step 1142 BSZ [I] is 1, goes to step 1143, otherwise go to step 1144;
Step 1143, calculating IL;
Step 1144, be 1 by BSL [I] assignment, go to step 1148;
Step 1145, calculating IR;
If step 1146 IR-I≤BW/2 and I-IL≤BW/2, go to step 1147, otherwise go to step 1144;
Step 1147, be 0 by BSL [I] assignment, go to step 1148;
Step 1148, by the value of I increase by 1;
If step 1149 I >=N-1, goes to step 1108, otherwise, go to step 1103;
Step 11410, postsearch screening judge to terminate.
Further, in above-mentioned steps 1143, the computation process of IL is as follows:
The value of step 11431, initialization J is I-1;
If step 11432 J >=0, goes to step 11433, otherwise goes to step 11434;
If step 11433 DL [J]≤1dB, goes to step 11436, otherwise goes to step 11435;
Step 11434, be 0 by IL assignment, go to step 11437;
Step 11435, the value of I is subtracted 1, go to step 11432;
Step 11436, be J by IL assignment;
Step 11437, calculating IL process terminate.
Further, in above-mentioned steps 1145, the computation process of IR is as follows:
The value of step 11451, initialization J is I+1;
If step 11452 J≤N-1, goes to step 11453, otherwise goes to step 11454;
If step 11453 DL [J]≤1dB, goes to step 11456, otherwise goes to step 11455;
Step 11454, be N-1 by IR assignment, go to step 11457;
Step 11455, the value of I is subtracted 1, go to step 11452;
Step 11456, be J by IR assignment;
Step 11457, calculating IL process terminate.
Tool of the present invention has the following advantages:
The inventive method is according to the generation principle of Brillouin's backscatter signal, noise grade judgement is carried out to Brillouin spectrum, and according to noise grade difference, the signal such as filtering, cross-correlation Enhancement Method is carried out to Brillouin spectrum signal, under reducing Low SNR, noise is on the impact of data, improves counting yield under high s/n ratio condition; The Brillouin spectrum data after filtering, enhancing are utilized to carry out first crest seeking, improve crest seeking precision, add crest seeking accuracy, and in conjunction with the relation of pulse width PW and half-peak bandwidth reference parameter BW, raw data is used to carry out twice screening, reduce further peak-seeking error, improve peak-seeking precision.The inventive method is under the prerequisite of not changing hardware, and accuracy when improving BOTDR product test multimodal optical fiber and large sstrain condition place an order the spatial resolution of peak fiber, expand the application of BOTDR product further.
Accompanying drawing explanation
Fig. 1 is the schematic flow sheet of a kind of BOTDR Brillouin spectrum Peak Search Method in the present invention;
Fig. 2 is Brillouin spectrum data noise grade SL decision method schematic diagram in the present invention;
Fig. 3 is that in the present invention, first peak-seeking judges schematic diagram;
Fig. 4 is that in the present invention, first screening judges schematic diagram;
Fig. 5 is that in the present invention, postsearch screening judges schematic diagram;
Fig. 6 is that in the present invention, IL calculates schematic diagram;
Fig. 7 is that in the present invention, IR calculates schematic diagram.
Embodiment
Below in conjunction with accompanying drawing and embodiment, the present invention is described in further detail:
As shown in Figure 1, a kind of BOTDR Brillouin spectrum Peak Search Method, comprises the steps:
Step 101, reading Brillouin spectrum data points N and Brillouin spectrum data D [0] ~ D [N-1] frequently.
Step 102, the maximal value DMAX calculated in Brillouin spectrum data D [0] ~ D [N-1].
Step 103, judge Brillouin spectrum data noise grade SL.
As shown in Figure 2, Brillouin spectrum data noise grade SL decision method is as follows:
If step 1031 DMAX≤50dB, goes to step 1032, otherwise goes to step 1033;
Step 1032, noise grade SL are 1, noise threshold NT is 40dB, goes to step 1038;
If step 1,033 50 < is DMAX≤60dB, goes to step 1034, otherwise go to step 1035;
Step 1034, noise grade SL are 2, noise threshold NT is 45dB, goes to step 1038;
If step 1,035 60 < is DMAX≤70dB, goes to step 1036, otherwise go to step 1037;
Step 1036, noise grade SL are 3, noise threshold NT is 50dB, goes to step 1038;
Step 1037, noise grade SL are 4, noise threshold NT is 55dB;
Step 1038, Brillouin spectrum data noise grade SL judge to terminate.
Step 104, points being less than noise threshold NT all in Brillouin spectrum data D are set to 0, generate interim Brillouin spectrum data DL [0] ~ DL [N-1].
If step 105 noise grade SL≤2, perform step 106, otherwise perform step 107.
Step 106, data enhancing is carried out to Brillouin spectrum data D [0] ~ D [N-1], the Brillouin spectrum that is enhanced data DZ [0] ~ DZ [N-1];
Operable signal Enhancement Method comprises the methods such as cross-correlation method, average filter, medium filtering, wavelet filtering.
Step 107, copy Brillouin spectrum data D [0] ~ D [N-1], the Brillouin spectrum that is enhanced data DZ [0] ~ DZ [N-1].
Step 108, foundation strengthen Brillouin spectrum crest flag data BSZ [0] ~ BSZ [N-1], interim Brillouin spectrum crest flag data BSL [0] ~ BSL [N-1], Brillouin spectrum crest flag data BS [0] ~ BS [N-1], and are all initialized as 0.
Step 109, the preliminary wave crest point found in enhancing Brillouin spectrum data DZ.
Shown in composition graphs 3, first peak-seeking decision principle is: traversal strengthens Brillouin spectrum data DZ, as the maximal value that data DZ [I] (0 < I < N-1) is in DZ [I-2 ~ I+2] five points, be 1 by BSZ [I] assignment, otherwise BSZ [I] assignment is 0.
Concrete decision method is as follows:
Step 1091, initialization I=1,0 < I < N-1;
If step 1092 DZ [I] is the maximal value in DZ [I-2] ~ DZ [I+2], goes to step 1093, otherwise go to step 1094;
Step 1093, be 1 by BSZ [I] assignment;
Step 1094, by the value of I increase by 1;
If step 1095 I >=N-1, goes to step 1096, otherwise, go to step 1092;
Step 1096, first secondary wave crest judge to terminate.
Step 110, enhancing Brillouin spectrum data DZ all points being labeled as crest to be screened for the first time.
As shown in Figure 4, first screening principle is: if BSZ [I] is 1, then judge interim Brillouin spectrum data DL [I] whether as the maximal value in DL [I-2] ~ DL [I+2] five points, if DL [I] is maximal value, then BSL [I] assignment is 1, otherwise BSL [I] assignment is 0.
Concrete screening technique is as follows:
Step 1101, initialization I=1,0 < I < N-1;
If step 1102 BSZ [I] is 1, goes to step 1103, otherwise go to step 1104;
If step 1103 DL [I] is maximal value in DL [I-2] ~ DL [I+2], goes to step 1105, otherwise go to step 1104;
Step 1104, be 0 by BSL [I] assignment, go to step 1106;
Step 1105, be 1 by BSL [I] assignment;
Step 1106, by the value of I increase by 1;
If step 1107 I >=N-1, goes to step 1108, otherwise, go to step 1103;
Step 1108, first screening judge to terminate.
Step 111, reading pulse width parameter PW.
Step 112, according to pulse width PW, calculate half-peak bandwidth reference parameter BW, half envelope bandwidth reference parameter BW calculating principle is in table 1.
Table 1 pulse width PW and half-peak bandwidth reference parameter BW relation table
| Pulse width PW | Half-peak bandwidth reference parameter BW |
| 5ns | 200MHz |
| 10ns | 100MHz |
| 20ns | 60MHz |
| 50ns and more than 50ns | 40MHz |
Step 113, reading frequency spacing parameter FG.
Step 114, postsearch screening is carried out to the point that all values in interim Brillouin spectrum crest flag data BSL is 1.
As shown in Figure 5, postsearch screening principle is as follows: if BSL [I] is 1, then find some D [IL] <D [I]-3 (IL<I), and the point that on the left of I, first value is less than D [I]-3 is IL; Find some D [IR] <D [I]-3 (IR>I), the point that on the right side of I, first value is less than D [I]-3 is IR, as IR-I<BW/2 and I-IR<BW/2, then BSL [I] assignment is 0.
Concrete screening step is as follows:
Step 1141, initialization I=1,0 < I < N-1;
If step 1142 BSZ [I] is 1, goes to step 1143, otherwise go to step 1144;
Step 1143, calculating IL;
Step 1144, be 1 by BSL [I] assignment, go to step 1148;
Step 1145, calculating IR;
If step 1146 IR-I≤BW/2 and I-IL≤BW/2, go to step 1147, otherwise go to step 1144;
Step 1147, be 0 by BSL [I] assignment, go to step 1148;
Step 1148, by the value of I increase by 1;
If step 1149 I >=N-1, goes to step 1108, otherwise, go to step 1103;
Step 11410, postsearch screening judge to terminate.
Wherein, in above-mentioned steps 1143, the computation process of IL as shown in Figure 6, comprises the steps:
The value of step 11431, initialization J is I-1;
If step 11432 J >=0, goes to step 11433, otherwise goes to step 11434;
If step 11433 DL [J]≤1dB, goes to step 11436, otherwise goes to step 11435;
Step 11434, be 0 by IL assignment, go to step 11437;
Step 11435, the value of I is subtracted 1, go to step 11432;
Step 11436, be J by IL assignment;
Step 11437, calculating IL process terminate.
In above-mentioned steps 1145, the computation process of IR as shown in Figure 7, comprises the steps:
The value of step 11451, initialization J is I+1;
If step 11452 J≤N-1, goes to step 11453, otherwise goes to step 11454;
If step 11453 DL [J]≤1dB, goes to step 11456, otherwise goes to step 11455;
Step 11454, be N-1 by IR assignment, go to step 11457;
Step 11455, the value of I is subtracted 1, go to step 11452;
Step 11456, be J by IR assignment;
Step 11457, calculating IL process terminate.
Step 115, the value of BSL [I] be corresponding in turn to be assigned to BS [I], 0 < I < N-1, export Brillouin spectrum crest flag data BS, BS data intermediate value be 1 point be crest.
Certainly; more than illustrate and be only preferred embodiment of the present invention; the present invention is not limited to enumerate above-described embodiment; should be noted that; any those of ordinary skill in the art are under the instruction of this instructions; made all equivalently to substitute, obvious form of distortion, within the essential scope all dropping on this instructions, protection of the present invention ought to be subject to.
Claims (7)
1. a BOTDR Brillouin spectrum Peak Search Method, is characterized in that, comprises the steps:
Step 101, reading Brillouin spectrum data points N and Brillouin spectrum data D [0] ~ D [N-1] frequently;
Step 102, the maximal value DMAX calculated in Brillouin spectrum data D [0] ~ D [N-1];
Step 103, judge Brillouin spectrum data noise grade SL;
Step 104, points being less than noise threshold NT all in Brillouin spectrum data D are set to 0, generate interim Brillouin spectrum data DL [0] ~ DL [N-1];
If step 105 noise grade SL≤2, perform step 106, otherwise perform step 107;
Step 106, data enhancing is carried out to Brillouin spectrum data D [0] ~ D [N-1], the Brillouin spectrum that is enhanced data DZ [0] ~ DZ [N-1];
Step 107, copy Brillouin spectrum data D [0] ~ D [N-1], the Brillouin spectrum that is enhanced data DZ [0] ~ DZ [N-1];
Step 108, foundation strengthen Brillouin spectrum crest flag data BSZ [0] ~ BSZ [N-1], interim Brillouin spectrum crest flag data BSL [0] ~ BSL [N-1], Brillouin spectrum crest flag data BS [0] ~ BS [N-1], and are all initialized as 0;
Step 109, the preliminary wave crest point found in enhancing Brillouin spectrum data DZ;
Step 110, enhancing Brillouin spectrum data DZ all points being labeled as crest to be screened for the first time;
Step 111, reading pulse width parameter PW;
Step 112, according to pulse width PW, calculate half-peak bandwidth reference parameter BW;
Step 113, reading frequency spacing parameter FG;
Step 114, postsearch screening is carried out to the point that all values in interim Brillouin spectrum crest flag data BSL is 1;
Step 115, the value of BSL [I] be corresponding in turn to be assigned to BS [I], 0 < I < N-1, export Brillouin spectrum crest flag data BS, BS data intermediate value be 1 point be crest.
2. a kind of BOTDR Brillouin spectrum Peak Search Method according to claim 1, is characterized in that, in described step 103, Brillouin spectrum data noise grade SL decision method is as follows:
If step 1031 DMAX≤50dB, goes to step 1032, otherwise goes to step 1033;
Step 1032, noise grade SL are 1, noise threshold NT is 40dB, goes to step 1038;
If step 1,033 50 < is DMAX≤60dB, goes to step 1034, otherwise go to step 1035;
Step 1034, noise grade SL are 2, noise threshold NT is 45dB, goes to step 1038;
If step 1,035 60 < is DMAX≤70dB, goes to step 1036, otherwise go to step 1037;
Step 1036, noise grade SL are 3, noise threshold NT is 50dB, goes to step 1038;
Step 1037, noise grade SL are 4, noise threshold NT is 55dB;
Step 1038, Brillouin spectrum data noise grade SL judge to terminate.
3. a kind of BOTDR Brillouin spectrum Peak Search Method according to claim 1, is characterized in that, in described step 109, first secondary wave crest decision method is as follows:
Step 1091, initialization I=1,0 < I < N-1;
If step 1092 DZ [I] is the maximal value in DZ [I-2] ~ DZ [I+2], goes to step 1093, otherwise go to step 1094;
Step 1093, be 1 by BSZ [I] assignment;
Step 1094, by the value of I increase by 1;
If step 1095 I >=N-1, goes to step 1096, otherwise, go to step 1092;
Step 1096, first secondary wave crest judge to terminate.
4. a kind of BOTDR Brillouin spectrum Peak Search Method according to claim 1, is characterized in that, in described step 110, first screening decision method is as follows:
Step 1101, initialization I=1,0 < I < N-1;
If step 1102 BSZ [I] is 1, goes to step 1103, otherwise go to step 1104;
If step 1103 DL [I] is maximal value in DL [I-2] ~ DL [I+2], goes to step 1105, otherwise go to step 1104;
Step 1104, be 0 by BSL [I] assignment, go to step 1106;
Step 1105, be 1 by BSL [I] assignment;
Step 1106, by the value of I increase by 1;
If step 1107 I >=N-1, goes to step 1108, otherwise, go to step 1103;
Step 1108, first screening judge to terminate.
5. a kind of BOTDR Brillouin spectrum Peak Search Method according to claim 1, it is characterized in that, in described step 114, postsearch screening decision method is as follows:
Step 1141, initialization I=1,0 < I < N-1;
If step 1142 BSZ [I] is 1, goes to step 1143, otherwise go to step 1144;
Step 1143, calculating IL;
Step 1144, be 1 by BSL [I] assignment, go to step 1148;
Step 1145, calculating IR;
If step 1146 IR-I≤BW/2 and I-IL≤BW/2, go to step 1147, otherwise go to step 1144;
Step 1147, be 0 by BSL [I] assignment, go to step 1148;
Step 1148, by the value of I increase by 1;
If step 1149 I >=N-1, goes to step 1108, otherwise, go to step 1103;
Step 11410, postsearch screening judge to terminate.
6. a kind of BOTDR Brillouin spectrum Peak Search Method according to claim 5, it is characterized in that, in described step 1143, the computation process of IL is as follows:
The value of step 11431, initialization J is I-1;
If step 11432 J >=0, goes to step 11433, otherwise goes to step 11434;
If step 11433 DL [J]≤1dB, goes to step 11436, otherwise goes to step 11435;
Step 11434, be 0 by IL assignment, go to step 11437;
Step 11435, the value of I is subtracted 1, go to step 11432;
Step 11436, be J by IL assignment;
Step 11437, calculating IL process terminate.
7. a kind of BOTDR Brillouin spectrum Peak Search Method according to claim 5, it is characterized in that, in described step 1145, the computation process of IR is as follows:
The value of step 11451, initialization J is I+1;
If step 11452 J≤N-1, goes to step 11453, otherwise goes to step 11454;
If step 11453 DL [J]≤1dB, goes to step 11456, otherwise goes to step 11455;
Step 11454, be N-1 by IR assignment, go to step 11457;
Step 11455, the value of I is subtracted 1, go to step 11452;
Step 11456, be J by IR assignment;
Step 11457, calculating IL process terminate.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510178877.4A CN104864894B (en) | 2015-04-16 | 2015-04-16 | A kind of BOTDR Brillouin spectrums Peak Search Method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510178877.4A CN104864894B (en) | 2015-04-16 | 2015-04-16 | A kind of BOTDR Brillouin spectrums Peak Search Method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104864894A true CN104864894A (en) | 2015-08-26 |
| CN104864894B CN104864894B (en) | 2017-10-31 |
Family
ID=53910891
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510178877.4A Active CN104864894B (en) | 2015-04-16 | 2015-04-16 | A kind of BOTDR Brillouin spectrums Peak Search Method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104864894B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106017522A (en) * | 2016-05-11 | 2016-10-12 | 武汉理工大学 | Rapid high-precision signal demodulation method of fiber F-P sensor |
| CN107560645A (en) * | 2017-08-29 | 2018-01-09 | 北京航空航天大学 | A kind of fiber Bragg grating sensor Wavelength demodulation Peak Search Method |
| CN107907507A (en) * | 2017-09-28 | 2018-04-13 | 南昌航空大学 | A kind of lubricating oil elasticity modulus detection method based on Brillouin scattering spectrum |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS63191950A (en) * | 1987-02-04 | 1988-08-09 | Nippon X-Ray Kk | Peak position detecting method for spectrum signal |
| JP2009281792A (en) * | 2008-05-20 | 2009-12-03 | Dai Ichi High Frequency Co Ltd | Fbg light spectrum analyzer |
| CN102098101A (en) * | 2011-01-04 | 2011-06-15 | 武汉理工大学 | Crest searching method and system |
| CN102435214A (en) * | 2011-09-20 | 2012-05-02 | 华北电力大学（保定） | High-precision and fast optical fiber grating spectrum peak-searching method |
| CN102706372A (en) * | 2012-03-28 | 2012-10-03 | 中国航空工业集团公司北京长城计量测试技术研究所 | Optical fiber wavelength demodulating light spectrum peak value positioning method |
| CN102944255A (en) * | 2012-11-20 | 2013-02-27 | 武汉理工光科股份有限公司 | Instrument and method for demodulating digital fiber bragg grating |
| CN103487074A (en) * | 2013-10-12 | 2014-01-01 | 重庆邮电大学 | Method for processing FBG (fiber bragg grating) sensing signal by utilizing three-point peek-seeking algorithm |
| CN104316089A (en) * | 2014-10-27 | 2015-01-28 | 中国电子科技集团公司第四十一研究所 | Method for enhancing BOTDR (Brillouin Optical Time Domain Reflectometer) spatial resolution |
-
2015
- 2015-04-16 CN CN201510178877.4A patent/CN104864894B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS63191950A (en) * | 1987-02-04 | 1988-08-09 | Nippon X-Ray Kk | Peak position detecting method for spectrum signal |
| JP2009281792A (en) * | 2008-05-20 | 2009-12-03 | Dai Ichi High Frequency Co Ltd | Fbg light spectrum analyzer |
| CN102098101A (en) * | 2011-01-04 | 2011-06-15 | 武汉理工大学 | Crest searching method and system |
| CN102435214A (en) * | 2011-09-20 | 2012-05-02 | 华北电力大学（保定） | High-precision and fast optical fiber grating spectrum peak-searching method |
| CN102706372A (en) * | 2012-03-28 | 2012-10-03 | 中国航空工业集团公司北京长城计量测试技术研究所 | Optical fiber wavelength demodulating light spectrum peak value positioning method |
| CN102944255A (en) * | 2012-11-20 | 2013-02-27 | 武汉理工光科股份有限公司 | Instrument and method for demodulating digital fiber bragg grating |
| CN103487074A (en) * | 2013-10-12 | 2014-01-01 | 重庆邮电大学 | Method for processing FBG (fiber bragg grating) sensing signal by utilizing three-point peek-seeking algorithm |
| CN104316089A (en) * | 2014-10-27 | 2015-01-28 | 中国电子科技集团公司第四十一研究所 | Method for enhancing BOTDR (Brillouin Optical Time Domain Reflectometer) spatial resolution |
Non-Patent Citations (2)
| Title |
|---|
| 刘泉等: "高速光纤光栅解调的寻峰算法研究", 《武汉理工大学学报》 * |
| 尚秋峰等: "光纤Bragg光栅传感***典型寻峰算法的比较分析", 《电测与仪表》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106017522A (en) * | 2016-05-11 | 2016-10-12 | 武汉理工大学 | Rapid high-precision signal demodulation method of fiber F-P sensor |
| CN107560645A (en) * | 2017-08-29 | 2018-01-09 | 北京航空航天大学 | A kind of fiber Bragg grating sensor Wavelength demodulation Peak Search Method |
| CN107907507A (en) * | 2017-09-28 | 2018-04-13 | 南昌航空大学 | A kind of lubricating oil elasticity modulus detection method based on Brillouin scattering spectrum |
| CN107907507B (en) * | 2017-09-28 | 2019-11-15 | 南昌航空大学 | A kind of lubricating oil elasticity modulus detection method based on Brillouin scattering spectrum |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104864894B (en) | 2017-10-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111751392B (en) | Steel bar corrosion detection method based on dual-polarization ground penetrating radar | |
| CN103063737B (en) | Coiled tubing magnetic detection method | |
| CN106525096B (en) | A kind of brillouin distributed optical fiber sensing device and reduce gain spectral line width method | |
| CN103344989B (en) | The analytical approach of impulse noise interference in vibroseis seismologic record | |
| CN104864894A (en) | Brillouin scattering spectrum peak searching method for BOTDR (Brillouin optical time domain reflectometer) | |
| CN107678064B (en) | Real-time extraction method for sound wave time difference | |
| CN105303181A (en) | Stochastic resonance weak impact feature enhancement extraction method on the basis of sliding window | |
| CN106932337A (en) | The data processing method of multi fiber grating steel reinforcement corrosion sensor system | |
| CN106546275A (en) | A kind of preparation method of Brillouin spectrum characteristic initial value | |
| CN113343402B (en) | Pipeline corrosion grade evaluation method based on multilayer convolution sparse coding | |
| CN106199349A (en) | Partial discharge of transformer region decision method based on calibration square-wave signal transmission ratio | |
| CN114674352A (en) | Distributed disturbance sensing and demodulation method based on Rayleigh scattering spectrum non-similarity | |
| CN112798142B (en) | Brillouin optical fiber sensor strain and temperature two-stage rapid identification method based on Bayesian updating and random simulation | |
| CN104296673A (en) | Brillouin spectrum signal quality improving method | |
| CN110333506B (en) | Method for extracting inhaul cable position parameters of cable force measurement radar | |
| CN112113148A (en) | Pipe network side leakage detection method based on pressure and flow for rail transit | |
| CN101033617A (en) | Dynamic detecting method for basic structure testing signal | |
| CN110455438B (en) | Automatic calibration and automatic test method for optical fiber length of optical fiber temperature distribution tester | |
| CN116839865A (en) | Long-distance distributed optical fiber performance rapid detection device and method based on distributed amplification and optical chirp chain | |
| CN102313858B (en) | Method for identifying traveling wave in initial reversed polarity direction | |
| CN106771928A (en) | A kind of online pick-up method of partial discharge pulse's initial time | |
| CN102518422A (en) | Method for detecting and identifying current stress of oil field downhole casing | |
| CN116088050A (en) | Rock burst prediction method based on microseismic monitoring of space-time strong parameters of fracture source | |
| KR102376920B1 (en) | Diagnosis method for wedge of a stator of a generator | |
| CN103291279B (en) | The optimization method of the micro-abnormal signal of a kind of Gas Logging Value |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20190227 Address after: 266555 Xiangjiang 98, Huangdao District, Qingdao City, Shandong Province Patentee after: China Electronics Technology Instrument and Meter Co., Ltd. Address before: 266555 No. 98 Xiangjiang Road, Qingdao economic and Technological Development Zone, Shandong Patentee before: The 41st Institute of CETC |
|
| CP01 | Change in the name or title of a patent holder | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 266555 Xiangjiang 98, Huangdao District, Qingdao City, Shandong Province Patentee after: CLP kesiyi Technology Co.,Ltd. Address before: 266555 Xiangjiang 98, Huangdao District, Qingdao City, Shandong Province Patentee before: CHINA ELECTRONIC TECHNOLOGY INSTRUMENTS Co.,Ltd. |