CN111665047A

CN111665047A - Method, device, equipment and medium for identifying vibration signal of wind generating set

Info

Publication number: CN111665047A
Application number: CN201910173392.4A
Authority: CN
Inventors: 叶月光
Original assignee: Beijing Goldwind Smart Energy Service Co Ltd
Current assignee: Beijing Goldwind Smart Energy Service Co Ltd
Priority date: 2019-03-07
Filing date: 2019-03-07
Publication date: 2020-09-15
Anticipated expiration: 2039-03-07
Also published as: CN111665047B

Abstract

The embodiment of the application provides a method, a device, equipment and a medium for identifying vibration signals of a wind generating set. The method comprises the following steps: acquiring real-time operation data of the wind generating set; the real-time operating data includes: the method comprises the following steps of (1) time tagging, real-time vibration data of a cabin, real-time rotating speed of a generator and real-time wind speed; the real-time vibration data comprises first direction vibration data and second direction vibration data; the first direction is the axis direction of the engine room, and the second direction is vertical to the first direction; determining a standard deviation of vibration data in a first direction and a standard deviation of vibration data in a second direction of the engine room according to the real-time operation data; and identifying a signal representing the non-stationary vibration of the generator according to the standard deviation of the vibration data in the first direction and the standard deviation of the vibration data in the second direction. The method can realize the identification of the non-stable vibration signal of the wind generating set in the bearing running process through the real-time running data of the wind generating set.

Description

Method, device, equipment and medium for identifying vibration signal of wind generating set

Technical Field

The application relates to the technical field of vibration detection, in particular to a method and a device for identifying vibration signals of a wind generating set, electronic equipment and a computer readable storage medium.

Background

The generator in the wind generating set is an important power generation equipment component in the wind power plant, and in the long-term operation of the actual wind power plant, the bearing of the generator inevitably has faults such as abrasion, cracking and the like, and most of the generator bearings are not provided with a bearing vibration monitoring system, and related monitoring data records do not exist. There is therefore a safety risk for the operation of such critical mechanical rotating parts.

The existing generator bearing problem in the wind power plant is mainly realized by replacing after-accident faults after triggering faults such as regular maintenance, regular bearing clearance test, generator running bearing blocking faults and the like. The preventive scheme adopted in the industry is also to monitor the vibration by additionally adding monitoring equipment, and the scheme increases the equipment operation cost, so that the scheme is not popularized and applied in wind power plants.

More schemes in the abnormal operation vibration monitoring of the direct drive generator bearing in the current industry adopt manual inspection, periodic maintenance and additional installation of a vibration monitoring device to realize external abnormal monitoring of the unit, so that the manual maintenance cost is increased, and real-time monitoring cannot be realized. The vibration fault early warning of the abnormal operation of the generator bearing in the current industry is mainly realized by additionally adding a new bearing vibration monitoring system to the generator set to detect an abnormal vibration signal, the wind generator set needs additional technical improvement, the technical improvement work greatly increases the cost, and the generator bearing is not popularized in engineering application actually.

In conclusion, the prior art has the defects that the vibration of the engine cannot be monitored in real time or additional detection equipment is required.

Disclosure of Invention

The application aims at the defects of the existing mode and provides a method and a device for identifying a vibration signal of a wind generating set, electronic equipment and a computer readable storage medium, which are used for solving the technical problem that the vibration of an engine cannot be monitored in real time or extra detection equipment needs to be added.

In a first aspect, the application provides a method for identifying a vibration signal of a wind generating set, including:

acquiring real-time operation data of the wind generating set; the real-time operating data includes: the method comprises the following steps of (1) time tagging, real-time vibration data of a cabin, real-time rotating speed of a generator and real-time wind speed; the real-time vibration data comprises first direction vibration data and second direction vibration data; the first direction is the axis direction of the engine room, and the second direction is vertical to the first direction;

determining a standard deviation of vibration data in a first direction and a standard deviation of vibration data in a second direction of the engine room according to the real-time operation data;

and identifying a signal representing the non-stationary vibration of the generator according to the standard deviation of the vibration data in the first direction and the standard deviation of the vibration data in the second direction.

In a second aspect, the present application provides an apparatus for recognizing a vibration signal of a wind turbine generator system, including:

the first processing module is used for acquiring real-time operation data of the wind generating set; the real-time operating data includes: the method comprises the following steps of (1) time tagging, real-time vibration data of a cabin, real-time rotating speed of a generator and real-time wind speed; the real-time vibration data comprises first direction vibration data and second direction vibration data; the first direction is the axis direction of the engine room, and the second direction is vertical to the first direction;

the second processing module is used for determining the standard deviation of the vibration data in the first direction and the standard deviation of the vibration data in the second direction of the cabin according to the real-time operation data;

and the third processing module is used for identifying a signal representing the non-stationary vibration of the generator according to the standard deviation of the vibration data in the first direction and the standard deviation of the vibration data in the second direction.

In a third aspect, the present application provides an electronic device, comprising: a processor, a memory, and a bus;

a bus for connecting the processor and the memory;

a memory for storing operating instructions;

and the processor is used for executing the identification method of the vibration signal of the wind generating set in the first aspect of the application by calling the operation instruction.

In a fourth aspect, the present application provides a computer-readable storage medium storing a computer program for executing the method for identifying a wind turbine generator set vibration signal of the first aspect of the present application.

The technical scheme provided by the embodiment of the application at least has the following beneficial effects:

1) acquiring real-time operation data of the wind generating set; determining a standard deviation of vibration data in a first direction and a standard deviation of vibration data in a second direction of the engine room according to the real-time operation data; identifying a signal representing non-stationary vibration of the generator according to the standard deviation of the vibration data in the first direction and the standard deviation of the vibration data in the second direction, wherein the whole process of the identification method in the embodiment of the application is automatically executed by starting in real time; that is to say, the embodiment of the application realizes the real-time identification of the non-stationary vibration signal in the running process of the generator (including the bearing of the generator) of the wind generating set, thereby realizing the early discovery, early processing and early solution of the vibration problem of the generator.

2) The method and the device can monitor the non-steady vibration signal of the running of the generator (including the bearing) of the wind generating set without adding extra detection equipment, reduce the detection cost of detecting the non-steady vibration of the generator, and have wider application range.

Additional aspects and advantages of the present application will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the present application.

Drawings

The foregoing and/or additional aspects and advantages of the present application will become apparent and readily appreciated from the following description of the embodiments, taken in conjunction with the accompanying drawings of which:

fig. 1 is a schematic flow chart of a method for identifying a vibration signal of a wind turbine generator system according to an embodiment of the present application;

FIG. 2 is a schematic flow chart of another method for identifying a vibration signal of a wind turbine generator system according to an embodiment of the present application;

fig. 3 is a schematic flow chart of a further method for identifying a vibration signal of a wind turbine generator system according to an embodiment of the present application;

fig. 4 is a schematic view of fluctuation rates (fluctuation is relatively gentle) of two vibration directions of a normal unit changing with a rotation speed according to an embodiment of the present application;

FIG. 5 is a schematic diagram of the variation rate of the abnormal unit in two vibration directions with the variation of the rotation speed (the variation rate of the unit vibration with the worn 3MW bearing) provided by the embodiment of the present application;

FIG. 6 is a schematic view of a first-order mode of a direct-drive wind turbine generator showing an arch distribution characteristic in a specific rotation speed domain according to an embodiment of the present application;

FIG. 7 is a schematic illustration of a 1.5Mw wind turbine bearing wear characteristic provided by an embodiment of the present application;

FIG. 8 is a schematic illustration of a normal characterization of a 1.5Mw direct drive wind turbine bearing provided by an embodiment of the present application;

fig. 9 is a schematic structural diagram of an identification device for a vibration signal of a wind generating set according to an embodiment of the present application;

fig. 10 is a schematic structural diagram of an electronic device according to an embodiment of the present application.

Detailed Description

Reference will now be made in detail to the present application, examples of which are illustrated in the accompanying drawings, wherein like reference numerals refer to the same or similar parts or parts having the same or similar functions throughout. In addition, if a detailed description of the known art is not necessary for illustrating the features of the present application, it is omitted. The embodiments described below with reference to the drawings are exemplary only for the purpose of explaining the present application and are not to be construed as limiting the present application.

It will be understood by those within the art that, unless otherwise defined, all terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs. It will be further understood that terms, such as those defined in commonly used dictionaries, should be interpreted as having a meaning that is consistent with their meaning in the context of the prior art and will not be interpreted in an idealized or overly formal sense unless expressly so defined herein.

