Fan tower health management method based on modal model and storage medium

Abstract

The invention relates to the field of wind power equipment structure health monitoring, and discloses a modal model-based fan tower health management method. A sensor array is deployed to collect vibration mode, environment and operation data, three working conditions of shutdown / low wind speed / rated wind speed are divided in an initial operation period, and a reference mode frequency and a vibration mode matrix are generated by adopting wind speed weighted average to establish a health baseline. And in the operation stage, single working condition data is screened through a dynamic time window, a real-time modal frequency deviation ratio beta and a vibration mode matrix MAC value are calculated, and a fault risk level is output in combination with a multi-level early warning threshold. The problems that a traditional method is high in false alarm rate and inaccurate in damage positioning are solved, and the safety monitoring precision and efficiency of the fan tower drum structure are remarkably improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of structural health monitoring of electrical equipment, and in particular to a fan tower health state dynamic evaluation and early warning method based on vibration modal analysis. BACKGROUND

[0002] As a core representative of clean and renewable energy, wind power has developed rapidly in recent years. With the increase in the number of wind turbines and the accumulation of operating time, the structural safety detection of wind turbine towers has become increasingly prominent. Traditional detection methods (such as manual visual inspection or regular shutdown detection) have significant drawbacks: detection personnel need to climb high-altitude towers, which is high-risk and low-efficiency; minor defects (such as cracks and corrosion) are easily overlooked, leading to potential failure escalation; the detection process relies on subjective experience and lacks real-time and quantitative accuracy. In addition, wind turbine operation is affected by dynamic environmental factors (such as wind speed and temperature), and traditional methods cannot establish a working condition adaptive model, making it difficult to achieve early warning. SUMMARY

[0003] In view of the technical defects and technical disadvantages in the prior art, the embodiments of the present application provide a wind turbine health management method based on a modal model to overcome the above problems or at least partially solve the above problems, and the specific solutions are as follows: As a first aspect of the present application, a wind turbine health management method based on a modal model is provided, comprising the following steps: Step S1: deploying a sensor array at a predetermined height position of the wind turbine tower, and collecting vibration modal data, environmental parameters and operating state data in real time; wherein: the vibration modal data includes a sway frequency and a mode shape matrix; the environmental parameters include wind speed, wind direction and temperature; the operating state data includes power generation power; Step S2: based on the data collected in step S1, during the initial operating period after the installation or overhaul of the wind turbine, the data is divided into three working conditions of shutdown, low wind speed and rated wind speed according to the mapping relationship between wind speed and power generation power; the data of each working condition is processed by weighted average to generate a reference modal frequency and a reference mode shape matrix, and a health baseline database of working condition-modal mapping is established; Step S3: in the operation stage of the wind turbine, based on the health baseline database of step S2, real-time data is collected through a dynamic time window, and a data set meeting the continuity of a single working condition is selected; the deviation ratio of the real-time modal frequency relative to the reference value is calculated, as well as the correlation coefficient of the real-time mode shape matrix and the reference mode shape matrix, i.e. the MAC value; Step S4: according to the calculation results of the deviation ratio and the MAC value in step S3, match the preset multi-level early warning threshold table, and output the alarm signal corresponding to the fault risk level.

[0004] Further, step S2 includes: Step S2.1: Based on the power generation and wind speed data collected in step S1, perform dynamic operating condition classification: When the power generation is zero and the wind speed is less than level 3, it is marked as a shutdown state; When the power generation is greater than zero but less than the rated power, it is marked as a low wind speed state; When the power generation capacity equals the rated power, it is marked as the rated wind speed state; Specifically, when dividing the operating conditions, the jump points where the power generation power changes from 0 to non-0 and their adjacent data are removed; Step S2.2: For each type of operating condition data, assign a weight related to wind speed, where the closer the wind speed is to the median of the typical wind speed range for the corresponding operating condition, the larger the weight. The weights are used to calculate the weighted average of the sway frequencies to generate the reference modal frequencies for the corresponding working conditions; the weights are also used to calculate the weighted average of the column vectors of the mode shape matrix to generate the reference mode shape matrix for the corresponding working conditions. Step S2.3: Store the reference modal frequencies and reference mode shape matrices of the three types of working conditions as a queryable mapping database according to the working condition labels, which will serve as the comparison benchmark for subsequent health status assessment.

[0005] Furthermore, S2.2 includes: S2.21, Determine the median reference wind speed for each type of operating condition based on historical data or design parameters. v mid ; S2.22, For the sway frequency data of each type of operating condition, the reference modal frequency is generated according to the following rules: Calculate the weights of each data point Where v is the real-time wind speed and σ is an adjustable standard deviation parameter; the reference modal frequency for this type of operating condition is output by weighting the swaying frequency with the aforementioned weights. S2.23, mode shape matrix for each type of working condition Φ n×m Generate the reference mode shape matrix according to the following rules: The column vector of the mode shape for each mode (j=1,2, The reference mode shape matrix is ​​output by individually weighting m, using the same weights w as in S2.22, and averaging by column vector weights. Φ n×m ; S2.24, store the reference modal frequencies and reference mode shape matrices into the database according to the working condition labels.

[0006] Furthermore, the correlation coefficient (MAC) between the real-time mode shape matrix and the reference mode shape matrix is ​​calculated, including: Extract the reference mode matrix for the current operating condition from the healthy baseline database. Φ n×mAnd obtain the real-time mode shape matrix of the real-time dynamic window. Φ' n×m ; calculate Φ' n×m and Φ n×m Projection matrix: P = Φ' n×m × Φ n×m T ,in P Let T be an m×m matrix, and T denote the transpose calculation; calculate Φ' n×m Autocorrelation diagonal vector E B = diag ( Φ' n×m × Φ' n×m T )and Φ n×m autocorrelation diagonal vector E A = diag ( Φ n×m T × Φ n×m ); pass E A and E B The outer product generates the denominator matrix Thus, matrix D is an m×n matrix; Squaring the projection matrix P element by element: P squared = P. × P, Where .× indicates element-wise division or multiplication; Dividing the numerator matrix by the denominator matrix yields M = P squared . / D ,in . / indicates element-wise division; The resulting M matrix is ​​an m×m matrix. Taking the values ​​of its diagonal elements gives the MAC value, where MAC = {M...} 11 M 22 ,...,M mm}, M ii M is the correlation calculation result between the baseline i-th mode and the real-time i-th mode. ii The closer to 1, the higher the correlation; equal to 1, they are completely identical.

