WO2022259437A1 - 振動測定器及び振動測定方法 - Google Patents
振動測定器及び振動測定方法 Download PDFInfo
- Publication number
- WO2022259437A1 WO2022259437A1 PCT/JP2021/021991 JP2021021991W WO2022259437A1 WO 2022259437 A1 WO2022259437 A1 WO 2022259437A1 JP 2021021991 W JP2021021991 W JP 2021021991W WO 2022259437 A1 WO2022259437 A1 WO 2022259437A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- optical
- vibration
- phase
- frequency
- pulse
- Prior art date
Links
- 238000000691 measurement method Methods 0.000 title claims description 6
- 238000005259 measurement Methods 0.000 title abstract description 21
- 230000003287 optical effect Effects 0.000 claims abstract description 146
- 239000013307 optical fiber Substances 0.000 claims abstract description 50
- 230000003595 spectral effect Effects 0.000 claims description 5
- 238000001514 detection method Methods 0.000 abstract description 14
- 238000001228 spectrum Methods 0.000 abstract 1
- 239000013598 vector Substances 0.000 description 34
- 238000000034 method Methods 0.000 description 23
- 238000012545 processing Methods 0.000 description 23
- 238000005070 sampling Methods 0.000 description 14
- 239000000835 fiber Substances 0.000 description 11
- 238000012360 testing method Methods 0.000 description 10
- 238000000253 optical time-domain reflectometry Methods 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000005562 fading Methods 0.000 description 5
- 230000001427 coherent effect Effects 0.000 description 4
- 239000000523 sample Substances 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 238000012937 correction Methods 0.000 description 3
- 230000003111 delayed effect Effects 0.000 description 3
- 230000005684 electric field Effects 0.000 description 3
- 230000002123 temporal effect Effects 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000001174 ascending effect Effects 0.000 description 1
- 230000008033 biological extinction Effects 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 1
- 230000002040 relaxant effect Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
- G01H9/004—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means using fibre optic sensors
Definitions
- the present disclosure relates to a vibration measuring instrument and a vibration measuring method for performing DAS-P.
- Non-Patent Document 1 DAS (Distributed Acoustic Sensing) that detects backscattered light due to Rayleigh scattering by injecting pulsed test light into the optical fiber under test as a means of measuring the physical vibration applied to the optical fiber in a distributed manner in the longitudinal direction of the optical fiber ) is known (Non-Patent Document 1).
- the DAS captures the change in the optical path length of the optical fiber due to the physical vibration applied to the optical fiber, and senses the vibration. By detecting the vibration, it is possible to detect the movement of an object around the optical fiber to be measured.
- DAS-I DAS-intensity
- DAS-P DAS-phase
- DAS-phase DAS-phase
- the phase changes linearly with respect to changes in the optical path length of the fiber due to vibration, and the rate of change is the same in the longitudinal direction of the optical fiber.
- the vibration can be measured quantitatively, and the vibration applied to the optical fiber to be measured can be faithfully reproduced (for example, Non-Patent Document 2).
- the phase of the scattered light at the time t when the pulsed light is injected is distributed in the longitudinal direction of the optical fiber. That is, the phase ⁇ (l, t) of the scattered light is measured as the distance l from the incident end of the optical fiber.
- the time at which the point at distance l is measured is delayed from the time at which the pulse is incident by the time it takes for the pulsed light to propagate from the incident end to the distance l. Furthermore, it should be noted that the measuring time is delayed by the time required for the scattered light to return to the incident end.
- the magnitude of the physical vibration applied to the section from the distance l to the distance l + ⁇ l at each time nT is the phase ⁇ (l + ⁇ l, nT) at the distance l + ⁇ l, the phase ⁇ (l, nT) at the distance l, and is proportional to the difference .delta..theta.(l, nT). That is, if time zero is used as a reference, the following expression is satisfied.
- the device configuration for detecting the phase of the scattered light includes a direct detection configuration in which the backscattered light from the optical fiber under test is directly detected by a photodiode, etc., and a coherent detection configuration in which the light is combined with a separately prepared reference light There is a configuration using detection (for example, Non-Patent Document 1).
- the mechanism that performs coherent detection and calculates the phase is subdivided into two: a software-based processing mechanism using the Hilbert transform and a hardware-based processing mechanism using a 90-degree optical hybrid. Also in this method, the in-phase component I(l, nT) and the quadrature component Q(l, nT) of the scattered light are obtained, and the phase is calculated by the following equation.
- the output value by the four-quadrant arctangent operator Arctan is in the range of (- ⁇ , ⁇ ] in radian units, and m is an arbitrary integer, and 2m ⁇ + ⁇ (l, nT) are all the same vector direction on the xy plane.
- phase unwrapping if the unwrapped phase is ⁇ cal unwrap , for example, when processing is performed in ascending order of time, ⁇ cal unwrap is the same as ⁇ cal at the start point of phase unwrapping.
- ⁇ cal unwrap if any integer p, calculate ⁇ cal unwrap (l, (p+1)T) from ⁇ cal unwrap ( l , pT) as follows.
- phase unwrapping processing is often performed on the calculated difference.
- the direction of the vector when there is no noise corresponds to the phase to be measured.
- the influence of noise is large, there is a probability that the direction of the vector will be in the opposite direction, and the actual measured phase value will differ by ⁇ radians from the ideal phase value in the absence of noise. growing.
- the magnitude of vibration from equation (1) it leads to an erroneous recognition that a large physical force is applied to the optical fiber.
- the influence of the noise of the measuring instrument increases is that the intensity of the scattered light itself decreases when the noise of the measuring instrument can be considered to be the same at each point and time. If it can be increased at points and at each time, it will be possible to reduce the influence of noise in the measuring instrument.
- the reason why the scattered light intensity itself decreases is not only the loss due to absorption and scattering that occurs as the probe pulse light propagates through the optical fiber to be measured. Since pulsed light with a finite time width is incident on the optical fiber under test and the scattering of the pulsed light is detected, scattering from a large number of very finely distributed scatterers on the optical fiber Light interference occurs. As a result of the interference, there occurs a point where the intensity of the scattered light becomes small according to the distribution of the scatterers in the longitudinal direction of the optical fiber under measurement at each time. This phenomenon is called fading (Non-Patent Document 3).
- Patent Document 1 when measuring the phase of scattered light in DAS-P, a phase measurement method and a signal processing device that can reduce the influence of noise in the measuring instrument without increasing the peak intensity of the incident light pulse has been proposed.
- pulsed light which is wavelength-multiplexed by arranging pulses of different optical frequency components at time intervals at which changes in the state of the fiber due to vibration can be ignored, is incident on the optical fiber to be measured.
- the scattered light vector obtained by plotting the scattered light at each wavelength from the optical fiber on a two-dimensional plane with the in-phase component on the horizontal axis and the orthogonal component on the vertical axis is created.
- the direction is matched by rotating each wavelength at each point on the fiber, and a new vector is generated by averaging the vectors whose directions are matched, and the in-phase and quadrature components of the generated new vector value is used to calculate the phase.
- Non-Patent Document 4 has been proposed as a solution to the above problem.
- pulses of different optical frequency components are arranged at equal intervals in time and wavelength-multiplexed pulsed light is incident on the optical fiber to be measured, and each wavelength from the optical fiber to be measured Scattered light vector obtained by plotting the scattered light at , on a two-dimensional plane in which the in-phase component is plotted on the horizontal axis and the orthogonal component is plotted on the vertical axis, is created.
- the phase is calculated using the scattered light vector obtained.
- the upper limit of the sampling rate can be set to N ⁇ fs by N-wave multiplexing. Note that the wavelength multiplexing number "N" is an arbitrary natural number.
- Non-Patent Document 4 when performing the frequency multiplexing method as described in Non-Patent Document 4, the angle of the scattered light vector obtained at each optical frequency is simply connected without correcting the angle difference between each optical frequency. If the phase change is calculated by , a problem arises that the calculated phase change is distorted with respect to the actual phase change, and an accurate vibration waveform cannot be measured.
- the single frequency after first calculating the temporal phase difference of each optical frequency, by concatenating the calculated phase differences of each optical frequency, the single frequency Proposed a method for correctly estimating the frequency even for a signal with a vibration frequency exceeding the Nyquist frequency f v in the case of . That is, up to frequency N ⁇ fv can be estimated without aliasing.
- this proposal since the angular difference between the optical frequencies is not obtained, there is a problem that the vibration waveform cannot be measured.
- Non-Patent Document 5 by correcting the angular difference between the respective optical frequencies using the correction frequency, it is possible to measure the vibration waveform under the condition that the upper limit of the sampling rate is increased to N ⁇ fs .
- the compensating optical frequency is used separately from the main optical frequency to improve the sampling rate, and the main frequency component and the compensating optical frequency component are periodically applied to the fiber under test at the same timing.
- the incident probe pulse train is used to correct the angular difference between the main optical frequency components and the compensating optical frequency components, thereby correcting said angular difference between the main optical frequencies.
- Non-Patent Document 5 can measure the vibration waveform, it is also effective in relaxing such restrictions.
- Non-Patent Document 5 a frequency multiplexing method in which pulses of different optical frequencies are incident at different times in order to improve the sampling rate and a frequency multiplexing method described in Patent Document 1 for fading countermeasures are used. , also proposes a method for constructing optical frequency pulses and a method for processing the received signal for simultaneous implementation.
- Non-Patent Document 7 The relationship between the magnitude of the phase change and the amount of strain applied to the fiber due to vibration is described in Non-Patent Document 7, for example.
- Non-Patent Document 7 when a fiber with a total length l is elongated by ⁇ 1 due to a strain amount ⁇ , the amount of increase ⁇ in phase change when light passes due to the amount of expansion by ⁇ 1 is given by the following equation.
- n is the effective refractive index of the fiber
- ⁇ p is Poisson's ratio
- p11 and p12 are the strain-optic tensor components.
- Non-Patent Document 5 pulses of different optical frequencies are injected at different times in order to improve the sampling rate.
- the band occupied by the scattered light of each optical frequency component has a sinc function shape, and side lobes exist. Due to the side lobes, crosstalk occurs in which the bands occupied by the scattered light of each optical frequency component overlap each other. Due to this crosstalk, the signal of the scattered light of a certain optical frequency component fm contains the components in the band near the optical frequency fm among the sidelobes of the scattered light of another optical frequency component fn incident at another time. They are added together and cannot be separated by a digital bandpass filter or the like.
- the above two components are added in a vector state on a plane with the in-phase component of the scattered light on the horizontal axis and the orthogonal component on the vertical axis. Since addition is performed in a vector state, there is a problem that erroneous detection of vibration occurs in sections other than the location where the vibration actually occurred. In addition, a phenomenon occurs in which the vibration waveform at the vibrating portion is distorted with respect to the actual waveform, and there is an effect that the angle difference cannot be completely corrected using the compensating light frequency.
- the present disclosure aims to reduce crosstalk between different optical frequency components, and reduce false detection of vibration and distortion of vibration waveforms.
- the present disclosure suppresses the extent to which the bands occupied by the scattered light of each optical frequency component overlap each other by making the optical pulse profile incident on the optical fiber under test a waveform with small side lobes in the frequency domain, thereby suppressing the crosstalk. To reduce problems such as erroneous detection of vibration in a section other than the location where the vibration actually occurred.
- a waveform such as a Raised-Cosine wave, which has a smaller spectral sidelobe than a rectangular wave, is used as the pulse shape of each optical frequency component.
- the vibration measuring instrument and vibration measuring method of the present disclosure include: In a measuring instrument that performs DAS-P by repeatedly injecting a plurality of optical pulses with different optical frequencies into an optical fiber, As the plurality of optical pulses, the optical pulses having waveforms with smaller spectral side lobes than rectangular waves are used.
- crosstalk between different optical frequency components can be reduced, and erroneous vibration detection and vibration waveform distortion can be reduced.
- FIG. 4 is a diagram illustrating a vibration detection system that detects vibrations in the DAS-P of this embodiment; 4 shows a configuration example of a pulse pattern. An example of crosstalk that occurs when the pulse profile of each optical frequency component is a rectangular waveform is shown. An example of side lobes included in the scattered light signal is shown. An example of a main signal component and a crosstalk component is shown. 4 illustrates an example temporal waveform of each optical frequency pulse according to the present disclosure; 1 shows an example of a pulse train used in the present disclosure; 4 shows an example of optical pulse generation using an SSB modulator.
- FIG. 1 is a diagram for explaining a vibration detection system that detects vibrations with the DAS-P of this embodiment.
- This vibration detection system comprises: a light source for injecting a frequency-multiplexed optical pulse train into one end of an optical fiber to be measured; and a signal processing unit that observes the vibration of the measurement optical fiber as a time change of the phase component of the scattered light.
- the vibration measuring instrument 31 includes a CW light source 1, a coupler 2, an optical modulator 3, a 90-degree optical hybrid 7, and balance detectors (13, 14).
- a CW light source 1, a coupler 2, and an optical modulator 3 correspond to the light source.
- the 90-degree optical hybrid 7 and balance detectors (13, 14) correspond to the light receivers.
- the receiver performs coherent detection using a 90 degree optical hybrid 7 .
- the signal processing device 17 corresponds to the signal processing section. However, it is not always necessary to use a 90-degree optical hybrid for the receiving system, and another device or signal processing may be used as long as the in-phase component and quadrature component of the scattered light can be measured.
- the signal processing device 17 of the present disclosure can also be realized by a computer and a program, and the program can be recorded on a recording medium or provided through a network.
- the vibration measuring device 31 measures the scattered light from the optical fiber 6 to be measured as follows.
- a CW light source 1 emits continuous light of a single wavelength with an optical frequency of f 0 and is split by a coupler 2 into reference light and probe light.
- the probe light is shaped into a wavelength-multiplexed optical pulse 4 by an optical modulator 3 .
- As the optical pulse 4 multiple pulses using compensation optical frequencies as described in Non-Patent Document 5 can be used.
- a configuration example of the optical pulse 4 is shown in FIG.
- N+1 groups arranged in order are prepared.
- N(N+1) delimiting pulse pairs are generated for every M (M is an arbitrary natural number) like 202 in the whole sequence.
- a pulse pair 203 is generated by adding a compensating optical frequency f NM+1 to 202 every N+1 pulse pairs.
- an actually incident light pulse pair train is constructed as indicated by 204.
- FIG. As a result, a pulse pattern is generated in which N(N+1) pulse pairs are arranged in a constant temporal cycle.
- compensating optical frequency f4 is added to the pulse pair with optical frequency f1
- optical frequency f4 is added to a pulse pair with optical frequency f3.
- the limit on the minimum value of how much TN can be reduced by the length of the optical fiber 6 under test is 1/N times that in the case of using a single optical frequency pulse. only mitigated. Also, in the configuration of FIG. 2, fading noise can be reduced using the M pulses present in each pulse pair.
- the type of the optical modulator 3 is not specifically specified as long as it can generate the optical pulse 4, and there may be multiple types.
- an SSB (Single Side Band) modulator or a frequency-variable AO (Acousto-Optics) modulator may be used, and in order to increase the extinction ratio in pulsing, an SOA (Semiconductor Optical Amplifier) or the like may be used. Modulation may be performed.
- the pulse of each optical frequency component shown in 204 has a rectangular wave shape, it is also possible to use a waveform other than a rectangular wave.
- the light pulse 4 is incident on the optical fiber 6 to be measured via the circulator 5 .
- Light scattered at each point in the longitudinal direction of the optical fiber 6 to be measured returns to the circulator 5 as backscattered light and enters one input portion of the 90-degree optical hybrid 7 .
- the reference light split by the coupler 2 enters the other input portion of the 90-degree optical hybrid 7 .
- the internal configuration of the 90-degree optical hybrid 7 can be anything as long as it has the function of a 90-degree optical hybrid.
- a configuration example is shown in FIG.
- the backscattered light is incident on the coupler 8 with a branching ratio of 50:50, and the scattered light split into two is incident on the input portions of the coupler 12 with a branching ratio of 50:50 and the coupler 11 with a branching ratio of 50:50.
- the reference light is incident on the coupler 9 with a branching ratio of 50:50, one of the two-branched reference light is incident on the input part of the coupler 11, and the other is phase-shifted by ⁇ /2 by the phase shifter 10. and is incident on the input of the coupler 12 .
- the two outputs of coupler 11 are detected by balance detector 13 and output is electrical signal 15 which is the analog in-phase component I analog .
- the two outputs of coupler 12 are detected by balance detector 14 and output is electrical signal 16 which is the analog quadrature component Q analog .
- the electric signal 15 and the electric signal 16 are sent to a signal processing device 17 having an AD (Analog to Digital) conversion element 17a and an AD conversion element 17b capable of sampling the frequency band of the signal without aliasing.
- NM+1 can be used.
- a calculation method may be considered in which I digital and Q digital are passed through a band-pass filter whose center frequency is f 0 +f i to compensate for phase delay. If the pulse width of each optical frequency component is W, the passband can be set to 2/W.
- the in-phase component and the quadrature component in the state of an analog electric signal may be separated into respective frequency components by an analog electric filter, and then AD-converted by the AD conversion elements 17a and 17b.
- the signal processing unit 17d calculates the phase. First, a complex vector r i is created on the xy plane with the in-phase component on the x-axis (real number axis) and the quadrature component on the y-axis (imaginary number axis).
- k ⁇ T N +n ⁇ N ⁇ T N (where n is an arbitrary integer) be the time when the head of pulse pair k is incident.
- (1-1) in M different optical frequency bands excluding the compensating optical frequencies constituting the pulse pair By averaging the calculated vectors, the phase at the position of the distance z from the incident end is calculated.
- the state of the optical fiber 6 to be measured at the position of the distance z from the incident end in the longitudinal direction on the optical fiber 6 to be measured is time k ⁇ T N +n ⁇ N ⁇ T N +z/ It is measured by ⁇ (n is any integer).
- ⁇ is the speed of light in the optical fiber 6 to be measured. Furthermore, considering the time for the scattered scattered light to propagate and return to the incident end, the measurement time at the vibration measuring device 31 is k ⁇ T N +n ⁇ N ⁇ T N +2z/ ⁇ (n is an arbitrary integer). becomes. Therefore, the phase calculated at the point of distance z is explicitly represented by the measurement time of the vibration measuring device 31, and
- the phase change due to the vibration applied to the section from the distance z 1 to the distance z 2 on the optical fiber 6 to be measured is the difference between the equations (1-3a) and (1-3b), that is, the equation (1- 3c).
- the instant time when the state of the optical fiber 6 to be measured does not include the time required for the scattered light to return to the incident end as described above, so at the point at the distance z1, the time mT N +z1 / ⁇ , at the point of distance z 2 , the time is mT N +z 2 / ⁇ , and there is a time difference of (z 1 ⁇ z 2 )/ ⁇ .
- the difference in distance between z 1 and z 2 is about the same as the spatial resolution, and is usually set to several meters to several tens of meters. ns, which is very short compared to the scale of time change of normal vibration to be measured, so the difference in time when the state of the optical fiber 6 to be measured is measured can be ignored. Therefore, it is possible to correctly measure the vibration applied to the corresponding section.
- Non-Patent Document 5 proposes a method for correcting the angular difference using a compensating optical frequency. In order to completely correct the angular difference between different optical frequencies, it is necessary to correct the angular difference of the leading optical frequencies of any two pulse pairs.
- the leading optical frequency of pulse pair j is f j pf
- the optical frequency of pulse pair i is f i pf .
- z, f j pf , f i pf ) can be expanded using f NM+1 as follows.
- i, j are arbitrary positive integers. However, i ⁇ j.
- the first pulse pair includes optical frequencies f1 and f4
- the fifth pulse pair includes optical frequencies f2 and f4
- the ninth pulse pair includes optical frequencies f2 and f4 .
- the pulse pair includes optical frequency f3 and optical frequency f4 .
- Equation (1-4) each term on the right side of Equation (1-4) can be calculated based on the same principle as the technique of Patent Document 1.
- the final phase is calculated from ⁇ (z, mT N +2z/ ⁇ ).
- Crosstalk 403 occurs between the side lobes of the scattered light due to the light pulse of the frequency f1 component and the main lobe of the scattered light due to the light pulse of the light frequency f3 component. Since the crosstalk 403 exists, the signal actually separated as the optical frequency component f3 is a superposition of the crosstalk component and the main signal component, as indicated by 305 in FIG. 3 , and the two cannot be separated.
- the value of the distance from the incident end of the crosstalk component in equation (1-5) is l', which is different from l.
- the reason for this is that the crosstalk component is received at a different timing than the main signal component, so the vector corresponding to the position 308 of the main signal component overlaps the vector corresponding to the different position 309 of the crosstalk component. be.
- Position 308 is in front of vibrating position 306, while position 309 is behind vibrating position 306, as shown in positional relationship 310 on the fiber in FIG. As shown in FIG. 5, considering the change from the state 501 at time T to the state 502 at time T', the position 308 is before the vibration position 306, so the vector of the main signal component does not change.
- the crosstalk component undergoes an angular change due to the vibration.
- the observation vector also changes in angle. If the angle of the observation vector r 3 (l, T) is ⁇ 3 (l, T), the amount of change in the angle is given by Equation (1-6).
- the amount of change in formula (1-6) depends not only on the magnitude of the vibration occurring at the vibration position 306, but also on the lengths of the main signal component vector and the crosstalk component vector, and the angular difference between the two vectors. do. Considering that the length and angle of the vector change for each position due to the fading effect, the change amount of the angle in Equation (1-6) changes depending on the distance l. For example, an angular change at a position .delta.l behind the position 308 can be similarly considered as long as it is within the overlapping section 307. FIG.
- This embodiment proposes a method of solving the above-mentioned problem due to the above-mentioned crosstalk by making the pulse shape of each optical frequency component have a smoothly changing outline.
- 610 in FIG. 6 shows the time waveform of each optical frequency pulse used in this proposal. It shows the outline of the electric field envelope before being multiplied by the carrier frequency.
- Waveform 601 is the shape of a conventional square wave, also used in Non-Patent Document 5, and waveform 602 represents the shape of the waveform used in this disclosure.
- the waveform shape is not limited to a specific shape as long as it is a waveform capable of suppressing crosstalk.
- a raised cosine waveform can be used.
- a waveform 602 shows a raised cosine waveform as an example.
- the spatial resolution of the vibration measuring apparatus can be maintained at the same level and deterioration of the SN ratio of the signal can be prevented.
- Power density in the frequency domain is shown at 620 .
- 605 in the case of a raised cosine waveform suppresses side lobes, and it can be seen that the magnitude of the crosstalk can be reduced. Therefore, the pulse train when using the raised cosine waveform has outline 701 in FIG. By using the pulse train shown in 701, it is possible to reduce the vector length of the crosstalk component.
- vibrations can be observed at locations other than the vibrations caused by the crosstalk without changing the signal processing method for phase calculation. It is possible to reduce the problem that the vibration waveform at the vibration location is distorted with respect to the actual waveform, and the influence that the angle difference cannot be completely corrected using the compensation light frequency.
- FIG. 8 An example of an optical pulse generation method using an SSB modulator will be described with reference to FIG. This is an example in which a modulation signal 802 is given to an SSB modulator 801 to modulate continuous light to create an optical pulse pair 803 . If the modulation signal 802 to be applied is set so that the envelope of each modulation frequency component has the contour of a raised cosine waveform, the magnitude of the electric field of each optical frequency component can also have a contour that changes smoothly.
- the generated pulse pair 803 has an optical frequency of f 0 +f 1 where f 0 is the frequency of the original continuous light.
- f 0 + f 2 and f 0 + f 3 can be approximately regarded as a pulse pair in which f 0 + f 2 and f 0 + f 3 are arranged, and the signal processing unit performs signal processing by coherently detecting the scattered light generated by the incident pulse pair. It is also sufficient to set the center frequencies of the bandpass filters for extracting the respective optical frequency components at 17c to f1 , f2 and f3, respectively .
- the FWHM (Half Width at Half Maximum) of the raised cosine waveform which is the envelope of each modulation frequency component of the modulated signal 802, is approximately the FWHM of each optical frequency component of the pulse pair electric field magnitude 803. It can be set assuming that they match.
- the passband width of the bandpass filter in the signal processing unit 17c is determined by the FWHM of the raised cosine waveform.
- the resolution can also be determined by adjusting the FWHM of the raised cosine waveform in consideration of the characteristics of SSB.
- each main optical frequency in frequency multiplexing if possible, the position of the side lobe trough of each main optical frequency shown in 606 in FIG.
- the influence of crosstalk due to the side lobes can be reduced.
- optical pulse pairs applicable to the present disclosure are not limited to the configuration shown in FIG. 2, and can be used for general techniques that improve the sampling rate using frequency multiplexing.
- the pulse shape of each optical frequency component is a Raised-Cosine wave instead of a conventional rectangular wave.
- This disclosure can be applied to the information and communications industry.
- CW light source 2 coupler 3: optical modulator 4: optical pulse 5: circulator 6: optical fiber to be measured 7: 90 degree optical hybrid 8, 9: coupler 10: phase shifter 11, 12: coupler 13, 14: balance Detector 15: analog electrical signal of in-phase component 16: electrical signal of analog quadrature component 17: signal processors 17a, 17b: AD conversion elements 17c, 17d: signal processor 31: vibration measuring instrument
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
Description
光周波数の異なる複数の光パルスを光ファイバに繰り返し入射し、DAS-Pを行う測定器において、
前記複数の光パルスとして、矩形波と比較してスペクトルサイドローブの小さい波形を有する前記光パルスを用いる。
図1は、本実施形態のDAS-Pで振動検出を行う振動検出システムを説明する図である。本振動検出システムは、周波数多重した光パルス列を被測定光ファイバの一端に入射する光源と、前記被測定光ファイバの前記一端に戻ってきた各波長の散乱光を受光する受光器と、前記被測定光ファイバの振動を前記散乱光の位相成分の時間変化として観測する信号処理部と、を備える。
本開示によれば、非特許文献5で記載される異なる時刻に異なる光周波数成分を入射する光周波数多重型位相OTDRにおいて、各光周波数成分のパルス外形として従来の矩形波に代わりRaised-Cosine波に代表されるスペクトルサイドローブの小さい波形を使用することで、異なる光周波数成分の間のクロストークを小さくし、振動箇所以外での振動の誤検知を抑えつつ振動箇所における振動波形が実際の振動波形に対して歪む現象の発生、補償光周波数を用いた前記角度差の補正が完全に行えない影響を低減する。
2:カプラ
3:光変調器
4:光パルス
5:サーキュレータ
6:被測定光ファイバ
7:90度光ハイブリッド
8、9:カプラ
10:位相シフタ
11、12:カプラ
13、14:バランス検出器
15:アナログの同相成分の電気信号
16:アナログの直交成分の電気信号
17:信号処理装置
17a、17b:AD変換素子
17c、17d:信号処理部
31:振動測定器
Claims (7)
- 光周波数の異なる複数の光パルスを光ファイバに繰り返し入射し、DAS-P(Distributed Acoustic Sensing-phase)を行う振動測定器において、
前記複数の光パルスとして、矩形波と比較してスペクトルサイドローブの小さい波形を有する前記光パルスを用いる、
振動測定器。 - 前記複数の光パルスは、メインの光周波数の成分と補償光周波数の成分を含む複数のパルス対からなり、
各メインの光周波数のサイドローブ成分の谷に位置する周波数位置に、その他のパルス対に含まれるメインの光周波数のメインローブの中心が配置されている、
請求項1に記載の振動測定器。 - 連続光を変調信号で変調することによって前記複数の光パルスを生成する、シングルサイドバンド変調器を備える、
請求項1又は2に記載の振動測定器。 - 前記変調信号の各変調周波数成分の包絡線がRaised cosine波形を有する、
請求項3に記載の振動測定器。 - 前記変調信号の変調周波数成分を用いて、前記複数の光パルスが前記光ファイバで散乱されて発生した散乱光を、前記複数の光パルスごとに分離する、
請求項3又は4に記載の振動測定器。 - 前記複数の光パルスごとに分離するバンドパスフィルタにおける通過帯域幅がRaised cosine波形のFWHM(Half Width at Half Maximum)で定められている、
請求項3から5のいずれかに記載の振動測定器。 - 光周波数の異なる複数の光パルスを光ファイバに繰り返し入射し、DAS-P(Distributed Acoustic Sensing-phase)を行う振動測定方法において、
前記複数の光パルスとして、矩形波と比較してスペクトルサイドローブの小さい波形を有する前記光パルスを用いる、
振動測定方法。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2021/021991 WO2022259437A1 (ja) | 2021-06-09 | 2021-06-09 | 振動測定器及び振動測定方法 |
EP21945113.5A EP4354100A1 (en) | 2021-06-09 | 2021-06-09 | Vibration measurement device and vibration measurement method |
JP2023526737A JPWO2022259437A1 (ja) | 2021-06-09 | 2021-06-09 | |
CN202180098778.9A CN117480364A (zh) | 2021-06-09 | 2021-06-09 | 振动测量器及振动测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2021/021991 WO2022259437A1 (ja) | 2021-06-09 | 2021-06-09 | 振動測定器及び振動測定方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2022259437A1 true WO2022259437A1 (ja) | 2022-12-15 |
Family
ID=84425962
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/JP2021/021991 WO2022259437A1 (ja) | 2021-06-09 | 2021-06-09 | 振動測定器及び振動測定方法 |
Country Status (4)
Country | Link |
---|---|
EP (1) | EP4354100A1 (ja) |
JP (1) | JPWO2022259437A1 (ja) |
CN (1) | CN117480364A (ja) |
WO (1) | WO2022259437A1 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117392796A (zh) * | 2023-12-11 | 2024-01-12 | 云南保利天同水下装备科技有限公司 | 分区探测方法、分区探测***及其防御探测组件 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014157134A (ja) * | 2013-02-18 | 2014-08-28 | Nippon Telegr & Teleph Corp <Ntt> | 光パルス試験装置及び光パルス試験方法 |
JP2020169904A (ja) | 2019-04-03 | 2020-10-15 | 日本電信電話株式会社 | 位相測定方法及び信号処理装置 |
WO2021075015A1 (ja) * | 2019-10-17 | 2021-04-22 | 日本電信電話株式会社 | 光パルス試験方法及び光パルス試験装置 |
-
2021
- 2021-06-09 EP EP21945113.5A patent/EP4354100A1/en active Pending
- 2021-06-09 WO PCT/JP2021/021991 patent/WO2022259437A1/ja active Application Filing
- 2021-06-09 CN CN202180098778.9A patent/CN117480364A/zh active Pending
- 2021-06-09 JP JP2023526737A patent/JPWO2022259437A1/ja active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014157134A (ja) * | 2013-02-18 | 2014-08-28 | Nippon Telegr & Teleph Corp <Ntt> | 光パルス試験装置及び光パルス試験方法 |
JP2020169904A (ja) | 2019-04-03 | 2020-10-15 | 日本電信電話株式会社 | 位相測定方法及び信号処理装置 |
WO2021075015A1 (ja) * | 2019-10-17 | 2021-04-22 | 日本電信電話株式会社 | 光パルス試験方法及び光パルス試験装置 |
Non-Patent Citations (10)
Title |
---|
A. E. ALEKSEEV ET AL.: "Fidelity of the dual-pulse phase-OTDR response to spatially distributed external perturbation", LASER PHYS., vol. 29, no. 055106 |
ALI. MASOUDIT. P. NEWSON: "Contributed Review: Distributed optical fibre dynamic strain sensing", REVIEW OF SCIENTIFIC INSTRUMENTS, vol. 87, 2016, pages 011501 |
C. D. BUTTERG. B. HOCKER: "Fiber optics strain gauge", APPL. OPT., vol. 17, 1978, pages 2867 - 2869 |
D. IIDAK. TOGET. MANABE: "Distributed measurement of acoustic vibration location with frequency multiplexed phase-OTDR", OPT. FIBER TECHNOL., vol. 36, 2017, pages 19 - 25, XP085114643, DOI: 10.1016/j.yofte.2017.02.005 |
G. YANG ET AL.: "Long-Range Distributed Vibration Sensing Based on Phase Extraction from Phase-Sensitive OTDR", IEEE PHOTONICS JOURNAL, vol. 8, no. 3, 2016, XP011607611, DOI: 10.1109/JPHOT.2016.2552820 |
IIDA, DAISUKE ET AL.: "Distributed measurement of acoustic vibration location with frequency multiplexed phase-OTDR", OPTICAL FIBER TECHNOLOGY, vol. 36, 21 February 2017 (2017-02-21), pages 19 - 25, XP085114643, DOI: 10.1016/j.yofte.2017.02.005 * |
KEN'ICHI NISHIGUCHILI CHE-HSIENARTUR GUZIKMITSUNORI YOKOYAMAKINZO KISHIDA: "Fabrication of Fiber-Optic Distributed Acoustic Sensor and Its Signal Processing", IEICE TECHNICAL REPORT, vol. 115, no. 202, 2015, pages 29 - 34 |
MARIA ROSARIO FERNANDEZ-RUIZHUGO F. MARTINS: "Steady-Sensitivity Distributed Acoustic Sensors", J. LIGHTWAVE TECHNOL., vol. 36, 2018, pages 5690 - 5696, XP011697725, DOI: 10.1109/JLT.2018.2877849 |
WAKISAKA, YOSHIFUMI ET AL.: "Distortion-suppressed sampling rate enhancement in phase-OTDR vibration sensing with newly designed FDM pulse sequence for correctly monitoring various waveforms", 2020 OPTICAL FIBER COMMUNICATIONS CONFERENCE AND EXHIBITION(OFC, 8 March 2020 (2020-03-08), pages 1 - 3, XP033767147, DOI: 10.1364/OFC.2020.Th3F.5 * |
Y. WAKISAKAD. IIDAH. OSHIDA: "Distortion-Suppressed Sampling Rate Enhancement in Phase-OTDR Vibration Sensing with Newly Designed FDM Pulse Sequence for Correctly Monitoring Various Waveforms", 2020 OPTICAL FIBER COMMUNICATIONS CONFERENCE AND EXHIBITION (OFC, 2020, pages 1 - 3, XP033767147, DOI: 10.1364/OFC.2020.Th3F.5 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117392796A (zh) * | 2023-12-11 | 2024-01-12 | 云南保利天同水下装备科技有限公司 | 分区探测方法、分区探测***及其防御探测组件 |
CN117392796B (zh) * | 2023-12-11 | 2024-03-22 | 云南保利天同水下装备科技有限公司 | 分区探测方法、分区探测***及其防御探测组件 |
Also Published As
Publication number | Publication date |
---|---|
CN117480364A (zh) | 2024-01-30 |
JPWO2022259437A1 (ja) | 2022-12-15 |
EP4354100A1 (en) | 2024-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Muanenda | Recent advances in distributed acoustic sensing based on phase‐sensitive optical time domain reflectometry | |
Xiong et al. | Single-shot COTDR using sub-chirped-pulse extraction algorithm for distributed strain sensing | |
US11867555B2 (en) | Phase measurement method and signal processing device | |
JP6893137B2 (ja) | 光ファイバ振動検知センサおよびその方法 | |
JP7435160B2 (ja) | 光ファイバ振動検知装置及び振動検知方法 | |
US11725965B2 (en) | Method and device for reconstructing a backscattered electromagnetic vector wave | |
JPWO2018083732A1 (ja) | ブリルアン散乱測定方法およびブリルアン散乱測定装置 | |
JP7298706B2 (ja) | 光パルス試験方法及び光パルス試験装置 | |
CN114543973A (zh) | 一种分布式超高频振动信号测量方法及光纤传感器 | |
JP6751379B2 (ja) | 光時間領域反射測定方法および光時間領域反射測定装置 | |
WO2022259437A1 (ja) | 振動測定器及び振動測定方法 | |
US11522606B2 (en) | Phase measurement method, signal processing device, and program | |
EP3922964B1 (en) | Vibration detection method, signal processing device, and program | |
WO2020194856A1 (ja) | 光コヒーレントセンサ及び光コヒーレントセンシング方法 | |
WO2022259436A1 (ja) | 信号処理装置、振動検出システム及び信号処理方法 | |
WO2023053263A1 (ja) | 光パルス試験方法及び光パルス試験装置 | |
Wu et al. | High performance distributed acoustic sensor based on ultra-weak FBG array | |
JP5927079B2 (ja) | レーザ光特性測定方法及び測定装置 | |
JP6751378B2 (ja) | 光時間領域反射測定方法および光時間領域反射測定装置 | |
WO2022059050A1 (ja) | 光周波数領域反射計測装置及び方法 | |
JP7435772B2 (ja) | 歪変化計測装置及び歪変化計測方法 | |
EP4047333B1 (en) | Optical pulse testing method and optical pulse testing device | |
US20240053172A1 (en) | Optical fiber sensor and brillouin frequency shift measurement method | |
CN117242320A (zh) | 信号处理方法及信号处理装置 | |
CN114199514A (zh) | 基于光频域反射分布式传感的假峰消除方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 21945113 Country of ref document: EP Kind code of ref document: A1 |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2023526737 Country of ref document: JP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 202180098778.9 Country of ref document: CN |
|
WWE | Wipo information: entry into national phase |
Ref document number: 18567606 Country of ref document: US |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2021945113 Country of ref document: EP |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
ENP | Entry into the national phase |
Ref document number: 2021945113 Country of ref document: EP Effective date: 20240109 |