CN112629838B - Wind turbine blade fault monitoring method - Google Patents

Wind turbine blade fault monitoring method Download PDF

Info

Publication number
CN112629838B
CN112629838B CN202011269371.1A CN202011269371A CN112629838B CN 112629838 B CN112629838 B CN 112629838B CN 202011269371 A CN202011269371 A CN 202011269371A CN 112629838 B CN112629838 B CN 112629838B
Authority
CN
China
Prior art keywords
wind turbine
turbine blade
blade
wind
time
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011269371.1A
Other languages
Chinese (zh)
Other versions
CN112629838A (en
Inventor
黄力
代朝阳
刘钢
刘鹏
谢黄海
石青松
廖梓榕
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Three Gorges University CTGU
Original Assignee
China Three Gorges University CTGU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Three Gorges University CTGU filed Critical China Three Gorges University CTGU
Priority to CN202011269371.1A priority Critical patent/CN112629838B/en
Publication of CN112629838A publication Critical patent/CN112629838A/en
Application granted granted Critical
Publication of CN112629838B publication Critical patent/CN112629838B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M13/00Testing of machine parts
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/886Radar or analogous systems specially adapted for specific applications for alarm systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Wind Motors (AREA)

Abstract

A wind turbine blade fault monitoring method comprises the steps of establishing a wind turbine blade model according to actually detected wind turbine blades; determining the number of sampling points for solving the backward scattering electric field of the wind turbine blade; solving a backscattering electric field of the wind turbine blade; obtaining a radar echo signal of the wind turbine blade according to the backscattering electric field of the wind turbine blade; carrying out short-time Fourier transform according to the radar echo signals of the wind turbine blade to obtain a time-frequency graph of the radar echo signals of the wind turbine blade; monitoring the wind turbine blade which actually runs in real time by using a radar, and acquiring a radar echo of the monitored wind turbine blade; carrying out short-time Fourier transform on the obtained radar echo signals of the monitored wind turbine blade to obtain a time-frequency diagram of the radar echo signals of the monitored wind turbine blade; and comparing the obtained time-frequency graphs of the radar echo signals of the wind motor, and judging whether the blades of the wind motor are in failure. The invention solves the technical problem that the running state of the blades of the wind motor cannot be accurately acquired in real time in the running process of the wind motor.

Description

Wind turbine blade fault monitoring method
Technical Field
The invention belongs to the technical field of wind turbine blade fault monitoring, and particularly relates to a wind turbine blade fault monitoring method.
Background
With the initiation of national energy internet engineering targeting 'transmission of clean energy as a main network frame by using an extra-high voltage power grid and transportation of the clean energy as a leading factor and global interconnection of ubiquitous and strong smart power grids', the energy internet engineering serves as the basis of an energy strategy in China and is also an important component of clean and renewable energy, and the development and utilization of wind power energy are more and more emphasized by the nation. A large number of large and medium-sized wind power plant projects built or regularly planned in the construction are generated in areas rich in wind energy resources in China, the wind power produced in 2017 reaches 1503 ten thousand kW, and the accumulated grid-connected installed capacity reaches 1.64 hundred million kW.
With the construction of large and medium-sized wind power plants, the number of wind motors put into use is increased very rapidly, and it is very important to acquire the operating state of the blades of the wind motors in real time.
Disclosure of Invention
In view of the technical problems in the background art, the method for monitoring the faults of the blades of the wind turbine provided by the invention solves the technical problem that the running state of the blades of the wind turbine cannot be accurately acquired in real time in the running process of the wind turbine.
In order to solve the technical problems, the invention adopts the following technical scheme to realize:
a wind turbine blade fault monitoring method comprises the following steps:
the method comprises the following steps: establishing a wind turbine blade model according to the actual monitored wind turbine blade;
according to the wind turbine blade fault monitoring method provided by the invention, the micro Doppler characteristic of the radar echo when the monitored wind turbine operates normally needs to be obtained, so that a wind turbine model needs to be established according to the actual monitored wind turbine for obtaining the micro Doppler characteristic of the wind turbine blade in the normal operation state.
Step two: and determining the number of sampling points for solving the backward scattering electric field of the wind turbine blade according to the Nyquist sampling theorem.
The "Nyquist sampling theorem" requires a sampling frequency fcSatisfy fc>2fdmax,fdmaxIs the maximum Doppler frequency, f, generated during the operation of the wind turbine bladedmax2f ω R/c; ω is the angular velocity of the wind turbine blade rotation; r is the wind turbine blade length. When sampling frequency fcWhen the electric field is determined, the number of sampling points for solving the backward scattering electric field of the wind turbine blade in a corresponding period is 360fc/ω。
Step three: and solving the backscattering electric field of the wind turbine blade by using a mixed algorithm formed by a physical optics method (PO) and a moment method (MOM).
The hybrid algorithm adopts a triangular surface element to disperse the wind turbine blade, and in PO areas such as the front surface of the wind turbine blade and the front surface of an engine room, induced current generated by radar signals on the triangular surface element on the surface of the wind turbine blade is used as a secondary scattering source to replace the wind turbine blade by a physical optical method (PO), and then the induced current on the surface is integrated to calculate the scattering electric field of the wind turbine blade; and establishing a matrix equation to solve the induced current of the MOM region in the MOM region at the back of the blade, the back of the engine room, the discontinuous part of the surface of the wind turbine blade and some complex edge parts, and further solving the scattering electric field of the MOM region. And finally, obtaining the scattering electric field of the wind turbine blade through vector superposition. And solving the scattering electric field of the wind motor blade in a certain fixed attitude, and combining a quasi-static method to obtain the dynamic scattering electric field of the wind motor blade.
Step four: obtaining radar echo signals of the wind turbine blade according to the solved backward scattering electric field of the wind turbine blade;
the echo at the scattering point is expressed as: s (t) exp (j2 pi ft) · σ (t) exp (-j4 pi r (t)/λ), where f and λ are radar signal frequency and wavelength, respectively; r (t) and sigma (t) are respectively the distance from a scattering point to the radar at the moment t and the scattering coefficient of the scattering point, wherein sigma (t) is generally subjected to homogenization treatment. Constant carrier term exp (j2 pi f)ct) is the carrier of the signal, the useful information of the echo is stored in the baseband signal exp (-j4 π R (t)/λ);
the expression in time-harmonic electromagnetic fields is: es(r,t)=Es(r)cos(2πfct + φ (r)), using a complex number to take the real part of this equation, we can:
Figure GDA0003476858730000021
in the formula
Figure GDA0003476858730000022
Is EsA complex amplitude vector of (r, t), wherein Es(r) is the amplitude of the electric field, and phi (r) is the phase of the electric field; comparing the echo expression of the scattering point with the formula, the baseband signal in the echo corresponds to the complex amplitude vector in the scattering electric field of the target,therefore, the obtained scattering electric field of the target can be converted into an echo signal of the target.
Step five: carrying out short-time Fourier transform (STFT) on the obtained radar echo signal of the wind turbine blade, and obtaining a time-frequency diagram of the radar echo signal of the wind turbine blade;
and (2) to acquire a Time-frequency diagram of the radar echo signals of the wind turbine blade, processing the radar echo signals of the wind turbine blade obtained in the step three by adopting Short Time Fourier Transform (STFT), wherein the calculation formula is as follows:
Figure GDA0003476858730000023
wherein x (k) is a radar echo signal of the wind turbine blade; n is the maximum sampling point number; k. m and n are 0,1,2,3L N-1; Δ t is the sampling interval of the time variable; ω (t) is a window function; the superscript "+" indicates the complex conjugate.
Step six: and monitoring the wind motor blade in actual operation in real time by using a radar, and acquiring a radar echo of the monitored wind motor blade. And carrying out short-time Fourier transform (STFT) on the obtained radar echo signals of the monitored wind turbine blade to obtain a time-frequency diagram of the radar echo signals of the monitored wind turbine blade.
Step seven: comparing the time-frequency graphs of the radar echo signals of the wind motor obtained in the fifth step with the time-frequency graphs of the radar echo signals of the wind motor obtained in the sixth step, and judging whether the blades of the wind motor break down or not;
when the wind motor runs, because the blades of the wind motor rotate periodically, a time-frequency diagram of a radar echo signal of the wind motor presents a unique micro Doppler characteristic. And fifthly, obtaining an echo time-frequency diagram of the blade of the wind turbine in normal operation, and obtaining an echo time-frequency diagram of the blade of the monitored wind turbine in the sixth step. Comparing the two signals, and if the micro Doppler characteristics of the two signals are the same, indicating that the blade of the monitored wind turbine is not in fault; if the micro Doppler characteristics of the two are different, the monitored wind turbine blade is judged to be in fault.
Step eight: after the operation state of the wind outlet motor blade is judged according to the seventh step, if the wind motor blade fails, an alarm is sent out in time; and if the wind motor blade normally runs, continuing to monitor the running state of the wind motor blade.
This patent can reach following beneficial effect:
1. the method solves the problem that the running state of the blades of the wind motor cannot be accurately acquired in real time in the running process of the wind motor, and can quickly and accurately judge the running state of the blades of the wind motor by comparing the micro Doppler characteristics of the echoes of the blades of the wind motor.
2. The invention monitors the running state of the wind turbine blade in real time, can send out an alarm in time when the blade is in fault, transmits the information to workers, reminds the workers to process in time and reduces loss.
Drawings
The invention is further illustrated by the following examples in conjunction with the accompanying drawings:
FIG. 1 is a flow chart of a wind turbine fault monitoring method of the present invention;
FIG. 2 is a wind turbine blade model constructed according to the present invention from a Vestas-V82 wind turbine;
FIG. 3 is a time domain diagram of the blade echo of the Vestas-V82 wind turbine of the present invention during normal operation;
FIG. 4 is a time-frequency diagram of the blade echo of the Vestas-V82 wind turbine in normal operation;
FIG. 5 is a wind turbine blade model after a blade fracture fault occurs in a Vestas-V82 wind turbine blade simulated by the method;
FIG. 6 is a time domain diagram of blade echoes of a Vestas-V82 wind turbine blade after the blade is broken and then fails;
FIG. 7 is a time-frequency diagram of blade echoes of a Vestas-V82 wind turbine blade after the blade is broken and then fails.
Detailed Description
The preferable scheme is as shown in fig. 1 to 7, and the wind motor fault monitoring method is used for judging the operation state of the wind motor blade based on the micro doppler characteristics of the blade echo signal of the wind motor during operation. The method comprises the following steps:
the method comprises the following steps: establishing a wind turbine blade model according to the actual monitored wind turbine blade;
according to the wind turbine blade fault monitoring method provided by the invention, the micro Doppler characteristic of the radar echo when the monitored wind turbine operates normally needs to be obtained, so that a wind turbine model needs to be established according to the actual monitored wind turbine for obtaining the micro Doppler characteristic of the wind turbine blade in the normal operation state.
Step two: and determining the number of sampling points for solving the backward scattering electric field of the wind turbine blade according to the Nyquist sampling theorem.
The "Nyquist sampling theorem" requires a sampling frequency fcSatisfy fc>2fdmax,fdmaxIs the maximum Doppler frequency, f, generated during the operation of the wind turbine bladedmax2f ω R/c; ω is the angular velocity of the wind turbine blade rotation; r is the wind turbine blade length. When sampling frequency fcWhen the electric field is determined, the number of sampling points for solving the backward scattering electric field of the wind turbine blade in a corresponding period is 360fc/ω。
Step three: and solving the backscattering electric field of the wind turbine blade by using a mixed algorithm formed by a physical optics method (PO) and a moment method (MOM).
The hybrid algorithm adopts a triangular surface element to disperse the wind turbine blade, and in PO areas such as the front surface of the wind turbine blade and the front surface of an engine room, induced current generated by radar signals on the triangular surface element on the surface of the wind turbine blade is used as a secondary scattering source to replace the wind turbine blade by a physical optical method (PO), and then the induced current on the surface is integrated to calculate the scattering electric field of the wind turbine blade; and establishing a matrix equation to solve the induced current of the MOM region in the MOM region at the back of the blade, the back of the engine room, the discontinuous part of the surface of the wind turbine blade and some complex edge parts, and further solving the scattering electric field of the MOM region. And finally, obtaining the scattering electric field of the wind turbine blade through vector superposition. And solving the scattering electric field of the wind motor blade in a certain fixed attitude, and combining a quasi-static method to obtain the dynamic scattering electric field of the wind motor blade.
Step four: obtaining radar echo signals of the wind turbine blade according to the solved backward scattering electric field of the wind turbine blade;
the echo at the scattering point is expressed as: s (t) ═ exp (j2 π ft). σ (t) exp (-j4 π R (t)/λ), itF and lambda are radar signal frequency and wavelength respectively; r (t) and sigma (t) are respectively the distance from a scattering point to the radar at the moment t and the scattering coefficient of the scattering point, wherein sigma (t) is generally subjected to homogenization treatment. Constant carrier term exp (j2 pi f)ct) is the carrier of the signal, the useful information of the echo is stored in the baseband signal exp (-j4 π R (t)/λ);
the expression in time-harmonic electromagnetic fields is: es(r,t)=Es(r)cos(2πfct + φ (r)), using a complex number to take the real part of this equation, we can:
Figure GDA0003476858730000051
in the formula
Figure GDA0003476858730000052
Is EsA complex amplitude vector of (r, t), wherein Es(r) is the amplitude of the electric field, and phi (r) is the phase of the electric field; comparing the echo expression of the scattering point with the formula, the baseband signal in the echo corresponds to the complex amplitude vector in the scattering electric field of the target, so that the obtained scattering electric field of the target can be converted into the echo signal of the target.
Step five: carrying out short-time Fourier transform (STFT) on the obtained radar echo signal of the wind turbine blade, and obtaining a time-frequency diagram of the radar echo signal of the wind turbine blade;
and (2) to acquire a Time-frequency diagram of the radar echo signals of the wind turbine blade, processing the radar echo signals of the wind turbine blade obtained in the step three by adopting Short Time Fourier Transform (STFT), wherein the calculation formula is as follows:
Figure GDA0003476858730000053
wherein x (k) is a radar echo signal of the wind turbine blade; n is the maximum sampling point number; k. m and n are 0,1,2,3L N-1; Δ t is the sampling interval of the time variable; ω (t) is a window function; the superscript "+" indicates the complex conjugate.
Step six: and monitoring the wind motor blade in actual operation in real time by using a radar, and acquiring a radar echo of the monitored wind motor blade. And carrying out short-time Fourier transform (STFT) on the obtained radar echo signals of the monitored wind turbine blade to obtain a time-frequency diagram of the radar echo signals of the monitored wind turbine blade.
Step seven: comparing the time-frequency graphs of the radar echo signals of the wind motor obtained in the fifth step with the time-frequency graphs of the radar echo signals of the wind motor obtained in the sixth step, and judging whether the blades of the wind motor break down or not;
when the wind motor runs, because the blades of the wind motor rotate periodically, a time-frequency diagram of a radar echo signal of the wind motor presents a unique micro Doppler characteristic. And fifthly, obtaining an echo time-frequency diagram of the blade of the wind turbine in normal operation, and obtaining an echo time-frequency diagram of the blade of the monitored wind turbine in the sixth step. Comparing the two signals, and if the micro Doppler characteristics of the two signals are the same, indicating that the blade of the monitored wind turbine is not in fault; if the micro Doppler characteristics of the two are different, the monitored wind turbine blade is judged to be in fault.
Step eight: after the operation state of the wind outlet motor blade is judged according to the seventh step, if the wind motor blade fails, an alarm is sent out in time; and if the wind motor blade normally runs, continuing to monitor the running state of the wind motor blade.
The following description will be made in detail by taking a Vestas-V82 wind motor as an example:
a simulation model of a wind turbine blade is established according to a Vestas-V82 wind turbine, the wind turbine blade model is shown in figure 2, the number of blades is 3, the length of the blade is 44 meters, and it is assumed that a radar signal is incident on a yoz plane, and an incident angle theta and the rotation direction of the blade are shown in figure 2. The transmitting frequency of the radar signal is 1GHZ, the number of sampling points for solving a backward scattering electric field of the wind turbine blade is 6000 according to the Nyquist sampling theorem, the polarization direction of electromagnetic waves is vertical polarization, and a time domain graph and a time frequency graph of blade echoes of the Vestas-V82 wind turbine during normal operation are shown in FIGS. 3-4.
The micro Doppler characteristic of the wind motor blade echo mainly comprises a zero frequency band, a sinusoidal envelope curve and time-frequency flicker. The "zero band" is due to the nacelle being relatively stationary with the radar at the wind turbine blade axis. The blade tips are equivalent to circular motion, and the "sinusoidal envelope" corresponds to the echoes of the blade tips. When the radar sight is perpendicular to the wind turbine blade, the blade is completely irradiated by a radar signal, the echo reaches the peak value, the time-frequency domain waveform generates a flickering phenomenon, the blade presents positive time-frequency flickering when being close to the radar, and the blade presents negative time-frequency flickering when being far away from the radar. The blade is completely irradiated, and when the echo reaches the peak value, the time domain waveform can generate a wave peak, so that the time when the time-frequency flicker occurs corresponds to the time when the wave peak of the time domain waveform occurs. When positive Doppler flicker occurs, the radar irradiates the front edge of the wind turbine blade, the front edge structure is thick, the energy of an echo is large, the energy of the positive Doppler flicker is higher than that of the negative Doppler flicker, and the complex shape of the wind turbine blade enables the time-frequency domain waveform to have the phenomenon of time-frequency flicker bending.
FIG. 5 is a wind turbine blade model after blade breakage failure of a simulated Vestas-V82 wind turbine blade, and a time domain graph and a time frequency graph of blade echo of the Vestas-V82 wind turbine blade when blade breakage failure occurs are shown in FIGS. 6-7.
When the radar sight is perpendicular to the blades of the wind motor, the time-frequency domain waveform has a flickering phenomenon, and the Vestas-V82 wind motor has three blades, so that the time-frequency domain waveform has three groups of time-frequency flickering, and the maximum Doppler frequency is equal. However, when a certain blade breaks down, the time-frequency domain waveform of the radar echo changes to a certain extent, so that the blade of the wind turbine can be judged to break down.
In conclusion, the wind motor fault monitoring method can effectively solve the technical problem that the running state of the wind motor blade cannot be accurately acquired in real time in the running process of the wind motor, and can remind workers to process the wind motor blade in time when the wind motor blade fails, so that loss is reduced.
The above-described embodiments are merely preferred embodiments of the present invention, and should not be construed as limiting the present invention, and the scope of the present invention is defined by the claims, and equivalents including technical features described in the claims. I.e., equivalent alterations and modifications within the scope hereof, are also intended to be within the scope of the invention.