[0007] Further, according to claim 1, the method is characterized in that step S3 includes: Based on the operating characteristics and frequency of changes in operating conditions of the wind turbine, the time window parameters are adaptively set to ensure that the window covers at least one complete operating cycle. Modal parameter calculation is triggered if and only if both the wind speed range and the power generation range within the window are within the preset operating condition range threshold. Real-time detection of events switching between zero and non-zero power generation values; removal of event points and their adjacent data to form a clean dataset. Synchronous execution based on the clean dataset: Calculate the relative ratio of the weighted average real-time modal frequency to the reference frequency; calculate the MAC value between the real-time mode shape matrix and the reference mode shape matrix; The output relative proportions and MAC values ​​are stored in a time-series database, and statistical feature vectors are generated in response to periodic analysis commands. (The preset operating condition interval thresholds are dynamically generated through historical operating condition clustering, and the periodic analysis commands are triggered by an external maintenance system or a fixed time period.)

[0008] Furthermore, the method also includes: performing temperature compensation on the real-time modal frequencies, specifically including: A mapping model between temperature and modal frequency is established based on historical data, and the temperature response coefficient is extracted. When the real-time temperature deviates from the reference temperature used when the health baseline was established, the compensation process is activated. The modal frequencies acquired in real time are corrected according to the following rules:

[0009] f raw The original frequency measurement value is represented by α, the temperature response coefficient is represented by α, and the normalized difference between the real-time temperature and the reference temperature is represented by ΔT.

[0010] Furthermore, the hierarchical early warning decision in step S4 is executed through the following conditional trigger chain: Receive the modal frequency relative ratio β and MAC value output in step S3; Level 1 warning trigger: When β is between the baseline preset lower limit and upper limit of fluctuation and the MAC value is not lower than the morphological correlation threshold, a level 1 warning signal is output, indicating that there is no abnormality in the structure; Level 2 warning trigger: When β is lower than the baseline preset fluctuation lower limit but does not reach the moderate degradation threshold, and the MAC value is lower than the morphological correlation threshold, a level 2 warning signal is output to indicate the risk of connector degradation. Level 3 warning trigger: When β reaches the moderate degradation threshold but not the severe degradation threshold, and the MAC value is lower than the morphological stability threshold, a level 3 warning signal is output to indicate the initiation of structural damage. Level 4 warning trigger: When β exceeds the severe deterioration threshold and the dynamic amplitude exceeds the design margin, a Level 4 warning signal is output to indicate the risk of structural instability.

[0011] Furthermore, the method further includes: achieving compatibility with multi-height tower barrels through a dimension-independent modal analysis framework, specifically including: S8.1: According to the distribution density of the tower barrel flange and the target height range, dynamically configure the installation positions and quantities of the sensor arrays to make the sensor topology cover the key modal deformation areas of the tower barrel; S8.2: Package the real-time collected vibration modal data into a dynamic dimension mode shape matrix Φ' n×m , which matrix serves as the real-time mode shape matrix in step S3. Here, n is the number of effective sensors of the current tower barrel, determined by the deployment rules of S8.1; m is the standard modal order, bound to the physical properties of the tower barrel and independent of height; S8.3: In the health baseline establishment stage, generate a reference mode shape matrix according to the same rules as S8.2 to ensure the dimension alignment between the dynamic dimension mode shape matrix and the reference mode shape matrix; S8.4: When executing step S3, apply a dimension-independent MAC value calculation core to the dynamic dimension mode shape matrix and the reference mode shape matrix, so that the output index only depends on modal correlation rather than physical size.

[0012] Furthermore, the method further includes: enhancing the fault location ability through a local-global collaborative diagnosis mechanism, specifically including: S9.1: Based on the tower barrel structure stress distribution model, identify the high-damage risk sections as the local diagnosis target areas; S9.2: In the data acquisition of step S1, synchronously activate the sensor subset corresponding to the local diagnosis target area to generate a target area mode shape matrix Ψ k×m , where: k is the number of effective sensors in the target area, k < n; m is the same modal order as the global matrix; S9.3: When performing the health assessment in step S3, synchronously calculate: The global MAC value is based on the complete modal analysis of the full-tower mode shape matrix Φ n×m ; The local diagnosis target area MAC value is based on the same modal calculation process analysis of the target area mode shape matrix Ψ k×m ; S9.4: When the deterioration rate of the local diagnosis target area MAC value relative to the full-tower MAC value exceeds the preset sensitivity threshold, generate a target area damage location signal; S9.5: Perform a logical OR operation on the target area damage location signal and the warning level in step S4 to trigger a targeted maintenance instruction.

