CN114069599A - Polymerization temperature control load improvement model prediction control under preset performance condition - Google Patents
Polymerization temperature control load improvement model prediction control under preset performance condition Download PDFInfo
- Publication number
- CN114069599A CN114069599A CN202010754960.2A CN202010754960A CN114069599A CN 114069599 A CN114069599 A CN 114069599A CN 202010754960 A CN202010754960 A CN 202010754960A CN 114069599 A CN114069599 A CN 114069599A
- Authority
- CN
- China
- Prior art keywords
- temperature
- model
- bilinear
- temperature control
- load
- 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
- 238000006116 polymerization reaction Methods 0.000 title claims abstract description 30
- 230000006872 improvement Effects 0.000 title claims description 6
- 238000000034 method Methods 0.000 claims abstract description 26
- 239000011159 matrix material Substances 0.000 claims description 27
- 239000000463 material Substances 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 8
- 239000000126 substance Substances 0.000 claims description 7
- 230000007704 transition Effects 0.000 claims description 3
- 230000002776 aggregation Effects 0.000 abstract description 9
- 238000004220 aggregation Methods 0.000 abstract description 9
- 238000011217 control strategy Methods 0.000 abstract description 7
- 238000004088 simulation Methods 0.000 abstract description 7
- 238000011426 transformation method Methods 0.000 abstract 1
- 238000004422 calculation algorithm Methods 0.000 description 5
- 238000010248 power generation Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 230000001052 transient effect Effects 0.000 description 3
- 230000033228 biological regulation Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000178 monomer Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000000087 stabilizing effect Effects 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000004931 aggregating effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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/003—Load forecast, e.g. methods or systems for forecasting future load demand
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
-
- 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]
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S50/00—Market activities related to the operation of systems integrating technologies related to power network operation or related to communication or information technologies
- Y04S50/16—Energy services, e.g. dispersed generation or demand or load or energy savings aggregation
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Power Engineering (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Feedback Control In General (AREA)
- Air Conditioning Control Device (AREA)
Abstract
The method for providing the load tracking auxiliary service by using the aggregated temperature control load on the demand side to stabilize the power fluctuation of the new energy is an economic and effective method. The invention discloses an improved model predictive control strategy with a preset performance based on aggregation temperature control load balancing new energy power fluctuation. Firstly, expanding an original temperature control load bilinear polymerization model through a second-order equivalent thermodynamic parameter model to obtain an improved temperature control load bilinear polymerization model so as to optimize the influence of the accumulated error of the original bilinear polymerization model; then, in order to ensure the specified performance of the output tracking error, an error transformation method is utilized, a preset performance controller is designed for the improved temperature control load bilinear polymerization model based on model prediction control of a Lyapunov function, and the stability of the control strategy is proved; finally, simulation results show that the improved temperature control load bilinear model has higher precision, the proposed control method has high convergence rate, and the tracking error can be controlled within the limit of the prediction performance function.
Description
Technical Field
The invention belongs to the field of power grid auxiliary service and demand response control, and particularly relates to an improved model predictive control strategy with preset performance based on aggregated temperature control load balancing new energy power fluctuation.
Background
The large-scale grid connection of the intermittent new energy sources brings challenges to the maintenance of power balance between a power source and a load of a power system. Usually, a controllable power supply such as thermal power supplies provides auxiliary services, but on one hand, the operation efficiency of the system is reduced, and on the other hand, the operation cost is increased due to deep peak shaving of the thermal power unit. Previous researches show that Thermal Controlled Loads (TCLs) of demand side resources such as air conditioners, electric water heaters and the like have large regulation potential after a large amount of polymerization. In addition, the Aggregation Thermal Controlled Load (ATCL) has the advantages of large energy storage capacity, quick response and the like, and becomes an ideal choice for providing auxiliary service for the power grid on the demand side.
The establishment of an accurate and practical temperature control load aggregation model is the basis for researching the participation of the temperature control load in the regulation and control of the power system. The current temperature control load aggregation model for control mainly comprises a state space model, a state sequence model, a Fokker-Planck equation model, a bilinear model and the like. The bilinear model which adopts the finite difference method to disperse the temperature in the finite range is a method which has higher precision and is easy to combine with the dispatching of the power system. The output power is adjusted by changing the temperature set point. Most of the temperature-controlled load bilinear polymerization models in the prior art are based on a first-order Equivalent Thermal Parameter (ETP) model to simulate the thermodynamic process, and the material temperature is ignored. Obviously, due to the coupled thermodynamic characteristics of the material temperature and the indoor temperature, the first-order equivalent thermodynamic parameter model is not accurate enough in practical application. Therefore, it is necessary to introduce the material temperature to obtain an aggregation model which accurately reflects the transient and steady-state response of the cluster temperature control load.
At present, many documents research control methods based on temperature-controlled load polymerization models. Such as back-thrust control, sliding mode control, Model Predictive Control (MPC), and the like. MPC is widely used due to its robustness and control accuracy. However, the conventional MPC is optimized repeatedly each time sampling, which requires a large amount of calculation, and causes a certain execution delay, which may reduce the control performance. Moreover, the above-mentioned aggregate temperature-controlled load power control method hardly enables the tracking error to be converged quickly according to a preset performance range. With respect to predetermined performance control, this method has been widely used in other fields. Here, predetermined properties are applied to the polymerization temperature-controlled load control.
Disclosure of Invention
In order to solve the problems, the invention improves a primary TCL Bilinear Aggregation Model (TBAM) and provides an improved MPC control strategy with preset performance. The problem that the traditional MPC is large in calculation amount is solved, and the tracking error is made to be converged in a pre-designed performance range.
The innovation points are embodied in the following two aspects: (1) an original temperature control load bilinear model is improved by introducing the material temperature, and a more accurate improved temperature control load bilinear aggregation model (ITBAM) in practical application is obtained; (2) based on an improved temperature control load bilinear polymerization model, an MPC control scheme with a preset performance and based on a Lyapunov function is designed, and transient and steady-state performances required by a tracking error are guaranteed.
In order to achieve the above object, the present invention provides a method for predictive control of a polymerization temperature-controlled load improvement model under a preset performance condition, comprising the steps of:
(1) expanding the original temperature control load bilinear polymerization model through a second-order equivalent thermodynamic parameter model to obtain an improved temperature control load bilinear polymerization model;
(2) aiming at the improved temperature control load bilinear polymerization model, a preset performance controller is designed based on model prediction control of a Lyapunov function.
According to the step (1), the original TCL bilinear polymerization model is expanded through a second-order ETP model to obtain an improved TCL bilinear polymerization model:
wherein,x (t) is the number of temperature control loads in each temperature interval, A is the state matrix, B is the input matrix, C is the output matrix, U (t) is the control input, CaIs the heat capacity of indoor air, RaIs the thermal resistance of indoor air, CmIs the heat capacity of indoor material, RmIs the thermal resistance of the indoor substance,as a derivative of the temperature set point, TsetIs a temperature set value, and is used as a temperature setting value,to a desired temperature set point, Tm0Is the initial room temperature, PT(T) Total output Power of the aggregate temperature controlled load,. DELTA.TsetFor variation of the temperature setting, Δ TmIs the change of the temperature of the indoor material u1(t) and u2(t) respectively reflecting the real-time influence and the accumulated influence of the temperature set value change on the dynamic process of the temperature control load, u3(t) reflects the influence of heterogeneity on the dynamic process of temperature-controlled load,as an estimate of the temperature of the room mass, Tm(t) is the room temperature.
Due to the temperature T of the indoor materialm(t) is not easily obtained in practical applications, and its estimated value can be usedInstead of this, the user can,the calculation can be made using the following formula:
the step (2) is based on an improved temperature control load bilinear polymerization model, and a preset performance controller is designed based on model prediction control of a Lyapunov function:
wherein, ρ(t)=(ρ0-ρ∞)e-rt+ρ∞,-δ1ρ(t)<ey(t)<δ2ρ(t),ρ0p (0), u (t) is the resulting optimal control input, Pref(t) is a given power reference value, s is a converted error variable, α is a virtual control quantity, ξ is an intermediate variable, ey(t) is the tracking error, X (t) is the number of temperature control loads in each temperature interval, A is the state matrix, B is the input matrix, C is the output matrix, xi is the intermediate variable and xi > 0, epsilon is a constant greater than 0, rho (t) is the selected predetermined performance function, rho (0) is the value of rho (t) at zero time, k is the value of1、δ1And delta2Are predetermined performance parameters and are all greater than 0, r is a time coefficient, rho0For a predetermined performance function initial value, p∞Is a steady state value of a predetermined performance function, and p0>ρ∞,Are the transition variables.
Drawings
FIG. 1 is a comparison of Monte Carlo simulations of a second order ETP model and a first order ETP model;
FIG. 2 is a comparison of simulation results of TBAM, ITBAM and second order ETP Monte Carlo models;
FIG. 3 is a schematic view of an embodiment;
FIG. 4 is a diagram of an improved MPC algorithm power tracking error with predetermined performance;
FIG. 5 is an MPC algorithm power tracking error based on a Lyapunov function;
FIG. 6 is a system power generation and utilization curve;
fig. 7 is a tracking error curve for different control algorithms.
Detailed Description
The invention is realized by the following technical scheme:
modeling of TCL bilinear polymerization model
1.1 thermodynamic model of monomer TCL
Most literature often describes thermodynamic processes using a first-order ETP model, which can be described as:
because the second-order dynamic effect caused by the temperature of the substance is neglected, the first-order ETP model cannot accurately describe the real TCL dynamic process caused by the change of the temperature set value. Therefore, TCL is described using a second order ETP model, also known as a two-mass model, which takes changes in air temperature and mass temperature as observed two state variables. The specific differentiation process is as follows:
wherein T is the indoor temperature, C is the indoor average heat capacity, R is the indoor average heat resistance, TaIs the temperature of the indoor air, TmIs the temperature of the indoor material, T∞Is the outdoor ambient temperature, P is the TCL rated power of the monomer, TmaxAnd TminRespectively, the upper and lower limits of the temperature limit, TsetIs the temperature set point, δdbM (t) is a temperature dead zone, m (t) is a switching variable, Δ t is a time interval, CaFor heat capacity of indoor air,RaIs the thermal resistance of indoor air, CmIs the heat capacity of indoor material, RmIs the thermal resistance of the indoor substance,
1.2 improved bilinear temperature control load aggregation model
A control-oriented temperature control load bilinear polymerization model (TBAM) established based on a first-order ETP model is as follows:
wherein X (t) is the number of temperature control loads in each temperature interval, A is a state matrix, B is an input matrix, C is an output matrix, u (t) is a control input, PT(t) is the total output power of the aggregate temperature control load.
The TBAM model is only suitable for a scene with a small change range of the load temperature set value, and once the temperature set value deviates from the expected value to be too large, the bilinear polymerization model fails under long-time operation. Therefore, the TBAM model is extended by using a second-order ETP model, and a more accurate improved temperature-controlled load bilinear aggregation model (ITBAM) is obtained as follows:
wherein,x (t) is the number of temperature control loads in each temperature interval, A is the state matrix, B is the input matrix, C is the output matrix, U (t) is the control input, CaIs the heat capacity of indoor air, RaIs the thermal resistance of indoor air, CmIs the heat capacity of indoor material, RmIs the thermal resistance of the indoor substance,as a derivative of the temperature set point, TsetIs a temperature set value, and is used as a temperature setting value,to a desired temperature set point, Tm0Is the initial room temperature, PT(T) Total output Power of the aggregate temperature controlled load,. DELTA.TsetFor variation of the temperature setting, Δ TmIs the change of the temperature of the indoor material u1(t) and u2(t) respectively reflecting the real-time influence and the accumulated influence of the temperature set value change on the dynamic process of the temperature control load, u3(t) reflects the influence of heterogeneity on the dynamic process of temperature-controlled load,as an estimate of the temperature of the room mass, Tm(t) is the room temperature.
2. Polymerization temperature control load improvement model prediction control under preset performance condition
2.1 control problem description
The control object of the invention is to adjust the ATCL output power to track the given reference power by using the proposed MPC control scheme based on Lyapunov function with the predetermined performance based on the improved temperature controlled load bilinear polymerization model.
2.2 error conversion
Defining a tracking error ey(t)=PT(t)-Pref(t), the predetermined property of which can be described as:
-δ1ρ(t)<ey(t)<δ2ρ(t) (5)
wherein, PT(t) Total output Power of the aggregate temperature controlled load, PrefIs the reference power. ρ (t) ═ p (ρ0-ρ∞)e-rt+ρ∞For a selected performance function, δ1,δ2A predetermined performance parameter greater than 0; rho0ρ (0) is a value at time ρ (t) 0. r is a time coefficient, ρ0For a predetermined performance function initial value, p∞Is a steady state value of a predetermined performance function, and p0>ρ∞。
The constrained tracking error is converted to an unconstrained error variable. In particular, we can define:
ey(t)=ρ(t)φ(s(t)) (6)
wherein s is a conversion error; phi(s) belongs to (-delta)1,δ2) Is a smooth monotonically increasing function.
From φ(s), s can be defined as:
As long as s (t) is bounded, ey(t) satisfies the predetermined property described in the above formula (6).
2.3 Lyapunov function based MPC controller design with predetermined Performance
Using the state space model of equation (4) as the prediction model, it is rewritten as:
wherein X (t) is the number of temperature control loads in each temperature interval, A is a state matrix, B is an input matrix, C is an output matrix, U (t) is a control input, PT(t) is the total output power of the aggregate temperature control load.
From equations (8) and (9), it is possible to obtain:
The derivative of V is:
in the formula, s is an error variable,for the derivative of the error variable, X (t) is the number of temperature control loads in each temperature interval, A is the state matrix, B is the input matrix, C is the output matrix, U (t) is the control input, Pref(t) is reference power, ey(t) is the tracking error, ρ (t) is the performance function,is the derivative of the performance function.
The optimal control law is designed as follows:
by using the theorem 1, s ξ cx (t) u (t) in the formula (11) is represented as:
sξCX(t)U(t)≤|sξCBX(t)U(t)|≤ε-sα。 (13)
where ε is a constant greater than 0, s is an error variable, X (t) is the number of temperature control loads in each temperature interval, A is a state matrix, B is an input matrix, C is an output matrix, U (t) is a control input,is the derivative of the reference power, ey(t) is the tracking error, ρ (t) is the performance function,is the derivative of the performance function, xi is the intermediate variable, k1Is a parameter greater than 0.
The theory 1 is as follows:
for any x ∈ R and a constant ε > 0, the following holds:
substituting the above formula (12) into the above formula (11) can obtain:
in the formula, k1Is a parameter greater than 0, η ═ k1ε is a constant greater than 0 and s is the error variable.
As can be obtained from equation (15) above, s is bounded, thus ey(t) satisfies the above-specified properties.
3. Improved TCL bilinear polymerization model accuracy analysis
In order to verify the accuracy of ITBAM, 1000 TCL devices are selected, and a first-order ETP model, a second-order ETP model, TBAM and ITBAM are simulated and compared through Monte Carlo. The simulation parameters are shown in table 1. The simulation results are shown in fig. 1-2.
As shown in fig. 1, when the temperature setting value is changed, the second-order ETP model has a longer response period and a larger fluctuation amplitude than the first-order ETP model under the same parameters, which is caused by the thermal capacitance and thermal resistance properties exhibited by the introduced indoor substance. The coupling characteristic of the indoor air temperature and the indoor substance temperature is considered, and the second-order ETP model can reflect the real use condition better.
As shown in fig. 2, the ITBAM effectively improves the accumulated error of the bilinear model under long-time operation, improves the accuracy of model description, and proves that the proposed ITBAM can more accurately describe the dynamic process of aggregating TCLs.
4. Simulation analysis of control strategies
To verify the effectiveness of the proposed control strategy, a schematic diagram of an example of a multi-level hierarchical architecture is adopted, as shown in fig. 3. Including conventional generator sets, renewable energy sources such as wind farms and photovoltaic generation and TCL aggregates containing 10000 devices. The dispatching plan is directly issued by the dispatching center of the upper layer; the distributed reference power is adjusted and tracked by the TCL aggregate of the middle layer; the lower layer aggregates a large number of TCLs to participate in the load tracking service using the proposed controller.
The effectiveness of the control strategy for stabilizing the power fluctuation provided by the invention is further analyzed and verified in two calculation examples implemented below by using MPC based on Lyapunov function as a comparison algorithm.
Example 1: tracking of a given power value
The goal was to provide 40MW of power in 30 minutes, with 6 initial states e selected separatelyy(0)=[0.0175,0.0451,0.0947,-0.0239,-0.0515,-0.0929]With a design parameter of rho0=0.1009,δ1=1,δ21, r 0.9 and ρ∞0.009. The simulation results are shown in fig. 4-5. As can be seen from fig. 4-5, the control method proposed herein can adjust the tracking error to zero and can adapt to the predetermined performance requirements at different initial conditions. In particular, compared with a model prediction control algorithm based on the Lyapunov function, the control method has faster transient convergence and can control the tracking error within the limit of a prediction performance function.
Example 2: stabilizing new energy power fluctuation
The ATCMs are utilized to carry out power fluctuation stabilization, and reference power (expected ATCMs power) is provided by a dispatching center and can be expressed as the total generated power containing new energy minus rigid load required power. The electricity generation and utilization in the system within two hours is shown in fig. 6.
As can be seen from fig. 6, the total load power demand of the system in the first two hours is greater than the total power generation power, because the uncontrollable load power demand is greater, even if the conventional energy generator set is adjusted, the load demand cannot be met due to the randomness of new energy power generation, and at this time, the ATCLs needs to be reduced to achieve the balance of the power consumption; and in the last two hours, the total generated power is larger than the total load power, because the new energy power generation is suddenly increased and is far larger than the load power demand, the new energy power generation can be consumed by increasing the ATCL power so as to meet the demand and supply balance. The desired ATCLs power is expressed as the total generated power containing new energy minus the demanded power without TCL loads.
To track the reference power, a parameter ρ is selected0=0.12,δ1=1,δ21, r 0.9 and ρ∞0.007. The simulation results are shown in fig. 7.
As can be seen from fig. 7, the proposed control method has a smaller steady-state error (about 0.002p.u) and faster convergence performance compared to the Lyapunov function-based MPC method with a steady-state error of about 0.005 p.u.
Claims (3)
1. The method is characterized by comprising the following steps of:
(1) expanding the original temperature control load bilinear polymerization model through a second-order equivalent thermodynamic parameter model to obtain an improved temperature control load bilinear polymerization model;
(2) aiming at the improved temperature control load bilinear polymerization model, a preset performance controller is designed based on model prediction control of a Lyapunov function.
2. The method for predicting and controlling the aggregated temperature-controlled load improvement model under the preset performance condition of claim 1, wherein: in the step (1), the original bilinear polymerization model is expanded, and the improved temperature control load bilinear polymerization model is obtained by:
wherein,x (t) is the number of temperature control loads in each temperature interval, A is the state matrix, B is the input matrix, C is the output matrix, U (t) is the control input, CaIs the heat capacity of indoor air, RaIs the thermal resistance of indoor air, CmIs the heat capacity of indoor material, RmIs the thermal resistance of the indoor substance,as a derivative of the temperature set point, TsetIs a temperature set value, and is used as a temperature setting value,to a desired temperature set point, Tm0Is the initial room temperature, PT(T) Total output Power of the aggregate temperature controlled load,. DELTA.TsetFor variation of the temperature setting, Δ TmIs the change of the temperature of the indoor material u1(t) and u2(t) respectively reflecting the real-time influence and the accumulated influence of the temperature set value change on the dynamic process of the temperature control load, u3(t) reflects the influence of heterogeneity on the dynamic process of temperature-controlled load,as an estimate of the temperature of the room mass, Tm(t) is the room temperature.
Due to the temperature T of the indoor materialm(t) is not easily obtained in practical applications, and its estimated value can be usedInstead of this, the user can,can be calculated by:
3. The method for predicting and controlling the aggregated temperature-controlled load improvement model under the preset performance condition of claim 1, wherein: the step (2) is based on an improved temperature control load bilinear polymerization model, and a preset performance controller is designed based on model prediction control of a Lyapunov function:
wherein, ρ(t)=(ρ0-ρ∞)e-rt+ρ∞,-δ1ρ(t)<ey(t)<δ2ρ(t),ρ0p (0), u (t) is the resulting optimal control input, Pref(t) is a given power reference value, s is a converted error variable, α is a virtual control quantity, ξ is an intermediate variable, ey(t) is the tracking error, X (t) is the number of temperature control loads in each temperature interval, A is the state matrix, B is the input matrix, C is the output matrix, xi is the intermediate variable and xi > 0, epsilon is a constant greater than 0, rho (t) is the selected predetermined performance function, rho (0) is the value of rho (t) at zero time, k is the value of1、δ1And delta2Are predetermined performance parameters and are all greater than 0, r is a time coefficient, rho0For a predetermined performance function initial value, p∞Is a steady state value of a predetermined performance function, and p0>ρ∞,Are the transition variables.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010754960.2A CN114069599B (en) | 2020-07-29 | 2020-07-29 | Polymerization temperature control load improvement model predictive control under preset performance condition |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010754960.2A CN114069599B (en) | 2020-07-29 | 2020-07-29 | Polymerization temperature control load improvement model predictive control under preset performance condition |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114069599A true CN114069599A (en) | 2022-02-18 |
CN114069599B CN114069599B (en) | 2023-12-19 |
Family
ID=80227884
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010754960.2A Active CN114069599B (en) | 2020-07-29 | 2020-07-29 | Polymerization temperature control load improvement model predictive control under preset performance condition |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114069599B (en) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107800157A (en) * | 2017-11-14 | 2018-03-13 | 武汉大学 | The virtual power plant dual-layer optimization dispatching method of the temperature control load containing polymerization and new energy |
CN108490791A (en) * | 2018-05-10 | 2018-09-04 | 燕山大学 | Temperature control load Cost Controlling Policy |
CN109787249A (en) * | 2018-12-21 | 2019-05-21 | 华北电力大学(保定) | The association system model of a kind of polymerization temperature control load and generator excitation and its counter push away control method for coordinating |
CN109787251A (en) * | 2019-01-16 | 2019-05-21 | 华北电力大学(保定) | A kind of cluster temperature control Load aggregation model, Parameter identification and counter push away control method |
-
2020
- 2020-07-29 CN CN202010754960.2A patent/CN114069599B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107800157A (en) * | 2017-11-14 | 2018-03-13 | 武汉大学 | The virtual power plant dual-layer optimization dispatching method of the temperature control load containing polymerization and new energy |
CN108490791A (en) * | 2018-05-10 | 2018-09-04 | 燕山大学 | Temperature control load Cost Controlling Policy |
CN109787249A (en) * | 2018-12-21 | 2019-05-21 | 华北电力大学(保定) | The association system model of a kind of polymerization temperature control load and generator excitation and its counter push away control method for coordinating |
CN109787251A (en) * | 2019-01-16 | 2019-05-21 | 华北电力大学(保定) | A kind of cluster temperature control Load aggregation model, Parameter identification and counter push away control method |
Non-Patent Citations (1)
Title |
---|
余洋 等: "基于反推控制的聚合温控负荷与发电机励磁协调控制", 《电力***自动化》, pages 101 - 124 * |
Also Published As
Publication number | Publication date |
---|---|
CN114069599B (en) | 2023-12-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109799851B (en) | User side load response method based on load cluster temperature regulation and control | |
US11277009B2 (en) | Method and apparatus for smoothing link-line power of electrothermal microgrid using thermal storage heat pump | |
CN105042800A (en) | Variable-frequency air conditioner load modeling and operation controlling method based on demand responses | |
CN107860057B (en) | Heat load economic optimization scheduling method for cogeneration heating system | |
CN113757931B (en) | Air conditioner control method and system | |
CN110739699A (en) | distribution network voltage regulation method based on air conditioner | |
CN113159983A (en) | Coordination control method for off-grid phosgene cogeneration comprehensive energy system | |
CN108490791B (en) | Temperature controlled load cost control strategy | |
CN113078629A (en) | Aggregate power distribution model for cluster temperature control load aggregate power regulation and control and distributed consistency control method | |
CN112994036B (en) | Temperature control load participation micro-grid regulation and control method and system based on model prediction | |
CN113028604A (en) | Temperature regulation control and scheduling method based on aggregated air conditioner | |
CN112838580B (en) | Improved heterogeneous temperature control load aggregation model and multi-target coordination control method | |
CN114069599A (en) | Polymerization temperature control load improvement model prediction control under preset performance condition | |
CN111242412A (en) | Thermal control load cluster cooperative management and control method based on demand response | |
CN108224692B (en) | Consider the air-conditioning flexible control responding ability prediction technique of outside air temperature prediction error | |
CN111193261A (en) | Day-ahead optimization method of multi-energy flow system based on building equivalent heat energy storage | |
CN115289518A (en) | Heating system thermal and hydraulic balance control method and system | |
CN114079293B (en) | New energy power fluctuation smoothing method of electric-gas interconnection system considering temperature control load and micro gas turbine | |
CN109787249A (en) | The association system model of a kind of polymerization temperature control load and generator excitation and its counter push away control method for coordinating | |
CN115377963B (en) | Novel temperature set value control method for heterogeneous cluster air conditioner load | |
CN113188239B (en) | Air conditioner load cluster modeling method based on deep belief network | |
CN111082442A (en) | Energy storage capacity optimal configuration method based on improved FPA | |
Cao et al. | Online optimization of heat pump systems for building heating based on Lyapunov method | |
CN117592861A (en) | Load demand response multi-state target quantitative evaluation device, evaluation and verification method | |
CN114498651A (en) | Hybrid load cluster management and control method and 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 |