CN105470957B - Power grid load modeling method for production simulation - Google Patents

Power grid load modeling method for production simulation Download PDF

Info

Publication number
CN105470957B
CN105470957B CN201511004804.XA CN201511004804A CN105470957B CN 105470957 B CN105470957 B CN 105470957B CN 201511004804 A CN201511004804 A CN 201511004804A CN 105470957 B CN105470957 B CN 105470957B
Authority
CN
China
Prior art keywords
load
power
unit
time
establishing
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
Application number
CN201511004804.XA
Other languages
Chinese (zh)
Other versions
CN105470957A (en
Inventor
刘纯
黄越辉
礼晓飞
曹阳
王跃峰
刘德伟
高云峰
许晓艳
马烁
许彦平
张楠
杨硕
李驰
王晶
潘霄峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
Original Assignee
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI filed Critical State Grid Corp of China SGCC
Priority to CN201511004804.XA priority Critical patent/CN105470957B/en
Publication of CN105470957A publication Critical patent/CN105470957A/en
Application granted granted Critical
Publication of CN105470957B publication Critical patent/CN105470957B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications

Abstract

The invention provides a power grid load modeling method for production simulation, which comprises the steps of establishing a load transfer power upper limit model of each time section; establishing a load transfer electric quantity upper limit model in each period; establishing different load clearing mode models; establishing a regional load balance model based on load side peak regulation; and establishing a load side peak regulation model for improving the wind power absorption capacity, and optimizing the load side peak regulation electric quantity of each time section. According to the method, on the basis of ensuring the calculation efficiency, the wind curtailment electric quantity of the wind power is effectively reduced through effective management of the load, the optimization result provides guidance and suggestions for power grid operators, the active power balance difficulty of a power system is effectively reduced, and the efficient and stable operation of the power grid is further ensured.

Description

Power grid load modeling method for production simulation
Technical Field
The invention relates to the field of new energy power generation, in particular to a power grid load modeling method for improving wind power absorption capacity.
Background
With the continuous and high-speed development of economy, gradual exhaustion of fossil fuel and continuous aggravation of pollution to the environment and greenhouse effect of fossil fuel in China, the government of China highly pays attention to the development and utilization of new energy based on a sustainable development strategy, develops and utilizes the new energy as an important measure for treating atmospheric pollution, adjusting energy structure and changing economic development, and takes wind power generation and photovoltaic power generation as one of main modes for developing and utilizing the new energy.
Because wind energy resources are thin and low in space-time energy density, and cannot be enriched, transported and stored, the wind energy resources must be directly converted into electric energy, so that wind power intermittence and fluctuation and space-time non-adjustability of power generation are brought, in addition, the peak regulation capacity of a power grid in the three-north area with rapid wind power development in China is insufficient, the transmission section of a local area is limited, the active power balance difficulty of a power system is increased, and the phenomenon of wind abandonment in the area is serious.
Based on this, the power grid operating personnel need optimize the load side peak regulation, and reasonable transfer is carried out to the load to improve power grid wind-powered electricity generation reception ability, reduce and abandon wind electric quantity.
Disclosure of Invention
In view of the above, according to the power grid load modeling method for improving the wind power absorption capacity provided by the invention, on the basis of ensuring the calculation efficiency, the abandoned wind power quantity of the wind power is effectively reduced through effective management of the load, the optimization result provides guidance and suggestions for power grid operators, the active power balance difficulty of a power system is effectively reduced, and the efficient and stable operation of the power grid is further ensured.
The purpose of the invention is realized by the following technical scheme:
a method of modeling grid load for production simulation, the method comprising the steps of:
step 1, establishing a load transfer power upper limit model of each time section;
step 2, establishing a load transfer electric quantity upper limit model in each period;
step 3, establishing different load clearing mode models;
step 4, establishing a regional load balance model based on load side peak regulation;
and 5, establishing a load side peak regulation model for improving the wind power absorption capacity, and optimizing the load side peak regulation electric quantity of each time section.
Preferably, the step 1 comprises:
establishing a load transfer power upper limit model of each time section according to the load increasing and decreasing power values of each time section:
Figure GDA0002467689700000021
in the formula (1), the reaction mixture is,
Figure GDA0002467689700000024
increasing a power value for the n load in the time t region;
Figure GDA0002467689700000025
power reduction for n loads in time t regionValue of, and
Figure GDA0002467689700000026
and
Figure GDA0002467689700000027
are all positive variables;
Figure GDA0002467689700000028
the upper limit of transferable power for the region n load at time t.
Preferably, the step 2 comprises:
according to the load transfer power upper limit model of each time section, establishing a load transfer electric quantity upper limit model in each period:
Figure GDA0002467689700000022
in the formula (2), T is a total scheduling period; qnThe total amount of power is transferred for the load.
Preferably, the step 3 comprises:
establishing different load clearing mode models according to the load transfer power upper limit models of the time sections:
Figure GDA0002467689700000023
preferably, the step 4 comprises:
establishing a regional load balance model based on load side peak regulation:
Figure GDA0002467689700000031
in the formula (4), the reaction mixture is,
Figure GDA0002467689700000032
the sum of the total power of all the conventional units at the t moment;
Figure GDA0002467689700000033
is the power load at the t-th moment;
Figure GDA0002467689700000034
for the tie line power values between time t, region n and region nn, and
Figure GDA0002467689700000035
when the value is positive, the current flowing area is in the positive direction;
Figure GDA0002467689700000036
when the value is negative, the current inflow region is in the negative direction;
Figure GDA0002467689700000037
the wind power received by the region n at the time t is used as the wind power.
Preferably, the step 5 comprises:
5-1, establishing a load side peak regulation model for improving the wind power absorption capacity according to the models in the steps 1 to 4;
and 5-2, optimizing the load side peak shaving electric quantity of each time section according to the load side peak shaving model.
Preferably, according to the models in the steps 1 to 4, a constraint condition and an objective function of a load side peak regulation model for improving the wind power absorption capacity are established:
a. and (3) unit optimization power constraint:
Figure GDA0002467689700000038
in the formula (5), the reaction mixture is,
Figure GDA0002467689700000039
a binary variable of the unit j at the time t; pj,max，Pj,minRespectively setting the upper output limit and the lower output limit of the jth unit; pj(t) optimizing power for the unit;
b. minimum on-off time constraint:
Figure GDA00024676897000000310
in the formula (6), the reaction mixture is,
Figure GDA00024676897000000311
respectively representing binary variables of the starting and stopping states of the unit j at the moment t,
Figure GDA00024676897000000312
a "1" indicates that the unit is starting,
Figure GDA0002467689700000041
a "0" indicates that the unit is not in the start-up state,
Figure GDA0002467689700000042
a "1" indicates that the unit is shutting down,
Figure GDA0002467689700000043
a "0" indicates that the unit is not in a shutdown state; k is a radical ofonThe minimum starting time of the unit is set; k is a radical ofoffMinimum down time for the unit; i is a calculation variable;
c. and (3) output constraint of a heat supply unit in a heat supply period:
Figure GDA0002467689700000044
in the formula (7), the reaction mixture is,
Figure GDA0002467689700000045
the output of the back pressure unit is large or small;
Figure GDA0002467689700000046
the output of the air extractor set is large or small;
Figure GDA0002467689700000047
the thermal load at time t;
Figure GDA0002467689700000048
are all the thermoelectric coupling coefficients of the heating unit;
d. start-stop logic state constraint:
Figure GDA0002467689700000049
in the formula (8), the reaction mixture is,
Figure GDA00024676897000000410
the binary variable of the unit j at the time of t-1;
e. and (3) unit climbing rate constraint:
Figure GDA00024676897000000411
in the formula (9), the reaction mixture is,
Figure GDA00024676897000000412
the maximum upward climbing speed and the maximum downward climbing speed of the unit j are respectively set;
Figure GDA00024676897000000413
is the power of the unit j at the time t;
Figure GDA0002467689700000051
the power of the unit j at the moment of t + 1;
f. rotating standby constraint:
Figure GDA0002467689700000052
in the formula (10), J is the total number of the units; pre and Nre are respectively positive rotation standby and negative rotation standby;
Figure GDA0002467689700000053
the electric power at the t-th moment;
g. inter-area line transfer capacity constraints:
Figure GDA0002467689700000054
in the formula (11), the reaction mixture is,
Figure GDA0002467689700000055
the upper limit of the transmission power of the connecting line between the region n and the region nn at the time t,
Figure GDA0002467689700000056
a lower limit of transmission power of a connecting line between the region n and the region nn at the time t;
h. wind power constraint:
Figure GDA0002467689700000057
in the formula (12), Pw,n(t) is wind power;
Figure GDA0002467689700000058
the theoretical output of wind power is obtained;
i. an objective function:
Figure GDA0002467689700000059
in the formula (13), the reaction mixture is,Nis the total number of regions.
According to the technical scheme, the invention provides the power grid load modeling method for production simulation, which comprises the steps of establishing a load transfer power upper limit model of each time section; establishing a load transfer electric quantity upper limit model in each period; establishing different load clearing mode models; establishing a regional load balance model based on load side peak regulation; and establishing a load side peak regulation model for improving the wind power absorption capacity, and optimizing the load side peak regulation electric quantity of each time section. According to the method, on the basis of ensuring the calculation efficiency, the wind curtailment electric quantity of the wind power is effectively reduced through effective management of the load, the optimization result provides guidance and suggestions for power grid operators, the active power balance difficulty of a power system is effectively reduced, and the efficient and stable operation of the power grid is further ensured.
Compared with the closest prior art, the technical scheme provided by the invention has the following excellent effects:
1. according to the technical scheme provided by the invention, the schedulability of the load is optimized and modeled, and the optimization of the load model is realized through load transfer on the premise of ensuring that the load electric quantity is not changed in the scheduling period, so that the wind curtailment electric quantity of the wind power is effectively reduced.
2. According to the technical scheme provided by the invention, the wind power output characteristic, the load time sequence characteristic, the peak regulation characteristic of the unit, the power grid output capacity and other factors of the power grid can be comprehensively considered, and the wind power-containing power balance of the whole power grid can be optimized. The calculation result is more consistent with the actual power system scheduling condition, and the most intuitive judgment basis can be provided for a dispatcher.
3. According to the technical scheme provided by the invention, on the basis of ensuring the calculation efficiency, the wind power abandon electric quantity is effectively reduced through effective management of the load, the optimization result provides guidance and suggestions for power grid operators, the active power balance difficulty of a power system is effectively reduced, and the efficient and stable operation of the power grid is further ensured.
4. The technical scheme provided by the invention has wide application and obvious social benefit and economic benefit.
Drawings
FIG. 1 is a flow chart of a method of grid load modeling for production simulation of the present invention;
FIG. 2 is a flow chart illustrating step 5 of the grid load modeling method of the present invention;
fig. 3 is a schematic flow chart of a specific application example of the power grid load modeling method for production simulation of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments of the present invention without making any creative effort, shall fall within the protection scope of the present invention.
As shown in fig. 1, the present invention provides a power grid load modeling method for production simulation, comprising the following steps:
step 1, establishing a load transfer power upper limit model of each time section;
step 2, establishing a load transfer electric quantity upper limit model in each period;
step 3, establishing different load clearing mode models;
step 4, establishing a regional load balance model based on load side peak regulation;
and 5, establishing a load side peak regulation model for improving the wind power absorption capacity, and optimizing the load side peak regulation electric quantity of each time section.
Wherein, step 1 includes:
establishing a load transfer power upper limit model of each time section according to the load increasing and decreasing power values of each time section:
Figure GDA0002467689700000071
in the formula (1), the reaction mixture is,
Figure GDA0002467689700000072
increasing a power value for the n load in the time t region;
Figure GDA0002467689700000073
for a time t region n load reduction power value, and
Figure GDA0002467689700000074
and
Figure GDA0002467689700000075
are all positive variables;
Figure GDA0002467689700000076
the upper limit of transferable power for the region n load at time t.
Wherein, step 2 includes:
according to the load transfer power upper limit model of each time section, establishing a load transfer electric quantity upper limit model in each period:
Figure GDA0002467689700000081
in the formula (2), T is a total scheduling period; qnThe total amount of power is transferred for the load.
Wherein, step 3 includes:
according to the load transfer power upper limit model of each time section, establishing different load clearing mode models:
Figure GDA0002467689700000082
wherein, step 4 includes:
establishing a regional load balance model based on load side peak regulation:
Figure GDA0002467689700000083
in the formula (4), the reaction mixture is,
Figure GDA0002467689700000084
the sum of the total power of all the conventional units at the t moment;
Figure GDA0002467689700000085
is the power load at the t-th moment;
Figure GDA0002467689700000086
for the tie line power values between time t, region n and region nn, and
Figure GDA0002467689700000087
the value is positive, the current inflow region is positiveDirection;
Figure GDA0002467689700000088
when the value is negative, the current inflow region is in the negative direction;
Figure GDA0002467689700000089
the wind power received by the region n at the time t is used as the wind power.
As shown in fig. 2, step 5 includes:
5-1, establishing a load side peak regulation model for improving the wind power absorption capacity according to the models in the steps 1 to 4;
and 5-2, optimizing the load side peak shaving electric quantity of each time section according to the load side peak shaving model.
Wherein, 5-1 comprises: according to the models in the steps 1 to 4, establishing a constraint condition and an objective function of a load side peak regulation model for improving the wind power absorption capacity:
a. and (3) unit optimization power constraint:
Figure GDA0002467689700000091
in the formula (5), the reaction mixture is,
Figure GDA0002467689700000092
a binary variable of the unit j at the time t; pj,max，Pj,minRespectively setting the upper output limit and the lower output limit of the jth unit; pj(t) optimizing power for the unit;
b. minimum on-off time constraint:
Figure GDA0002467689700000093
in the formula (6), the reaction mixture is,
Figure GDA0002467689700000094
respectively representing binary variables of the starting and stopping states of the unit j at the moment t,
Figure GDA0002467689700000095
a "1" indicates that the unit is starting,
Figure GDA0002467689700000096
a "0" indicates that the unit is not in the start-up state,
Figure GDA0002467689700000097
a "1" indicates that the unit is shutting down,
Figure GDA0002467689700000098
a "0" indicates that the unit is not in a shutdown state; k is a radical ofonThe minimum starting time of the unit is set; k is a radical ofoffMinimum down time for the unit; i is a calculation variable;
c. and (3) output constraint of a heat supply unit in a heat supply period:
Figure GDA0002467689700000099
in the formula (7), the reaction mixture is,
Figure GDA00024676897000000910
the output of the back pressure unit is large or small;
Figure GDA00024676897000000911
the output of the air extractor set is large or small;
Figure GDA00024676897000000912
the thermal load at time t;
Figure GDA00024676897000000913
are all the thermoelectric coupling coefficients of the heating unit;
d. start-stop logic state constraint:
Figure GDA0002467689700000101
in the formula (8), the reaction mixture is,
Figure GDA0002467689700000102
the binary variable of the unit j at the time of t-1;
e. and (3) unit climbing rate constraint:
Figure GDA0002467689700000103
in the formula (9), the reaction mixture is,
Figure GDA0002467689700000104
the maximum upward climbing speed and the maximum downward climbing speed of the unit j are respectively set;
Figure GDA0002467689700000105
is the power of the unit j at the time t;
Figure GDA0002467689700000106
the power of the unit j at the moment of t + 1;
f. rotating standby constraint:
Figure GDA0002467689700000107
in the formula (10), J is the total number of the units; pre and Nre are respectively positive rotation standby and negative rotation standby;
Figure GDA0002467689700000108
the electric power at the t-th moment;
g. inter-area line transfer capacity constraints:
Figure GDA0002467689700000109
in the formula (11), the reaction mixture is,
Figure GDA0002467689700000111
the upper limit of the transmission power of the connecting line between the region n and the region nn at the time t,
Figure GDA0002467689700000112
a lower limit of transmission power of a connecting line between the region n and the region nn at the time t;
h. wind power constraint:
Figure GDA0002467689700000113
in the formula (12), Pw,n(t) is wind power;
Figure GDA0002467689700000114
the theoretical output of wind power is obtained;
i. an objective function:
Figure GDA0002467689700000115
in the formula (13), the reaction mixture is,Nis the total number of regions.
As shown in fig. 3, the present invention provides a specific application example of a power grid load modeling method for production simulation, which is as follows:
in the first step, the upper limit of the load transfer power of each time section is modeled.
(1) Load transfer power cap modeling
Figure GDA0002467689700000116
In the formula (I), the compound is shown in the specification,
Figure GDA0002467689700000117
and
Figure GDA0002467689700000118
respectively representing the n load increasing power of the t time zone and the n load decreasing power of the t time zone, which are positive variables,
Figure GDA0002467689700000119
transferable Power for n loads in time t regionAn upper limit. This constraint limits the upper limit of the load transfer power at time t.
And secondly, modeling the upper limit of the load transfer electric quantity in each period.
(2) Load transfer capacity upper limit modeling
Figure GDA00024676897000001110
In the formula, T is a scheduling period of the load model, taking a simulation time step of 1 hour as an example, if the simulation time step is a load transfer daily electric quantity constraint, then T is 24, if the simulation time step is a load transfer weekly electric quantity constraint, then T is 168, if the simulation time step is a load transfer monthly electric quantity constraint (a month is calculated by 30 days), then T is 720. Q is the total electric quantity of load transfer, and the value can be determined according to the value situation of T. This constraint limits the upper limit of the load transfer capacity over the entire scheduling period.
And thirdly, modeling different load clearing modes.
(3) Modeling in different load clearing modes
Figure GDA0002467689700000121
The constraint means that the sum of the power increasing upwards and the sum of the power decreasing downwards of all time section loads are the same in the total scheduling period of T, namely the power consumption of the loads in the scheduling period is ensured to be unchanged. The value of T can be selected according to the load management mode adopted by the power grid dispatching operator, if the power grid dispatching operator performs load management by adopting a load day clearing mode, then T is 24, if the power grid dispatching operator performs load management by adopting a week clearing mode, then T is 168, and if the power grid dispatching operator performs load management by adopting a month clearing mode (30 days in a month), then T is 720.
And fourthly, establishing a regional load balance model based on the load model.
(4) Regional load balancing constraints
Figure GDA0002467689700000122
In the formula (I), the compound is shown in the specification,
Figure GDA0002467689700000123
the sum of the total power of all conventional units in the t-th period;
Figure GDA0002467689700000124
it represents the electric load of the t-th period. In the formula
Figure GDA0002467689700000125
The magnitude of the tie line power between the region n and the region nn at the time t. Setting the current reference direction as follows: the inflow region is in the positive direction and the outflow region is in the negative direction. Therefore, it is not only easy to use
Figure GDA0002467689700000126
The values can be positive and negative, which represent the direction of power transfer.
Figure GDA0002467689700000127
The wind power received for the t period of the n region.
And fifthly, comprehensively considering factors such as wind power output time sequence characteristics, load time sequence characteristics, peak regulation characteristics of the unit, power grid output capacity and the like of the power grid, establishing a provincial power grid time sequence simulation model, and optimizing the wind power-containing power balance of the whole power grid. The provincial power grid time sequence simulation model is consistent with a method for making a wind power annual plan based on time sequence simulation (volume 38, page 11 and page 13 of power system automation), which is described briefly herein.
(5) Optimized power constraint of unit
Figure GDA0002467689700000131
(6) Minimum on-off time constraint
Figure GDA0002467689700000132
Figure GDA0002467689700000133
(7) Output constraint of heat supply unit in heat supply period
Figure GDA0002467689700000134
Figure GDA0002467689700000135
(8) Start-stop logic state constraints
Figure GDA0002467689700000136
(9) Unit ramp rate constraint
Figure GDA0002467689700000137
Figure GDA0002467689700000138
(10) Rotational back-up restraint
Figure GDA0002467689700000139
(11) Inter-area line transfer capacity constraints
Figure GDA0002467689700000141
(12) Wind power constraint
Figure GDA0002467689700000142
(13) Objective function
Figure GDA0002467689700000143
In the formula, Pj,max，Pj,minRespectively is the upper limit and the lower limit of the output of the jth unit.
Figure GDA0002467689700000144
Respectively representing binary variables of the starting and stopping states of the unit j in the time period t,
Figure GDA0002467689700000145
a "1" indicates that the unit is starting, a "0" indicates that the unit is not in a starting state,
Figure GDA0002467689700000146
a value of "1" indicates that the unit is in a shutdown state, and a value of "0" indicates that the unit is not in a shutdown state; k is a radical ofonThe minimum starting time of the unit is set; k is a radical ofoffMinimum down time for the unit; it reflects the minimum startup or shutdown time length, with different types of units having different startup and shutdown time parameters.
Figure GDA0002467689700000147
The output of the back pressure unit is large or small;
Figure GDA0002467689700000148
the output of the air extractor set is large or small;
Figure GDA0002467689700000149
the thermal load is the t time period;
Figure GDA00024676897000001410
the thermoelectric coupling coefficient of the heating unit is shown.
Figure GDA00024676897000001411
Are respectively provided withThe maximum upward climbing rate and the maximum downward climbing rate of the unit j. Pre and Nre are positive spinning standby and negative spinning standby, respectively.
Figure GDA00024676897000001412
The upper limit of the transmission power of the connecting line between the region n and the region nn at the time t,
Figure GDA00024676897000001413
the lower limit of the transmission power of the tie line between the region n and the region nn at the time t.
Figure GDA00024676897000001414
The theoretical output of wind power is obtained.
And sixthly, optimizing the load electric quantity of each time section by adopting the mathematical model established by the new method, and optimizing the starting mode of the unit by transferring the load power, so that the abandoned wind electric quantity of the wind power is greatly reduced, and the optimization result can provide guidance and suggestion for power grid dispatching personnel.
Although the present invention has been described in detail with reference to the above embodiments, those skilled in the art can make modifications and equivalents to the embodiments of the present invention without departing from the spirit and scope of the present invention, which is set forth in the claims of the present application.

Claims (7)

1. A method for modeling grid load for production simulation, the method comprising the steps of:
step 1, establishing a load transfer power upper limit model of each time section;
step 2, establishing a load transfer electric quantity upper limit model in each period;
step 3, establishing different load clearing mode models;
step 4, establishing a regional load balance model based on load side peak regulation;
and 5, establishing a load side peak regulation model for improving the wind power absorption capacity, and optimizing the load side peak regulation electric quantity of each time section.
2. The method of claim 1, wherein step 1 comprises:
establishing a load transfer power upper limit model of each time section according to the load increasing and decreasing power values of each time section:
Figure FDA0000893561720000011
in the formula (1), the reaction mixture is,
Figure FDA0000893561720000012
increasing a power value for the n load in the time t region;
Figure FDA0000893561720000013
for a time t region n load reduction power value, and
Figure FDA0000893561720000014
and
Figure FDA0000893561720000015
are all positive variables;
Figure FDA0000893561720000016
the upper limit of transferable power for the region n load at time t.
3. The method of claim 2, wherein step 2 comprises:
according to the load transfer power upper limit model of each time section, establishing a load transfer electric quantity upper limit model in each period:
Figure FDA0000893561720000017
in the formula (2), T is a total scheduling period; qnThe total amount of power is transferred for the load.
4. The method of claim 3, wherein step 3 comprises:
establishing different load clearing mode models according to the load transfer power upper limit models of the time sections:
Figure FDA0000893561720000021
5. the method of claim 4, wherein step 4 comprises:
establishing a regional load balance model based on load side peak regulation:
Figure FDA0000893561720000022
in the formula (4), the reaction mixture is,
Figure FDA0000893561720000023
the sum of the total power of all the conventional units at the t moment;
Figure FDA0000893561720000024
is the power load at the t-th moment;
Figure FDA0000893561720000025
for the tie line power values between time t, region n and region nn, and
Figure FDA0000893561720000026
when the value is positive, the current flowing area is in the positive direction;
Figure FDA0000893561720000027
when the value is negative, the current inflow region is in the negative direction;
Figure FDA0000893561720000028
the wind power received by the region n at the time t is used as the wind power.
6. The method of claim 5, wherein the step 5 comprises:
5-1, establishing a load side peak regulation model for improving the wind power absorption capacity according to the models in the steps 1 to 4;
and 5-2, optimizing the load side peak shaving electric quantity of each time section according to the load side peak shaving model.
7. The method of claim 6, wherein the 5-1 comprises: according to the models in the steps 1 to 4, establishing a constraint condition and an objective function of a load side peak regulation model for improving the wind power absorption capacity:
a. and (3) unit optimization power constraint:
Figure FDA0000893561720000029
in the formula (5), the reaction mixture is,
Figure FDA00008935617200000210
a binary variable of the unit j at the time t; pj,max，Pj,minRespectively setting the upper output limit and the lower output limit of the jth unit; pj(t) optimizing power for the unit;
b. minimum on-off time constraint:
Figure FDA0000893561720000031
in the formula (6), the reaction mixture is,
Figure FDA0000893561720000032
respectively representing binary variables of the starting and stopping states of the unit j at the moment t,
Figure FDA0000893561720000033
a "1" indicates that the unit is starting,
Figure FDA0000893561720000034
a "0" indicates that the unit is not in the start-up state,
Figure FDA0000893561720000035
a "1" indicates that the unit is shutting down,
Figure FDA0000893561720000036
a "0" indicates that the unit is not in a shutdown state; k is a radical ofonThe minimum starting time of the unit is set; k is a radical ofoffMinimum down time for the unit; i is a calculation variable;
c. and (3) output constraint of a heat supply unit in a heat supply period:
Figure FDA0000893561720000037
in the formula (7), the reaction mixture is,
Figure FDA0000893561720000038
the output of the back pressure unit is large or small;
Figure FDA0000893561720000039
the output of the air extractor set is large or small;
Figure FDA00008935617200000310
the thermal load at time t;
Figure FDA00008935617200000311
are all the thermoelectric coupling coefficients of the heating unit;
d. start-stop logic state constraint:
Figure FDA00008935617200000312
in the formula (8), the reaction mixture is,
Figure FDA00008935617200000313
the binary variable of the unit j at the time of t-1;
e. and (3) unit climbing rate constraint:
Figure FDA0000893561720000041
in the formula (9), the reaction mixture is,
Figure FDA0000893561720000042
the maximum upward climbing speed and the maximum downward climbing speed of the unit j are respectively set;
Figure FDA0000893561720000043
is the power of the unit j at the time t;
Figure FDA0000893561720000044
the power of the unit j at the moment of t + 1;
f. rotating standby constraint:
Figure FDA0000893561720000045
in the formula (10), J is the total number of the units; pre and Nre are respectively positive rotation standby and negative rotation standby;
Figure FDA0000893561720000046
the electric power at the t-th moment;
g. inter-area line transfer capacity constraints:
Figure FDA0000893561720000047
in the formula (11), the reaction mixture is,
Figure FDA0000893561720000048
the upper limit of the transmission power of the connecting line between the region n and the region nn at the time t,
Figure FDA0000893561720000049
a lower limit of transmission power of a connecting line between the region n and the region nn at the time t;
h. wind power constraint:
Figure FDA00008935617200000410
in the formula (12), Pw,n(t) is wind power;
Figure FDA00008935617200000411
the theoretical output of wind power is obtained;
i. an objective function:
Figure FDA0000893561720000051
in formula (13), N is the total number of regions.
CN201511004804.XA 2015-12-29 2015-12-29 Power grid load modeling method for production simulation Active CN105470957B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201511004804.XA CN105470957B (en) 2015-12-29 2015-12-29 Power grid load modeling method for production simulation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201511004804.XA CN105470957B (en) 2015-12-29 2015-12-29 Power grid load modeling method for production simulation

Publications (2)

Publication Number Publication Date
CN105470957A CN105470957A (en) 2016-04-06
CN105470957B true CN105470957B (en) 2021-01-15

Family

ID=55608402

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201511004804.XA Active CN105470957B (en) 2015-12-29 2015-12-29 Power grid load modeling method for production simulation

Country Status (1)

Country Link
CN (1) CN105470957B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105932711B (en) * 2016-04-26 2018-06-29 国网电力科学研究院 The assessment measuring method of wind is abandoned in a kind of wind-power electricity generation-load real time coordination consumption
CN107681673B (en) * 2017-09-18 2019-09-03 中国电力工程顾问集团西北电力设计院有限公司 The selection method of the tower photo-thermal Generating Set Peak shaving Mode of fused salt based on power producing characteristics
CN109038589B (en) * 2018-08-02 2020-05-05 中国电力工程顾问集团西北电力设计院有限公司 Multi-provincial power grid coordinated operation production simulation method
CN111489009B (en) * 2019-06-06 2023-07-25 国网辽宁省电力有限公司 Optimization calculation method and device for operation mode of electric vehicle charging station

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101789598A (en) * 2010-03-05 2010-07-28 湖北省电力试验研究院 Power system load modelling method
CN103296679A (en) * 2013-05-20 2013-09-11 国家电网公司 Modeling method for medium and long-term wind power output model of power system capable of optimally running for medium and long terms
CN103762589A (en) * 2014-01-08 2014-04-30 河海大学 Method for optimizing new energy capacity ratio in layers in power grid
CN104092213A (en) * 2014-07-30 2014-10-08 东南大学 Power analyzing method for indeterminate power flow branches based on optimization method
CN104993523A (en) * 2015-07-16 2015-10-21 国家电网公司 Pumped storage power station characteristic accurate simulation method for optimized operation of wind power contained power grid system

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5557801B2 (en) * 2010-06-15 2014-07-23 株式会社日立製作所 System state calculation device, system control device, system state calculation system, distribution system power flow simulation device, system state calculation method, system control method, power distribution system power flow simulation method, and program thereof
ES2535059B1 (en) * 2013-10-31 2016-02-09 Control Techniques Iberia S.A. Method and system to control a power supply to a load

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101789598A (en) * 2010-03-05 2010-07-28 湖北省电力试验研究院 Power system load modelling method
CN103296679A (en) * 2013-05-20 2013-09-11 国家电网公司 Modeling method for medium and long-term wind power output model of power system capable of optimally running for medium and long terms
CN103762589A (en) * 2014-01-08 2014-04-30 河海大学 Method for optimizing new energy capacity ratio in layers in power grid
CN104092213A (en) * 2014-07-30 2014-10-08 东南大学 Power analyzing method for indeterminate power flow branches based on optimization method
CN104993523A (en) * 2015-07-16 2015-10-21 国家电网公司 Pumped storage power station characteristic accurate simulation method for optimized operation of wind power contained power grid system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
A Novel Annual Wind Power Planning Method of Provincial Power Grid Based on Time Sequential Simulations;Cao Yang et al.;《2014 International Conference on Power System Technology (POWERCON 2014)》;20141022;第2693-2699页 *
基于时序仿真的风电年度计划制定方法;刘纯 等;《电力系统自动化》;20140610;第38卷(第11期);第13-19页 *

Also Published As

Publication number Publication date
CN105470957A (en) 2016-04-06

Similar Documents

Publication Publication Date Title
WO2022100091A1 (en) Centralized control method for scheduling of generalized source storage system
CN108154309B (en) Energy internet economic dispatching method considering multi-load dynamic response of cold, heat and electricity
CN103066626B (en) Multi-source coordinating day-ahead generation scheduling method
CN106300336A (en) A kind of meter and the virtual plant Multiobjective Optimal Operation method of load side and mains side
CN103728881B (en) A kind of optimizing operation method of many edifice control system system
CN104795846A (en) Optimized operation method of pumped-storage power station and wind power combined system
WO2019128012A1 (en) Robust optimal coordinated dispatching method for alternating-current and direct-current hybrid micro-grid
CN105337415A (en) Regional power grid dispatching system and method based on prediction control
CN105470957B (en) Power grid load modeling method for production simulation
CN107276122B (en) Peak-shaving resource calling decision method suitable for large-scale renewable energy grid connection
CN103296682A (en) Multiple spatial and temporal scale gradually-advancing load dispatching mode designing method
CN108628176B (en) Thermoelectric combined optimization scheduling method considering heat supply network flexibility
CN112039056A (en) Two-stage optimal scheduling method for new energy
CN104779611A (en) Economic dispatch method for micro grid based on centralized and distributed double-layer optimization strategy
CN108629449B (en) Distributed robust type optimal scheduling method for alternating current-direct current hybrid microgrid
CN104882905A (en) New energy available capability assessment method considering transient security constraints
CN104361416A (en) Power-grid double-layer optimized dispatching method considering large-scale electric automobile access
CN103219751A (en) Control method of active power of clustered wind power plants
CN109993345B (en) Garden-oriented dynamic economic dispatching method for multi-energy complementary system for island operation
CN111049192A (en) Power generation control method considering renewable energy resource bidding on same station
CN104124713B (en) A kind of Wind turbines optimizes progress control method
CN105186583B (en) Energy router and its energy dispatching method based on multi-agent modeling
CN110854933A (en) Monthly unit combination optimization method utilizing flexible resources
CN104239960A (en) Electricity generating schedule optimizing method considering pump storage unit
CN106712105A (en) New energy acceptance space discrete probability sequence calculation method

Legal Events

Date Code Title Description
C06 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