[0013] As a second aspect of the present invention, a computer-readable storage medium is provided, wherein a computer program is stored therein, and when executed by a computer, the computer program causes the computer to perform the wind turbine tower health management method based on modal model as described above.

[0014] The present invention has the following beneficial effects: 1. Improves detection accuracy and real-time performance, enabling the identification of minute structural changes (such as loosening or corrosion) and reducing the false negative rate; 2. Reduce manual intervention and downtime costs, and improve wind turbine operation safety through automated early warning systems; 3. The system is highly robust and suitable for towers of different heights and environments, supporting long-term structural health management; 4. Enhance data comparability and reduce false alarms by using temperature compensation and wind speed weighting. Attached Figure Description

[0015] Figure 1 This is a flowchart illustrating a wind turbine tower health management method based on modal modeling, provided as an embodiment of the present invention. Detailed Implementation

[0016] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only a part of the present invention, and not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0017] To enable those skilled in the art to better understand the technical solutions of the present invention, exemplary embodiments of the present invention are described below in conjunction with the accompanying drawings, including various details of the embodiments of the present invention to aid understanding. These should be considered merely exemplary. Therefore, those skilled in the art should recognize that various changes and modifications can be made to the embodiments described herein without departing from the scope and spirit of the present invention. Similarly, for clarity and brevity, descriptions of well-known functions and structures are omitted in the following description.

[0018] Where there is no conflict, the various embodiments of the present invention and the features thereof may be combined with each other.

[0019] As used herein, the term “and / or” includes any and all combinations of one or more related enumerated entries.

[0020] The terminology used herein is for the purpose of describing particular embodiments only and is not intended to limit the invention. As used herein, the singular forms “a” and “the” are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will also be understood that when the terms “comprising” and / or “made of” are used in this specification, the presence of the stated feature, integral, step, operation, element, and / or component is specified, but the presence or addition of one or more other features, integrals, steps, operations, elements, components, and / or groups thereof is not excluded. Terms such as “connected” or “linked” are not limited to physical or mechanical connections but can include electrical connections, whether direct or indirect.

[0021] Unless otherwise specified, all terms used herein (including technical and scientific terms) have the same meaning as commonly understood by one of ordinary skill in the art. It will also be understood that terms such as those defined in commonly used dictionaries should be interpreted as having the meaning consistent with their meaning in the context of the relevant art and the invention, and will not be interpreted as having an idealized or overly formal meaning unless expressly so defined herein.

[0022] In the technical solution of this invention, the collection, storage, use, processing, transmission, provision, and disclosure of user personal information all comply with relevant laws and regulations and do not violate public order and good morals. The use of user data in this technical solution follows relevant national laws and regulations (e.g., the "Information Security Technology - Personal Information Security Specification"). For example: appropriate measures are taken for personal information access control; restrictions are imposed on the display of personal information; the purpose of using personal information does not exceed the scope of direct or reasonable association; and explicit identity targeting is eliminated when using personal information to avoid precisely locating a specific individual.

[0023] To address at least one of the technical problems existing in the aforementioned related technologies, the present invention provides a wind turbine tower health management method based on modal model. Figure 1 A flowchart illustrating a wind turbine tower health management method based on modal modeling, provided in this embodiment of the invention, includes the following steps: Step S1: Deploy a sensor array at a preset height position on the wind turbine tower to collect vibration modal data, environmental parameters, and operating status data in real time; wherein: the vibration modal data includes sway frequency and mode shape matrix; the environmental parameters include wind speed, wind direction, and temperature; and the operating status data includes power generation. Step S2: Based on the data collected in Step S1, during the initial operating cycle after the wind turbine is newly installed or overhauled, the data is divided into three operating conditions: shutdown, low wind speed, and rated wind speed, according to the mapping relationship between wind speed and power generation. The data of each operating condition is weighted and averaged to generate the reference modal frequency and reference mode shape matrix, and a healthy baseline database of operating condition-modal mapping is established. Step S3: During the wind turbine operation phase, based on the health baseline database of step S2, real-time data is collected through a dynamic time window, and datasets that meet the continuity of a single operating condition are selected; the deviation ratio of the real-time modal frequency relative to the reference value and the correlation coefficient between the real-time mode matrix and the reference mode matrix, i.e., the MAC value, are calculated. Step S4: Based on the deviation ratio and MAC value calculation results from Step S3, match the preset multi-level early warning threshold table and output the alarm signal corresponding to the fault risk level.

[0024] This invention integrates vibration mode (frequency + mode shape), environmental parameters (wind speed / temperature), and operating status (power) data to overcome the limitations of traditional methods that rely on a single vibration signal. It establishes a benchmark database according to operating conditions (shutdown / low wind speed / rated wind speed) to eliminate misjudgments caused by environmental interference and improve the accuracy of health assessment. Finally, it uses the modal frequency deviation ratio (quantifying stiffness degradation) and MAC value (mode shape correlation) for collaborative judgment, which can identify minor damage earlier than traditional methods that only monitor frequency.