As used herein, the singular forms "a", "an", "the" and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms "comprises" and/or "comprising," when used in this specification, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof. It will be understood that when an element is referred to as being "connected" or "coupled" to another element, it can be directly connected or coupled to the other element or intervening elements may also be present. Further, "connected" or "coupled" as used herein may include wirelessly connected or wirelessly coupled. As used herein, the term "and/or" includes all or any element and all combinations of one or more of the associated listed items.

The basic principle of the embodiments of the present application is first described as follows:

standard Deviation, also commonly known as mean square error in the Chinese context, is the square root of the arithmetic mean of the squared mean squared Deviation, expressed as σ_XAnd (4) showing. The standard deviation is the arithmetic square root of the variance. The standard deviation can reflect the degree of dispersion of a data set. The standard deviation is not necessarily the same for two sets of data with the same mean. The standard deviation is applied to investment and can be used as an index for measuring the stability of the return. The larger the standard deviation value, the more the return is from the past average value, and the more unstable the return is, the higher the risk. Conversely, a smaller standard deviation value indicates a more stable response and a lower risk. The standard deviation is affected by the extreme values. Smaller standard deviations indicate more data aggregation; the larger the standard deviation, the more discrete the data. The degree of dispersion of a set of data is reflected by the standard deviation, which is also an important index representing the accuracy. The standard deviation is calculated as follows:

in the formula, X_iIs a real number, the mean (arithmetic mean) is μ_XStandard deviation of σ_X。

Kurtosis (kurtosis), also known as the kurtosis coefficient, is a characteristic number that characterizes the height of the peak of a probability density distribution curve at the mean, and reflects the sharpness of the peak. The kurtosis of a sample is a statistic compared to a normal distribution, and generally if the kurtosis is greater than three, the shape of the peak is sharper than the normal distribution peak, and vice versa. Kurtosis is the kurtosis for measuring the probability distribution of real random variables, and the calculation formula of the kurtosis is as follows:

in the formula, X_iFor real numbers, n is the sample size and StDev is the standard deviation estimated based on the samples.

First-order mode: the vibration mode of the object is called a first-order mode or a main mode when the excitation frequency of the external force is equal to the natural frequency of the object. The mode is the natural vibration characteristic of the mechanical structure, and each order mode corresponds to a specific natural frequency, a mode shape and a damping ratio. The concept of mode is generally used in the field of vibration, and can be understood as a vibration state, each object has its own natural frequency, and under the excitation action of an external force, the object can show different vibration characteristics. The first-order mode appears when the excitation frequency of the external force is equal to the natural frequency (first order) of the object, and the vibration form of the object is called a first-order mode or a main mode; the second-order mode appears when the excitation frequency of the external force is twice of the natural frequency (second order) of the object, and the vibration form of the object is called a second-order mode shape at the moment, and so on.

Correlation coefficient: the correlation coefficient is a measure of the degree of linear correlation between the study variables. The correlation coefficient has several definitions according to different research objects.

The simple correlation coefficient, also called correlation coefficient or linear correlation coefficient, is generally denoted by the letter r and is used to measure the linear relationship between two variables. The correlation coefficient is calculated as follows:

in the formula, Cov (X, Y) is the covariance of X and Y, Var [ X ] is the variance of X, and Var [ Y ] is the variance of Y.

The inventor of the application finds that the unevenness of the vibration of the unit in the power generation process can be represented by using the standard deviation distribution characteristics of real-time vibration data in the direction of the cabin acceleration X, Y of the unit at a specific rotating speed. The method mainly comprises the steps of representing a specific small rotating speed area and a large rotating speed area, wherein the small rotating speed area represents the vibration distribution which is in an arched shape and appears in a specific part of rotating speed area of a generator in the speed increasing process, the large rotating speed area represents the vibration distribution which is in a linear increasing shape and appears in the specific part of rotating speed area of a generator bearing in the speed increasing process, optionally, the rotating speed in the rotating speed increasing process is 8-17 rpm, the rotating speed which is increased from 8rpm to 17rpm is the specific part of rotating speed, and rpm (revolutions per minute) represents the rotating times per minute of equipment. The description will be made in two cases according to different characteristic expressions.

1) The characteristic of the arc vibration instability in the generating process of the unit is mainly the vibration characteristic aiming at the specific first-order bending torque.

Specifically, by a structural mode analysis method, the vibration characteristics of each order mode of the mechanical structure in a susceptible frequency range and the vibration response results of the mechanical structure in the frequency range and under the excitation action of various internal or external vibration sources can be obtained. The stator and the rotor of the permanent magnet direct drive generator are of a disc structure, and the modal vibration type characteristics of the stator and the rotor have certain similarity. The low frequency modes are first-order bending, first-order axial oscillation, second-order bending and first-order torsion. The frequency is the modal frequency of the stator and rotor structure, and the frequency can be coupled and changed under the action of electromagnetic force and impeller load. Due to the influence of uneven air gaps and impeller loads of the generator, the generator is easy to generate bending moment load excitation around a central shaft, and current fluctuation or torque unbalance is easy to generate excitation in a torsional direction. Therefore, the first-order bending and first-order torsion modes of the stator and the rotor are easy to cause forced vibration and resonance under electromagnetic excitation, and different impeller configurations and loads have the effect of inhibiting or amplifying the vibration. In the result of the acceleration frequency domain analysis, the vibration near 8.xHz is related to the fundamental excitation and the first order bending mode of the stator and rotor, 8.x represents 8.0, 8.1, etc., x represents a positive integer and zero, and the torsional vibration around 14rpm is related to the first order torsional mode of the rotor. The vibration amplitude is related to electromagnetic excitation, current imbalance and structural damping.

The first-order bending vibration of the generator may cause the cracking problem of a generator stator support, the high-frequency abnormal vibration of a generator bearing and other abnormalities, and due to the influence of current harmonic waves and wind turbulence, the electromagnetic torque of the generator generates a large amount of pulsation, which inevitably causes remarkable influence on the service life of mechanical parts such as a bearing of a transmission shaft system.

2) The linear non-stationarity linearity of the generator bearing during operation mainly shows the problem of bearing abrasion of the large direct-drive wind driven generator.

Specifically, when wear occurs in a bearing of a unit, and the sliding surface inside the bearing is not smooth enough, an unstable phenomenon of vibration inevitably occurs during rotation of the unit. The greater the change in the rotational speed, the more pronounced the jerky behavior. This is because if a peeling pit is formed in the raceway due to wear or peeling, and the peeling pit is located in the load zone of the low-speed rolling bearing, the raceway is locally unloaded and the pipe is locally overloaded at the instant when the raceway is locally unloaded and the center of contact occurs between the roller and the raceway which are in relative circular motion with the raceway, and then returns to normal, and then the stress of the raceway where the peeling pit is located and the vicinity thereof fluctuates intermittently. If this is the case to some extent, the nearby stress, strain and contact stress between the contact surfaces will all change regularly to a large extent. Irregular movement caused by metal friction or deformation of rough surfaces and movement of cracks and swarf will result in the generation of stress waves. The PLC is an electronic system using digital operation in industrial environment, and it uses the instructions designed by user stored in the Programmable memory to implement Logic operation, sequential operation, timing, counting and arithmetic operation and control various types of machinery and processes by digital or analog input and output. In the second-level data sampled randomly in the operation of the unit, the characteristic of the periodic signal is almost completely lost, and the vibration performance of the unit can be described by the turbulence of the vibration.

If the inner ring of the bearing is not coincident with the axial lead, the bearing is easy to generate alternating axial force action once when rotating for one circle, and a certain angle exists between the elevation angle and the horizontal direction of the unit, so the vibration value of the unit in the X direction in the rotating process is more obvious in the rotating speed change process, and if the contact surface of the bearing is worn, and the temperature rise caused by peeling, pockmarks and clearance increase or the gap between the inner sleeve of the bearing and the shaft neck is uneven or scars appear, the running vibration in the rotating plane direction can influence the represented stability.

In the assembly process of a bearing and a main shaft of a generator, the temperature of a bearing ring is not completely recovered to the ambient temperature for assembly, the inner ring of an NJ bearing (the NJ is a bearing type structure code, two letters of N and J are combined together to represent the single flange of the inner ring, and the cylindrical roller bearing with two flanges of the outer ring) is subjected to abnormal motion, so that the position of the bearing assembly is deviated, the outer ring and part of rolling bodies are damaged, the axial displacement of the inner ring of the NG (No Good, bad and abnormal) bearing is found after disassembly, and a large gap exists between the bearing and a retainer.