Claims (10)

1. A wind turbine blade fault monitoring method is characterized by comprising the following steps:
the method comprises the following steps: establishing a wind motor blade model according to the actually detected wind motor blade;
step two: determining the number of sampling points for solving the backward scattering electric field of the wind turbine blade according to the Nyquist sampling theorem;
step three: solving a backscattering electric field of the wind turbine blade by using a hybrid algorithm formed by a physical optical method and a moment method;
step four: obtaining radar echo signals of the wind turbine blade according to the solved backward scattering electric field of the wind turbine blade;
step five: carrying out short-time Fourier transform on the obtained radar echo signal of the wind turbine blade to obtain a time-frequency graph of the radar echo signal of the wind turbine blade;
step six: monitoring the wind turbine blade which actually runs in real time by using a radar, and acquiring a radar echo of the monitored wind turbine blade; carrying out short-time Fourier transform on the obtained radar echo signals of the monitored wind turbine blade to obtain a time-frequency diagram of the radar echo signals of the monitored wind turbine blade;
step seven: comparing the time-frequency graphs of the radar echo signals of the wind motor obtained in the fifth step with the time-frequency graphs of the radar echo signals of the wind motor obtained in the sixth step, and judging whether the blades of the wind motor break down or not;
step eight: after the running state of the wind outlet motor blade is judged, if the wind motor blade breaks down, an alarm is sent out in time; and if the wind motor blade normally runs, continuing to monitor the running state of the wind motor blade.
2. The wind turbine blade fault monitoring method according to claim 1, characterized in that: the step of obtaining the radar echo signal of the wind motor blade by four-way passing through the backward scattering electric field of the wind motor blade is that a baseband signal in the echo corresponds to a complex amplitude vector of a target scattering electric field.
3. The wind turbine blade fault monitoring method according to claim 2, characterized in that: fifthly, obtaining an echo time-frequency diagram of the wind turbine blade in normal operation, and sixthly obtaining an echo time-frequency diagram of the monitored wind turbine blade; when the wind turbine runs, because the blades of the wind turbine rotate periodically, a time-frequency diagram of a radar echo signal of the wind turbine can present a unique micro Doppler characteristic; and fifthly, acquiring the micro Doppler characteristics of the radar echo signals of the wind turbine blade, namely performing short-time Fourier transform on the radar echo signals of the wind turbine blade obtained in the fourth step.
4. The wind turbine blade fault monitoring method according to claim 3, characterized in that: and sixthly, obtaining the radar echo and the micro Doppler characteristic of the monitored wind turbine blade through radar real-time monitoring.
5. The wind turbine blade fault monitoring method according to claim 1, characterized in that: and judging whether the wind motor blade has a fault or not by comparing whether the blade micro Doppler characteristics obtained in the fifth step and the sixth step are the same or not.
6. The wind turbine blade fault monitoring method according to claim 1, characterized in that: after the running state of the wind outlet motor blade is judged, if the wind motor blade breaks down, an alarm is sent out in time; and if the wind motor blade normally runs, continuing to monitor the running state of the wind motor blade.
7. The wind turbine blade fault monitoring method according to claim 1, characterized in that: in the second step, the number of sampling points for solving the backward scattering electric field of the wind turbine blade is determined according to the Nquist sampling theorem, and the method comprises the following steps:
the "Nquist sampling theorem" requires a sampling frequency fcSatisfy fc>2fdmax,fdmaxIs the maximum Doppler frequency, f, generated during the operation of the wind turbine bladedmax=2fωR/c;
ω is the angular velocity of the wind turbine blade rotation; r is the wind turbine blade length; when sampling frequency fcWhen determined, thenThe number of sampling points for solving the backward scattering electric field of the wind turbine blade in a corresponding period is 360fc/ω;
f is the radar signal frequency.
8. The wind turbine blade fault monitoring method according to claim 7, wherein: in the third step, a hybrid algorithm formed by a physical optical method and a moment method is used for solving the backscattering electric field of the wind turbine blade, and the method comprises the following steps:
in the hybrid algorithm, a triangular surface element is adopted to disperse the wind turbine blade, induced current generated by radar signals on the triangular surface element on the surface of the wind turbine blade is used as a secondary scattering source to replace the wind turbine blade in PO areas such as the front surface of the wind turbine blade and the front surface of an engine room by a physical optical method, and then the surface induced current is integrated to calculate the scattering electric field of the wind turbine blade; establishing a matrix equation in MOM areas such as the back of a blade, the back of an engine room, the surface discontinuity of a wind turbine blade and some complex edge parts to solve induced currents of the MOM areas so as to obtain a scattering electric field of the MOM areas; finally, a scattering electric field of the wind motor blade can be obtained through vector superposition; and solving the scattering electric field of the wind motor blade in a certain fixed attitude, and combining a quasi-static method to obtain the dynamic scattering electric field of the wind motor blade.
9. The wind turbine blade fault monitoring method according to claim 8, wherein: obtaining radar echo signals of the wind turbine blade according to the solved backward scattering electric field of the wind turbine blade;
the echo at the scattering point is expressed as: s (t) exp (j2 pi ft) · σ (t) exp (-j4 pi r (t)/λ), where f and λ are radar signal frequency and wavelength, respectively; r (t) and sigma (t) are respectively the distance from a scattering point to the radar at the moment t and the scattering coefficient of the scattering point, wherein sigma (t) is generally subjected to homogenization treatment; constant carrier term exp (j2 pi f)ct) is the carrier of the signal, the useful information of the echo is stored in the baseband signal exp (-j4 π R (t)/λ);
the expression in time-harmonic electromagnetic fields is: es(r,t)=Es(r)cos(2πfct + φ (r)), using a complex number to take the real part of this equation, we can:
Figure FDA0003476858720000031
in the formula
Figure FDA0003476858720000032
Is EsA complex amplitude vector of (r, t), wherein Es(r) is the amplitude of the electric field, and phi (r) is the phase of the electric field; comparing the echo expression of the scattering point with the formula, the baseband signal in the echo corresponds to the complex amplitude vector in the scattering electric field of the target, so that the obtained scattering electric field of the target can be converted into the echo signal of the target.
10. The wind turbine blade fault monitoring method of claim 9, wherein: carrying out short-time Fourier transform on the obtained radar echo signal of the wind turbine blade to obtain a time-frequency graph of the radar echo signal of the wind turbine blade;
and (3) acquiring a time-frequency diagram of the radar echo signals of the wind turbine blade, wherein the radar echo signals of the wind turbine blade acquired in the third short-time Fourier transform processing step need to be processed, and the calculation formula is as follows:
Figure FDA0003476858720000033
wherein x (k) is a radar echo signal of the wind turbine blade; n is the maximum sampling point number; k. m and N are 0,1,2,3 … N-1; Δ t is the sampling interval of the time variable; ω (t) is a window function; the superscript "+" indicates the complex conjugate.
CN202011269371.1A 2020-11-13 2020-11-13 Wind turbine blade fault monitoring method Active CN112629838B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011269371.1A CN112629838B (en) 2020-11-13 2020-11-13 Wind turbine blade fault monitoring method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011269371.1A CN112629838B (en) 2020-11-13 2020-11-13 Wind turbine blade fault monitoring method