[0025] In some embodiments, step S2 includes: Step S2.1: Based on the power generation and wind speed data collected in step S1, perform dynamic operating condition classification: When the power generation is zero and the wind speed is less than level 3, it is marked as a shutdown state; When the power generation is greater than zero but less than the rated power, it is marked as a low wind speed state; When the power generation capacity equals the rated power, it is marked as the rated wind speed state; Specifically, when dividing the operating conditions, the jump points where the power generation power changes from 0 to non-0 and their adjacent data are removed; Step S2.2: For each type of operating condition data, assign a weight related to wind speed, where the closer the wind speed is to the median of the typical wind speed range for the corresponding operating condition, the larger the weight. The weights are used to calculate the weighted average of the sway frequencies to generate the reference modal frequencies for the corresponding working conditions; the weights are also used to calculate the weighted average of the column vectors of the mode shape matrix to generate the reference mode shape matrix for the corresponding working conditions. Step S2.3: Store the reference modal frequencies and reference mode shape matrices of the three types of working conditions as a queryable mapping database according to the working condition labels, which will serve as the comparison benchmark for subsequent health status assessment.

[0026] In the above embodiments, power jump point filtering (such as power from 0 to non-0) is used to avoid polluting the benchmark library with transitional operating condition data and improve data purity; data close to the median of typical wind speed are given higher weights so that the benchmark modal parameters are closer to the actual operating conditions (such as tower flexibility is more significant at low wind speeds).

[0027] In some embodiments, S2.2 includes: S2.21, Determine the median reference wind speed for each type of operating condition based on historical data or design parameters. v mid ; S2.22, For the sway frequency data of each type of operating condition, the reference modal frequency is generated according to the following rules: Calculate the weights of each data point Where v is the real-time wind speed and σ is an adjustable standard deviation parameter; the reference modal frequency for this type of operating condition is output by weighting the swaying frequency with the aforementioned weights. S2.23, mode shape matrix for each type of working condition Φ n×m Generate the reference mode shape matrix according to the following rules: The column vector of the mode shape for each mode (j=1,2, The reference mode shape matrix is ​​output by individually weighting m, using the same weights w as in S2.22, and averaging by column vector weights. Φ n×m ; S2.24, store the reference modal frequencies and reference mode shape matrices into the database according to the working condition labels.

[0028] In the above embodiments, an exponential weighting function is used. Strengthen the contribution of typical wind speed data, suppress the noise impact of extreme wind speeds (such as gusts), improve the stability of the reference frequency, independently weight the column vector of each mode, preserve the physical characteristics of the mode shape (such as the shape of the first bending mode), and avoid the shape distortion caused by overall averaging.

[0029] In some embodiments, calculating the correlation coefficient (MAC) between the real-time mode shape matrix and the reference mode shape matrix includes: Extract the reference mode matrix for the current operating condition from the healthy baseline database. Φ n×m And obtain the real-time mode shape matrix of the real-time dynamic window. Φ' n×m ; calculate Φ' n×m and Φ n×m Projection matrix: P = Φ' n×m ×Φ n×m T ,in P Let T be an m×m matrix, and T denote the transpose calculation; calculate Φ' n×m Autocorrelation diagonal vector E B = diag ( Φ' n×m × Φ' n×m T )and Φ n×m autocorrelation diagonal vector E A = diag ( Φ n×m T × Φ n×m ); pass E A and E B The outer product generates the denominator matrix Thus, matrix D is an m×n matrix; Squaring the projection matrix P element by element: P squared = P. × P, Where .× indicates element-wise division or multiplication; Dividing the numerator matrix by the denominator matrix yields M = P squared . / D ,in . / indicates element-wise division; The resulting M matrix is ​​an m×m matrix. Taking the values ​​of its diagonal elements gives the MAC value, where MAC = {M...} 11 M 22 ,...,M mm}, M ii M is the correlation calculation result between the baseline i-th mode and the real-time i-th mode. ii The closer to 1, the higher the correlation; equal to 1, they are completely identical.

[0030] In the above embodiments, the projection matrix P and the autocorrelation vector are used... E A / E BThe element-wise operation decouples the MAC value calculation from the number of sensors n, supports dynamic sensor deployment (such as when a sensor is damaged), and only requires that the modal order m be consistent (such as the first 3 modes inherent in the tower), thus solving the problem of data incomparability caused by the difference in the number of sensors for towers of different heights.

[0031] In some embodiments, step S3 includes: Based on the operating characteristics and frequency of changes in operating conditions of the wind turbine, the time window parameters are adaptively set to ensure that the window covers at least one complete operating cycle. Modal parameter calculation is triggered if and only if both the wind speed range and the power generation range within the window are within the preset operating condition range threshold. Real-time detection of events switching between zero and non-zero power generation values; removal of event points and their adjacent data to form a clean dataset. Synchronous execution based on the clean dataset: Calculate the relative ratio of the weighted average real-time modal frequency to the reference frequency, and calculate the MAC value of the real-time mode shape matrix and the reference mode shape matrix; The output relative proportions and MAC values ​​are stored in a time series database, and statistical feature vectors are generated in response to periodic analysis commands.

[0032] The preset operating condition range threshold is dynamically generated through historical operating condition clustering, and the periodic analysis command is triggered by an external maintenance system or a fixed time period.

[0033] In the above embodiments, by using an adaptive window to cover the complete operating cycle (such as a 10-minute low wind speed steady state), misjudgments caused by mixed operating data are avoided, power switching event points (such as shutdown → startup) are removed in real time, transient vibration interference is eliminated, and the dataset is purified to reduce MAC calculation errors.

