CN112636366B - Wind power plant dynamic frequency control method based on control process data fitting - Google Patents
Wind power plant dynamic frequency control method based on control process data fitting Download PDFInfo
- Publication number
- CN112636366B CN112636366B CN202011398967.1A CN202011398967A CN112636366B CN 112636366 B CN112636366 B CN 112636366B CN 202011398967 A CN202011398967 A CN 202011398967A CN 112636366 B CN112636366 B CN 112636366B
- Authority
- CN
- China
- Prior art keywords
- fan
- moment
- matrix
- dynamic
- control
- 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
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/24—Arrangements for preventing or reducing oscillations of power in networks
- H02J3/241—The oscillation concerning frequency
-
- 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/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
-
- 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/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- 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)
- Power Engineering (AREA)
- Feedback Control In General (AREA)
Abstract
The invention provides a data-driven dynamic optimization control method for the frequency stability control of a high-proportion wind power independent power system, so that the wind power plant is stabilized and optimized in the internal transient process while participating in system frequency modulation. According to the wind field frequency optimization control method, an initial data set is constructed, on-line dynamic modeling and centralized wind field control modeling are carried out, and wind field frequency optimization control is carried out according to a fan on-line dynamic modeling result. According to the invention, the online equivalent dynamic model of the fan is obtained by fitting the historical data generated in the operation process of the fan, and the online equivalent dynamic model has the characteristic of pure data driving, so that the online equivalent dynamic model can be suitable for operation control of fans of different manufacturers under different working conditions, meanwhile, the algorithm has a pure linear form, the solving process is simple, the calculation load is small, the complex nonlinear optimization control problem is not required to be solved, the real-time operation requirement is met, and the balance of dynamic response precision and speed can be considered.
Description
Technical Field
The invention belongs to the technical field of electric power, and particularly relates to a wind power plant dynamic frequency control method based on control process data fitting.
Background
In recent decades, with the increasing energy crisis and environmental problems caused by the rapid consumption of global fossil energy, renewable energy sources such as wind power and photovoltaics keep a rapid increase trend in the installed capacity around the world. Taking wind power as an example, by the end of 2019, the total installed amount of wind power worldwide reaches 594GW, and the average annual growth from 2010 is 46GW. Because the variable speed fan can adapt to a wider wind speed working range, most of newly-added fans adopt variable speed fans. However, while the wind energy utilization rate is improved, the active output and the frequency of the variable speed fan are decoupled, and as the permeability of the power grid of the variable speed fan is continuously improved, the power grid has the trend and the characteristics of power electronization, so that the problem of the frequency stability of the whole power grid becomes a great challenge for large-scale access of the fan to the power grid.
The existing fan participating in the frequency modulation technical route comprises two types of local control and coordinated control, wherein the local control introduces frequency response characteristics by modifying a local control loop of the fan, but because the in-field fan lacks coordinated interaction, the overall response is difficult to ensure optimality under the condition of uneven in-field wind speed distribution, and meanwhile, the control parameters of a plurality of devices in the field are required to be set one by one at regular intervals, so that the flexibility of the overall response is greatly reduced. And the coordination control is performed by relying on a dynamic time domain model of the fan to perform the coordination optimization control of the whole wind field. At present, a dynamic modeling method for a fan mainly depends on experience cognition and parameter measurement of a physical model of the fan, on one hand, the method is difficult to fully consider actual processes such as each time lag of a complete control closed loop in the process of describing the physical model, and on the other hand, the accuracy and even the feasibility of parameter measurement are often limited by actual working conditions, so that the actual approximation effect of a final dynamic model is difficult to ensure, and the final optimal control result is greatly influenced.
Therefore, a data-driven wind field frequency optimization control method needs to be researched, and the method is suitable for the actual requirements of wind fields for rapidly and flexibly participating in dynamic response of the power grid frequency.
Disclosure of Invention
Aiming at the problems, the wind farm dynamic frequency control method based on the control process data fitting is suitable for the actual requirements of the wind farm for rapidly and flexibly participating in the power grid frequency dynamic response.
A wind farm dynamic frequency control method based on control process data fitting, the control method comprising:
on-line dynamic modeling of the fan;
and optimally controlling the wind field frequency according to the online dynamic modeling result of the fan.
Further, on-line dynamic modeling of fans is performed by:
constructing an initial data set;
on-line dynamic modeling;
and (5) modeling the concentrated wind field control.
Further, the initial dataset is constructed by:
establishing a state equation of a single fan or a power generation unit:
ω k+1 =f(ω k ,u k ) (1.1)
f is a state transition relation function, w k The fan rotation speed at the moment k is u k Is the control system input variable at time k;
P k is an active instruction at time k, v k The wind speed at time k;
is the generalized state variable at time k, +.>Is a generalized state variable at time k+1, and N pairs of data (x k ,y k ) The initial data set is obtained for the matrix arrangement:
X=[x 1 x 2 … x N ],Y=[y 1 y 2 … y N ] (1.3)。
further, online dynamic modeling is performed by:
the linear model characterizes the current dynamic characteristics of the fan:
y k =Ax k (1.4)
solving an optimization problem by algebraic operation and obtaining a fan dynamic model:
is the generalized state variable at time k, +.>Is a generalized state variable at time k+1, and N pairs of data (x k ,y k ) The initial data set is obtained for the matrix arrangement:
X=[x 1 x 2 … x N ],Y=[y 1 y 2 … y N ] (1.3)
matrix a is the state transition matrix corresponding to the generalized state variable,is a matrix pseudo-inverse operation.
Further, a concentrated wind field control model is established by:
dynamic model of each power generation monomer in wind field:
is the fan rotating speed of the ith fan at the moment k+1, A i Is the state transition matrix corresponding to the generalized state variable of the ith fan,/the system is provided with a plurality of control modules>Is the fan rotating speed of the ith fan at the moment k, B i Is the input matrix of the ith fan, +.>Is the input variable of the control system of the ith fan at the moment k, and M is the total number of fans;
state variable χ of concentrated wind field control model at k moment k :
Input vector u of current wind speed and active command of each power generation unit at k moment k :
Obtaining a control model of the concentrated state vector from equation (1.6):
χ k+1 =Aχ k +Bu k (1.9)
constructing a state transition matrix of each power generation unit according to a diagonal form:
A 1 is a state transition matrix corresponding to the generalized state variable of the 1 st fan, A M Is a state transition matrix corresponding to the generalized state variable of the Mth fan, B 1 Is the input matrix of the 1 st fan, B M An input matrix of the Mth fan.
Further, online dynamic optimization is performed by:
wind field dynamic optimization control algorithm:
t is the length of a prediction interval of a model prediction control algorithm, and J is an objective function of the control algorithm:
Q k a semi-positive target coefficient matrix for the state variable at time k, R k A semi-positive target coefficient matrix for the input variable at time k,target coefficient vector for state variable at time k, < ->A target coefficient vector for the input variable at time k, E k Constraint coefficient matrix for state variable boundary at k moment, F k For the input variable at time kBoundary constraint coefficient matrix, b k The design of the coefficient matrix and the coefficient vector depends on a wind field dynamic optimization control target for the boundary constraint coefficient vector at the moment k;
wind power frequency modulation optimization target:
Active command of ith fan at k momentControl command at time k in maximum power tracking mode with respect to local controller of ith blower>If the wind farm adopts a load shedding mode, the active adjustment amount +.>
For the active instruction of the ith fan at the moment k under the load shedding working mode, R d To subtract the amplitude coefficient, the frequency offset is:
Δf=f meas -f ref (1.17)
the frequency f of the parallel point measurement meas Relative to the reference frequency f ref Deviation amount, K between df For the sag coefficient of the wind farm to the external power frequency characteristic curve, Q x Weight coefficients for balancing two optimization objectives;
time sequence constraint, input constraint and initial state constraint of concentrated state vector:
χ k+1 =Aχ k +Bu k ,k=0,1,…,T-1 (1.18)
χ k for the state variable of the concentrated wind field control model at the moment k, u k For the input variables of the concentrated wind field control model at time k, the matrices A and B are the state transition matrix and the input matrix of the concentrated wind field control model,for the fan rotating speed omega of the ith fan at the moment k min ,ω max Respectively a lower limit and an upper limit of the rotating speed of the fan, < ->For the active instruction of the ith fan at the moment k, P min ,P max Is the lower and upper bounds of the active instruction.
The invention provides a data-driven dynamic optimization control method for the frequency stability control of a high-proportion wind power independent power system, so that the wind power plant is stabilized and optimized in the internal transient process while participating in system frequency modulation. According to the wind field frequency optimization control method, an initial data set is constructed, on-line dynamic modeling and centralized wind field control modeling are carried out, and wind field frequency optimization control is carried out according to a fan on-line dynamic modeling result. According to the invention, the online equivalent dynamic model of the fan is obtained by fitting the historical data generated in the operation process of the fan, and the online equivalent dynamic model has the characteristic of pure data driving, so that the online equivalent dynamic model can be suitable for operation control of fans of different manufacturers under different working conditions, meanwhile, the algorithm has a pure linear form, the solving process is simple, the calculation load is small, the complex nonlinear optimization control problem is not required to be solved, the real-time operation requirement is met, and the balance of dynamic response precision and speed can be considered.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention may be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions of the prior art, the following description will briefly explain the drawings used in the embodiments or the description of the prior art, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic diagram of a wind farm frequency optimization control module according to an embodiment of the present invention.
Fig. 2 is a schematic diagram of a wind farm frequency optimization control implementation procedure according to an embodiment of the present invention.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present invention more apparent, the technical solutions of 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, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments of the present invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
A wind power plant dynamic frequency control method based on control process data fitting is shown in figure 1, and is based on a data-driven wind power plant frequency optimization control module relationship, and mainly comprises two modules, namely a fan online dynamic modeling module and a wind power plant frequency optimization control module; as shown in fig. 2, the wind field frequency optimization control implementation step is based on data driving.
The object of the method is a wind driven generator in a large wind farm group, and because the doubly-fed asynchronous fans in the wind power market at present occupy a larger market share, the fans are assumed to be doubly-fed asynchronous fans. On the other hand, because fans with similar in-site operation conditions and equipment parameters can be polymerized into one power generation monomer, the control strategy provided by the system aims at the active power output of each power generation monomer in the wind field.
For convenience of description, it is assumed that the converters of the doubly-fed asynchronous fan are all three-phase grid-connected, and single-phase or two-phase access can be performed according to requirements in actual engineering, so that implementation of a control method is not affected. And it is assumed that the power control loop of the converter has decoupled the active and reactive control by dq decomposition of the rotor current, so that the active output of the fan can be quickly adjusted following an external active command. Assuming that M fans generate electricity in total, the active instruction of the ith generating electricity monomer at the moment k is thatThe corresponding real-time wind speed is->The active command and the real-time wind speed form an input variable of a state equation, and the rotating speed of the fan is +.>Form the state variable of the fan.
(1) Construction of an initial dataset
The state equation of a single fan or a power generation unit is established from the frequency modulation control point of view as follows:
ω k+1 =f(ω k ,u k ) (1.1)
wherein f represents a state transition relation function, ω k The fan rotation speed at the moment k is u k Is the control system input variable at time k:
P k is an active instruction at time k, v k The wind speed at time k;
in order to establish a pure data driving method, only the accumulated data pairs (x k ,y k ),Is the generalized state variable at time k, +.>Is a generalized state variable at the time of k+1, comprises a current state and a current input, and forms a data pair which is accumulated to a preset value (when the step length is 0.1s, the preset value is larger than 6000) and comprises a time sequence characteristic of fan state transition, and the accumulated data pairs (x k ,y k ) The initial data set is obtained for the matrix arrangement as follows:
X=[x 1 x 2 … x N ],Y=[y 1 y 2 … y N ] (1.3)
(2) Online dynamic modeling
From the data driving perspective, the dominant dynamic equation can be deduced from an empirical model formula, and then nonlinear fitting is carried out on model parameters according to data formed by state tracks, but the method generally needs to measure a fan to some extent, otherwise, the fitting effect of the situation with more model parameters is often not ideal. Therefore, on-line dynamic modeling of the fan dynamic characteristics is considered through a pure data driving mode, specifically, N data pairs forward at the current moment form a data window, N depends on the size of the control step T, and typically, the product of NT is in a time scale of about one minute. The data set formed in this way not only contains the dynamic characteristics of the current running working point of the fan, but also avoids the nonlinear characteristics of the fan state under various wind speed levels, so that the current dynamic characteristics of the fan can be represented by a linear model:
y k =Ax k (1.4)
in order to restore the state transition matrix from a least squares perspective through the dataset, it is necessary to find matrix A such that AX-Y|| 2 At minimum, according to the theory of linear algebra, the optimization problem can be solved by the following algebraic operation and a fan dynamic model is obtained:
wherein the method comprises the steps ofRepresenting the pseudo-inverse of the matrix. It should be noted that, the matrix a here is a state transition matrix corresponding to a generalized state quantity, and for a general form required by a control model, sub-block interception is performed according to dimensions of the state quantity and the input quantity.
(3) Centralized wind field control model
According to the online dynamic modeling method, a dynamic model of each power generation monomer in the wind field can be obtained:
is the fan rotating speed of the ith fan at the moment k+1, A i Is the state transition matrix corresponding to the generalized state variable of the ith fan,/the system is provided with a plurality of control modules>Is the fan rotating speed of the ith fan at the moment k, B i Is the input matrix of the ith fan, +.>Is the input variable of the control system of the ith fan at the moment k, and M is the total number of fans;
on the basis, defining state vector χ in the k-moment concentrated wind field control model k :
Simultaneously defining an input vector u containing the current wind speed and the active command of each power generation unit at k moment k :
From equation (1.6) a control model corresponding to the concentrated state vector can be given:
χ k+1 =Aχ k +Bu k (1.9)
wherein the state transition matrix of each power generation unit is constructed in the following diagonal form:
a1 is a state transition matrix corresponding to a generalized state variable of the 1 st fan, AM is a state transition matrix corresponding to a generalized state variable of the M-th fan, B1 is an input matrix of the 1 st fan, and BM is an input matrix of the M-th fan.
It can be seen that the matrix A and the matrix B constructed in the mode have special sparse structures, and provide more convenient conditions for rapidly solving the online dynamic optimization control problem.
(4) Online dynamic optimization
On the basis of obtaining a linear dynamic model of the fan by a data driving method, according to the general form of a model predictive control architecture, a wind field dynamic optimization control algorithm can be expressed as follows:
where T is the prediction interval length of the model predictive control algorithm and the objective function J is of the form:
Q k a semi-positive target coefficient matrix for the state variable at time k, R k A semi-positive target coefficient matrix for the input variable at time k,target coefficient vector for state variable at time k, < ->A target coefficient vector for the input variable at time k, E k Constraint coefficient matrix for state variable boundary at k moment, F k Constraint coefficient matrix for boundary of input variable at k moment, b k And (3) for the boundary constraint coefficient vector at the moment k, the design of the coefficient matrix and the coefficient vector depends on a wind field dynamic optimization control target.
In the problem of wind power frequency modulation, in view of the fact that a control target needs to be optimized in two aspects of wind power plant active frequency modulation instruction following effect and rotational speed transient stability of a power generation unit, the transient stability degree of the wind power plant active frequency modulation instruction is measured by the fluctuation degree of the rotational speed of a fan, and the following optimization targets are given:
Active command of ith fan at k momentControl command +_at time k in maximum power tracking Mode (MPPT) with respect to local controller of ith blower>If the wind farm adopts load shedding mode (load), the active regulation is +.>/>
Wherein the method comprises the steps ofFor the active instruction of the ith fan at the moment k under the load shedding working mode, R d To subtract the amplitude coefficient, the frequency offset is defined as:
Δf=f meas -f ref (1.17)
representing the frequency f of measurement of the grid-connected point meas Relative to the reference frequency f ref Offset between and parameter K df For the sag coefficient of the wind farm to the external power frequency characteristic curve, Q x To balance the weight coefficients of the two optimization objectives.
The first term of the optimization target in the formula (1.14) enables the distribution of the frequency modulation tasks of the wind turbines to be carried out on the basis of considering the respective power generation level by adjusting the active command on the basis of the local control command, so that the overall distribution of the active output force corresponds to the local working condition of the wind turbines, and the wind power plant group provides frequency modulation service and simultaneously guarantees the overall wind energy conversion benefit of the wind power plant group. The second item of the optimization target enables the running state of the fan to be more stable by optimizing the fluctuation degree of the rotating speed of the fan, so that the mechanical fatigue of the fan, which is increased to fragile mechanical parts such as a gearbox and the like due to participation in a frequency modulation response process, is reduced, and the service life of the fan is prolonged.
On the basis of an optimization target, timing constraints, input constraints and initial state constraints of the concentrated state vector are given:
χ k+1 =Aχ k +Bu k ,k=0,1,…,T-1 (1.18)
χ k for the state variable of the concentrated wind field control model at the moment k, u k For the input variables of the concentrated wind field control model at time k, the matrices A and B are the state transition matrix and the input matrix of the concentrated wind field control model,for the fan rotating speed omega of the ith fan at the moment k min ,ω max Respectively a lower limit and an upper limit of the rotating speed of the fan, < ->For the active instruction of the ith fan at the moment k, P min ,P max Is the lower and upper bounds of the active instruction.
The complete wind field frequency dynamic optimization control method is formed by the formulas (1.14) to (1.21), and the constraint condition and the objective function have convexity, so that the whole model forms a convex optimized QP problem, and the existing optimization solver can accurately and rapidly solve the problem.
The invention provides a data-driven dynamic optimization control method for the frequency stability control of a high-proportion wind power independent power system, so that the wind power plant is stabilized and optimized in the internal transient process while participating in system frequency modulation. According to the wind field frequency optimization control method, an initial data set is constructed, on-line dynamic modeling and centralized wind field control modeling are carried out, and wind field frequency optimization control is carried out according to a fan on-line dynamic modeling result. According to the invention, the online equivalent dynamic model of the fan is obtained by fitting the historical data generated in the operation process of the fan, and the online equivalent dynamic model has the characteristic of pure data driving, so that the online equivalent dynamic model can be suitable for operation control of fans of different manufacturers under different working conditions, meanwhile, the algorithm has a pure linear form, the solving process is simple, the calculation load is small, the complex nonlinear optimization control problem is not required to be solved, the real-time operation requirement is met, and the balance of dynamic response precision and speed can be considered.
Although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention.
Claims (2)
1. The wind farm dynamic frequency control method based on control process data fitting is characterized by comprising the following steps of:
on-line dynamic modeling of the fan, wherein the on-line dynamic modeling of the fan is performed in the following manner: an initial dataset is constructed, wherein the initial dataset is constructed by:
establishing a state equation of a single fan or a power generation unit:
ω k+1 =f(ω k ,u k ) (1.1)
f is a state transition relation function, ω k The fan rotation speed at the moment k is u k Is the control system input variable at time k;
P k is an active instruction at time k, v k The wind speed at time k;
is the generalized state variable at time k, +.>Is a generalized state variable at time k+1, and N pairs of data (x k ,y k ) The initial data set is obtained for the matrix arrangement:
X=[x 1 x 2 … x N ],Y=[y 1 y 2 … y N ] (1.3);
online dynamic modeling, wherein online dynamic modeling is performed by:
the linear model characterizes the current dynamic characteristics of the fan:
y k =Ax k (1.4)
solving an optimization problem by algebraic operation and obtaining a fan dynamic model:
is the generalized state variable at time k, +.>Is a generalized state variable at time k+1, and N pairs of data (x k ,y k ) The initial data set is obtained for the matrix arrangement:
X=[x 1 x 2 … x N ],Y=[y 1 y 2 … y N ] (1.3)
matrix a is the state transition matrix corresponding to the generalized state variable,is a matrix pseudo-inverse operation;
establishing a concentrated wind field control model, wherein the concentrated wind field control model is established by the following steps:
dynamic model of each power generation monomer in wind field:
is the fan rotating speed of the ith fan at the moment k+1, A i Is the state transition matrix corresponding to the generalized state variable of the ith fan,/the system is provided with a plurality of control modules>Is the fan rotating speed of the ith fan at the moment k, B i Is the input matrix of the ith fan, +.>Is the input variable of the control system of the ith fan at the moment k, and M is the total number of fans;
state variable χ of concentrated wind field control model at k moment k :
Input vector u of current wind speed and active command of each power generation unit at k moment k :
Obtaining a control model of the concentrated state vector from equation (1.6):
χ k+1 =Aχ k +Bu k (1.9)
constructing a state transition matrix of each power generation unit according to a diagonal form:
A 1 is a state transition matrix corresponding to the generalized state variable of the 1 st fan, A M Is a state transition matrix corresponding to the generalized state variable of the Mth fan, B 1 Is the input matrix of the 1 st fan, B M An input matrix of the M-th fan;
and optimally controlling the wind field frequency according to the online dynamic modeling result of the fan.
2. The control method according to claim 1, characterized in that the on-line dynamic optimization is performed by:
wind field dynamic optimization control algorithm:
t is the length of a prediction interval of a model prediction control algorithm, and J is an objective function of the control algorithm:
Q k a semi-positive target coefficient matrix for the state variable at time k, R k A semi-positive target coefficient matrix for the input variable at time k,target coefficient vector for state variable at time k, < ->A target coefficient vector for the input variable at time k, E k Constraint coefficient matrix for state variable boundary at k moment, F k Constraint coefficient matrix for boundary of input variable at k moment, b k The design of the coefficient matrix and the coefficient vector depends on a wind field dynamic optimization control target for the boundary constraint coefficient vector at the moment k;
wind power frequency modulation optimization target:
Active command of ith fan at k momentRelative toControl command of local controller of ith fan at k moment in maximum power tracking mode +.>If the wind farm adopts a load shedding mode, the active adjustment amount +.>
For the active instruction of the ith fan at the moment k under the load shedding working mode, R d To subtract the amplitude coefficient, the frequency offset is:
Δf=f meas -f ref (1.17)
the frequency f of the parallel point measurement meas Relative to the reference frequency f ref Deviation amount, K between df For the sag coefficient of the wind farm to the external power frequency characteristic curve, Q x Weight coefficients for balancing two optimization objectives;
time sequence constraint, input constraint and initial state constraint of concentrated state vector:
χ k+1 =Aχ k +Bu k ,k=0,1,…,T-1 (1.18)
χ k for the state variable of the concentrated wind field control model at the moment k, u k For the input variables of the concentrated wind field control model at time k, the matrices A and B are the state transition matrix and the input matrix of the concentrated wind field control model,for the fan rotating speed omega of the ith fan at the moment k min ,ω max Respectively a lower limit and an upper limit of the rotating speed of the fan, < ->For the active instruction of the ith fan at the moment k, P min ,P max Is the lower and upper bounds of the active instruction. />
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011398967.1A CN112636366B (en) | 2020-12-01 | 2020-12-01 | Wind power plant dynamic frequency control method based on control process data fitting |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011398967.1A CN112636366B (en) | 2020-12-01 | 2020-12-01 | Wind power plant dynamic frequency control method based on control process data fitting |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112636366A CN112636366A (en) | 2021-04-09 |
CN112636366B true CN112636366B (en) | 2023-05-16 |
Family
ID=75307559
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011398967.1A Active CN112636366B (en) | 2020-12-01 | 2020-12-01 | Wind power plant dynamic frequency control method based on control process data fitting |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112636366B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114188991B (en) * | 2021-11-03 | 2024-03-12 | 国网吉林省电力有限公司 | Wind farm primary frequency modulation model identification method and device, electronic equipment and storage medium |
CN114583767B (en) * | 2022-03-10 | 2023-03-17 | 中国电力科学研究院有限公司 | Data-driven wind power plant frequency modulation response characteristic modeling method and system |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103887815B (en) * | 2014-02-21 | 2016-03-02 | 华南理工大学 | Based on wind energy turbine set parameter identification and the Dynamic Equivalence of service data |
EP3096004A1 (en) * | 2015-05-18 | 2016-11-23 | ABB Technology AG | Wind farm inertial response |
CN109085755A (en) * | 2018-08-14 | 2018-12-25 | 中国电力科学研究院有限公司 | The advanced control method for frequency of wind-powered electricity generation group system and system based on distributed model predictive control |
-
2020
- 2020-12-01 CN CN202011398967.1A patent/CN112636366B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN112636366A (en) | 2021-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Abdelbaky et al. | Design and implementation of partial offline fuzzy model-predictive pitch controller for large-scale wind-turbines | |
Zhao et al. | Fatigue load sensitivity-based optimal active power dispatch for wind farms | |
Mansour et al. | A grid-connected variable-speed wind generator driving a fuzzy-controlled PMSG and associated to a flywheel energy storage system | |
CN112632774B (en) | Data-driven wind power plant frequency control method based on dynamic modal decomposition | |
CN106786677A (en) | A kind of interconnected electric power system distributed dynamic matrix frequency control method of meter and wind-powered electricity generation | |
CN112636366B (en) | Wind power plant dynamic frequency control method based on control process data fitting | |
Zamzoum et al. | Performance analysis of a robust adaptive fuzzy logic controller for wind turbine power limitation | |
CN109296500A (en) | Maximal wind-energy capture method based on robust control theory | |
Barrera-Cardenas et al. | Optimal LQG controller for variable speed wind turbine based on genetic algorithms | |
Hamid et al. | Optimal MPPT and BES control for grid-tied DFIG-based wind energy conversion system | |
Geng et al. | Robust pitch controller for output power levelling of variable-speed variable-pitch wind turbine generator systems | |
Song et al. | An overview of renewable wind energy conversion system modeling and control | |
Hosseini et al. | Control of pitch angle in wind turbine based on doubly fed induction generator using fuzzy logic method | |
Chojaa et al. | Optimization of DFIG wind turbine power quality through adaptive fuzzy control | |
Chakib et al. | Comparative study of active disturbance rejection control with RST control for variable wind speed turbine based on doubly fed induction generator connected to the grid | |
Bernard et al. | Ant-based optimal tuning of PID controllers for load frequency control in power systems | |
Tidjani et al. | Augmented robust TS fuzzy control based PMSG wind turbine improved with H∞ performance | |
CN109611274A (en) | A kind of wind power generating set high wind speed area LQG optimal control method | |
Manna et al. | A review of control techniques for wind energy conversion system | |
Wang et al. | Wind power compound model-free adaptive predictive control based on full wind speed | |
Yao et al. | Pitch regulated LQG controller design for variable speed wind turbine | |
Puangdownreong | Optimal PIλDµ controller design based on spiritual search for wind turbine systems | |
CN115642613A (en) | Wind power plant robust frequency coordination control method, system and computer storage medium | |
CN112072693B (en) | High-dimensional multi-fractional order optimization method for high-dimensional multi-fractional order doubly-fed fan control | |
CN113922384A (en) | Wind power plant distributed reactive voltage optimization coordination control 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 |