The retainer ring of the NJ bearing inner ring raceway part of the generator is damaged, the outer ring rolling element retainer is riveted and broken, and the rolling elements are scattered after disassembly. The retainer is loosened, influences are generated on the rolling body, and when a large-area retainer rivet joint is broken, the retainer generates friction on the rolling body, the bearing is scratched, abrasion of the rolling body is aggravated, and a bearing gap is enlarged.

The following describes the technical solutions of the present application and how to solve the above technical problems with specific embodiments. The following several specific embodiments may be combined with each other, and details of the same or similar concepts or processes may not be repeated in some embodiments. Embodiments of the present application will be described below with reference to the accompanying drawings.

Example one

The embodiment of the application provides a method for identifying a vibration signal of a wind generating set, a flow schematic diagram of the method is shown in fig. 1, and the method comprises the following steps:

s101, acquiring real-time operation data of the wind generating set; the real-time operating data includes: the method comprises the following steps of (1) time tagging, real-time vibration data of a cabin, real-time rotating speed of a generator and real-time wind speed; the real-time vibration data comprises first direction vibration data and second direction vibration data; the first direction is an axial direction of the nacelle, and the second direction is perpendicular to the first direction.

And S102, determining the standard deviation of the vibration data of the first direction of the cabin and the standard deviation of the vibration data of the Y direction according to the real-time operation data.

Optionally, before determining the standard deviation of the first direction vibration data and the standard deviation of the second direction vibration data of the nacelle, the method further includes: judging whether the quantity of the real-time vibration data meets a first threshold value and whether the real-time vibration data meets a preset data quality condition; when the quantity of the real-time vibration data meets a first threshold value and the real-time vibration data meets a data quality condition, screening the real-time vibration data to obtain screened real-time vibration data; and determining the standard deviation of the vibration data in the first direction and the standard deviation of the vibration data in the second direction based on the rotating speed according to the screened real-time vibration data and the real-time rotating speed.

Optionally, the determining whether the real-time vibration data meets a preset data quality condition includes: and judging whether the quantity of the real-time vibration data of which the absolute value of the vibration data in the second direction is greater than 1 is smaller than a second threshold value or not, wherein the vibration data comprises acceleration.

Optionally, the screening the real-time vibration data to obtain the screened real-time vibration data includes: according to the first direction vibration data and the second direction vibration data in the real-time vibration data, the real-time vibration data meeting the following first data retention conditions are determined and serve as the real-time vibration data after the first screening: the absolute value of the first-direction vibration data and the absolute value of the second-direction vibration data are both smaller than the high boundary value, and the absolute value of the second-direction vibration data is larger than the low boundary value or the absolute value of the first-direction vibration data is larger than the low boundary value.

Optionally, the screening the real-time vibration data to obtain the screened real-time vibration data, further includes: determining real-time vibration data meeting the following second data retention conditions from the real-time vibration data after the first screening as the real-time vibration data after the second screening: the number of unique values of the first direction vibration data or the second direction vibration data is not less than a third threshold value, the maximum value of the real-time wind speed unique value is at least greater than a cut-in operation wind speed point, the number of real-time wind speeds greater than the cut-in operation wind speed point is at least greater than a fourth threshold value or the duration is not less than a first time period, and the real-time operation data of real-time wind speeds greater than the threshold wind speed point is not less than a fifth threshold value.

And S103, identifying a signal representing the non-stationary vibration of the generator according to the standard deviation of the vibration data in the first direction and the standard deviation of the vibration data in the second direction.

Optionally, identifying a signal indicative of non-stationary vibration of the generator based on the standard deviation of the first direction vibration data and the standard deviation of the second direction vibration data comprises: determining a first rotating speed interval according to the standard deviation of the vibration data in the first direction and the standard deviation of the vibration data in the second direction; determining a second rotating speed interval according to the first rotating speed interval, the rotating speed center position and the rotating speed center deviation position; judging whether a first arch condition and a second arch condition are met or not according to the standard deviation of the second-direction vibration data in the first rotating speed interval, the standard deviation of the second-direction vibration data in the second rotating speed interval, the rotating speed center position, the rotating speed center deviation position and the vibration standard deviation threshold parameter; and when the first and second camber conditions are met, determining real-time vibration data meeting the first and second camber conditions as a camber signal representing non-stationary vibration of the generator.

Optionally, the determining whether the first camber condition and the second camber condition are satisfied according to the standard deviation of the second-direction vibration data in the first rotation speed interval, the standard deviation of the second-direction vibration data in the second rotation speed interval, the rotation speed center position, the rotation speed center deviation position, and the vibration standard deviation threshold parameter includes: determining that a first dome condition is satisfied when: the coefficient of variation of the standard deviation of the second-direction vibration data of the second rotation speed interval is greater than or equal to a sixth threshold value, the polar distance value of the standard deviation of the second-direction vibration data of the second rotation speed interval is greater than or equal to a seventh threshold value, the number of the standard deviation of the second-direction vibration data of the second rotation speed interval is greater than or equal to a vibration standard deviation threshold parameter is greater than or equal to an eighth threshold value, the maximum value of the standard deviation data of the second-direction vibration of the second rotation speed interval is greater than or equal to a ninth threshold value, the number of the standard deviation of the second-direction vibration data of the first rotation speed interval is greater than or equal to a vibration standard deviation threshold parameter is greater than a tenth threshold value, the absolute value of the kurtosis of the standard deviation of the second-direction vibration data of the second rotation speed interval is greater than an eleventh threshold value, and the position where the standard deviation of the second-direction vibration data of the first rotation speed interval, and the position where the standard deviation of the second direction vibration data in the first rotating speed interval is the maximum value is smaller than the second position value; carrying out differential calculation on the standard deviation of the vibration data in the second direction in the first rotating speed interval to obtain a differential calculation value; taking the maximum value in the absolute values of the difference calculation values as the maximum variation amplitude of the operation of the adjacent points of the arch distribution characteristics; and when the maximum variation amplitude of the operation of the adjacent points of the arch-shaped distribution characteristics is smaller than a twelfth threshold value, or the ratio of the maximum variation amplitude of the operation of the adjacent points of the arch-shaped distribution characteristics to the variation coefficient is smaller than or equal to a thirteenth threshold value, determining that the second arch condition is satisfied.

Optionally, identifying a signal indicative of non-stationary vibration of the generator based on the standard deviation of the first direction vibration data and the standard deviation of the second direction vibration data comprises: determining a first rotating speed interval according to the standard deviation of the vibration data in the first direction and the standard deviation of the vibration data in the second direction; determining a relative vibration difference parameter between the first direction and the second direction according to the standard deviation of the first direction vibration data and the standard deviation of the second direction vibration data in the first rotation speed interval, wherein the non-steady vibration in the second direction is linearly increased based on the whole correlation coefficient of the rotation speed when the first direction vibration is stable, the multiple of the rotation speed domain range occupied by the disturbance degree of the second direction vibration relative to the first direction vibration, the rotation speed domain integral value of the non-steady and large second direction vibration and the rotation speed interval proportion of the large second direction vibration; judging whether a first linear condition is met or not according to the relative vibration difference parameter in the first direction and the second direction, the correlation coefficient, the multiple of the range of the rotating speed domain, the integral value of the rotating speed domain and the proportion of the rotating speed interval; and when the first linear condition is met, determining the real-time vibration data meeting the first linear condition as a linear signal representing the non-stationary vibration of the generator.

Optionally, the determining whether the first linear condition is satisfied according to the relative vibration difference parameter in the first direction and the second direction, the correlation coefficient, the multiple of the range of the rotation speed domain, the integral value of the rotation speed domain, and the ratio of the rotation speed interval includes: determining that a first line condition is satisfied when: the difference parameter of the relative vibration of the first direction and the second direction is greater than or equal to a fourteenth threshold value, the first correlation coefficient is greater than a fifteenth threshold value, when the vibration of the first direction is stable, the correlation coefficient of the overall linear increase of the non-stable vibration of the second direction based on the rotating speed comprises the first correlation coefficient, the second correlation coefficient and the third correlation coefficient, the second correlation coefficient is less than or equal to a sixteenth threshold value, the third correlation coefficient is greater than a seventeenth threshold value, the multiple of the rotating speed domain range occupied by the disturbance degree of the vibration of the second direction relative to the vibration of the first direction is greater than an eighteenth threshold value, the rotating speed domain integral value of the vibration of the second direction which is not stable and greater than a nineteenth threshold value, and the rotating speed interval proportion of the vibration of the second direction which is increased is greater than a twentieth threshold value. By applying the first embodiment of the present application, at least the following beneficial effects can be achieved:

1) the method for identifying the non-stationary vibration information number obtains real-time operation data of the wind generating set; determining a standard deviation of vibration data in a first direction and a standard deviation of vibration data in a second direction of the engine room according to the real-time operation data; identifying a signal representing non-stationary vibration of the generator according to the standard deviation of the vibration data in the first direction and the standard deviation of the vibration data in the second direction, wherein the whole process of the identification method in the embodiment of the application is automatically executed by starting in real time; that is to say, the embodiment of the application realizes the real-time identification of the non-stationary vibration signal in the running process of the generator (including the bearing of the generator) of the wind generating set, thereby realizing the early discovery, early processing and early solution of the vibration problem of the generator.

2) In addition, extra detection equipment is not required to be added in the embodiment of the application, the monitoring of the non-steady vibration signal of the running of the generator (including the bearing) of the wind generating set can be realized, the detection cost for detecting the non-steady vibration of the generator is reduced, and the application range of the embodiment of the application is wider.

Example two

Based on the same inventive concept of the first embodiment, the present embodiment explains the technical solution of the first embodiment by specific examples.

Referring to fig. 2, fig. 2 is a schematic flow chart of another method for identifying a vibration signal of a wind turbine generator system according to an embodiment of the present application, and as shown in fig. 2, the method specifically includes:

s201, acquiring real-time operation data of the wind generating set.

Specifically, setting the period of the calculation data as 10 days, namely selecting the real-time operation data of the wind generating set for 10 days; the bearing of the generator does not belong to an easily-consumed product, and the monitoring period does not need high-frequency monitoring.

Optionally, the fields of the real-time operation data include the following corresponding items: the real-time vibration data of the engine room comprises X-direction vibration data and Y-direction vibration data, the X-direction vibration data is X-direction engine room acceleration, and the Y-direction vibration data is Y-direction engine room acceleration.

Optionally, in this embodiment of the application, the first direction is specifically an X direction, and the second direction is specifically a Y direction; the X direction is the axis direction of the engine room, the Y direction is perpendicular to the X direction, and the straight line of the Y direction is coplanar with the straight line of the X direction. The method comprises the steps of utilizing real-time data generated by a wind turbine generator PLC or second-level vibration data recorded by a wind power plant central monitoring System (SCADA), wherein the sampling frequency is generally acquired according to the second level, And the sampling period can be 1-10 seconds. And taking the data in the running state as calculation data. The higher the sampling frequency, the better the data quality.

Optionally, the real-time operation data field is standardized, and the time tag, the X-direction nacelle acceleration, the Y-direction nacelle acceleration, the real-time rotation speed of the generator, and the real-time wind speed are renamed to respectively define: time, x _ ace, y _ ace, rpm, and wind; or renaming field names according to user requirements for program compatibility and uniformity of program parameters of the model.

S202, judging whether the quantity of the real-time vibration data is larger than or equal to a first threshold value or not; and executing S203 when the number of the real-time vibration data is judged to be larger than or equal to the first threshold, and ending the process when the number of the real-time vibration data is judged to be smaller than the first threshold.

Optionally, the first threshold is 10000 stripes.

S203, judging whether the real-time vibration data meet the data quality condition; when the real-time vibration data meets the data quality condition, executing S204; and when the real-time vibration data does not meet the data quality condition, executing S211 and ending the process.

Specifically, when the number of the real-time vibration data with the absolute value of the Y-direction vibration data Y _ ace greater than 1 is determined to be smaller than the second threshold, S204 is executed, and when the number of the real-time vibration data with the absolute value of the Y-direction vibration data greater than 1 is determined to be not smaller than the second threshold, the process is ended. Optionally, the second threshold is 10.

And S204, carrying out primary data screening to obtain real-time vibration data after the primary screening.

Optionally, the real-time vibration data after the first screening satisfies: the absolute value of the X-direction vibration data and the absolute value of the Y-direction vibration data are both smaller than the high boundary value, and the absolute value of the Y-direction vibration data is larger than the low boundary value or the absolute value of the X-direction vibration data is larger than the low boundary value.

Specifically, a low boundary value low _ set is set, a high boundary value high _ set is set, both the low boundary value and the high boundary value are non-negative real numbers, the low boundary value is smaller than the high boundary value, and real-time vibration data of which the absolute value of the vibration data in the X direction and the absolute value of the vibration data in the Y direction of the unit are not between (low _ set and high _ set) (the boundary low _ set is not included) are rejected, so that the remaining real-time vibration data meet the following logic conditions: (abs (x _ ace) < high _ set & abs (y _ ace) < high _ set) & (abs (x _ ace) > low _ set | abs (y _ ace) > low _ set)).

And S205, performing secondary data screening to obtain real-time vibration data subjected to secondary screening.

Optionally, the real-time vibration data after the second screening satisfies: the number of unique values of the X-direction vibration data or the Y-direction vibration data is not less than a third threshold value, the maximum value of the real-time wind speed unique value is at least greater than a cut-in operation wind speed point, the number of real-time wind speeds greater than the cut-in operation wind speed point is at least greater than a fourth threshold value or the duration time is not less than a first time period, and the real-time operation data of real-time wind speeds greater than the threshold wind speed point is not less than a fifth threshold value.

Specifically, the inventor of the present application considers that real-time vibration data has a certain wind speed interval and vibration requirements, and cannot always vibrate or wind speed cannot represent the situation of always no wind. Therefore, in the embodiment of the present application, the number of the unique values x _ ace or y _ ace is not less than 3 (belonging to the third threshold), the maximum value of the unique values of the wind speed is at least greater than 3 (the cut-in operating wind speed point of the unit), the data volume of the wind speed greater than 3 (the cut-in operating wind speed point of the unit) is at least greater than the fourth threshold 80 (or the time is not less than the first time period 10 minutes), and the data volume of the wind speed greater than the threshold wind speed point 5 is not less than the fifth threshold 30. It should be noted that, for example, the value of x _ ace is 0.1, 0.2, and 0.3, if there are multiple identical values of 0.1, only one unique value of 0.1 is calculated here, and there are 3 unique values, which are 0.1, 0.2, and 0.3, respectively.

S206, judging whether the non-empty real-time vibration data quantity meets the aggregation calculation condition; when the non-empty real-time vibration data amount meets the aggregation calculation condition, executing S207; and when the non-empty real-time vibration data amount does not meet the aggregation calculation condition, executing S211 and ending the process.

Specifically, when the amount of the non-empty real-time vibration data is greater than 10, S207 is executed, otherwise, the process ends.

And S207, performing aggregation calculation on the real-time vibration data subjected to the secondary screening to obtain the standard deviation of the vibration data in the X direction and the standard deviation of the vibration data in the Y direction.

Specifically, aggregation calculation is performed according to the rotating speed (accurate to decimal point 1 bit) of the generator of the (wind power generation) unit, the standard deviation sd _ xacc of the X-direction vibration data and the standard deviation sd _ yacc of the Y-direction vibration data at the same rotating speed are obtained, and the vibration stability of the unit at different (generator) rotating speeds is expressed through the standard deviations.

When the pattern recognition is performed on the camber uneven vibration signal at a specific rotation speed of the generator, the following S208-S210 are used for recognition.

And S208, setting a rotating speed interval.

Specifically, a specific rotating speed interval in the aggregation calculation result is selected, that is, a rotating speed area of the normal generator is selected.

For example, for a 1.5MW (megawatt) model, aggregate data with a rotation speed interval of 8.5-16.5 is selected. In the rotation speed section, the value number of the unique rotation speed value range of the power generation rotation speed region in the aggregation calculation result is judged to be not lower than 20 values, and the rotation speed span (the rotation speed interval after the rotation speed is rounded) is not lower than 4 rotation speed intervals. If the condition is not met, the calculation is not continued, and the data condition is not met. If the conditions are met, selecting a data mark of a first rotating speed interval of 10-14 rotating speed intervals as rpm _ xyacce; setting parameters: the rotational speed center position (a position within the selected rotational speed interval is taken as a point in the abnormal vibration region where the occurrence of the abnormality is likely) center _ point is 12.5, the rotational speed center deviation position div _ diff is 0.5, and the vibration standard deviation threshold parameter ypeak _ value is 0.01.