[0034] In some embodiments, the method further includes: performing temperature compensation on the real-time modal frequencies, specifically including: A mapping model between temperature and modal frequency is established based on historical data, and the temperature response coefficient is extracted. When the real-time temperature deviates from the reference temperature used when the health baseline was established, the compensation process is activated. The modal frequencies acquired in real time are corrected according to the following rules:

[0035] f raw The original frequency measurement value is represented by α, the temperature response coefficient is represented by α, and the normalized difference between the real-time temperature and the reference temperature is represented by ΔT. The corrected f adj Enter step S3 and verify whether the frequency deviation converges after compensation.

[0036] In the above embodiments, the temperature response coefficient is used. α (e.g., for steel, α≈-0.03% / ℃) Correct modal frequencies to eliminate measurement bias caused by temperature differences and reduce seasonal false alarm rates. Trigger only when the temperature difference ΔT exceeds the limit (e.g., >5℃) to avoid unnecessary computational overhead.

[0037] In some embodiments, the hierarchical early warning decision in step S4 is executed through the following conditional trigger chain: Receive the modal frequency relative ratio β and MAC value output in step S3; Level 1 warning trigger: When β is between the baseline preset lower limit and upper limit of fluctuation and the MAC value is not lower than the morphological correlation threshold, a level 1 warning signal is output, indicating that there is no abnormality in the structure; Level 2 warning trigger: When β is lower than the baseline preset fluctuation lower limit but does not reach the moderate degradation threshold, and the MAC value is lower than the morphological correlation threshold, a level 2 warning signal is output to indicate the risk of connector degradation. Level 3 warning trigger: When β reaches the moderate degradation threshold but not the severe degradation threshold, and the MAC value is lower than the morphological stability threshold, a level 3 warning signal is output to indicate the initiation of structural damage. Level 4 warning trigger: When β exceeds the severe degradation threshold and the dynamic amplitude exceeds the design margin, a Level 4 warning signal is output to indicate the risk of structural instability.

[0038] The table below shows a specific form of early warning at various levels provided in the embodiments of the present invention.

[0039]

[0040] In the above embodiments, the second-level warning indicates the risk of connector deterioration, such as loose bolts (improving sensitivity); the third-level warning indicates the initiation of structural damage, such as crack initiation (improving early detection rate); and the fourth-level warning indicates the risk of structural instability. By combining frequency deviation (stiffness deterioration) and MAC value (morphological distortion), the damage to connectors and the damage to the structural body are distinguished, guiding targeted maintenance.

[0041] In some embodiments, the method further includes: achieving multi-height tower compatibility through a dimension-independent modal analysis architecture, specifically including: S8.1: Based on the distribution density of tower flanges and the target height range, dynamically configure the installation position and number of sensor arrays to ensure that the sensor topology covers the key modal deformation zone of the tower. S8.2: Encapsulate the real-time acquired vibration modal data into a dynamic dimensional mode shape matrix. Φ n×mThis matrix serves as the real-time mode shape matrix for step S3, where n is the number of effective sensors in the current tower, determined by the deployment rules in S8.1; and m is the standard mode order, which is bound to the physical properties of the tower and is independent of its height. S8.3: During the healthy baseline establishment phase, generate the reference mode matrix according to the same rules as in S8.2 to ensure that the dynamic dimension mode matrix is ​​aligned with the dimension of the reference mode matrix; S8.4: When performing step S3, adjust the dynamic dimension mode matrix. Φ n×m The application of dimension-independent MAC value calculation core with the reference mode matrix makes the output index depend only on modal correlation rather than physical dimensions.

[0042] The number of sensors (n) varies depending on the tower height, resulting in different mode shape matrices. Φ n×m The row dimension changes, which traditional matrix operations cannot directly compare. The core of this invention uses mathematical reconstruction to make the MAC value calculation depend only on the shape similarity of the mode shape, and is independent of the number of sensors n (i.e., the tower height). MAC is calculated as follows: ; For each mode i, the scalar MAC value is calculated independently, freeing it from the constraint of matrix dimension n. Here, MACi represents the criterion value of the i-th mode, i.e., the calculated output, reflecting the similarity between the real-time mode shape and the reference mode shape. Let be the reference mode shape vector, representing the standardized vibration mode of the i-th mode in the healthy baseline database, i.e. Φ n×m The i-th column, The real-time mode shape vector represents the standardized vibration mode of the i-th mode within the current dynamic window, i.e. Φ' n×m The i-th column.

[0043] In the above embodiments, the number of rows n in the mode shape matrix adapts to the number of sensors (e.g., n=8 for a 70m tower, n=12 for a 140m tower), while m is fixed (e.g., analyzing the first 3 modes), enabling towers of different heights to share the same analysis kernel. Core Formula By eliminating the dependence on n, the MAC results for 70m / 140m towers can be directly compared.

[0044] In some embodiments, the method further includes: enhancing fault location capabilities through a local-global collaborative diagnostic mechanism, specifically including: S9.1: Based on the stress distribution model of the tower structure, identify high-damage-risk sections as local diagnostic target areas; S9.2: In the data acquisition step S1, the sensor subset corresponding to the target area is activated synchronously to generate the target area vibration mode matrix.Ψ k×m , where: k is the number of effective sensors in the target area, k < n; m is the modal order same as the global matrix; S9.3: When performing the health assessment in step S3, synchronously calculate: The global MAC value is based on the complete modal analysis of the full-tower mode shape matrix Φ n×m ; The local MAC value is based on the analysis of the same modal calculation process of the target area mode shape matrix Ψ k×m ; S9.4: When the deterioration rate of the target area MAC value relative to the full-tower MAC value exceeds the preset sensitivity threshold, generate a target area damage location signal; S9.5: Perform a logical OR operation on the target area damage location signal and the warning level in step S4 to trigger a targeted maintenance instruction.

