CN102269333B - Method for eliminating pipe blockage acoustic signal strong interference by utilizing frequency domain self-adaptive filtering - Google Patents

Abstract

The invention relates to a method for eliminating pipe blockage acoustic signal strong interference by utilizing frequency domain self-adaptive filtering, The method comprises the following steps: (1) a main sensor is arranged in a pipeline, a reference sensor is arranged on the part close to a pump, and simultaneously a sound wave signal is collected; (2) the sound wave signals of the main sensor and the reference sensor sound wave are inputted a monitoring system; (3) a frequency domain spectrum of the inputted low-frequency sound wave signal is extracted; (4) the estimation is carried out on frequency domain spectrum of the signal, and the threshold value is set; (5) when the spectrum is higher than the threshold value, and an adaptive filter is started to eliminate the strong interference spectrum; and (6) the inverse fourier transform is utilized, thereby obtaining low-frequency blockage echo signal after the strong interference signal is eliminated. The system adopts a disturb-counteracted filter, the frequency domain spectrum adopts an adaptive algorithm to gradually iterative update, the required data-storing volume is smaller, thereby greatly improving the response time of the system; the system provided by the invention has the advantages that the construction is convenient, the reliability of the system is strong, the detection speed is rapid, therefore, the normal production run of the pipeline can not be influenced, and the blockage state of the oil gas pipeline is monitored in on-line manner.

Description

A kind of frequency domain adaptive filtering is eliminated the strongly disturbing method of line clogging acoustic signal

Technical field

The present invention relates to the removing method that oil and gas pipes stops up high reject signal in monitoring equipment.

Background technique

Crude oil that China produces is the higher waxy crude oil of condensation point and sticky heavy crude more than 80%, exists and the relevant potential safety hazard that flows so it is carried always.If the wax component in crude oil contacts with cryogenic object, will separate out and condense, can blocking pipe when serious.

Therefore, find that in time it is the key that ensures the submarine pipeline safe operation that the mobile sexual abnormality of submarine pipeline prevents trouble before it happens.Submarine pipeline stops up monitoring technology and is conducive in time to find and processes the flow abnormalities such as the leakage of submarine pipeline generation and obstruction, is the critical technological means of guarantee submarine pipeline safe operation.At present, submarine pipeline obstruction monitoring technology means relatively lack.

Low-frequency sound wave has the advantages that the decay brief biography is broadcast distance, is the important means of carrying out flowing state monitoring in pipeline.Monitoring pipeline blocking system based on sound wave utilizes acoustic emission apparatus emission low-frequency sound wave, and use sonic sensor that the low frequency echo signal of line clogging point reflection is converted to electrical signal, be transferred to oil gas pipe network blockage positioning server by data transmission network.The flowing state such as obstructions occurs when abnormal when pipeline, system can send warning automatically according to signal processing results, and calculates to stop up and put a position.

In practical operation, the acoustic signals impact that is interfered of interference acoustic signals that produce power is higher during the pump machine start and stop at oil and gas pipes two ends work, the obstruction reflection echo signal that sonic sensor collects has greatly increased the rate of false alarm of system.

Summary of the invention

For the problems referred to above, the invention provides a kind of adaptive frequency domain filter that utilizes and eliminate the method that the oil and gas pipes low frequency stops up interference acoustic signals in echo signal.This method can increase the reliability of system, reduces rate of false alarm.

For achieving the above object, the present invention takes following technological scheme: a kind of frequency domain adaptive filtering is eliminated the strongly disturbing method of line clogging acoustic signal, it is characterized in that, comprises the following steps:

1) master reference is set in pipeline, near pump machine place, reference sensor being set, gathers simultaneously acoustic signals;

2) to monitoring system input master reference and reference sensor acoustic signals;

3) its frequency domain spectra of low-frequency sound wave signal extraction to inputting;

4) frequency domain spectra of signal is estimated setting threshold;

5), start sef-adapting filter and eliminate the strong jamming frequency spectrum during higher than threshold value when frequency spectrum;

