CN105470957A - Power network load modeling method for production analogue simulation - Google Patents

Abstract

The invention provides a power network load modeling method for production analogue simulation. The method comprises the following steps: building a load transfer power upper-limit models for various time sections; building a load transfer electricity quantity upper-limit model within various periods; building different load settlement mode models; building a regional load balance model on the basis of load side peak regulation; building a load side peak regulation model for improving the wind power absorptive capability; and optimizing the load side peak regulation electricity quantity of various time sections. According to the method provided by the invention, the abandoned electricity quantity of wind power is effectively reduced through effective management of loads on the basis of ensuring the computational efficiency; the optimization result provides guidance and advice for power network dispatching personnel; and the active balance difficulty of a power system is effectively lowered, so that the condition that the power network efficiently and stably runs is ensured.

Description

A kind of network load modeling method for the production of analog simulation

Technical field

The present invention relates to field of new energy generation, be specifically related to a kind of network load modeling method for the production of analog simulation.

Background technology

Along with China's sustained and rapid development of economy, fossil fuel peter out and it is to the continuous aggravation of the pollution of environment and greenhouse effect, based on the strategy of sustainable development, the Chinese government pays much attention to develop new forms of energy, new energy development is utilized as control atmosphere pollution, readjusted the energy structure and the major action of Transformation of Economic Growth, and using wind power generation and photovoltaic generation one of the major way as new energy development and utilization.

Due to the thin property of wind energy resources, space time energy metric density is low, cannot enrichment, transport and storage, directly must be converted to electric energy, bring intermittence and the fluctuation of wind-powered electricity generation thus, and the space-time non-scheduling of generating, " three Norths " area power grid peak modulation capacity that China's Wind Power Development is swift and violent is in addition not enough and the restriction of local area transmission section, increase electric power system active balance difficulty, it is serious that wind phenomenon is abandoned in such area.

Based on this, operation of power networks personnel need be optimized load side peak regulation, by reasonably shifting load, thus improving power grid wind receiving ability, reducing and abandoning wind-powered electricity generation amount.

Summary of the invention

In view of this, a kind of network load modeling method improving wind electricity digestion capability provided by the invention, the method is on the basis ensureing computational efficiency, by the effective management to load, effective minimizing wind-powered electricity generation abandon wind-powered electricity generation amount, optimum results for operation of power networks personnel provide instruct and suggestion, effectively reduce the active balance difficulty of electric power system, so ensure that electrical network efficiently and stably run.

The object of the invention is to be achieved through the following technical solutions:

For the production of a network load modeling method for analog simulation, described method comprises the steps:

The load transfer plan power upper limit model of discontinuity surface when step 1. is set up each;

Step 2. sets up the load transfer plan electricity Upper-Bound Model in each cycle;

Step 3. sets up different load cleancut models;

Step 4. sets up the region load balancing model based on load side peak regulation;

Step 5. sets up the load side peak regulation model improving wind electricity digestion capability, the load side peak regulation power of discontinuity surface when optimizing each.

Preferably, described step 1 comprises:

According to the performance number that the load of discontinuity surface time each increases and declines, the load transfer plan power upper limit model of discontinuity surface when setting up each:

