DISTRIBUTED FIBER-OPTIC ACOUSTIC SENSING SYSTEM AND SIGNAL PROCESSING METHOD USING THE SAME

Abstract

A distributed fiber-optic acoustic sensing system and a signal processing method. The distributed fiber-optic acoustic sensing system is based on a high spatial resolution distributed fiber-optic acoustic sensor. The interval between adjacent sensing units is centimeter or millimeter level. Through specific digital signal processing, signal enhancement can be realized, noise in the system and environment are suppressed, at the same time, problems such as interference fading is solved, and the sensor signal-to-noise ratio of subunits can be increased by two to three orders of magnitude. Each subunit can serve as an independent high-sensitivity sensor for sensing. The multiple subunits can form one or more new sensor arrays. The azimuth estimation and spatial orientation of signal sources can be realized by the array signal processing method.

Claims

1. A distributed fiber-optic acoustic sensing system, comprising a distributed fiber-optic acoustic sensor comprising a plurality of sensing units along an axial direction, and a sensing optical fiber, wherein the distributed fiber-optic acoustic sensor is connected to the sensing optical fiber and is used for emitting a detection light pulse to the sensing optical fiber, collecting and demodulating sound fields of signal sources sensed along the sensing optical fiber; a total number of the sensing units in the axial direction of the distributed fiber-optic acoustic sensor is L, the sensing units are divided into M subunits in the axial direction, and the number of the sensing units in each subunit is N, L=M×N, and the sensing units in each subunit are used for spatially coherent combination of sound field signals obtained by detection; and spatial resolution of the distributed fiber-optic acoustic sensor is centimeter-level or higher.

2. The distributed fiber-optic acoustic sensing system according to claim 1, wherein the distributed fiber-optic acoustic sensor comprises a distributed fiber-optic sensor based on a sweep frequency pulse and a pulse compression technology, an optical frequency domain reflectometer, and one or more of phase-sensitive optical time-domain reflectometers based on an ultra-narrow pulse or coherent detection.

3. The distributed fiber-optic acoustic sensing system according to claim 1, wherein the signal sources comprise one or more of objects that vibrate in the air, soil, or water to produce sound waves.

4. The distributed fiber-optic acoustic sensing system according to claim 1, wherein the number M of the subunits 10≤M≤10,000.

5. The distributed fiber-optic acoustic sensing system according to claim 1, wherein the number N of the sensing units in each subunit is 10≤N≤10,000.

6. The distributed fiber-optic acoustic sensing system according to claim 1, wherein the spatial scale of the subunits is meter or sub-meter level.

7. A signal processing method using the distributed fiber-optic acoustic sensing system as described in claim 1, comprising step (1) collecting distributed sound field signal, wherein the distributed sound field signals collected by the distributed fiber-optic acoustic sensor are S(l,t), l represents one-dimensional axial space of the sensing optical fiber and t represents time; step (2) dividing the subunits, wherein the total number of the sensing units in the axial direction of the distributed fiber-optic acoustic sensor is L, the sensing units are divided into M subunits in the axial direction, and the number of the sensing units in each subunit is N, and L=M×N is met, the sound field signals collected by the N sensing units in the mth subunit are: [S.sub.m(l.sub.0+(m−1).Math.N.Math.Δz,t), S.sub.m(l.sub.0+(m−1).Math.N.Math.Δz+Δz,t), K S.sub.m(l.sub.0+(m−1).Math.N.Math.Δz+(N−1)Δz,t)] that is abbreviated as [S.sub.m,1, S.sub.m,2, K S.sub.m,N]; m meets 1≤m≤M; l.sub.0 represents a starting position of the optical fiber; Δz represents the distance between the two adjacent sensing units; and the minimum value of Δz is determined by the spatial resolution of the distributed fiber-optic acoustic sensor; step (3) conducting spatially coherent combination in the subunits and realizing signal-to-noise ratio and sensitivity enhancement, wherein the sound field signals [S.sub.m(l.sub.0+(m−1).Math.N.Math.Δz,t), S.sub.m(l.sub.0+(m−1).Math.N.Math.Δz+Δz,t), K S.sub.m(l.sub.0+(m−1).Math.N.Math.Δz+(N−1)Δz,t)] collected by the N sensing units in the mth subunit are enhanced by a spatially coherent combination algorithm to generate a signal S′.sub.m(l.sub.0+(m−1).Math.N.Math.Δz+Δz, t) that is abbreviated as S′.sub.m; and after the M subunits are enhanced by array signal processing, obtained new sound field signals are [S′.sub.1, S′.sub.2, K S′.sub.m, K S′.sub.M]; and step (4) realizing azimuth estimation and space orientation by the subunits serving as independent high-sensitivity signal-to-noise ratio sensors or a new array formed by the subunits, wherein each subunit can serve as the independent high-sensitivity sensor for sensing, and at the same time, M.sub.x subunits which are adjacent or at regular intervals can also be selected from the M subunits to serve as the new sensor array; the M subunits can form one or more sensing arrays; the newly formed sensor array performs the azimuth estimation or space orientation of the signal sources by a directional array signal processing method; and the process similar to the step (3) can be repeated again to achieve signal enhancement.