6) utilize Fourier inversion, the time domain low frequency after the high reject signal that is eliminated stops up echo signal.

In described step 3), the low frequency echo signal of frequency domain spectra extracting method stop up to(for) the oil and gas pipes that contains high reject signal is:

1. define oil and gas pipes obstruction echo signal model as follows:

y(n)＝s(n)+c(n) （1）

Wherein, y (n) is the acoustic signals that master reference collects, and s (n) is that low frequency stops up echo signal, and c (n) is the high reject signal that pipeline start and stop pump power traction rises;

2. input signal y (n) being divided into length is the time window that contains M sampling point, to the sampling point of the M in each time window, extracts frequency domain spectra with Fourier transformation, if M is not 2 integral number power, to N, N is 2 integral number power, that is: with this M sampling point zero padding

$Y\left(k\right)=\underset{n=0}{\overset{N-1}{\mathrm{\Σ}}}y\left(n\right)\·e^{-j\frac{2\mathrm{\π}}{N}\mathrm{kn}}---\left(2\right)$

K=1 wherein, 2 ..., N-1, y (n) expression master reference time-domain signal sequence, Y (k) is exactly the frequency domain spectra of current input signal, also referred to as Fourier transform spectrum;

3. wushu (1) transforms to frequency domain spectra:

Y(k)＝S(k)+C(k) （3）

Wherein, Y (k) is the frequency domain spectra of master reference signal, and S (k) is the frequency domain spectra of stopping up echo signal, and C (k) is the frequency domain spectra of high reject signal.

In step 5), utilize Adaptive Interference Cancelling Filter elimination strong jamming spectral method to be:

1. the reference undesired signal frequency domain spectra that gathers of computing reference sensor:

$C_1\left(k\right)=\underset{n=0}{\overset{N-1}{\mathrm{\Σ}}}{c}_1\left(n\right)\·e^{-j\frac{2\mathrm{\π}}{N}\mathrm{kn}}---\left(4\right)$

2. as reference sensor frequency domain spectra C ₁During (k) higher than threshold epsilon, according to pipeline field situation, C ₁(k) through the adjustable Interference Cancellation wave filter w (k) of parameter, the estimation of output undesired signal frequency domain spectra:

$\hat{C}\left(k\right)=C_1\left(k\right)w\left(k\right)---\left(5\right)$

Wherein w (k) is the Interference Cancellation wave filter;

3. deduct the interfering noise signal frequency domain spectra of estimating from master reference signal frequency domain spectrum Y (k) Be exactly system output signal frequency domain spectra Z (k):

$Z\left(k\right)=Y\left(k\right)-\hat{C}\left(k\right)=Y\left(k\right)-C_1\left(k\right)w\left(k\right)---\left(6\right)$

Z (k) approaches the frequency domain spectra S (k) of target signal.

Wherein, the step of adjusting Interference Cancellation filter parameter is as follows:

1. initialization, selected Interference Cancellation wave filter initial weight: w (k);

2. calculate t constantly, the output of Interference Cancellation wave filter:

3. Interference Cancellation: $Z\left(k\right)=Y\left(k\right)-\hat{C}\left(k\right)=S\left(k\right)+C\left(k\right)-C_1\left(k\right)w\left(k\right);$

4. the Interference Cancellation filter parameter upgrades: w (k+1)=w (k)+2 μ Z (k) Y (k), and adaptive step parameter μ adjusts according to actual conditions, and the μ value is less than 1 usually;

5. enter next constantly: t=t+1, jump to step 2., repeat 2.～4. iteration;

6. the output Z of system (k) is namely the frequency domain spectra of having eliminated the low frequency obstruction echo signal after the high reject signal.

Method that time domain low frequency after the high reject signal stops up echo signal is to utilize Fourier inversion to be eliminated:

$\hat{z}\left(n\right)=\underset{k=0}{\overset{N-1}{\mathrm{\Σ}}}Z\left(k\right)\·e^{j\frac{2\mathrm{\π}}{N}k\·n}---\left(7\right)$

Namely to have eliminated the time domain low frequency after the high reject signal to stop up echo signal.