[0045] In the above embodiments, through the deterioration rate of the target area MAC (such as the Ψ k×m MAC decline rate in the flange area is faster than the global), the fault location accuracy is improved from "tower barrel level" to "flange level", and the local sensor subset focuses on the high-stress area (such as the bottom of the tower barrel), making the change amplitude of the MAC value of the micro-crack magnified by more than several times.

[0046] An embodiment of the present invention also provides a computer-readable medium. A computer program is stored on the computer-readable medium, where when the program is executed by a processor, it implements the steps in any one of the above embodiments of the modal model-based fan tower barrel health management method. The computer-readable storage medium can be a volatile or non-volatile computer-readable storage medium.

[0047] An embodiment of the present invention also provides a computer program product, including computer-readable code, or a non-volatile computer-readable storage medium carrying the computer-readable code. When the computer-readable code runs in the processor of an electronic device, the processor in the electronic device executes the above modal model-based fan tower barrel health management method.

[0048] Those skilled in the art will understand that all or some of the steps, systems, and apparatuses disclosed above, and their functional modules / units, can be implemented as software, firmware, hardware, or suitable combinations thereof. In hardware implementations, the division between functional modules / units mentioned above does not necessarily correspond to the division of physical components; for example, a physical component may have multiple functions, or a function or step may be performed collaboratively by several physical components. Some or all physical components may be implemented as software executed by a processor, such as a central processing unit, digital signal processor, or microprocessor, or as hardware, or as an integrated circuit, such as an application-specific integrated circuit (ASIC). Such software can be distributed on a computer-readable storage medium, which may include computer storage media (or non-transitory media) and communication media (or transient media).

[0049] As is known to those skilled in the art, the term computer storage medium includes volatile and non-volatile, removable and non-removable media implemented in any method or technology for storing information, such as computer-readable program instructions, data structures, program modules, or other data. Computer storage media includes, but is not limited to, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM), static random access memory (SRAM), flash memory or other memory technologies, portable compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical disc storage, magnetic cartridges, magnetic tape, disk storage or other magnetic storage devices, or any other medium that can be used to store desired information and is accessible to a computer. Furthermore, it is known to those skilled in the art that communication media typically contain computer-readable program instructions, data structures, program modules, or other data in modulated data signals such as carrier waves or other transmission mechanisms, and may include any information delivery medium.

[0050] The computer-readable program instructions described herein can be downloaded from computer-readable storage media to various computing / processing devices, or downloaded via a network, such as the Internet, local area network, wide area network, and / or wireless network, to an external computer or external storage device. The network may include copper transmission cables, fiber optic transmission, wireless transmission, routers, firewalls, switches, gateway computers, and / or edge servers. A network adapter card or network interface in each computing / processing device receives the computer-readable program instructions from the network and forwards them to the computer-readable storage media in the respective computing / processing device.

[0051] The computer program instructions used to perform the operations of this invention may be assembly instructions, instruction set architecture (ISA) instructions, machine instructions, machine-dependent instructions, microcode, firmware instructions, state setting data, or source code or object code written in any combination of one or more programming languages, including object-oriented programming languages ​​such as Smalltalk, C++, etc., and conventional procedural programming languages ​​such as the "C" language or similar programming languages. The computer-readable program instructions may be executed entirely on the user's computer, partially on the user's computer, as a standalone software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In cases involving a remote computer, the remote computer may be connected to the user's computer via any type of network—including a local area network (LAN) or a wide area network (WAN)—or may be connected to an external computer (e.g., via the Internet using an Internet service provider). In some embodiments, electronic circuitry, such as programmable logic circuitry, field-programmable gate arrays (FPGAs), or programmable logic arrays (PLAs), is personalized by utilizing state information from the computer-readable program instructions. This electronic circuitry can execute the computer-readable program instructions to implement various aspects of the invention.

[0052] The computer program product described herein can be implemented specifically through hardware, software, or a combination thereof. In one alternative embodiment, the computer program product is specifically embodied in a computer storage medium; in another alternative embodiment, the computer program product is specifically embodied in a software product, such as a software development kit (SDK), etc.

[0053] Various aspects of the present invention are described herein with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It should be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer-readable program instructions.

[0054] These computer-readable program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing apparatus to produce a machine such that, when executed by the processor of the computer or other programmable data processing apparatus, they create means for implementing the functions / actions specified in one or more blocks of the flowchart and / or block diagram. These computer-readable program instructions can also be stored in a computer-readable storage medium that causes a computer, programmable data processing apparatus, and / or other device to operate in a particular manner; thus, the computer-readable medium storing the instructions comprises an article of manufacture that includes instructions for implementing aspects of the functions / actions specified in one or more blocks of the flowchart and / or block diagram.

[0055] Computer-readable program instructions may also be loaded onto a computer, other programmable data processing apparatus, or other device to cause a series of operational steps to be performed on the computer, other programmable data processing apparatus, or other device to produce a computer-implemented process, thereby causing the instructions executed on the computer, other programmable data processing apparatus, or other device to perform the functions / actions specified in one or more boxes of a flowchart and / or block diagram.