8. The signal processing method according to claim 7, wherein the spatially coherent combination in the step (3) comprises a vector superposition or beamforming array signal processing method.

9. The signal processing method according to claim 7, wherein the array signal processing method performing azimuth estimation or space orientation in the step (4) comprises one of a conventional beamforming algorithm, a spatial spectrum estimation algorithm, or a direction-of-arrival estimation algorithm.

10. The signal processing method according to claim 7, wherein the array signal processing method implemented in the step (3) is a self-adaptive spatial filtering method, a delay-and-sum method, or both.

11. The signal processing method according to claim 10, wherein the self-adaptive spatial filtering method is a method adopting a minimum variance distortionless response beamformer (MVDR), a method adopting a linearly constrained minimum variance beamformer (LCMV), or a method adopting a generalized sidelobe cancellation beamformer (GSC).

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0066] FIG. 1 shows the basic structure of the distributed fiber-optic acoustic sensing system of the present invention.

[0067] FIGS. 2A to 2D show the signal processing method of the present invention, where FIG. 2A shows distributed sound field signal acquisition; FIG. 2B shows division of subunits; FIG. 2C shows spatial coherent combination in subunits and realization of signal-to-noise ratio and sensitivity enhancement; and FIG. 2D shows realization of azimuth estimation and spatial orientation by subunits serving as independent high-sensitivity and high-signal-to-noise ratio sensors or a new array formed by the subunits.

[0068] FIG. 3 shows the distributed fiber-optic sensor based on a sweep frequency pulse and a pulse compression technology of the present invention.

[0069] FIG. 4 shows the optical frequency domain reflectometer with a Michelson interferometer structure in the present invention.

[0070] FIG. 5 shows the phase-sensitive optical time domain reflectometer based on an ultra-narrow pulse and coherent detection in the present invention

DETAILED DESCRIPTION OF THE INVENTION

[0071] The present invention is described in detail in connection with the figures and the following embodiments, however, they do not serve to limit the scope of the present invention.

[0072] One of skilled in the art may modify or change the present invention as shown in the embodiments without departing from the scope of the present invention.

Embodiment I

[0073] In the first embodiment of the present invention as shown in FIG. 1, the distributed fiber-optic acoustic sensing system comprises a high spatial resolution distributed fiber-optic acoustic sensor 1 and a sensing optical fiber 2. The total number of sensing units in the axial direction of the distributed fiber-optic acoustic sensor is L. The sensing units are divided into M subunits in the axial direction, and the number of the sensing units in each subunit is N, and the condition that L=M×N is met. The distributed fiber-optic acoustic sensor is connected to the sensing optical fiber and is used for emitting a detection light pulse to the sensing optical fiber, collecting and demodulating the sound fields of signal sources 3 sensed along the sensing optical fiber. The distributed fiber-optic acoustic sensor is based on a sweep frequency pulse and a pulse compression technology.