The present invention is owing to taking above technological scheme, and it has the following advantages: 1, system adopts master reference, reference sensor collection signal.Reference sensor mainly gathers near the pump machine undesired signal that start and stop pump machine produces.Master reference mainly gathers the line clogging reflection echo signal, and effect of signals also can be interfered.So when the master reference signal was transformed to frequency domain spectra, its frequency spectrum also comprised obstruction reflection echo signal and undesired signal.The present invention can eliminate the undesired signal that contains in master reference discriminatively.2, system adopts the Interference Cancellation wave filter, and frequency domain spectra is adopted adaptive algorithm, and progressively iteration is upgraded, and the data volume of the required storage of system is less, has greatly improved the reaction time of system.3, system construction is convenient, and reliability is strong, and detection speed is fast, and on not impact of the normal production run of pipeline, stopping state that can the on-line monitoring oil and gas pipes.

Description of drawings

Fig. 1 is basic flow sheet of the present invention;

Fig. 2 adopts sef-adapting filter to eliminate the schematic diagram of pipeline undesired signal.

Embodiment

Below in conjunction with drawings and Examples, the present invention is described in detail.

As shown in Figure 1, the basic step of the inventive method is:

1) master reference is set in pipeline, near pump machine place, reference sensor being set, gathers simultaneously acoustic signals;

2) acoustic signals of input master reference and reference sensor in the monitoring system;

3) system extracts its frequency domain spectra to the acoustic signals of input;

4) frequency domain spectra of signal is estimated setting threshold;

5) when frequency spectrum during higher than threshold value, start sef-adapting filter and eliminate and disturb sound wave spectrum;

6) utilize Fourier inversion, the time domain low frequency after the high reject signal that is eliminated stops up echo signal.

In step 3), contain low frequency in the low-frequency sound wave signal due to the master reference input and stop up echo signal and strong jamming acoustic signals, so its frequency domain spectra extracting method is:

1. define oil and gas pipes obstruction echo signal model as follows:

y(n)＝s(n)+c(n) （1）

Wherein, y (n) is the signal that master reference collects, and usually, the signal that collects comprises that low frequency stops up echo signal and disturbs acoustic signals; S (n) is target signal, and namely low frequency stops up echo signal; C (n) is the high reject signal that the normal start and stop pump of pipeline power traction rises.

2. master reference input signal y (n) being divided into length is the time window that contains M sampling point, to the sampling point of the M in each time window, extracts frequency domain spectra with Fourier transformation, if M is not 2 integral number power, to N, N is 2 integral number power, that is: with this M sampling point zero padding

$Y\left(k\right)=\underset{n=0}{\overset{N-1}{\mathrm{\Σ}}}y\left(n\right)\·e^{-j\frac{2\mathrm{\π}}{N}\mathrm{kn}}---\left(2\right)$

K=1 wherein, 2 ..., N-1, y (n) expression master reference time-domain signal sequence, Y (k) is the frequency domain spectra that the master reference signal obtains by Fourier transformation.Frequency domain spectra Y (k) is analyzed, both can improve and screen low frequency obstruction echo signal, can reduce computational processing again.

3. wushu (1) transforms to frequency domain spectra:

Y(k)＝S(k)+C(k) （3）

Wherein, Y (k) is the frequency domain spectra of master reference acoustic signals; S (k) is the frequency domain spectra that master reference stops up echo signal; C (k) is the frequency domain spectra of master reference high reject signal.

In step 5), very near the pump machine, so the signal that collects of reference sensor is mainly the high reject signal that start and stop pump machine produces due to reference sensor.The reference undesired signal frequency domain spectra that the computing reference sensor gathers:

$C_1\left(k\right)=\underset{n=0}{\overset{N-1}{\mathrm{\Σ}}}{c}_1\left(n\right)\·e^{-j\frac{2\mathrm{\π}}{N}\mathrm{kn}}---\left(4\right)$