[0056] The flowcharts and block diagrams in the accompanying drawings illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer program products according to various embodiments of the present invention. In this regard, each block in a flowchart or block diagram may represent a module, segment, or portion of an instruction, which contains one or more executable instructions for implementing a specified logical function. In some alternative implementations, the functions marked in the blocks may occur in a different order than those shown in the drawings. For example, two consecutive blocks may actually be executed substantially in parallel, and they may sometimes be executed in reverse order, depending on the functions involved. It should also be noted that each block in the block diagrams and / or flowcharts, and combinations of blocks in the block diagrams and / or flowcharts, may be implemented using a dedicated hardware-based system that performs the specified function or action, or using a combination of dedicated hardware and computer instructions.

[0057] Example embodiments have been disclosed herein, and while specific terminology has been used, it is for illustrative purposes only and should be construed as such, and is not intended to be limiting. In some instances, it will be apparent to those skilled in the art that features, characteristics, and / or elements described in conjunction with particular embodiments may be used alone, or in combination with features, characteristics, and / or elements described in conjunction with other embodiments, unless otherwise expressly indicated. Therefore, those skilled in the art will understand that various changes in form and detail may be made without departing from the scope of the invention as set forth in the appended claims.

[0058] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the protection scope of the present invention.

Claims

1. A wind turbine health management method based on modal models, characterized in that, Includes the following steps: Step S1: Deploy a sensor array at a preset height position on the wind turbine tower to collect vibration modal data, environmental parameters, and operating status data in real time; wherein: the vibration modal data includes sway frequency and mode shape matrix; the environmental parameters include wind speed, wind direction, and temperature; and the operating status data includes power generation. Step S2: Based on the data collected in Step S1, during the initial operating cycle after the wind turbine is newly installed or overhauled, the data is divided into three operating conditions: shutdown, low wind speed, and rated wind speed, according to the mapping relationship between wind speed and power generation. The data of each operating condition is weighted and averaged to generate the reference modal frequency and reference mode shape matrix, and a healthy baseline database of operating condition-modal mapping is established. Step S3: During the wind turbine operation phase, based on the health baseline database of step S2, real-time data is collected through a dynamic time window, and datasets that meet the continuity of a single operating condition are selected; the deviation ratio of the real-time modal frequency relative to the reference value and the correlation coefficient between the real-time mode matrix and the reference mode matrix, i.e., the MAC value, are calculated. Step S4: Based on the deviation ratio and MAC value calculation results from Step S3, match the preset multi-level early warning threshold table and output the alarm signal corresponding to the fault risk level.

2. The wind turbine tower health management method based on modal model according to claim 1, characterized in that, Step S2 includes: Step S2.1: Based on the power generation and wind speed data collected in step S1, perform dynamic operating condition classification: When the power generation is zero and the wind speed is less than level 3, it is marked as a shutdown state; When the power generation is greater than zero but less than the rated power, it is marked as a low wind speed state; When the power generation capacity equals the rated power, it is marked as the rated wind speed state; Specifically, when dividing the operating conditions, the jump points where the power generation power changes from 0 to non-0 and their adjacent data are removed; Step S2.2: For each type of operating condition data, assign a weight related to wind speed, where the closer the wind speed is to the median of the typical wind speed range for the corresponding operating condition, the larger the weight. The weights are used to calculate the weighted average of the sway frequencies to generate the reference modal frequencies for the corresponding working conditions; the weights are also used to calculate the weighted average of the column vectors of the mode shape matrix to generate the reference mode shape matrix for the corresponding working conditions. Step S2.3: Store the reference modal frequencies and reference mode shape matrices of the three types of working conditions as a queryable mapping database according to the working condition labels, which will serve as the comparison benchmark for subsequent health status assessment.

3. The wind turbine tower health management method based on modal model according to claim 2, characterized in that, S2.2 includes: S2.21, Based on historical data or design parameters, determine the median reference wind speed v for each type of operating condition. mid ; S2.22, For the sway frequency data of each type of operating condition, the reference modal frequency is generated according to the following rules: Calculate the weights of each data point Where v is the real-time wind speed and σ is an adjustable standard deviation parameter; the reference modal frequency for this type of operating condition is output by weighting the swaying frequency with the aforementioned weights. S2.23, the mode shape matrix Φ for each type of working condition n×m Generate the reference mode shape matrix according to the following rules: The column vector of the mode shape for each mode (j=1,2, The reference mode shape matrix Φ is output by individually weighting m, using the same weights w as in S2.22, and averaging by column vector weights. n×m ; S2.24, store the reference modal frequencies and reference mode shape matrices into the database according to the working condition labels.

4. The wind turbine tower health management method based on modal model according to claim 2, characterized in that, The correlation coefficient (MAC) between the real-time mode shape matrix and the reference mode shape matrix is ​​calculated as follows: Extract the reference mode matrix Φ for the current operating condition from the health baseline database. n×m And obtain the real-time mode shape matrix Φ′ of the real-time dynamic window. n×m ; Calculate Φ′ n×m With Φ n×m Projection matrix: P=Φ′ n×m ×Φ n×m T , where P is an m×m projection matrix and T represents transpose calculation; Calculate Φ′ n×m The autocorrelation diagonal vector E B =diag(Φ′ n×m ×Φ′ n×m T ) and Φ n×m The autocorrelation diagonal vector E A =diag(Φ n×m T ×Φ n×m ); via E A With E B The outer product generates the denominator matrix Thus, matrix D is an m×n matrix; Squaring the projection matrix P element by element: P squared =P.×P, where .× indicates element-wise division and multiplication; Dividing the numerator matrix by the denominator matrix, we get M=P squared . / D, where . / represents element-wise division; The resulting M matrix is ​​an m×m matrix. Taking the values ​​of its diagonal elements gives the MAC value, where MAC = {M...} 11 M 22 ,...,M mm }, M ii M is the correlation calculation result between the baseline i-th mode and the real-time i-th mode. ii The closer to 1, the higher the correlation; equal to 1, they are completely identical.