[0074] As shown in FIG. 3, the range of the sweep frequency pulse is 10 GHz, the spatial resolution is 1 cm, and the sensing optical fiber is a SM-28 ordinary single-mode optical fiber/cable. The distributed fiber-optic acoustic sensor based on the sweep frequency pulse and the pulse compression technology comprises a single-frequency narrow linewidth laser 6-1, a first fiber-optic coupler 6-2, a linear frequency-modulation pulse modulation module 6-3, a circulator 6-4, a sensing optical fiber interface 101, a second fiber-optic coupler 6-5, a partial polarization coherent detection unit 6-6, a data acquisition unit 6-7, and a digital signal processing unit 6-8. Single-frequency narrow linewidth laser light output by the single-frequency narrow linewidth laser 6-1 is output to a second port 202 and a third port 203 of the first fiber-optic coupler 6-2 through a first port 201 of the first fiber-optic coupler 6-2. The second port 202 of the first fiber-optic coupler 6-2 is connected to the linear frequency-modulation pulse modulation module 6-3; the linear frequency-modulation pulse modulation module 6-3 is used for modulating continuous light into a linear frequency modulation light pulse and is connected to a first port 401 of the circulator 6-4. A second port 402 of the circulator 6-4 is connected to the sensing optical fiber interface 101. A third port 403 of the circulator 6-4 and the third port 203 of the first fiber-optic coupler 6-2 are respectively connected to a first port 501 and a third port 503 of the second fiber-optic coupler 6-5. A second port 502 and a fourth port 504 of the second fiber-optic coupler 6-5 are connected to an input port of the partial polarization coherent detection unit 6-6. An output port of the partial polarization coherent detection unit 6-6 is connected to an input port of the data acquisition unit 6-7. The data acquisition unit 6-7 is connected to the digital signal processing unit 6-8, and the digital signal processing unit 6-8 is used for realizing pulse compression processing, signal demodulation etc.

[0075] The spatially coherent combination method in the present invention adopts a delay-and-sum beamforming algorithm. After the subunits are enhanced, the subunits serve as high-sensitivity distributed fiber-optic acoustic sensors and do not perform azimuth estimation and space orientation of signal sources; and the implementation steps are as follows as shown in FIGS. 2A to 2D:

[0076] step (1): as shown in FIG. 2A, the high spatial resolution distributed fiber-optic acoustic sensor is connected to the sensing optical fiber and is used for emitting the detection light pulse to the sensing optical fiber, collecting and quantitatively demodulating the sound fields of the signal sources sensed along the sensing optical fiber to obtain distributed sound field signals S(l,t), wherein l represents one-dimensional axial space of the sensing optical fiber, and t represents time;

[0077] step (2): as shown in FIG. 2B, the total number of the sensing units in the axial direction of the high spatial resolution distributed fiber-optic acoustic sensor is L, the sensing units are divided into M subunits in the axial direction, the number of the sensing units in each subunit is N, and L=M×N is met, wherein the sound field signals collected by the N sensing units in the mth subunit are: [S.sub.m(l.sub.0+(m−1).Math.N.Math.Δz,t), S.sub.m (l.sub.0+(m−1).Math.N.Math.Δz+Δz,t), K S.sub.m(l.sub.0+(m−1).Math.N.Math.Δz+(N−1)Δz,t)], which is abbreviated as [S.sub.m,1, S.sub.m,2, K S.sub.m,N]; m is greater than or equal to 1 and less than or equal to M; l.sub.0 represents the starting position of the optical fiber; Δz represents the distance between the two adjacent sensing units; the minimum value of Δz is determined by the spatial resolution of the distributed fiber-optic acoustic sensor; the spatial resolution of the distributed fiber-optic acoustic sensor is centimeter level or higher; N is about tens or hundreds; and the spatial scale of each subunit is meter or submeter level;

[0078] step (3): calculating the time delay vector of the signal source with respect to each sensing unit in the subunit according to the spatial position information of each sensing unit in each subunit, for example, the time delay vector of the signal source with respect to each sensing unit in the mth subunit is a.sub.p (Δτ)=[Δτ.sub.m,1Δτ.sub.m,2 L Δτ.sub.m,N], or the phase retardation is a.sub.p (Δφ)=[Δφ.sub.m,1 Δφ.sub.m,2 L Δφ.sub.m,N], wherein p=1, 2, L, P, P represents the number of signal sources;

[0079] step (4): calculating delay compensation weight W.sub.p(Δτ)=[−Δτ.sub.m,1 −Δτ.sub.m,2 L −Δτ.sub.m,N] according to the time delay vector a.sub.p (Δτ)=[Δτ.sub.m,1 Δτ.sub.m,2 L Δτ.sub.m,N] of the signal source with respect to each sensing unit in the mth subunit;