As reference sensor frequency domain spectra C ₁During (k) higher than threshold epsilon, system utilizes frequency domain adaptive filtering that the high reject signal of start and stop pump machine is effectively suppressed.Utilize processing method that sef-adapting filter suppresses undesired signal referring to Fig. 2:

According to the pipeline field situation, reference sensor undesired signal frequency domain spectra C ₁(k) after the adjustable Interference Cancellation wave filter w (k) of parameter, obtain the estimation of master reference undesired signal frequency domain spectra C (k)

$\hat{C}\left(k\right)=C_1\left(k\right)w\left(k\right)---\left(5\right)$

Wherein w (k) is the Interference Cancellation wave filter.Utilize adaptive algorithm to regulate the parameter of Interference Cancellation wave filter w (k), make the output signal frequency domain spectra

Approach master reference interfering noise signal frequency domain spectra C (k).

Then, deduct the interfering noise signal frequency domain spectra of estimating from master reference signal frequency domain spectrum Y (k)

System output signal frequency domain spectra Z (k):

Due to Interference Cancellation wave filter output frequency domain spectra

Approach the interfering noise signal frequency domain spectra C (k) of master reference signal superposition, i.e. C (k)-C ₁(k) w (k) levels off to 0, so the frequency domain spectra Z (k) of system's output will approach target signal frequency domain spectra S (k).

The parameter estimation algorithm of Interference Cancellation wave filter w (k) adopts Minimum Mean Square Error error criterion self adaption to estimate to obtain.The adaptive algorithm step is as follows:

1. initialization, selected Interference Cancellation wave filter initial weight: w (k);

2. calculate t constantly, the output of Interference Cancellation wave filter:

3. Interference Cancellation: $Z\left(k\right)=Y\left(k\right)-\hat{C}\left(k\right)=S\left(k\right)+C\left(k\right)-C_1\left(k\right)w\left(k\right);$

4. the Interference Cancellation filter parameter upgrades: w (k+1)=w (k)+2 μ Z (k) Y (k), and adaptive step parameter μ adjusts according to actual conditions, and the μ value is less than 1 usually.

5. enter next constantly: t=t+1, jump to step 2., repeat 2.～4. iteration;

6. the output Z of system (k) is namely the frequency domain spectra of having eliminated the low frequency obstruction echo signal after the high reject signal.

The time domain low frequency that utilizes the Fourier inversion algorithm to be eliminated after the high reject signal stops up echo signal:

$\hat{z}\left(n\right)=\underset{k=0}{\overset{N-1}{\mathrm{\Σ}}}Z\left(k\right)\·e^{j\frac{2\mathrm{\π}}{N}k\·n}---\left(7\right)$

Namely to have eliminated the time domain low frequency after the high reject signal to stop up echo signal.

Claims (4)

1. a frequency domain adaptive filtering is eliminated the strongly disturbing method of line clogging acoustic signal, it is characterized in that, comprises the following steps:

1) master reference is set in pipeline, near pump machine place, reference sensor being set, gathers simultaneously acoustic signals;

2) to monitoring system input master reference and reference sensor acoustic signals;

3) its frequency domain spectra of low-frequency sound wave signal extraction to inputting; The low frequency echo signal of frequency domain spectra extracting method stop up to(for) the oil and gas pipes that contains high reject signal is:

1. define oil and gas pipes obstruction echo signal model as follows:

y(n)＝s(n)+c(n) （1）

Wherein, y (n) is the acoustic signals that master reference collects, and s (n) is that low frequency stops up echo signal, and c (n) is the high reject signal that pipeline start and stop pump power traction rises;

2. input signal y (n) being divided into length is the time window that contains M sampling point, to the sampling point of the M in each time window, extracts frequency domain spectra with Fourier transformation, if M is not 2 integral number power, to N, N is 2 integral number power, that is: with this M sampling point zero padding

$Y\left(k\right)=\underset{n=0}{\overset{N-1}{\mathrm{\Σ}}}y\left(n\right)\·e^{-j\frac{2\mathrm{\π}}{N}\mathrm{kn}}---\left(2\right)$