The second rotational speed interval is an interval between the first position value (center _ point-div _ diff) and the second position value (center _ point + div _ diff). The result that the rotating speed is in the second rotating speed interval is marked as rpm _ cenacce; the number of data points in rpm _ cenacce that satisfy sd _ yacc > ═ ypeak _ value is calculated and marked as warp. The coefficient of variation of sd _ yacc in rpm _ cenacce was recalculated: i.e., the ratio of the standard deviation and the mean of sd _ yacc, the result is labeled ycv.

The kurtosis of sd _ yacc in rpm _ xyacce was calculated and the result was labeled peak _ py.

Calculating the number of data points in sd _ yacc > ypeak _ value in rpm _ xyacce, and recording the result as welen _ py; and the maximum, minimum, and polar distance (difference between maximum and minimum) among sd _ yacc in rpm _ xyacce were calculated and recorded as yacc _ max, yacc _ min, yrange, respectively.

Calculating a rotation speed point of maximum vibration instability: the label is yrpm _ pointmax, i.e., the rotation point at which yrpm _ pointmax is the maximum value of sd _ yacc in rpm _ xyacce.

Calculating the maximum allowable variation amplitude in the steady ascending and descending variation in the vibration instability point arch distribution characteristics: labeled yvmax _ diff; that is, the difference calculation is performed on sd _ yacc in rpm _ xyacce, and the maximum value in the absolute value of the calculation result is taken as the maximum variation amplitude of the operation of the neighboring point of the arch distribution characteristic, so that the jump characteristic, such as the jump phenomenon caused by unstable vibration of the rotating speed point, can be avoided.

S209, judging whether the real-time vibration data simultaneously meet a first arch condition and a second arch condition; executing S210 when the real-time vibration data simultaneously meets the first arch condition and the second arch condition; when the real-time vibration data does not satisfy the first and second arching conditions at the same time, S211 is executed, and the process ends.

Specifically, real-time vibration data satisfying a first camber condition and a second camber condition is determined as a camber signal characterizing non-stationary vibration of the generator. And executing S210 when the real-time vibration data meets the first arch condition and the second arch condition, otherwise, ending the process.

Determining that a first vaulting condition war _ xlwtj is satisfied when: the coefficient of variation ycv of the standard deviation of the Y-direction vibration data of the second rotation speed interval is greater than or equal to a sixth threshold, optionally, the sixth threshold is 0.14; and the polar distance value yrange of the standard deviation of the vibration data in the Y direction in the second rotation speed interval is greater than or equal to a seventh threshold value, optionally, the seventh threshold value is 0.004; the number weleny of the standard deviation of the vibration data in the Y direction in the second rotating speed interval is greater than or equal to the eighth threshold value, and optionally, the eighth threshold value is 1; and the maximum value yacc _ max of the standard deviation data of the Y-direction vibration of the second rotation speed section is equal to or greater than a ninth threshold value, optionally, the ninth threshold value is 0.01; the quantity warlen _ py of the standard deviation of the vibration data in the Y direction in the first rotating speed interval is greater than or equal to the tenth threshold value, and optionally, the tenth threshold value is 1; and the absolute value abs (peak _ py) of the kurtosis of the standard deviation of the Y-direction vibration data of the second rotation speed section is greater than an eleventh threshold, optionally, the eleventh threshold is 0.25; and the position yrpm _ pointmax where the standard deviation of the vibration data in the Y direction in the first rotation speed interval is the maximum value is larger than the first position value, and the position yrpm _ pointmax where the standard deviation of the vibration data in the Y direction in the first rotation speed interval is the maximum value is smaller than the second position value.

For example, the first dome condition may be expressed by the expression:

war_xlwtj＝ycv>＝0.14&yrange>＝0.004&warleny>＝1&yacc_max>＝0.01&warlen_py>1&abs(peak_py)>0.25&yrpm_pointmax>(center_point-div_diff)&yrpm_pointmax<(center_point+div_diff)。

the first camber condition is that unstable vibration occurs in a specific rotating speed interval, and the unstable vibration phenomenon is presented in an arch distribution.

Optionally, performing difference calculation on the standard deviation of the vibration data in the Y direction in the first rotation speed interval to obtain a difference calculation value; taking the maximum value in the absolute values of the difference calculated values as the maximum variation amplitude yvmax _ diff of the adjacent point operation of the arch distribution characteristic; when the maximum variation amplitude yvmax _ diff of the adjacent points of the arch distribution characteristics is smaller than a twelfth threshold value, optionally, the twelfth threshold value is 0.003; or the ratio of the maximum variation amplitude yvmax _ diff of the neighboring points of the arch distribution feature to the variation coefficient ycv is less than or equal to a thirteenth threshold value, which is optionally 0.015, and it is determined that the second arch condition is satisfied. Specifically, the second camber condition is war _ xlwtj2 ═ (yvmax _ diff <0.003| yvmax _ diff/ycv ≦ 0.015), which is an allowable local fluctuation amplitude that exhibits the camber characteristic.

S210, triggering and alarming 'abnormal vibration of the generator in the non-steady speed-up process'.

And S211, ending the process.

Referring to fig. 3, fig. 3 is a schematic flow chart of a further method for identifying a vibration signal of a wind turbine generator system according to an embodiment of the present application, and as shown in fig. 3, the method specifically includes:

s301, acquiring real-time operation data of the wind generating set.

The specific method in this step is the same as the specific method in step S201, and is not described herein again.

S302, judging whether the quantity of the real-time vibration data is larger than or equal to a first threshold value or not; if the number of the real-time vibration data is greater than or equal to the first threshold, S303 is executed, and if the number of the real-time vibration data is less than the first threshold, S314 is executed to end the process.

The specific method in this step is the same as the specific method in step S202, and is not described here again.

S303, judging whether the real-time vibration data meet the data quality condition; when the real-time vibration data meets the data quality condition, executing S304; and when the real-time vibration data does not meet the data quality condition, executing the S314 process and ending the process.

The specific method in this step is the same as the specific method in step S203, and is not described herein again.

And S304, performing first-time data screening to obtain real-time vibration data after the first-time screening.

The specific method in this step is the same as the specific method in step S204, and is not described herein again.

And S305, performing secondary data screening to obtain real-time vibration data subjected to secondary screening.

The specific method in this step is the same as the specific method in step S205, and is not described herein again.

S306, judging whether the non-empty real-time vibration data quantity meets the aggregation calculation condition; when the non-empty real-time vibration data amount meets the aggregation calculation condition, executing S307; and when the non-empty real-time vibration data amount does not meet the aggregation calculation condition, executing the S314 and ending the process.

The specific method in this step is the same as the specific method in step S206, and is not described here again.

And S307, performing aggregation calculation on the real-time vibration data subjected to the secondary screening to obtain the standard deviation of the vibration data in the X direction and the standard deviation of the vibration data in the Y direction.

The specific method in this step is the same as the specific method in step S207, and is not described here again.

When the pattern recognition is performed for the linear unsteady state of the generator bearing wear, the recognition is performed using the following S308-S310.

And S308, setting a rotating speed interval.

Specifically, a specific rotating speed area in the polymerization calculation result is selected, namely a rotating speed area of a normal generator is selected, and a first rotating speed interval of 8.5-13 rotating speeds is selected by taking 2.5MW and 3MW machine types as examples; calculating the vibration condition of the vibration rotating speed value range data in the direction of X, Y by the following method: the value number of the unique rotating speed value range in the power generation rotating speed region is not lower than 20 values, and the rotating speed span (rotating speed interval after the rotating speed is integrated) is not lower than 4 rotating speed intervals. If the condition is not met, the calculation is not continued, and the data condition is not met. If the conditions are satisfied, calculating a deviation value sd _ yjx of the stability of the Y-direction vibration to the X-direction vibration; for convenience of statistics, the label with sd _ yjx greater than 0 is marked as red, otherwise, the label is marked as green, and then the label quantity of red and green is counted respectively and marked as tab _ red and tab _ green; calculating the multiple mark of tab _ red as the multiple of tab _ green, namely the multiple of the rotation speed domain range occupied by the disturbance degree of the vibration in the Y direction relative to the vibration in the X direction; if tab _ red >0 and tab _ green >0, then multiple _ rg equals tab _ red/tab _ green; otherwise, judging whether tab _ green is 0 and tab _ red is greater than 0, if so, setting multiple _ rg to 100, otherwise, setting multiple _ rg to 0. The Y-direction vibration fluctuates more obviously irregularly along with the increase of the rotating speed. Calculating the rotating speed domain integral value multiple _ th with unstable and larger vibration in the Y direction as follows: if tab _ red >0 and tab _ green >0, then the sum of all records sd _ yjx of sd _ yjx >0.0002 is calculated and marked as green _ th, the sum of all records sd _ yjx < -0.0002 is calculated and the negative mark is taken as less _ th, multiple _ th is marked as green _ th/less _ th, when tab _ green is 0 and tab _ red >0, the multiple _ th is set to 100, otherwise, 0 is set.