Publications (2)

Publication Number Publication Date
CN112629838A CN112629838A (en) 2021-04-09
CN112629838B true CN112629838B (en) 2022-03-08

Family

ID=75303244

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011269371.1A Active CN112629838B (en) 2020-11-13 2020-11-13 Wind turbine blade fault monitoring method

Country Status (1)

Country Link
CN (1) CN112629838B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114030490B (en) * 2022-01-12 2022-04-26 深圳市永达电子信息股份有限公司 Collision determination method in operation of movable platform door and computer storage medium

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3806929A (en) * 1971-06-24 1974-04-23 Us Navy Method for the detection of radar targets
CN106291482A (en) * 2016-09-29 2017-01-04 三峡大学 A kind of wind turbine blade radar echo signal Doppler frequency spectrum method for solving
CN108256488A (en) * 2018-01-19 2018-07-06 中国人民解放军陆军装甲兵学院 A kind of radar target identification method based on micro-Doppler feature extraction and deep learning
CN108708871A (en) * 2018-04-19 2018-10-26 浙江大学 A kind of axial flow blower vibration signal subconstiuent extracting method based on inverse Short Time Fourier Transform
CN109959930A (en) * 2017-12-22 2019-07-02 英飞凌科技股份有限公司 Use the system and method for millimetre-wave radar sensor monitoring of structures object
CN110658519A (en) * 2019-09-18 2020-01-07 武汉大学 Wind turbine generator blade state monitoring method based on external radiation source radar
EP3663732A1 (en) * 2018-12-05 2020-06-10 Tata Consultancy Services Limited Method and system for monitoring machine health using radar based segregation for induced machine vibrations

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3806929A (en) * 1971-06-24 1974-04-23 Us Navy Method for the detection of radar targets
CN106291482A (en) * 2016-09-29 2017-01-04 三峡大学 A kind of wind turbine blade radar echo signal Doppler frequency spectrum method for solving
CN109959930A (en) * 2017-12-22 2019-07-02 英飞凌科技股份有限公司 Use the system and method for millimetre-wave radar sensor monitoring of structures object
CN108256488A (en) * 2018-01-19 2018-07-06 中国人民解放军陆军装甲兵学院 A kind of radar target identification method based on micro-Doppler feature extraction and deep learning
CN108708871A (en) * 2018-04-19 2018-10-26 浙江大学 A kind of axial flow blower vibration signal subconstiuent extracting method based on inverse Short Time Fourier Transform
EP3663732A1 (en) * 2018-12-05 2020-06-10 Tata Consultancy Services Limited Method and system for monitoring machine health using radar based segregation for induced machine vibrations
CN110658519A (en) * 2019-09-18 2020-01-07 武汉大学 Wind turbine generator blade state monitoring method based on external radiation source radar

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于动态RCS的风电机叶片多普勒特性;唐波 等;《高电压技术》;20171031;第43卷(第10期);第3435-3442页 *
直升机旋翼桨叶外形对雷达特征信号的影响;蒋相闻 等;《航空学报》;20141125;第35卷(第11期);第3123-3136页 *