5. The wind turbine tower health management method based on modal model according to claim 1, characterized in that, According to the method of claim 1, step S3 includes: Based on the operating characteristics and frequency of changes in operating conditions of the wind turbine, the time window parameters are adaptively set to ensure that the window covers at least one complete operating cycle. Modal parameter calculation is triggered if and only if both the wind speed range and the power generation range within the window are within the preset operating condition range threshold. Real-time detection of events switching between zero and non-zero power generation values; removal of event points and their adjacent data to form a clean dataset. Synchronous execution based on the clean dataset: Calculate the relative ratio of the weighted average real-time modal frequency to the reference frequency; calculate the MAC value between the real-time mode shape matrix and the reference mode shape matrix; The output relative proportions and MAC values ​​are stored in a time-series database, and statistical feature vectors are generated in response to periodic analysis commands. (The preset operating condition interval thresholds are dynamically generated through historical operating condition clustering, and the periodic analysis commands are triggered by an external maintenance system or a fixed time period.) 6. The wind turbine tower health management method based on modal model according to claim 1, characterized in that, The method further includes: temperature compensation for real-time modal frequencies, specifically including: A mapping model between temperature and modal frequency is established based on historical data, and the temperature response coefficient is extracted. When the real-time temperature deviates from the reference temperature used when the health baseline was established, the compensation process is activated. The modal frequencies acquired in real time are corrected according to the following rules: ； f raw The original frequency measurement value is represented by α, the temperature response coefficient is represented by α, and the normalized difference between the real-time temperature and the reference temperature is represented by ΔT.

7. The wind turbine tower health management method based on modal model according to claim 1, characterized in that, The hierarchical early warning decision in step S4 is executed through the following condition trigger chain: Receive the modal frequency relative ratio β and MAC value output in step S3; Level 1 warning trigger: When β is between the baseline preset lower limit and upper limit of fluctuation and the MAC value is not lower than the morphological correlation threshold, a level 1 warning signal is output, indicating that there is no abnormality in the structure; Level 2 warning trigger: When β is lower than the baseline preset fluctuation lower limit but does not reach the moderate degradation threshold, and the MAC value is lower than the morphological correlation threshold, a level 2 warning signal is output to indicate the risk of connector degradation. Level 3 warning trigger: When β reaches the moderate degradation threshold but not the severe degradation threshold, and the MAC value is lower than the morphological stability threshold, a level 3 warning signal is output to indicate the initiation of structural damage. Level 4 warning trigger: When β exceeds the severe degradation threshold and the dynamic amplitude exceeds the design margin, a Level 4 warning signal is output to indicate the risk of structural instability.

8. The wind turbine tower health management method based on modal model according to claim 1, characterized in that, The method further includes: achieving multi-height tower compatibility through a dimension-independent modal analysis architecture, specifically including: S8.1: Based on the distribution density of tower flanges and the target height range, dynamically configure the installation position and number of sensor arrays to ensure that the sensor topology covers the key modal deformation zone of the tower. S8.2: Encapsulate the real-time acquired vibration modal data into a dynamic dimensional mode shape matrix Φ′ n×m This matrix serves as the real-time mode shape matrix for step S3, where n is the number of effective sensors in the current tower and m is the standard mode order. S8.3: During the healthy baseline establishment phase, generate the reference mode matrix according to the same rules as in S8.2 to ensure that the dynamic dimension mode matrix is ​​aligned with the dimension of the reference mode matrix; S8.4: When performing step S3, apply a dimension-independent MAC value calculation core to the dynamic dimension mode matrix and the reference mode matrix, so that the output index depends only on modal correlation rather than physical size.

9. The wind turbine tower health management method based on modal model according to claim 1, characterized in that, The method further includes: enhancing fault location capabilities through a local-global collaborative diagnostic mechanism, specifically including: S9.1: Based on the stress distribution model of the tower structure, identify high-damage-risk sections as local diagnostic target areas; S9.2: In the data acquisition of step S1, synchronously activate the sensor subset corresponding to the target area to generate the target area mode shape matrix Ψ k×m , where: k is the number of effective sensors in the target area, k < n; m is the modal order the same as the global matrix; S9.3: During the health assessment in step S3, calculate simultaneously: The global MAC value is based on the full tower mode matrix Φ n×m Complete modal analysis; Local MAC values ​​are based on the target area mode shape matrix Ψ k×m Analysis of the same modal calculation process; S9.4: When the degradation rate of the target area MAC value relative to the total tower MAC value exceeds the preset sensitivity threshold, a target area damage localization signal is generated. S9.5: Perform a logical OR operation between the target area damage location signal and the warning level in step S4 to trigger a targeted maintenance command.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program, which, when executed by a computer, causes the computer to perform the wind turbine tower health management method based on a modal model as described in any one of claims 1 to 9.

Images