S309, determining relative vibration difference parameters in the X direction and the Y direction.

Specifically, X, Y marks the difference parameter of relative vibration in two directions as ncahr _ yggtth, and when sd _ yacc > sd _ xacc, the relative difference parameter of vibration is 1, otherwise it is 0; and calculating the frequency of 1 appearing in the rotating speed interval (8.5-13.0) set by the vector, marking the frequency as ncahr _ yggth, and expressing the number of points in the rotating speed domain of the relatively unstable vibration rotating speed interval in two directions.

S310, determining a correlation coefficient of the unstable vibration in the Y direction which is linearly increased on the basis of the rotating speed when the vibration in the X direction is stable.

Specifically, when the vibration in the X direction is stable, the correlation coefficients of the unstable vibration in the Y direction which linearly increase on the basis of the rotating speed are marked as cor _ coef, cor _ coefx and cor _ coefy, and the linear proportionality coefficients and the correlation coefficients of the power generation rotating speed point in the first rotating speed interval of 8.5-13 and sd _ yjx are respectively marked as lm _ coef and cor _ coef; correlation coefficients of the rotation speed point and sd _ xacc and the rotation speed point and sd _ yacc are respectively marked as cor _ coefx and cor _ coefy, and linear proportionality coefficients adopt a linear regression fitting method.

And S311, determining the rotating speed interval proportion of the increased vibration in the Y direction.

Specifically, the rotation speed interval proportion of the increase of the Y-direction vibration is marked as rpm _ increase, and when tab _ red >0 and tab _ green > is 0, rpm _ increase is tab _ red/(tab _ red + tab _ green), otherwise rpm _ increase is 0.

S312, judging whether the real-time vibration data meet a first linear condition; when the real-time vibration data satisfies the first linear condition, executing S313; when the real-time vibration data does not satisfy the first linear condition, S314 is executed and the process ends.

Optionally, the first linear condition reatult is determined to be satisfied when: the X-direction and Y-direction relative vibration difference parameter nchar _ ygth is greater than or equal to a fourteenth threshold value, and the first correlation coefficient cor _ coef is greater than a fifteenth threshold value, the correlation coefficient in which the non-stationary vibration in the Y-direction linearly increases as a whole based on the rotation speed when the vibration in the X-direction is stable includes the first correlation coefficient cor _ coef, the second correlation coefficient cor _ coefx, and the third correlation coefficient cor _ coefy, and the second correlation coefficient cor _ coefx is less than or equal to a sixteenth threshold value, and the third correlation coefficient cor _ coefy is greater than a seventeenth threshold value, and the multiple _ rg of the rotation speed domain range occupied by the disturbance of the degree of vibration in the Y-direction with respect to the vibration in the X-direction is greater than an eighteenth threshold value, and the rotation speed domain integrated value multiple _ th in which the vibration in the Y-direction is not stably large is greater than a nineteenth threshold value, and the rotation speed interval proportion rpm _ initial in which the vibration in the Y-direction is increased is greater than a twentieth threshold value.

Specifically, the first linearity condition is satisfied that the measured is 25& cor _ coef >0.6& cor _ coeff > <0.75 & cor _ coefy >0.9& multiple _ rg >1.5& multiple _ th >2& rpm _ increment > -0.6, that is, the fact that the instability of the vibration appearing as the rotation speed increases is expressed by the coefficient parameter of the correlation and the linear fitting coefficient, and the instability of the vibration in the Y direction is more significant and the speed of change is faster and the instability of the vibration in the Y direction exceeds more than half of the instability of the vibration in the X direction when the correlation coefficient in the X direction is set lower and the correlation coefficient in the Y direction is set higher and the rotation speed is relatively stable in the X direction is expressed.

And S313, triggering the abnormal abrasion and vibration of the generator bearing.

And S314, ending the process.

It should be noted that the data parameters related to the present embodiment are expressions of first-order torsion and first-order bending. The method can be used for calculating the operation data of 5-7 seconds in a common sampling period in the current wind power plant once every 7-14 days, and the data volume of 7 days or more is adopted once. In practical application, the data quantity and the data sampling period can be adjusted according to the effect.

The visual graphs of the early warning results of the bearing wear of the 3MW wind turbine output by the identification method of the non-stationary vibration signals provided by the embodiment are shown in fig. 4 and 5. Fig. 4 is a schematic diagram of fluctuation rates (fluctuation is relatively smooth) of two vibration directions of the normal unit changing with the rotation speed according to the embodiment of the application, and as shown in fig. 4, fluctuations of the X-direction cabin acceleration X _ ace and the Y-direction cabin acceleration Y _ ace of the normal unit changing with the rotation speed are relatively smooth. Fig. 5 is a schematic diagram of the fluctuation rate of the abnormal unit in two vibration directions changing with the rotation speed (unit vibration fluctuation rate change of 3MW bearing wear) provided in the embodiment of the present application, and as shown in fig. 5, the unit vibration fluctuation rate change of 3MW bearing wear shows a linear growth characteristic that the instability of the Y-direction cabin acceleration Y _ ace increases with the rotation speed.

For example, the visualized graphs of the non-stable vibration signal monitoring result in the power generation process of the wind generating set output by the non-stable vibration signal identification method provided by the above embodiment are shown in fig. 6, fig. 7 and fig. 8.

Specifically, fig. 6 is a schematic view of a first-order mode of the direct-drive wind turbine generator showing an arch distribution characteristic in a specific rotation speed domain, as shown in fig. 6, the first-order mode of the direct-drive wind turbine generator showing the arch distribution characteristic in the specific rotation speed domain (the mode around 12 turns showing the arch characteristic), and the first-order mode or the first-order bending of the direct-drive wind turbine generator easily causes fatigue damage to a stator and rotor support of the generator, and finally causes a severe influence of mechanical cracking.

Fig. 7 is a schematic view of a wear characteristic of a bearing of a 1.5Mw wind turbine provided in an embodiment of the present application, and as shown in fig. 7, the wear characteristic of the bearing of the 5Mw wind turbine (an arch characteristic appears around 12.5 revolutions), a discrete difference is 0.0028, a discrete variation ratio is 0.007, and the wear of the bearing of the 1.5Mw wind turbine can face abnormal performances such as large bearing clearance and bearing bush indentation of a wind turbine in actual engineering.

FIG. 8 is a schematic representation of the normal characteristic behavior of a 1.5Mw direct drive wind turbine bearing provided by an embodiment of the present application, as shown in FIG. 8, with no camber characteristics occurring between 11-14 turns. In fig. 6, 7 and 8, sd-acid represents the standard deviation sd _ xacc of the vibration data in the X direction or the standard deviation sd _ yacc of the vibration data in the Y direction, rpm (revolutions per minute) represents the number of revolutions per minute of the direct-drive wind turbine generator, X-acid represents the acceleration of the cabin in the X direction, Y-acid represents the acceleration of the cabin in the Y direction, and max-acid represents the instantaneous maximum acceleration.

By applying the third embodiment of the present application, at least the following beneficial effects can be achieved:

in the embodiment of the application, extra detection equipment is not required to be added, the monitoring of the non-stationary vibration signal of the running of the generator (including the bearing) of the wind generating set can be realized, the detection cost for detecting the non-stationary vibration of the generator is reduced, and the application range of the embodiment of the application is wider.

EXAMPLE III

Based on the same inventive concept, the embodiment of the present application further provides an apparatus for identifying a vibration signal of a wind generating set, a schematic structural diagram of the apparatus for identifying a vibration signal of a wind generating set is shown in fig. 9, and an apparatus 1300 for identifying a vibration signal of a wind generating set includes a first processing module 1301, a second processing module 1302, and a third processing module 1303.