Also Published As

Publication number Publication date
CN112629838A (en) 2021-04-09

Similar Documents

Publication Publication Date Title
Li et al. Experimental investigation of solidity and other characteristics on dual vertical axis wind turbines in an urban environment
Naqvi et al. Investigation of Doppler features from wind turbine scattering
Doolan et al. WIND TURBINE NOISE MECHANISMS AND SOME CONCEPTS FOR ITS CONTROL.
CN102954857B (en) Vane unbalance fault diagnosis method of wind turbine generator set based on current signal
Oerlemans et al. Experimental demonstration of wind turbine noise reduction through optimized airfoil shape and trailing-edge serrations
CN103216383B (en) Control system and control method for wind generating set
Laratro et al. A discussion of wind turbine interaction and stall contributions to wind farm noise
Lee et al. Numerical and experimental study of aerodynamic noise by a small wind turbine
CN112629838B (en) Wind turbine blade fault monitoring method
Zhao et al. Remote structural health monitoring for industrial wind turbines using short-range Doppler radar
Xie et al. An integration fault detection method using stator voltage for marine current turbines
Burlibaşa et al. Rotationally sampled spectrum approach for simulation of wind speed turbulence in large wind turbines
Hu et al. An experimental investigation on the effects of turbine rotation directions on the wake interference of wind turbines
Kurniawati et al. Experimental investigation on performance of crossflow wind turbine as effect of blades number
CN112018784A (en) Subsynchronous resonance tracing method based on synchronous phasor measurement data
CN113572156B (en) Power spectral density-based power system equivalent inertia evaluation method
CN108278184A (en) Impeller of wind turbine set imbalance monitoring method based on empirical mode decomposition
Wei et al. A VMD denoising-based imbalance fault detection method for marine current Turbine
Viterna The NASA-LeRC wind turbine sound prediction code
Bakar et al. Experimental Study of a Hybrid Turbine for Hydrokinetic Applications on Small Rivers in Malaysia
Pearson et al. Experimental characterization of vertical-axis wind turbine noise
Choy et al. Real-time hardware simulator for grid-tied PMSG wind power system
Toshio et al. Failure forecast diagnosis of small wind turbine using acoustic emission sensor
CN110608133A (en) Offshore wind power generation control system and method
CN113357082B (en) Wind turbine generator protection method

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
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20210409

Assignee: Yichang Yizhixing Technology Co.,Ltd.

Assignor: CHINA THREE GORGES University

Contract record no.: X2023980034895

Denomination of invention: Fault Monitoring Method for Wind Turbine Blades

Granted publication date: 20220308

License type: Common License

Record date: 20230426

EE01 Entry into force of recordation of patent licensing contract