CN116399379B - Distributed optical fiber acoustic wave sensing system and measuring method thereof - Google Patents

Distributed optical fiber acoustic wave sensing system and measuring method thereof Download PDF

Info

Publication number
CN116399379B
CN116399379B CN202310664799.3A CN202310664799A CN116399379B CN 116399379 B CN116399379 B CN 116399379B CN 202310664799 A CN202310664799 A CN 202310664799A CN 116399379 B CN116399379 B CN 116399379B
Authority
CN
China
Prior art keywords
information
optical fiber
dimensional image
dimensional
demodulation
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
Application number
CN202310664799.3A
Other languages
Chinese (zh)
Other versions
CN116399379A (en
Inventor
渠帅
王伟涛
尚盈
王晨
黄胜
李常
曹冰
赵文安
倪家升
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Laser Institute of Shandong Academy of Science
Original Assignee
Laser Institute of Shandong Academy of Science
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Laser Institute of Shandong Academy of Science filed Critical Laser Institute of Shandong Academy of Science
Priority to CN202310664799.3A priority Critical patent/CN116399379B/en
Publication of CN116399379A publication Critical patent/CN116399379A/en
Application granted granted Critical
Publication of CN116399379B publication Critical patent/CN116399379B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING 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/00Mechanical 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/26Mechanical 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/32Mechanical 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/34Mechanical 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/353Mechanical 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/35338Mechanical 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/35354Sensor working in reflection
    • G01D5/35358Sensor working in reflection using backscattering to detect the measured quantity
    • G01D5/35361Sensor working in reflection using backscattering to detect the measured quantity using elastic backscattering to detect the measured quantity, e.g. using Rayleigh backscattering
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING 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/00Mechanical 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/26Mechanical 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/32Mechanical 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/34Mechanical 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/353Mechanical 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/3537Optical fibre sensor using a particular arrangement of the optical fibre itself

Abstract

The embodiment of the application provides a distributed optical fiber acoustic wave sensing system and a measuring method thereof, which relate to the technical field of distributed optical fiber sensing detection and comprise the steps of collecting backward Rayleigh scattering signals generated by vibration signals in a sensing optical fiber along a line in a preset time period; demodulating the collected backward Rayleigh scattering signal to obtain first demodulation information; constructing an original two-dimensional image through the first demodulation information, wherein the original two-dimensional image comprises an image of a time-distance domain formed by the first demodulation information corresponding to each position of the sensing optical fiber along the line in a preset time period; performing two-dimensional wavelet transformation image processing on the original two-dimensional image to eliminate noise information and obtain a noise-removed two-dimensional image; and reconstructing second demodulation information according to the denoising two-dimensional image, wherein the second demodulation information comprises demodulation information after the denoising information of each position along the sensing optical fiber, so that measurement errors are reduced, and the measurement accuracy of the distributed optical fiber acoustic wave sensing system is improved.

Description

Distributed optical fiber acoustic wave sensing system and measuring method thereof
Technical Field
The application relates to the technical field of distributed optical fiber sensing detection, in particular to a distributed optical fiber acoustic wave sensing system and a measuring method thereof.
Background
The principle of the distributed optical fiber acoustic wave sensing technology is that the change of physical quantity (such as sound, vibration and the like) required to be detected is obtained by detecting the phase change of Rayleigh scattered light along the sensing optical fiber. The distributed optical fiber acoustic wave sensing system has the advantages of light weight, small volume, high sensitivity, strong electromagnetic interference resistance and the like, and has extremely high application value in the fields of perimeter security, oil and gas exploration, pipeline monitoring and the like.
In the related art, in the measuring process of the distributed optical fiber acoustic wave sensing system, incident light can interfere with backward scattered light at a certain point along the sensing optical fiber, and the change of physical quantities such as acoustic wave or vibration at the point can cause the change of the phase of interference light, so that the change of acoustic wave or vibration at the point can be determined by demodulating the phase of interference light. Because the sensing optical fibers are continuously distributed in the space, the distributed optical fiber acoustic wave sensing system can quantitatively detect the change of the physical quantity of any point in the space, thereby realizing distributed sensing.
However, in the distributed optical fiber acoustic wave sensing system, when the information of each position along the transmission fiber is demodulated under the influence of random noise fluctuation such as laser phase noise, polarization noise, environmental noise, and the like, a plurality of abnormal information are contained in the obtained demodulation result. In addition, since the demodulation result is obtained from the noisy information, the time-domain waveform and amplitude fluctuation of the demodulation result are also large, thereby bringing about measurement errors.
Disclosure of Invention
The embodiment of the application provides a measuring method of a distributed optical fiber acoustic wave sensor system, which aims to solve the technical problem that a distributed optical fiber acoustic wave sensor in the prior art has larger measuring error.
The embodiment of the application provides a measuring method of a distributed optical fiber acoustic wave sensing system, which comprises the steps of collecting backward Rayleigh scattering signals generated by vibration signals in a sensing optical fiber along a line in a preset time period;
demodulating the collected backward Rayleigh scattering signals to obtain first demodulation information, wherein the first demodulation information comprises demodulation information containing noise information at each position along the sensing optical fiber;
constructing an original two-dimensional image through the first demodulation information, wherein the original two-dimensional image comprises an image of a time-distance domain formed by the first demodulation information corresponding to each position of the sensing optical fiber along the line in a preset time period;
performing two-dimensional wavelet transformation image processing on the original two-dimensional image to eliminate noise information and obtain a noise-removed two-dimensional image;
reconstructing second demodulation information according to the denoising two-dimensional image, wherein the second demodulation information comprises demodulation information after the denoising information of each position along the sensing optical fiber.
In one possible implementation, the original two-dimensional image is subjected to two-dimensional wavelet transform image processing, including,
and carrying out wavelet decomposition on the original two-dimensional image through a discrete wavelet transform function to obtain a plurality of sub-images corresponding to the original two-dimensional image.
In one possible implementation, the two-dimensional wavelet transform image processing is performed on an original two-dimensional image, further comprising,
the level of wavelet decomposition is adjusted to obtain a corresponding number of sub-images, and the high-frequency component and the low-frequency component of the sub-images are acquired and processed to eliminate the abnormal value on the sub-images so as to obtain the denoising sub-images.
In one possible implementation, there is an adjustment threshold for the level of wavelet decomposition, within which the higher the level of wavelet decomposition, the less outliers the noise-removed sub-image is obtained.
In one possible implementation, the two-dimensional wavelet transform image processing is performed on the original two-dimensional image, and the method further comprises the step of performing wavelet reconstruction on the high-frequency component and the low-frequency component of the denoising sub-image by adopting an inverse discrete wavelet transform function to obtain a denoising two-dimensional image.
In one possible implementation, after reconstructing the second demodulation information from the denoised two-dimensional image, extracting time domain information of the vibration signal in the second demodulation information is included.
In one possible implementation, the original two-dimensional image and the noise-removed two-dimensional image each include two-dimensional graphical information of a time-distance domain or two-dimensional matrix information of a time-distance domain.
The embodiment of the application also provides a distributed optical fiber acoustic wave sensing system, which comprises a demodulation unit and a controller, wherein the demodulation unit is connected with the controller through signals;
the demodulation unit is configured to collect backward Rayleigh scattering signals within a preset time period, wherein the backward Rayleigh scattering signals comprise signals generated by a vibration signal in the distributed optical fiber acoustic wave sensing system;
the demodulation unit is further configured to demodulate the collected backward Rayleigh scattering signals to obtain first demodulation information, wherein the first demodulation information comprises demodulation information containing noise information at each position along the sensing optical fiber;
the controller is configured to construct an original two-dimensional image through the first demodulation information, wherein the two-dimensional image comprises a time-distance domain image formed by the first demodulation information corresponding to each position of the sensing optical fiber along the line in a preset time period;
the controller is further configured to perform two-dimensional wavelet transform image processing on the original two-dimensional image to eliminate noise information to obtain a noise-removed two-dimensional image;
the controller is further configured to reconstruct second demodulation information from the denoised two-dimensional image, the second demodulation information including the noise-cancelled demodulation information for each location along the sensing fiber.
The embodiment of the application provides a measuring method of a distributed optical fiber acoustic wave sensing system, which is characterized in that first demodulation information is obtained by demodulating a backward Rayleigh scattering signal, an original two-dimensional image is constructed by the first demodulation information, a two-dimensional wavelet transformation image processing is carried out on the original two-dimensional image to obtain a denoising two-dimensional image, a second demodulation information is obtained by reconstructing the denoising two-dimensional image, demodulation information after noise information is eliminated is obtained, abnormal information in the demodulation information is eliminated, time domain waveform and amplitude fluctuation of a demodulation result are reduced, measuring errors are reduced, and measuring precision of the distributed optical fiber acoustic wave sensing system is improved.
The embodiment of the application also provides a distributed optical fiber acoustic wave sensing system, which adopts the measuring method in any one of the technical schemes, so that the distributed optical fiber acoustic wave sensing system has all the beneficial effects of the measuring method in any one of the technical schemes, and is not repeated here.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this specification, illustrate embodiments of the application and together with the description serve to explain the application and do not constitute a limitation on the application. In the drawings:
FIG. 1 is a schematic diagram of a distributed optical fiber acoustic wave sensor system according to an embodiment of the present application;
FIG. 2 is a diagram showing steps of a method for measuring a distributed optical fiber acoustic wave sensor system according to an embodiment of the present application;
FIG. 3 is a diagram of an implementation step of S400 in FIG. 2;
FIG. 4 is a two-dimensional image corresponding to demodulation information of vibration signals measured by a measuring method of a distributed optical fiber acoustic wave sensing system according to the related art;
FIG. 5 is a two-dimensional image corresponding to demodulation information of a vibration signal measured by a measuring method of a distributed optical fiber acoustic wave sensing system according to an embodiment of the present application;
FIG. 6 is a time domain diagram of vibration positions corresponding to the two-dimensional image of FIG. 4;
fig. 7 is a time domain diagram of vibration positions corresponding to the two-dimensional image in fig. 5.
Reference numerals illustrate:
1-a laser emitting unit; 2-a circulator; a 3-sensing unit; a 4-demodulation unit;
a 101-laser; 102-an isolator; 103-an acousto-optic modulator; 1031-a function generator; 104-a first erbium-doped fiber amplifier; 105-a first filter;
301-sensing optical fibers; 302-a first piezoelectric ceramic;
401-a second erbium-doped fiber amplifier; 402-a second filter; 403-couplers; 404-a first faraday rotator mirror; 405-a second faraday rotator mirror; 406-a second piezoelectric ceramic; 407-detector; 408-acquisition card.
Detailed Description
The principle of the distributed optical fiber acoustic wave sensing technology is to obtain the change of the physical quantity (such as sound, vibration, etc.) to be detected by detecting the phase change of the rayleigh scattered light along the sensing optical fiber 301. The distributed optical fiber acoustic wave sensing system has the advantages of light weight, small volume, high sensitivity, strong electromagnetic interference resistance and the like, and has extremely high application value in the fields of perimeter security, oil and gas exploration, pipeline monitoring and the like.
In the related art, in the measurement process of the distributed optical fiber acoustic wave sensing system, incident light interferes with back scattered light at a point along the sensing optical fiber 301, and the change of physical quantity such as acoustic wave or vibration at the point causes the change of the phase of interference light, so that the change of acoustic wave or vibration at the point can be determined by demodulating the phase of interference light. Since the sensing optical fibers 301 are continuously distributed in space, the distributed optical fiber acoustic wave sensing system can quantitatively detect the change of the physical quantity at any point in space, thereby realizing distributed sensing.
However, in the distributed optical fiber acoustic wave sensing system, when the information of each position along the transmission fiber is demodulated under the influence of random noise fluctuation such as laser phase noise, polarization noise, environmental noise, and the like, a plurality of abnormal information are contained in the obtained demodulation result. In addition, since the demodulation result is obtained from the noisy information, the time-domain waveform and amplitude fluctuation of the demodulation result are also large, thereby bringing about measurement errors.
Therefore, the embodiment of the application provides a distributed optical fiber acoustic wave sensing system and a measuring method thereof, which are used for solving the technical problem that the measuring error is caused by the fact that the obtained demodulation result contains a plurality of abnormal information when the distributed optical fiber acoustic wave sensing system demodulates the information of all positions along the transmission optical fiber in the prior art.
Fig. 1 is a schematic structural diagram of a distributed optical fiber acoustic wave sensor system according to an embodiment of the present application.
In some examples, referring to fig. 1, the distributed optical fiber acoustic wave sensing system generally includes a laser emitting unit 1, a circulator 2, a sensing unit 3, a demodulation unit 4, and a controller, where an output terminal of the laser emitting unit 1 is connected to a first input terminal of the circulator 2, a first output terminal of the circulator 2 is connected to an output terminal of the sensing unit 3, a second output terminal of the circulator 2 is connected to the demodulation unit 4, and an output terminal of the demodulation unit 4 is connected to the controller.
Illustratively, the laser emitting unit 1 emits a laser signal, the laser signal is output to the sensing unit 3 via the first input end and the first output end of the circulator 2, the laser signal propagates along the sensing optical fiber 301 of the sensing unit 3 to detect a vibration signal in an environment to be detected, and generates backward rayleigh scattered light, the backward rayleigh scattered light is output to the demodulating unit 4 via the second output end of the circulator 2, the demodulating unit 4 demodulates the backward rayleigh scattered light, and the demodulated information is output to the controller.
Fig. 2 is a diagram of implementation steps of a measurement method of a distributed optical fiber acoustic wave sensor system according to an embodiment of the present application. Fig. 3 is a diagram of an implementation step of S400 in fig. 2.
The embodiment of the application provides a measuring method of a distributed optical fiber acoustic wave sensing system, referring to fig. 2 and 3, comprising the following steps:
collecting backward Rayleigh scattering signals generated by vibration signals in the line of the sensing optical fiber 301 in a preset time period;
demodulating the collected backward Rayleigh scattering signals to obtain first demodulation information, wherein the first demodulation information comprises demodulation information containing noise information at each position along the sensing optical fiber 301;
constructing an original two-dimensional image by the first demodulation information, wherein the original two-dimensional image comprises an image of a time-distance domain formed by the first demodulation information corresponding to each position of the sensing optical fiber 301 along the line in a preset time period;
performing two-dimensional wavelet transformation image processing on the original two-dimensional image to eliminate noise information and obtain a noise-removed two-dimensional image;
the second demodulation information is reconstructed from the denoised two-dimensional image, and includes the demodulation information after the noise information is removed at each position along the sensing fiber 301.
The embodiment of the application provides a measuring method of a distributed optical fiber acoustic wave sensing system, which is characterized in that first demodulation information is obtained by demodulating a backward Rayleigh scattering signal, an original two-dimensional image is constructed by the first demodulation information, a two-dimensional wavelet transformation image processing is carried out on the original two-dimensional image to obtain a denoising two-dimensional image, a second demodulation information is obtained by reconstructing the denoising two-dimensional image, demodulation information after noise information is eliminated is obtained, abnormal information in the demodulation information is eliminated, time domain waveform and amplitude fluctuation of a demodulation result are reduced, measuring errors are reduced, and measuring precision of the distributed optical fiber acoustic wave sensing system is improved.
In some examples, performing two-dimensional wavelet transform image processing on the original two-dimensional image includes performing wavelet decomposition on the original two-dimensional image by a discrete wavelet transform function to obtain a plurality of sub-images corresponding to the original two-dimensional image.
For example, the wavelet decomposition may be calculated by referring to the following formula:
wherein, the liquid crystal display device comprises a liquid crystal display device,represents a scaling factor->Represents a translation factor->Representing vibration signal +.>The function of the mother wavelet is represented,is->Complex conjugate of->Representing a wavelet decomposition function.
The first sub-image is obtained after the primary wavelet decomposition of the original two-dimensional image, wherein the first sub-image is averagely divided into a low-frequency component and three high-frequency components, the three high-frequency components comprise three parts of horizontal high frequency, vertical high frequency and diagonal high frequency, the secondary decomposition is carried out on the original two-dimensional image, namely, the area of the low-frequency component of the first sub-image is continuously averagely divided into the low-frequency component and the three high-frequency components, and the multi-stage decomposition is carried out on the original two-dimensional image, so that the method is just like.
The low frequency component after wavelet decomposition represents contour information of an image and can be regarded as effective information.
According to the embodiment of the application, the original two-dimensional image is subjected to wavelet decomposition through the discrete wavelet transform function, so that a plurality of sub-images corresponding to the original two-dimensional image can be obtained, and further, the high-frequency component and the low-frequency component of each sub-image are obtained, so that the sub-images can be subjected to noise removal.
The method includes the steps of performing two-dimensional wavelet transform image processing on an original two-dimensional image, and adjusting the level of wavelet decomposition to obtain a corresponding number of sub-images, and acquiring and processing high-frequency components and low-frequency components of the sub-images to eliminate abnormal values on the sub-images and remove noise of the sub-images.
Illustratively, the level of wavelet decomposition is adjusted to obtain a low-frequency component and a high-frequency component of the sub-image, the pixel value of the low-frequency component is increased according to a nonlinear stretching rule, and noise reduction processing is performed on the high-frequency component to obtain a noise-removed sub-image.
Illustratively, there is an adjustment threshold for the level of wavelet decomposition, within which the higher the level of wavelet decomposition, the less outliers the noise-removed sub-image is obtained.
When the level of wavelet decomposition increases, the number of low-frequency components and the number of high-frequency components increase, and noise of the sub-image increases due to the increase of the high-frequency components, so that noise reduction processing is required for more high-frequency components, and there is a residual outlier in the high-frequency components after the noise reduction processing, and therefore, when the high-frequency components increase to a certain value, the residual outlier increases accordingly, and therefore, within the adjustment threshold, the higher the level of wavelet decomposition, the less the outlier of the noise-removed sub-image is obtained.
According to the embodiment of the application, the sub-image in the adjustment threshold is obtained by adjusting the wavelet decomposition level, and the low-frequency component and the high-frequency component of the sub-image are processed, so that the noise interference in the sub-image is eliminated, and the noise-removed sub-image is obtained.
The method for processing the two-dimensional wavelet transform image of the original two-dimensional image further comprises the step of carrying out wavelet reconstruction on high-frequency components and low-frequency components of the denoising sub-image by adopting an inverse discrete wavelet transform function to obtain the denoising two-dimensional image.
The wavelet reconstruction is, for example, image reconstruction by using the low-frequency component with increased pixel value and the high-frequency component with noise reduction, so as to obtain a contrast enhanced image, i.e. an undenoised two-dimensional image.
According to the embodiment of the application, through wavelet reconstruction, the high-frequency component and the low-frequency component of the denoising sub-image are subjected to wavelet reconstruction by adopting an inverse discrete wavelet transform function, so that a denoising two-dimensional image with the noises eliminated is obtained, and the subsequent reconstruction of second demodulation information is facilitated.
In some examples, reconstructing the second demodulation information from the denoised two-dimensional image includes extracting time domain information of the vibration signal in the second demodulation information.
In some examples, the original two-dimensional image and the denoised two-dimensional image each include two-dimensional graphical information of a time-distance domain or two-dimensional matrix information of a time-distance domain.
It should be noted that, the time in the time-distance domain corresponds to a preset time period, and the distance corresponds to each position along the sensing optical fiber 301.
In some examples, referring to fig. 1, the laser emitting unit 1 includes a laser 101, an isolator 102, an acousto-optic modulator 103, a first erbium-doped amplifier, and a first filter 105, an output of the laser 101 is connected to an input of the isolator 102, an output of the isolator 102 is connected to an input of the acousto-optic modulator 103, an output of the acousto-optic modulator 103 is connected to an input of the erbium-doped amplifier, and an output of the erbium-doped amplifier is connected to an input of the first filter 105. The laser 101 is configured to emit a laser signal to an input of the isolator 102; the isolator 102 is configured such that the laser beam is output unidirectionally by the isolator 102 to the acousto-optic modulator 103. The acousto-optic modulator 103 is connected to a function generator 1031, the function generator 1031 being configured to drive the acousto-optic modulator 103 to modulate the pulse width between pulses of the laser signal and the number of frequencies between pulses. The first erbium-doped fiber amplifier 104 is configured to amplify the laser signal output by the output terminal of the acousto-optic modulator 103. The first filter 105 is configured to filter the disturbance light of the laser signal.
The sensing unit 3 comprises a sensing optical fiber 301, the first output end of the circulator 2 is connected with the sensing optical fiber 301, and the sensing optical fiber 301 is configured to transmit laser signals and interfere with vibration signals along the sensing optical fiber 301 to generate backward Rayleigh scattering signals.
The demodulation unit 4 comprises a second erbium doped fiber amplifier 401, a second filter 402, a detector 407, an acquisition card 408, a coupler 403, a first faraday rotator mirror 404 and a second faraday rotator mirror 405. The input end of the second erbium-doped fiber amplifier 401 is connected with the second output end of the circulator 2, the output end of the second erbium-doped fiber amplifier 401 is connected with the input end of the filter, the output end of the filter is connected with the first end of the coupler 403, the second end of the coupler 403 is connected with the input end of the detector 407, the output end of the detector 407 is connected with the input end of the acquisition card 408, the third end of the coupler 403 is connected with the first faraday rotator mirror 404, and the fourth end of the coupler 403 is connected with the second faraday rotator mirror 405.
The laser signal is input by the laser 101, sequentially passes through the isolator 102, the acousto-optic modulator 103, the erbium-doped fiber amplifier and the filter, and then is output to the sensing fiber 301 through the first output end of the circulator 2, the vibration signal along the sensing fiber 301 generates a backward rayleigh scattering signal within a preset time, and the backward rayleigh scattering signal is output to the demodulation unit 4 through the first input end and the second output end of the circulator 2 in sequence.
The backward Rayleigh scattering signal after entering the demodulation unit 4 enters the first Faraday rotator mirror 404 through the erbium-doped fiber amplifier and the filter, a part of the first end of the coupler 403 and the third end of the coupler 403 are reflected to the third end of the coupler 403 by the first Faraday rotator mirror 404, the other part of the backward Rayleigh scattering signal is output to the second Faraday rotator mirror 405 through the fourth end of the coupler 403, a second piezoelectric ceramic 406 is arranged at the second Faraday rotator mirror 405, a carrier signal is loaded through the second piezoelectric ceramic 406 to modulate the other part of the backward Rayleigh scattering signal, the modulated backward Rayleigh scattering signal is reflected to the fourth end of the coupler 403 by the second Faraday rotator mirror 405, the backward Rayleigh scattering signal of the fourth end of the coupler 403 and the backward Rayleigh scattering signal of the third end of the coupler 403 are all output to the second end of the coupler 403 and are mutually interfered and collected by the collecting card 408, the collecting card 408 is connected with the controller, the collecting card 408 outputs the collected signal to the controller, and the interference generated demodulation information is processed in the controller.
In some examples, referring to fig. 2 and 3, a method of measuring a distributed fiber optic acoustic wave sensing system includes the steps of:
s100: the backward Rayleigh scattering signal generated by the vibration signal along the line of the sensing optical fiber 301 in a preset time period is collected.
For example, in testing the distributed optical fiber acoustic wave sensing system, the distributed optical fiber acoustic wave sensing system may be set up, a vibration signal is applied to the first piezoelectric ceramic 302 on the sensing optical fiber 301 on the function generator 1031, and the distributed optical fiber acoustic wave sensing system is used to collect backward rayleigh scattering signals generated by the sensing optical fiber 301 along the line in a preset time period.
And S200, demodulating the collected backward Rayleigh scattering signals to obtain first demodulation information, wherein the first demodulation information comprises demodulation information containing noise information at each position along the sensing optical fiber 301.
The collected backward Rayleigh scattering signal is subjected to phase demodulation to obtain first demodulation information.
S300, constructing an original two-dimensional image through the first demodulation information, wherein the original two-dimensional image comprises an image of a time-distance domain formed by the first demodulation information corresponding to each position of the sensing optical fiber 301 along the line in a preset time period.
Illustratively, the original two-dimensional image includes two-dimensional graphic information of a time-distance domain or two-dimensional matrix information of a time-distance domain.
S400, performing two-dimensional wavelet transformation image processing on the original two-dimensional image to eliminate noise information and obtain a noise-removed two-dimensional image.
Illustratively, the noise-removed two-dimensional image includes two-dimensional graphic information of a time-distance domain or two-dimensional matrix information of a time-distance domain.
S410, carrying out wavelet decomposition on the original two-dimensional image through a discrete wavelet transform function to obtain a plurality of sub-images corresponding to the original two-dimensional image.
S420: the level of wavelet decomposition is adjusted to obtain a corresponding number of sub-images, and the high-frequency component and the low-frequency component of the sub-images are acquired and processed to eliminate the abnormal value on the sub-images so as to obtain the denoising sub-images.
S430: and carrying out wavelet reconstruction on the high-frequency component and the low-frequency component of the denoising sub-image by adopting an inverse discrete wavelet transform function to obtain a denoising two-dimensional image.
S500, reconstructing second demodulation information according to the denoising two-dimensional image, wherein the second demodulation information comprises demodulation information after the denoising information of each position along the sensing optical fiber 301.
And S510, extracting time domain information of the vibration signal in the second demodulation information.
Fig. 4 is a two-dimensional image corresponding to demodulation information of a vibration signal measured by a measuring method of a distributed optical fiber acoustic wave sensing system in the related art. Fig. 5 is a two-dimensional image corresponding to demodulation information of a vibration signal measured by a measuring method of a distributed optical fiber acoustic wave sensing system according to an embodiment of the present application. Fig. 6 is a time domain diagram of vibration positions corresponding to the two-dimensional image in fig. 4. Fig. 7 is a time domain diagram of vibration positions corresponding to the two-dimensional image in fig. 5.
For example, referring to fig. 4, the schematic diagram in fig. 4 is that the vibration signal applied to the sensing optical fiber 301 is measured by using the distributed optical fiber acoustic wave sensing system in the related art, and as can be seen from viewing fig. 3, there are multiple highlight interferences a in the area a and multiple interferences at the vibration signal analog diagram in the area B in fig. 3, so that the resolution of the vibration signal analog diagram in the area B is lower.
Referring to fig. 5, the schematic diagram in fig. 5 is that the vibration signal along the sensing optical fiber 301 is measured by using the measuring method of the distributed optical fiber acoustic wave sensing system in the embodiment of the present application, and it can be seen from observing fig. 5 that, compared with fig. 4, the highlight interference in the area A1 corresponding to the area a in fig. 5 is reduced and darkened, and similarly, the highlight interference in the area B1 corresponding to the area B in fig. 3 is also reduced in fig. 4, and the resolution of the vibration signal diagram in the area B1 is higher.
Similarly, referring to fig. 6, a vibration position time domain diagram corresponding to the two-dimensional image in fig. 4, a curve C is a vibration position time domain curve of the vibration signal, and referring to fig. 7, a vibration position time domain diagram corresponding to the two-dimensional image in fig. 5, a curve D is a vibration position time domain curve of the vibration signal. As can be seen from a comparison between fig. 6 and fig. 7, the vibration position time domain curve in fig. 7 has smaller fluctuation amplitude than the vibration position time domain curve in fig. 6.
As can be seen from comparison of fig. 4 and fig. 5, and comparison of fig. 6 and fig. 7, a measurement method of the distributed optical fiber acoustic wave sensing system is adopted to measure a certain vibration signal applied to the sensing optical fiber 301, so that interference of noise on the vibration signal is removed, and measurement accuracy is improved.
A distributed optical fiber acoustic wave sensing system comprises a demodulation unit 4 and a controller, wherein the demodulation unit 4 is connected with the controller through signals;
the demodulation unit 4 is configured to collect backward Rayleigh scattering signals within a preset time period, wherein the backward Rayleigh scattering signals comprise signals generated by a distributed optical fiber acoustic wave sensing system by vibration signals;
the demodulation unit 4 is further configured to demodulate the collected backward rayleigh scattering signals to obtain first demodulation information, where the first demodulation information includes demodulation information including noise information at each position along the sensing optical fiber 301;
the controller is configured to construct an original two-dimensional image by the first demodulation information, wherein the two-dimensional image comprises a time-distance domain image formed by the first demodulation information corresponding to each position of the sensing optical fiber 301 along the line in a preset time period;
the controller is further configured to perform two-dimensional wavelet transform image processing on the original two-dimensional image to eliminate noise information to obtain a noise-removed two-dimensional image;
the controller is further configured to reconstruct second demodulation information from the denoised two-dimensional image, the second demodulation information including the noise-cancelled demodulation information for each location along the sensing fiber 301.
The embodiment of the application also provides a distributed optical fiber acoustic wave sensing system, which adopts the measuring method in any one of the technical schemes, so that the distributed optical fiber acoustic wave sensing system has all the beneficial effects of the measuring method in any one of the technical schemes, and is not repeated here.

Claims (7)

1. A method for measuring a distributed optical fiber acoustic wave sensing system, comprising:
collecting backward Rayleigh scattering signals generated by vibration signals along the sensing optical fiber within a preset time period;
demodulating the collected backward Rayleigh scattering signals to obtain first demodulation information, wherein the first demodulation information comprises demodulation information of time-distance domains containing noise information of all positions along the sensing optical fiber;
constructing an original two-dimensional image of a time-distance domain through the first demodulation information, wherein the original two-dimensional image comprises an image of the time-distance domain formed by the first demodulation information corresponding to each position of the sensing optical fiber along the line in the preset time period, and the original two-dimensional image comprises two-dimensional graphic information of the time-distance domain or two-dimensional matrix information of the time-distance domain;
performing two-dimensional wavelet transformation image processing on the original two-dimensional image to eliminate noise information to obtain a noise-removed two-dimensional image, wherein the noise-removed two-dimensional image comprises two-dimensional graphic information of a time-distance domain or two-dimensional matrix information of a time-distance domain;
reconstructing second demodulation information according to the denoising two-dimensional image, wherein the second demodulation information comprises demodulation information after the denoising information of each position along the sensing optical fiber.
2. The method for measuring a distributed optical fiber acoustic wave sensor system according to claim 1, wherein said performing two-dimensional wavelet transform image processing on said original two-dimensional image comprises,
and carrying out wavelet decomposition on the original two-dimensional image through a discrete wavelet transform function to obtain a plurality of sub-images corresponding to the original two-dimensional image.
3. The method for measuring a distributed optical fiber acoustic wave sensor system according to claim 2, wherein said performing two-dimensional wavelet transform image processing on said original two-dimensional image further comprises,
and adjusting the wavelet decomposition level to obtain a corresponding number of sub-images, and acquiring and processing the high-frequency components and the low-frequency components of the sub-images to eliminate abnormal values on the sub-images and remove noise of the sub-images.
4. A method of measuring a distributed optical fiber acoustic wave sensing system according to claim 3 wherein the level of wavelet decomposition has an adjustment threshold within which the higher the level of wavelet decomposition, the less outliers the noise-canceling sub-image.
5. The method for measuring a distributed optical fiber acoustic wave sensor system according to claim 3, wherein the processing of the two-dimensional wavelet transform image on the original two-dimensional image further comprises performing wavelet reconstruction on the high-frequency component and the low-frequency component of the denoising sub-image by using an inverse discrete wavelet transform function to obtain the denoising two-dimensional image.
6. The method according to any one of claims 1 to 5, wherein after reconstructing the second demodulation information from the noise-removed two-dimensional image, extracting time domain information of the vibration signal in the second demodulation information.
7. The distributed optical fiber acoustic wave sensing system is characterized by comprising a demodulation unit and a controller, wherein the demodulation unit is connected with the controller through signals;
the demodulation unit is configured to collect backward Rayleigh scattering signals within a preset time period, wherein the backward Rayleigh scattering signals comprise signals generated by vibration signals in the distributed optical fiber acoustic wave sensing system;
the demodulation unit is further configured to demodulate the collected backward Rayleigh scattering signals to obtain first demodulation information, wherein the first demodulation information comprises demodulation information of time-distance domains containing noise information of all positions along the sensing optical fiber;
the controller is configured to construct an original two-dimensional image of a time-distance domain through the first demodulation information, wherein the original two-dimensional image comprises a time-distance domain image formed by the first demodulation information corresponding to each position of the sensing optical fiber along the line in the preset time period, and the original two-dimensional image comprises two-dimensional graphic information of the time-distance domain or two-dimensional matrix information of the time-distance domain;
the controller is further configured to perform two-dimensional wavelet transform image processing on the original two-dimensional image to eliminate noise information to obtain a noise-removed two-dimensional image, wherein the noise-removed two-dimensional image comprises two-dimensional graphic information of a time-distance domain or two-dimensional matrix information of a time-distance domain;
the controller is further configured to reconstruct second demodulation information from the denoised two-dimensional image, the second demodulation information including noise information removed demodulation information for each location along the sensing fiber.
CN202310664799.3A 2023-06-07 2023-06-07 Distributed optical fiber acoustic wave sensing system and measuring method thereof Active CN116399379B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310664799.3A CN116399379B (en) 2023-06-07 2023-06-07 Distributed optical fiber acoustic wave sensing system and measuring method thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310664799.3A CN116399379B (en) 2023-06-07 2023-06-07 Distributed optical fiber acoustic wave sensing system and measuring method thereof

Publications (2)

Publication Number Publication Date
CN116399379A CN116399379A (en) 2023-07-07
CN116399379B true CN116399379B (en) 2023-11-03

Family

ID=87018294

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310664799.3A Active CN116399379B (en) 2023-06-07 2023-06-07 Distributed optical fiber acoustic wave sensing system and measuring method thereof

Country Status (1)

Country Link
CN (1) CN116399379B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116907677B (en) * 2023-09-15 2023-11-21 山东省科学院激光研究所 Distributed optical fiber temperature sensing system for concrete structure and measuring method thereof

Citations (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5844415A (en) * 1994-02-03 1998-12-01 Massachusetts Institute Of Technology Method for three-dimensional positions, orientation and mass distribution
KR20050032134A (en) * 2003-10-01 2005-04-07 한국표준과학연구원 An intrusion sensing device using optical time domain reflectometry detecting structure
CN103411660A (en) * 2013-08-29 2013-11-27 山东省科学院激光研究所 Optical fiber distributed type sound wave monitor system
CN105092014A (en) * 2015-05-12 2015-11-25 天津大学 Distributed fiber sound wave detection apparatus and method based on wave beam formation
CN106441386A (en) * 2016-09-29 2017-02-22 西南交通大学 Data processing method and device based on distributed fiber sensing system
CN107270952A (en) * 2017-07-27 2017-10-20 天津求实飞博科技有限公司 Based on the distributed destabilization sensing signal processing method of probe beam deflation long-distance optical fiber
CN107436158A (en) * 2017-07-27 2017-12-05 天津求实飞博科技有限公司 A kind of Demodulation System for Fiber Bragg
CN108225385A (en) * 2017-11-23 2018-06-29 上海电力学院 A kind of overlapping FBG transducing signal peak value localization methods
CN108253999A (en) * 2016-12-29 2018-07-06 中国科学院半导体研究所 A kind of noise-reduction method for distribution type fiber-optic sound sensor-based system
CN108732614A (en) * 2018-05-24 2018-11-02 华中科技大学 A kind of online railway monitoring system and method based on distributed acoustic sensor
CN108875684A (en) * 2018-06-29 2018-11-23 电子科技大学 Train operation state method for parameter estimation based on optical fiber sensing time-frequency figure processing
CN109919929A (en) * 2019-03-06 2019-06-21 电子科技大学 A kind of fissuring of tongue feature extracting method based on wavelet transformation
CN110160569A (en) * 2019-04-24 2019-08-23 国网浙江省电力有限公司信息通信分公司 For the noise-reduction method of distributing optical fiber sensing signal, system and storage medium
CN110608760A (en) * 2019-07-05 2019-12-24 中国地质大学(武汉) Method for improving disturbance detection signal-to-noise ratio of phase sensitive optical time domain reflection sensing system
CN110751073A (en) * 2019-10-12 2020-02-04 武汉理工大学 Pipeline early damage mode identification method based on distributed optical fiber sensing and deep learning
CN110768715A (en) * 2019-12-26 2020-02-07 华中科技大学 Polarized light time domain reflectometer based on time division multiplexing of three polarization states and detection method
CN112033522A (en) * 2020-08-10 2020-12-04 太原理工大学 Sound signal detection system and method of distributed optical fiber sensor
WO2021004076A1 (en) * 2019-07-05 2021-01-14 山东大学 Ai chip-based conformal wearable biological information monitoring device and system
CN113188461A (en) * 2021-05-06 2021-07-30 山东大学 OFDR large strain measurement method under high spatial resolution
CN113447110A (en) * 2021-06-10 2021-09-28 天津大学 Distributed optical fiber vibration sensing system and phase carrier demodulation method thereof
CN113499059A (en) * 2021-06-01 2021-10-15 武汉理工大学 BCG signal processing system and method based on optical fiber sensing non-contact
CN113553997A (en) * 2021-09-17 2021-10-26 南京信息工程大学 Signal processing method, storage medium and system for jointly improving wavelet threshold
CN114353844A (en) * 2021-12-03 2022-04-15 湖北民族大学 Fiber grating signal demodulation method
CN114461975A (en) * 2022-01-27 2022-05-10 无锡科晟光子科技有限公司 Optical fiber sound measuring method based on wavelet transformation denoising
CN114623921A (en) * 2022-03-17 2022-06-14 中国科学技术大学 Demodulation method of distributed optical fiber vibration sensing
CN115096427A (en) * 2022-06-09 2022-09-23 电子科技大学 Distributed optical fiber array sensing sound wave detection system and method
CN115265750A (en) * 2022-07-19 2022-11-01 温州市质量技术检测科学研究院 Optical fiber distributed acoustic wave sensing system and method
CN115307715A (en) * 2022-08-09 2022-11-08 南昌航空大学 Improved wavelet denoising method based on sagnac optical fiber acoustic sensing system

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6859275B2 (en) * 1999-04-09 2005-02-22 Plain Sight Systems, Inc. System and method for encoded spatio-spectral information processing
US7260272B2 (en) * 2003-07-10 2007-08-21 Samsung Electronics Co.. Ltd. Method and apparatus for noise reduction using discrete wavelet transform
JP5948035B2 (en) * 2011-10-05 2016-07-06 ニューブレクス株式会社 Distributed optical fiber acoustic wave detector
GB201408125D0 (en) * 2014-05-08 2014-06-25 Optasense Holdings Ltd Fibre optic distributed sensing
WO2016010550A1 (en) * 2014-07-18 2016-01-21 Halliburton Energy Services, Inc. Distributed sensing systems and methods with efficient energy spectrum analysis
US20180045542A1 (en) * 2015-06-22 2018-02-15 Omnisens Sa A method for reducing noise in measurements taken by a distributed sensor
US10018714B2 (en) * 2015-06-24 2018-07-10 The United States Of America As Represented By The Secretary Of The Army Two-dimensional RF harmonic imaging system and algorithm
WO2019148024A1 (en) * 2018-01-26 2019-08-01 Park Jong Kang Systems and methods to reduce scattering in temporal focusing multiphoton microscopy
CN108981767B (en) * 2018-08-17 2020-10-16 山东省科学院激光研究所 Optical fiber distributed monitoring system and method
US11468667B2 (en) * 2019-06-19 2022-10-11 Nec Corporation Distributed intelligent traffic informatics using fiber sensing
FR3099245B1 (en) * 2019-07-26 2021-07-30 Febus Optics METHOD AND DEVICE FOR RECONSTRUCTING A RETRODUCED VECTOR ELECTROMAGNETIC WAVE
US20210356332A1 (en) * 2020-05-18 2021-11-18 Nec Laboratories America, Inc Joint wavelet denoising for distributed temperature sensing

Patent Citations (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5844415A (en) * 1994-02-03 1998-12-01 Massachusetts Institute Of Technology Method for three-dimensional positions, orientation and mass distribution
KR20050032134A (en) * 2003-10-01 2005-04-07 한국표준과학연구원 An intrusion sensing device using optical time domain reflectometry detecting structure
CN103411660A (en) * 2013-08-29 2013-11-27 山东省科学院激光研究所 Optical fiber distributed type sound wave monitor system
CN105092014A (en) * 2015-05-12 2015-11-25 天津大学 Distributed fiber sound wave detection apparatus and method based on wave beam formation
CN106441386A (en) * 2016-09-29 2017-02-22 西南交通大学 Data processing method and device based on distributed fiber sensing system
CN108253999A (en) * 2016-12-29 2018-07-06 中国科学院半导体研究所 A kind of noise-reduction method for distribution type fiber-optic sound sensor-based system
CN107436158A (en) * 2017-07-27 2017-12-05 天津求实飞博科技有限公司 A kind of Demodulation System for Fiber Bragg
CN107270952A (en) * 2017-07-27 2017-10-20 天津求实飞博科技有限公司 Based on the distributed destabilization sensing signal processing method of probe beam deflation long-distance optical fiber
CN108225385A (en) * 2017-11-23 2018-06-29 上海电力学院 A kind of overlapping FBG transducing signal peak value localization methods
CN108732614A (en) * 2018-05-24 2018-11-02 华中科技大学 A kind of online railway monitoring system and method based on distributed acoustic sensor
CN108875684A (en) * 2018-06-29 2018-11-23 电子科技大学 Train operation state method for parameter estimation based on optical fiber sensing time-frequency figure processing
CN109919929A (en) * 2019-03-06 2019-06-21 电子科技大学 A kind of fissuring of tongue feature extracting method based on wavelet transformation
CN110160569A (en) * 2019-04-24 2019-08-23 国网浙江省电力有限公司信息通信分公司 For the noise-reduction method of distributing optical fiber sensing signal, system and storage medium
WO2021004076A1 (en) * 2019-07-05 2021-01-14 山东大学 Ai chip-based conformal wearable biological information monitoring device and system
CN110608760A (en) * 2019-07-05 2019-12-24 中国地质大学(武汉) Method for improving disturbance detection signal-to-noise ratio of phase sensitive optical time domain reflection sensing system
CN110751073A (en) * 2019-10-12 2020-02-04 武汉理工大学 Pipeline early damage mode identification method based on distributed optical fiber sensing and deep learning
CN110768715A (en) * 2019-12-26 2020-02-07 华中科技大学 Polarized light time domain reflectometer based on time division multiplexing of three polarization states and detection method
CN112033522A (en) * 2020-08-10 2020-12-04 太原理工大学 Sound signal detection system and method of distributed optical fiber sensor
CN113188461A (en) * 2021-05-06 2021-07-30 山东大学 OFDR large strain measurement method under high spatial resolution
CN113499059A (en) * 2021-06-01 2021-10-15 武汉理工大学 BCG signal processing system and method based on optical fiber sensing non-contact
CN113447110A (en) * 2021-06-10 2021-09-28 天津大学 Distributed optical fiber vibration sensing system and phase carrier demodulation method thereof
CN113553997A (en) * 2021-09-17 2021-10-26 南京信息工程大学 Signal processing method, storage medium and system for jointly improving wavelet threshold
CN114353844A (en) * 2021-12-03 2022-04-15 湖北民族大学 Fiber grating signal demodulation method
CN114461975A (en) * 2022-01-27 2022-05-10 无锡科晟光子科技有限公司 Optical fiber sound measuring method based on wavelet transformation denoising
CN114623921A (en) * 2022-03-17 2022-06-14 中国科学技术大学 Demodulation method of distributed optical fiber vibration sensing
CN115096427A (en) * 2022-06-09 2022-09-23 电子科技大学 Distributed optical fiber array sensing sound wave detection system and method
CN115265750A (en) * 2022-07-19 2022-11-01 温州市质量技术检测科学研究院 Optical fiber distributed acoustic wave sensing system and method
CN115307715A (en) * 2022-08-09 2022-11-08 南昌航空大学 Improved wavelet denoising method based on sagnac optical fiber acoustic sensing system

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
一种空域与小波域相结合的水印算法;季亚新;李鸿林;;应用科技(第03期);全文 *
光纤分布式振动系统中的信号处理算法研究;文仲寺;江毅;;仪表技术与传感器(第06期);全文 *
分布式光纤声波传感技术研究进展;蔡海文;叶青;王照勇;卢斌;曹玉龙;;应用科学学报(第01期);全文 *
基于均值滤波和小波变换的图像去噪;彭姝姝;;现代计算机(第12期);全文 *
基于实时小波消噪的光纤光栅解调系统研究;史振江;李志全;郑华;李靖;;压电与声光(第03期);全文 *
基于灰度图像的光纤相干光时域反射计信号处理研究;赵世琦;庞拂飞;贺梦婷;徐慧;王廷云;;中国激光(第03期);全文 *
基于自适应噪声完备经验模态分解的Φ-OTDR信号去噪算法;熊兴隆;魏永兴;张琬童;李猛;马愈昭;;半导体光电(第04期);全文 *
星载雷达弱目标小波去噪方法研究;侯颖辉;许荣庆;;黑龙江大学自然科学学报(第05期);全文 *

Also Published As

Publication number Publication date
CN116399379A (en) 2023-07-07

Similar Documents

Publication Publication Date Title
JP6342019B2 (en) Distributed optical fiber acoustic wave detector
CN116399379B (en) Distributed optical fiber acoustic wave sensing system and measuring method thereof
EP3364147B1 (en) Method and apparatus for motion compensation in interferometric sensing systems
CN107907151B (en) Phi OTDR and sensing method based on optical frequency modulation and direct detection
CN110132329B (en) Stress, temperature and vibration composite detection optical fiber sensor and signal processing method
CN106840222B (en) A kind of distributed optical fiber sensing system and its suppressing method inhibiting common-mode noise
CN110864797A (en) Differential COTDR distributed acoustic sensing device and method for heterogeneous double-sideband chirped pulses
CN111157097B (en) Vibration signal detection method based on phase-sensitive optical time domain reflection system
Tu et al. Enhancement of signal identification and extraction in a Φ-OTDR vibration sensor
CN116380140B (en) Distributed acoustic wave sensing system based on mean value filtering technology and measuring method thereof
Cunzheng et al. Large dynamic range optical fiber distributed acoustic sensing (DAS) with differential-unwrapping-integral algorithm
CN109391321A (en) Disturbance positioning method in a kind of phase sensitive OTDR sensing
CN108827447B (en) Different-frequency double-pulse COTDR sensing device and method
Qu et al. Phase sensitive optical time domain reflectometry based on compressive sensing
Chen et al. Denoising method of the Φ-OTDR system based on EMD-PCC
Chen et al. Denoising of distributed acoustic sensing seismic data using an integrated framework
Wu et al. NLM Parameter Optimization for $\varphi $-OTDR Signal
Bai et al. Noise reduction method of Φ-OTDR system based on EMD-TFPF algorithm
Mao et al. Denoising method based on VMD-PCC in φ-OTDR system
CN116818085B (en) Distributed optical fiber acoustic wave sensing system based on weak grating and measuring method thereof
Li et al. Denoising method based upon ECEMD and selective addition for phi-OTDR
CN115333913B (en) Phase demodulation method and device for phase sensitive optical time domain reflectometer
Li et al. Denoising algorithm of Φ-OTDR signal based on curvelet transform with adaptive threshold
Boutkedjirt et al. Reconstruction of ultrasonic fields by deconvolving the hydrophone aperture effects: I. Theory and simulation
Sun et al. Noise Suppression of Distributed Acoustic Sensing Based on f-x Deconvolution and Wavelet Transform

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