CN109870693A - A kind of wind turbine clutter detection method - Google Patents
A kind of wind turbine clutter detection method Download PDFInfo
- Publication number
- CN109870693A CN109870693A CN201910149161.XA CN201910149161A CN109870693A CN 109870693 A CN109870693 A CN 109870693A CN 201910149161 A CN201910149161 A CN 201910149161A CN 109870693 A CN109870693 A CN 109870693A
- Authority
- CN
- China
- Prior art keywords
- distance unit
- processed
- wind turbine
- echo
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
The present invention relates to a kind of wind turbine clutter detection methods, it is widely different on spectrum width based on land clutter, meteorologic signal, wind turbine clutter,, frequency spectrum wide this characteristic strong using wind turbine clutter power, weather radar wind turbine clutter is detected by the otherness of spectrum width, it can effectively detect the distance unit where wind turbine clutter, and operand is small, is suitble to engineer application.
Description
Technical field
The present invention relates to a kind of wind turbine clutter detection methods, belong to weather radar signal detection technique field.
Background technique
In recent years, with the continuous growth of wind-powered electricity generation number and scale, by wind turbine blade in high-speed rotation
In caused clutter have become influence meteorological radar sounding performance principal element, existing radar system is in land clutter, meteorology
In the case that signal, turbine clutter exist simultaneously, the clutter suppression method taken is mostly to assume that clutter is in static or slow
In slow motion state, and the blade clutter in quickly movement is had ignored, and turn of the different piece of the turbine blade in movement
Moving-wire speed is also different.Therefore, the micro-doppler frequency spectrum for studying turbine blade is returned to the radar as caused by blade clutter is reduced
The error detection of wave signal and tracking error are of great significance.
To solve the above-mentioned problems, Cornwall proposes to select, and there is the material for absorbing wave beam function to make turbine leaf
Piece, so that radar cross section (Radar Cross Square, RCS) is reduced to inhibit turbine clutter, but due to material
Cost of manufacture is high, it is difficult to promote；Gallardo B et al. proposes the revolving speed by radar return power Spectral Estimation turbine blade,
Then judge the revolving speed whether priori the wind turbine blade range of speeds, turbine clutter is detected with this, this method is excessively
Rely on priori number；Domestic He Weikun et al. proposes the turbine electromechanics magnetic scattering characteristic Simulation based on mixed model, the model
Influence of the factors such as the reflection coefficient of electromagnetic wave to echo is considered, but longer using the modeling time needed for FEKO simulation software,
Estimation result can not quickly be obtained.
Summary of the invention
Technical problem to be solved by the invention is to provide a kind of wind turbine clutter detection methods, utilize wind turbine
Machine clutter power is strong, this wide characteristic of frequency spectrum, can effectively detect the distance list where weather radar wind turbine clutter
Member, and operand is small, is suitble to engineer application.
In order to solve the above-mentioned technical problem the present invention uses following technical scheme: the present invention devises a kind of wind turbine
Clutter detection method realizes wind based on weather radar to wind turbine machine testing, received sequences of echo signals { x (n) }
The detection of power turbine clutter, n ∈ { 1 ..., K }, K indicate weather radar to wind turbine machine testing, received echo-signal
Sum, x (n) indicates n-th of echo-signal received by weather radar, and detection method includes the following steps:
It is single that step B. sequence is directed between weather radar and wind turbine all distances in echo-signal movement routine
Member carries out average division by default number of segment M, obtains sector unit group, and obtains corresponding to sector unit group difference
Power threshold；It is then based on each power threshold, is rejected there is only each distance unit of land clutter signal, is obtained remaining each
Distance unit as each distance unit to be processed, and enters step C；
Step C. is based on sequences of echo signals { x (n) }, and obtaining each distance unit to be processed, corresponding spectrum width is estimated respectively
Evaluation σl, subsequently into step D；Wherein, l ∈ { 1 ..., N }, N indicates the quantity of distance unit to be processed, σlIndicate first to
Handle spectrum width estimated value corresponding to distance unit；
Spectrum width estimation values sigma step D. corresponding by each distance unit difference to be processedl, single for institute's distance to be handled
Member carries out ascending order arrangement, for sequence front end, sequence tail end, rejects default d distance unit to be processed respectively；Then it is remained
The average value of spectrum width estimated value corresponding to remaining distance unit to be handled, and default spectrum width amplification coefficient γ is combined, it is composed
Wide detection threshold value σ0, and enter step E；
Step E. is directed to step B respectively and obtains each distance unit to be processed, judges corresponding to distance unit to be processed
Spectrum width estimation values sigmalWhether spectrum width detection threshold value σ is greater than0, be indicate the distance unit to be processed there are wind turbine clutter,
Otherwise indicating the distance unit to be processed, there is no wind turbine clutters.
As a preferred technical solution of the present invention: the detection method further includes step A as follows, and step is first carried out
A, subsequently into step B；
Step A. carries out offseting processing for sequences of echo signals { x (n) }, realizes the inhibition of noise signal over the ground.
As a preferred technical solution of the present invention: in the step A, cascading institute using by two two pulse cancellers
The three pulse canceller of composition carries out offseting processing for sequences of echo signals { x (n) }, realizes the inhibition of noise signal over the ground.
As a preferred technical solution of the present invention, in the step A, respectively in sequences of echo signals { x (n) }
Each echo-signal x (n), obtain following formulae results；
x(n)-2x(n+1)+x(n+2)
And corresponding echo-signal x (n) is directed to the result and is updated, it is respectively completed the update of each echo-signal x (n),
Sequences of echo signals { x (n) } is carried out to offset processing, realizes the inhibition of noise signal over the ground.
As a preferred technical solution of the present invention, the step B includes the following steps:
Step B-1. sequence is directed between weather radar and wind turbine all distances in echo-signal movement routine
Unit carries out average division by default number of segment M, obtains M sections of distance unit groups, respectively includes Z=in sector unit group
L/M distance unit, subsequently into step B-2；Wherein, L indicates that echo-signal is mobile between weather radar and wind turbine
The sum of distance unit on path；
Step B-2. is directed to sector unit group respectively, as follows:
Obtain the corresponding power threshold P of sector unit group differencedm, subsequently into step B-3；Wherein, m ∈
{ 1 ..., M }, PdmIndicate that power threshold corresponding to m sections of distance unit groups, G indicate predetermined amplitude coefficient, PmiIndicate m
The power of i-th of distance unit echo-signal in section distance unit group；
Step B-3. is directed to each distance unit respectively, judges whether the power of distance unit echo-signal is less than where it
The corresponding power threshold of distance unit group is to determine that there is only land clutter signals in the distance unit, and it is single to reject the distance
Member；Otherwise the distance unit is not handled, after all distance unit complete above-mentioned judgement operation, enters step B-4；
Step B-4. is by remaining each distance unit, as each distance unit to be processed, subsequently into step C.
As a preferred technical solution of the present invention, the step C includes the following steps:
Step C-1. obtains the corresponding R of each distance unit difference to be processed0lAnd R1l, subsequently into step C-2；Its
In, l ∈ { 1 ..., N }, N indicates the quantity of distance unit to be processed, R0lIndicate first of distance unit echo-signal sequence to be processed
Column zero-lag autocorrelation parameter, R1lIndicate first of distance unit sequences of echo signals single order autocorrelation parameter to be processed；
Step C-2. is directed to each distance unit to be processed respectively, as follows:
Obtain the corresponding spectrum width estimation values sigma of each distance unit difference to be processedl, subsequently into step D；Wherein, σlTable
Show spectrum width estimated value corresponding to first of distance unit to be processed, λ indicates weather radar wavelength, TrIndicate that front and back is twice adjacent
Time interval between wave signal.
As a preferred technical solution of the present invention, the step C-1 includes the following:
It is directed to each distance unit to be processed respectively, as follows:
Obtain distance unit sequences of echo signals zero-lag autocorrelation parameter R to be processed0l, wherein l ∈ { 1 ..., N }, N
Indicate the quantity of distance unit to be processed, R0lIndicate first of distance unit sequences of echo signals zero-lag auto-correlation ginseng to be processed
Number；yl(n) n-th of echo-signal of first of distance unit to be processed is indicated,Indicate yl(n) conjugation；
Meanwhile according to the following formula:
Obtain the sequences of echo signals single order autocorrelation parameter R of the distance unit to be processed1l, wherein R1lIt indicates first
The sequences of echo signals single order autocorrelation parameter of distance unit to be processed；yl(b+1) it indicates in first of distance unit to be processed
The b+1 echo-signal,Indicate yl(b) conjugation, yl(b) b-th of echo letter of first of distance unit to be processed is indicated
Number；
And then obtain the corresponding R of each distance unit difference to be processed0lAnd R1l, subsequently into step C-2.
As a preferred technical solution of the present invention: single for remaining institute's distance to be handled in the step D
The average value of spectrum width estimated value corresponding to member obtains spectrum width detection threshold by the product with the default spectrum width amplification coefficient γ
Value σ0。
A kind of wind turbine clutter detection method of the present invention, compared with the prior art by using the above technical solution,
It has following technical effect that
Wind turbine clutter detection method designed by the present invention is based on land clutter, meteorologic signal, wind turbine clutter
Widely different on spectrum width, frequency spectrum wide this characteristic strong using wind turbine clutter power is examined by the otherness of spectrum width
Weather radar wind turbine clutter is surveyed, can effectively detect the distance unit where wind turbine clutter, and operand
It is small, it is suitble to engineer application.
Detailed description of the invention
Fig. 1 is wind turbine and radar relational graph；
Fig. 2 is signal processing flow figure of the present invention；
Fig. 3 (a) is that three pulses offset preceding radar echo signal (land clutter+meteorologic signal) power spectrum；
Fig. 3 (b) is that three pulses offset rear radar echo signal (land clutter+meteorologic signal) power spectrum；
Fig. 4 (a) is that three pulses offset preceding radar echo signal (land clutter+WTC) power spectrum；
Fig. 4 (b) is that three pulses offset rear radar echo signal (land clutter+WTC) power spectrum；
Fig. 5 (a) is that three pulses offset preceding radar echo signal (meteorologic signal+land clutter+WTC) power spectrum；
Fig. 5 (b) is that three pulses offset rear radar echo signal (meteorologic signal+land clutter+WTC) power spectrum；
Fig. 6 is the estimation of meteorologic signal spectrum width；
Fig. 7 is the estimation of meteorologic signal+WTC spectrum width.
Specific embodiment
Specific embodiments of the present invention will be described in further detail with reference to the accompanying drawings of the specification.
The present invention devises a kind of wind turbine clutter detection method, based on weather radar to wind turbine machine testing,
Received sequences of echo signals { x (n) }, realizes the detection of wind turbine clutter, and n ∈ { 1 ..., K }, K indicate meteorological thunder
Up to the sum to wind turbine machine testing, received echo-signal, x (n) indicates n-th of echo received by weather radar
Signal.
For each echo-signal x (n) in sequences of echo signals { x (n) }, expression formula is x (n)=W (n)+C
(n)+S (n)+U (n) or x (n)=C (n)+S (n)+U (n), wherein W (n) indicates wind turbine noise signal, and C (n) is indicated
Land clutter signal, S (n) indicate that meteorologic signal, U (n) indicate noise.In practical application, that is, indicate to carry out echo-signal x (n)
It says, wherein may include wind turbine noise signal W (n), it is also possible to not include wind turbine noise signal W (n).
Land clutter signal C (n): using distance unit partitioning, and the width of each distance unit isfsTo adopt
Sample frequency, c are the light velocity.One distance unit is divided into 180 clutter scattering units, the clutter data of each distance unit is
For the superposition of this 180 clutter scattering units, the land clutter echo-signal of a-th of distance unit are as follows:
ω in formulat=2 π fd/frFor time domain angular frequency, RaFor the distance of a-th of distance unit, F is antenna gain, fdIt is more
General Le frequency, frFor pulse recurrence frequency, n ∈ { 1 ..., K }, K indicate weather radar to wind turbine machine testing, received
The sum of echo-signal, L indicate the sum of distance unit in echo-signal movement routine between weather radar and wind turbine.
Meteorologic signal S (n): meteorological target is distributed object, and echo-signal can be folded by the echo of multiple scattering units
Add to be formed.Coherent superposition, as meteorologic signal echo are carried out to the echo of all scattering units in a-th of distance unit.A-th
Distance unit, meteorologic signal can be indicated in n-th of echo signal sample are as follows:
Wherein, I is the sum of scattering unit in a-th of distance unit, AuFor the amplitude of meteorological intended particle u, by radar
Equation determines.ωtFor time domain angular frequency, can be obtained by following formula:
vuFor the radial velocity of meteorological intended particle u, λ is radar wavelength
Wind turbine signal: wind turbine is made of blade, mast and engine room, wind turbine clutter mainly by
Blade and mast clutter composition.Blade is rotation, mast be it is static, using scattering point Additive Model, wind turbine with
Radar relationship is as shown in Figure 1.
Blade signals: there are three blades for wind turbine, their interval is 120 degree, so between any two in angle
There is one 120 degree of lag, consider a blade first, then other two blades only differ 120 degree.Each blade is seen
It is formed at by multiple scattering points, as soon as the clutter of blade is exactly the superposition of these scattering points, the superposition of three blade scattering points
It is wind turbine blade clutter.The echo-signal of one scattering point p of wind turbine blade:
In formula, fcFor radar frequency,For time delay, RpFor the distance of scattering point p to radar:
RgFor the distance of wind turbine to radar, riIt is distance of the blade scattering point to blade axle center,It is with reference to blade
With the angle of radar line of sight:
α, β are respectively azimuth and pitch angle, and θ is the angle with reference to blade and Y-axis positive direction.
Remove carrier wave, baseband signal can be obtained are as follows:Then the total echo-signal of all blades is (wind-force whirlpool
Turbine blade number is generally 3):
P is the sum of blade scattering particles in formula.
Mast signal: mast be it is static, mast clutter is also considered as the superposition of multiple scattering points.Then mast scattering point q
Clutter are as follows:RqIt is distance of the scattering point q to radar:
hiIt is distance of the scattering point q to blade axle center.
Remove carrier wave, baseband signal can be obtained:Then the echo-signal of entire mast is all scattering points
Superposition:In conclusion wind turbine noise signal:
A kind of wind turbine clutter detection method designed by the present invention, in practical application, as shown in Fig. 2, specific packet
Include following steps.
Step A. carries out offseting processing for sequences of echo signals { x (n) }, realizes the inhibition of noise signal over the ground, then
Enter step B.
Above-mentioned steps A in practical applications, using cascading constituted three pulse canceller by two two pulse cancellers,
Respectively for each echo-signal x (n) in sequences of echo signals { x (n) }, following formulae results are obtained；
x(n)-2x(n+1)+x(n+2)
And corresponding echo-signal x (n) is directed to the result and is updated, it is respectively completed the update of each echo-signal x (n),
Sequences of echo signals { x (n) } is carried out to offset processing, realizes the inhibition of noise signal over the ground.
Since there is only the distance unit of land clutter after offseting by three pulses, corresponding performance number is largely
On be inhibited, according to power threshold reject there is only the distance unit of land clutter, by wind turbine clutter and meteorologic signal
Distance unit that may be present extracts, and then reduces the range of follow-up data processing, that is, continues to execute following steps.
The how far that step B. studies detective distance radar is different, the power of the radar return data of different distance unit
Value is also different, and distance is remoter, and performance number is smaller, on the contrary then bigger, calculates accuracy to improve, proposes segmental averaging method.Sequentially
For all distance unit between weather radar and wind turbine in echo-signal movement routine, carried out by default number of segment M
It is average to divide, sector unit group is obtained, and obtain the corresponding power threshold of sector unit group difference；It is then based on
Each power threshold is rejected there is only each distance unit of land clutter signal, and remaining each distance unit is obtained, as it is each to
Distance unit is handled, and enters step C.
In above-mentioned steps B practical application, following steps are specifically executed:
Step B-1. sequence is directed between weather radar and wind turbine all distances in echo-signal movement routine
Unit carries out average division by default number of segment M, obtains M sections of distance unit groups, respectively includes Z=in sector unit group
L/M distance unit, subsequently into step B-2；Wherein, L indicates that echo-signal is mobile between weather radar and wind turbine
The sum of distance unit on path.
Step B-2. is directed to sector unit group respectively, as follows:
Obtain the corresponding power threshold P of sector unit group differencedm, subsequently into step B-3；Wherein, m ∈
{ 1 ..., M }, PdmIndicate that power threshold corresponding to m sections of distance unit groups, G indicate predetermined amplitude coefficient, practical application is worked as
In, G=1.2, PmiIndicate the power of i-th of distance unit echo-signal in m sections of distance unit groups.
Step B-3. is directed to each distance unit respectively, judges whether the power of distance unit echo-signal is less than where it
The corresponding power threshold of distance unit group is to determine that there is only land clutter signals in the distance unit, and it is single to reject the distance
Member；Otherwise the distance unit is not handled, after all distance unit complete above-mentioned judgement operation, enters step B-4.
Step B-4. is by remaining each distance unit, as each distance unit to be processed, subsequently into step C.
Step C. is based on sequences of echo signals { x (n) }, and obtaining each distance unit to be processed, corresponding spectrum width is estimated respectively
Evaluation σl, subsequently into step D；Wherein, l ∈ { 1 ..., N }, N indicates the quantity of distance unit to be processed, σlIndicate first to
Handle spectrum width estimated value corresponding to distance unit.
In above-mentioned steps C practical application, following steps are specifically executed:
Step C-1. obtains the corresponding R of each distance unit difference to be processed0lAnd R1l, subsequently into step C-2；Its
In, l ∈ { 1 ..., N }, N indicates the quantity of distance unit to be processed, R0lIndicate first of distance unit echo-signal sequence to be processed
Column zero-lag autocorrelation parameter, R1lIndicate first of distance unit sequences of echo signals single order autocorrelation parameter to be processed.
In practical application, the corresponding R of each distance unit difference to be processed0lAnd R1l, obtained by following specific execute
?.
It is directed to each distance unit to be processed respectively, as follows:
Obtain distance unit sequences of echo signals zero-lag autocorrelation parameter R to be processed0l, wherein l ∈ { 1 ..., N }, N
Indicate the quantity of distance unit to be processed, R0lIndicate first of distance unit sequences of echo signals zero-lag auto-correlation ginseng to be processed
Number；yl(n) n-th of echo-signal of first of distance unit to be processed is indicated,Indicate yl(n) conjugation.
Meanwhile according to the following formula:
Obtain the sequences of echo signals single order autocorrelation parameter R of the distance unit to be processed1l, wherein R1lIt indicates first
The sequences of echo signals single order autocorrelation parameter of distance unit to be processed；yl(b+1) it indicates in first of distance unit to be processed
The b+1 echo-signal,Indicate yl(b) conjugation, yl(b) b-th of echo letter of first of distance unit to be processed is indicated
Number.
And then obtain the corresponding R of each distance unit difference to be processed0lAnd R1l。
Step C-2. is directed to each distance unit to be processed respectively, as follows:
Obtain the corresponding spectrum width estimation values sigma of each distance unit difference to be processedl, subsequently into step D；Wherein, σlTable
Show spectrum width estimated value corresponding to first of distance unit to be processed, λ indicates weather radar wavelength, TrIndicate that front and back is twice adjacent
Time interval between wave signal.
Spectrum width estimation values sigma step D. corresponding by each distance unit difference to be processedl, single for institute's distance to be handled
Member carries out ascending order arrangement, for sequence front end, sequence tail end, rejects default d distance unit to be processed respectively, practical application is worked as
In, it can design for sequence front end, sequence tail end, reject preset 2 distance unit to be processed respectively；Then remaining institute is obtained
The average value of spectrum width estimated value corresponding to distance unit to be handled is obtained by the product with default spectrum width amplification coefficient γ
Spectrum width detection threshold value σ0, and enter step E.
Step E. is directed to step B respectively and obtains each distance unit to be processed, judges corresponding to distance unit to be processed
Spectrum width estimation values sigmalWhether spectrum width detection threshold value σ is greater than0, be indicate the distance unit to be processed there are wind turbine clutter,
Otherwise indicating the distance unit to be processed, there is no wind turbine clutters.
It is applied to wind turbine clutter detection method designed by the present invention in specific embodiment, radar system is imitative
True parameter is as shown in table 1 below, and wind turbine system simulation parameter is as shown in table 2 below.
| Carrier frequency f0 | 5.5GHz |
| Pulse recurrence frequency fr | 1000Hz |
| Pulse number | 64 |
| Radar altitude (m) | 1000 |
Table 1
| Length of blade (m) | 26 |
| Mast height (m) | 70 |
| Distance (km) of the axle center to radar | 36 |
| Rotating vane revolving speed (r/min) | 15 |
| Surfaces of revolution angle | 73.5° |
| With reference to the initial angle of blade | 90° |
Table 2
In this emulation experiment, radar echo signal is offseted by three pulses first, filters off land clutter and wind turbine
Machine mast clutter.If there are meteorologic signals for the 1st to the 50th distance unit, each distance unit spectrum width is estimated respectively.?
Wind turbine clutter is added in 25 distance unit, estimates the distance unit spectrum width.Compare the 25th distance unit and other
The size relation of distance unit spectrum width calculates spectrum width detection threshold value according to the relationship of the two.
In this emulation experiment, land clutter is clutter, and Doppler frequency center is in 36Hz, spectrum width 20Hz.Due to three
Pulse canceller is merely able to effectively inhibit the clutter near zero-frequency, and place first has to carry out frequency before carrying out three pulses and offseting
Spectrum is moved, by the Doppler frequency central region of land clutter to zero-frequency.Before three pulses offset, radar return (land clutter+meteorology letter
Number) shown in power spectrum such as Fig. 3 (a), after three pulses offset, radar return (land clutter+meteorologic signal) power spectrum such as Fig. 3 (b) institute
Show.It can be obtained from the figure that going out, before three pulses offset, land clutter power is about 38dB, after three pulses offset, under land clutter power
- 20dB is dropped to, 58dB is had dropped, land clutter has obtained effective inhibition.Before three pulses offset, meteorologic signal power is about
48dB is offseted by three pulses, and meteorologic signal power rise has arrived 55dB, it can be seen that three pulse canceller to meteorologic signal without
Larger impact.Before three pulses offset, shown in radar return (land clutter+WTC) power spectrum such as Fig. 4 (a), after three pulses offset, thunder
Up to shown in echo (land clutter+WTC) power spectrum such as Fig. 4 (b).It can be obtained from the figure that before three pulses offset, the signal (wind-force at zero-frequency
Turbine mast clutter+land clutter) power is about 49dB, and after three pulses offset, 0dB is dropped to, has dropped 49dB.And
The clutter power of other frequency ranges offsets front and back in three pulses and is not much different.Fig. 5 (a) is that three pulses offset preceding radar return
(meteorologic signal+land clutter+WTC) power spectrum, Fig. 5 (b) are that three pulses offset rear radar return (meteorologic signal+land clutter+WTC)
Power spectrum.
After three pulse canceller, land clutter and wind turbine mast clutter have obtained effective inhibition, are left gas
Picture signals and wind turbine blade clutter.Meteorologic signal and wind turbine blade clutter have very big difference on spectrum width.It is first
First assume the 1 to 50th distance unit only exist meteorologic signal, spectrum width estimation is carried out to it, the result obtained as shown in fig. 6, by
Emulation experiment it is found that meteorologic signal spectrum width σsIn 0.36-0.54m/s.WTC is added in the 25th distance unit, to the distance
Unit carries out spectrum width estimation, as shown in fig. 7, obtaining the distance unit spectrum width is σd=5.0206m/s.Compare σdWith(For σsIt asks
Average gained) relationship, it calculatesSpectrum width amplification coefficient γ=6 can be set, spectrum width detection threshold value σ can be obtained0=6 × σs.When
Distance unit spectrum width σ to be detectedl> σ0When, there are wind turbine clutters for the distance unit.
Wind turbine clutter detection method designed by above-mentioned technical proposal is based on land clutter, meteorologic signal, wind turbine
Machine clutter is widely different on spectrum width, and, frequency spectrum wide this characteristic strong using wind turbine clutter power passes through the difference of spectrum width
Property detects weather radar wind turbine clutter, can effectively detect the distance unit where wind turbine clutter, and
Operand is small, is suitble to engineer application.
Embodiments of the present invention are explained in detail above in conjunction with attached drawing, but the present invention is not limited to above-mentioned implementations
Mode within the knowledge of a person skilled in the art can also be without departing from the purpose of the present invention
It makes a variety of changes.
Claims (8)
1. a kind of wind turbine clutter detection method, based on weather radar to wind turbine machine testing, received echo
Signal sequence { x (n) }, realizes the detection of wind turbine clutter, and n ∈ { 1 ..., K }, K indicate weather radar to wind turbine
It detects, the sum of received echo-signal, x (n) indicates that n-th of echo-signal received by weather radar, feature exist
In detection method includes the following steps:
Step B. sequence is directed between weather radar and wind turbine all distance unit in echo-signal movement routine, presses
Default number of segment M carries out average division, obtains sector unit group, and obtains the corresponding power of sector unit group difference
Threshold value；It is then based on each power threshold, is rejected there is only each distance unit of land clutter signal, it is single to obtain remaining each distance
Member as each distance unit to be processed, and enters step C；
Step C. is based on sequences of echo signals { x (n) }, obtains the corresponding spectrum width estimated value of each distance unit difference to be processed
σl, subsequently into step D；Wherein, l ∈ { 1 ..., N }, N indicates the quantity of distance unit to be processed, σlIndicate first it is to be processed
Spectrum width estimated value corresponding to distance unit；
Spectrum width estimation values sigma step D. corresponding by each distance unit difference to be processedl, for institute's distance unit to be handled into
The arrangement of row ascending order rejects default d distance unit to be processed for sequence front end, sequence tail end respectively；Then remaining institute is obtained
The average value of spectrum width estimated value corresponding to distance unit to be handled, and default spectrum width amplification coefficient γ is combined, obtain spectrum width inspection
Survey threshold value σ0, and enter step E；
Step E. is directed to step B respectively and obtains each distance unit to be processed, judges spectrum width corresponding to distance unit to be processed
Estimation values sigmalWhether spectrum width detection threshold value σ is greater than0, it is to indicate the distance unit to be processed there are wind turbine clutters, otherwise
Indicating the distance unit to be processed, there is no wind turbine clutters.
2. a kind of wind turbine clutter detection method according to claim 1, it is characterised in that: the detection method is also wrapped
It is as follows to include step A, step A is first carried out, subsequently into step B；
Step A. carries out offseting processing for sequences of echo signals { x (n) }, realizes the inhibition of noise signal over the ground.
3. a kind of wind turbine clutter detection method according to claim 2, it is characterised in that: in the step A, use
Constituted three pulse canceller is cascaded by two two pulse cancellers, carries out offseting place for sequences of echo signals { x (n) }
Reason realizes the inhibition of noise signal over the ground.
4. a kind of wind turbine clutter detection method according to claim 2, which is characterized in that in the step A, respectively
For each echo-signal x (n) in sequences of echo signals { x (n) }, following formulae results are obtained；
x(n)-2x(n+1)+x(n+2)
And corresponding echo-signal x (n) is directed to the result and is updated, pair it is respectively completed the update of each echo-signal x (n), i.e.,
Sequences of echo signals { x (n) } carries out offseting processing, realizes the inhibition of noise signal over the ground.
5. a kind of wind turbine clutter detection method according to claim 1, which is characterized in that the step B includes such as
Lower step:
Step B-1. sequence is directed between weather radar and wind turbine all distance unit in echo-signal movement routine,
Average division is carried out by default number of segment M, obtains M sections of distance unit groups, respectively includes Z=L/M in sector unit group
Distance unit, subsequently into step B-2；Wherein, L is indicated between weather radar and wind turbine in echo-signal movement routine
The sum of distance unit；
Step B-2. is directed to sector unit group respectively, as follows:
Obtain the corresponding power threshold P of sector unit group differencedm, subsequently into step B-3；Wherein, m ∈ 1 ...,
M }, PdmIndicate that power threshold corresponding to m sections of distance unit groups, G indicate predetermined amplitude coefficient, PmiIndicate that m sections of distances are single
The power of i-th of distance unit echo-signal in tuple；
Step B-3. is directed to each distance unit respectively, judges whether the power of distance unit echo-signal is less than distance where it
The corresponding power threshold of unit group is to determine that there is only land clutter signals in the distance unit, and reject the distance unit；It is no
The distance unit is not handled then, after all distance unit complete above-mentioned judgement operation, enters step B-4；
Step B-4. is by remaining each distance unit, as each distance unit to be processed, subsequently into step C.
6. a kind of wind turbine clutter detection method according to claim 1, which is characterized in that the step C includes such as
Lower step:
Step C-1. obtains the corresponding R of each distance unit difference to be processed0lAnd R1l, subsequently into step C-2；Wherein, l ∈
{ 1 ..., N }, N indicates the quantity of distance unit to be processed, R0lIndicate that first of distance unit sequences of echo signals zero to be processed is prolonged
Slow autocorrelation parameter, R1lIndicate first of distance unit sequences of echo signals single order autocorrelation parameter to be processed；
Step C-2. is directed to each distance unit to be processed respectively, as follows:
Obtain the corresponding spectrum width estimation values sigma of each distance unit difference to be processedl, subsequently into step D；Wherein, σlIndicate the
Spectrum width estimated value corresponding to l distance unit to be processed, λ indicate weather radar wavelength, TrIndicate the adjacent two echoes letter in front and back
Time interval between number.
7. a kind of wind turbine clutter detection method according to claim 6, which is characterized in that the step C-1 includes
It is as follows:
It is directed to each distance unit to be processed respectively, as follows:
Obtain distance unit sequences of echo signals zero-lag autocorrelation parameter R to be processed0l, wherein l ∈ { 1 ..., N }, N are indicated
The quantity of distance unit to be processed, R0lIndicate first of distance unit sequences of echo signals zero-lag autocorrelation parameter to be processed；yl
(n) n-th of echo-signal of first of distance unit to be processed is indicated,Indicate yl(n) conjugation；
Meanwhile according to the following formula:
Obtain the sequences of echo signals single order autocorrelation parameter R of the distance unit to be processed1l, wherein R1lIndicate first it is to be processed
The sequences of echo signals single order autocorrelation parameter of distance unit；yl(b+1) b+1 in first of distance unit to be processed is indicated
A echo-signal,Indicate yl(b) conjugation, yl(b) b-th of echo-signal of first of distance unit to be processed is indicated；Into
And obtain the corresponding R of each distance unit difference to be processed0lAnd R1l, subsequently into step C-2.
8. a kind of wind turbine clutter detection method according to claim 1, it is characterised in that: in the step D, for
The average value of spectrum width estimated value corresponding to remaining the distance unit to be handled, by with the default spectrum width amplification coefficient
The product of γ obtains spectrum width detection threshold value σ0。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910149161.XA CN109870693B (en) | 2019-02-28 | 2019-02-28 | Wind turbine clutter detection method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910149161.XA CN109870693B (en) | 2019-02-28 | 2019-02-28 | Wind turbine clutter detection method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109870693A true CN109870693A (en) | 2019-06-11 |
| CN109870693B CN109870693B (en) | 2022-06-24 |
Family
ID=66919432
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910149161.XA Active CN109870693B (en) | 2019-02-28 | 2019-02-28 | Wind turbine clutter detection method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109870693B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110261836A (en) * | 2019-06-27 | 2019-09-20 | 中国航空工业集团公司雷华电子技术研究所 | A method of inhibit to cause decoy by blades of large-scale wind driven generator |
| CN110297247A (en) * | 2019-07-22 | 2019-10-01 | 河海大学 | Weather radar wind power plant clutter suppression method based on the sparse recovery of low-rank matrix |
| CN111257885A (en) * | 2020-03-17 | 2020-06-09 | 河海大学 | Meteorological radar wind power plant clutter suppression method based on extreme learning machine |
| CN111624556A (en) * | 2020-06-08 | 2020-09-04 | 河海大学 | Meteorological radar WTC (wind turbine controller) inhibition method based on morphological component analysis |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106342236B (en) * | 2009-03-16 | 2012-12-19 | 中国航空工业集团公司雷华电子技术研究所 | Based on the airborne radar clutter suppression method of difference beam |
| CN102854496A (en) * | 2012-09-19 | 2013-01-02 | 中国民航大学 | Airborne meteorological radar ground clutter suppression method based on double threshold control |
| CN103383448A (en) * | 2013-06-25 | 2013-11-06 | 西安电子科技大学 | Clutter suppression method suitable for high pulse repetition frequency (HPRF) waveform airborne radar |
| CN104931949A (en) * | 2015-06-15 | 2015-09-23 | 中国民航大学 | Aircraft target detection method under radar scanning mode in wind turbine clutter background |
| CN109116326A (en) * | 2018-09-27 | 2019-01-01 | 中国科学院电子学研究所苏州研究院 | A kind of adaption radar ocean clutter cancellation method based on medion estimator |
-
2019
- 2019-02-28 CN CN201910149161.XA patent/CN109870693B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106342236B (en) * | 2009-03-16 | 2012-12-19 | 中国航空工业集团公司雷华电子技术研究所 | Based on the airborne radar clutter suppression method of difference beam |
| CN102854496A (en) * | 2012-09-19 | 2013-01-02 | 中国民航大学 | Airborne meteorological radar ground clutter suppression method based on double threshold control |
| CN103383448A (en) * | 2013-06-25 | 2013-11-06 | 西安电子科技大学 | Clutter suppression method suitable for high pulse repetition frequency (HPRF) waveform airborne radar |
| CN104931949A (en) * | 2015-06-15 | 2015-09-23 | 中国民航大学 | Aircraft target detection method under radar scanning mode in wind turbine clutter background |
| CN109116326A (en) * | 2018-09-27 | 2019-01-01 | 中国科学院电子学研究所苏州研究院 | A kind of adaption radar ocean clutter cancellation method based on medion estimator |
Non-Patent Citations (1)
| Title |
|---|
| 夏鹏等: "《基于形态分量分析的风力发电机杂波抑制方法》", 《空军预警学院学报》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110261836A (en) * | 2019-06-27 | 2019-09-20 | 中国航空工业集团公司雷华电子技术研究所 | A method of inhibit to cause decoy by blades of large-scale wind driven generator |
| CN110261836B (en) * | 2019-06-27 | 2023-03-03 | 中国航空工业集团公司雷华电子技术研究所 | Method for inhibiting false target caused by large-scale wind driven generator blade |
| CN110297247A (en) * | 2019-07-22 | 2019-10-01 | 河海大学 | Weather radar wind power plant clutter suppression method based on the sparse recovery of low-rank matrix |
| CN110297247B (en) * | 2019-07-22 | 2022-06-24 | 河海大学 | Meteorological radar wind power plant clutter suppression method based on low-rank matrix sparse recovery |
| CN111257885A (en) * | 2020-03-17 | 2020-06-09 | 河海大学 | Meteorological radar wind power plant clutter suppression method based on extreme learning machine |
| CN111257885B (en) * | 2020-03-17 | 2022-04-15 | 河海大学 | Meteorological radar wind power plant clutter suppression method based on extreme learning machine |
| CN111624556A (en) * | 2020-06-08 | 2020-09-04 | 河海大学 | Meteorological radar WTC (wind turbine controller) inhibition method based on morphological component analysis |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109870693B (en) | 2022-06-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109870693A (en) | A kind of wind turbine clutter detection method | |
| AU653812B2 (en) | Microburst precursor detector utilizing microwave radar | |
| CN104914415B (en) | Single-pulse radar coherent jamming method based on target range profile template matching | |
| CN111913158B (en) | Radar signal processing method for detecting low-speed small target under complex clutter background | |
| CN104931949B (en) | Wind turbine clutter background is got off the plane object detection method under radar scanning pattern | |
| CA2649887C (en) | A method for filtering sea clutter in a radar echo using a hydrographic model | |
| RU2419815C1 (en) | Method of tracking air target of "turbojet engine aircraft" class | |
| Janssen et al. | The shape of Doppler spectra from precipitation | |
| KR20130140693A (en) | Radar system | |
| CN108802722A (en) | It is a kind of based on tracking before the Faint target detection virtually composed | |
| CN103176187B (en) | A kind of airborne early warning radar ground surface high speed highway goal filtering method | |
| CN108181620B (en) | Three-coordinate radar trace point quality evaluation method | |
| Norin | A quantitative analysis of the impact of wind turbines on operational Doppler weather radar data | |
| CN111796288B (en) | Clutter frequency spectrum compensation technology-based three-coordinate radar moving target processing method | |
| US7548187B1 (en) | Adaptive clutter filtering to improve high sub-clutter visibility radar detection performance | |
| Molchanov et al. | Aerial target classification by micro-Doppler signatures and bicoherence-based features | |
| RU2617110C1 (en) | Method to support group air targets of "aircraft with turbojet" class in radar location station at exposure of rate interference | |
| CN107643515A (en) | Turbulent flow target velocity spectrum width method of estimation based on space-time covariance fitting | |
| CN110261840A (en) | A kind of quasi- high pulse repetition frequency radar waveform group and its signal processing method | |
| CN110244289A (en) | A kind of adaptive particle filter ground wave radar target integrative detection method | |
| CN115113208A (en) | Continuous wave radar clutter cancellation method based on accurate clutter feature recognition | |
| Sozen et al. | Scatter and Doppler effect of wind power plants to land radars | |
| Goncharenko et al. | Adaptive moving target indication in a windblown clutter environment | |
| CN111198358B (en) | Multi-rotor unmanned aerial vehicle echo accumulation method based on correlation detection | |
| Kulpa et al. | Parametric detection of the helicopter hub echo |
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 |