[0080] obtaining an enhanced signal of the pth signal source by calculating according to the following mode: S′.sub.m=W.sub.p.Math.[S.sub.m,1, S.sub.m,2, K S.sub.m,N].sup.T, wherein T represents transposition operation, other P−1 signal sources and environmental noise can be suppressed, and the signal suppression effect of the P−1 signal sources is limited by the number of the sensing units in the entire subunit. After the sound field signals collected by the N sensing units in the mth subunit are enhanced by the above-mentioned array signal processing, the environmental noise can be effectively suppressed, the noise floor of the sensing system is reduced, the signal-to-noise ratio and sensitivity of the sensing system are improved, and at the same time, the interference fading and other problems of the fiber-optic sensing system can be solved. Wherein, the sound field signals [S.sub.m(l.sub.0+(m−1).Math.N.Math.Δz,t), S.sub.m(l.sub.0+(m−1).Math.N.Math.Δz+Δz,t), K S.sub.m(l.sub.0+(m−1).Math.N.Math.Δz+(N−1)Δz,t)] collected by the N sensing units in the mth subunit are enhanced by array signal processing to generate a signal S′.sub.m (l.sub.0+(m−1).Math.N.Math.Δz+Δz,t), which is abbreviated as S′.sub.m; and after the M subunits are enhanced by array signal processing, obtained new sound field signals are [S′.sub.1, S′.sub.2, K S′.sub.m, K S′.sub.M], which is shown in FIG. 2C; and

[0081] step (5): as shown in FIG. 2D, each subunit can serve as an independent high-sensitivity sensor for sensing, and the process similar to the steps (3) and (4) can also be repeated again as needed to further realize signal enhancement.

Embodiment II

[0082] In the second embodiment of the present invention, the distributed fiber-optic acoustic sensor adopts an optical frequency domain reflectometer with a Michelson interferometer structure. As shown in FIG. 4, the sweep frequency range of a linear sweep frequency laser is 100 GHz, the spatial resolution of the sensor is about 1 mm, and the sensing optical fiber is a SM-28 ordinary single-mode optical fiber/cable. The optical frequency domain reflectometer is of the Michelson interferometer structure and comprises a frequency tunable laser 7-1, a fiber-optic coupler 7-2, an optical fiber reflecting end face 7-3, a photodetector 7-4, a data acquisition unit 7-5, a digital signal processing unit 7-6, and a sensing optical fiber interface 101. The frequency tunable laser 7-1 is used for outputting linear frequency modulation laser light which is output to a second port 202 and a fourth port 204 of the fiber-optic coupler 7-2 through a first port 201 of the fiber-optic coupler 7-2; the fourth port 204 of the fiber-optic coupler 7-2 is connected to the optical fiber reflecting end face 7-3; the second port 202 of the fiber-optic coupler 7-2 is connected to the sensing optical fiber interface 101; a third port 203 of the fiber-optic coupler 7-2 is connected to the photodetector 7-4; a signal output by the photodetector 7-4 is connected to the data acquisition unit 7-5; and the data acquisition unit 7-5 is connected to the digital signal processing unit 7-6.

[0083] The spatially coherent combination method of the present invention adopts the delay-and-sum beamforming algorithm; after the subunits are enhanced, a new array is formed by the subunits and is used for performing the two-dimensional spatial spectrum estimation of the signal sources. The implementation steps are as follows:

[0084] step (1): the high spatial resolution distributed fiber-optic acoustic sensor is connected to the sensing optical fiber, and is used for emitting the detection light pulse to the sensing optical fiber, collecting and quantitatively demodulating the sound fields of the signal sources sensed along the sensing optical fiber to obtain distributed sound field signals S(l,t), wherein l represents one-dimensional axial space of the sensing optical fiber, and t represents time;

[0085] step (2): the total number of the sensing units in the axial direction of the high spatial resolution distributed fiber-optic acoustic sensor is L, the sensing units are divided into M subunits in the axial direction, the number of the sensing units in each subunit is N, and L=M×N is met, wherein the sound field signals collected by the N sensing units in the mth subunit are:

[0086] [S.sub.m (l.sub.0+(m−1).Math.N.Math.Δz,t), S.sub.m(l.sub.0+(m−1).Math.N.Math.Δz+Δz,t), K S.sub.m (l.sub.0+(m−1).Math.N.Math.Δz+(N−1)Δz,t)], which is abbreviated as [S.sub.m,1, S.sub.m,2, K S.sub.m,N]; m is greater than or equal to 1 and less than or equal to M; l.sub.0 represents the starting position of the optical fiber; Δz represents the distance between the two adjacent sensing units; the minimum value of Δz is determined by the spatial resolution of the distributed fiber-optic acoustic sensor; the spatial resolution of the distributed fiber-optic acoustic sensor is centimeter level or higher; N is about tens or hundreds; and the spatial scale of each subunit is meter or submeter level;

[0087] step (3): calculating the time delay vector of the signal source with respect to each sensing unit in the subunit according to the spatial position information of each sensing unit in each subunit, for example, the time delay vector of the signal source with respect to each sensing unit in the mth subunit is a.sub.p (Δτ)=[Δτ.sub.m,1 Δτ.sub.m,2 L Δτ.sub.m,N], or the phase retardation is a.sub.p (Δφ)=[Δφ.sub.m,1Δφ.sub.m,2 L Δφ.sub.m,N], wherein p=1, 2, L, P, P represents the number of signal source;

[0088] step (4): calculating delay compensation weight W.sub.p(Δτ)=[−Δτ.sub.m,1 −Δτ.sub.m,2 L −Δτ.sub.m,N] according to the time delay vector a.sub.p (Δτ)=[Δτ.sub.m,1 Δτ.sub.m,2 L Δτ.sub.m,N] of the signal source with respect to each sensing unit in the mth subunit;

[0089] obtaining the enhanced signal of the pth signal source by calculating according to the following mode: S′.sub.m=W.sub.p.Math.[S.sub.m,1, S.sub.m,2, K S.sub.m,N].sup.T, wherein T represents transposition operation, other P−1 signal sources and environmental noise can be suppressed, and the signal suppression effect of the P−1 signal sources is limited by the number of the sensing units in the entire subunit. After the sound field signals collected by the N sensing units in the mth subunit are enhanced by the above-mentioned array signal processing, the environmental noise can be effectively suppressed, the noise floor of the sensing system is reduced, the signal-to-noise ratio and sensitivity of the sensing system are improved, and at the same time, the interference fading and other problems of the fiber-optic sensing system can be solved. Wherein the sound field signals [S.sub.m(l.sub.0+(m−1).Math.N.Math.Δz,t), S.sub.m(l.sub.0+(m−1).Math.N.Math.Δz+Δz,t), K S.sub.m(l.sub.0+(m−1).Math.N.Math.Δz+(N−1)Δz,t)] collected by the N sensing units in the mth subunit are enhanced by array signal processing to generate a signal S′.sub.m (l.sub.0+(m−1).Math.N.Math.Δz+Δz, t), which is abbreviated as S′.sub.m; and after the M subunits are enhanced by array signal processing, the obtained new sound field signals are [S′.sub.1, S′.sub.2, K S′.sub.m, K S′.sub.M]; and

[0090] step (5): the main steps of two-dimensional space spectrum estimation are as follows:

[0091] firstly, selecting new sound field signals [S′.sub.m, S′.sub.m+1, K S′.sub.m+M.sub.x.sub.−1] after the M.sub.x subunits are enhanced from the new sound field signals generated in the step (4), wherein the new sound field signals are distributed along the optical fiber;

[0092] secondly, setting up a coordinate system (x, y, z) in the three-dimensional space, and performing coordinate system conversion on [S′.sub.m, S′.sub.m+1, K S′.sub.m+M.sub.x.sub.−l] according to the specific sensing optical fiber laying method T [x, y, z], wherein S′.sub.m(l.sub.0+(m−1).Math.N.Math.Δz+Δz, t) is converted into S″.sub.m(x.sub.m, y.sub.m, z.sub.m, t), which is abbreviated as S″.sub.m, that is, [S′.sub.m, S′.sub.m+1, K S′.sub.m+M.sub.x.sub.−1] is converted into [S″.sub.m, S″.sub.m+1, K S″.sub.m+M.sub.x.sub.−1] after being subjected to coordinate conversion;

[0093] thirdly, calculating the covariance matrix of the above sound field signals;