K=1 wherein, 2 ..., N-1, y (n) expression master reference time-domain signal sequence, Y (k) is exactly the frequency domain spectra of current input signal, also referred to as Fourier transform spectrum;

3. wushu (1) transforms to frequency domain spectra:

Y(k)＝S(k)+C(k) （3）

Wherein, Y (k) is the frequency domain spectra of master reference signal, and S (k) is the frequency domain spectra of stopping up echo signal, and C (k) is the frequency domain spectra of high reject signal;

4) frequency domain spectra of signal is estimated setting threshold;

5), start sef-adapting filter and eliminate the strong jamming frequency spectrum during higher than threshold value when frequency spectrum;

6) utilize Fourier inversion, the time domain low frequency after the high reject signal that is eliminated stops up echo signal.

2. a kind of frequency domain adaptive filtering as claimed in claim 1 is eliminated the strongly disturbing method of line clogging acoustic signal, it is characterized in that, in described step 5), Adaptive Interference Cancelling Filter is eliminated the strong jamming spectral method and is:

1. the reference undesired signal frequency domain spectra that gathers of computing reference sensor:

$C_1\left(k\right)=\underset{n=0}{\overset{N-1}{\mathrm{\Σ}}}{c}_1\left(n\right)\·e^{-j\frac{2\mathrm{\π}}{N}\mathrm{kn}}---\left(4\right)$

c ₁(n) high reject signal that rises of the pipeline start and stop pump power traction of reference sensor collection;

2. as reference sensor frequency domain spectra C ₁During (k) higher than threshold epsilon, according to pipeline field situation, C ₁(k) through the adjustable Interference Cancellation wave filter w (k) of parameter, the estimation of output undesired signal frequency domain spectra:

$\hat{C}\left(k\right)=C_1\left(k\right)w\left(k\right)---\left(5\right)$

Wherein w (k) is the Interference Cancellation wave filter;

3. deduct the interfering noise signal frequency domain spectra of estimating from master reference signal frequency domain spectrum Y (k) Be exactly system output signal frequency domain spectra Z (k):

$Z\left(k\right)=Y\left(k\right)-\hat{C}\left(k\right)=Y\left(k\right)-C_1\left(k\right)w\left(k\right)---\left(6\right)$

Z (k) approaches the frequency domain spectra S (k) of target signal.

3. a kind of frequency domain adaptive filtering as claimed in claim 2 is eliminated the strongly disturbing method of line clogging acoustic signal, it is characterized in that, the step of regulating the Interference Cancellation filter parameter is as follows:

1. initialization, selected Interference Cancellation wave filter initial weight: w (k);

2. calculate t constantly, the output of Interference Cancellation wave filter:

3. Interference Cancellation: $Z\left(k\right)=Y\left(k\right)-\hat{C}\left(k\right)=S\left(k\right)+C\left(k\right)-C_1\left(k\right)w\left(k\right),$ C (k) is the frequency domain spectra of master reference high reject signal;

4. the Interference Cancellation filter parameter upgrades: w (k+1)=w (k)+2 μ Z (k) Y (k), and adaptive step parameter μ adjusts according to actual conditions, and the μ value is less than 1 usually;

5. enter next constantly: t=t+1, jump to step 2., repeat 2.～4. iteration;

6. the output Z of system (k) is namely the frequency domain spectra of having eliminated the low frequency obstruction echo signal after the high reject signal.

4. eliminate the strongly disturbing methods of line clogging acoustic signal as claim 1 or 2 or 3 described a kind of frequency domain adaptive filterings, it is characterized in that, method that the time domain low frequency after the high reject signal stops up echo signal is to utilize Fourier inversion to be eliminated:

$\hat{z}\left(n\right)=\underset{k=0}{\overset{N-1}{\mathrm{\Σ}}}Z\left(k\right)\·e^{j\frac{2\mathrm{\π}}{N}k\·n}---\left(7\right)$

Be namely to have eliminated the time domain low frequency after the high reject signal to stop up echo signal, Z (k) is the frequency domain spectra of having eliminated the low frequency obstruction echo signal after the high reject signal.

Images