CN114583767A - Data-driven wind power plant frequency modulation response characteristic modeling method and system - Google Patents
Data-driven wind power plant frequency modulation response characteristic modeling method and system Download PDFInfo
- Publication number
- CN114583767A CN114583767A CN202210229805.8A CN202210229805A CN114583767A CN 114583767 A CN114583767 A CN 114583767A CN 202210229805 A CN202210229805 A CN 202210229805A CN 114583767 A CN114583767 A CN 114583767A
- Authority
- CN
- China
- Prior art keywords
- data
- frequency modulation
- neural network
- transfer function
- working condition
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000004044 response Effects 0.000 title claims abstract description 94
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000012546 transfer Methods 0.000 claims abstract description 71
- 238000003062 neural network model Methods 0.000 claims abstract description 40
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 23
- 238000000691 measurement method Methods 0.000 claims abstract description 14
- 238000013528 artificial neural network Methods 0.000 claims description 22
- 238000004364 calculation method Methods 0.000 claims description 13
- 238000005259 measurement Methods 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 238000012360 testing method Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 5
- 238000011156 evaluation Methods 0.000 abstract description 4
- 238000007781 pre-processing Methods 0.000 abstract description 3
- 238000012549 training Methods 0.000 abstract description 2
- 238000011157 data evaluation Methods 0.000 abstract 1
- 230000006870 function Effects 0.000 description 58
- 230000008569 process Effects 0.000 description 13
- 238000010586 diagram Methods 0.000 description 10
- 230000008859 change Effects 0.000 description 7
- 239000011159 matrix material Substances 0.000 description 6
- 238000011478 gradient descent method Methods 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000021615 conjugation Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000014509 gene expression Effects 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 230000006386 memory function Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for AC mains or AC distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/48—Controlling the sharing of the in-phase component
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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
-
- 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/048—Activation functions
-
- 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for AC mains or AC distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/381—Dispersed generators
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
- H02J2300/20—The dispersed energy generation being of renewable origin
- H02J2300/28—The renewable source being wind energy
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E10/00—Energy generation through renewable energy sources
- Y02E10/70—Wind energy
- Y02E10/76—Power conversion electric or electronic aspects
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Data Mining & Analysis (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Artificial Intelligence (AREA)
- Computational Linguistics (AREA)
- Mathematical Optimization (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Power Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Molecular Biology (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- Medical Informatics (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Wind Motors (AREA)
Abstract
Description
技术领域technical field
本发明涉及风电场调频响应特性技术领域,特别是涉及一种数据驱动的风电场调频响应特性建模方法及系统。The invention relates to the technical field of frequency modulation response characteristics of wind farms, in particular to a data-driven method and system for modeling frequency modulation response characteristics of wind farms.
背景技术Background technique
风电场内频率稳定的研究,一个重要的问题就是建立风场调频响应特性的模型,一般来说,建模方法分为两种,一种为基于组件的建模,一种是基于数据的建模。基于组件的建模往往需要风电机组内部组件的参数和运行时的动力学数学模型,而实际上由于保密协议的存在,风机厂商不能够提供详细的风机信息,这使得组件建模有一定的局限性。基于风场运行数据的建模关注研究对象的输入和输出,避免了对风电机组组件的物理描述和获取相应的参数,且模型也能达到满意的精度。因此数据驱动的建模方法被越来越多的应用到风电机组及风电场建模过程中。In the research of frequency stability in wind farms, an important issue is to establish a model of the frequency modulation response characteristics of the wind farm. Generally speaking, there are two modeling methods, one is component-based modeling, and the other is data-based modeling. mold. Component-based modeling often requires the parameters of the internal components of the wind turbine and the dynamic mathematical model during operation. In fact, due to the existence of the confidentiality agreement, the wind turbine manufacturer cannot provide detailed wind turbine information, which makes the component modeling have certain limitations. sex. The modeling based on wind farm operation data focuses on the input and output of the research object, avoiding the physical description of the wind turbine components and obtaining the corresponding parameters, and the model can also achieve satisfactory accuracy. Therefore, data-driven modeling methods are more and more applied to the modeling process of wind turbines and wind farms.
带外部输入的非线性自回归神经网络模型能够很好的表征非线性时序系统。风电场并网点调频响应就是一个时序的非线性过程。根据现有可得的相关文献,与本发明最相近的现有风电场调频响应特性建模方法,传递函数方法往往是基于组件的,并且考虑单工况的较多,非线性自回归神经网络多用于风电场风速、功率的预测,应用风电场调频建模的极少。我们知道,当来自并网点的调频指令下达时,风电场不一定处于最大功率运行状态,若在其它有功功率运行状态下进行调频,单一模型可能无法表征。这可能导致针对单一模型进行的风电场频率响应评估不准确,调频效果不理想等问题。The nonlinear autoregressive neural network model with external input can well characterize the nonlinear time series system. The frequency modulation response of wind farm grid connection point is a nonlinear process of time series. According to the existing relevant literature, the most similar to the present invention is the existing wind farm frequency modulation response characteristic modeling method, the transfer function method is often based on components, and considers many single working conditions, nonlinear autoregressive neural network It is mostly used for wind speed and power prediction of wind farms, and there are very few applications of wind farm frequency modulation modeling. We know that when the frequency regulation command from the grid connection point is issued, the wind farm is not necessarily in the maximum power operating state. If the frequency regulation is performed in other active power operating states, a single model may not be able to characterize it. This may lead to inaccurate evaluation of the frequency response of the wind farm based on a single model, and unsatisfactory frequency modulation effect.
发明内容SUMMARY OF THE INVENTION
为了克服现有技术的不足,本发明的目的是提供一种数据驱动的风电场调频响应特性建模方法及系统。In order to overcome the deficiencies of the prior art, the purpose of the present invention is to provide a data-driven method and system for modeling the frequency modulation response characteristics of a wind farm.
为实现上述目的,本发明提供了如下方案:For achieving the above object, the present invention provides following scheme:
一种数据驱动的风电场调频响应特性建模方法,包括:A data-driven modeling method for frequency modulation response characteristics of wind farms, comprising:
基于阶跃响应动态性能指标求解算法对风电场调频响应特性实测数据进行分析和预处理,得到各个工况下处理后的数据;Based on the step response dynamic performance index solution algorithm, the measured data of the frequency modulation response characteristics of the wind farm are analyzed and preprocessed, and the processed data under each working condition are obtained;
根据所述处理后的数据对每个工况建立传递函数模型,并利用间隙测度方法测量各模型间的间隙值,以确定非线性自回归神经网络模型表征的工况区域;Establish a transfer function model for each working condition according to the processed data, and measure the gap value between the models by using the gap measurement method to determine the working condition area represented by the nonlinear autoregressive neural network model;
根据所述间隙值对工况的调频数据进行合并,并根据合并后的数据对所述非线性自回归神经网络模型进行训练,得到训练好的非线性自回归神经网络模型。The frequency modulation data of the working conditions are combined according to the gap value, and the nonlinear autoregressive neural network model is trained according to the combined data to obtain a trained nonlinear autoregressive neural network model.
优选地,所述基于阶跃响应动态性能指标求解算法对风电场调频响应特性实测数据进行分析和预处理,得到各个工况下处理后的数据,包括:Preferably, the step-response-based dynamic performance index solution algorithm analyzes and preprocesses the measured data of the frequency modulation response characteristics of the wind farm to obtain processed data under various operating conditions, including:
根据风电场调频时的初始有功功率对所述风电场调频响应特性实测数据进行划分,得到多个工况;Divide the measured data of the frequency modulation response characteristic of the wind farm according to the initial active power during the frequency modulation of the wind farm to obtain a plurality of working conditions;
获取所述风电场调频响应特性实测数据的数据集的功率曲线,并计算所述功率曲线的各个数据点的上峰值集合和下峰值集合;Obtaining the power curve of the data set of the measured data of the frequency modulation response characteristic of the wind farm, and calculating the upper peak set and the lower peak set of each data point of the power curve;
计算所述上峰值集合的最大值和所述下峰值集合的最小值;calculating the maximum value of the upper peak set and the minimum value of the lower peak set;
基于所述数据集中的数据点,计算误差带上界和误差带下界;based on the data points in the data set, calculating an upper error band and a lower error band;
判断所述误差带上界是否大于或者等于所述上峰值集合的最大值,且所述误差带下界是否大于或者等于所述下峰值集合的最小值,若判断结果均为是,则根据采样间隔确定调节时间;Judging whether the upper bound of the error band is greater than or equal to the maximum value of the upper peak set, and whether the lower bound of the error band is greater than or equal to the minimum value of the lower peak set, if the judgment results are all yes, according to the sampling interval determine the adjustment time;
判断所述风电场调频响应特性实测数据是否为频率阶跃下扰,若是,则根据所述上峰值集合的最大值确定超调量;Determine whether the measured data of the frequency modulation response characteristic of the wind farm is a frequency step underdisturbance, and if so, determine the overshoot amount according to the maximum value of the upper peak set;
根据所述数据集确定上升时间或下降时间。Rise time or fall time is determined from the data set.
优选地,所述根据所述处理后的数据对每个工况建立传递函数模型,包括:Preferably, establishing a transfer function model for each operating condition according to the processed data includes:
分别根据各个工况下的所述处理后的数据中的频率变化量为输入,功率变化量为输出,构建初始函数模型;According to the frequency variation in the processed data under each working condition as the input, and the power variation as the output, an initial function model is constructed;
设定传递函数的分子阶数和分母阶数;Set the numerator order and denominator order of the transfer function;
对所述初始传递函数模型进行模型辨识,并调整所述分子阶数和所述分母阶数,得到最优的所述传递函数模型。Perform model identification on the initial transfer function model, and adjust the numerator order and the denominator order to obtain the optimal transfer function model.
优选地,所述利用间隙测度方法测量各模型间的间隙值,以确定非线性自回归神经网络模型表征的工况区域,包括:Preferably, the use of the gap measurement method to measure the gap value between the models to determine the working condition area represented by the nonlinear autoregressive neural network model includes:
根据不同的所述传递函数模型的正交投影矩阵确定间隙度量公式;Determine the gap measurement formula according to the orthogonal projection matrices of the different transfer function models;
根据所述间隙度量公式计算所述间隙值;所述间隙值用于确定所述工况区域。The gap value is calculated according to the gap metric formula; the gap value is used to determine the operating condition area.
优选地,所述根据所述间隙值对工况的调频数据进行合并,包括:Preferably, the merging of the frequency modulation data of the working conditions according to the gap value includes:
计算不同工况的所述传递函数模型之间的间隙值与0之间的距离,得到第一距离;Calculate the distance between the gap value and 0 between the transfer function models of different working conditions to obtain the first distance;
计算不同工况的所述传递函数模型之间的间隙值与1之间的距离,得到第二距离;Calculate the distance between the gap value and 1 between the transfer function models of different working conditions to obtain the second distance;
判断所述第一距离是否小于或者等于所述第二距离,若是,则按时序合并这两个工况的调频数据,得到合并后的数据;若否,则选取任意一个工况的调频数据与已有的调频数据进行时序合并,得到合并后的数据。Determine whether the first distance is less than or equal to the second distance, and if so, combine the FM data of the two working conditions according to the time series to obtain the combined data; if not, select the FM data of any working condition and The existing FM data is time-series merged to obtain merged data.
优选地,所述根据合并后的数据对所述非线性自回归神经网络模型进行训练,得到训练好的非线性自回归神经网络模型,包括:Preferably, the nonlinear autoregressive neural network model is trained according to the combined data to obtain a trained nonlinear autoregressive neural network model, including:
构建初始神经网络;Build the initial neural network;
基于列文伯格-马尔夸克算法,根据所述合并后的数据对所述初始神经网络进行训练;Based on the Levenberg-Marquark algorithm, the initial neural network is trained according to the combined data;
利用每一个工况的调频数据对训练后的神经网络进行测试,并根据测试结果确定所述非线性自回归神经网络模型。The trained neural network is tested by using the frequency modulation data of each working condition, and the nonlinear autoregressive neural network model is determined according to the test results.
一种数据驱动的风电场调频响应特性建模系统,包括:A data-driven wind farm frequency modulation response characteristic modeling system, comprising:
数据处理模块,用于基于阶跃响应动态性能指标求解算法对风电场调频响应特性实测数据进行分析和预处理,得到各个工况下处理后的数据;The data processing module is used to analyze and preprocess the measured data of the frequency modulation response characteristics of the wind farm based on the step response dynamic performance index solution algorithm, and obtain the processed data under each working condition;
传递函数构建模块,用于根据所述处理后的数据对每个工况建立传递函数模型,并利用间隙测度方法测量各模型间的间隙值,以确定非线性自回归神经网络模型表征的工况区域;The transfer function building module is used to establish a transfer function model for each working condition according to the processed data, and use the gap measurement method to measure the gap value between the models to determine the working condition represented by the nonlinear autoregressive neural network model area;
神经网络建模模块,用于根据所述间隙值对工况的调频数据进行合并,并根据合并后的数据对所述非线性自回归神经网络模型进行训练,得到训练好的非线性自回归神经网络模型。A neural network modeling module is used to merge the frequency modulation data of the working condition according to the gap value, and train the nonlinear autoregressive neural network model according to the merged data to obtain a trained nonlinear autoregressive neural network. network model.
优选地,所述数据处理模块具体包括:Preferably, the data processing module specifically includes:
工况划分单元,用于根据风电场调频时的初始有功功率对所述风电场调频响应特性实测数据进行划分,得到多个工况;A working condition dividing unit, configured to divide the measured data of the frequency modulation response characteristic of the wind farm according to the initial active power during frequency modulation of the wind farm to obtain a plurality of working conditions;
峰值计算单元,用于获取所述风电场调频响应特性实测数据的数据集的功率曲线,并计算所述功率曲线的各个数据点的上峰值集合和下峰值集合;a peak value calculation unit, configured to obtain the power curve of the data set of the measured data of the frequency modulation response characteristic of the wind farm, and calculate the upper peak value set and the lower peak value set of each data point of the power curve;
最值计算单元,用于计算所述上峰值集合的最大值和所述下峰值集合的最小值;a maximum value calculation unit, configured to calculate the maximum value of the upper peak set and the minimum value of the lower peak set;
上下界计算单元,用于基于所述数据集中的数据点,计算误差带上界和误差带下界;an upper and lower bound calculation unit, configured to calculate the upper bound of the error band and the lower bound of the error band based on the data points in the data set;
第一判断单元,用于判断所述误差带上界是否大于或者等于所述上峰值集合的最大值,且所述误差带下界是否大于或者等于所述下峰值集合的最小值,若判断结果均为是,则根据采样间隔确定调节时间;A first judging unit, configured to judge whether the upper bound of the error band is greater than or equal to the maximum value of the upper peak set, and whether the lower bound of the error band is greater than or equal to the minimum value of the lower peak set, if the judgment results are all If yes, determine the adjustment time according to the sampling interval;
第二判断单元,用于判断所述风电场调频响应特性实测数据是否为频率阶跃下扰,若是,则根据所述上峰值集合的最大值确定超调量;a second judging unit, configured to judge whether the measured data of the frequency modulation response characteristic of the wind farm is a frequency step underdisturbance, and if so, determine the overshoot amount according to the maximum value of the upper peak set;
时间确定单元,用于根据所述数据集确定上升时间或下降时间。A time determination unit, configured to determine a rise time or a fall time according to the data set.
优选地,所述传递函数构建模块具体包括:Preferably, the transfer function building module specifically includes:
初始模型构建单元,用于分别根据各个工况下的所述处理后的数据中的频率变化量为输入,功率变化量为输出,构建初始函数模型;an initial model building unit, configured to construct an initial function model according to the frequency variation in the processed data under each working condition as the input and the power variation as the output;
阶数设定单元,用于设定传递函数的分子阶数和分母阶数;Order setting unit, used to set the numerator order and denominator order of the transfer function;
模型确定单元,用于对所述初始传递函数模型进行模型辨识,并调整所述分子阶数和所述分母阶数,得到最优的所述传递函数模型。A model determining unit, configured to perform model identification on the initial transfer function model, and adjust the numerator order and the denominator order to obtain the optimal transfer function model.
优选地,所述传递函数构建模块具体包括:Preferably, the transfer function building module specifically includes:
公式确定单元,用于根据不同的所述传递函数模型的正交投影矩阵确定间隙度量公式;a formula determination unit, used for determining a gap measurement formula according to the orthogonal projection matrices of the different transfer function models;
间隙计算单元,用于根据所述间隙度量公式计算所述间隙值;所述间隙值用于确定所述工况区域。A clearance calculation unit, configured to calculate the clearance value according to the clearance measurement formula; the clearance value is used to determine the working condition area.
根据本发明提供的具体实施例,本发明公开了以下技术效果:According to the specific embodiments provided by the present invention, the present invention discloses the following technical effects:
本发明提供了一种数据驱动的风电场调频响应特性建模方法及系统,所述方法包括:基于阶跃响应动态性能指标求解算法对风电场调频响应特性实测数据进行分析和预处理,得到各个工况下处理后的数据;根据所述处理后的数据对每个工况建立传递函数模型,并利用间隙测度方法测量各模型间的间隙值,以确定非线性自回归神经网络模型表征的工况区域;根据所述间隙值对工况的调频数据进行合并,并根据合并后的数据对所述非线性自回归神经网络模型进行训练,得到训练好的非线性自回归神经网络模型。本发明利用风电场实际调频响应数据,设计了能够良好表征不同工况下风电场并网点调频响应特性的建模方案,能够提高风电场频率响应评估的准确性和调频效果。The present invention provides a data-driven method and system for modeling frequency modulation response characteristics of wind farms. The method includes: analyzing and preprocessing the measured data of frequency modulation response characteristics of wind farms based on a step response dynamic performance index solution algorithm to obtain each The processed data under the working conditions; the transfer function model is established for each working condition according to the processed data, and the gap value between the models is measured by the gap measurement method to determine the working conditions represented by the nonlinear autoregressive neural network model. The frequency modulation data of the working condition is combined according to the gap value, and the nonlinear autoregressive neural network model is trained according to the combined data to obtain a trained nonlinear autoregressive neural network model. The invention uses the actual frequency modulation response data of the wind farm to design a modeling scheme that can well characterize the frequency modulation response characteristics of the grid connection point of the wind farm under different working conditions, and can improve the accuracy of the frequency response evaluation of the wind farm and the frequency modulation effect.
附图说明Description of drawings
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the accompanying drawings required in the embodiments will be briefly introduced below. Obviously, the drawings in the following description are only some of the present invention. In the embodiments, for those of ordinary skill in the art, other drawings can also be obtained according to these drawings without creative labor.
图1为本发明提供的实施例中的风电场调频响应特性建模方法的方法流程图;1 is a method flowchart of a method for modeling a frequency modulation response characteristic of a wind farm in an embodiment provided by the present invention;
图2为本发明提供的实施例中的传递函数模型的输出与实际数据的第一幅输出对比图;2 is a comparison diagram of the output of the transfer function model in the embodiment provided by the present invention and the first output of actual data;
图3为本发明提供的实施例中的传递函数模型的输出与实际数据的第二幅输出对比图;3 is a comparison diagram of the output of the transfer function model in the embodiment provided by the present invention and the second output of actual data;
图4为本发明提供的实施例中的传递函数模型的输出与实际数据的第三幅输出对比图;Fig. 4 is the output of the transfer function model in the embodiment provided by the present invention and the third output comparison diagram of actual data;
图5为本发明提供的实施例中的传递函数模型的输出与实际数据的第四幅输出对比图;Fig. 5 is the output of the transfer function model in the embodiment provided by the present invention and the fourth output comparison diagram of actual data;
图6为本发明提供的实施例中的带外部输入的非线性自回归神经网络结构图;6 is a structural diagram of a nonlinear autoregressive neural network with external input in an embodiment provided by the present invention;
图7为本发明提供的实施例中的多工况数据训练的非线性自回归神经网络模型被单工况的频率干扰数据测试的第一幅结果图;7 is a first result diagram of a nonlinear autoregressive neural network model trained with multi-operating condition data in an embodiment provided by the present invention tested by frequency interference data of a single operating condition;
图8为本发明提供的实施例中的多工况数据训练的非线性自回归神经网络模型被单工况的频率干扰数据测试的第二幅结果图;8 is a second result diagram of a nonlinear autoregressive neural network model trained with multi-operating condition data in an embodiment provided by the present invention tested by frequency interference data of a single operating condition;
图9为本发明提供的实施例中的多工况数据训练的非线性自回归神经网络模型被单工况的频率干扰数据测试的第三幅结果图;9 is a third result diagram of a nonlinear autoregressive neural network model trained on multi-operating condition data in an embodiment provided by the present invention tested by frequency interference data of a single operating condition;
图10为本发明提供的实施例中的多工况数据训练的非线性自回归神经网络模型被单工况的频率干扰数据测试的第四幅结果图。FIG. 10 is a fourth result diagram of a nonlinear autoregressive neural network model trained with data of multiple operating conditions tested by frequency interference data of a single operating condition in the embodiment provided by the present invention.
具体实施方式Detailed ways
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, but not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.
在本文中提及“实施例”意味着,结合实施例描述的特定特征、结构或特性可以包含在本申请的至少一个实施例中。在说明书中的各个位置出现该短语并不一定均是指相同的实施例,也不是与其它实施例互斥的独立的或备选的实施例。本领域技术人员显式地和隐式地理解的是,本文所描述的实施例可以与其它实施例相结合。Reference herein to an "embodiment" means that a particular feature, structure, or characteristic described in connection with the embodiment can be included in at least one embodiment of the present application. The appearances of the phrase in various places in the specification are not necessarily all referring to the same embodiment, nor a separate or alternative embodiment that is mutually exclusive of other embodiments. It is explicitly and implicitly understood by those skilled in the art that the embodiments described herein may be combined with other embodiments.
本申请的说明书和权利要求书及所述附图中的术语“第一”、“第二”、“第三”和“第四”等是用于区别不同对象,而不是用于描述特定顺序。此外,术语“包括”和“具有”以及它们任何变形,意图在于覆盖不排他的包含。例如包含了一系列步骤、过程、方法等没有限定于已列出的步骤,而是可选地还包括没有列出的步骤,或可选地还包括对于这些过程、方法、产品或设备固有的其它步骤元。The terms "first", "second", "third" and "fourth" in the description and claims of the present application and the drawings are used to distinguish different objects, rather than to describe a specific order . Furthermore, the terms "comprising" and "having", and any variations thereof, are intended to cover non-exclusive inclusion. For example, including a series of steps, processes, methods, etc. is not limited to the listed steps, but optionally also includes unlisted steps, or optionally also includes inherent to these processes, methods, products or devices. other steps.
本发明的目的是提供一种数据驱动的风电场调频响应特性建模方法及系统,能够提高风电场频率响应评估的准确性和调频效果。The purpose of the present invention is to provide a data-driven wind farm frequency modulation response characteristic modeling method and system, which can improve the accuracy of wind farm frequency response evaluation and the frequency modulation effect.
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。In order to make the above objects, features and advantages of the present invention more clearly understood, the present invention will be described in further detail below with reference to the accompanying drawings and specific embodiments.
图1为本发明提供的实施例中的风电场调频响应特性建模方法的方法流程图,如图1所示,本发明提供了一种数据驱动的风电场调频响应特性建模方法,包括:FIG. 1 is a method flowchart of a method for modeling frequency modulation response characteristics of a wind farm in an embodiment provided by the present invention. As shown in FIG. 1 , the present invention provides a data-driven method for modeling frequency modulation response characteristics of wind farms, including:
步骤100:基于阶跃响应动态性能指标求解算法对风电场调频响应特性实测数据进行分析和预处理,得到各个工况下处理后的数据;Step 100: analyze and preprocess the measured data of the frequency modulation response characteristic of the wind farm based on the step response dynamic performance index solution algorithm, and obtain the processed data under each working condition;
步骤200:根据所述处理后的数据对每个工况建立传递函数模型,并利用间隙测度方法测量各模型间的间隙值,以确定非线性自回归神经网络模型表征的工况区域;Step 200: establishing a transfer function model for each operating condition according to the processed data, and measuring the gap value between the models by using the gap measurement method to determine the operating condition area represented by the nonlinear autoregressive neural network model;
步骤300:根据所述间隙值对工况的调频数据进行合并,并根据合并后的数据对所述非线性自回归神经网络模型进行训练,得到训练好的非线性自回归神经网络模型。Step 300: Combine the frequency modulation data of the working conditions according to the gap value, and train the nonlinear autoregressive neural network model according to the combined data to obtain a trained nonlinear autoregressive neural network model.
优选地,所述步骤100包括:Preferably, the
根据风电场调频时的初始有功功率对所述风电场调频响应特性实测数据进行划分,得到多个工况;Divide the measured data of the frequency modulation response characteristic of the wind farm according to the initial active power during the frequency modulation of the wind farm to obtain a plurality of working conditions;
获取所述风电场调频响应特性实测数据的数据集的功率曲线,并计算所述功率曲线的各个数据点的上峰值集合和下峰值集合;Obtaining the power curve of the data set of the measured data of the frequency modulation response characteristic of the wind farm, and calculating the upper peak set and the lower peak set of each data point of the power curve;
计算所述上峰值集合的最大值和所述下峰值集合的最小值;calculating the maximum value of the upper peak set and the minimum value of the lower peak set;
基于所述数据集中的数据点,计算误差带上界和误差带下界;based on the data points in the data set, calculating an upper error band and a lower error band;
判断所述误差带上界是否大于或者等于所述上峰值集合的最大值,且所述误差带下界是否大于或者等于所述下峰值集合的最小值,若判断结果均为是,则根据采样间隔确定调节时间;Judging whether the upper bound of the error band is greater than or equal to the maximum value of the upper peak set, and whether the lower bound of the error band is greater than or equal to the minimum value of the lower peak set, if the judgment results are all yes, according to the sampling interval determine the adjustment time;
判断所述风电场调频响应特性实测数据是否为频率阶跃下扰,若是,则根据所述上峰值集合的最大值确定超调量;Determine whether the measured data of the frequency modulation response characteristic of the wind farm is a frequency step underdisturbance, and if so, determine the overshoot amount according to the maximum value of the upper peak set;
根据所述数据集确定上升时间或下降时间。Rise time or fall time is determined from the data set.
本实施例中的第一步是对风电场调频响应特性实测数据进行分析和预处理,在实际中,描述一个系统动态性能的指标一般为上升(下降)时间tr、调节时间ts以及超调量δ%表示。其中tr为响应值从零到达终值p∞所用的时间,ts为响应值从零到其稳定在人为规定的误差带内所需时间,δ%由以下公式求得The first step in this embodiment is to analyze and preprocess the measured data of the frequency modulation response characteristics of the wind farm. In practice, the indicators describing the dynamic performance of a system are generally the rise (fall) time tr , the adjustment time ts , and the overrun time. The adjustment amount is expressed as δ%. where t r is the time it takes for the response value to reach the final value p ∞ from zero, t s is the time it takes for the response value to change from zero to the time it is stabilized within the artificially specified error band, and δ% is obtained from the following formula
其中pex为调节时间内响应值的最大偏移量。where p ex is the maximum offset of the response value within the adjustment time.
由于风电场频率阶跃响应无法用简单的一阶或二阶的数学模型描述,故无法准确的用解析表达式求出响应的动态性能指标,为了解决这个问题,现提出基于实测数据的阶跃响应动态性能指标求解算法,为方便理解,该算法以伪代码形式体现,算法的步骤如下:Since the frequency step response of the wind farm cannot be described by a simple first-order or second-order mathematical model, the dynamic performance index of the response cannot be accurately obtained by analytical expressions. In order to solve this problem, a step based on the measured data is proposed The response dynamic performance index solution algorithm, for the convenience of understanding, the algorithm is embodied in the form of pseudo code, the steps of the algorithm are as follows:
选取指标中调节时间内的数据来进行建模,并对数据进行如下处理:以风电场频率干扰发生时刻t0为初始时刻,选取t0时刻对应的有功功率p0,则风电场调频过程中各时刻对应的有功功率变化量为The data in the adjustment time in the index is selected for modeling, and the data is processed as follows: take the time t 0 of the frequency disturbance of the wind farm as the initial time, and select the active power p 0 corresponding to the time t 0 , then the wind farm frequency adjustment process The active power change corresponding to each moment is
Δpi=pi-p0 (2) Δpi = pi -p 0 (2)
其中i=0,1,…,n,n为发生完整频率阶跃事件的数据个数。pi为ti时刻对应的实测数据中风电场的有功功率值。以50Hz为并网标准频率则频率阶跃过程中各时刻对应的频率偏差量为where i=0,1,...,n,n is the number of data with complete frequency step events. p i is the active power value of the wind farm in the measured data corresponding to time t i . Taking 50Hz as the grid-connected standard frequency, the frequency deviation corresponding to each moment in the frequency step process is:
Δfi=fi-50 (3)Δf i =f i -50 (3)
其中fi为ti时刻对应的实测数据中风电场的频率值。为方便计算,定义Δf=(Δf0,Δf1,...,Δfn)T,以Δp=(Δp0,Δp1,...,Δpn)T。where f i is the frequency value of the wind farm in the measured data corresponding to time t i . For the convenience of calculation, Δf=(Δf 0 ,Δf 1 ,...,Δf n ) T is defined, and Δp=(Δp 0 ,Δp 1 ,...,Δp n ) T .
进一步地,以频率上扰为例,以风电场调频时的初始有功功率划分不同的工况,根据风电场实测数据,具体分为,工况1:限功率25%下的频率阶跃上扰;工况2:限功率50%下的频率阶跃上扰;工况3:不限功率下的频率阶跃上扰;工况4:复合频率上扰;根据本发明提出的算法可求出不同工况下风电场调频响应特性指标,如表1所示,表1为风电场在各工况下频率阶跃上扰响应特性指标,根据表1中的下降时间,截取建模所需数据,值得注意的是,工况4,是复合频率干扰,不是典型的阶跃响应,故求取性能指标没有实际的价值,故表1没有给出工况4的指标,并且针对工况4,建模数据选取完整的调频响应时间内的数据。选定好建模所需数据之后,根据方案第一步中预处理方法对数据进行预处理。Further, taking the frequency up-disturbance as an example, the initial active power of the wind farm during frequency modulation is divided into different working conditions, and according to the measured data of the wind farm, it is specifically divided into, working condition 1: frequency step up-disturbance at 25% of the limited power. ; working condition 2: frequency step up disturbance under 50% power limit; working condition 3: frequency step up disturbance under unlimited power; working condition 4: compound frequency up disturbance; The frequency modulation response characteristic indexes of the wind farm under different working conditions are shown in Table 1. Table 1 is the frequency step up disturbance response characteristic index of the wind farm under various working conditions. According to the fall time in Table 1, the data required for modeling is intercepted , It is worth noting that working condition 4 is a composite frequency interference, not a typical step response, so the performance index has no practical value, so Table 1 does not give the index of working condition 4, and for working condition 4, Modeling data selects data within the full FM response time. After selecting the data required for modeling, preprocess the data according to the preprocessing method in the first step of the plan.
表1Table 1
优选地,所述步骤200包括:Preferably, the
分别根据各个工况下的所述处理后的数据中的频率变化量为输入,功率变化量为输出,构建初始函数模型;According to the frequency variation in the processed data under each working condition as the input, and the power variation as the output, an initial function model is constructed;
设定传递函数的分子阶数和分母阶数;Set the numerator order and denominator order of the transfer function;
对所述初始传递函数模型进行模型辨识,并调整所述分子阶数和所述分母阶数,得到最优的所述传递函数模型。Perform model identification on the initial transfer function model, and adjust the numerator order and the denominator order to obtain the optimal transfer function model.
优选地,所述步骤200还包括:Preferably, the
根据不同的所述传递函数模型的正交投影矩阵确定间隙度量公式;Determine the gap measurement formula according to the orthogonal projection matrices of the different transfer function models;
根据所述间隙度量公式计算所述间隙值;所述间隙值用于确定所述工况区域。The gap value is calculated according to the gap metric formula; the gap value is used to determine the operating condition area.
具体的,本实施例中第二步为根据各工况下处理好的数据对每个工况建立传递函数模型并利用间隙测度方法测量各模型间的间隙值,以确定建模工况区域,具体方法如下:Specifically, the second step in this embodiment is to establish a transfer function model for each working condition according to the processed data under each working condition, and use the gap measurement method to measure the gap value between the models, so as to determine the modeled working condition area, The specific method is as follows:
以频率变化量为输入,相应的功率变化量为输出,给出以下传递函数模型Taking the frequency change as the input and the corresponding power change as the output, the following transfer function model is given
其中m和k代表传递函数分子和分母的阶数,并假设k≥m,j=1,2…n,代表第j次频率波动事件,θj=[b1,j,…,bk,j,a0,j,…,am,j]T是离散传递函数参数组成的向量。where m and k represent the order of the numerator and denominator of the transfer function, and assuming k≥m, j=1, 2...n, represents the jth frequency fluctuation event, θ j =[b 1,j ,...,b k, j ,a 0,j ,…,am ,j ] T is a vector of discrete transfer function parameters.
接着(4)式可化为:Then (4) can be transformed into:
ΔPj(z)(1+b1,jz-1+…+bk,jz-k)=Δfj(z)(a0,j+a1,jz-1+…+am,jz-m) (5)ΔP j (z)(1+b 1,j z -1 +...+b k,j z -k )=Δf j (z)(a 0,j +a 1,j z -1 +...+am ,j z -m ) (5)
根据z-变换的性质,(5)可以直接转换为差分方程如下:According to the properties of the z-transform, (5) can be directly transformed into the difference equation as follows:
ΔPj(t)=-b1,jΔPj(t-1)-…-bk,jΔPj(t-k)+a0,jΔfj(t)+…+am,jΔfj(t-m) (6)ΔP j (t)=-b 1,j ΔP j (t-1)-…-b k,j ΔP j (tk)+a 0,j Δf j (t)+…+am ,j Δf j ( tm) (6)
假设第j次频率波动事件有N个实测样本点,并令Assume that the jth frequency fluctuation event has N measured sample points, and let
Yj=[ΔPj(t),…,ΔPj(N)]T,Y j =[ΔP j (t),...,ΔP j (N)] T ,
则(6)式可由上述矩阵表达为如下紧凑形式:Then equation (6) can be expressed by the above matrix in the following compact form:
Yj=Ξjθj (7)Y j = Ξ j θ j (7)
由此我们可以利用线性最小二乘法来求解θj,即From this, we can use the linear least squares method to solve θ j , namely
在建模过程中,可以假设传递函数分子和分母的阶数m和k已知,由传递函数的性质可知:m和k的值可以反应出输出向量曲线的极值和零点个数。本发明在风电场调频响应特性分析的过程中,首先根据调频事件过程中有功功率变化量的曲线,初步判断传递函数的极点和零点个数,给出m和k的值,然后进行模型参数的辨识,最后通过不断调整m和k的值,来得到最优模型。In the modeling process, it can be assumed that the order m and k of the numerator and denominator of the transfer function are known. From the properties of the transfer function, it can be known that the values of m and k can reflect the extreme value and the number of zero points of the output vector curve. In the process of analyzing the frequency regulation response characteristics of the wind farm, the present invention firstly judges the number of poles and zeros of the transfer function according to the curve of the active power change in the frequency regulation event process, and gives the values of m and k, and then calculates the model parameters. Identify, and finally obtain the optimal model by continuously adjusting the values of m and k.
在得到的传递传递函数模型之后,利用间隙度量的方法比较每个模型之间的动态特性差别,来达到用一个模型表征在一定工况范围内的调频响应的目的。若两个系统的传递函数分别为Q1和Q2则这两个系统间的间隙可以表示为:After the transfer function model is obtained, the gap measurement method is used to compare the dynamic characteristics difference between each model, so as to achieve the purpose of using a model to characterize the frequency modulation response within a certain range of operating conditions. If the transfer functions of the two systems are Q1 and Q2 , respectively, then the gap between the two systems can be expressed as:
式中,G(Qi)(i=1,2)表示Qi对应的系统的图,正交投影矩阵定义为:In the formula, G(Q i ) (i=1,2) represents the graph of the system corresponding to Q i , and the orthogonal projection matrix defined as:
其中,i=1,2,Ni,Di∈RH∞可由Qi规范右互质分解得到,即Qi=NiDi -1,并满足(I为单位矩阵,*表示共轭),由(9)、(10)可得两个系统的间隙度量公式为:Among them, i=1,2,N i ,D i ∈RH ∞ can be obtained from the normative right coprime decomposition of Qi i , that is, Qi i =N i D i -1 , and satisfy (I is the identity matrix, * means conjugation), the gap metric formula of the two systems can be obtained from (9) and (10):
其中,P为任意的希尔伯特矩阵,0≤g≤1,代表两个系统动态特性差别越大,g值越靠近0,则两个系统动态特性差别越小。g值越靠近1,动态特性差别越大。Among them, P is an arbitrary Hilbert matrix, 0≤g≤1, which means that the greater the difference in the dynamic characteristics of the two systems, the closer the g value is to 0, the smaller the difference in the dynamic characteristics of the two systems. The closer the g value is to 1, the greater the difference in dynamic characteristics.
进一步地,利用本实施例中的第二步,得到不同工况下的传递函数模型,其拟合图如图2至图5所示,由图2至图5可知,传递函数拟合度较好,然后根据间隙测度方法,得到各模型间的间隙值均接近于1,因此我们将所有工况的数据进行合并。与此同时,间隙值也说明了传递函数模型需要针对不同工况的模型分别建模,当工况增加时,工作量也会变得巨大。根据图5可知,传递函数模型对复杂的频率响应事件表征能力不强。Further, using the second step in this embodiment, the transfer function models under different working conditions are obtained, and the fitting diagrams thereof are shown in Figs. Well, then according to the gap measurement method, the gap values between the models are all close to 1, so we merge the data of all working conditions. At the same time, the gap value also shows that the transfer function model needs to be modeled separately for different working conditions. When the working conditions increase, the workload will also become huge. According to Fig. 5, it can be seen that the transfer function model has poor ability to characterize complex frequency response events.
优选地,所述步骤300包括:Preferably, the
计算不同工况的所述传递函数模型之间的间隙值与0之间的距离,得到第一距离;Calculate the distance between the gap value and 0 between the transfer function models of different working conditions to obtain the first distance;
计算不同工况的所述传递函数模型之间的间隙值与1之间的距离,得到第二距离;Calculate the distance between the gap value and 1 between the transfer function models of different working conditions to obtain the second distance;
判断所述第一距离是否小于或者等于所述第二距离,若是,则按时序合并这两个工况的调频数据,得到合并后的数据;若否,则选取任意一个工况的调频数据与已有的调频数据进行时序合并,得到合并后的数据。Determine whether the first distance is less than or equal to the second distance, and if so, combine the FM data of the two working conditions according to the time series to obtain the combined data; if not, select the FM data of any working condition and The existing FM data is time-series merged to obtain merged data.
优选地,所述步骤300还包括:Preferably, the
构建初始神经网络;Build the initial neural network;
基于列文伯格-马尔夸克算法,根据所述合并后的数据对所述初始神经网络进行训练;Based on the Levenberg-Marquark algorithm, the initial neural network is trained according to the combined data;
利用每一个工况的调频数据对训练后的神经网络进行测试,并根据测试结果确定所述非线性自回归神经网络模型。The trained neural network is tested by using the frequency modulation data of each working condition, and the nonlinear autoregressive neural network model is determined according to the test results.
具体的,本实施例中的第三步是根据各工况下的传递函数模型的间隙值,合并相应数据,进行非线性自回归神经网络模型的构建,具体方法如下:Specifically, the third step in this embodiment is to combine the corresponding data according to the gap value of the transfer function model under each working condition to construct the nonlinear autoregressive neural network model. The specific method is as follows:
由第二步得到各模型间的间隙值之后,若两个工况传递函数模型之间的间隙值接近1,则说明其动态特性差别很大,则按时序合并这两个工况的调频数据,得到数据集Du,若间隙值接近于0,则说明两个模型有很相似的动态特性,则选取任意一个工况的调频数据与已有Du进行时序合并,得到新的Du。利用合并完成的数据集Du对神经网络进行训练,下面介绍非线性自回归神经网络:After the gap value between the models is obtained from the second step, if the gap value between the transfer function models of the two working conditions is close to 1, it means that the dynamic characteristics are very different, and the FM data of the two working conditions are merged according to the time series. , to obtain the data set Du , if the gap value is close to 0, it means that the two models have very similar dynamic characteristics, then select the frequency modulation data of any working condition and merge the existing Du to obtain a new Du . Use the merged dataset Du to train the neural network. The following describes the nonlinear autoregressive neural network:
非线性自回归神经网络的输出不仅依赖于当前的输入,也依赖于过去的输入和输出,这使得其具有一定的记忆功能。风电场有调频过程是动态的且调频数据也是时序序列,故风电场调频过程的非线性自回归神经网络模型如下:The output of the nonlinear autoregressive neural network not only depends on the current input, but also on the past input and output, which makes it have a certain memory function. The frequency modulation process of the wind farm is dynamic and the frequency modulation data is also a time series, so the nonlinear autoregressive neural network model of the frequency modulation process of the wind farm is as follows:
ΔP(t)=h(ΔP(t-1),ΔP(t-1),…,ΔP(t-d),Δf(t),Δf(t-1),…,Δf(t-d)) (12)ΔP(t)=h(ΔP(t-1), ΔP(t-1), ..., ΔP(t-d), Δf(t), Δf(t-1), ..., Δf(t-d)) (12)
其中,ΔP(t)为t时刻对应的有功功率变化量即为网络的输出,Δf(t)为t时刻对应的频率变化量即为网络的输入,h(·)为网络构造的关于ΔP(t)和Δf(t)的非线性函数。d为延时数。其网络结构如图6所示。Among them, ΔP(t) is the active power change corresponding to time t, which is the output of the network, Δf(t) is the frequency change corresponding to time t, which is the input of the network, and h( ) is the network construction about ΔP( t) and non-linear functions of Δf(t). d is the delay number. Its network structure is shown in Figure 6.
本发明采用列文伯格-马尔夸克算法来训练非线性自回归神经网络,该算法继承了牛顿法和梯度下降法的优点,并且收敛速度快且稳定。其参数更新规则为:The present invention adopts the Levenberg-Marquark algorithm to train the nonlinear autoregressive neural network, the algorithm inherits the advantages of the Newton method and the gradient descent method, and the convergence speed is fast and stable. Its parameter update rules are:
式中,qn为第n次迭代的网络权重,e(·)为误差向量。H和J分别为雅可比矩阵和海森矩阵,该算法利用雅可比矩阵给出海森矩阵的近似:In the formula, q n is the network weight of the nth iteration, and e( ) is the error vector. H and J are the Jacobian and Hessian matrices, respectively. The algorithm uses the Jacobian to give an approximation of the Hessian matrix:
H≈ηI+JTJ (14)H≈ηI+J T J (14)
其中,η和I分别为学习系数和单位矩阵。where n and I are the learning coefficient and identity matrix, respectively.
列文伯格-马尔夸克算法训练过程基于调整η值,通过迭代寻找最小的误差,当误差值小于上次迭代误差则减小η值,反之,则增大η值。η接近于0时该算法具有近似于梯度下降法的局部收敛性,η很大时该算法具有近似于高斯牛顿法的全局收敛性。The training process of the Levenberg-Marquark algorithm is based on adjusting the value of η, and iteratively finds the smallest error. When the error value is smaller than the error of the previous iteration, the value of η is decreased, otherwise, the value of η is increased. When η is close to 0, the algorithm has the local convergence approximate to the gradient descent method, and when η is large, the algorithm has the global convergence approximate to the Gauss-Newton method.
进一步地,利用方案第三步,将相应工况的数据合并,对非线性自回归神经网络进行训练,并利用不用工况的数据去测试训练好的网络模型,得到图7至图10,由图7至图10可知相较于传递函数模型在非线性自回归神经网络模型可以同时表征多种工况,同时对非线性更强的调频响应同样有良好的表征能力。Further, using the third step of the scheme, the data of the corresponding working conditions are merged, the nonlinear autoregressive neural network is trained, and the trained network model is tested with the data of the different working conditions, and Figures 7 to 10 are obtained. Figures 7 to 10 show that, compared with the transfer function model, the nonlinear autoregressive neural network model can simultaneously represent multiple operating conditions, and also has a good ability to characterize the more nonlinear FM response.
本发明还提供了一种数据驱动的风电场调频响应特性建模系统,包括:The present invention also provides a data-driven wind farm frequency modulation response characteristic modeling system, comprising:
数据处理模块,用于基于阶跃响应动态性能指标求解算法对风电场调频响应特性实测数据进行分析和预处理,得到各个工况下处理后的数据;The data processing module is used to analyze and preprocess the measured data of the frequency modulation response characteristics of the wind farm based on the step response dynamic performance index solution algorithm, and obtain the processed data under each working condition;
传递函数构建模块,用于根据所述处理后的数据对每个工况建立传递函数模型,并利用间隙测度方法测量各模型间的间隙值,以确定非线性自回归神经网络模型表征的工况区域;The transfer function building module is used to establish a transfer function model for each working condition according to the processed data, and use the gap measurement method to measure the gap value between the models to determine the working condition represented by the nonlinear autoregressive neural network model area;
神经网络建模模块,用于根据所述间隙值对工况的调频数据进行合并,并根据合并后的数据对所述非线性自回归神经网络模型进行训练,得到训练好的非线性自回归神经网络模型。A neural network modeling module is used to merge the frequency modulation data of the working condition according to the gap value, and train the nonlinear autoregressive neural network model according to the merged data to obtain a trained nonlinear autoregressive neural network. network model.
优选地,所述数据处理模块具体包括:Preferably, the data processing module specifically includes:
工况划分单元,用于根据风电场调频时的初始有功功率对所述风电场调频响应特性实测数据进行划分,得到多个工况;A working condition dividing unit, configured to divide the measured data of the frequency modulation response characteristic of the wind farm according to the initial active power during frequency modulation of the wind farm to obtain a plurality of working conditions;
峰值计算单元,用于获取所述风电场调频响应特性实测数据的数据集的功率曲线,并计算所述功率曲线的各个数据点的上峰值集合和下峰值集合;a peak value calculation unit, configured to obtain the power curve of the data set of the measured data of the frequency modulation response characteristic of the wind farm, and calculate the upper peak value set and the lower peak value set of each data point of the power curve;
最值计算单元,用于计算所述上峰值集合的最大值和所述下峰值集合的最小值;a maximum value calculation unit, configured to calculate the maximum value of the upper peak set and the minimum value of the lower peak set;
上下界计算单元,用于基于所述数据集中的数据点,计算误差带上界和误差带下界;an upper and lower bound calculation unit, configured to calculate the upper bound of the error band and the lower bound of the error band based on the data points in the data set;
第一判断单元,用于判断所述误差带上界是否大于或者等于所述上峰值集合的最大值,且所述误差带下界是否大于或者等于所述下峰值集合的最小值,若判断结果均为是,则根据采样间隔确定调节时间;A first judging unit, configured to judge whether the upper bound of the error band is greater than or equal to the maximum value of the upper peak set, and whether the lower bound of the error band is greater than or equal to the minimum value of the lower peak set, if the judgment results are all If yes, determine the adjustment time according to the sampling interval;
第二判断单元,用于判断所述风电场调频响应特性实测数据是否为频率阶跃下扰,若是,则根据所述上峰值集合的最大值确定超调量;a second judging unit, configured to judge whether the measured data of the frequency modulation response characteristic of the wind farm is a frequency step underdisturbance, and if so, determine the overshoot amount according to the maximum value of the upper peak set;
时间确定单元,用于根据所述数据集确定上升时间或下降时间。A time determination unit, configured to determine a rise time or a fall time according to the data set.
优选地,所述传递函数构建模块具体包括:Preferably, the transfer function building module specifically includes:
初始模型构建单元,用于分别根据各个工况下的所述处理后的数据中的频率变化量为输入,功率变化量为输出,构建初始函数模型;an initial model building unit, configured to construct an initial function model according to the frequency variation in the processed data under each working condition as the input and the power variation as the output;
阶数设定单元,用于设定传递函数的分子阶数和分母阶数;Order setting unit, used to set the numerator order and denominator order of the transfer function;
模型确定单元,用于对所述初始传递函数模型进行模型辨识,并调整所述分子阶数和所述分母阶数,得到最优的所述传递函数模型。A model determining unit, configured to perform model identification on the initial transfer function model, and adjust the numerator order and the denominator order to obtain the optimal transfer function model.
优选地,所述传递函数构建模块具体包括:Preferably, the transfer function building module specifically includes:
公式确定单元,用于根据不同的所述传递函数模型的正交投影矩阵确定间隙度量公式;a formula determination unit, used for determining a gap measurement formula according to the orthogonal projection matrices of the different transfer function models;
间隙计算单元,用于根据所述间隙度量公式计算所述间隙值;所述间隙值用于确定所述工况区域。A gap calculation unit, configured to calculate the gap value according to the gap metric formula; the gap value is used to determine the working condition area.
本发明的有益效果如下:The beneficial effects of the present invention are as follows:
(1)发明提出了一种基于数据的频率阶跃响应特性指标求解算法,用来分析风电场调频响应特性以及处理建模所需数据。(1) The invention proposes a data-based frequency step response characteristic index solution algorithm, which is used to analyze the frequency modulation response characteristic of a wind farm and process the data required for modeling.
(2)本发明构建了不同工况下的风电场调频响应特性的传递函数模型,并利用一种模型间隙测度方法来衡量模型间的动态相似度,来确定工况区域。(2) The present invention constructs the transfer function model of the frequency modulation response characteristic of the wind farm under different working conditions, and uses a model gap measurement method to measure the dynamic similarity between the models to determine the working condition area.
(3)本发明构建了NARX神经网络模型来表征多工况下的风电场调频响应特性。此模型克服了传递函数模型表征工况范围能力和对非线性强的响应表征能力交差的缺点,经过仿真验证,此模型在风电场调频响应特性建模领域有着良好的应用前景。(3) The present invention constructs a NARX neural network model to characterize the frequency modulation response characteristics of the wind farm under multiple operating conditions. This model overcomes the disadvantage of the cross between the transfer function model's ability to characterize the range of operating conditions and the ability to characterize strong nonlinear responses. After simulation verification, this model has a good application prospect in the field of wind farm frequency modulation response characteristics modeling.
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。The various embodiments in this specification are described in a progressive manner, and each embodiment focuses on the differences from other embodiments, and the same and similar parts between the various embodiments can be referred to each other. For the system disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant part can be referred to the description of the method.
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。In this paper, specific examples are used to illustrate the principles and implementations of the present invention. The descriptions of the above embodiments are only used to help understand the methods and core ideas of the present invention; meanwhile, for those skilled in the art, according to the present invention There will be changes in the specific implementation and application scope. In conclusion, the contents of this specification should not be construed as limiting the present invention.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210229805.8A CN114583767B (en) | 2022-03-10 | 2022-03-10 | Data-driven wind power plant frequency modulation response characteristic modeling method and system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210229805.8A CN114583767B (en) | 2022-03-10 | 2022-03-10 | Data-driven wind power plant frequency modulation response characteristic modeling method and system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114583767A true CN114583767A (en) | 2022-06-03 |
CN114583767B CN114583767B (en) | 2023-03-17 |
Family
ID=81773321
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210229805.8A Active CN114583767B (en) | 2022-03-10 | 2022-03-10 | Data-driven wind power plant frequency modulation response characteristic modeling method and system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114583767B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115498660A (en) * | 2022-11-17 | 2022-12-20 | 华北电力大学 | Wind farm frequency modulation dynamic modeling method, device and electronic equipment |
CN117665824A (en) * | 2023-12-22 | 2024-03-08 | 中山大学 | A sea surface wind field reconstruction method and system |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109004687A (en) * | 2018-08-03 | 2018-12-14 | 山东大学 | The intelligent inertia response control mehtod and system of wind power plant participation power grid frequency modulation |
CN110426953A (en) * | 2019-07-18 | 2019-11-08 | 国网山东省电力公司电力科学研究院 | AGC method of evaluating performance based on fired power generating unit generation model |
CN110556842A (en) * | 2019-09-16 | 2019-12-10 | 湖南大学 | Control method of direct-drive wind power plant inductive weak grid-connected subsynchronous oscillation suppression device |
CN110765703A (en) * | 2019-11-04 | 2020-02-07 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | Wind power plant aggregation characteristic modeling method |
CN111525593A (en) * | 2020-04-10 | 2020-08-11 | 中国电力科学研究院有限公司 | A method and system for fitting primary frequency modulation parameters of a wind farm |
AU2020102245A4 (en) * | 2019-01-08 | 2020-10-29 | Nanjing Institute Of Technology | A grid hybrid rolling dispatching method considering congestion and energy storage tou price |
US20210079888A1 (en) * | 2016-12-26 | 2021-03-18 | Beijing Goldwind Science & Creation Windpower Equipment Co., Ltd. | Test and control apparatus, system and method for wind farm |
CN112560352A (en) * | 2020-12-24 | 2021-03-26 | 华北电力大学 | System frequency response model modeling method based on AM-LSTM neural network |
CN112636366A (en) * | 2020-12-01 | 2021-04-09 | 国家电网有限公司 | Wind power plant dynamic frequency control method based on control process data fitting |
CN113708389A (en) * | 2021-09-10 | 2021-11-26 | 国网湖南省电力有限公司 | Wind power plant primary frequency modulation model parameter identification method and system based on actual power response |
CN113837432A (en) * | 2021-08-12 | 2021-12-24 | 华北电力大学 | Power system frequency prediction method driven by physics-data combination |
CN113964884A (en) * | 2021-11-17 | 2022-01-21 | 国家电网有限公司华东分部 | Power grid active frequency regulation and control method based on deep reinforcement learning |
CN113988481A (en) * | 2021-12-23 | 2022-01-28 | 南京鼐威欣信息技术有限公司 | Wind power prediction method based on dynamic matrix prediction control |
WO2022021470A1 (en) * | 2020-07-28 | 2022-02-03 | 南京东博智慧能源研究院有限公司 | Prediction error distribution estimation method for frequency modulation potential of wind turbines |
CN114065460A (en) * | 2020-08-03 | 2022-02-18 | 北京国电智深控制技术有限公司 | Model processing method, storage medium and electronic device in thermal power generation system |
-
2022
- 2022-03-10 CN CN202210229805.8A patent/CN114583767B/en active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20210079888A1 (en) * | 2016-12-26 | 2021-03-18 | Beijing Goldwind Science & Creation Windpower Equipment Co., Ltd. | Test and control apparatus, system and method for wind farm |
CN109004687A (en) * | 2018-08-03 | 2018-12-14 | 山东大学 | The intelligent inertia response control mehtod and system of wind power plant participation power grid frequency modulation |
AU2020102245A4 (en) * | 2019-01-08 | 2020-10-29 | Nanjing Institute Of Technology | A grid hybrid rolling dispatching method considering congestion and energy storage tou price |
CN110426953A (en) * | 2019-07-18 | 2019-11-08 | 国网山东省电力公司电力科学研究院 | AGC method of evaluating performance based on fired power generating unit generation model |
CN110556842A (en) * | 2019-09-16 | 2019-12-10 | 湖南大学 | Control method of direct-drive wind power plant inductive weak grid-connected subsynchronous oscillation suppression device |
CN110765703A (en) * | 2019-11-04 | 2020-02-07 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | Wind power plant aggregation characteristic modeling method |
CN111525593A (en) * | 2020-04-10 | 2020-08-11 | 中国电力科学研究院有限公司 | A method and system for fitting primary frequency modulation parameters of a wind farm |
WO2022021470A1 (en) * | 2020-07-28 | 2022-02-03 | 南京东博智慧能源研究院有限公司 | Prediction error distribution estimation method for frequency modulation potential of wind turbines |
CN114065460A (en) * | 2020-08-03 | 2022-02-18 | 北京国电智深控制技术有限公司 | Model processing method, storage medium and electronic device in thermal power generation system |
CN112636366A (en) * | 2020-12-01 | 2021-04-09 | 国家电网有限公司 | Wind power plant dynamic frequency control method based on control process data fitting |
CN112560352A (en) * | 2020-12-24 | 2021-03-26 | 华北电力大学 | System frequency response model modeling method based on AM-LSTM neural network |
CN113837432A (en) * | 2021-08-12 | 2021-12-24 | 华北电力大学 | Power system frequency prediction method driven by physics-data combination |
CN113708389A (en) * | 2021-09-10 | 2021-11-26 | 国网湖南省电力有限公司 | Wind power plant primary frequency modulation model parameter identification method and system based on actual power response |
CN113964884A (en) * | 2021-11-17 | 2022-01-21 | 国家电网有限公司华东分部 | Power grid active frequency regulation and control method based on deep reinforcement learning |
CN113988481A (en) * | 2021-12-23 | 2022-01-28 | 南京鼐威欣信息技术有限公司 | Wind power prediction method based on dynamic matrix prediction control |
Non-Patent Citations (3)
Title |
---|
姚琦等: "含异步变速风机的风电场一次调频等值建模与仿真", 《电力系统自动化》 * |
宁夏电力调度控制中心等: "《新能源消纳调控技术》", 30 June 2021 * |
栾贻贺: "适用于频率控制的风电场建模与调频策略研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115498660A (en) * | 2022-11-17 | 2022-12-20 | 华北电力大学 | Wind farm frequency modulation dynamic modeling method, device and electronic equipment |
US11847530B1 (en) | 2022-11-17 | 2023-12-19 | North China Electric Power University | Frequency modulation dynamic modeling method and device for wind farm, and electronic device |
CN117665824A (en) * | 2023-12-22 | 2024-03-08 | 中山大学 | A sea surface wind field reconstruction method and system |
CN117665824B (en) * | 2023-12-22 | 2024-08-20 | 中山大学 | Sea surface wind field reconstruction method and system |
Also Published As
Publication number | Publication date |
---|---|
CN114583767B (en) | 2023-03-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106443285B (en) | Quantitative Analysis Method of Harmonic Responsibility for Multi-harmonic Sources Based on Total Least Squares | |
Chen et al. | Data-driven piecewise linearization for distribution three-phase stochastic power flow | |
CN106548021B (en) | A kind of electronic mutual inductor error state appraisal procedure based on higher-dimension random matrix | |
CN106408223A (en) | Short-term load prediction based on meteorological similar day and error correction | |
CN115455793B (en) | Mechanical analysis method of complex components of high-rise structures based on multi-scale model modification | |
CN110442911B (en) | An uncertainty analysis method for high-dimensional complex systems based on statistical machine learning | |
CN114583767B (en) | Data-driven wind power plant frequency modulation response characteristic modeling method and system | |
CN112485747A (en) | Error state evaluation method for capacitor voltage transformer | |
CN118311464B (en) | A battery health status estimation method and system | |
CN113139605A (en) | Power load prediction method based on principal component analysis and LSTM neural network | |
CN109060393B (en) | Constant-load response time domain fusion analysis method for bridge structure | |
CN114935688B (en) | Harmonic evaluation method and system for electric arc furnace power supply system based on power segmentation | |
CN102779208A (en) | Sequential accelerated degradation test optimal design method based on relative entropy | |
CN111175608A (en) | Power distribution network harmonic responsibility quantitative division method based on accelerated independent component analysis | |
CN114372344A (en) | A Quantitative Identification Method for Influencing Factors of Subsynchronous Oscillation Damping Characteristics | |
CN105406461A (en) | Adaptive dynamic load monitoring method for power distribution network power failure events | |
CN107025609A (en) | Electrical Power System Dynamic method for estimating state based on singular value decomposition CDKF | |
CN105808962A (en) | Assessment method considering voltage probabilities of multiple electric power systems with wind power output randomness | |
CN114626207B (en) | A Method of Constructing a General Probabilistic Model for Harmonic Emission Levels of Industrial Loads | |
CN118627925A (en) | A method for evaluating causes of short-term wind power forecasting errors based on principal component analysis | |
CN110309873B (en) | Integrated transient stability evaluation method and system based on Bayesian classification of power flow vectors | |
CN117236744A (en) | A harmonic impact assessment method considering multiple harmonic sources | |
JP5989571B2 (en) | Fluctuation estimation method, fluctuation estimation apparatus, and fluctuation estimation program for total power output of natural energy type distributed power supply group | |
CN116681158A (en) | Reference crop evapotranspiration prediction method based on integrated extreme learning machine | |
CN113125037B (en) | Cable conductor temperature estimation method based on distributed optical fiber on-line temperature measurement system |
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 |