CN116757101B - Cabin wind speed correction method and system based on mechanism model and neural network - Google Patents
Cabin wind speed correction method and system based on mechanism model and neural network Download PDFInfo
- Publication number
- CN116757101B CN116757101B CN202311047504.4A CN202311047504A CN116757101B CN 116757101 B CN116757101 B CN 116757101B CN 202311047504 A CN202311047504 A CN 202311047504A CN 116757101 B CN116757101 B CN 116757101B
- Authority
- CN
- China
- Prior art keywords
- wind speed
- data
- frequency
- low
- scada
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D17/00—Monitoring or testing of wind motors, e.g. diagnostics
- F03D17/005—Monitoring or testing of wind motors, e.g. diagnostics using computation methods, e.g. neural networks
- F03D17/006—Estimation methods
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
- F03D—WIND MOTORS
- F03D17/00—Monitoring or testing of wind motors, e.g. diagnostics
- F03D17/009—Monitoring or testing of wind motors, e.g. diagnostics characterised by the purpose
- F03D17/025—Monitoring or testing of wind motors, e.g. diagnostics characterised by the purpose for calibrating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/044—Recurrent networks, e.g. Hopfield networks
- G06N3/0442—Recurrent networks, e.g. Hopfield networks characterised by memory or gating, e.g. long short-term memory [LSTM] or gated recurrent units [GRU]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/06—Wind turbines or wind farms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Life Sciences & Earth Sciences (AREA)
- Software Systems (AREA)
- Geometry (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- General Health & Medical Sciences (AREA)
- Biomedical Technology (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Health & Medical Sciences (AREA)
- Computer Hardware Design (AREA)
- Mechanical Engineering (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Sustainable Energy (AREA)
- Chemical & Material Sciences (AREA)
- Medical Informatics (AREA)
- Combustion & Propulsion (AREA)
- Sustainable Development (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Wind Motors (AREA)
Abstract
Description
技术领域Technical field
本申请涉及数据处理和数据传输领域,更具体的,涉及一种基于机理模型和神经网络的机舱风速修正方法和系统。This application relates to the fields of data processing and data transmission, and more specifically, to a cabin wind speed correction method and system based on a mechanism model and a neural network.
背景技术Background technique
随着风电技术的发展,在风场进行经济效益评估、发电量评估、机群布置、运行控制时,需要获取准确的风速数据。由于真实风速不易直接测量,风机评估计算时通常将数据采集与监视控制系统(Supervisory ControI And Data AcquiSition System ,简称SCADA)记录的机舱风速作为分析的重要参数。但是,由于受到叶片转动和机舱引起的气流畸变等因素影响,SCADA测得的风速值往往不能直接反映叶轮的真实风速,必须经过修正后才能使用,这种对风速的修正称为机舱传递函数。With the development of wind power technology, it is necessary to obtain accurate wind speed data when conducting economic benefit assessment, power generation assessment, fleet layout, and operation control in wind farms. Since the real wind speed is not easy to measure directly, the nacelle wind speed recorded by the Supervisory Control And Data AcquiSition System (SCADA) is usually used as an important parameter for analysis when calculating wind turbines. However, due to factors such as blade rotation and airflow distortion caused by the nacelle, the wind speed value measured by SCADA often cannot directly reflect the true wind speed of the impeller and must be corrected before use. This correction to the wind speed is called the nacelle transfer function.
目前,对SCADA风速进行修正的方法主要分为理论计算修正法与函数拟合法。其中,理论计算修正法主要基于空气动力学理论,利用风力发电机组的多项参数和SCADA记录的运行数据通过理论计算来实现机舱风速修正,但是由于其参数数量较多且内部机理较为复杂,此类方法计算精度不足且无法覆盖不同机型;函数拟合法通过拟合测风塔风速与机舱风速直接的函数关系对机舱风速进行修正,该方法简单易用,但是严重依赖于测风塔数据,而测风塔的配置场地要求苛刻、费用高昂,通常只选取其中一台典型风电机组为其配备测风塔进行特性评估,无法满足全场所有风机需求,因此,这类方法应用成本高、应用场景有限。At present, the methods for correcting SCADA wind speed are mainly divided into theoretical calculation correction methods and function fitting methods. Among them, the theoretical calculation correction method is mainly based on aerodynamic theory, using multiple parameters of the wind turbine and the operating data recorded by SCADA to achieve cabin wind speed correction through theoretical calculations. However, due to its large number of parameters and complex internal mechanism, this method The calculation accuracy of the class method is insufficient and cannot cover different aircraft types; the function fitting method corrects the cabin wind speed by fitting the direct functional relationship between the wind speed of the wind tower and the wind speed of the cabin. This method is simple and easy to use, but it relies heavily on the wind tower data. However, the configuration site of the wind measurement tower is demanding and expensive. Usually, only one typical wind turbine is selected to equip it with a wind measurement tower for characteristic evaluation, which cannot meet the needs of all wind turbines in the entire site. Therefore, the application cost of this method is high and the application cost is high. Scenes are limited.
发明内容Contents of the invention
为了克服了现有技术的不足,本发明提供了一种运用神经网络、理论特性、模型评估与修正等技术并且考虑风电运行工况的机舱风速修正方法和系统,具体提出了一种基于机理模型和神经网络的机舱风速修正方法和系统。本发明通过在叶轮处安装激光雷达,测量真实风速,之后把物理机理与数据模型结合,具有计算速度快,精度高的特点且泛化性好,可广泛用于不同机型的风力发电机风速修正。In order to overcome the shortcomings of the existing technology, the present invention provides a nacelle wind speed correction method and system that uses neural network, theoretical characteristics, model evaluation and correction and other technologies and considers wind power operating conditions. Specifically, a method and system based on a mechanism model are proposed and neural network cabin wind speed correction method and system. This invention measures the real wind speed by installing a laser radar at the impeller, and then combines the physical mechanism with the data model. It has the characteristics of fast calculation speed, high accuracy and good generalization, and can be widely used to measure the wind speed of different types of wind turbines. Correction.
第一方面,本发明提供了一种基于机理模型和神经网络的机舱风速修正方法,包括如下步骤:In the first aspect, the present invention provides a cabin wind speed correction method based on a mechanism model and a neural network, which includes the following steps:
选取SCADA系统中测得的风速数据和功率数据,将所述风速数据和功率数据根据预设规则进行划分,得到变桨工况数据集和未变桨工况数据集;Select the wind speed data and power data measured in the SCADA system, divide the wind speed data and power data according to preset rules, and obtain a pitch-changing operating condition data set and a non-pitch operating condition data set;
其中,所述根据预设规则进行划分包括:判断所述风速数据是否小于额定风速数据,若是,则判定所述风速数据为未变桨工况时风速;若否,则判定所述风速数据为变桨工况时风速;Wherein, the dividing according to the preset rules includes: determining whether the wind speed data is less than the rated wind speed data; if so, determining that the wind speed data is the wind speed under non-pitch operating conditions; if not, determining that the wind speed data is Wind speed during pitch operation;
所述SCADA系统中的风速数据和功率数据分辨率为1s级;The resolution of wind speed data and power data in the SCADA system is 1s level;
将所述变桨工况数据集和未变桨工况数据集平均转化成分辨率为30s级的数据作为特征数据集;Convert the pitch-changing operating condition data set and the non-pitch operating condition data set into average data with a resolution of 30s as a feature data set;
获取空气密度信息、扫风面积信息,将所述空气密度信息、所述扫风面积信息、所述风速数据和所述功率数据根据预设的经验公式进行计算,得到对应工况下的理论风速;Obtain the air density information and sweep area information, calculate the air density information, sweep area information, wind speed data and power data according to the preset empirical formula to obtain the theoretical wind speed under the corresponding working conditions. ;
其中,所述预设的经验公式为:Among them, the preset empirical formula is:
; ;
其中是功率数据,/>是风速数据,/>是空气密度信息,/>是扫风面积信息;in is the power data,/> is the wind speed data,/> is the air density information,/> It is the sweep area information;
获取真实风速,之后通过小波变换将所述理论风速/>和所述真实风速/>分别进行分解得到对应的高频数据和低频数据,求得相应的高频残差/>和低频残差/>,将所述高频残差/>和所述低频残差/>作为目标数据集,所述风速数据为特征数据集;Get true wind speed , and then use wavelet transform to convert the theoretical wind speed/> and the actual wind speed/> Decompose them separately to obtain the corresponding high-frequency data and low-frequency data, and obtain the corresponding high-frequency residuals/> and low-frequency residual/> , convert the high-frequency residual/> and the low frequency residual/> As the target data set, the wind speed data is a feature data set;
其中,通过激光雷达测得的真实风速,再通过小波变换将真实风速/>和理论风速/>分别进行分解,得到对应的高频数据/>、/>和低频数据/>、/>,然后将/>和/>相减得到真实风速高频数据与理论风速高频数据的高频残差/>,将/>和/>相减得到真实风速低频数据与理论风速低频数据的低频残差/>,具体的,得到真实风速数据后通过小波变换将真实风速数据分解成高频数据和低频数据,并记录对应时间点,然后根据真实风速的时间点将理论风速分解成高频数据和低频数据;Among them, the true wind speed measured by lidar , and then use wavelet transform to convert the real wind speed/> and theoretical wind speed/> Decompose them separately to obtain the corresponding high-frequency data/> ,/> and low frequency data/> ,/> , then // and/> Subtract the high-frequency residual of the real wind speed high-frequency data and the theoretical wind speed high-frequency data/> , will/> and/> Subtract to obtain the low-frequency residual of the real wind speed low-frequency data and the theoretical wind speed low-frequency data/> , Specifically, after obtaining the real wind speed data, the real wind speed data is decomposed into high-frequency data and low-frequency data through wavelet transformation, and the corresponding time points are recorded, and then the theoretical wind speed is decomposed into high-frequency data and low-frequency data according to the time point of the real wind speed;
将所述特征数据集作为输入,所述目标数据集作为输出,使用神经网络进行训练,用优化方法确定时间步长,根据工况和高低频数据,建立风速-残差预测模型;Use the characteristic data set as input and the target data set as output, use a neural network for training, use an optimization method to determine the time step, and establish a wind speed-residual prediction model based on operating conditions and high and low frequency data;
将一定时间步长的SCADA风速数据输入至训练好的风速-残差预测模型中,得到对应时间步长的高频残差预测值数据和低频残差预测值数据;Input SCADA wind speed data of a certain time step into the trained wind speed-residual prediction model, and obtain high-frequency residual prediction value data and low-frequency residual prediction value data corresponding to the time step;
将所述SCADA风速数据与高频残差预测值数据和低频残差预测值数据进行线性相加,得到最终的真实风速修正值。The SCADA wind speed data is linearly added to the high-frequency residual prediction value data and the low-frequency residual prediction value data to obtain the final true wind speed correction value.
可选的,其中,通过激光雷达测得风机轮毂前方的真实风速;Optional, where the real wind speed in front of the wind turbine hub is measured through lidar ;
通过SCADA系统测得机舱风速。Measuring cabin wind speed through SCADA system .
可选的,其中,所述将所述特征数据集作为输入,所述目标数据集作为输出,使用神经网络进行训练,用优化方法确定时间步长,根据工况和高低频数据,建立风速-残差预测模型包括:Optionally, the feature data set is used as input, the target data set is used as output, a neural network is used for training, an optimization method is used to determine the time step, and the wind speed- Residual prediction models include:
使用神经网络进行训练,建立两种工况下的风速-残差预测模型,共计四个:一是变桨工况下SCADA风速作为输入,高频残差作为输出的高频残差预测模型;二是变桨工况下SCADA风速作为输入,低频残差作为输出的低频残差预测模型;三是未变桨工况下SCADA风速作为输入,高频残差作为输出的高频残差预测模型;四是未变桨工况下SCADA风速作为输入,低频残差作为输出的低频残差预测模型。Use neural networks for training to establish wind speed-residual prediction models under two working conditions, a total of four: one is a high-frequency residual prediction model in which SCADA wind speed is used as input and high-frequency residuals are output under pitch working conditions; The second is a low-frequency residual prediction model in which SCADA wind speed is used as input under pitch changing conditions and low-frequency residuals are used as output. The third is a high-frequency residual prediction model in which SCADA wind speed is used as input under non-pitching conditions and high-frequency residuals are used as output. ; The fourth is a low-frequency residual prediction model using SCADA wind speed as input and low-frequency residual as output under non-pitch operating conditions.
可选的,其中,在所述建立风速-残差预测模型过程中,使用伽马检测技术对训练时间步长和预测时间步长进行优化求解,采用逐层随机跳动方法对各层参数进行优化组合求解。Optionally, in the process of establishing the wind speed-residual prediction model, gamma detection technology is used to optimize the training time step and prediction time step, and the layer-by-layer random jumping method is used to optimize the parameters of each layer. Combination solution.
可选的,其中,所述逐层随机跳动方法包括如下步骤:Optionally, the layer-by-layer random jumping method includes the following steps:
1)确定随机跳动的参数、跳动范围和跳动步长:1) Determine the parameters, jump range and jump step size of random jumps:
确定需要寻优的参数为第二层LSTM单元数、dropout值、激活函数;第四层全连接层单元数、激活函数;设定各参数的范围和跳动步长;激活函数为sigmod函数、tanh函数和relu函数随机选择,参数值计算公式如下:Determine the parameters that need to be optimized as the second layer LSTM unit number, dropout value, and activation function; the fourth layer fully connected layer unit number and activation function; set the range and jump step size of each parameter; the activation function is the sigmod function, tanh The function and relu function are randomly selected, and the parameter value calculation formula is as follows:
; ;
为当前轮计算值,/>为上一轮计算值,/>为范围最小值,/>为范围最大值,为跳动步长,/>为在给定范围内随机选取整数的函数,/>; Calculates the value for the current round,/> Calculated value for the previous round,/> is the minimum value of the range,/> is the maximum value of the range, is the jump step size,/> is a function that randomly selects integers within a given range,/> ;
2)第一轮时,在给定范围内随机初始化所选层的单元数、dropout值、映射维度和激活函数,然后将一定量的风速和残差数据集分成训练集和测试集进行训练测试并计算精度,得到对应参数的精度值,记录该参数组合和结果;然后将所选层的单元数和dropout值加上一个随机正整数或负整数乘上跳动步长得到下一组的参数值,其中,所述参数值不能超出给定范围,激活函数也随机选择,重复执行五轮;2) In the first round, the number of units, dropout value, mapping dimension and activation function of the selected layer are randomly initialized within a given range, and then a certain amount of wind speed and residual data sets are divided into training sets and test sets for training and testing. And calculate the accuracy, get the accuracy value of the corresponding parameter, record the parameter combination and result; then add a random positive or negative integer to the number of units and dropout value of the selected layer multiplied by the jump step size to get the next set of parameter values. , where the parameter values cannot exceed the given range, and the activation function is also randomly selected and executed repeatedly for five rounds;
3)计算五组随机参数的精度,选择精度最高的参数组合,然后将第四层的参数按步骤2)执行,同样计算五组参数的精度后选择精度最高的参数组合,最后得到最优的参数组合。3) Calculate the accuracy of five groups of random parameters, select the parameter combination with the highest accuracy, and then execute the parameters of the fourth layer according to step 2). Also calculate the accuracy of the five groups of parameters and select the parameter combination with the highest accuracy, and finally get the optimal Parameter combination.
可选的,其中,获取SCADA风速数据并且根据该SCADA风速数据判断工况后,输入至对应的高频残差预测模型和低频残差预测模型,得到对应时间步长的高频残差和低频残差/>,最后将机舱风速/>加上所述对应时间步长的高频残差/>和低频残差/>得到修正的真实风速/>,即Optionally, after obtaining the SCADA wind speed data and judging the working conditions based on the SCADA wind speed data, input it into the corresponding high-frequency residual prediction model and low-frequency residual prediction model to obtain the high-frequency residual of the corresponding time step. and low-frequency residual/> , and finally the cabin wind speed/> Add the high-frequency residuals of the corresponding time steps/> and low-frequency residual/> Get corrected true wind speed/> ,Right now
。 .
第二方面,本发明提供了一种基于机理模型和神经网络的机舱风速修正系统,用于实现第一方面所述的一种基于机理模型和神经网络的机舱风速修正方法,包括:In a second aspect, the present invention provides a cabin wind speed correction system based on a mechanism model and a neural network, which is used to implement a cabin wind speed correction method based on a mechanism model and a neural network described in the first aspect, including:
划分模块,用于选取SCADA系统中测得的风速数据和功率数据,将所述风速数据和功率数据根据预设规则进行划分,得到变桨工况数据集和未变桨工况数据集;A division module, used to select wind speed data and power data measured in the SCADA system, divide the wind speed data and power data according to preset rules, and obtain a pitch-changing operating condition data set and a non-pitch operating condition data set;
计算模块,用于获取空气密度信息、扫风面积信息,将所述空气密度信息、所述扫风面积信息、所述风速数据和所述功率数据根据预设的经验公式进行计算,得到对应工况下的理论风速;The calculation module is used to obtain air density information and sweep area information, and calculate the air density information, sweep area information, wind speed data and power data according to a preset empirical formula to obtain the corresponding work The theoretical wind speed under ;
分解模块,用于获取真实风速,之后通过小波变换将所述理论风速/>和所述真实风速/>分别进行分解得到对应的高频数据和低频数据,求得相应的高频残差/>和低频残差/>,将所述高频残差/>和所述低频残差/>作为目标数据集,所述风速数据为特征数据集;Decomposition module for obtaining true wind speed , and then use wavelet transform to convert the theoretical wind speed/> and the actual wind speed/> Decompose them separately to obtain the corresponding high-frequency data and low-frequency data, and obtain the corresponding high-frequency residuals/> and low-frequency residual/> , convert the high-frequency residual/> and the low frequency residual/> As the target data set, the wind speed data is a feature data set;
建立模块,用于将所述特征数据集作为输入,所述目标数据集作为输出,使用神经网络进行训练,用优化方法确定时间步长,根据工况和高低频数据,建立风速-残差预测模型;Establish a module for taking the characteristic data set as input and the target data set as output, using a neural network for training, using an optimization method to determine the time step, and establishing wind speed-residual predictions based on working conditions and high and low frequency data. Model;
输入模块,用于将一定时间步长的SCADA风速数据输入至训练好的风速-残差预测模型中,得到对应时间步长的高频残差预测值数据和低频残差预测值数据;The input module is used to input SCADA wind speed data of a certain time step into the trained wind speed-residual prediction model to obtain high-frequency residual prediction value data and low-frequency residual prediction value data corresponding to the time step;
输出模块,用于将所述SCADA风速数据与高频残差预测值数据和低频残差预测值数据进行线性相加,得到最终的真实风速修正值。The output module is used to linearly add the SCADA wind speed data to the high-frequency residual prediction value data and the low-frequency residual prediction value data to obtain the final true wind speed correction value.
本发明公开的一种基于机理模型和神经网络的机舱风速修正方法和系统,选取SCADA系统中测得的风速数据和功率数据,将SCADA风速按额定风速划分为不同工况风速数据集,通过经验公式计算得到理论风速;利用小波变换计算真实风速与理论风速的高低频残差,采用神经网络建立SCADA风速与高低频残差的关系,将来自实际运行风电场的SCADA风速数据输入已训练的神经网络,获取对应高低频残差;将SCADA风速数据与残差数据一一对应,线性相加得到真实风速修正值。本发明具有计算速度快,精度高的特点且泛化性好,可广泛用于不同机型的风力发电机风速修正。The invention discloses a cabin wind speed correction method and system based on a mechanism model and a neural network. The wind speed data and power data measured in the SCADA system are selected, and the SCADA wind speed is divided into different working condition wind speed data sets according to the rated wind speed. Through experience The theoretical wind speed is calculated using the formula; wavelet transform is used to calculate the high- and low-frequency residuals between the real wind speed and the theoretical wind speed, and a neural network is used to establish the relationship between SCADA wind speed and high- and low-frequency residuals. The SCADA wind speed data from the actual operating wind farm is input into the trained neural network. network to obtain the corresponding high and low frequency residuals; correspond the SCADA wind speed data and the residual data one-to-one, and linearly add them to obtain the true wind speed correction value. The invention has the characteristics of fast calculation speed, high accuracy and good generalization, and can be widely used for wind speed correction of different types of wind turbines.
附图说明Description of the drawings
图1示出了本发明一种基于机理模型和神经网络的机舱风速修正方法步骤的流程示意图;Figure 1 shows a schematic flow chart of the steps of a cabin wind speed correction method based on a mechanism model and a neural network according to the present invention;
图2示出了本发明一种基于机理模型和神经网络的机舱风速修正方法的架构图;Figure 2 shows the architecture diagram of a cabin wind speed correction method based on a mechanism model and neural network according to the present invention;
图3示出了本发明一种神经网络结构图;Figure 3 shows a neural network structure diagram of the present invention;
图4示出了本发明一种基于机理模型和神经网络的机舱风速修正系统的模块示意图。Figure 4 shows a schematic module diagram of a cabin wind speed correction system based on a mechanism model and a neural network according to the present invention.
具体实施方式Detailed ways
为了能够更清楚地理解本发明的上述目的、特征和优点,下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是,在不冲突的情况下,本申请的实施例及实施例中的特征可以相互组合。In order to more clearly understand the above objects, features and advantages of the present invention, the present invention will be described in further detail below in conjunction with the accompanying drawings and specific embodiments. It should be noted that, as long as there is no conflict, the embodiments of the present application and the features in the embodiments can be combined with each other.
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是,本发明还可以采用其他不同于在此描述的其他方式来实施,因此,本发明的保护范围并不受下面公开的具体实施例的限制。Many specific details are set forth in the following description in order to fully understand the present invention. However, the present invention can also be implemented in other ways different from those described here. Therefore, the protection scope of the present invention is not limited by the specific details disclosed below. Limitations of Examples.
如图1和图2示出了本发明一种基于机理模型和神经网络的机舱风速修正方法步骤的流程示意图和架构图,一种基于机理模型和神经网络的机舱风速修正方法100,包括:Figures 1 and 2 show a schematic flowchart and architecture diagram of the steps of a cabin wind speed correction method based on a mechanism model and a neural network according to the present invention. A cabin wind speed correction method 100 based on a mechanism model and a neural network includes:
S102:选取SCADA系统中测得的风速数据和功率数据,将所述风速数据和功率数据根据预设规则进行划分,得到变桨工况数据集和未变桨工况数据集;S102: Select the wind speed data and power data measured in the SCADA system, divide the wind speed data and power data according to the preset rules, and obtain a pitch-changing operating condition data set and a non-pitch operating condition data set;
该步骤用于建立原始数据集,选取SCADA系统中的风速和功率数据,根据风速是否达到额定风速将风速和功率数据分为变桨工况风速数据集和未变桨工况风速数据集,具体的:This step is used to establish the original data set, select the wind speed and power data in the SCADA system, and divide the wind speed and power data into the wind speed data set under changing pitch conditions and the wind speed data set under unchanging pitch conditions according to whether the wind speed reaches the rated wind speed. Specifically, of:
在将风速数据分工况时,首先需要确定风机型号,然后确定达到额定功率时对应的额定风速值,之后将额定风速作为判断依据。以某风场2MW直驱机型风机为例,确定额定功率为2000kw、额定风速为12m/s。When dividing wind speed data into working conditions, you first need to determine the fan model, then determine the corresponding rated wind speed value when the rated power is reached, and then use the rated wind speed as the basis for judgment. Taking a 2MW direct-drive wind turbine in a wind farm as an example, it is determined that the rated power is 2000kw and the rated wind speed is 12m/s.
可选的,通过激光雷达测得风机轮毂前方的真实风速;Optionally, use lidar to measure the true wind speed in front of the turbine hub. ;
通过SCADA系统测得机舱风速,Measuring cabin wind speed through SCADA system ,
进一步地,将SCADA系统记录的该风机的风速以额定风速为界,判断所述风速数据是否小于额定风速数据;Further, the wind speed of the fan recorded by the SCADA system Taking the rated wind speed as the boundary, determine whether the wind speed data is less than the rated wind speed data;
若是,则判定所述风速数据为未变桨工况时风速;If so, it is determined that the wind speed data is the wind speed under non-pitch operating conditions;
若否,则判定所述风速数据为变桨工况时风速。If not, it is determined that the wind speed data is the wind speed under the pitching condition.
可选的,其中,所述SCADA系统中的风速数据和功率数据分辨率为1s级;将所述变桨工况数据集和未变桨工况数据集平均转化成分辨率为30s级的数据作为特征数据集。Optionally, the resolution of the wind speed data and power data in the SCADA system is 1s level; the pitch operating condition data set and the non-pitch operating condition data set are averagely converted into data with a resolution of 30s level. as a feature dataset.
S104:获取空气密度信息、扫风面积信息,将所述空气密度信息、所述扫风面积信息、所述风速数据和所述功率数据根据预设的经验公式进行计算,得到对应工况下的理论风速;S104: Obtain the air density information and sweep area information, and calculate the air density information, sweep area information, wind speed data and power data according to the preset empirical formula to obtain the Theoretical wind speed ;
该步骤用于对数据集进行处理,将SCADA风速、功率以及环境和机型参数输入经验公式得到对应的理论风速。具体,获取空气密度信息、扫风面积信息,将所述空气密度信息、所述扫风面积信息、所述风速数据和所述功率数据根据预设的经验公式进行计算,得到对应的理论风速/>;This step is used to process the data set and input SCADA wind speed, power, environmental and aircraft model parameters into the empirical formula to obtain the corresponding theoretical wind speed. . Specifically, air density information and sweep area information are obtained, and the air density information, sweep area information, wind speed data and power data are calculated according to a preset empirical formula to obtain the corresponding theoretical wind speed/ > ;
其中,在计算理论风速时,所使用的经验公式如下:Among them, when calculating the theoretical wind speed When , the empirical formula used is as follows:
; ;
其中是发电机输出功率,/>是SCADA记录的风速,/>是空气密度,/>是扫风面积。in is the generator output power,/> is the wind speed recorded by SCADA,/> is the air density,/> is the swept area.
该经验公式的推导过程如下:The derivation process of this empirical formula is as follows:
根据空气动力学理论,当空气质量以速度流过区域/>时,空气运动的功率表示为:According to aerodynamic theory, when the air mass moves at a speed Flow through area/> When , the power of air movement is expressed as:
(1) (1)
公式(1)中,是空气密度,/>是扫风面积。In formula (1), is the air density,/> is the swept area.
功率系数通常定义为:Power coefficient Usually defined as:
(2) (2)
公式(2)中,是风力发电机的机械功率。In formula (2), is the mechanical power of the wind turbine.
根据等式(2),风机机械功率和功率系数之间的关系可以写成:according to Equation (2), the relationship between fan mechanical power and power coefficient can be written as:
(3) (3)
其中: ,in: ,
为叶尖速比,/>为叶片半径,/>为发电机转速。 is the tip speed ratio,/> is the blade radius,/> is the generator speed.
考虑到风机机械功率难以直接测量,实际测量的功率是SCADA系统中的发电机功率,上式(3)可以改写为:Considering that the mechanical power of the wind turbine is difficult to measure directly, the actual measured power is the generator power in the SCADA system. The above equation (3) can be rewritten as:
(4) (4)
其中,为能量转换率。in, is the energy conversion rate.
根据空气动力学理论,通过风机叶片的气流可以分为三类,一类表示叶片前方的风速;一类表示通过叶片时的风速/>;一类表示叶片后方的风速/>。由于风速仪通常安装在风机叶片附近,因此可将SCADA风速视为/>。根据空气动力学理论,风机叶片吸收的空气动能可以表示为:According to aerodynamic theory, the airflow passing through the fan blades can be divided into three categories. One category represents the wind speed in front of the blades. ; One type represents the wind speed when passing through the blades/> ; One type represents the wind speed behind the blade/> . Since anemometers are usually installed near wind turbine blades, SCADA wind speed can be considered/> . According to aerodynamic theory, the aerodynamic energy absorbed by fan blades can be expressed as:
(5) (5)
(6) (6)
将公式(6)代入上式(5)得到:Substituting formula (6) into the above formula (5) we get:
(7) (7)
忽略风机的机械功率和发电机功率之间的差异,并将式(4)代入式(7)得:Ignoring the difference between the mechanical power of the fan and the power of the generator, and substituting equation (4) into equation (7) we get:
(8) (8)
由式(4)和式(8)可得:From formula (4) and formula (8), we can get:
(9) (9)
最终得到真实风速与SCADA风速的经验公式如下所示:The empirical formulas for the actual wind speed and SCADA wind speed are finally obtained as follows:
(10) (10)
其中是发电机输出功率,/>是SCADA记录的风速。in is the generator output power,/> is the wind speed recorded by SCADA.
S106:获取真实风速,之后通过小波变换将所述理论风速/>和所述真实风速/>分别进行分解得到对应的高频数据和低频数据,求得相应的高频残差/>和低频残差/>,将所述高频残差/>和所述低频残差/>作为目标数据集,所述风速数据为特征数据集;S106: Get the real wind speed , and then use wavelet transform to convert the theoretical wind speed/> and the actual wind speed/> Decompose them separately to obtain the corresponding high-frequency data and low-frequency data, and obtain the corresponding high-frequency residuals/> and low-frequency residual/> , convert the high-frequency residual/> and the low frequency residual/> As the target data set, the wind speed data is a feature data set;
此步骤中,通过激光雷达测得的真实风速,再通过小波变换将真实风速/>和理论风速/>分别进行分解,得到对应的高频数据/>、/>和低频数据/>、/>,然后将/>和/>相减得到真实风速高频数据与理论风速高频数据的高频残差/>,将/>和/>相减得到真实风速低频数据与理论风速低频数据的低频残差/>。In this step, the true wind speed measured by lidar , and then use wavelet transform to convert the real wind speed/> and theoretical wind speed/> Decompose them separately to obtain the corresponding high-frequency data/> ,/> and low frequency data/> ,/> , then // and/> Subtract the high-frequency residual of the real wind speed high-frequency data and the theoretical wind speed high-frequency data/> , will/> and/> Subtract to obtain the low-frequency residual of the real wind speed low-frequency data and the theoretical wind speed low-frequency data/> .
具体的,得到真实风速数据后通过小波变换将真实风速数据分解成高频数据和低频数据,并记录对应时间点,然后根据真实风速的时间点将理论风速分解成高频数据和低频数据,小波变换公式如公式(11)所示:Specifically, after obtaining the real wind speed data, the real wind speed data is decomposed into high-frequency data and low-frequency data through wavelet transformation, and the corresponding time points are recorded. Then the theoretical wind speed is decomposed into high-frequency data and low-frequency data according to the time point of the real wind speed. Wavelet The transformation formula is shown in formula (11):
(11) (11)
公式(11)中,表示母小波函数,/>分别表示伸缩因子和平移因子。/>必须满足两个条件:In formula (11), represents the mother wavelet function,/> represent the stretching factor and the translation factor respectively. /> Two conditions must be met:
第一:,/>;First: ,/> ;
第二:;second: ;
其中是/>的傅里叶变换。这时信号/>的连续小波变换由公式(12)计算:in Yes/> The Fourier transform of . Signal at this time/> The continuous wavelet transform of is calculated by formula (12):
(12) (12)
其中,为小波系数,/>是/>的共轭。in, is the wavelet coefficient,/> Yes/> of conjugation.
为了描述时间序列的小波变换,将连续小波变换离散化。令公式(11)中的/>,这时(11)式变为:To describe a time series The wavelet transform discretizes the continuous wavelet transform. Let //> in formula (11) , then equation (11) becomes:
(13) (13)
时间序列的离散小波系数计算公式为:sequentially The calculation formula of discrete wavelet coefficient is:
(14) (14)
从而,由小波系数重构原始信号的逆变换公式为:Therefore, the inverse transformation formula for reconstructing the original signal from the wavelet coefficients is:
(15) (15)
S108:将所述特征数据集作为输入,所述目标数据集作为输出,使用神经网络进行训练,用优化方法确定时间步长,根据工况和高低频数据,建立风速-残差预测模型;S108: Use the feature data set as input and the target data set as output, use neural networks for training, use optimization methods to determine the time step, and establish a wind speed-residual prediction model based on working conditions and high and low frequency data;
图3示出了本发明一种神经网络结构图,所使用的神经网络为基于自注意力机制的长短期记忆循环神经网络(Self-Attention-LSTM),该神经网络为五层结构,第一层为输入层,输入特征数据集,第二层为LSTM层,对第一层的数据依次按批数通过LSTM神经元计算得到对应的隐状态作为下一层的输入;第三层为自注意力层,将第二层输入的隐状态使用自注意力机制进行映射重组,Q、K、V映射后的维度为10;第四层为全连接层(Dense);第五层为输出层,输出模型最终计算得到的数值。此神经网络权重和偏置参数即为风速-残差之间的关系模型,将风速输入此神经网络即可得到残差值。Figure 3 shows a neural network structure diagram of the present invention. The neural network used is a long short-term memory recurrent neural network (Self-Attention-LSTM) based on the self-attention mechanism. The neural network has a five-layer structure. The first The layer is the input layer, which inputs the feature data set. The second layer is the LSTM layer. The data of the first layer is calculated through LSTM neurons in batches to obtain the corresponding hidden state as the input of the next layer; the third layer is self-attention. In the force layer, the hidden state of the second layer input is mapped and reorganized using the self-attention mechanism. The dimension after Q, K, and V mapping is 10; the fourth layer is the fully connected layer (Dense); the fifth layer is the output layer. Output the value finally calculated by the model. The weights and bias parameters of this neural network are the relationship model between wind speed and residual. If the wind speed is input into this neural network, the residual value can be obtained.
可选的,在一实施例中,所述将所述特征数据集作为输入,所述目标数据集作为输出,使用神经网络进行训练,用优化方法确定时间步长,根据工况和高低频数据,建立风速-残差预测模型包括:Optionally, in one embodiment, the feature data set is used as input and the target data set is used as output, a neural network is used for training, and an optimization method is used to determine the time step, based on working conditions and high and low frequency data. , establishing a wind speed-residual prediction model includes:
使用神经网络进行训练,建立两种工况下的风速-残差预测模型,共计四个:一是变桨工况下SCADA风速作为输入,高频残差作为输出的高频残差预测模型;二是变桨工况下SCADA风速作为输入,低频残差作为输出的低频残差预测模型;三是未变桨工况下SCADA风速作为输入,高频残差作为输出的高频残差预测模型;四是未变桨工况下SCADA风速作为输入,低频残差作为输出的低频残差预测模型。Use neural networks for training to establish wind speed-residual prediction models under two working conditions, a total of four: one is a high-frequency residual prediction model in which SCADA wind speed is used as input and high-frequency residuals are output under pitch working conditions; The second is a low-frequency residual prediction model in which SCADA wind speed is used as input under pitch changing conditions and low-frequency residuals are used as output. The third is a high-frequency residual prediction model in which SCADA wind speed is used as input under non-pitching conditions and high-frequency residuals are used as output. ; The fourth is a low-frequency residual prediction model using SCADA wind speed as input and low-frequency residual as output under non-pitch operating conditions.
在所述建立风速-残差预测模型过程中,使用伽马检测技术对训练时间步长和预测时间步长进行优化求解,采用逐层随机跳动方法对各层参数进行优化组合求解。In the process of establishing the wind speed-residual prediction model, gamma detection technology is used to optimize the training time step and prediction time step, and the layer-by-layer random jumping method is used to optimize and combine the parameters of each layer.
具体地,针对不同的风机,输入输出的时间步长通过伽马检测算法进行确定,各层的参数使用逐层随机跳动方法进行确定,但神经网络的总体结构不会改变。Specifically, for different wind turbines, the time step of the input and output is determined through the gamma detection algorithm, and the parameters of each layer are determined using the layer-by-layer random jumping method, but the overall structure of the neural network will not change.
伽马检测直接通过观测数据估计非线性模型中存在的最小均方误差,可以通过这种方式决定输入数据的参数。Gamma detection directly estimates the minimum mean square error present in the nonlinear model through observation data, and in this way the parameters of the input data can be determined.
设样本数据的形式为,/>表示第/>行输入,每个/>可能包含多个特征,用/>表示第几个特征,/>表示第/>个输出。假设向量包含影响输出的有用因素,可以把与之间的关系表示为:Suppose the form of sample data is ,/> Indicates the first/> Line input, each/> May contain multiple features, use/> Indicates which feature,/> Indicates the first/> output. Assuming that the vector contains useful factors that affect the output, the relationship between and can be expressed as:
(16) (16)
是光滑未知的回归函数,/>表示噪声随机变量。因为任何常数偏差都可以包含在未知函数中,期望/>,方差/>。 is a smooth unknown regression function,/> represents a noisy random variable. Since any constant deviation can be contained in the unknown function, the expectation/> , variance/> .
如果两个点和/>(/>)在输入空间中距离很接近,那么它们相应的输出/>和/>在输出空间中的距离也应该很接近,如果不满足,那么我们认为这种差异是由噪声/>造成的。/>函数计算输入空间中每个点的临近距离的平方值的均值,如式(17)所示:If two points and/> (/> ) are very close in the input space, then their corresponding output/> and/> The distance in the output space should also be very close, if not, then we think this difference is caused by noise/> Caused. /> The function calculates the mean of the squared value of the proximity distance of each point in the input space, as shown in Equation (17):
(17) (17)
其中,为/>的第/>临近点,即/>为最靠近/>的第/>个点,/>为/>的第/>临近点,即/>为最靠近/>的第/>个点,/>表示欧几里得距离,在输出空间中,对应的/>函数如下式:in, for/> of/> Close to the point, i.e./> is closest/> of/> points,/> for/> of/> Close to the point, i.e./> is closest/> of/> points,/> Represents the Euclidean distance, in the output space, the corresponding /> The function is as follows:
(18) (18)
伽马检测通过在点上作线性回归线来计算伽马统计量,计算公式如下:Gamma detection calculates gamma statistics by drawing a linear regression line on the points ,Calculated as follows:
(19) (19)
在式(19)中,垂直截距(当为零)描述了由/>表示的伽马统计量的值。如果/>很小,则回归函数/>存在,输出值/>在很大程度上取决于输入变量/>,输入和输出之间存在紧密的相关关系。如果/>很大,表明输入与输出无关。In equation (19), the vertical intercept (when is zero) described by/> The value of the gamma statistic represented. if/> is very small, then the regression function/> Exists, output value/> Depends a lot on the input variables/> , there is a close correlation between input and output. if/> Very large, indicating that the input has nothing to do with the output.
由于伽马统计量的值受样本数值大小的影响,考虑统计量使结果标准化。统计量/>定义为:Since the value of the gamma statistic is affected by the numerical size of the sample, consider the statistic Standardize the results. Statistics/> defined as:
(20) (20)
公式(20)中,表示输出/>的方差。接近0的/>值表明输出/>具有高度的可预测性,接近1的/>值表明输出/>是随机的,与输入变量/>无关。In formula (20), Indicates output/> Variance. Close to 0/> The value indicates the output/> Highly predictable, close to 1/> The value indicates the output/> is random, with the input variable/> Nothing to do.
伽马检测选择嵌入维数的相关说明:为嵌入维数取M时对应的伽马统计量的值,越小的伽马统计量对应的嵌入维数越合理,也就是说统计量/>的值越接近0,对应的嵌入维数越合理。因此可以按嵌入维数递增,分别进行伽马检测,使用统计量/>进行判断,选择合适的嵌入维数。Relevant instructions for selecting embedding dimensions for gamma detection: The value of the corresponding gamma statistic when M is chosen for the embedding dimension. The smaller the gamma statistic, the more reasonable the corresponding embedding dimension, that is to say, the statistic/> The closer the value is to 0, the more reasonable the corresponding embedding dimension is. Therefore, gamma detection can be performed separately according to the increasing embedding dimension, using statistics/> Make judgments and choose appropriate embedding dimensions.
各层的参数使用逐层随机跳动方法进行确定过程中,所述逐层随机跳动方法包括如下步骤:In the process of determining the parameters of each layer using the layer-by-layer random jumping method, the layer-by-layer random jumping method includes the following steps:
1)确定随机跳动的参数、跳动范围和跳动步长:1) Determine the parameters, jump range and jump step size of random jumps:
确定需要寻优的参数为第二层LSTM单元数、dropout值、激活函数;第四层全连接层单元数、激活函数;设定各参数的范围和跳动步长;激活函数为sigmod函数、tanh函数和relu函数随机选择,参数值计算公式如下:Determine the parameters that need to be optimized as the second layer LSTM unit number, dropout value, and activation function; the fourth layer fully connected layer unit number and activation function; set the range and jump step size of each parameter; the activation function is the sigmod function, tanh The function and relu function are randomly selected, and the parameter value calculation formula is as follows:
(21) (twenty one)
式中,为当前轮计算值,/>为上一轮计算值,/>为范围最小值,/>为范围最大值,/>为跳动步长,/>为在给定范围内随机选取整数的函数,/>;In the formula, Calculates the value for the current round,/> Calculated value for the previous round,/> is the minimum value of the range,/> is the maximum value of the range,/> is the jump step size,/> is a function that randomly selects integers within a given range,/> ;
2)第一轮时,在给定范围内随机初始化所选层的单元数、dropout值、映射维度和激活函数,然后将一定量的风速和残差数据集分成训练集和测试集进行训练测试并计算精度,得到对应参数的精度值,记录该参数组合和结果;然后将所选层的单元数和dropout值加上一个随机正整数或负整数乘上跳动步长得到下一组的参数值,其中,所述参数值不能超出给定范围,激活函数也随机选择,重复执行五轮;2) In the first round, the number of units, dropout value, mapping dimension and activation function of the selected layer are randomly initialized within a given range, and then a certain amount of wind speed and residual data sets are divided into training sets and test sets for training and testing. And calculate the accuracy, get the accuracy value of the corresponding parameter, record the parameter combination and result; then add a random positive or negative integer to the number of units and dropout value of the selected layer multiplied by the jump step size to get the next set of parameter values. , where the parameter values cannot exceed the given range, and the activation function is also randomly selected and executed repeatedly for five rounds;
3)计算五组随机参数的精度,选择精度最高的参数组合,然后将第四层的参数按步骤2)执行,同样计算五组参数的精度后选择精度最高的参数组合,最后得到最优的参数组合。3) Calculate the accuracy of five groups of random parameters, select the parameter combination with the highest accuracy, and then execute the parameters of the fourth layer according to step 2). Also calculate the accuracy of the five groups of parameters and select the parameter combination with the highest accuracy, and finally get the optimal Parameter combination.
S110:将一定时间步长的SCADA风速数据输入至训练好的风速-残差预测模型中,得到对应时间步长的高频残差预测值数据和低频残差预测值数据;S110: Input SCADA wind speed data of a certain time step into the trained wind speed-residual prediction model, and obtain high-frequency residual prediction value data and low-frequency residual prediction value data corresponding to the time step;
根据训练时的输入时间步长与分辨率,每隔10至20倍分辨率乘时间步长的时间对SCADA秒级风速进行一次提取。将SCADA风速根据其是否达到额定风速划分为变桨工况分速和未变桨工况风速,然后根据设定的分辨率求平均,将降低分辨率后的变桨工况风速输入到训练好的变桨工况高低频残差预测模型,得到变桨工况风速高频残差预测值和低频残差预测值;将降低分辨率后的未变桨工况风速输入到训练好的未变桨工况高低频残差预测模型,得到未变桨工况风速高频残差预测值和低频残差预测值。According to the input time step and resolution during training, the SCADA second-level wind speed is extracted every 10 to 20 times the resolution times the time step. The SCADA wind speed is divided into the pitch speed and the non-pitch wind speed according to whether it reaches the rated wind speed, and then the average is calculated according to the set resolution, and the reduced pitch wind speed is input into the training system. The high-frequency and low-frequency residual prediction model of the pitch operating condition is used to obtain the high-frequency residual prediction value and the low-frequency residual prediction value of the wind speed under the pitch operating condition; the wind speed of the unvaried pitch condition after the resolution is reduced is input to the trained unvaried wind speed. The high-frequency and low-frequency residual prediction model of the propeller operating condition is used to obtain the high-frequency residual prediction value and the low-frequency residual prediction value of the wind speed under the unchanging propeller operating condition.
S112:将所述SCADA风速数据与高频残差预测值数据和低频残差预测值数据进行线性相加,得到最终的真实风速修正值。S112: Linearly add the SCADA wind speed data, high-frequency residual prediction value data and low-frequency residual prediction value data to obtain the final true wind speed correction value.
在分别得到变桨工况风速高频残差预测值和低频残差预测值/>,以及未变桨工况风速高频残差预测值和低频残差预测值后,然后将SCADA变桨工况风速、变桨工况风速高频残差预测值和变桨工况风速低频残差预测值线性相加得到最终的真实变桨风速修正值;将SCADA未变桨工况风速、未变桨工况风速高频残差预测值和未变桨工况风速低频残差预测值线性相加得到最终的真实未变桨风速修正值。The high-frequency residual prediction values of wind speed under pitching conditions are obtained respectively. and low-frequency residual prediction value/> , and the high-frequency residual prediction value and low-frequency residual prediction value of the wind speed in the un-pitch operating condition, and then the SCADA wind speed in the pitch operating condition, the high-frequency residual prediction value of the wind speed in the pitch operating condition and the low-frequency residual prediction value of the wind speed in the pitch operating condition. The difference prediction values are linearly added to obtain the final true pitch wind speed correction value; the SCADA un-pitch operating condition wind speed, the un-pitch operating condition wind speed high-frequency residual prediction value and the un-pitch operating condition wind speed low-frequency residual prediction value are linearly added Add together to get the final true unpitched wind speed correction value.
也就是,将机舱风速加上所述对应时间步长的高频残差/>和低频残差/>得到修正的真实风速/>,即That is, the cabin wind speed Add the high-frequency residuals of the corresponding time steps/> and low-frequency residual/> Get corrected true wind speed/> ,Right now
。 (22) . (twenty two)
根据本实施例,本发明公开的一种基于机理模型和神经网络的机舱风速修正方法,选取SCADA系统中测得的风速数据和功率数据,将SCADA风速按额定风速划分为不同工况风速数据集,通过经验公式计算得到理论风速;利用小波变换计算真实风速与理论风速的高低频残差,采用神经网络建立SCADA风速与高低频残差的关系,将来自实际运行风电场的SCADA风速数据输入已训练的神经网络,获取对应高低频残差;将SCADA风速数据与残差数据一一对应,线性相加得到真实风速修正值。本发明具有计算速度快,精度高的特点且泛化性好,可广泛用于不同机型的风力发电机风速修正。According to this embodiment, the invention discloses a cabin wind speed correction method based on a mechanism model and a neural network. The wind speed data and power data measured in the SCADA system are selected, and the SCADA wind speed is divided into different working condition wind speed data sets according to the rated wind speed. , the theoretical wind speed is calculated through empirical formulas; the wavelet transform is used to calculate the high and low frequency residuals of the real wind speed and the theoretical wind speed, the neural network is used to establish the relationship between SCADA wind speed and the high and low frequency residuals, and the SCADA wind speed data from the actual operating wind farm is input. The trained neural network obtains the corresponding high and low frequency residuals; the SCADA wind speed data and the residual data are corresponded one-to-one, and the true wind speed correction value is obtained by linear addition. The invention has the characteristics of fast calculation speed, high accuracy and good generalization, and can be widely used for wind speed correction of different types of wind turbines.
图4示出了本发明一种基于机理模型和神经网络的机舱风速修正系统的模块示意图。一种基于机理模型和神经网络的机舱风速修正系统200,用于实现前述的一种基于机理模型和神经网络的机舱风速修正方法,包括:Figure 4 shows a schematic module diagram of a cabin wind speed correction system based on a mechanism model and a neural network according to the present invention. A cabin wind speed correction system 200 based on a mechanism model and a neural network is used to implement the aforementioned cabin wind speed correction method based on a mechanism model and a neural network, including:
划分模块210,用于选取SCADA系统中测得的风速数据和功率数据,将所述风速数据和功率数据根据预设规则进行划分,得到变桨工况数据集和未变桨工况数据集;The dividing module 210 is used to select the wind speed data and power data measured in the SCADA system, divide the wind speed data and power data according to the preset rules, and obtain a pitch-changing operating condition data set and a non-pitch operating condition data set;
计算模块220,用于获取空气密度信息、扫风面积信息,将所述空气密度信息、所述扫风面积信息、所述风速数据和所述功率数据根据预设的经验公式进行计算,得到对应工况下的理论风速;The calculation module 220 is used to obtain air density information and sweep area information, and calculate the air density information, sweep area information, wind speed data and power data according to preset empirical formulas to obtain the corresponding Theoretical wind speed under working conditions ;
分解模块230,用于获取真实风速,之后通过小波变换将所述理论风速/>和所述真实风速/>分别进行分解得到对应的高频数据和低频数据,求得相应的高频残差/>和低频残差/>,将所述高频残差/>和所述低频残差/>作为目标数据集,所述风速数据为特征数据集;Decomposition module 230, used to obtain the true wind speed , and then use wavelet transform to convert the theoretical wind speed/> and the actual wind speed/> Decompose them separately to obtain the corresponding high-frequency data and low-frequency data, and obtain the corresponding high-frequency residuals/> and low-frequency residual/> , convert the high-frequency residual/> and the low frequency residual/> As the target data set, the wind speed data is a feature data set;
建立模块240,用于将所述特征数据集作为输入,所述目标数据集作为输出,使用神经网络进行训练,用优化方法确定时间步长,根据工况和高低频数据,建立风速-残差预测模型;The establishment module 240 is used to take the characteristic data set as input and the target data set as output, use neural network for training, use optimization method to determine time step, and establish wind speed-residual error according to working conditions and high and low frequency data. Predictive models;
输入模块250,用于将一定时间步长的SCADA风速数据输入至训练好的风速-残差预测模型中,得到对应时间步长的高频残差预测值数据和低频残差预测值数据;The input module 250 is used to input SCADA wind speed data of a certain time step into the trained wind speed-residual prediction model to obtain high-frequency residual prediction value data and low-frequency residual prediction value data corresponding to the time step;
输出模块260,用于将所述SCADA风速数据与高频残差预测值数据和低频残差预测值数据进行线性相加,得到最终的真实风速修正值。The output module 260 is used to linearly add the SCADA wind speed data to the high-frequency residual prediction value data and the low-frequency residual prediction value data to obtain the final true wind speed correction value.
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,本发明该电子设备的具体工作过程,可以参考前述各实施例中的对应过程,以及与前述实施例相同的部分及有益效果,在此不再赘述。Those skilled in the art can clearly understand that for the convenience and simplicity of description, the specific working process of the electronic device of the present invention can be referred to the corresponding processes in the foregoing embodiments, as well as the same parts and beneficial effects as the foregoing embodiments. , which will not be described in detail here.
根据本实施例,本发明公开的一种基于机理模型和神经网络的机舱风速修正系统,选取SCADA系统中测得的风速数据和功率数据,将SCADA风速按额定风速划分为不同工况风速数据集,通过经验公式计算得到理论风速;利用小波变换计算真实风速与理论风速的高低频残差,采用神经网络建立SCADA风速与高低频残差的关系,将来自实际运行风电场的SCADA风速数据输入已训练的神经网络,获取对应高低频残差;将SCADA风速数据与残差数据一一对应,线性相加得到真实风速修正值。本发明具有计算速度快,精度高的特点且泛化性好,可广泛用于不同机型的风力发电机风速修正。According to this embodiment, the invention discloses a cabin wind speed correction system based on a mechanism model and a neural network. The wind speed data and power data measured in the SCADA system are selected, and the SCADA wind speed is divided into different working condition wind speed data sets according to the rated wind speed. , the theoretical wind speed is calculated through empirical formulas; the wavelet transform is used to calculate the high and low frequency residuals of the real wind speed and the theoretical wind speed, the neural network is used to establish the relationship between SCADA wind speed and the high and low frequency residuals, and the SCADA wind speed data from the actual operating wind farm is input. The trained neural network obtains the corresponding high and low frequency residuals; the SCADA wind speed data and the residual data are corresponded one-to-one, and the true wind speed correction value is obtained by linear addition. The invention has the characteristics of fast calculation speed, high accuracy and good generalization, and can be widely used for wind speed correction of different types of wind turbines.
在本申请所提供的几个实施例中,应该理解到,所揭露的方法和系统,可以通过其它的方式实现。以上所描述的设备实施例仅仅是示意性的,例如,所述模块的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,如:多个模块或组件可以结合,或可以集成到另一个系统,或一些特征可以忽略,或不执行。另外,所显示或讨论的各组成部分相互之间的耦合、或直接耦合、或通信连接可以是通过一些接口,设备或单元的间接耦合或通信连接,可以是电性的、机械的或其它形式的。In the several embodiments provided in this application, it should be understood that the disclosed methods and systems can be implemented in other ways. The device embodiments described above are only illustrative. For example, the division of modules is only a logical function division. In actual implementation, there may be other division methods, such as: multiple modules or components can be combined, or can be integrated into another system, or some features can be ignored, or not implemented. In addition, the coupling, direct coupling, or communication connection between the components shown or discussed may be through some interfaces, and the indirect coupling or communication connection of the devices or units may be electrical, mechanical, or other forms. of.
上述作为分离部件说明的模块可以是、或也可以不是物理上分开的,作为模块显示的部件可以是、或也可以不是物理单元;既可以位于一个地方,也可以分布到多个网络单元上;可以根据实际的需要选择其中的部分或全部单元来实现本实施例方案的目的。The modules described above as separate components may or may not be physically separated, and the components shown as modules may or may not be physical units; they may be located in one place or distributed to multiple network units; Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of this embodiment.
另外,在本发明各实施例中的各功能单元可以全部集成在一个处理单元中,也可以是各单元分别单独作为一个单元,也可以两个或两个以上单元集成在一个单元中;上述集成的单元既可以采用硬件的形式实现,也可以采用硬件加软件功能单元的形式实现。In addition, each functional unit in each embodiment of the present invention can be all integrated into one processing unit, or each unit can be separately used as a unit, or two or more units can be integrated into one unit; the above-mentioned integration The unit can be implemented in the form of hardware or in the form of hardware plus software functional units.
本领域普通技术人员可以理解:实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成,前述的程序可以存储于计算机可读取存储介质中,该程序在执行时,执行包括上述方法实施例的步骤;而前述的存储介质包括:移动存储设备、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。Those of ordinary skill in the art can understand that all or part of the steps to implement the above method embodiments can be completed through hardware related to program instructions. The aforementioned program can be stored in a computer-readable storage medium. When the program is executed, the execution includes: The steps of the above method embodiment; and the aforementioned storage media include: mobile storage devices, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic disks or optical disks, etc. The medium on which program code is stored.
或者,本发明上述集成的单元如果以软件功能模块的形式实现并作为独立的产品销售或使用时,也可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明实施例的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机、服务器、或者网络设备等)执行本发明各个实施例所述方法的全部或部分。而前述的存储介质包括:移动存储设备、ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。Alternatively, if the above-mentioned integrated unit of the present invention is implemented in the form of a software function module and sold or used as an independent product, it can also be stored in a computer-readable storage medium. Based on this understanding, the technical solutions of the embodiments of the present invention can be embodied in the form of software products in essence or those that contribute to the existing technology. The computer software products are stored in a storage medium and include a number of instructions to A computer device (which may be a personal computer, a server, a network device, etc.) is caused to execute all or part of the methods described in various embodiments of the present invention. The aforementioned storage media include: mobile storage devices, ROM, RAM, magnetic disks or optical disks and other media that can store program codes.
虽然已经通过示例对本发明的一些特定实施例进行了详细说明,但是本领域的技术人员应该理解,以上示例仅是为了进行说明,而不是为了限制本发明的范围。本领域的技术人员还应理解,可以对实施例进行多种修改而不脱离本发明的范围和精神。本发明开的范围由所附权利要求来限定。Although some specific embodiments of the invention have been described in detail by way of examples, those skilled in the art will understand that the above examples are for illustration only and are not intended to limit the scope of the invention. It will also be understood by those skilled in the art that various modifications may be made to the embodiments without departing from the scope and spirit of the invention. The scope of the invention is defined by the appended claims.
Claims (7)
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311047504.4A CN116757101B (en) | 2023-08-21 | 2023-08-21 | Cabin wind speed correction method and system based on mechanism model and neural network |
NL2036572A NL2036572A (en) | 2023-08-21 | 2023-12-18 | Nacelle wind speed correction method and system based on mechanism model and neural network |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311047504.4A CN116757101B (en) | 2023-08-21 | 2023-08-21 | Cabin wind speed correction method and system based on mechanism model and neural network |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116757101A CN116757101A (en) | 2023-09-15 |
CN116757101B true CN116757101B (en) | 2023-11-07 |
Family
ID=87959518
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311047504.4A Active CN116757101B (en) | 2023-08-21 | 2023-08-21 | Cabin wind speed correction method and system based on mechanism model and neural network |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN116757101B (en) |
NL (1) | NL2036572A (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102562469A (en) * | 2011-12-27 | 2012-07-11 | 华北电力大学 | Short-term wind driven generator output power predicting method based on correction algorithm |
CN110414045A (en) * | 2019-06-18 | 2019-11-05 | 东华大学 | Short-term wind speed prediction method based on VMD-GRU |
CN113537582A (en) * | 2021-07-02 | 2021-10-22 | 东北电力大学 | An ultra-short-term prediction method of photovoltaic power based on short-wave radiation correction |
CN114091377A (en) * | 2022-01-21 | 2022-02-25 | 南京理工大学 | Calculation method of wake wind speed of dynamic double Gaussian wind turbine based on spatial variation |
CN115374710A (en) * | 2022-09-19 | 2022-11-22 | 兰州理工大学 | Wind speed prediction method based on improved seagull optimization algorithm and multi-core extreme learning machine |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2017527059A (en) * | 2014-08-29 | 2017-09-14 | ズーナム エアロ,インコーポレイテッド | System and method for implementing a regional air traffic network using a hybrid electric airplane |
CN110261394B (en) * | 2019-06-24 | 2022-09-16 | 内蒙古工业大学 | Online real-time diagnosis system and method for damage of fan blade |
CN110566406B (en) * | 2019-10-16 | 2020-08-04 | 上海海事大学 | Robust control system and method for real-time pitch pitch of wind turbine based on reinforcement learning |
KR20210006874A (en) * | 2020-12-30 | 2021-01-19 | 정인우 | Kalman Filter and Deep Reinforcement Learning based Wind Turbine Yaw Misalignmnet Control Method |
CN113177187A (en) * | 2021-03-04 | 2021-07-27 | 明阳智慧能源集团股份公司 | Equivalent laser radar wind speed calculation method based on long-term and short-term memory neural network |
CN116362144A (en) * | 2021-12-22 | 2023-06-30 | 北京金风科创风电设备有限公司 | Laser radar wind speed correction method, device and equipment |
CN115345076B (en) * | 2022-08-19 | 2023-09-12 | 清华大学 | Wind speed correction processing method and device |
-
2023
- 2023-08-21 CN CN202311047504.4A patent/CN116757101B/en active Active
- 2023-12-18 NL NL2036572A patent/NL2036572A/en unknown
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102562469A (en) * | 2011-12-27 | 2012-07-11 | 华北电力大学 | Short-term wind driven generator output power predicting method based on correction algorithm |
CN110414045A (en) * | 2019-06-18 | 2019-11-05 | 东华大学 | Short-term wind speed prediction method based on VMD-GRU |
CN113537582A (en) * | 2021-07-02 | 2021-10-22 | 东北电力大学 | An ultra-short-term prediction method of photovoltaic power based on short-wave radiation correction |
CN114091377A (en) * | 2022-01-21 | 2022-02-25 | 南京理工大学 | Calculation method of wake wind speed of dynamic double Gaussian wind turbine based on spatial variation |
CN115374710A (en) * | 2022-09-19 | 2022-11-22 | 兰州理工大学 | Wind speed prediction method based on improved seagull optimization algorithm and multi-core extreme learning machine |
Non-Patent Citations (2)
Title |
---|
Power Data Preprocessing Method of Mountain Wind Farm Based on POT-DBSCAN;Zhao Xiao等;《Energy Engineering》;第549-563页 * |
基于运行规律和TICC算法的风电SCADA高维时序数据聚类方法;肖钊;《机械工程学报》;第58卷(第23期);第196-206页 * |
Also Published As
Publication number | Publication date |
---|---|
CN116757101A (en) | 2023-09-15 |
NL2036572A (en) | 2024-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kisvari et al. | Wind power forecasting–A data-driven method along with gated recurrent neural network | |
CN108334672B (en) | Maximum wind energy capture method for variable speed wind turbines based on effective wind speed estimation | |
CN102562469B (en) | Short-term wind driven generator output power predicting method based on correction algorithm | |
CN107507097A (en) | A kind of short-term wind power prediction method | |
CN111881572A (en) | Wind power plant collaborative yaw intelligent control method based on multi-objective optimization | |
Xu et al. | Correlation based neuro-fuzzy Wiener type wind power forecasting model by using special separate signals | |
CN106979126A (en) | Wind power generating set high wind speed section effective wind speed method of estimation based on SVR | |
Gao et al. | Data-driven yaw misalignment correction for utility-scale wind turbines | |
CN111931983B (en) | Precipitation prediction method and system | |
Dimitrov et al. | From SCADA to lifetime assessment and performance optimization: how to use models and machine learning to extract useful insights from limited data | |
CN115526376A (en) | Generative Adversarial Network Ultra-short-term Wind Power Prediction Method Based on Multi-feature Fusion | |
CN115013258A (en) | Method, device and server for intelligent soft measurement of wind speed in front of wind turbine | |
CN115660232A (en) | Ultra-short-term prediction method, device and system for wind power | |
CN113610285A (en) | Power prediction method for distributed wind power | |
Andersen et al. | Global trends in the performance of large wind farms based on high-fidelity simulations | |
CN114320770B (en) | Method and device for measuring effective wind speed of wind turbine generator | |
CN111598313A (en) | Wind power prediction method | |
CN113343562B (en) | Fan power prediction method and system based on hybrid modeling strategy | |
CN116757101B (en) | Cabin wind speed correction method and system based on mechanism model and neural network | |
Kim et al. | Accuracy of a short-term wind power forecasting model based on deep learning using LiDAR-SCADA integration: A case study of the 400-MW Anholt offshore wind farm. | |
CN117892099B (en) | Wind turbine service effectiveness evaluation method and system based on probabilistic neural network | |
CN112561149B (en) | Wind power prediction method based on dynamic downscaling data fusion and statistical model | |
CN116882269A (en) | Wind power plant error correction model construction method and device and wind power plant control system | |
Bertelè et al. | Brief communication: Wind inflow observation from load harmonics–wind tunnel validation of the rotationally symmetric formulation | |
CN113312836B (en) | Short-term wind speed prediction method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
OL01 | Intention to license declared | ||
OL01 | Intention to license declared |