The first processing module 1301 is used for acquiring real-time operation data of the wind generating set; the real-time operating data includes: the method comprises the following steps of (1) time tagging, real-time vibration data of a cabin, real-time rotating speed of a generator and real-time wind speed; the real-time vibration data comprises first direction vibration data and second direction vibration data; the first direction is the axis direction of the engine room, and the second direction is vertical to the first direction;

the second processing module 1302 is configured to determine a standard deviation of the first-direction vibration data of the nacelle and a standard deviation of the Y-direction vibration data according to the real-time operation data;

and the third processing module 1303 is configured to identify a signal representing non-stationary vibration of the generator according to the standard deviation of the first-direction vibration data and the standard deviation of the second-direction vibration data.

Optionally, the apparatus 1300 for identifying a vibration signal of a wind turbine generator system further includes a screening module.

Optionally, the screening module is configured to determine whether the number of the real-time vibration data meets a first threshold and whether the real-time vibration data meets a preset data quality condition; when the quantity of the real-time vibration data meets a first threshold value and the real-time vibration data meets a data quality condition, screening the real-time vibration data to obtain screened real-time vibration data; and determining the standard deviation of the vibration data in the first direction and the standard deviation of the vibration data in the second direction based on the rotating speed according to the screened real-time vibration data and the real-time rotating speed.

Optionally, the screening module is specifically configured to determine whether the number of the real-time vibration data of which the absolute value of the vibration data in the second direction is greater than 1 is smaller than a second threshold, where the vibration data includes an acceleration.

Optionally, the screening module is specifically configured to determine, according to the first direction vibration data and the second direction vibration data in the real-time vibration data, real-time vibration data that satisfies the following first data retention condition, as real-time vibration data after the first screening: the absolute value of the first-direction vibration data and the absolute value of the second-direction vibration data are both smaller than the high boundary value, and the absolute value of the second-direction vibration data is larger than the low boundary value or the absolute value of the first-direction vibration data is larger than the low boundary value.

Optionally, the screening module is specifically configured to determine, from the real-time vibration data after the first screening, real-time vibration data that satisfies the following second data retention condition, as the real-time vibration data after the second screening: the number of unique values of the first direction vibration data or the second direction vibration data is not less than a third threshold value, the maximum value of the real-time wind speed unique value is at least greater than a cut-in operation wind speed point, the number of real-time wind speeds greater than the cut-in operation wind speed point is at least greater than a fourth threshold value or the duration is not less than a first time period, and the real-time operation data of real-time wind speeds greater than the threshold wind speed point is not less than a fifth threshold value.

Optionally, the third processing module 1303 is specifically configured to determine a first rotation speed interval according to the standard deviation of the first-direction vibration data and the standard deviation of the second-direction vibration data; determining a second rotating speed interval according to the first rotating speed interval, the rotating speed center position and the rotating speed center deviation position; judging whether a first arch condition and a second arch condition are met or not according to the standard deviation of the second-direction vibration data in the first rotating speed interval, the standard deviation of the second-direction vibration data in the second rotating speed interval, the rotating speed center position, the rotating speed center deviation position and the vibration standard deviation threshold parameter; and when the first and second camber conditions are met, determining real-time vibration data meeting the first and second camber conditions as a camber signal representing non-stationary vibration of the generator.

Optionally, the third processing module 1303 is specifically configured to determine that the first arch condition is satisfied when: the coefficient of variation of the standard deviation of the second-direction vibration data of the second rotation speed interval is greater than or equal to a sixth threshold value, the polar distance value of the standard deviation of the second-direction vibration data of the second rotation speed interval is greater than or equal to a seventh threshold value, the number of the standard deviation of the second-direction vibration data of the second rotation speed interval is greater than or equal to a vibration standard deviation threshold parameter is greater than or equal to an eighth threshold value, the maximum value of the standard deviation data of the second-direction vibration of the second rotation speed interval is greater than or equal to a ninth threshold value, the number of the standard deviation of the second-direction vibration data of the first rotation speed interval is greater than or equal to a vibration standard deviation threshold parameter is greater than a tenth threshold value, the absolute value of the kurtosis of the standard deviation of the second-direction vibration data of the second rotation speed interval is greater than an eleventh threshold value, and the position where the standard deviation of the second-direction vibration data of the first rotation speed interval, and the position where the standard deviation of the second direction vibration data in the first rotating speed interval is the maximum value is smaller than the second position value; carrying out differential calculation on the standard deviation of the vibration data in the second direction in the first rotating speed interval to obtain a differential calculation value; taking the maximum value in the absolute values of the difference calculation values as the maximum variation amplitude of the operation of the adjacent points of the arch distribution characteristics; and when the maximum variation amplitude of the operation of the adjacent points of the arch-shaped distribution characteristics is smaller than a twelfth threshold value, or the ratio of the maximum variation amplitude of the operation of the adjacent points of the arch-shaped distribution characteristics to the variation coefficient is smaller than or equal to a thirteenth threshold value, determining that the second arch condition is satisfied.

Optionally, the third processing module 1303 is specifically configured to determine a first rotation speed interval according to the standard deviation of the first-direction vibration data and the standard deviation of the second-direction vibration data; determining a relative vibration difference parameter between the first direction and the second direction according to the standard deviation of the first direction vibration data and the standard deviation of the second direction vibration data in the first rotation speed interval, a correlation coefficient of overall linear increase of the non-stationary vibration in the Y direction based on the rotation speed when the first direction vibration is stable, a multiple of a rotation speed domain range occupied by the disturbance degree of the second direction vibration relative to the first direction vibration, a rotation speed domain integral value of unstable and larger second direction vibration and a rotation speed interval proportion of larger second direction vibration; judging whether a first linear condition is met or not according to the relative vibration difference parameter in the first direction and the second direction, the correlation coefficient, the multiple of the range of the rotating speed domain and the rotating speed interval proportion; and when the first linear condition is met, determining the real-time vibration data meeting the first linear condition as a linear signal representing the non-stationary vibration of the generator.

The content that does not detail in the recognition device of wind generating set vibration signal that this application embodiment provided can refer to the recognition method of above-mentioned wind generating set vibration signal, and the beneficial effect that the recognition device of wind generating set vibration signal that this application embodiment provided can reach is the same with the recognition method of above-mentioned wind generating set vibration signal, and no longer repeated here.

Based on the same inventive concept, the embodiment of the present application further provides an electronic device, a schematic structural diagram of the electronic device is shown in fig. 10, the electronic device 1400 includes at least one processor 1401, a memory 1402 and a bus 1403, and the at least one processor 1401 is electrically connected to the memory 1402; the memory 1402 is configured to store at least one computer executable instruction, and the processor 1401 is configured to execute the at least one computer executable instruction, so as to execute the steps of any method for identifying a wind turbine generator set vibration signal, as provided in any one of the first to second embodiments or any one alternative embodiment of the present application.

Further, the processor 1403 may be an FPGA (Field-Programmable Gate Array) or other device with logic processing capability, such as an MCU (micro controller Unit) and a CPU (Central processing Unit).

Based on the same inventive concept, the embodiment of the present application provides a computer-readable storage medium, on which a computer program is stored, and the computer program, when executed by a processor, implements the method for identifying the vibration signal of the wind generating set provided by the embodiment of the present application.

The computer readable medium includes, but is not limited to, any type of disk including floppy disks, hard disks, optical disks, CD-ROMs, and magnetic-optical disks, ROMs, RAMs, EPROMs (Erasable Programmable Read-Only memories), EEPROMs, flash memories, magnetic cards, or optical cards. That is, a readable medium includes any medium that stores or transmits information in a form readable by a device (e.g., a computer).

The computer-readable storage medium provided by the present embodiment has the same inventive concept and the same advantages as the foregoing embodiments, and is not repeated herein.

The application of the embodiment of the application has at least the following beneficial effects:

3) On the real-time second-level vibration data conventionally collected in a wind generating set of a wind power plant, mode recognition can be achieved for vibration signals of mechanical damage caused by bending, torsion and the like in a first-order structural mode of a permanent magnet direct-drive wind driven generator, and meanwhile mode recognition can be achieved for similar Bush indentations of a bearing in the running process of the generator.

Those of skill in the art will appreciate that the various operations, methods, steps in the processes, acts, or solutions discussed in this application can be interchanged, modified, combined, or eliminated. Further, other steps, measures, or schemes in various operations, methods, or flows that have been discussed in this application can be alternated, altered, rearranged, broken down, combined, or deleted. Further, steps, measures, schemes in the prior art having various operations, methods, procedures disclosed in the present application may also be alternated, modified, rearranged, decomposed, combined, or deleted.