[00007]$R^{″}=\underset{i=1}{\overset{K}{.Math.}}{\left[{S^{″}}_{m+1},{S^{″}}_{m+2},K{S^{″}}_{m+m_x}\right]}^T.Math.\left[{S^{″}}_{m+1},{S^{″}}_{m+2},K{S^{″}}_{m+M_x}\right]/K,$

K is the number of repeated measurements of the high spatial resolution distributed fiber-optic acoustic sensor;

[0094] fourthly, performing eigenvalue decomposition on the covariance matrix R″ to obtain noise subspace E.sub.n, and calculating signal subspace E.sub.s(θ.sub.i,r.sub.i) corresponding to the sensor array by the sensing optical fiber laying method, wherein θ.sub.i represents the azimuth of the ith signal source, and r.sub.i represents the linear distance from the ith signal source to the central position of the sensor array; and calculating the loss function according to the following relationship:

[00008]$\Gamma \left(\theta ,r\right)={\frac{1}{\left[E_s\left(\theta ,r\right)E_n\right]}}^H{E}_n{E}_s\left(\theta ,r\right),$

wherein H represents conjugate transposition, and the orientation vector (θ.sub.i, r.sub.i) of the signal source with respect to the sensor array is obtained by searching the maximum value of the loss function and calculating; and

[0095] fifthly, calculating the two-dimensional space coordinates of the ith signal source according to the following relationship: x.sub.i=r.sub.i sin θ.sub.i, y.sub.i=r.sub.i cos θ.sub.i.

Embodiment III

[0096] In the third embodiment of the present invention, the distributed fiber-optic acoustic sensor adopts a phase-sensitive optical time domain reflectometer based on an ultra-narrow pulse and coherent detection, and as shown in FIG. 5, the sensing optical fiber is a SM-28 ordinary single-mode optical fiber/cable. The phase-sensitive optical time domain reflectometer based on the ultra-narrow pulse and coherent detection comprises a single-frequency narrow linewidth laser 8-1, a first fiber-optic coupler 8-2, a second fiber-optic coupler 8-6, an ultra-narrow pulse modulator 8-3, an optical amplifier 8-4, a fiber-optic circulator 8-5, a double-balanced photodetector 8-7, a data acquisition unit 8-8, a digital signal processing unit 8-9, and a sensing optical fiber interface 101, wherein laser light of the single-frequency narrow linewidth laser 8-1 is input to a first port 201 of the first fiber-optic coupler 8-2, and then is connected to the ultra-narrow pulse modulator 8-3 and a third port 603 of the second fiber-optic coupler 8-6 through a second port 202 and a third port 203 of the first fiber-optic coupler 8-2; the output of the ultra-narrow pulse modulator 8-3 is connected to the optical amplifier 8-4; the output of the optical amplifier 8-4 is connected to a first port 501 of the fiber-optic circulator 8-5, is output by a second port 502 and is connected to the sensing optical fiber interface 101; a third port 503 of the fiber-optic circulator 8-5 is connected to a first port 601 of the second fiber-optic coupler 8-6; an output signal of the second fiber-optic coupler 8-6 is output and connected to the double-balanced photodetector 8-7 through a second port 602 and a fourth port 604; and an output signal of the double-balanced photodetector 8-7 is connected to the data acquisition unit 8-8 to be amplified, and then is input to a digital signal processing unit 12-8.

[0097] The spatially coherent combination method of the present invention adopts a delay-and-sum beamforming algorithm; and after the subunits are enhanced, the subunits serve as the high sensitivity distributed fiber-optic sensors and do not perform the azimuth estimation and the spatial orientation of the signal sources. The implementation steps are as follows:

[0098] step (1): the high spatial resolution distributed fiber-optic acoustic sensor is connected to the sensing optical fiber, and is used for emitting the detection light pulse to the sensing optical fiber, collecting and quantitatively demodulating the sound fields of the signal sources sensed along the sensing optical fiber to obtain distributed sound field signals S(l,t), wherein l represents one-dimensional axial space of the sensing optical fiber, and t represents time;

[0099] step (2): the total number of the sensing units in the axial direction of the high spatial resolution distributed fiber-optic acoustic sensor is L, the sensing units are divided into M subunits in the axial direction, the number of the sensing units in each subunit is N, and L=M×N is met, wherein the sound field signals collected by the N sensing units in the mth subunit are:

[0100] [S.sub.m(l.sub.0+(m−1).Math.N.Math.Δz, t), S.sub.m(l.sub.0+(m−1).Math.N.Math.Δz+Δz, t), K S.sub.m(l.sub.0+(m−1).Math.N.Math.Δz+(N−1)Δz, t)], which is abbreviated as [S.sub.m,1, S.sub.m,2, K S.sub.m,N]; m is greater than or equal to 1 and less than or equal to M; l.sub.0 represents the starting position of the optical fiber; Δz represents the distance between the two adjacent sensing units; the minimum value of Δz is determined by the spatial resolution of the distributed fiber-optic acoustic sensor; the spatial resolution of the distributed fiber-optic acoustic sensor is centimeter level or higher; N is about tens or hundreds; and the spatial scale of each subunit is meter or submeter level;

[0101] step (3): calculating the time delay vector of the signal source with respect to each sensing unit in the subunit according to the spatial position information of each sensing unit in each subunit, for example, the time delay vector of the signal source with respect to each sensing unit in the mth subunit is a.sub.p (Δτ)=[Δτ.sub.m,1 Δτ.sub.m,2 L Δτ.sub.m,N] or the phase retardation is a.sub.p (Δφ)=[Δφ.sub.m,1 Δφ.sub.m,2 L Δφ.sub.m,N], wherein p=1, 2, L, P, P represents the number of the signal source;

[0102] step (4): calculating delay compensation weight W.sub.p(Δτ)=[−Δτ.sub.m,1 −Δτ.sub.m,2 L −Δτ.sub.m,N] according to the time delay vector a.sub.p(Δτ)=[Δτ.sub.m,1 Δτ.sub.m,2 L Δτ.sub.m,N] of the signal source with respect to each sensing unit in the mth subunit; obtaining the enhanced signal of the pth signal source by calculating according to the following mode: S′.sub.m=W.sub.p.Math.[S.sub.m,1, S.sub.m,2, K S.sub.m,N].sup.T, wherein T represents transposition operation, other P−1 signal sources and environmental noise can be suppressed, and the signal suppression effect of the P−1 signal sources is limited by the number of the sensing units in the entire subunit. After the sound field signals collected by the N sensing units in the mth subunit are enhanced by the above-mentioned array signal processing, the environmental noise can be effectively suppressed, the noise floor of the sensing system is reduced, the signal-to-noise ratio and sensitivity of the sensing system are improved, and at the same time, the interference fading and other problems of the fiber-optic sensing system can be solved, wherein the sound field signals S.sub.m(l.sub.0+(m−1).Math.N.Math.Δz,t), S.sub.m(l.sub.0+(m−1).Math.N.Math.Δz+Δz,t), K S.sub.m(+(m−1).Math.N.Math.Δz+(N−1)Δz,t) collected by the N sensing units in the mth subunit are enhanced by array signal processing to generate a signal S.sub.m (l.sub.0+(m−1).Math.N.Math.Δz+Δz, t), which is abbreviated as S′.sub.m; and after the M subunits are enhanced by array signal processing, obtained new sound field signals are [S′.sub.1, S′.sub.2, K S′.sub.m, K S′.sub.M]; and

[0103] step (5): each subunit can serve as the independent high-sensitivity sensor for sensing, and the process similar to the steps (3) and (4) can also be repeated again as needed to further realize signal enhancement.

[0104] The present invention is simple to implement, and can solve the problems of interference fading and limited signal-to-noise ratio of the distributed fiber-optic sensor; according to the present invention, a high signal-to-noise ratio sensor array can be constructed without changing the optical fiber structure, and the azimuth estimation and spatial orientation of the signal sources are realized; and the present invention is of great significance to the fields such as marine underwater acoustic detection and oil and gas exploration.

[0105] Although the present invention is described in detail with reference to the above-mentioned specific embodiments, it shall be understood that the present invention is not limited to the disclosed embodiments and examples, and for those skilled in the art, various changes can be made to the form and details of the present invention. For example, the working waveband of the laser can be replaced with other wavebands, and the circulator can be replaced with the coupler. It shall be understood that the above are only specific embodiments of the present invention and are not used to limit the present invention, and any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be comprised in the protection scope of the present invention.