$\left\{\begin{array}{c}{\mathrm{\Δl}}_{n,u}^t\≤{\mathrm{\Δl}}_{n,max}^t\\ {\mathrm{\Δl}}_{n,d}^t\≤{\mathrm{\Δl}}_{n,\max }^t\end{array}\right.---\left(1\right)$

In formula (1), for t region n load increases performance number; for t region n load decline performance number, and with be positive variable; for the transferable power upper limit of t region n load.

Preferably, described step 2 comprises:

According to the described load transfer plan power upper limit model of discontinuity surface time each, set up the load transfer plan electricity Upper-Bound Model in each cycle:

$\underset{t=1}{\overset{T}{\Σ}}({\mathrm{\Δl}}_{n,u}^t+{\mathrm{\Δl}}_{n,d}^t)/2\≤Q_n---\left(2\right)$

In formula (2), T is the total activation cycle; Q _nfor the total electricity of load transfer plan.

Preferably, described step 3 comprises:

According to the described load transfer plan power upper limit model of discontinuity surface time each, set up different load cleancut models:

$\underset{t=1}{\overset{T}{\Σ}}{\mathrm{\Δl}}_{n,u}^t-\underset{t=1}{\overset{T}{\Σ}}{\mathrm{\Δl}}_{n,d}^t=0---\left(3\right).$

Preferably, described step 4 comprises:

Set up the region load balancing model based on load side peak regulation:

$P_{all.n}^t+P_w^{t,n}+\underset{nn=1}{\overset{N}{\Σ}}{L}_{n,nn}^t=P_{l.n}^t+{\mathrm{\Δl}}_{n,u}^t-{\mathrm{\Δl}}_{n,d}^t---\left(4\right)$

In formula (4), it is the gross power sum of all conventional power unit of t; it is the electric load of t; for t, dominant eigenvalues value between region n and region nn, and value is timing, and electric current inflow region is positive direction; value is for time negative, and electric current inflow region is negative direction; be the wind power generation power of t, region n receiving.

Preferably, described step 5 comprises:

5-1., according to the model in step 1 to 4, sets up the load side peak regulation model improving wind electricity digestion capability;

5-2. according to described load side peak regulation model, the load side peak regulation power of discontinuity surface when optimizing each.

Preferably, according to the model in step 1 to 4, set up constraints and the target function of the load side peak regulation model improving wind electricity digestion capability:

A. set optimization power constraint:

$X_j^t\·P_{j,min}\≤P_j\left(t\right)\≤P_{j,max}\·X_j^t---\left(5\right)$

In formula (5), for unit j is at the binary variable of t; P _{j, max}, P _{j, min}be respectively the exert oneself upper limit and the lower limit of exerting oneself of jth platform unit; P _jt () is set optimization power;

B. minimum start and stop time-constrain:

$\left\{\begin{array}{c}{Y}_j^t+\underset{i=1}{\overset{k_{on}}{\Σ}}{Z}_j^{t+i}\≤1\\ Z_j^t+\underset{i=1}{\overset{k_{off}}{\Σ}}{Y}_j^{t+i}\≤1\end{array}\right.---\left(6\right)$

In formula (6), be respectively and represent that unit j starts in t, the binary variable of stopped status, for " 1 " represents that unit starts, for " 0 " represents that unit is not at starting state, for " 1 " represents that unit is shut down, for " 0 " represents that unit is not in stopped status; k _onthe machine time is opened for unit is minimum; k _offfor unit minimum downtime; I is for calculating variable;

C. heat supply phase thermal power plant unit units limits:

$\left\{\begin{array}{c}{P}_{j,BY}^t=C_j^b\·H_j^t\\ H_j^t\·C_j^b\≤P_{j,CQ}^t\≤P_{j,\max }-H_j^t\·C_j^v\end{array}\right.---\left(7\right)$

In formula (7), for back pressure unit is exerted oneself size; for unit output size of bleeding; for t load of heat; be thermal power plant unit coupled thermomechanics coefficient;

D. start and stop logic state constraint:

$\left\{\begin{array}{c}{X}_j^t-X_j^{t-1}-Y_j^t+Z_j^t=0\\ -X_j^t-X_j^{t-1}+Y_j^t\≤0\\ X_j^t+X_j^{t-1}+Y_j^t\≤2\\ -X_j^t-X_j^{t-1}+Z_j^t\≤0\\ X_j^t+X_j^{t-1}+Z_j^t\≤2\end{array}\right.---\left(8\right)$

In formula (8), for unit j is at the binary variable in t-1 moment;

E. unit climbing rate constraint:

$\left\{\begin{array}{c}{P}_j^{t+1}-P_j^t\≤{\mathrm{\ΔP}}_{j,up}\\ P_j^t-P_j^{t+1}\≤{\mathrm{\ΔP}}_{j,down}\end{array}\right.---\left(9\right)$

In formula (9), be respectively the maximum upper creep speed of unit j and lower creep speed; for t unit j power; for t+1 moment unit j power;

F. spinning reserve constraint:

$\left\{\begin{array}{c}-\underset{j=1}{\overset{J}{\Σ}}{P}_{j,max}\·X_j^t\≤-P_l^t-\mathrm{Pr}e\\ \underset{j=1}{\overset{J}{\Σ}}{P}_{j,min}\·X_j^t\≤P_l^t-Nre\end{array}\right.---\left(10\right)$

In formula (10), J is unit sum; Pre and Nre is respectively positive rotation spinning reserve for subsequent use and negative; it is the electric power of t;

G. interregional line transmission capacity-constrained:

$L_{n,nn}^{\min }\≤L_{n,nn}^t\≤L_{n,nn}^{\max }---\left(11\right)$

In formula (11), for the tie-line power transmission upper limit between t region n and region nn, for tie-line power transmission lower limit between t region n and region nn;

H. wind power constraint:

$0\≤P_{w,n}\left(t\right)\≤P_{w,n}^{*}\left(t\right)---\left(12\right)$

In formula (12), P _w,nt () is wind power; for wind-powered electricity generation theory is exerted oneself;

I. target function:

$max\underset{t=1}{\overset{T}{\Σ}}\underset{n=1}{\overset{N}{\Σ}}{P}_{w,n}\left(t\right)---\left(13\right)$

In formula (13), N is the sum in region.

As can be seen from above-mentioned technical scheme, the invention provides a kind of network load modeling method for the production of analog simulation, the load transfer plan power upper limit model of discontinuity surface during by setting up each; Set up the load transfer plan electricity Upper-Bound Model in each cycle; Set up different load cleancut models; Set up the region load balancing model based on load side peak regulation; Set up the load side peak regulation model improving wind electricity digestion capability, the load side peak regulation power of discontinuity surface when optimizing each.The method that the present invention proposes is on the basis ensureing computational efficiency, by the effective management to load, effective minimizing wind-powered electricity generation abandon wind-powered electricity generation amount, optimum results instructs and suggestion for operation of power networks personnel provide, effectively reduce the active balance difficulty of electric power system, so ensure that electrical network efficiently and stably run.

With immediate prior art ratio, technical scheme provided by the invention has following excellent effect:

1, in technical scheme provided by the present invention, be optimized modeling to the schedulability of load, under the prerequisite that guarantee internal loading electricity dispatching cycle is constant, by load transfer plan, realize the optimization to load model, that effectively reduces wind-powered electricity generation abandons wind-powered electricity generation amount.

2, technical scheme provided by the present invention, can consider the wind power output characteristic of electrical network, load temporal characteristics, peak load regulation characteristic, electrical network send the factors such as ability, optimize the power balance of the whole network containing wind-powered electricity generation.The more realistic power system dispatching situation of result of calculation, can provide basis for estimation the most intuitively for dispatcher.

3, technical scheme provided by the present invention, on the basis ensureing computational efficiency, by the effective management to load, effective minimizing wind-powered electricity generation abandon wind-powered electricity generation amount, optimum results instructs and suggestion for operation of power networks personnel provide, effectively reduce the active balance difficulty of electric power system, so ensure that electrical network efficiently and stably run.

4, technical scheme provided by the invention, is widely used, and has significant Social benefit and economic benefit.

Accompanying drawing explanation

Fig. 1 is the flow chart of a kind of network load modeling method for the production of analog simulation of the present invention;

Fig. 2 is the schematic flow sheet of the step 5 of network load modeling method of the present invention;

Fig. 3 is the schematic flow sheet of the embody rule example of a kind of network load modeling method for the production of analog simulation of the present invention.

Embodiment

Below in conjunction with the accompanying drawing in the embodiment of the present invention, be clearly and completely described the technical scheme in the embodiment of the present invention, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiments.Based on embodiments of the invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong to the scope of protection of the invention.

As shown in Figure 1, the invention provides a kind of network load modeling method for the production of analog simulation, comprise the steps:

The load transfer plan power upper limit model of discontinuity surface when step 1. is set up each;

Step 2. sets up the load transfer plan electricity Upper-Bound Model in each cycle;

Step 3. sets up different load cleancut models;

Step 4. sets up the region load balancing model based on load side peak regulation;

Step 5. sets up the load side peak regulation model improving wind electricity digestion capability, the load side peak regulation power of discontinuity surface when optimizing each.

Wherein, step 1 comprises:

According to the performance number that the load of discontinuity surface time each increases and declines, the load transfer plan power upper limit model of discontinuity surface when setting up each:

$\left\{\begin{array}{c}{\mathrm{\Δl}}_{n,u}^t\≤{\mathrm{\Δl}}_{n,max}^t\\ {\mathrm{\Δl}}_{n,d}^t\≤{\mathrm{\Δl}}_{n,\max }^t\end{array}\right.---\left(1\right)$

In formula (1), for t region n load increases performance number; for t region n load decline performance number, and with be positive variable; for the transferable power upper limit of t region n load.

Wherein, step 2 comprises:

According to the load transfer plan power upper limit model of discontinuity surface time each, set up the load transfer plan electricity Upper-Bound Model in each cycle:

$\underset{t=1}{\overset{T}{\Σ}}({\mathrm{\Δl}}_{n,u}^t+{\mathrm{\Δl}}_{n,d}^t)/2\≤Q_n---\left(2\right)$

In formula (2), T is the total activation cycle; Q _nfor the total electricity of load transfer plan.

Wherein, step 3 comprises:

According to the load transfer plan power upper limit model of discontinuity surface time each, set up different load cleancut models:

$\underset{t=1}{\overset{T}{\Σ}}{\mathrm{\Δl}}_{n,u}^t-\underset{t=1}{\overset{T}{\Σ}}{\mathrm{\Δl}}_{n,d}^t=0---\left(3\right).$

Wherein, step 4 comprises:

Set up the region load balancing model based on load side peak regulation:

$P_{all.n}^t+P_w^{t,n}+\underset{nn=1}{\overset{N}{\Σ}}{L}_{n,nn}^t=P_{l.n}^t+{\mathrm{\Δl}}_{n,u}^t-{\mathrm{\Δl}}_{n,d}^t---\left(4\right)$

In formula (4), it is the gross power sum of all conventional power unit of t; it is the electric load of t; for t, dominant eigenvalues value between region n and region nn, and value is timing, and electric current inflow region is positive direction; value is for time negative, and electric current inflow region is negative direction; be the wind power generation power of t, region n receiving.

As shown in Figure 2, step 5 comprises:

5-1., according to the model in step 1 to 4, sets up the load side peak regulation model improving wind electricity digestion capability;

5-2. according to load side peak regulation model, the load side peak regulation power of discontinuity surface when optimizing each.

Wherein, 5-1 comprises: according to the model in step 1 to 4, sets up constraints and the target function of the load side peak regulation model improving wind electricity digestion capability:

A. set optimization power constraint:

$X_j^t\·P_{j,min}\≤P_j\left(t\right)\≤P_{j,max}\·X_j^t---\left(5\right)$

In formula (5), for unit j is at the binary variable of t; P _{j, max}, P _{j, min}be respectively the exert oneself upper limit and the lower limit of exerting oneself of jth platform unit; P _jt () is set optimization power;

B. minimum start and stop time-constrain:

$\left\{\begin{array}{c}{Y}_j^t+\underset{i=1}{\overset{k_{on}}{\Σ}}{Z}_j^{t+i}\≤1\\ Z_j^t+\underset{i=1}{\overset{k_{off}}{\Σ}}{Y}_j^{t+i}\≤1\end{array}\right.---\left(6\right)$

In formula (6), be respectively and represent that unit j starts in t, the binary variable of stopped status, for " 1 " represents that unit starts, for " 0 " represents that unit is not at starting state, for " 1 " represents that unit is shut down, for " 0 " represents that unit is not in stopped status; k _onthe machine time is opened for unit is minimum; k _offfor unit minimum downtime; I is for calculating variable;

C. heat supply phase thermal power plant unit units limits:

$\left\{\begin{array}{c}{P}_{j,BY}^t=C_j^b\·H_j^t\\ H_j^t\·C_j^b\≤P_{j,CQ}^t\≤P_{j,\max }-H_j^t\·C_j^v\end{array}\right.---\left(7\right)$

In formula (7), for back pressure unit is exerted oneself size; for unit output size of bleeding; for t load of heat; be thermal power plant unit coupled thermomechanics coefficient;

D. start and stop logic state constraint:

$\left\{\begin{array}{c}{X}_j^t-X_j^{t-1}-Y_j^t+Z_j^t=0\\ -X_j^t-X_j^{t-1}+Y_j^t\≤0\\ X_j^t+X_j^{t-1}+Y_j^t\≤2\\ -X_j^t-X_j^{t-1}+Z_j^t\≤0\\ X_j^t+X_j^{t-1}+Z_j^t\≤2\end{array}\right.---\left(8\right)$

In formula (8), for unit j is at the binary variable in t-1 moment;

E. unit climbing rate constraint:

$\left\{\begin{array}{c}{P}_j^{t+1}-P_j^t\≤{\mathrm{\ΔP}}_{j,up}\\ P_j^t-P_j^{t+1}\≤{\mathrm{\ΔP}}_{j,down}\end{array}\right.---\left(9\right)$

In formula (9), be respectively the maximum upper creep speed of unit j and lower creep speed; for t unit j power; for t+1 moment unit j power;

F. spinning reserve constraint:

$\left\{\begin{array}{c}-\underset{j=1}{\overset{J}{\Σ}}{P}_{j,max}\·X_j^t\≤-P_l^t-\mathrm{Pr}e\\ \underset{j=1}{\overset{J}{\Σ}}{P}_{j,min}\·X_j^t\≤P_l^t-Nre\end{array}\right.---\left(10\right)$

In formula (10), J is unit sum; Pre and Nre is respectively positive rotation spinning reserve for subsequent use and negative; it is the electric power of t;

G. interregional line transmission capacity-constrained:

$L_{n,nn}^{\min }\≤L_{n,nn}^t\≤L_{n,nn}^{\max }---\left(11\right)$

In formula (11), for the tie-line power transmission upper limit between t region n and region nn, for tie-line power transmission lower limit between t region n and region nn;

H. wind power constraint:

$0\≤P_{w,n}\left(t\right)\≤P_{w,n}^{*}\left(t\right)---\left(12\right)$

In formula (12), P _w,nt () is wind power; for wind-powered electricity generation theory is exerted oneself;

I. target function:

$max\underset{t=1}{\overset{T}{\Σ}}\underset{n=1}{\overset{N}{\Σ}}{P}_{w,n}\left(t\right)---\left(13\right)$

In formula (13), _nfor the sum in region.

As shown in Figure 3, the invention provides a kind of embody rule example of the network load modeling method for the production of analog simulation, as follows:

The first step, carries out modeling to the load transfer plan power upper limit of discontinuity surface time each.

(1) load transfer plan power upper limit modeling

$\begin{array}{c}{\mathrm{\Δl}}_{n,u}^t\≤{\mathrm{\Δl}}_{n,max}^t\\ {\mathrm{\Δl}}_{n,d}^t\≤{\mathrm{\Δl}}_{n,\max }^t\end{array}---\left(1\right)$

In formula, with represent t region n load respectively and increase watt level, and t region n load decline watt level, be positive variable, for the transferable power upper limit of t region n load.This constrained upper limit of t load transfer plan power.

Second step, carries out modeling to the load transfer plan electricity upper limit in each cycle.

(2) load transfer plan electricity upper limit modeling

$\underset{t=1}{\overset{T}{\Σ}}({\mathrm{\Δl}}_{n,u}^t+{\mathrm{\Δl}}_{n,d}^t)/2\≤Q_n---\left(2\right)$

In formula, T is the dispatching cycle of load model, for simulation time step-length 1 hour, as being load transfer plan day Constraint, then now T=24, as being load transfer plan week Constraint, then now T=168, as load transfer plan moon Constraint (January by 30 days calculate), then now T=720.Q is the total electricity of load transfer plan, and its value can be determined according to the value condition of T.In the whole scheduling slot of this constrained, the upper limit of load transfer plan electricity.

3rd step, to different load cleancut modelings.

(3) different load cleancut modeling

$\underset{t=1}{\overset{T}{\Σ}}{\mathrm{\Δl}}_{n,u}^t-\underset{t=1}{\overset{T}{\Σ}}{\mathrm{\Δl}}_{n,d}^t=0---\left(3\right)$

This constraint refers to, in the total activation cycle of T, the power summation that upwards increases of free section load should be identical with the power summation reduced downwards, namely ensure that the load power consumption in dispatching cycle is constant.The load management mode that the value of T can adopt according to dispatching of power netwoks operations staff is chosen, if the mode that dispatching of power netwoks operations staff adopts load day to close carries out load management, then now T=24, load management is carried out according to the mode closed in week, then now T=168, load management (January by 30 days calculate) is carried out, then now T=720 according to the mode closed by the moon.

4th step, sets up the region load balancing model based on load model.

(4) region account load balancing constraints

$P_{all.n}^t+P_w^{t,n}+\underset{nn=1}{\overset{N}{\Σ}}{L}_{n,nn}^t=P_{l.n}^t+{\mathrm{\Δl}}_{n,u}^t-{\mathrm{\Δl}}_{n,d}^t---\left(4\right)$

In formula, the gross power sum of t period all conventional power unit; then represent the electric load of t period.In formula for dominant eigenvalues size between t region n and region nn.Setting current reference direction is: inflow region is positive direction, and outflow region is negative direction.So can positive and negative values be got, positive and negative, represent the direction of power delivery. for the wind power generation power that the n region t period receives.

5th step, considers the wind power output temporal characteristics of electrical network, load temporal characteristics, peak load regulation characteristic, electrical network sends the factors such as ability, set up provincial power network time stimulatiom model, optimize the power balance of the whole network containing wind-powered electricity generation.Provincial power network time stimulatiom model is consistent with document one " the wind-powered electricity generation annual plan formulating method based on time stimulatiom " (Automation of Electric Systems the 38th volume o. 11th the 13rd page), only does simple introduction here.

(5) set optimization power constraint

$X_j^t\·P_{j,\min }\≤P_j\left(t\right)\≤P_{j,\max }\·X_j^t---\left(5\right)$

(6) minimum start and stop time-constrain

$Y_j^t+\underset{i=1}{\overset{k_{on}}{\Σ}}{Z}_j^{t+i}\≤1---\left(6\right)$

$Z_j^t+\underset{i=1}{\overset{k_{off}}{\Σ}}{Y}_j^{t+i}\≤1---\left(7\right)$

(7) heat supply phase thermal power plant unit units limits

$P_{j,BY}^t=C_j^b\·H_j^t---\left(8\right)$

$H_j^t\·C_j^b\≤P_{j,CQ}^t\≤P_{j,max}-H_j^t\·C_j^v---\left(9\right)$

(8) start and stop logic state constraint

$\left\{\begin{array}{c}{X}_j^t-X_j^{t-1}-Y_j^t+Z_j^t=0\\ -X_j^t-X_j^{t-1}+Y_j^t\≤0\\ X_j^t+X_j^{t-1}+Y_j^t\≤2\\ -X_j^t-X_j^{t-1}+Z_j^t\≤0\\ X_j^t+X_j^{t-1}+Z_j^t\≤2\end{array}\right.---\left(10\right)$

(9) unit climbing rate constraint

$P_j^{t+1}-P_j^t\≤{\mathrm{\ΔP}}_{j,up}---\left(11\right)$

$P_j^t-P_j^{t+1}\≤{\mathrm{\ΔP}}_{j,down}---\left(12\right)$

(10) spinning reserve constraint

$\begin{array}{c}-\underset{j=1}{\overset{J}{\Σ}}{P}_{j,max}\·X_j^t\≤-P_l^t-\mathrm{Pr}e\\ \underset{j=1}{\overset{J}{\Σ}}{P}_{j,min}\·X_j^t\≤P_l^t-Nre\end{array}---\left(13\right)$

(11) interregional line transmission capacity-constrained

$L_{n,nn}^{\min }\≤L_{n,nn}^t\≤L_{n,nn}^{\max }---\left(14\right)$

(12) wind power constraint

$0\≤P_{w,n}\left(t\right)\≤P_{w,n}^{*}\left(t\right)---\left(15\right)$

(13) target function

$max\underset{t=1}{\overset{T}{\Σ}}\underset{n=1}{\overset{N}{\Σ}}{P}_{w,n}\left(t\right)---\left(16\right)$

In formula, P _{j, max}, P _{j, min}be respectively the exert oneself upper limit and the lower limit of exerting oneself of jth platform unit. be respectively and represent that unit j starts at period t, the binary variable of stopped status, for " 1 " represents that unit starts, for " 0 " represents that unit is not at starting state, for " 1 " represents that unit is shut down, for " 0 " represents that unit is not in stopped status; k _onthe machine time is opened for unit is minimum; k _offfor unit minimum downtime; That reflects the minimum time span opening machine or shutdown, dissimilar Unit Commitment machine time parameter is different. for back pressure unit is exerted oneself size; for unit output size of bleeding; for t period load of heat; for thermal power plant unit coupled thermomechanics coefficient. be respectively the maximum upper creep speed of unit j and lower creep speed.Pre and Nre is respectively positive rotation spinning reserve for subsequent use and negative. for the tie-line power transmission upper limit between t region n and region nn, for tie-line power transmission lower limit between t region n and region nn. for wind-powered electricity generation theory is exerted oneself.

6th step, finally by the Mathematical Modeling adopting said new method to set up, the power load of discontinuity surface when optimizing each, by the transfer of load power, unit start-up mode can be optimized, thus significantly reduce wind-powered electricity generation abandon wind-powered electricity generation amount, optimum results can be dispatching of power netwoks personnel provide instruct and suggestion.

Above embodiment is only in order to illustrate that technical scheme of the present invention is not intended to limit; although with reference to above-described embodiment to invention has been detailed description; those of ordinary skill in the field still can modify to the specific embodiment of the present invention or equivalent replacement; and these do not depart from any amendment of spirit and scope of the invention or equivalent replacement, it is all being applied within the claims of the present invention awaited the reply.

Claims (7)

1. for the production of a network load modeling method for analog simulation, it is characterized in that, described method comprises the steps:

The load transfer plan power upper limit model of discontinuity surface when step 1. is set up each;

Step 2. sets up the load transfer plan electricity Upper-Bound Model in each cycle;

Step 3. sets up different load cleancut models;

Step 4. sets up the region load balancing model based on load side peak regulation;

Step 5. sets up the load side peak regulation model improving wind electricity digestion capability, the load side peak regulation power of discontinuity surface when optimizing each.

2. the method for claim 1, is characterized in that, described step 1 comprises:

According to the performance number that the load of discontinuity surface time each increases and declines, the load transfer plan power upper limit model of discontinuity surface when setting up each:

$\left\{\begin{array}{c}{\mathrm{\Δl}}_{n,u}^t\≤{\mathrm{\Δl}}_{n,max}^t\\ {\mathrm{\Δl}}_{n,d}^t\≤{\mathrm{\Δl}}_{n,max}^t\end{array}\right.---\left(1\right)$

In formula (1), for t region n load increases performance number; for t region n load decline performance number, and with be positive variable; for the transferable power upper limit of t region n load.

3. method as claimed in claim 2, it is characterized in that, described step 2 comprises:

According to the described load transfer plan power upper limit model of discontinuity surface time each, set up the load transfer plan electricity Upper-Bound Model in each cycle:

$\underset{t=1}{\overset{T}{\Σ}}({\mathrm{\Δl}}_{n,u}^t+{\mathrm{\Δl}}_{n,d}^t)/2\≤Q_n---\left(2\right)$

In formula (2), T is the total activation cycle; Q _nfor the total electricity of load transfer plan.

4. method as claimed in claim 3, it is characterized in that, described step 3 comprises:

According to the described load transfer plan power upper limit model of discontinuity surface time each, set up different load cleancut models:

$\underset{t=1}{\overset{T}{\Σ}}{\mathrm{\Δl}}_{n,u}^t-\underset{t=1}{\overset{T}{\Σ}}{\mathrm{\Δl}}_{n,d}^t=0---\left(3\right).$

5. method as claimed in claim 4, it is characterized in that, described step 4 comprises:

Set up the region load balancing model based on load side peak regulation:

$P_{all.n}^t+P_w^{t,n}+\underset{nn=1}{\overset{N}{\Σ}}{L}_{n,nn}^t=P_{l.n}^t+{\mathrm{\Δl}}_{n,u}^t-{\mathrm{\Δl}}_{n,d}^t---\left(4\right)$

In formula (4), it is the gross power sum of all conventional power unit of t; it is the electric load of t; for t, dominant eigenvalues value between region n and region nn, and value is timing, and electric current inflow region is positive direction; value is for time negative, and electric current inflow region is negative direction; be the wind power generation power of t, region n receiving.

6. method as claimed in claim 5, it is characterized in that, described step 5 comprises:

5-1., according to the model in step 1 to 4, sets up the load side peak regulation model improving wind electricity digestion capability;

5-2. according to described load side peak regulation model, the load side peak regulation power of discontinuity surface when optimizing each.

7. method as claimed in claim 6, it is characterized in that, described 5-1 comprises: according to the model in step 1 to 4, sets up constraints and the target function of the load side peak regulation model improving wind electricity digestion capability:

A. set optimization power constraint:

$X_j^t\·P_{j,min}\≤P_j\left(t\right)\≤P_{j,max}\·X_j^t---\left(5\right)$

In formula (5), for unit j is at the binary variable of t; P _{j, max}, P _{j, min}be respectively the exert oneself upper limit and the lower limit of exerting oneself of jth platform unit; P _jt () is set optimization power;

B. minimum start and stop time-constrain:

$\left\{\begin{array}{c}{Y}_j^t+\underset{i=1}{\overset{k_{on}}{\Σ}}{Z}_j^{t+i}\≤1\\ Z_j^t+\underset{i=1}{\overset{k_{off}}{\Σ}}{Y}_j^{t+i}\≤1\end{array}\right.---\left(6\right)$

In formula (6), be respectively and represent that unit j starts in t, the binary variable of stopped status, for " 1 " represents that unit starts, for " 0 " represents that unit is not at starting state, for " 1 " represents that unit is shut down, for " 0 " represents that unit is not in stopped status; k _onthe machine time is opened for unit is minimum; k _offfor unit minimum downtime; I is for calculating variable;

C. heat supply phase thermal power plant unit units limits:

$\left\{\begin{array}{c}{P}_{j,BY}^t=C_j^b\·H_j^t\\ H_j^t\·C_j^b\≤P_{j,CQ}^t\≤P_{j,\max }-H_j^t\·C_j^v\end{array}\right.---\left(7\right)$

In formula (7), for back pressure unit is exerted oneself size; for unit output size of bleeding; for t load of heat; be thermal power plant unit coupled thermomechanics coefficient;

D. start and stop logic state constraint:

$\left\{\begin{array}{c}{X}_j^t-X_j^{t-1}-Y_j^t+Z_j^t=0\\ -X_j^t-X_j^{t-1}+Y_j^t\≤0\\ X_j^t+X_j^{t-1}+Y_j^t\≤2\\ -X_j^t-X_j^{t-1}+Z_j^t\≤0\\ X_j^t+X_j^{t-1}+Z_j^t\≤2\end{array}\right.---\left(8\right)$

In formula (8), for unit j is at the binary variable in t-1 moment;

E. unit climbing rate constraint:

$\left\{\begin{array}{c}{P}_j^{t+1}-P_j^t\≤{\mathrm{\ΔP}}_{j,up}\\ P_j^t-P_j^{t+1}\≤{\mathrm{\ΔP}}_{j,down}\end{array}\right.---\left(9\right)$

In formula (9), be respectively the maximum upper creep speed of unit j and lower creep speed; for t unit j power; for t+1 moment unit j power;

F. spinning reserve constraint:

$\left\{\begin{array}{c}-\underset{j=1}{\overset{J}{\Σ}}{P}_{j,max}\·X_j^t\≤-P_l^t-\mathrm{Pr}e\\ \underset{j=1}{\overset{J}{\Σ}}{P}_{j,min}\·X_j^t\≤P_l^t-Nre\end{array}\right.---\left(10\right)$

In formula (10), J is unit sum; Pre and Nre is respectively positive rotation spinning reserve for subsequent use and negative; it is the electric power of t;

G. interregional line transmission capacity-constrained:

$L_{n,nn}^{min}\≤L_{n,nn}^t\≤L_{n,nn}^{\max }---\left(11\right)$

In formula (11), for the tie-line power transmission upper limit between t region n and region nn, for tie-line power transmission lower limit between t region n and region nn;

H. wind power constraint:

$0\≤P_{w,n}\left(t\right)\≤P_{w,n}^{*}\left(t\right)---\left(12\right)$

In formula (12), P _w,nt () is wind power; for wind-powered electricity generation theory is exerted oneself;

I. target function:

$max\underset{t=1}{\overset{T}{\Σ}}\underset{n=1}{\overset{N}{\Σ}}{P}_{w,n}\left(t\right)---\left(13\right)$

In formula (13), N is the sum in region.