It should be understood that, although the steps in the flowcharts of the figures are shown in order as indicated by the arrows, the steps are not necessarily performed in order as indicated by the arrows. The steps are not performed in the exact order shown and may be performed in other orders unless explicitly stated herein. Moreover, at least a portion of the steps in the flow chart of the figure may include multiple sub-steps or multiple stages, which are not necessarily performed at the same time, but may be performed at different times, which are not necessarily performed in sequence, but may be performed alternately or alternately with other steps or at least a portion of the sub-steps or stages of other steps.

The foregoing is only a partial embodiment of the present application, and it should be noted that, for those skilled in the art, several modifications and decorations can be made without departing from the principle of the present application, and these modifications and decorations should also be regarded as the protection scope of the present application.

Claims

1. A method for identifying vibration signals of a wind generating set is characterized by comprising the following steps:

acquiring real-time operation data of the wind generating set; the real-time operating data comprises: the method comprises the following steps of (1) time tagging, real-time vibration data of a cabin, real-time rotating speed of a generator and real-time wind speed; the real-time vibration data comprises first direction vibration data and second direction vibration data; the first direction is an axial direction of the nacelle, and the second direction is perpendicular to the first direction;

and identifying a signal representing the non-stationary vibration of the generator according to the standard deviation of the first direction vibration data and the standard deviation of the second direction vibration data.

2. The method of claim 1, further comprising, prior to the determining the standard deviation of the first direction vibration data and the standard deviation of the second direction vibration data of the nacelle:

judging whether the quantity of the real-time vibration data meets a first threshold value and whether the real-time vibration data meets a preset data quality condition;

when the quantity of the real-time vibration data meets a first threshold value and the real-time vibration data meets the data quality condition,

screening the real-time vibration data to obtain screened real-time vibration data;

and determining the standard deviation of the vibration data in the first direction and the standard deviation of the vibration data in the second direction based on the rotating speed according to the screened real-time vibration data and the real-time rotating speed.

3. The method of claim 2, wherein the determining whether the real-time vibration data meets a preset data quality condition comprises:

and judging whether the quantity of the real-time vibration data of which the absolute value is greater than 1 in the second direction is smaller than a second threshold value or not, wherein the vibration data comprise acceleration.

4. The method of claim 3, wherein the screening the real-time vibration data to obtain the screened real-time vibration data comprises:

according to the first direction vibration data and the second direction vibration data in the real-time vibration data, determining the real-time vibration data meeting the following first data retention conditions as the real-time vibration data after the first screening: the absolute value of the first-direction vibration data and the absolute value of the second-direction vibration data are both smaller than the high boundary value, and the absolute value of the second-direction vibration data is larger than the low boundary value or the absolute value of the first-direction vibration data is larger than the low boundary value.

5. The method of claim 4, wherein the screening the real-time vibration data to obtain the screened real-time vibration data further comprises:

determining real-time vibration data meeting the following second data retention conditions from the real-time vibration data after the first screening as the real-time vibration data after the second screening: the number of unique values of the first direction vibration data or the second direction vibration data is not less than a third threshold value, the maximum value of the real-time wind speed unique value is at least greater than a cut-in operation wind speed point, the number of real-time wind speeds greater than the cut-in operation wind speed point is at least greater than a fourth threshold value or the duration is not less than a first time period, and the real-time operation data of the real-time wind speed greater than the threshold wind speed point is not less than a fifth threshold value.

6. The method of claim 1, wherein identifying a signal indicative of non-stationary vibration of the generator based on the standard deviation of the first direction vibration data and the standard deviation of the second direction vibration data comprises:

determining a first rotation speed interval according to the standard deviation of the first direction vibration data and the standard deviation of the second direction vibration data;

determining a second rotating speed interval according to the first rotating speed interval, the rotating speed center position and the rotating speed center deviation position;

judging whether a first arch condition and a second arch condition are met according to the standard deviation of the second-direction vibration data in the first rotating speed interval, the standard deviation of the second-direction vibration data in the second rotating speed interval, the rotating speed center position, the rotating speed center deviation position and a vibration standard deviation threshold parameter;

and when the first and second camber conditions are met, determining the real-time vibration data meeting the first and second camber conditions as a camber signal representing the non-stable vibration of the generator.

7. The method of claim 6, wherein the determining whether the first camber condition and the second camber condition are satisfied based on the standard deviation of the second directional vibration data in the first speed interval, the standard deviation of the second directional vibration data in the second speed interval, the speed center position, the speed center deviation position, and a vibration standard deviation threshold parameter comprises:

determining that the first dome condition is satisfied when: the variation coefficient of the standard deviation of the second direction vibration data of the second rotation speed interval is greater than or equal to a sixth threshold value, the polar distance value of the standard deviation of the second direction vibration data of the second rotation speed interval is greater than or equal to a seventh threshold value, the number of the standard deviation of the second direction vibration data in the second rotation speed interval which is greater than or equal to the vibration standard deviation threshold parameter is greater than or equal to an eighth threshold value, the maximum value of the standard deviation data of the second direction vibration in the second rotation speed interval is greater than or equal to a ninth threshold value, the number of the standard deviation data of the second direction vibration data in the first rotation speed interval which is greater than or equal to the vibration standard deviation threshold parameter is greater than a tenth threshold value, the absolute value of the kurtosis of the standard deviation of the second direction vibration data in the second rotation speed interval is greater than an eleventh threshold value, and the position where the standard deviation of the second direction vibration data in the first rotation speed interval is the maximum value is greater than a first position, and the position where the standard deviation of the second direction vibration data in the first rotating speed interval is the maximum value is smaller than the second position value;

carrying out differential calculation on the standard deviation of the vibration data in the second direction in the first rotating speed interval to obtain a differential calculation value;

taking the maximum value in the absolute values of the difference calculation values as the maximum variation amplitude of the operation of the neighboring points of the arch distribution characteristics;

determining that the second arch condition is satisfied when the maximum variation amplitude of the neighboring points of the arch distribution feature is smaller than a twelfth threshold, or the ratio of the maximum variation amplitude of the neighboring points of the arch distribution feature to the coefficient of variation is smaller than or equal to a thirteenth threshold.

8. The method of claim 1, wherein identifying a signal indicative of non-stationary vibration of the generator based on the standard deviation of the first direction vibration data and the standard deviation of the second direction vibration data comprises:

determining a relative vibration difference parameter between the first direction and the second direction according to the standard deviation of the first direction vibration data and the standard deviation of the second direction vibration data in the first rotation speed interval, a correlation coefficient of overall linear increase of the non-stationary vibration in the second direction based on the rotation speed when the first direction vibration is stable, a multiple of a rotation speed domain range occupied by the disturbance degree of the second direction vibration relative to the first direction vibration, a rotation speed domain integral value of unstable and large second direction vibration and a rotation speed interval proportion of large second direction vibration;

judging whether a first linear condition is met or not according to the relative vibration difference parameter in the first direction and the second direction, the correlation coefficient, the multiple of the rotating speed domain range, the rotating speed domain integral value and the rotating speed interval proportion;

and when the first linear condition is met, determining the real-time vibration data meeting the first linear condition as a linear signal representing the non-stationary vibration of the generator.

9. An identification device for vibration signals of a wind generating set, characterized in that the device comprises:

the first processing module is used for acquiring real-time operation data of the wind generating set; the real-time operating data comprises: the method comprises the following steps of (1) time tagging, real-time vibration data of a cabin, real-time rotating speed of a generator and real-time wind speed; the real-time vibration data comprises first direction vibration data and second direction vibration data; the first direction is an axial direction of the nacelle, and the second direction is perpendicular to the first direction;

and the third processing module is used for identifying a signal representing the non-stationary vibration of the generator according to the standard deviation of the first-direction vibration data and the standard deviation of the second-direction vibration data.

10. An electronic device, comprising: a processor, a memory, and a bus;

the bus is used for connecting the processor and the memory;

the memory is used for storing operation instructions;

the processor is used for executing the method for identifying the vibration signal of the wind generating set according to any one of the claims 1 to 8 by calling the operation instruction.

11. A computer-readable storage medium, characterized in that a computer program is stored which, when being executed by a processor, is adapted to carry out the method for identifying a wind park vibration signal according to any one of claims 1 to 8.