CN113639652B - Automatic testing and analyzing method for strain distribution of optical fiber hydrophone - Google Patents
Automatic testing and analyzing method for strain distribution of optical fiber hydrophone Download PDFInfo
- Publication number
- CN113639652B CN113639652B CN202111015471.6A CN202111015471A CN113639652B CN 113639652 B CN113639652 B CN 113639652B CN 202111015471 A CN202111015471 A CN 202111015471A CN 113639652 B CN113639652 B CN 113639652B
- Authority
- CN
- China
- Prior art keywords
- executing
- data
- interference arm
- judging whether
- optical fiber
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000013307 optical fiber Substances 0.000 title claims abstract description 123
- 238000000034 method Methods 0.000 title claims abstract description 55
- 238000012360 testing method Methods 0.000 title claims abstract description 43
- 230000003287 optical effect Effects 0.000 claims abstract description 30
- 230000005540 biological transmission Effects 0.000 claims abstract description 14
- 238000001228 spectrum Methods 0.000 claims description 74
- 230000008569 process Effects 0.000 claims description 31
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 239000000835 fiber Substances 0.000 description 20
- 238000004458 analytical method Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 10
- 238000004804 winding Methods 0.000 description 9
- 238000005516 engineering process Methods 0.000 description 7
- 238000011161 development Methods 0.000 description 6
- 238000001514 detection method Methods 0.000 description 4
- 230000006872 improvement Effects 0.000 description 4
- 238000007689 inspection Methods 0.000 description 4
- 230000003595 spectral effect Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000007405 data analysis Methods 0.000 description 2
- 230000002706 hydrostatic effect Effects 0.000 description 2
- 239000003550 marker Substances 0.000 description 2
- 230000010287 polarization Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
- G01B11/161—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge by interferometric means
-
- 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/002—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means for representing acoustic field distribution
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Testing Of Optical Devices Or Fibers (AREA)
Abstract
The disclosure provides an automatic testing and analyzing method for strain distribution of an optical fiber hydrophone, which comprises the following steps: acquiring the length of a long interference arm of an optical fiber hydrophone, the length of a short interference arm of the optical fiber hydrophone and the length of a transmission part optical fiber in each optical switch channel; initializing a counting variable of an optical switch channel, and remotely setting parameters of an optical fiber strain distribution tester; strain distribution test is carried out on the optical fiber hydrophone to obtain strain distribution test data; rapidly analyzing original Brillouin distribution data in strain distribution test data to obtain the total optical fiber length, the interference arm reflector position and the long interference arm reflector position of the strain data; judging rules, and directly proceeding to the next step when the judging rules are met; analyzing the original Brillouin distribution data and the original strain data to obtain long interference arm strain data and short interference arm strain data; and judging whether the current optical switch channel counting variable is not smaller than the optical switch channel number, if so, directly outputting all channel long interference arm strain data and short interference arm strain data.
Description
Technical Field
The disclosure belongs to the technical field of underwater sound detection, and particularly relates to an automatic testing and analyzing method for strain distribution of an optical fiber hydrophone.
Background
The statements in this section merely provide background information related to the present disclosure and may not necessarily constitute prior art.
The fiber hydrophone has wide territory in the territory of the ocean and abundant ocean resources, is used as a novel underwater sound detection technology developed in recent years, and has important significance for maintaining ocean safety. Compared with the traditional piezoelectric hydrophone, the optical fiber hydrophone has the advantages of high sensitivity, full light at the wet end, good stability and frequency response characteristics, hydrostatic pressure resistance, corrosion resistance, long transmission distance, easiness in large-scale multiplexing and the like. The sensing principle is that the information such as the frequency, the intensity and the like of the sound wave is obtained by utilizing parameters such as the intensity, the polarization state, the phase and the like of the light wave in the sound wave modulation optical fiber. The phase interference type optical fiber hydrophone is the most mature technology at present and has the most research at home and abroad. The Michelson interference type optical fiber hydrophone converts the underwater acoustic signal into the phase change of the optical signal, extracts the underwater acoustic information from the phase change of the optical signal through an optical coherence detection technology and a signal processing system, and has the advantages of high sensitivity, convenience in multiplexing and the like.
The optical fiber hydrophone needs to work in a severe marine environment with high reliability for a long time, which puts high demands on the manufacturing process of the optical fiber hydrophone. The optical fiber coil is one of key devices in the optical fiber hydrophone, when the optical fiber hydrophone is disturbed by external sound pressure, the optical fiber coil deforms, so that the intensity, the polarization state, the frequency or the phase of optical waves in the optical fiber are modulated, and then a corresponding signal demodulation method is adopted to demodulate detection signals through the rear end of the hydrophone, so that external sound field information is obtained. Thus, the performance of the fiber optic coil will directly affect the optical and acoustic performance and long-term reliability of the fiber optic hydrophone.
The optical fiber coil winding process is to wind the optical fiber on the elastic body according to a certain winding method to obtain the optical fiber coil, wherein the winding tension is an important process parameter when the optical fiber coil is wound, and the winding tension is usually determined according to the performance parameters of the wound optical fiber and the elastic body material, the working water depth hydrostatic pressure of the optical fiber hydrophone and other factors. In the improvement process of the winding process of the optical fiber coil, the development process and the inspection process of the optical fiber hydrophone, the winding tension of the optical fiber coil can be accurately measured, and the method is particularly important for the improvement of the winding process, the development and the inspection of the optical fiber coil.
According to the knowledge of the inventor, the optical fiber strain distribution tester is mainly used for testing the strain distribution state of the optical fiber coil at present, but because the optical fiber coil is wound by using two interference arms of the Michelson optical fiber hydrophone, and the lengths of the two interference arms are different, when the optical fiber strain tester is used for testing, the strain data are difficult to accurately analyze due to mutual interference caused by superposition of self-Brillouin scattering signals of the two interference arms, and only the excessive length part of the longer interference arm compared with the shorter interference arm of the two interference arms can obtain more accurate strain data because the interference of the self-published Brillouin scattering signal of the other interference arm is avoided. Therefore, strain data of the overlapping part of the two interference arms cannot be obtained by adopting the prior art, only strain data of the residual length part of the longer interference arm can be obtained, the strain distribution state of the whole optical fiber coil is difficult to accurately reflect, the accurate judgment of the performance of the optical fiber coil in the improvement process of the optical fiber hydrophone winding process and the development and inspection process of the optical fiber hydrophone is influenced, and particularly, the judgment capability of the deformation state of the shorter interference arm is greatly influenced, so that the rapid development of the optical fiber hydrophone technology is severely limited.
Disclosure of Invention
In order to solve the problems, the disclosure provides an automatic testing and analyzing method for strain distribution of an optical fiber hydrophone, and solves the problem that the traditional technology cannot test the strain distribution of two interference arms of the Michelson optical fiber hydrophone, so that comprehensive distributed testing of strain distribution data of two interference optical fiber arms in an optical fiber coil is realized.
According to some embodiments, the scheme of the disclosure provides an automatic testing and analyzing method for strain distribution of an optical fiber hydrophone, which adopts the following technical scheme:
an automatic testing and analyzing method for strain distribution of an optical fiber hydrophone comprises the following steps:
step S01: acquiring the length of a long interference arm of an optical fiber hydrophone, the length of a short interference arm of the optical fiber hydrophone and the length of a transmission part optical fiber in each optical switch channel;
step S02: initializing a counting variable of an optical switch channel, and remotely setting parameters of an optical fiber strain distribution tester;
step S03: the optical fiber strain distribution tester performs strain distribution test of the optical fiber hydrophone to obtain strain distribution test data;
step S04: rapidly analyzing original Brillouin distribution data in strain distribution test data to respectively obtain the total optical fiber length, the interference arm reflector position and the long interference arm reflector position of the strain data;
Step S05: judging whether the total optical fiber length, the interference arm reflector position and the long interference arm reflector position of the obtained strain data meet the judging rule, and directly performing step S06 when the judging rule is met; otherwise, continuing to adjust parameters of the optical fiber strain distribution tester, and turning to step S03;
step S06: analyzing the original Brillouin distribution data and the original strain data to obtain long interference arm strain data and short interference arm strain data; and judging whether the current optical switch channel counting variable is not smaller than the optical switch channel number, if yes, directly outputting all the long interference arm strain data and the short interference arm strain data of the channels, and if not, adding 1 to the current optical switch channel counting variable, and turning to the step S02.
Compared with the prior art, the beneficial effects of the present disclosure are:
the method solves the problem that the traditional technology cannot test the strain distribution of the interference double arms of the Michelson optical fiber hydrophone, and realizes comprehensive distributed test of the strain distribution data of the two interference optical fiber arms in the optical fiber coil; the development and improvement process of the optical fiber coil winding process, the development and inspection process efficiency of the optical fiber hydrophone and the test accuracy are improved; further expands the application range of the optical fiber strain distribution tester and increases the application scene for BOTDR.
Drawings
The accompanying drawings, which are included to provide a further understanding of the disclosure, illustrate and explain the exemplary embodiments of the disclosure and together with the description serve to explain the disclosure, and do not constitute an undue limitation on the disclosure.
FIG. 1 is a schematic diagram of an automated test and analysis method for strain distribution in an optical fiber hydrophone in an embodiment of the disclosure;
FIG. 2 is a schematic diagram of an automatic test and analysis flow of strain distribution of a dual interference arm of an optical fiber hydrophone in an embodiment of the disclosure;
FIG. 3 is a flow chart of a fast fiber length analysis method based on fiber Brillouin scattering spectrum in an embodiment of the disclosure;
FIG. 4 is a schematic diagram of a temporary Brillouin spectral data calculation flow for removing a substrate according to an embodiment of the present disclosure;
FIG. 5 is a schematic diagram of an effective marker data analysis flow based on spectral maximum distribution data and spectral mean distribution data in an embodiment of the disclosure;
FIG. 6 is a schematic diagram of LSM and LLM analysis flow based on valid marker data in an embodiment of the present disclosure;
FIG. 7 is a schematic diagram of a flow chart of strain data analysis based on a long interference arm and a short interference arm in an embodiment of the disclosure;
FIG. 8 is a schematic diagram of an analysis flow of a redundant single reflection mirror brillouin scattering spectrum reconstruction algorithm according to an embodiment of the present disclosure;
FIG. 9 is a schematic diagram of an analysis flow of a reconstruction algorithm of a redundancy-free single reflection mirror image Brillouin scattering spectrum in an embodiment of the present disclosure;
FIG. 10 is a schematic diagram of an analysis flow (first part) of a multiple reflection mirror Brillouin scattering spectrum reconstruction algorithm according to an embodiment of the present disclosure;
fig. 11 is a schematic diagram of an analysis flow (second portion) of a multiple reflection mirror brillouin scattering spectrum reconstruction algorithm in an embodiment of the present disclosure.
The specific embodiment is as follows:
the disclosure is further described below with reference to the drawings and examples.
It should be noted that the following detailed description is illustrative and is intended to provide further explanation of the present disclosure. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this disclosure belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of exemplary embodiments in accordance with the present disclosure. As used herein, the singular is also intended to include the plural unless the context clearly indicates otherwise, and furthermore, it is to be understood that the terms "comprises" and/or "comprising" when used in this specification are taken to specify the presence of stated features, steps, operations, devices, components, and/or combinations thereof.
Embodiments of the present disclosure and features of embodiments may be combined with each other without conflict.
According to the method for automatically testing and analyzing the strain distribution of the optical fiber hydrophone, an optical fiber strain distribution tester (Brillouin optical time domain reflectometer, BOTDR for short) is adopted as testing equipment, a Michelson interference type optical fiber hydrophone optical fiber coil is tested, interference double-arm strain distribution data obtained through testing are automatically analyzed, and accordingly the respective optical fiber strain distribution data of two interference arms are obtained, and the typical embodiment is shown in fig. 1, and the specific composition is as follows:
(1) And the control computer is used for controlling the optical fiber strain distribution tester to test strain distribution and controlling the optical switch to switch the tested optical fiber coil.
(2) The optical fiber strain distribution tester is typically a brillouin optical time domain reflectometer based on a self-Brillouin scattering technology.
(3) The 1 XN optical switch is controlled by a control computer to switch the measured light path.
(4) The measured Michelson fiber optic hydrophone.
The method for automatically testing and analyzing the strain distribution of the optical fiber hydrophone shown in fig. 2 comprises the following specific steps:
step 101: starting a fiber coil strain distribution test, and turning to step 102;
Step 102: the control computer 1 reads the number N of the 1 XN optical switch channels and goes to step 103;
step 103: the user inputs the length LL 1-N of the long interference arm of the optical fiber hydrophone of each channel, the length LS 1-N of the short interference arm of each channel and the length L0-N of the optical fiber of the transmission part of each channel, and the step 104 is transferred;
step 104: initializing a current channel count variable osi=1 for the 1 xn optical switch, proceeding to step 105;
step 105: control computer 1 remotely controls the 1 xn optical switch to the OSI channel, proceeding to step 106
Step 106: controlling the computer 1 to remotely set parameters of the optical fiber strain distribution tester 2, wherein the spatial resolution SR is a settable minimum value, the measuring range value RAN is a minimum settable value exceeding LL+100deg.m, the sampling resolution DR is a settable minimum value, the starting frequency FS is a settable minimum value, the frequency interval FI is a settable minimum value, the ending frequency FE is a settable maximum value, and the step 107 is carried out;
step 107: control computer 1 remotely initiates a strain distribution test, go to step 108;
step 108: after the optical fiber strain distribution tester 2 completes the strain distribution test, the control computer 1 reads the strain data points SN, the frequency data points FM, the original Brillouin distribution data BORIDATA [ OSI ] [ 1-SN ] [ 1-FM ], the original strain data SORIDATA [ OSI ] [ 1-SN ], the distance data LDATA [ OSI ] [ 1-SN ], and then goes to step 109;
Step 109: analyzing BORIDATA [ OSI ] [ 1-SN ] [ 1-FM ] to obtain the total fiber length SL of the strain data, the position LSM of the interference arm reflector and the position LLM of the long interference arm reflector, turning to step 110 in detail, see step 201;
step 110: judging whether the SL, LSM and LLM obtained by analysis meet the basic data requirement, wherein LSTH is an analysis length difference judgment threshold, if the difference between SL and (LL [ OSI ] +L0) x 2 exceeds LSTH, namely SL- (LL [ OSI ] +L0) x 2 > LSTH or LSM=0 or LLM=0, judging that the SL, LSM and LLM obtained by analysis are invalid data, and replacing test parameters, and turning to step 111, otherwise, turning to step 112;
the LSTH typical value taking method comprises the following steps:
step 111: the control computer 1 increases the spatial resolution/pulse width SR set value by 1 gear step remotely, go to step 107;
step 112: analyzing the original Brillouin distribution data BORIDATA [ OSI ] [ 1-SN ] [ 1-FM ], the original strain data SORIDATA [ OSI ] [ 1-SN ] to obtain long interference arm strain data SSDATA [ OSI ] [ 1-SN ], and short interference arm strain data SSDATA [ OSI ] [ 1-SN ], see step 401 for details, and turning to step 113;
step 113: judging whether OSI is greater than or equal to N, if so, turning to step 115, otherwise turning to step 114;
Step 114: assign OSI to OSI +1, go to step 105;
step 115: all channel long interference arm strain data SSDATA [ 1-N ] [ 1-SN ] and short interference arm strain data SSDATA [ 1-N ] [ 1-SN ] are output.
As shown in fig. 3, the original brillouin spectrum data is rapidly analyzed, so as to obtain an optical fiber length SL, a short interference arm reflector position LSM, and a long interference arm reflector position LLM, which specifically include the following steps:
step 201: reading original Brillouin distribution data BORIDATA [ OSI ] [ 1-SN ] [ 1-FM ], distance data LDATA [ OSI ] [ 1-SN ], long interference arm length LL [ OSI ] of an optical fiber hydrophone, transmission part optical fiber length L0[ OSI ], spatial resolution SR, sampling resolution DR, starting frequency FS, frequency interval FI and ending frequency FE, and turning to step 202;
step 202: initializing counting variables II=1 and JJ=1, initializing temporary Brillouin data BDATA [ 1-SN ] [ 1-FM ], central frequency data CFDATA [ 1-SN ] to be all 0, effective marks BGSPU [ 1-SN ] to be all 0, brillouin spectrum intensity data BGSPDATA [ 1-FM ] to be all 0, brillouin spectrum frequency data BGSFDATA [ 1-FM ] to be all 0, noise intensity NOISEAVE=0, and turning to step 203;
step 203: calculating temporary Brillouin data BDATA [ 1-SN ] [ 1-FM ] according to the original Brillouin distribution data BORIDATA [ OSI ] [ 1-SN ] [ 1-FM ], see step 20301 for details, and turning to step 204;
Step 204: calculating effective marks BGSPU [ 1-SN ] of Brillouin spectra according to BDATA [ 1-SN ] [ 1-FM ], and turning to step 205, wherein the effective marks are shown in step 20401;
step 205: calculating a determination threshold slth=rounded up (SR/DR), i.e. the value of SR/DR is rounded up, if SLTH < 3, the value of SLTH is 3, go to step 206;
step 206: initializing a count variable ii=sn, initializing a valid brillouin spectrum count bgbouum=0, and turning to step 207;
step 207: judging whether the Brillouin spectrum effective mark BGSPU [ II ] is larger than the optical fiber end judging threshold value FETH, if so, turning to step 213, otherwise turning to step 208;
the typical value taking method of the optical fiber end judging threshold value FETH is as follows:
step 208: bgbouum=0, go to step 214;
step 209: judging that II is less than 2, if yes, turning to step 211, otherwise turning to step 210;
step 210: ii=ii-1, go to step 207;
step 211: fiber length subscript sln=0, go to step 212;
step 212: a short interference arm mirror position LSM value of 0 (refer to start position lsms=0, end position lsme=0) and a long interference arm mirror position LLM value of 0 (refer to start position llms=0, end position llme=0), go to step 217;
step 213: bgbouum=bgbouum+1, go to step 214;
Step 214: bgbouum > SLTH, go to step 215;
step 215: fiber length subscript sln=ii, go to step 216;
step 216: calculating a short interference arm reflector position LSM and a long interference arm reflector position LLM according to BGSPU [ 1-SN ], and turning to step 217;
step 217: fiber length sl=ldata [ OSI ] [ LSN ], go to step 218;
step 218: output fiber length SL, short interference arm mirror position LSM, and long interference arm mirror position LLM.
As shown in fig. 4, the noise-removed base calculation is performed on the original brillouin distributed data, so as to obtain temporary brillouin data BDATA, which includes the following specific steps:
step 20301: reading original brillouin distribution data boridinata [ OSI ] [ 1-SN ] [ 1-FM ], sampling resolution DR, initializing counting variables ii=1, jj=1, initializing temporary brillouin data BDATA [ 1-SN ] [ 1-FM ], noise intensity noise=0, and turning to step 20302;
step 20302: noise count noienum=100 meters/DR, step 20303;
step 20303: NOISEAVE = noiseav + borideta [ OSI ] [ SN-ii+1] [ JJ ], step 20304;
step 20304: judging whether JJ > FM is true, if so, turning to step 20306, otherwise turning to step 20305;
step 20305: jj=jj+1, go to step 20303;
Step 20306: jj=1, go to step 20307;
step 20307: judging whether II > NOISENUM is established, if so, turning to step 20309, and if not, turning to step 20308;
step 20308: ii=ii+1, go to step 20303;
step 20309: noiseal = noiseal/noiseal, step 203010;
step 20310: ii=1, jj=1, go to step 20311;
step 20311: BDATA [ II ] [ JJ ] = BORIDATA [ OSI ] [ II ] [ JJ ] -NOISEAVE, step 20312;
step 20312: judging whether JJ > FM is true, if true, turning to step 20314, and if false, turning to step 20313;
step 20313: jj=jj+1, go to step 20311;
step 20314: jj=1, go to step 20315;
step 20315: judging whether II > SN is true, if true, turning to step 20317, and if false, turning to step 20316;
step 20316: ii=ii+1, go to step 20311;
step 20317: and outputting the temporary Brillouin data BDATA [ 1-SN ] [ 1-FM ].
As shown in fig. 5, the temporary brillouin spectrum data is further analyzed to obtain effective tag data, which includes the following specific steps:
step 20401: reading temporary Brillouin data BDATA [ 1-SN ] [ 1-FM ] and initializing counting variables II=1 and JJ=1, wherein all initialized VMAX [ 1-SN ] are 0, all VAVE [ 1-SN ] are 0, and turning to step 20401;
Step 20402: VMAX [ II ] =bdata [ II ] [ JJ ], VAVE [ II ] =0, step 20403;
step 20403: judging whether VMAX [ II ] < BDATA [ II ] [ JJ ] is true, turning to step 20404 if true, turning to step 20405 if false;
step 20404: VMAX [ II ] =bdata [ II ] [ JJ ], go to step 20405;
step 20405: VAVE [ II ] =vave [ II ] +bdata [ II ] [ JJ ], go to step 20406;
step 20406: judging whether JJ > FM is true or not, if true, turning to step 20408, and if false, turning to step 20407;
step 20407: jj=jj+1, go to step 20403;
step 20408: VAVE [ II ] =vave [ II ]/FM, go to step 20409;
step 20409: judging whether II > SN is true, turning to step 20411, and turning to step 20410 is not true;
step 20410: ii=ii+1, go to step 20402;
step 20411: initializing temporary variable maxvmax=0, mavvave=0, ii=1, go to step 20412;
step 20412: initializing a temporary variable maxvmax=vmax [ II ], maxvave=vave [ II ], and turning to step 20413;
step 20413: judging whether MAXVMAX is less than VMAX [ II ] or not, if yes, turning to step 20414, otherwise turning to step 20415;
step 20414: maxvmax=vmax [ II ], step 20415;
step 20415: judging whether MAXVAVE < VAVE [ II ] is true, if so, turning to step 20416, otherwise, turning to step 20417;
Step 20416: maxvave=vave [ II ], step 20417;
step 20417: judging whether II > SN is true, if so, turning to step 20419, otherwise turning to step 20418;
step 20418: ii=ii+1, go to step 20413;
step 20419: initializing effective marks BGSPU [ 1-SN ] to be 0, II=1, and turning to step 20420;
step 20420:
BGSPU [ II ] =vmax [ II ]/MAXVMAX x VAVE [ II ]/MAXVAVE, step 20421;
step 20421: judging whether II > SN is true, if so, turning to step 20423, otherwise turning to step 20422;
step 20422: ii=ii+1, go to step 20420;
step 20423: and outputting BGSPU [ 1-SN ].
As shown in fig. 6, the effective mark data is further analyzed to obtain a short interference arm mirror position LSM and a long interference arm mirror position LLM, which specifically include the following steps:
step 301: reading BGSPU [ 1-SN ], distance data LDATA [ OSI ] [ 1-SN ], optical fiber length SL, initializing counting variables II=1 and PTEMP=1, initializing a short interference arm reflector position LSM, wherein the LSM comprises a starting position LSMS, a terminating position LSME and a long interference arm reflector position LLM, the LLM comprises a starting position LLMS and a terminating position LLME, and turning to step 302;
step 302: judging whether BGSPU [ II ] is more than or equal to 0.9 or not, if so, turning to step 305, otherwise turning to step 303;
Step 303: judging whether II > SN is true, turning to step 339, otherwise turning to step 304;
step 304: ii=ii+1, go to step 302;
step 305: initializing jj=0, lnum=0, and proceeding to step 306;
step 306: judging whether BGSPU [ II+JJ ] is not less than 0.9, if yes, turning to step 307, otherwise turning to step 310;
step 307: lnum=lnum+1, step 308;
step 308: judging whether II+JJ > SN is true, turning to step 310, otherwise turning to step 309;
step 309: jj=jj+1, go to step 306;
step 310: judging whether LNUM > SLTH is satisfied, if yes, turning to step 311, otherwise turning to step 308;
step 311: ptemp=ii+jj, go to step 312;
step 312: judging whether BGSPU [ II ] -BGSPU [ II-1] < 0 is established, establishing a step 315, otherwise, establishing a step 313; step 313: judging whether II > 2 is true, turning to step 314, otherwise turning to step 339;
step 314: ii=ii-1, go to step 312;
step 315: lsms=ii, go to step 315;
step 316: II = PTEMP, go to step 317;
step 317: judging whether BGSPU [ II ] -BGSPU [ II-1] < 0 is true, turning to step 320, otherwise turning to step 318; step 318: judging whether II < SN-1 is true, turning to step 319, and turning to step 339;
Step 319: ii=ii+1, go to step 317;
step 320: lsme=ii, go to step 321;
step 321: judging whether BGSPU [ II ] is more than or equal to 0.9, if so, turning to step 324, otherwise turning to step 322;
step 322: judging whether II > SN is true, turning to step 339, otherwise turning to step 323;
step 323: ii=ii+1, go to step 321;
step 324: initializing jj=0, lnum=0, and proceeding to step 325;
step 325: judging whether BGSPU [ II+JJ ] is not less than 0.9, if so, turning to step 326, otherwise turning to step 329;
step 326: lnum=lnum+1, go to step 327;
step 327: judging whether II+JJ > SN is true, turning to step 329, otherwise turning to step 328;
step 328: jj=jj+1, go to step 325;
step 329: judging whether LNUM > SLTH is satisfied, turning to step 330, otherwise turning to step 327;
step 330: ptemp=ii+jj, go to step 331;
step 331: judging whether BGSPU [ II ] -BGSPU [ II-1] < 0 is true, turning to step 334, otherwise turning to step 332; step 332: judging whether II > 2 is true, turning to step 333, otherwise turning to step 339;
step 333: ii=ii-1, go to step 331;
step 334: llms=ii, go to step 334;
step 335: II = PTEMP, go to step 336;
Step 336: judging whether BGSPU [ II ] -BGSPU [ II-1] < 0 is established, establishing a step 340, otherwise, establishing a step 337;
step 337: judging whether II < SN-1 is true, if so, turning to step 338, otherwise turning to step 339;
step 338: ii=ii+1, go to step 336;
step 339: calculating data errors, wherein the LSM value is 0 (refer to LSMS=0 and LSME=0), and the LLM value is 0 (refer to LLMS=0 and LLME=0);
step 340: lsme=ii, go to step 341;
step 341: outputting LSM and LLM.
As shown in fig. 7, the long interference arm strain data SSDATA [ OSI ] [1 to SN ], the short interference arm strain data SSDATA [ OSI ] [1 to SN ], are calculated as follows:
step 401: reading BORIDATA [ OSI ] [1 to SN ] [1 to FM ], SORIDATA [ OSI ] [1 to SN ], long interference arm length LL [ OSI ], short interference arm length LS [ OSI ], transmission section fiber length L0[ OSI ], fiber length SL, interference arm mirror position LSM and long interference arm mirror position LLM;
step 402: judging whether the LSMS is not more than 2 times LLMS, if yes, turning to step 403, otherwise turning to step 404;
step 403: redundant single reflection mirror image Brillouin scattering spectrum reconstruction is carried out on BORIDATA [ OSI ] [ 1-SN ] [ 1-FM ], the optical fiber length SL, the short interference arm reflector position LSM and the long interference arm reflector position LLM, short interference arm Brillouin spectrum data BSDATA [ 1-SN ] [ 1-FM ] and long interference arm Brillouin spectrum data BLDATA [ 1-SN ] [ 1-FM ] are constructed, and the step 501 is seen, and the step 407 is repeated;
Step 404: judging whether 3 x LSMS is less than or equal to 2 x LLMS is established, establishing a transition step 405, otherwise, switching to a step 406;
step 405: carrying out redundancy-free single reflection mirror image Brillouin scattering spectrum reconstruction on BORIDATA [ OSI ] [ 1-SN ] [ 1-FM ], the optical fiber length SL, the short interference arm reflector position LSM and the long interference arm reflector position LLM, and constructing short interference arm Brillouin spectrum data BSDATA [ 1-SN ] [ 1-FM ], and long interference arm Brillouin spectrum data BLDATA [ 1-SN ] [ 1-FM ], wherein the step 601 is seen in detail, and the step 407 is carried out;
step 406: carrying out multiple reflection mirror image Brillouin scattering spectrum reconstruction on BORIDATA [ OSI ] [ 1-SN ] [ 1-FM ], the optical fiber length SL, the short interference arm reflector position LSM and the long interference arm reflector position LLM to construct short interference arm Brillouin spectrum data BSDATA [ 1-SN ] [ 1-FM ], and long interference arm Brillouin spectrum data BLDATA [ 1-SN ] [ 1-FM ], wherein the steps are shown in step 701 and 407;
step 407: analyzing the brillouin spectrum data BSDATA [ 1-SN ] [ 1-FM ] of the short interference arm according to the brillouin scattering spectrum strain demodulation algorithm to obtain short interference arm strain data SSDATA [ OSI ] [ 1-SN ], and turning to step 408;
step 408: analyzing the brillouin spectrum data BLDATA [ 1-SN ] [ 1-FM ] of the long interference arm according to the brillouin scattering spectrum strain demodulation algorithm to obtain strain data SSDATA [ OSI ] [ 1-SN ] of the long interference arm, and turning to step 409;
Step 409: the long interference arm strain data SSDATA [ OSI ] [ 1-SN ] and the short interference arm strain data SSDATA [ OSI ] [ 1-SN ] are outputted.
As shown in fig. 8, for the case that 2×lsms is less than or equal to LLMS, the redundant single reflection mirror image brillouin scattering spectrum reconstruction algorithm is used to construct short interference arm brillouin spectrum data BSDATA [1 to SN ] [1 to FM ], and long interference arm brillouin spectrum data BLDATA [1 to SN ] [1 to FM ], which specifically comprises the following steps:
step 501: reading BORIDATA [ OSI ] [1 to SN ] [1 to FM ], long interference arm length LL [ OSI ], short interference arm length LS [ OSI ], transmission section fiber length L0[ OSI ], fiber length SL, interference arm mirror position LSM and long interference arm mirror position LLM, turning to step 502;
step 502: initializing all of short interference arm Brillouin spectrum data BSDATA [ 1-SN ] [ 1-FM ] to be 0, and all of long interference arm Brillouin spectrum data BLDATA [ 1-SN ] [ 1-FM ] to be 0, wherein II=1, turning to step 503;
step 503: jj=1, go to step 504;
step 504:
BSDATA[LSMS+1-II][JJ]＝BORIDATA[OSI][LSMS+1-II][JJ]-BORIDATA[OSI]
[ 2X LLMS-LSMS+II-1] [ JJ ], step 505.
Step 505: judging whether JJ > FM is true, if so, turning to step 507, otherwise turning to step 506;
step 506: jj=jj+1, go to step 504;
step 507: judging whether II > LSMS is established, if so, turning to step 509, otherwise turning to step 508;
Step 508: ii=ii+1, step 503;
step 509: ii=1, go to step 510;
step 510: jj=1, go to step 511;
step 511:
BLDATA [ lsme+1-II ] [ JJ ] =borideta [ OSI ] [2 x llms-lsme+ii-1] [ JJ ], go to step 512;
step 512: judging whether JJ > FM is true, if so, turning to step 514, otherwise turning to step 513;
step 513: jj=jj+1, go to step 511;
step 514: judging whether II > LSME is true, if so, turning to step 516, otherwise turning to step 515;
step 515: ii=ii+1, go to step 510;
step 516: ii=1, go to step 517;
step 517: jj=1, go to step 518;
step 518:
BLDATA[LSME+II][JJ]＝BORIDATA[OSI][LSME+II][JJ]-BSDATA[OSI]
[ 2. LSMS-LSME-II ] [ JJ ], step 519;
step 519: judging whether JJ > FM is true, if so, turning to step 521, otherwise turning to step 520;
step 520: jj=jj+1, go to step 518;
step 521: judging whether II is more than 2LSMS-LSME-1, if so, turning to step 523, otherwise turning to step 522;
step 522: ii=ii+1, go to step 517;
step 523: ii=1, go to step 524;
step 524: jj=1, go to step 525;
step 525: BLDATA [2 x lsme+ii ] [ JJ ] =borideta [ OSI ] [2 x lsme+ii ] [ JJ ], go to step 526;
step 526: determining whether JJ > FM is true, if so, turning to step 528, otherwise turning to step 527;
Step 527: jj=jj+1, go to step 525;
step 528: judging whether II is larger than LLMS-2LSMS, if so, turning to step 530, otherwise turning to step 529;
step 529: ii=ii+1, go to step 524;
step 530: short interference arm Brillouin spectrum data BSDATA [ 1-SN ] [ 1-FM ], long interference arm Brillouin spectrum data BLDATA [ 1-SN ] [ 1-FM ] are output.
As shown in fig. 9, for the case that 2/3×llms is greater than or equal to LSMS > 0.5×llms, the non-redundant single reflection mirror image brillouin scattering spectrum reconstruction algorithm is adopted to construct short interference arm brillouin spectrum data BSDATA [1 to SN ] [1 to FM ], and long interference arm brillouin spectrum data BLDATA [1 to SN ] [1 to FM ], which specifically comprises the following steps:
step 601: reading BORIDATA [ OSI ] [1 to SN ] [1 to FM ], long interference arm length LL [ OSI ], short interference arm length LS [ OSI ], transmission section fiber length L0[ OSI ], fiber length SL, interference arm mirror position LSM and long interference arm mirror position LLM, turning to step 602;
step 602: initializing all of short interference arm Brillouin spectrum data BSDATA [ 1-SN ] [ 1-FM ] to be 0, and all of long interference arm Brillouin spectrum data BLDATA [ 1-SN ] [ 1-FM ] to be 0, wherein II=1, turning to step 603;
step 603: jj=1, go to step 604;
step 604:
BSDATA[LSMS+1-II][JJ]＝BORIDATA[OSI][LSMS+1-II][JJ]-BORIDATA[OSI]
[2×LLMS-LSMS+II-1] [ JJ ], go to step 605;
Step 605: judging whether JJ > FM is true, turning to step 607, otherwise turning to step 606;
step 606: jj=jj+1, go to step 604;
step 607: judging whether II > LSMS is established, establishing a transition step 609, otherwise, transitioning to step 608;
step 608: ii=ii+1, go to step 603;
step 609: ii=1, go to step 610;
step 610: jj=1, go to step 611;
step 611:
BLDATA [ lsme+1-II ] [ JJ ] =borideta [ OSI ] [2 x llms-lsme+ii-1] [ JJ ], step 612;
step 612: judging whether JJ > FM is true, if so, turning to step 614, otherwise turning to step 613;
step 613: jj=jj+1, go to step 611;
step 614: judging whether II > LSME is true, if so, turning to step 616, otherwise turning to step 615;
step 615: ii=ii+1, go to step 610;
step 616: ii=1, go to step 617;
step 617: jj=1, go to step 618;
step 618:
BLDATA[LSME+II][JJ]＝BORIDATA[OSI][LSME+II][JJ]-BSDATA[OSI]
[ 2. LSMS-LSME-II ] [ JJ ], step 619;
step 619: judging whether JJ > FM is true, if so, turning to step 621, otherwise turning to step 620;
step 620: jj=jj+1, go to step 618;
step 621: judging whether II is more than LLMS-LSME-1, if so, turning to step 623, otherwise turning to step 622;
step 622: ii=ii+1, go to step 617;
Step 623: short interference arm Brillouin spectrum data BSDATA [ 1-SN ] [ 1-FM ], long interference arm Brillouin spectrum data BLDATA [ 1-SN ] [ 1-FM ] are output.
As shown in fig. 10 and 11, for the case of not less than LSMS2/3×llms, the short interference arm brillouin spectrum data BSDATA [1 to SN ] [1 to FM ], the long interference arm brillouin spectrum data BLDATA [1 to SN ] [1 to FM ] are constructed by adopting a multiple reflection mirror image brillouin scattering spectrum reconstruction algorithm, and the specific steps are as follows:
step 701: reading BORIDATA [ OSI ] [1 to SN ] [1 to FM ], long interference arm length LL [ OSI ], short interference arm length LS [ OSI ], transmission section fiber length L0[ OSI ], fiber length SL, interference arm mirror position LSM and long interference arm mirror position LLM, turning to step 702;
step 702: initializing temporary Brillouin data TBORIDATA [ 1-SN ] [ 1-FM ] = BORIDATA [ OSI ] [ 1-SN ] [ 1-FM ], short interference arm Brillouin spectrum data BSDATA [ 1-SN ] [ 1-FM ] are all 0, long interference arm Brillouin spectrum data BLDATA [ 1-SN ] [ 1-FM ] are all 0, turning to step 703;
step 703: calculating Δlm=llms-LSMS, initializing le=0, ls=2×Δlm, re=llms×2, rs=2×llms-2×Δlm, ii=1, turning to step 704;
step 704: judging whether LS is less than or equal to LSMS or not, if so, turning to step 705, otherwise, turning to step 719;
Step 705: jj=1, go to step 706;
step 706:
BSDATA [ LS+1-II ] [ JJ ] = TBORIDATA [ LS+1-II ] [ JJ ] -TBORIDATA [ RS+II-1] [ JJ ], BLDATA [ LS+1-II ] [ JJ ] = TBORIDATA [ RS+II-1] [ JJ ], step 707;
step 707: judging whether JJ > FM is true, if so, turning to step 709, otherwise turning to step 708;
step 708: jj=jj+1, go to step 706;
step 709: judging whether II is more than 2 DeltaLM is true, if so, turning to step 711, otherwise turning to step 710;
step 710: ii=ii+1, go to step 705;
step 711: re=re-2 Δlm, rs=rs-2 Δlm, ii=1, go to step 712;
step 712: jj=1, go to step 713;
step 713:
TBORIDATA [ RS-1+II ] [ JJ ] = TBORIDATA [ RS-1+II ] [ JJ ] -BSDATA [ LS-II+1] [ JJ ], go to step 714;
step 714: judging whether JJ > FM is true, if so, turning to step 716, otherwise turning to step 715;
step 715: jj=jj+1, go to step 713;
step 716: judging whether II is more than 2 DeltaLM is true, if so, turning to step 718, otherwise turning to step 717;
step 717: ii=ii+1, go to step 712;
step 718: le=le+2 Δlm, ls=ls+2 Δlm, ii=1, step 704;
step 719: ls=lsms, rs=llms+Δlm, step 720;
Step 720: jj=1, go to step 721;
step 721:
BSDATA [ LS+1-II ] [ JJ ] = TBORIDATA [ LS+1-II ] [ JJ ] -TBORIDATA [ RS+II-1] [ JJ ], BLDATA [ LS+1-II ] [ JJ ] = TBORIDATA [ RS+II-1] [ JJ ], go to step 722;
step 722: judging whether JJ > FM is true, if so, turning to step 724, otherwise turning to step 723;
step 723: jj=jj+1, go to step 721;
step 724: judging whether II > LS-LE is true, if so, turning to step 726, otherwise turning to step 725;
step 725: ii=ii+1, go to step 720;
step 726: re=re-2 Δlm, rs=llms, ii=1, go to step 727;
step 727: jj=1, go to step 728;
step 728:
BLDATA [ RS-1+II ] [ JJ ] = TBORIDATA [ RS-1+II ] [ JJ ] -BSDATA [ LS-II+1] [ JJ ], go to step 729;
step 729: judging whether JJ > FM is true, if so, turning to step 731, otherwise turning to step 730;
step 730: jj=jj+1, go to step 728;
step 731: judging whether II > LS-LE is true, if so, turning to step 733, otherwise turning to step 732;
step 732: ii=ii+1, go to step 727;
step 733: ls=lsme, le=lsms, re=lsms+ in the case of a ≡lm, rs=re+lsms-LSME, ii=1, step 734;
step 734: jj=1, go to step 735;
step 735: BLDATA [ LS+1-II ] [ JJ ] = TBORIDATA [ RS-1+II ] [ JJ ], step 736;
Step 736: judging whether JJ > FM is true, if so, turning to step 738, otherwise turning to step 737;
step 737: jj=jj+1, go to step 735;
step 738: judging whether II > LS-LE is true, if so, turning to step 740, otherwise turning to step 739;
step 739: ii=ii+1, step 734;
step 740: short interference arm Brillouin spectrum data BSDATA [ 1-SN ] [ 1-FM ], long interference arm Brillouin spectrum data BLDATA [ 1-SN ] [ 1-FM ] are output.
While the specific embodiments of the present disclosure have been described above with reference to the drawings, it should be understood that the present disclosure is not limited to the embodiments, and that various modifications and changes can be made by one skilled in the art without inventive effort on the basis of the technical solutions of the present disclosure while remaining within the scope of the present disclosure.
Claims (9)
1. The automatic testing and analyzing method for the strain distribution of the optical fiber hydrophone is characterized by comprising the following steps of:
step S01: acquiring the length of a long interference arm of an optical fiber hydrophone, the length of a short interference arm of the optical fiber hydrophone and the length of a transmission part optical fiber in each optical switch channel;
step S02: initializing a counting variable of an optical switch channel, and remotely setting parameters of an optical fiber strain distribution tester;
Step S03: the optical fiber strain distribution tester performs strain distribution test of the optical fiber hydrophone to obtain strain distribution test data;
step S04: rapidly analyzing original Brillouin distribution data in strain distribution test data to respectively obtain the total optical fiber length, the short interference arm reflector position and the long interference arm reflector position of the strain data;
step S05: judging whether the total optical fiber length, the short interference arm reflector position and the long interference arm reflector position of the obtained strain data meet a judging rule, and directly performing step S06 when the judging rule is met; otherwise, continuing to adjust parameters of the optical fiber strain distribution tester, and turning to step S03;
step S06: analyzing the original Brillouin distribution data and the original strain data to obtain long interference arm strain data and short interference arm strain data; judging whether the current optical switch channel counting variable is not smaller than the optical switch channel number, if yes, directly outputting all channel long interference arm strain data and short interference arm strain data, and if not, adding 1 to the current optical switch channel counting variable and turning to step S02;
in step S06, the specific process of obtaining the long interference arm strain data and the short interference arm strain data is as follows:
Step S601: taking original Brillouin distribution data BORIDATA [ OSI ] [ 1-SN ] [ 1-FM ], original strain data SORIDATA [ OSI ] [ 1-SN ], long interference arm length LL [ OSI ], short interference arm length LS [ OSI ], transmission part optical fiber length L0[ OSI ], optical fiber length SL, short interference arm reflector position LSM and long interference arm reflector position LLM; the number of the strain data points is SN, the number of the frequency data points is FM, and the counting variable of the optical switch channel is OSI;
step S602: judging whether the LSMS is not more than 2 times LLMS is established, if so, executing a step S603, otherwise, executing a step S604;
step S603: performing redundant single reflection mirror image Brillouin scattering spectrum reconstruction on BORIDATA [ OSI ] [ 1-SN ] [ 1-FM ], the optical fiber length SL, the short interference arm reflector position LSM and the long interference arm reflector position LLM, constructing short interference arm Brillouin spectrum data BSDATA [ 1-SN ] [ 1-FM ], and long interference arm Brillouin spectrum data BLDATA [ 1-SN ] [ 1-FM ], and executing step S607;
step S604: judging whether 3 x LSMS is less than or equal to 2 x LLMS is true, if so, executing step S605, otherwise, executing step S606;
step S605: carrying out redundancy-free single reflection mirror image Brillouin scattering spectrum reconstruction on BORIDATA [ OSI ] [ 1-SN ] [ 1-FM ], the optical fiber length SL, the short interference arm reflector position LSM and the long interference arm reflector position LLM, constructing short interference arm Brillouin spectrum data BSDATA [ 1-SN ] [ 1-FM ], and long interference arm Brillouin spectrum data BLDATA [ 1-SN ] [ 1-FM ], and turning to step S607, wherein the step is described in detail;
Step S606: carrying out multiple reflection mirror image Brillouin scattering spectrum reconstruction on BORIDATA [ OSI ] [ 1-SN ] [ 1-FM ], the optical fiber length SL, the short interference arm reflector position LSM and the long interference arm reflector position LLM, constructing short interference arm Brillouin spectrum data BSDATA [ 1-SN ] [ 1-FM ], and long interference arm Brillouin spectrum data BLDATA [ 1-SN ] [ 1-FM ], and executing step S607;
step S607: analyzing the brillouin spectrum data BSDATA [ 1-SN ] [ 1-FM ] of the short interference arm according to the brillouin scattering spectrum strain demodulation algorithm to obtain short interference arm strain data SSDATA [ OSI ] [ 1-SN ], and executing step S608;
step S608: analyzing the long interference arm Brillouin spectrum data BLDATA [ 1-SN ] [ 1-FM ] according to the Brillouin scattering spectrum strain demodulation algorithm to obtain long interference arm strain data SSDATA [ OSI ] [ 1-SN ], and executing step S609;
step S609: the long interference arm strain data SSDATA [ OSI ] [ 1-SN ] and the short interference arm strain data SSDATA [ OSI ] [ 1-SN ] are outputted.
2. The method for automatically testing and analyzing strain distribution of optical fiber hydrophone according to claim 1, wherein in step S02, in the process of remotely setting parameters of the optical fiber strain distribution tester, spatial resolution is set to a minimum value, a range value is set to a minimum value larger than a sum of a length of a long interference arm of the optical fiber hydrophone and 100m, sampling resolution is set to a minimum value, a start frequency is set to a minimum value, a frequency interval is set to a minimum value, and a stop frequency is set to a maximum value.
3. The method for automatically testing and analyzing strain distribution of optical fiber hydrophone as recited in claim 2, wherein in step S03, the strain distribution test data further includes strain data points, frequency data points, raw strain data, and distance data.
4. An automatic testing and analyzing method for strain distribution of optical fiber hydrophone as in claim 3, wherein in step S04, the total optical fiber length, the short interference arm reflector position and the long interference arm reflector position of the obtained strain data are as follows:
step S401: reading original Brillouin distribution data, distance data, long interference arm length of an optical fiber hydrophone, optical fiber length of a transmission part, spatial resolution, sampling resolution, starting frequency, frequency interval and ending frequency, setting initial values of counting variables to be 1, setting initial values of temporary Brillouin data, center frequency data, effective marks, brillouin spectrum intensity data, brillouin spectrum frequency data and noise intensity to be 0, setting effective Brillouin spectrum count to BGSPUNUM, and setting counting variables to be II;
step S402: calculating temporary Brillouin data according to the original Brillouin distribution data;
Step S403, calculating effective marks according to the temporary Brillouin data;
step S404: the calculation judgment threshold is rounded upwards, and if the calculation judgment threshold is smaller than 3, the calculation judgment threshold is taken as 3; setting the value of the counting variable as the number of the strain data points, and setting the initial value of BGSPUNUM as 0;
step S405: judging whether the Brillouin spectrum effective mark BGSPU [ II ] is larger than the optical fiber end judging threshold, if so, directly entering a step S409, otherwise, entering a step S406;
step S406, when BGSPUNUM is 0, turning to step S4010;
step S407, judging whether the counting variable is smaller than 2, if yes, directly marking the index of the optical fiber length as 0, otherwise, executing step S405 after the counting variable is reduced by 1;
step S408: assigning the LSM value of the short interference arm reflector to 0 and the LLM value of the long interference arm reflector to 0, and entering step S4011;
step S409: let bgspunum=bgspunum+1;
step S4010, setting BGSPUNUM to be larger than a calculation judgment threshold;
step S4011, when the index of the optical fiber length is equal to the counting variable, calculating the position of the short interference arm reflector and the position of the long interference arm reflector according to the effective mark, and outputting the optical fiber length, the position of the short interference arm reflector and the position of the long interference arm reflector.
5. The method for automatically testing and analyzing strain distribution of optical fiber hydrophone as recited in claim 4, wherein the specific process of step S402 is as follows:
step S40201: reading original brillouin distribution data and sampling resolution, setting an initial value of a counting variable to be 1, setting an initial value of noise intensity to be 0, initializing temporary brillouin data, and counting noise by NOISENUM=100 meters/DR; the counting variable of the optical switch channel is OSI, the number of the strain data points is SN, and the counting variables are II and JJ;
step S40202: calculating average noise
NOISEAVE = NOISEAVE + original brillouin profile borida [ OSI ] [ SN-ii+1] [ JJ ];
step S40203: judging whether the counting variable JJ is larger than the frequency data points, if so, outputting the counting variable JJ as 1, then executing the step S40204, otherwise, adding 1 to the counting variable JJ, and then executing the step S40202;
step S40204: judging whether the counting variable II is larger than noise counting, if so, outputting average noise NOISEAVE=NOISEAVE/NOISENUM, and then executing step S40205, otherwise, adding 1 to the counting variable II, and then turning to step S40202;
step S40205: design number variables II and JJ are both 1;
step S40206: temporary brillouin data BDATA [ II ] [ JJ ] =original brillouin data borida [ OSI ] [ II ] [ JJ ] -average noise;
Step S40207: judging whether the counting variable JJ is larger than the frequency data points, if so, outputting the counting variable JJ as 1, then executing the step S40208, otherwise, adding 1 to the counting variable JJ, and then executing the step S40206;
step S40208: and judging whether the counting variable II is larger than the number of the strain data points, if so, directly outputting temporary Brillouin data, otherwise, adding 1 to the counting variable II, and then turning to a step S40206.
6. The method for automatically testing and analyzing strain distribution of optical fiber hydrophone as recited in claim 4, wherein the specific process of step S403 is as follows:
step S40301: reading temporary Brillouin data, setting initial values of counting variables to be 1, setting initial values of maximum value data VMAX [ 1-SN ] and average value data VAVE [ 1-SN ] to be 0, setting frequency data points to be FM, setting strain data points to be SN, and setting counting variables to be II and JJ;
step S40302 sets VMAX [ II ] =temporary brillouin data BDATA [ II ] [ JJ ], VAVE [ II ] =0;
step S40303: judging whether the maximum value data VMAX [ II ] is larger than BDATA [ II ] [ JJ ], if yes, firstly setting VMAX [ II ] =BDATA [ II ] [ JJ ], then executing step S40304, otherwise, directly executing step S40304;
step S40304: setting average value data VAVE [ II ] =vave [ II ] +bdata [ II ] [ JJ ], judging whether the count variable JJ is greater than the frequency data point number, if yes, executing step S40305, otherwise, adding 1 to the count variable JJ, and then turning to step S40303;
Step S40305: setting VAVE [ II ] =VAVE [ II ]/FM, judging whether the counting variable II is larger than the number of the strain data points, if so, directly executing the step S40306, otherwise, adding 1 to the counting variable II, and then turning to the step S40302;
step S40306: initializing a temporary variable maxvmax=0, mavvave=0, ii=1, initializing a temporary variable maxvmax=vmax [ II ], maxvave=vave [ II ];
step S40307: judging whether the MAXVMAX is larger than VMAX [ II ], if so, setting MAXVMAX=VMAX [ II ] to execute the step S40308, otherwise, directly executing the step S40308;
step S40308: judging whether MAXVAVE is larger than VAVE [ II ], if yes, setting MAXVAVE=VAVE [ II ] and executing step S40309, otherwise, directly executing step S40309;
step S40309: judging whether the counting variable II is larger than the number of the strain data points, if yes, executing the step S40310, otherwise, adding 1 to the counting variable II, and then turning to the step S40307;
step S40310: initializing effective marks BGSPU [ 1-SN ] to be 0, wherein II=1;
step S40311: BGSPU [ II ] =vmax [ II ]/MAXVMAX x VAVE [ II ]/MAXVAVE;
step S40312: and judging whether the counting variable II is larger than the number of the strain data points, if so, directly outputting an effective mark, otherwise, adding 1 to the counting variable II, and then turning to a step S40311.
7. The method for automatically testing and analyzing strain distribution of optical fiber hydrophone as recited in claim 4, wherein in step S4011, the calculating the position LSM of the short interference arm reflector and the position LLM of the long interference arm reflector according to the effective marks comprises the following steps:
step S401101: reading effective marks BGSPU [ 1-SN ], distance data LDATA [ OSI ] [ 1-SN ] and optical fiber length SL of a Brillouin spectrum, initializing a short interference arm reflector position LSM, wherein the LSM comprises a start position LSMS, a stop position LSME and a long interference arm reflector position LLM, and the LLM comprises a start position LLMS and a stop position LLME; the number of the strain data points is SN, the judgment threshold is SLTH, the counting variables are II, JJ and PTEMP, and the judgment and counting of the threshold is LNUM;
step S401102: judging whether the Brillouin spectrum effective mark BGSPU [ II ] is more than or equal to 0.9 or not, if so, executing a step S401105, and if not, executing a step S401103;
step S401103: judging whether II > SN is true, if so, executing step S401139, otherwise, executing step S401104;
step S401104: ii=ii+1, step S401102 is performed;
step S401105: initializing jj=0, lnum=0;
Step S401106: judging whether the Brillouin spectrum effective mark BGSPU [ II+JJ ] is more than or equal to 0.9 is met, if so, executing a step S401107, otherwise, executing a step S401110;
step S401107: lnum=lnum+1;
step S401108: judging whether II+JJ > SN is true, if so, executing step S40110, otherwise, executing step S40109;
step S401109: jj=jj+1, and the process proceeds to step S401106;
step S401110: judging whether LNUM > SLTH is satisfied, if so, executing step S401111, otherwise, turning to step S401108;
step S401111: ptemp=ii+jj;
step S401112: judging whether BGSPU [ II ] -BGSPU [ II-1] < 0 is established, if yes, executing a step S401115, otherwise, executing a step S401113;
step S401113: judging whether II > 2 is true, if so, executing step S401114, otherwise, executing step S401139;
step S401114: ii=ii-1, and the process proceeds to step S401112;
step S401115: lsms=ii;
step S401116: II = PTEMP;
step S401117: judging whether BGSPU [ II ] -BGSPU [ II-1] < 0 is established, if yes, executing a step S401120, otherwise, executing a step S401118;
step S401118: judging whether II is smaller than SN-1, if so, executing step S401119, otherwise, executing step S401139;
Step S401119: ii=ii+1;
step S401120: lsme=ii;
step S401121: judging whether BGSPU [ II ] is more than or equal to 0.9 or not, if so, executing a step S401124, and if not, executing a step S401122;
step S401122: judging whether II > SN is true, if so, executing step S401139, otherwise, executing step S401123;
step S401123: ii=ii+1, and the process proceeds to step S401121;
step S401124: initializing jj=0, lnum=0;
step S401125: judging whether BGSPU [ II+JJ ] is more than or equal to 0.9, if so, executing a step S401126, and if not, executing a step S401129;
step S401126: lnum=lnum+1
Step S401127: judging whether II+JJ > SN is true, if so, executing step S401129, otherwise, executing step S401128;
step S401128: jj=jj+1;
step S401129: judging whether LNUM > SLTH is satisfied, if so, executing step S401130, otherwise, executing step S401127;
step S401130: ptemp=ii+jj;
step S401131: judging whether BGSPU [ II ] -BGSPU [ II-1] < 0 is established, if yes, executing a step S401134, otherwise, executing a step S401132;
step S401132: judging whether II > 2 is true, if so, executing step S401133, otherwise, executing step S401139;
Step S401133: ii=ii-1;
step S401134: llms=ii;
step S401135: II = PTEMP;
step S401136: judging whether BGSPU [ II ] -BGSPU [ II-1] < 0 is established, if yes, executing a step S401140, otherwise, executing a step S401137;
step S401137: judging whether II is smaller than SN-1, if so, executing step S401138, otherwise, executing step S401139;
step S401138: ii=ii+1;
step S401139: calculating data errors, wherein an LSM value of 0 means LSMS=0 and LSME=0; LLM value of 0 means llms=0, llme=0;
step S401140: lsme=ii;
step S401141: outputting LSM and LLM.
8. The method for automatically testing and analyzing strain distribution of optical fiber hydrophone according to claim 1, wherein the specific process of step S603 is as follows:
step S60301: reading original Brillouin distribution data BORIDATA [ OSI ] [ 1-SN ] [ 1-FM ], long interference arm length LL [ OSI ], short interference arm length LS [ OSI ], transmission part optical fiber length L0[ OSI ], optical fiber length SL, interference arm reflector position LSM and long interference arm reflector position LLM; the number of the strain data points is SN, the number of the frequency data points is FM, the counting variable of the optical switch channel is OSI, and the counting variable is II and JJ;
Step S60302: initializing all of short interference arm Brillouin spectrum data BSDATA [ 1-SN ] [ 1-FM ] to be 0, and all of long interference arm Brillouin spectrum data BLDATA [ 1-SN ] [ 1-FM ] to be 0, wherein II=1;
step S60303: jj=1;
step S60304:
BSDATA[LSMS+1-II][JJ]＝BORIDATA[OSI][LSMS+1-II][JJ]-BORIDATA[OSI]
[2*LLMS-LSMS+II-1][JJ]；
step S60305: judging whether JJ > FM is true, if so, executing step S60307, otherwise, executing step S60306;
step S60306: jj=jj+1;
step S60307: judging whether II > LSMS is established, if so, executing step S60309, otherwise, executing step S60308;
step S60308: ii=ii+1, and the process proceeds to step S60303;
step S60309: ii=1;
step S60310: jj=1;
step S60311:
BLDATA[LSME+1-II][JJ]＝BORIDATA[OSI][2*LLMS-LSME+II-1][JJ]；
step S60312: judging whether JJ > FM is true, if so, executing step S60314, otherwise, executing step S60313;
step S60313: jj=jj+1;
step S60314: judging whether II > LSME is true, if so, executing step S60316, otherwise, executing step S60315;
step S60315: ii=ii+1;
step S60316: ii=1;
step S60317: jj=1;
step S60318:
BLDATA[LSME+II][JJ]＝BORIDATA[OSI][LSME+II][JJ]-BSDATA[OSI]
[2*LSMS-LSME-II][JJ]；
step S60319: judging whether JJ > FM is true, if so, executing step S60321, otherwise, executing step S60320;
step S60320: jj=jj+1, and the process proceeds to step S60318;
Step S60321: judging whether II is more than 2LSMS-LSME-1 is true, if so, executing step S60323, otherwise, executing step S60322;
step S60322: ii=ii+1, and the process proceeds to step S60317;
step S60323: ii=1;
step S60324: jj=1;
step S60325: BLDATA [2 x lsme+ii ] [ JJ ] =borideta [ OSI ] [2 x lsme+ii ] [ JJ ];
step S60326: judging whether JJ > FM is true, if so, executing step S60328, otherwise, executing step S60327;
step S60327: jj=jj+1, and the process proceeds to step S60325;
step S60328: judging whether II is more than LLMS-2LSMS is established, if so, executing a step S60330, otherwise, executing a step S60329;
step S60329: ii=ii+1, and the process proceeds to step S60324;
step S60330: short interference arm Brillouin spectrum data BSDATA [ 1-SN ] [ 1-FM ], long interference arm Brillouin spectrum data BLDATA [ 1-SN ] [ 1-FM ] are output.
9. The method for automatically testing and analyzing strain distribution of optical fiber hydrophone as recited in claim 1, wherein the specific process of step S605 is as follows:
step S60501: reading original Brillouin distribution data BORIDATA [ OSI ] [ 1-SN ] [ 1-FM ], long interference arm length LL [ OSI ], short interference arm length LS [ OSI ], transmission part optical fiber length L0[ OSI ], optical fiber length SL, interference arm reflector position LSM and long interference arm reflector position LLM; the number of the strain data points is SN, the number of the frequency data points is FM, the counting variable of the optical switch channel is OSI, and the counting variable is II and JJ;
Step S60502: initializing all of short interference arm Brillouin spectrum data BSDATA [ 1-SN ] [ 1-FM ] to be 0, and all of long interference arm Brillouin spectrum data BLDATA [ 1-SN ] [ 1-FM ] to be 0, wherein II=1;
step S60503: jj=1;
step S60504:
BSDATA[LSMS+1-II][JJ]＝BORIDATA[OSI][LSMS+1-II][JJ]-BORIDATA[OSI]
[2*LLMS-LSMS+II-1][JJ]；
step S60505: judging whether JJ > FM is true, if so, executing step S60507, otherwise, executing step S60506:
step S60506: jj=jj+1, and the process proceeds to step S60504;
step S60507: judging whether II > LSMS is established, if so, executing step S60509, otherwise, executing step S60508;
step S60508: ii=ii+1, proceeding to step S60503;
step S60509: ii=1;
step S60510: jj=1;
step S60511:
BLDATA[LSME+1-II][JJ]＝BORIDATA[OSI][2*LLMS-LSME+II-1][JJ]；
step S60512: judging whether JJ > FM is true, if so, executing step S60514, otherwise, executing step S60513;
step S60513: jj=jj+1, and the process proceeds to step S60511;
step S60514: judging whether II > LSME is true, if so, executing step S60516, otherwise, executing step S60515;
step S60515: ii=ii+1, and the process proceeds to step S60510;
step S60516: ii=1;
step S60517: jj=1;
step S60518:
BLDATA[LSME+II][JJ]＝BORIDATA[OSI][LSME+II][JJ]-BSDATA[OSI]
[2*LSMS-LSME-II][JJ]；
step S60519: judging whether JJ > FM is true, if so, executing step S60521, otherwise, executing step S60520;
Step S60520: jj=jj+1, and the process proceeds to step S60518;
step S60521: judging whether II is more than LLMS-LSME-1, if so, executing a step S60523, otherwise, executing a step S60522;
step S60522: ii=ii+1, proceeding to step S60517;
step S60523: short interference arm Brillouin spectrum data BSDATA [ 1-SN ] [ 1-FM ], long interference arm Brillouin spectrum data BLDATA [ 1-SN ] [ 1-FM ] are output.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111015471.6A CN113639652B (en) | 2021-08-31 | 2021-08-31 | Automatic testing and analyzing method for strain distribution of optical fiber hydrophone |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111015471.6A CN113639652B (en) | 2021-08-31 | 2021-08-31 | Automatic testing and analyzing method for strain distribution of optical fiber hydrophone |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113639652A CN113639652A (en) | 2021-11-12 |
| CN113639652B true CN113639652B (en) | 2024-04-12 |
Family
ID=78424659
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111015471.6A Active CN113639652B (en) | 2021-08-31 | 2021-08-31 | Automatic testing and analyzing method for strain distribution of optical fiber hydrophone |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113639652B (en) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5268741A (en) * | 1992-01-31 | 1993-12-07 | Hewlett-Packard Company | Method and apparatus for calibrating a polarization independent optical coherence domain reflectometer |
| WO2010009951A1 (en) * | 2008-07-25 | 2010-01-28 | Thales | Self-referenced optical fibre sensor with stimulated brillouin scattering |
| CN103090813A (en) * | 2013-01-15 | 2013-05-08 | 电子科技大学 | High resolution sensing system measuring beat-length and strain of polarization maintaining optical fiber based on optical frequency domain reflectometer (OFDR) system |
| CN105333834A (en) * | 2015-12-11 | 2016-02-17 | 中国电子科技集团公司第四十一研究所 | System and method for automatic calibration of optical fiber strain coefficient |
| CN106482658A (en) * | 2016-11-17 | 2017-03-08 | 中国电子科技集团公司第四十研究所 | A kind of fibre strain coefficient automatic calibration method |
| CN107687939A (en) * | 2017-08-29 | 2018-02-13 | 鞍山睿科光电技术有限公司 | A kind of interference type optical fiber hydrophone pickup arm optical fiber detector and method |
| CN108168728A (en) * | 2017-12-11 | 2018-06-15 | 哈尔滨工程大学 | Non-equilibrium polarization maintaining optical fibre dual interferometer temperature strain simultaneous measuring apparatus and method |
| CN111103051A (en) * | 2019-12-30 | 2020-05-05 | 珠海任驰光电科技有限公司 | Optical fiber interference type hydrophone detection system and method |
| CN112525372A (en) * | 2020-11-10 | 2021-03-19 | 广东工业大学 | Device and method for simultaneously measuring strain temperature based on polarization maintaining optical fiber double-arm hetero-axial interferometer |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6728431B2 (en) * | 2001-10-15 | 2004-04-27 | The United States Of America As Represented By The Secretary Of The Navy | Fiber optic curvature sensor for towed hydrophone arrays |
-
2021
- 2021-08-31 CN CN202111015471.6A patent/CN113639652B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5268741A (en) * | 1992-01-31 | 1993-12-07 | Hewlett-Packard Company | Method and apparatus for calibrating a polarization independent optical coherence domain reflectometer |
| WO2010009951A1 (en) * | 2008-07-25 | 2010-01-28 | Thales | Self-referenced optical fibre sensor with stimulated brillouin scattering |
| CN103090813A (en) * | 2013-01-15 | 2013-05-08 | 电子科技大学 | High resolution sensing system measuring beat-length and strain of polarization maintaining optical fiber based on optical frequency domain reflectometer (OFDR) system |
| CN105333834A (en) * | 2015-12-11 | 2016-02-17 | 中国电子科技集团公司第四十一研究所 | System and method for automatic calibration of optical fiber strain coefficient |
| CN106482658A (en) * | 2016-11-17 | 2017-03-08 | 中国电子科技集团公司第四十研究所 | A kind of fibre strain coefficient automatic calibration method |
| CN107687939A (en) * | 2017-08-29 | 2018-02-13 | 鞍山睿科光电技术有限公司 | A kind of interference type optical fiber hydrophone pickup arm optical fiber detector and method |
| CN108168728A (en) * | 2017-12-11 | 2018-06-15 | 哈尔滨工程大学 | Non-equilibrium polarization maintaining optical fibre dual interferometer temperature strain simultaneous measuring apparatus and method |
| CN111103051A (en) * | 2019-12-30 | 2020-05-05 | 珠海任驰光电科技有限公司 | Optical fiber interference type hydrophone detection system and method |
| CN112525372A (en) * | 2020-11-10 | 2021-03-19 | 广东工业大学 | Device and method for simultaneously measuring strain temperature based on polarization maintaining optical fiber double-arm hetero-axial interferometer |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113639652A (en) | 2021-11-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108534811B (en) | Cavity length demodulation algorithm of short-cavity optical fiber Fabry-Perot sensor | |
| CN109167633A (en) | Optical-fiber laser sensing network multi-channel parallel demodulating system based on FPGA | |
| CN109282837B (en) | Demodulation method of Bragg fiber grating staggered spectrum based on LSTM network | |
| CN108896274B (en) | Distributed optical fiber strain demodulation method based on subset window length optimization algorithm | |
| CN101520509B (en) | Method for comparatively measuring range of interframe of chaotic light time domain waveform | |
| CN113639652B (en) | Automatic testing and analyzing method for strain distribution of optical fiber hydrophone | |
| CN111982000A (en) | Optical fiber shape reconstruction method and device based on Beta frame | |
| CN110424362B (en) | Optical fiber type temperature self-compensating static sounding sensor | |
| CN115597642A (en) | Optical frequency domain distributed optical fiber sensing method based on Kalman filtering | |
| CN103454034B (en) | Optical fiber micrometric displacement air pressure measuring apparatus | |
| CN111189556A (en) | Real-time multichannel fiber grating temperature measurement system based on AWG | |
| CN110940445A (en) | Optical fiber ring type residual stress test system and residual stress test method | |
| US20220390224A1 (en) | Special optical fiber for measuring 3d curved shape, manufacturing method thereof, and system for measuring 3d curved shape by using special optical fiber | |
| CN205027314U (en) | Novel experiment of fiber grating displacement sensing device | |
| CN210036956U (en) | Fiber grating double-edge correction demodulation vibration measuring device | |
| CN110071759B (en) | Optical cable fault positioning device and method based on polarized white light interference | |
| CN109844485A (en) | Method for estimating modal bandwidth spectral dependency | |
| CN116428977A (en) | Main cable strand positioning and adjusting method, system and storage medium based on laser array | |
| CN107144338B (en) | Equivalent noise pressure spectrum level calibration system of optical fiber acoustic sensor based on UPD low-noise photoelectric detector | |
| CN208833396U (en) | Silicon optical chip temperature measuring device based on OFDR | |
| CN108917636B (en) | A kind of distributive fiber optic strain demodulation method based on subset window adaption algorithm | |
| CN113405566B (en) | Optical fiber sensitive ring optical performance test method based on distributed polarization crosstalk | |
| CN111707303A (en) | Variable-step-length rapid high-precision signal demodulation method for optical fiber F-P sensor | |
| CN114279619B (en) | Anti-permeation high-sensitivity graphene hydraulic sensor and processing technology thereof | |
| CN113819931B (en) | Brillouin frequency shift extraction method for BOTDR and BOTDA fusion |
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 |