CN111220189A - Brillouin optical time domain analysis sensing device and non-local effect compensation method - Google Patents
Brillouin optical time domain analysis sensing device and non-local effect compensation method Download PDFInfo
- Publication number
- CN111220189A CN111220189A CN202010052968.4A CN202010052968A CN111220189A CN 111220189 A CN111220189 A CN 111220189A CN 202010052968 A CN202010052968 A CN 202010052968A CN 111220189 A CN111220189 A CN 111220189A
- Authority
- CN
- China
- Prior art keywords
- light
- pulse
- brillouin
- optical
- sensing
- 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
- 230000003287 optical effect Effects 0.000 title claims abstract description 56
- 230000000694 effects Effects 0.000 title claims abstract description 26
- 238000000034 method Methods 0.000 title claims abstract description 16
- 239000013307 optical fiber Substances 0.000 claims abstract description 29
- 238000001514 detection method Methods 0.000 claims abstract description 22
- 230000010287 polarization Effects 0.000 claims abstract description 10
- 239000000523 sample Substances 0.000 claims description 34
- 239000000835 fiber Substances 0.000 claims description 23
- 238000012546 transfer Methods 0.000 claims description 19
- 230000003993 interaction Effects 0.000 claims description 7
- 230000000630 rising effect Effects 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 238000001228 spectrum Methods 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D5/00—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
- G01D5/26—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light
- G01D5/32—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light
- G01D5/34—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells
- G01D5/353—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre
- G01D5/35338—Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable characterised by optical transfer means, i.e. using infrared, visible, or ultraviolet light with attenuation or whole or partial obturation of beams of light the beams of light being detected by photocells influencing the transmission properties of an optical fibre using other arrangements than interferometer arrangements
- G01D5/35354—Sensor working in reflection
- G01D5/35358—Sensor working in reflection using backscattering to detect the measured quantity
- G01D5/35364—Sensor working in reflection using backscattering to detect the measured quantity using inelastic backscattering to detect the measured quantity, e.g. using Brillouin or Raman backscattering
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Optical Transform (AREA)
Abstract
The invention discloses a Brillouin optical time domain analysis sensing device and a non-local effect compensation method.A continuous light output by a tunable laser is divided into an upper branch and a lower branch after passing through a coupler; continuous light in the upper branch enters the sensing optical fiber after sequentially passing through the polarization controller, the intensity modulator, the first erbium-doped optical fiber amplifier and the first optical circulator; continuous light in the lower branch passes through an acousto-optic modulator to realize optical pulse modulation, and then sequentially passes through a second erbium-doped optical fiber amplifier, a polarization scrambler and a second optical circulator and then enters a sensing optical fiber; the detection light enters a third optical circulator and a first photoelectric detector through a third port of the second optical circulator and then is sent to an upper computer through a data acquisition card; pulse light enters the adjustable attenuator through a third port of the first optical circulator, then passes through a second photoelectric detector and then is sent to an upper computer through a data acquisition card; the invention has simple structure, can effectively improve the performance of the sensing device and improve the compensation quality and precision.
Description
Technical Field
The invention relates to a distributed optical fiber sensing technology, in particular to a Brillouin optical time domain analysis sensing device and a non-local effect compensation method.
Background
In the last thirty years, with the rapid development of oil and gas pipelines, subways (high-speed railways), large buildings and the like, the safety of the pipelines is more and more emphasized by people. The distributed optical fiber sensing technology can use a common communication optical cable as a sensing medium, has the advantages of long sensing distance, more sensing nodes and the like, and is widely applied to infrastructure safety detection. The Brillouin optical time domain analysis technology is an important component of the distributed optical fiber sensing technology. The method is mainly applied to the fields of oil and gas pipelines, bridge structure health monitoring and the like at present. The brillouin optical time domain analysis technique obtains brillouin frequency shift (the parameter is in a linear relation with the temperature stress of the optical fiber) by measuring brillouin gain spectrum. The overall performance of the brillouin optical time domain analysis system is directly related to the signal-to-noise ratio of the brillouin gain. The signal-to-noise ratio is in turn directly proportional to the power of the continuous probe light input to the system. To mitigate non-local effects, the power of the continuous light in conventional single sideband gain configured systems is limited to-13 dBm. The current techniques for overcoming non-local effects are mainly classified into the following categories: 1. symmetric double-sideband detection light, which raises the detection light power to about-6 dBm; 2. the power of the detection light is raised to +5dBm by the symmetrical double-sideband detection light-scanning pulse light; 3. scanning the wavelength of a light source, and lifting the detection light to about-6.3 dBm; 4. modulating the frequency of the detection light, and increasing the power of the detection light to about-6 dBm; 5. the average Brillouin gain spectrum and the attenuation spectrum are used for increasing the detection light to 0 dBm; 6. the reflected rayleigh signal can be used to boost the probe light to about +2 dBm; the method 1 has a simple structure, but can only improve the detection optical power to-6 dBm; the methods 2 to 6 have complex structure configuration, and the signal-to-noise ratio of the compensation signal of the method 6 is low.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides a simple and high-performance Brillouin optical time domain analysis sensing device for compensating the non-local effect and a non-local effect compensation method.
The technical scheme adopted by the invention is as follows:
a Brillouin optical time domain analysis sensing device comprises an upper branch circuit and a lower branch circuit, wherein the upper branch circuit is used for generating continuous probe light, and the lower branch circuit is used for generating pulsed light; continuous light output by the tunable laser is divided into an upper branch and a lower branch after passing through the coupler; continuous light in the upper branch enters the sensing optical fiber after sequentially passing through the polarization controller, the intensity modulator, the first erbium-doped optical fiber amplifier and the first optical circulator; continuous light in the lower branch passes through an acousto-optic modulator to realize optical pulse modulation, and then sequentially passes through a second erbium-doped optical fiber amplifier, a polarization scrambler and a second optical circulator and then enters a sensing optical fiber; the intensity modulator is driven by a microwave source, and the acousto-optic modulator is driven by a pulse generator; after the stimulated Brillouin effect is carried out on the sensing light, the detection light enters a third optical circulator through a third port of the second optical circulator, then enters a first photoelectric detector and then is sent to an upper computer through a data acquisition card; after the stimulated Brillouin effect is performed on the sensing light, the pulse light enters the adjustable attenuator through the third port of the first optical circulator, then passes through the second photoelectric detector and then is sent to the upper computer through the data acquisition card.
Further, the third optical circulator is connected with an optical filter.
Further, the optical filter filtering window corresponds to an upper sideband of the double sidebands generated by the intensity modulator.
Further, the bandwidth of the first photoelectric detector is larger than the lowest sampling rate; the bandwidth of the second photodetector is greater than twice the frequency of the corresponding rising or falling edge of the pulse.
A non-local effect compensation method of a Brillouin optical time domain analysis sensing device comprises the following steps:
step 1: the probe light and the pulse light after undergoing the brillouin effect are expressed as formulas (1) and (3), and the probe light and the pulse light without the brillouin effect are expressed as formulas (2) and (4):
PSo(v,z)=Pise-α(L-z)e-αz(2)
Pp(v,z)=PiPe-αz(1-d(v,z)) (3)
Ppo(v,z)=PiPe-αz(4)
in the formula: pS(v, z) is the power of the probe light with energy transfer, PSo(v, z) is the power of the probe light without energy transfer; pp(v, z) is the power of the pulsed light with energy transfer, Ppo(v, z) is the power of the pulse light without energy transfer, v is the difference between the frequencies of the probe light and the pulse light, z is the position of the interaction between the probe light and the pulse light, Δ z is the length of the position of the interaction between the probe light and the pulse light, L is the total length of the sensing light, α is the light loss coefficient, P is the total length of the sensing light, andisfor detecting the power of light input into the sensing fibre, PiPFor pulsed light input into the sensing fiber, gB(v, z) is the linear gain of Brillouin, AeffA nonlinear effective area that is a light sensing mode; d (v, z) is the pulse consumption due to energy transfer from the 0 to z position due to the stimulated brillouin scattering effect;
step 2: probe light in z1The power gained by the location is as follows:
and step 3: assuming pulsed light in z1The position gives the probe light energy, which is carried to the z position by the pulsed light and also due to the propagation z-z1Distance and loss of energy; further assuming that the pulsed light does not impart probe light energy at the 0 to z positions, then these energies areEnergy E consumed by the pulsed(v,z):
Namely:
Ed(v,z)=(PPo(v,z)-PP(v,z))Δz=d(v,z)PiPΔze-αz(7)
and 4, step 4: the following formula can be obtained from formulas (6) and (7):
namely:
in the formula: pS(v,z)、PSo(v, z) and z are obtained according to the first photoelectric detector (14), and k is a consumption coefficient of pulsed light measured by the second photoelectric detector at the tail end of the optical fiber;
and 5: the power of the pulse at the fiber pigtail with and without energy transfer is as follows:
PP(v,z=L)=LssPiPe-αL(1-d(v,z=L)) (10)
PPo(v,z=L)=LssPiPe-αL(11)
in the formula: l isssIs a fixed loss factor; pP(v, z ═ L) is the power of the pulse at the fiber pigtail end with energy transfer, PPo(v, z ═ L) is the power of the pulse at the fiber pigtail end without energy transfer;
step 6: solving according to equations (10) and (11) yields d (v, z ═ L):
and 7: d (v, z) is obtained by substituting formula (12) for formula (9):
and 8: the original Brillouin gain can be obtained according to the formulao(v, z) and the compensated Brillouin gainc(v,z):
The invention has the beneficial effects that:
(1) the invention has simple structure and good compensation performance;
(2) the invention has high compensation signal quality, can effectively ensure the compensated precision, and has high precision of the measurement result;
(3) the invention has obviously improved performance, and can improve the detection light from-13 dBm to +2dBm under the condition of avoiding Brillouin error.
Drawings
FIG. 1 is a schematic view of the structure of the present invention.
Fig. 2 is a schematic structural diagram of the principle of the present invention.
Fig. 3 is a schematic diagram of a measurement result of pulsed light at the tail end of an optical fiber according to an embodiment of the present invention.
FIG. 4 is a d (v, z) distribution diagram in example.
FIG. 5 is a graphical comparison of Brillouin gain without compensation and compensated by the method of the present invention; a is a pulse consumption graph, b is an original Brillouin gain spectrum, and c is the compensated Brillouin gain spectrum by the method.
Fig. 6 is a comparison graph of the hot spot test results at 45 ℃, wherein a is an original brillouin frequency shift result graph, and b is a compensated brillouin frequency shift result graph.
In the figure: 1-tunable laser, 2-coupler, 3-polarization controller, 4-intensity modulator, 5-first erbium-doped fiber amplifier, 6-first optical circulator, 7-sensing fiber, 8-second optical circulator, 9-polarization scrambler, 10-acousto-optic modulator, 11-second erbium-doped fiber amplifier, 12-optical filter, 13-third optical circulator, 14-first photoelectric detector, 15-data acquisition card, 16-upper computer, 17-second photoelectric detector and 18-adjustable attenuator.
Detailed Description
The invention is further described with reference to the following figures and specific embodiments.
As shown in fig. 1, in a brillouin optical time domain analysis sensing device, continuous light output by a tunable laser 1 is divided into an upper path and a lower path through a coupler 2; continuous light in the upper branch sequentially passes through a polarization controller 3 to adjust the light polarization state, then is subjected to carrier suppression double-sideband modulation through an intensity modulator 4 driven by a microwave source, is improved in light power through a first erbium-doped fiber amplifier 5, and enters a sensing fiber 7 through a first optical circulator 6 to serve as detection light.
Continuous light in the lower branch is subjected to light pulse modulation through an acousto-optic modulator 10 driven by a pulse generator, then is amplified through a second erbium-doped fiber amplifier 11 in sequence, then is disturbed in the polarization state of the pulse light through a polarization scrambler 9, and is injected into a sensing fiber 7 through a second optical circulator 8 to serve as interaction of pulse light pre-detection light.
The detection light after the stimulated Brillouin effect in the sensing optical fiber 7 enters a third optical circulator 13 through a third port of a second optical circulator 8, the upper side band of the light is filtered by an optical filter 12, then the detection light enters a first photoelectric detector 14 for photoelectric conversion, and then the detection light is sent to an upper computer 16 through a data acquisition card 15 for data processing.
After the stimulated Brillouin effect in the sensing light 7, the pulse light enters the adjustable attenuator 18 through the third port of the first optical circulator 6 to adjust the light intensity, then is subjected to photoelectric conversion through the second photoelectric detector 17, and then is sent to the upper computer 20 through the data acquisition card 15 to be subjected to data processing.
The filter window of the optical filter 12 corresponds to the upper one of the double side bands generated by the intensity modulator 4. The bandwidth of the first photoelectric detector 14 is greater than the lowest sampling rate (the lowest sampling rate is 100MHz corresponding to the optical fiber spatial resolution 1m corresponding to the pulse light width 10 ns); the second photodetector 17 has a bandwidth greater than twice the frequency of the rising or falling edges of the pulses (10 ns for a rising edge corresponds to a frequency of 200 MHz).
The driving frequency of the intensity modulator 4 of the upper branch of the invention is the frequency sweep frequency fsThe frequency sweep range of the fiber needs to cover the Brillouin gain region of the sensing fiber; the acousto-optic modulator 10 of the lower branch not only generates pulsed light but also introduces optical frequency shift, and the optical frequency shift enables the lower sideband of the detection light to avoid the Brillouin loss region of the sensing optical fiber. And describing the distribution condition of pulse consumption in the optical fiber by adopting a continuous optical signal, and using the pulse at the tail end of the optical fiber to assist in calculating the pulse consumption coefficient.
A non-local effect compensation method of a Brillouin optical time domain analysis sensing device comprises the following steps:
step 1: the probe light and the pulse light after undergoing the brillouin effect are expressed as formulas (1) and (3), and the probe light and the pulse light without the brillouin effect are expressed as formulas (2) and (4):
PSo(v,z)=Pise-α(L-z)e-αz(2)
Pp(v,z)=PiPe-αz(1-d(v,z)) (3)
Ppo(v,z)=PiPe-αz(4)
in the formula: pS(v, z) is the power of the probe light with energy transfer, PSo(v, z) is the power of the probe light without energy transfer; pp(v, z) is the power of the pulsed light with energy transfer, Ppo(v, z) is the power of the pulse light without energy transfer, v is the difference between the frequencies of the probe light and the pulse light, z is the position of the interaction between the probe light and the pulse light, Δ z is the length of the position of the interaction between the probe light and the pulse light, L is the total length of the sensing light 7, α is the light loss coefficient, P is the total length of the sensing light 7, andisfor detecting the power, P, of light input into the sensing fibre 7iPFor pulsed light input into the power of the sensing fiber 7, gB(v, z) is the linear gain of Brillouin, AeffA nonlinear effective area that is a light sensing mode; d (v, z) isPulse consumption due to energy transfer from the 0 to z position by the stimulated brillouin scattering effect;
step 2: probe light in z1The power gained by the location is as follows:
and step 3: assuming pulsed light in z1The position gives the probe light energy, which is carried to the z position by the pulsed light and also due to the propagation z-z1Distance and loss of energy; further assuming that the pulse light does not impart probe light energy at the 0 to z positions, then these energies, i.e. the energy E consumed by the pulsed(v,z):
Namely:
Ed(v,z)=(PPo(v,z)-PP(v,z))Δz=d(v,z)PiPΔze-αz(7)
and 4, step 4: the following formula can be obtained from formulas (6) and (7):
namely:
in the formula: pS(v,z)、PSo(v, z) and z are obtained from the first photodetector 14, k is the consumption coefficient of the pulsed light at the end of the optical fiber measured by the second photodetector 17, z is ct/2, c is the propagation speed of the light in the optical fiber, the parameter to be determined in equation (9) is only k, and the loss coefficient of the optical fiber is α-1=22×103m。
And 5: the power of the pulse at the fiber pigtail with and without energy transfer is as follows:
PP(v,z=L)=LssPiPe-αL(1-d(v,z=L))(10)
PPo(v,z=L)=LssPiPe-αL(11)
in the formula: l isssA fixed loss factor introduced for the first optical circulator 6 and the adjustable attenuator 18; pP(v, z ═ L) is the power of the pulse at the fiber pigtail end with energy transfer, PPo(v, z ═ L) is the power of the pulse at the fiber pigtail end without energy transfer;
step 6: solving according to equations (10) and (11) yields d (v, z ═ L):
and 7: d (v, z) is obtained by substituting formula (12) for formula (9):
and 8: the original Brillouin gain can be obtained according to the formulao(v, z) and the compensated Brillouin gainc(v,z):
It can be seen from the above equation that the non-local effects of the pulse loss can be eliminated.
In practice, the frequency range needs to cover the fiber brillouin shift by the microwave source sweeping the frequency. FIGS. 3 and 4 are graphs of the results of tests using the apparatus of the present invention, using 39.1km of fiber. Fig. 3 is a graph showing the measurement result of the pulse light at the end of the optical fiber, and fig. 4 is a d (v, z) distribution graph. Pulse consumption plots at the proximal end of 6km and at the end of the fibre of 39.105km are also shown. It can be seen that the d (v, z) signal-to-noise ratio is high, thus ensuring the post-compensation measurement accuracy.
FIG. 5 is a graph of Brillouin gain versus Brillouin gain without compensation in accordance with the present invention; under the condition that the detection light is +2 dBm; a is the pulse consumption graph, b is the original brillouin gain spectrum, and c is the compensated brillouin gain spectrum. Comparing the graphs b and c can show that the line width of the gain spectrum before compensation is inconsistent with the Brillouin frequency shift, and the line width after compensation is highly consistent with the Brillouin frequency shift.
Fig. 6 is a comparison graph of the hot spot test results at 45 ℃, wherein a is an original brillouin frequency shift result graph, and b is a compensated brillouin frequency shift result graph. It can be seen from the graph that the original brillouin frequency shift at the hot spot is not consistent, and the compensated brillouin frequency shift is basically consistent.
The invention uses continuous detecting light to draw the distribution condition of pulse consumption in the optical fiber, and uses the pulse consumption value to compensate the non-local effect, thereby improving the performance and accuracy of the whole system. The system has a simple structure, only needs the configuration of the most basic Brillouin analysis system, and effectively improves the detection light from-13 dBm to +2dBm under the condition of avoiding the Brillouin error, thereby improving the performance of the sensing system; and the quality compensation signal can be improved, thereby ensuring the compensated precision.
Claims (5)
1. A Brillouin optical time domain analysis sensing device is characterized by comprising an upper branch circuit and a lower branch circuit, wherein the upper branch circuit is used for generating continuous probe light, and the lower branch circuit is used for generating pulsed light; continuous light output by the tunable laser (1) is divided into an upper branch and a lower branch after passing through the coupler (2); continuous light in the upper branch enters a sensing optical fiber (7) after sequentially passing through a polarization controller (3), an intensity modulator (4), a first erbium-doped optical fiber amplifier (5) and a first optical circulator (6); continuous light in the lower branch passes through an acousto-optic modulator (10) to realize optical pulse modulation, then sequentially passes through a second erbium-doped optical fiber amplifier (11), a polarization scrambler (9) and a second optical circulator (8) and then enters a sensing optical fiber (7); the intensity modulator (4) is driven by a microwave source, and the acousto-optic modulator (10) is driven by a pulse generator; the detection light enters a third optical circulator (13) through a third port of a second optical circulator (8) after being subjected to stimulated Brillouin action in the sensing light (7), then enters a first photoelectric detector (14), and then is sent to an upper computer (16) through a data acquisition card (15); after being subjected to the stimulated Brillouin action in the sensing light (7), the pulse light enters the adjustable attenuator (18) through the third port of the first optical circulator (6), then passes through the second photoelectric detector (17) and then is sent to the upper computer (20) through the data acquisition card (15).
2. A brillouin optical time domain analysis sensing apparatus according to claim 1, wherein said third optical circulator (13) is connected to an optical filter (12).
3. A brillouin optical time domain analysis sensing apparatus according to claim 2, wherein the filter window of the optical filter (12) corresponds to an upper one of the double sidebands produced by the intensity modulator (4).
4. A brillouin optical time domain analysis sensing apparatus according to claim 1, wherein said first photodetector (14) has a bandwidth greater than a lowest sampling rate; the bandwidth of the second photodetector (17) is greater than twice the corresponding frequency of the rising or falling edge of the pulse.
5. The non-local effect compensation method adopting the Brillouin optical time domain analysis sensing device as claimed in any one of claims 1 to 4, characterized by comprising the following steps:
step 1: the probe light and the pulse light after undergoing the brillouin effect are expressed as formulas (1) and (3), and the probe light and the pulse light without the brillouin effect are expressed as formulas (2) and (4):
PSo(v,z)=Pise-α(L-z)e-αz(2)
Pp(v,z)=PiPe-αz(1-d(v,z)) (3)
Ppo(v,z)=PiPe-αz(4)
in the formula: pS(v, z) is probe light with energy transferPower of PSo(v, z) is the power of the probe light without energy transfer; pp(v, z) is the power of the pulsed light with energy transfer, Ppo(v, z) is the power of the pulse light without energy transfer, v is the difference between the frequencies of the probe light and the pulse light, z is the position of the interaction between the probe light and the pulse light, Δ z is the length of the position of the interaction between the probe light and the pulse light, L is the total length of the sensing light (7), α is the light loss coefficient, P is the total length of the sensing light (7)isFor detecting the power, P, of light input into the sensing fibre (7)iPFor the pulsed light input into the power of the sensing fiber (7), gB(v, z) is the linear gain of Brillouin, AeffA nonlinear effective area that is a light sensing mode; d (v, z) is the pulse consumption due to energy transfer from the 0 to z position due to the stimulated brillouin scattering effect;
step 2: probe light in z1The power gained by the location is as follows:
and step 3: assuming pulsed light in z1The position gives the probe light energy, which is carried to the z position by the pulsed light and also due to the propagation z-z1Distance and loss of energy; further assuming that the pulse light does not impart probe light energy at the 0 to z positions, then these energies, i.e. the energy E consumed by the pulsed(v,z):
Namely:
Ed(v,z)=(PPo(v,z)-PP(v,z))Δz=d(v,z)PiPΔze-αz(7)
and 4, step 4: the following formula can be obtained from formulas (6) and (7):
namely:
in the formula: pS(v,z)、PSo(v, z) and z are obtained according to the first photoelectric detector (14), and k is a consumption coefficient of pulsed light measured by the second photoelectric detector (17) at the tail end of the optical fiber;
and 5: the power of the pulse at the fiber pigtail with and without energy transfer is as follows:
PP(v,z=L)=LssPiPe-αL(1-d(v,z=L)) (10)
PPo(v,z=L)=LssPiPe-αL(11)
in the formula: l isssIs a fixed loss factor; pP(v, z ═ L) is the power of the pulse at the fiber pigtail end with energy transfer, PPo(v, z ═ L) is the power of the pulse at the fiber pigtail end without energy transfer;
step 6: solving according to equations (10) and (11) yields d (v, z ═ L):
and 7: d (v, z) is obtained by substituting formula (12) for formula (9):
and 8: the original Brillouin gain can be obtained according to the formulao(v, z) and the compensated Brillouin gainc(v,z):
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010052968.4A CN111220189B (en) | 2020-01-17 | 2020-01-17 | Brillouin optical time domain analysis sensing device and non-local effect compensation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010052968.4A CN111220189B (en) | 2020-01-17 | 2020-01-17 | Brillouin optical time domain analysis sensing device and non-local effect compensation method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111220189A true CN111220189A (en) | 2020-06-02 |
CN111220189B CN111220189B (en) | 2021-04-13 |
Family
ID=70829659
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010052968.4A Active CN111220189B (en) | 2020-01-17 | 2020-01-17 | Brillouin optical time domain analysis sensing device and non-local effect compensation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111220189B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113595638A (en) * | 2021-07-21 | 2021-11-02 | 华南师范大学 | BOTDA system based on four-frequency-division driving |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050213869A1 (en) * | 2004-03-26 | 2005-09-29 | Anthony Brown | System and method for resolution enhancement of a distributed sensor |
CN102636217A (en) * | 2012-01-09 | 2012-08-15 | 宋牟平 | Sensing device based on joint detection of Brillouin optical time domain analysis and Mach-Zehnder interference |
CN102998025A (en) * | 2012-12-18 | 2013-03-27 | 华北电力大学(保定) | Measuring method for pulse pre-pump rayleigh BOTDA (Brilouin optical time domain analysis) temperature and strain |
CN103364106A (en) * | 2013-07-16 | 2013-10-23 | 扬州森斯光电科技有限公司 | Brillouin optical time domain analyzing and sensing system based on ultra-long annular laser pumping |
CN103913185A (en) * | 2014-03-31 | 2014-07-09 | 广西师范大学 | Brillouin optical fiber sensing system and method |
CN104567960A (en) * | 2015-01-04 | 2015-04-29 | 西南交通大学 | Coherent Brillouin optical time-domain analysis sensing system based on phase modulation probe light |
CN105675031A (en) * | 2016-01-23 | 2016-06-15 | 中国人民解放军国防科学技术大学 | Pre-pumped pulse and Gray code based BOTDA (Brillouin Optical Time Domain Analysis) instrument |
CN107091698A (en) * | 2017-06-16 | 2017-08-25 | 苏州光格设备有限公司 | Brillouin optical time domain analysis system and method |
CN109974760A (en) * | 2019-01-24 | 2019-07-05 | 西南交通大学 | A kind of Brillouin optical time domain analysis system based on Brillouin's phase shift demodulation |
-
2020
- 2020-01-17 CN CN202010052968.4A patent/CN111220189B/en active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050213869A1 (en) * | 2004-03-26 | 2005-09-29 | Anthony Brown | System and method for resolution enhancement of a distributed sensor |
CN102636217A (en) * | 2012-01-09 | 2012-08-15 | 宋牟平 | Sensing device based on joint detection of Brillouin optical time domain analysis and Mach-Zehnder interference |
CN102998025A (en) * | 2012-12-18 | 2013-03-27 | 华北电力大学(保定) | Measuring method for pulse pre-pump rayleigh BOTDA (Brilouin optical time domain analysis) temperature and strain |
CN103364106A (en) * | 2013-07-16 | 2013-10-23 | 扬州森斯光电科技有限公司 | Brillouin optical time domain analyzing and sensing system based on ultra-long annular laser pumping |
CN103913185A (en) * | 2014-03-31 | 2014-07-09 | 广西师范大学 | Brillouin optical fiber sensing system and method |
CN104567960A (en) * | 2015-01-04 | 2015-04-29 | 西南交通大学 | Coherent Brillouin optical time-domain analysis sensing system based on phase modulation probe light |
CN105675031A (en) * | 2016-01-23 | 2016-06-15 | 中国人民解放军国防科学技术大学 | Pre-pumped pulse and Gray code based BOTDA (Brillouin Optical Time Domain Analysis) instrument |
CN107091698A (en) * | 2017-06-16 | 2017-08-25 | 苏州光格设备有限公司 | Brillouin optical time domain analysis system and method |
CN109974760A (en) * | 2019-01-24 | 2019-07-05 | 西南交通大学 | A kind of Brillouin optical time domain analysis system based on Brillouin's phase shift demodulation |
Non-Patent Citations (1)
Title |
---|
梁家伟: "基于光梳放大的布里渊光时域分析***优化研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113595638A (en) * | 2021-07-21 | 2021-11-02 | 华南师范大学 | BOTDA system based on four-frequency-division driving |
CN113595638B (en) * | 2021-07-21 | 2022-05-17 | 华南师范大学 | BOTDA system based on four-frequency-division driving |
Also Published As
Publication number | Publication date |
---|---|
CN111220189B (en) | 2021-04-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111307054B (en) | High-precision dynamic strain monitoring device and method based on time-delay-free chaotic laser | |
CN110440851B (en) | Long-distance multi-parameter measuring device and method based on Brillouin and Raman scattering | |
CN110375800B (en) | Sensing device and method based on super-continuum spectrum Brillouin optical time domain analyzer | |
CN101762290A (en) | Distributed Raman amplification-based Brillouin optical time domain analysis system | |
CN111721338B (en) | Brillouin optical time domain analysis system for alternately modulating high frequency and low frequency of pump light | |
CN103115695B (en) | Double-sideband distributed type optical fiber sensing system parameter measuring device | |
CN103913185A (en) | Brillouin optical fiber sensing system and method | |
CN114088124B (en) | Brillouin optical time domain reflectometer based on double-sideband modulation | |
CN113049014B (en) | Time-frequency multiplexing BOTDA system based on pumping light frequency sweep and sensing method | |
CN112697181A (en) | Phase sensitive optical time domain reflection device and method based on frequency modulation | |
CN108801305B (en) | Method and device of Brillouin optical time domain reflectometer based on step pulse self-amplification | |
CN109163748A (en) | Single-ended dynamic distributed Brillouin's reflection unit and method based on frequency agility technology | |
CN111854918B (en) | Subband processing method for improving spatial resolution of distributed optical fiber detection | |
CN104483289A (en) | Birefringence detection device and birefringence detection method based on sweep frequency optical coherence tomography technology | |
CN111707299A (en) | BOTDA device for obtaining gain loss spectrum based on pumping light frequency sweep | |
CN111220189B (en) | Brillouin optical time domain analysis sensing device and non-local effect compensation method | |
CN203758532U (en) | Brillouin fiber-optic sensing system | |
CN103376124A (en) | Brillouin optical time domain analyzer | |
CN111141414B (en) | Temperature and strain simultaneous measurement device and method based on chaos BOCDA | |
CN110375960B (en) | Device and method based on super-continuum spectrum light source OTDR | |
CN116105079A (en) | Distributed petroleum pipeline monitoring method and system based on forward stimulated Brillouin scattering | |
CN111637910A (en) | Time domain differential high-speed chaotic Brillouin optical coherent domain monitoring device and method | |
CN114414993B (en) | Method for testing frequency response of chip of electro-optical intensity modulator | |
CN109959627B (en) | Rapid gas absorption spectrum measuring device and method based on optical frequency agility | |
Zheng et al. | Single shot OCC-BOTDA based on polarization diversity and image denoising |
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 |