CN101710702B - Method for realizing dynamic energy-saving distribution of electrical power system - Google Patents

Abstract

The invention discloses a method for realizing dynamic energy-saving distribution of an electrical power system, which is implemented by collecting data from the EMS systems of each node of a generator in certain element power grid; substituting the data into an objective function of multi-period dynamic energy-saving distribution; equaling the connecting lines for connecting other element power grids with the element power grid to obtain a new generation scheme; and finally feeding the new generation scheme back to the EMS systems of each node of the generator in the certain element power grid. The method can convert grid loss into standard coal consumption, and formulates a generation scheme with minimum multi-period generation coal consumption according to the change of the load of the electrical power system while satisfying the generator set security constraint and grid power transmission security constraint, namely, comprehensive coal consumption totaling the generation coal consumption and the equivalent coal consumption of grid loss is secured.

Description

Realize the method for dynamic energy-saving distribution of electrical power system

Technical field

The invention belongs to the dispatching automation of electric power systems technical field, relate in particular to a kind of method that realizes dynamic energy-saving distribution of electrical power system.

Background technology

Energy-saving and emission-reduction are important measures of alleviating energy supply contradiction and environmental constraints contradiction, are the important channels with benefit of improving the quality of economic growth, and are national fundamental state policies, are directly connected to the national sustainable development implementation.Implementing energy-saving distribution in the power-management centre is the effective means of implementing national energy-saving and emission-reduction policy.

The energy-saving distribution target is after the via net loss conversion is standard coal consumption, realizes satisfying the minimum optimization generation schedule of the whole network total consumption of coal of security constraint, ensures that promptly the comprehensive total consumption of coal of the equivalent coal consumption sum of gross coal consumption rate and via net loss is minimum.Because load changes in time, generation schedule also should change according to load variations.

In the prior art, the optimized calculation method in the energy-saving distribution only considers that generally the generating set coal consumption is minimum in a certain calculation interval, can't consider the demand of multi-period optimization; Or in multi-period calculating, can not take into account constraint of generating set security constraint, power grid security and via net loss simultaneously.

Summary of the invention

Above-mentioned deficiency to the prior art existence; The purpose of this invention is to provide a kind of when satisfying generating set security constraint, network transmission of electricity security constraint; Press the variation of power system load, formulate the method for realization dynamic energy-saving distribution of electrical power system of the generation schedule of multi-period minimum gross coal consumption rate.

The method of realization dynamic energy-saving distribution of electrical power system provided by the invention comprises the steps:

1) the burden with power P of acquisition node in the EMS system of each generator node from certain unit electrical network _D, branch road ij constraints P _Ijmax, the last one actual P that exerts oneself of generator bus i constantly _Gi ⁰, this moment generator bus i the actual P that exerts oneself _Gi, generator bus i the minimum P that exerts oneself _Gimin, generator bus i EIAJ P _Gimax, branch road quantity NT and generator quantity NG;

2) data of gathering are brought into the target function of multi-period dynamic energy-saving distribution, the calculating generator adjustment amount of exerting oneself:

minF(t)＝h ₁F ₁(t)+h ₂F ₂(t)

Constraints: $\underset{i\∈\mathrm{NG}}{\mathrm{\Σ}}{P}_{\mathrm{Gi}}\left(t\right)=\underset{k\∈\mathrm{ND}}{\mathrm{\Σ}}{P}_{\mathrm{Dk}}\left(t\right)+P_L\left(t\right)$

P _Gi(t)-P _Gi(t-1)≤ΔP _GRCimax i∈NG

P _Gi(t-1)-P _Gi(t)≤ΔP _GRCimax i∈NG

|P _ij(t)|≤P _ijmax ij∈NT

Wherein, h ₁: the weight factor of energy consumption in the very first time section target function

h ₂: the weight factor of generator output adjustment in the second time period target function

F ₁(t) be the target function of minimum comprehensive coal consumption of first period:

min $F_{}$ ${}_1=\underset{i=1}{\overset{\mathrm{NG}}{\mathrm{\Σ}}}{F}_i\left(P_{\mathrm{Gi}}\right)+\mathrm{\γ}{P}_L$

F ₂(t) be the minimum target function of the second period generator output adjustment amount:

min $F_{}$ ${}_2=\underset{i=1}{\overset{\mathrm{NG}}{\mathrm{\Σ}}}{(P_{\mathrm{Gi}}-P_{\mathrm{Gi}}^0)}^2$

Constraints comprises meritorious balance, generator output restriction and the constraint of circuit trend:

$\underset{i\∈\mathrm{NG}}{\mathrm{\Σ}}{P}_{\mathrm{Gi}}=\underset{k\∈\mathrm{ND}}{\mathrm{\Σ}}{P}_{\mathrm{Dk}}+P_L$

P _Gimin≤P _Gi≤P _Gimax i∈NG

|P _ij|≤P _ijmax ij∈NT

Wherein, P _D: the burden with power of node

system's total load

P _Ij: the trend on the branch road ij, carry out trend through Niu Lafa and calculate

P _Ijmax: the constraints of branch road ij, the heap(ed) capacity of circuit and transformer

P _GI ⁰: last one actual the exerting oneself of generator bus i constantly

P _Gi: actual the exerting oneself of this moment generator bus i

P _Gimin: the minimum of generator bus i is exerted oneself

P _Gimax: the EIAJ of generator bus i

The total consumption of coal of all generating sets of system

P _L: via net loss

F _i: the fuel consumption of generator unit i

NT: branch road quantity

NG: generator quantity

γ: active loss is converted to the conversion coefficient of coal consumption

3) application load distribution factor method is assigned to system's total load the load bus of each power plant;

4) with other unit electrical network through carrying out equivalence, the generation schedule that must make new advances with the joining interconnection of this unit electrical network;

5) new generation schedule is fed back in the EMS system of each generator node in this unit electrical network, as next step implementation plan.

Further, said h ₁Value be 1, h ₂Value be 1.25;

Further, in said step 3), the generator output adjustment amount that calculates goes out the sharing of load factor through being connected with the SCADA database to check with corrected Calculation;

Further, in said step 4), said equivalence is to be duty value or equivalent power supplys such as this unit electrical network is borderline with the Power Exchange between unit electrical network and other unit electrical network is equivalent.

Beneficial effect of the present invention: the present invention is through the data of the EMS system of each generator node in the collecting unit electrical network; And the target function that data are brought multi-period dynamic energy-saving distribution into calculated and through the equivalent generation schedule that must make new advances; New generation schedule is fed back in the EMS system of each generator node in this unit electrical network, implement target as next step.This method can be standard coal consumption with the via net loss conversion; When satisfying generating set security constraint, network transmission of electricity security constraint; Press the variation of power system load; Formulate the generation schedule of multi-period minimum gross coal consumption rate, ensure that promptly the comprehensive total consumption of coal of the equivalent coal consumption sum of gross coal consumption rate and via net loss is minimum.

Description of drawings

Fig. 1 is an overview flow chart of the present invention;

Fig. 2 is the topological structure sketch map of power system network;

Fig. 3 is the external system isopleth map of electrical network;

Fig. 4 is the whole network load chart;

Fig. 5 is the total system total consumption of coal comparison diagram before and after optimizing.

Embodiment

Below in conjunction with accompanying drawing and embodiment the present invention is done to describe in further detail.

Fig. 1 is an overview flow chart of the present invention, and Fig. 2 is the topological structure sketch map of power system network, and Fig. 3 is the external system isopleth map of electrical network, and is as shown in the figure.The concrete steps of the method for realization dynamic energy-saving distribution of electrical power system are following:

1), the burden with power P of acquisition node in the EMS system of each generator node from certain unit electrical network _D, branch road ij constraints P _Ijmax, the last one actual P that exerts oneself of generator bus i constantly _Gi ⁰, this moment generator bus i the actual P that exerts oneself _Gi, generator bus i the minimum P that exerts oneself _Gimin, generator bus i EIAJ P _Gimax, branch road quantity NT and generator quantity NG.In the present embodiment, certain unit electrical network refers to certain provincial power network, and other unit electrical network refers to other provinces electrical network.

2) in the present embodiment, be divided into 96 operating points (per 15 minutes operating points) with one day, one-time calculation goes out the generation schedule optimized distribution plan of 96 operating points.

Different with static energy-saving distribution, when formulating multi-period generation schedule and dynamic electricity generation plan, the meritorious adjusting constraint Δ P of generator _GRCimaxShould take in.The product of adjusting time and governing speed (speed of promptly climbing) is depended in the meritorious adjusting constraint of generator, that is:

ΔP _GRCimax＝ΔP _GRi×Δt i∈NG (1)

Wherein,

Δ P _GRCimax: take into account the meritorious adjusting constraint of the generator i of generator climbing speed

Δ P _GRi: the climbing speed of generator i

Δ t: running time

The meritorious adjusting constraint expression formula of considering the generator i of generator climbing speed does

Generator rises to gain merit and regulates the constraint expression formula

$P_{\mathrm{Gi}}\left(t\right)-P_{\mathrm{Gi}}(t-1)\≤\mathrm{\Δ}{P}_{\mathrm{GRCi}\max }^{\mathrm{up}}$ i∈NG (2)

Like some moment t:

The generator node number	P Gi(t)	P Gi(t-1)	ΔP GRCimax up
				85	45	0	45
86	90	45	45
				87	135	90	45

Generator descends to gaining merit and regulates the constraint expression formula

$P_{\mathrm{Gi}}(t-1)-P_{\mathrm{Gi}}\left(t\right)\≤\mathrm{\Δ}{P}_{\mathrm{GRCi}\max }^{\mathrm{down}}$ i∈NG (3)

Like some moment t:

The generator node number	P Gi(t)	P Gi(t-1)	ΔP GRCimax up
				85	0	45	45
86	45	90	45
				87	90	135	45

Wherein,

Δ P _GRCimax ^Up: take into account the meritorious adjusting constraint of the generator i of generator climbing speed (rising)

Δ P _GRCimax ^Down: take into account the meritorious adjusting constraint of the generator i of generator climbing speed (decline)

T: running time section

P _Gi(t): generator i is in the meritorious output of section running time (t)

P _Gi(t-1): generator i is in the meritorious output of section running time (t-1)

Ordinary circumstance can consider that the generator rate of climb and decrease speed are the same, so the meritorious adjusting constraint that generator rises and descends is equal, promptly

$\mathrm{\Δ}{P}_{\mathrm{GRCi}\max }^{\mathrm{up}}=\mathrm{\Δ}{P}_{\mathrm{GRCi}\max }^{\mathrm{down}}=\mathrm{\Δ}{P}_{\mathrm{GRCi}\max }$ i∈NG (4)

So adjacent period generator rising meritorious adjusting constraint expression formula and the meritorious adjusting of generator decline constraint expression formula can be expressed as:

P _Gi(t)-P _Gi(t-1)≤ΔP _GRCimax i∈NG (5)

P _Gi(t-1)-P _Gi(t)≤ΔP _GRCimax i∈NG (6)

The target function of minimum comprehensive coal consumption of first period:

${\mathrm{<figure-callout\; id="1"\; label="min\; F"\; filenames="H2009101917055E00022.png"\; state="\{\{state\}\}">min</figure-callout>}F}_1=\underset{i=1}{\overset{\mathrm{NG}}{\mathrm{\&Sigma;}}}{F}_i\left(P_{\mathrm{Gi}}\right)+\mathrm{\&gamma;}{P}_L---\left(7\right)$

The mathematical expression of generator consumption characteristic is:

F _i(P _Gi)＝a _iP _Gi ²+b _iP _Gi+c _i

Wherein,

a _i: the characteristic quadratic term of generator consumption

b _i: the characteristic once item of generator consumption

c _i: the characteristic constant term of generator consumption

Like some moment t:

The generator node number	a	b	c
				85	0.00065	0.01390	41.84502
86	0.00065	0.01390	41.84502
				87	0.00065	0.01390	41.84502

The second period generation schedule will be with changing load, and the adjusting target that the present invention proposes is to make the generator output adjustment amount minimum, to improve the stability of unit.Target function is:

${\mathrm{<figure-callout\; id="2"\; label="min\; F"\; filenames="H2009101917055E00022.png"\; state="\{\{state\}\}">min</figure-callout>}F}_2=\underset{i=1}{\overset{\mathrm{NG}}{\mathrm{\&Sigma;}}}{(P_{\mathrm{Gi}}-P_{\mathrm{Gi}}^0)}^2---\left(8\right)$

As being directed against a certain generating set:

Period	P Gi 0	P Gi
			0:00	100	100
0:15	110	108
			0:30	110	108

Constraints comprises meritorious balance, generator output restriction and the constraint of circuit trend (being network transmission of electricity security constraint), that is:

$\underset{i\&Element;\mathrm{NG}}{\mathrm{\&Sigma;}}{P}_{\mathrm{Gi}}=\underset{k\&Element;\mathrm{ND}}{\mathrm{\&Sigma;}}{P}_{\mathrm{Dk}}+P_L---\left(9\right)$

Like some moment t:

P _Gimin≤P _Gi≤P _Gimax i∈NG (10)

Like some moment t:

The generator node number	P Gimin	P Gi	P Gimax
				85	150	170	300
86	150	200	300
				87	150	250	300

|P _ij|≤P _ijmax ij∈NT (11)

As:

Branch road top node number	Branch road end segment period	P ij	P ijmax
				1	20	3600	4091.85
98	1	1600	2045.90
				38	2	-240	400

Wherein:

P _D: the burden with power of node

system's total load

P _Ij: the trend on the branch road ij, carry out trend through Niu Lafa and calculate

P _Ijmax: the constraints of branch road ij, the heap(ed) capacity of branch road (circuit and transformer)

P _Gi ⁰: last one actual the exerting oneself of generator bus i constantly

P _Gi: actual the exerting oneself of this moment generator bus i

P _Gimin: the minimum of generator bus i is exerted oneself

P _GimaxThe EIAJ of generator bus i

The total consumption of coal of

all generating sets of system

F _i: the fuel consumption of generator unit i

NT: branch road quantity

NG: generator quantity

γ: active loss is converted to the conversion coefficient of coal consumption

Therefore, the Mathematical Modeling of multi-period dynamic energy-saving distribution can continuous iterative computation be accomplished on the basis that following two period generation schedules are optimized:

minF(t)＝h ₁F ₁(t)+h ₂F ₂(t) (12)

Constraints:

$\underset{i\&Element;\mathrm{NG}}{\mathrm{\&Sigma;}}{P}_{\mathrm{Gi}}\left(t\right)=\underset{k\&Element;\mathrm{ND}}{\mathrm{\&Sigma;}}{P}_{\mathrm{Dk}}\left(t\right)+P_L\left(t\right)---\left(13\right)$

P _Gi(t)-P _Gi(t-1)≤ΔP _GRCimax i∈NG (14)

P _Gi(t-1)-P _Gi(t)≤ΔP _GRCimax i∈NG (15)

|P _ij(t)|≤P _ijmax ij∈NT (16)

Wherein:

H1: the weight factor of energy consumption in the very first time section target function

H2: the weight factor of generator output adjustment in the second time period target function

Usually the value of h2 is greater than the value of h1, the h of present embodiment ₁Value be 1, h ₂Value be 1.25.

3) the sharing of load factorization method in the dynamic energy-saving distribution

No matter it is that system safety is checked that energy-saving distribution calculates, still one day 96 economic dispatch is calculated, and all needs the load data of each load bus of using system.The load of predicting each node accurately is to calculate one of key.

The present invention is assigned to the system prediction total load method of each concrete load bus:

Prognoses system total load PD at first; Carry out match according to each concrete load bus in the historical load data of each period then; But recursion obtains the load value of each concrete load bus in each period; With its sharing of load factor, combine system's total load of prediction to carry out the load calculating of each concrete load bus again as this load bus.

If system's total load that the t period is predicted is PD (t), the sharing of load factor of load point k is PDFk (t), and so, load point k can go out through computes at the load value of period t

$P_{\mathrm{Dk}}\left(t\right)=\frac{P_{\mathrm{DFk}}\left(t\right)}{\underset{k=1}{\overset{\mathrm{ND}}{\mathrm{\&Sigma;}}}{P}_{\mathrm{DFk}}\left(t\right)}\&times;P_D\left(t\right)$ i∈ND

(17)

During actual the use through be connected the sharing of load factor that can check to go out with the EMS system database with corrected Calculation; That is: with the same day in the real time data each concrete load point at the load value of each period as the sharing of load factor, the load that combines system's total load of prediction to carry out each concrete load point again calculates.

If system's total load that the t period is predicted is PD (t), the load of the load point k that obtains through real time data is PRTDk (t), and so, load point k can go out through computes at the load value of period t during energy-saving distribution calculated

$P_{\mathrm{Dk}}\left(t\right)=\frac{P_{\mathrm{RTDk}}\left(t\right)}{\underset{k=1}{\overset{\mathrm{ND}}{\mathrm{\&Sigma;}}}{P}_{\mathrm{RTDk}}\left(t\right)}\&times;P_D\left(t\right)$ i∈ND

(18)

4) external system in the dynamic energy-saving distribution is equivalent

Provincial power network is the important component part of large regional grid company, and each is economized, and electrical network is coupled to each other between net, support each other, but is being again relatively independent aspect energy scheduling management, the benefit accounting.Provincial power network saving energy and decreasing loss, energy source optimization dispatching management serve as main carrying out with himself electrical network.Provincial power network is referred to as built-in system, is referred to as external system with the joining other provinces of this electrical network electrical network.Since in the net each to economize between net electrical network be to be coupled to each other, need to carry out equivalence to external system.

Concrete grammar is each province Netcom to be crossed with the joining interconnection of this provincial power network carry out equivalence; According to Power Exchange between this certain provincial power network and other provinces electrical network and power supply agreement; With other provinces network equivalence is borderline duty value or the equivalent power supply of waiting of certain provincial power network; As shown in Figure 3, P is an active power among the figure, and Q is a reactive power.

Equivalent through external province electrical network, can under the prerequisite that keeps certain provincial power network total system topology and characteristic, simplify network configuration, reduce amount of calculation.

Example:

Required system's basic model and the relevant parameter of dynamic energy-saving distribution is following:

1, system's basic parameter

A. branch road parameter

IFROM，ITO，R，X，B，TAP，SMAX，SMAXC，TAPMAX，TAPMIN

Wherein:

IFROM-branch road top node number,

ITO-branch road end segment period

R-branch road resistance

X-branch road reactance

B-line charging power (1/2)

TAP-transformer branch road no-load voltage ratio

SMAX-branch road maximum power output constraint

SMAXC-branch trouble maximum power output constraint

The maximum no-load voltage ratio of TAPMAX-transformer branch road

The minimum no-load voltage ratio of TAPMIN-transformer branch road

The real system instance is economized by China:

As:

IFROM，ITO，R，X，B，TAP，SMAX，SMAXC，TAPMAX，TAPMIN

1 20 0.00002 0.00053 0.1020 0.000 4091.85 4091.85 0.00 0.00

98 1 0.00060 0.00830 0.4688 0.000 2045.90 2045.90 0.00 0.00

38 2 0.00000 0.00010 0.0000 0.000 400.00 400.00 0.00 0.00

40 38 0.00050 0.00640 0.3408 0.000 2119.58 2119.58 0.00 0.00

39 37 0.00049 0.00646 1.4766 0.000 4239.16 4239.16 0.00 0.00

36 38 0.00021 0.00355 0.7878 0.000 4239.16 4239.16 0.00 0.00

B. node parameter:

BUSTYP，IBUS，PGEN，QGEN，PLOAD，QLOAD，PMAX1，PMIN1，QMAX1，QMIN1，VOLT，ANGLE，Vnup，Vnlo，Vcup，Vclo，Qgcup，Qgclo，Ownregn

Explain:

BUSTYP-node type

IBUS-node number

PGEN-generator is meritorious exerts oneself

QGEN-generator reactive is exerted oneself

PLOAD-node burden with power

QLOAD-node load or burden without work

PMAX-generator EIAJ of gaining merit

The meritorious minimum of PMIN-generator is exerted oneself

QMAX-generator reactive EIAJ

QMIN-generator reactive minimum is exerted oneself

VOLT-node initial voltage

ANGLE-node initial voltage phase angle

Vnup-node voltage the upper limit

Vnlo-node voltage lower limit:

Vcup-node failure upper voltage limit

Vclo-node failure lower voltage limit

The idle EIAJ of Qgcup-generator failure

The idle minimum of Qgclo-generator failure is exerted oneself

Ownregn-node location

The real system instance is economized by China:

Illustrate:

BUSTYP，IBUS，PGEN，QGEN，PLOAD，QLOAD，PMAX1，PMIN1，QMAX1，QMIN1，VOLT，ANGLE，Vnup，Vnlo，Vcup，Vclo，Qgcup，Qgclo，Ownregn

1 1 0.00 0.00 1800.60 45.00 0.00 0.00 0.00 0.00 1.050 0.00

1.100 0.900 1.100 0.900 0.00 0.00 1

1 2 0.00 0.00 1800.60 111.00 0.00 0.00 0.00 0.00 1.050

0.00 1.100 0.900 1.100 0.900 0.00 0.00 1

2 3 400.00 0.00 42.00 21.00 600.00 0.00 290.00 -90.00

1.000 0.00 1.100 0.900 1.100 0.900 290.00 -90.00 1

2 4 400.00 0.00 42.00 21.00 600.00 0.00 290.00 -90.00

1.000 0.00 1.100 0.900 1.100 0.900 290.00 -90.00 1

2 5 400.00 0.00 42.00 21.00 600.00 0.00 290.00 -90.00

1.000 0.00 1.100 0.900 1.100 0.900 290.00 -90.00 1

2 6 420.00 0.00 42.00 21.00 600.00 0.00 290.00 -90.00

1.000 0.00 1.100 0.900 1.100 0.900 290.00 -90.00 1

C. system network architecture parameter:

TIME，IFROMB，ITOB，R，X，B，TAP，SMAX，SMAXC，TAPMAX，TAPMIN，BRSTAT

Explain:

System network architecture parameter and the difference of front a. branch road parameter: front a. branch road parameter is the basic model parameter of every day, is static; The parameter here is dynamic, promptly only imports the parameter of localized variation.

There are following four kinds of situation to embody the network topology circuit and change (transformer that has four kinds of situation to embody localized network equally changes):

(1), the network parameter of certain time point and the same day basic model just the same, do not import any circuit-switched data this moment;

(2), the network parameter of certain time point with the same day basic model compare, have the parameter of several branch roads to change;

(3), the network parameter of certain time point with the same day basic model compare, have branch road to stop transport;

(4), the network parameter of certain time point with the same day basic model compare, have new branch road to put into operation;

First kind of situation, this time point is not imported any circuit-switched data;

For reaching (2)-(4) purpose, one line state (BRSTAT) of circuit-switched data end input

BRSTAT=0 representes that this branch road cut-offs (corresponding situation (3)) at this time point

BRSTAT=1 representes that this branch road is in this time point parameter update (corresponding situation (2))

BRSTAT=2 representes that this branch road increases newly at this time point, does not originally have (corresponding situation (4)) in the basic model

Wherein:

TIME-time period

IFROMB-branch road top node number

ITOB-branch road end segment period

R-branch road resistance

X-branch road reactance

B-line charging power (1/2)

TAP-transformer branch road no-load voltage ratio

SMAX-branch road maximum power output constraint

SMAXC-branch trouble maximum power output constraint

The maximum no-load voltage ratio of TAPMAX-transformer branch road

The minimum no-load voltage ratio of TAPMIN-transformer branch road

BRSTAT-Zhi line state

The real system instance is economized by China:

Illustrate:

TIME，IFROMB，ITOB，R，X，B，TAP，SMAX，SMAXC，TAPMAX，TAPMIN，BRSTAT

1 154 184 0.00860 0.04040 0.0659 0.000 406.33 406.33

0.00 0.00 0

2 154 184 0.00860 0.04040 0.0659 0.000 406.33 406.33

0.00 0.00 0

3 154 184 0.00860 0.04040 0.0659 0.000 406.33 406.33

0.00 0.00 0

4 154 184 0.00860 0.04040 0.0659 0.000 406.33 406.33

0.00 0.00 0

5 154 184 0.00860 0.04040 0.0659 0.000 406.33 406.33

0.00 0.00 0

6 154 184 0.00860 0.04040 0.0659 0.000 406.33 406.33

0.00 0.00 0

2, economic dispatch parameter

PGNO，ALPHA，BETA，GAMMA，PGMIN，PGMAX，AGCSTATUS，PGCTG，PGGRC

Explain:

PGNO-generator node number

ALPHA-generator (factory) coal consumption curve constant term

BETA-generator (factory) coal consumption curve is item once

GAMMA-generator (factory) coal consumption curve quadratic term

The meritorious minimum of PGMI N-generator is exerted oneself

The PGMAX-generator EIAJ of gaining merit

AGCSTATUS-generator AGC state:

AGCSTATUS=1, expression AGC state, promptly adjustable unit

AGCSTATUS=0 representes non-AGC state, promptly non-adjustable unit

The PGCTG-generator is stopped transport or the trouble hunting state:

PGCTG=1, the expression generator puts into operation

PGCTG=0, the expression generator is stopped transport or trouble hunting

PGGRC-generator climbing speed

The real system instance is economized by China:

Illustrate:

PGNO，ALPHA，BETA，GAMMA，PGMI N，PGMAX，AGCSTATUS，PGCTG，PGGRC

85 41.84502 0.01390 0.00065 150.000 300.000 1 1 45.0

86 41.84502 0.01390 0.00065 150.000 300.000 1 1 45.0

87 41.84502 0.01390 0.00065 150.000 300.000 1 1 45.0

88 41.84502 0.01390 0.00065 0.000 200.000 1 1 60.0

90 39.81035 0.01322 0.00062 0.000 200.000 1 1 60.0

91 39.81035 0.01322 0.00062 0.000 200.000 1 1 60.0

3, system's total load and generator parameter:

A. the given initial generation schedule parameter of generator (factory):

TIME，PGNO，PG，AGCST，UNCTGST，V0

Explain:

The TIME-time period

PGNO-generator node number

The given generation schedule of PG-

AGCST-this time period generator AGC state

AGCST=1, expression AGC state, promptly adjustable unit

AGCST=0 representes non-AGC state, promptly non-adjustable unit

UNCTGST-this time period generator is stopped transport or the trouble hunting state

UNCTGST=1, the expression generator puts into operation

UNCTGST=0, the expression generator is stopped transport or trouble hunting

V0-this time period generator voltage

The real system instance is economized by China:

Illustrate:

TIME，PGNO，PG，AGCST，UNCTGST，V0

1 3 0.000 0 1 1.000

1 4 375.000 0 0 1.000

1 5 375.000 0 0 1.000

1 6 0.000 0 1 1.000

1 7 0.000 0 1 0.970

1 8 400.000 0 0 1.015

B. system's total load parameter

TIME，TotalLoad

Explain:

The TIME-time period

TotalLoad-this time period system total load

The real system instance is economized by China:

Illustrate:

TIME，TotalLoad

1 10352.000

2 10325.000

3 10264.000

4 10052.000

5 10044.000

6 9782.000

Practical power systems is economized by China be optimized, comprising 520 load buses, 679 branch roads, 104 generator nodes (with reference to accompanying drawing 2,4 load buses that only draw among the figure, 6 branch road ij, 2 generator nodes, all the other do not draw).

min

min $F_{}$ ${}_2=\underset{i=1}{\overset{\mathrm{NG}}{\mathrm{\&Sigma;}}}{(P_{\mathrm{Gi}}-P_{\mathrm{Gi}}^0)}^2$

min $F\left(t\right)=h_1{F}_1\left(t\right)+h_2{F}_2\left(t\right)=h_1(\underset{i=1}{\overset{\mathrm{NG}}{\mathrm{\&Sigma;}}}{F}_i\left(P_{\mathrm{Gi}}\right)+\mathrm{\&gamma;}{P}_L)+h_2\left(\underset{i=1}{\overset{\mathrm{NG}}{\mathrm{\&Sigma;}}}{(P_{\mathrm{Gi}}-P_{\mathrm{Gi}}^0)}^2\right)$

Get h in the present embodiment ₁=1, h ₂=1.25.

γ: coal consumption conversion factor (being the coal tonnage that per 1000 kilowatt hours consume) value 0.53

F _i(P _Gi)＝a _iP _Gi ²+b _iP _Gi+c _i

Wherein:

a _i: the characteristic quadratic term of generator consumption

b _i: the characteristic once item of generator consumption

c _i: the characteristic constant term of generator consumption

The generator node number	a	b	c
				85	0.00065	0.01390	41.84502
86	0.00065	0.01390	41.84502
				87	0.00065	0.01390	41.84502

Got a bit in per 15 minutes, whole day 24 hours was exactly 96 points so.

To a certain generating set:

Period	P Gi 0	P Gi
			0:00	100	100
0:15	110	108
			0:30	110	108

After the optimization before optimizing

System's saving ratio of economizing on coal after system optimizes before the network optimizationization of the net electrical network of the total negative electricity net of the whole network

The total consumption of lotus (MW) damage (MW) damage (MW) system coal is always united, and the consumption coal always consumes coal

Period (ton) (ton) (ton)

0:00 10092 185.975 184.376 881.403 857.2365 24.16654 2.74％

0:15 10065 185.975 184.376 881.298 857.1367 24.16133 2.74％

0:30 10004 185.975 184.376 880.899 856.92 23.97905 2.72％

0:45 9792 185.975 184.376 880.388 856.4532 23.93475 2.72％

1:00 9784 185.218 183.626 877.687 854.3294 23.35761 2.66％

1:15 9522 173.328 171.751 840.666 821.775 18.89098 2.25％

1:30 9522 169.903 168.326 831.59 813.1909 18.3991 2.21％

1:45 9377 169.153 167.576 828.549 810.7939 17.75513 2.14％

2:00 9320 169.153 167.576 828.338 810.5963 17.7417 2.14％

2:15 9095 167.902 166.326 823.684 806.8153 16.8687 2.05％

2:30 8992 165.905 164.326 816.67 801.6251 15.0449 1.84％

2:45 8852 157.73 156.151 793.979 781.9872 11.99175 1.51％

3:00 8652 151.08 149.501 777.224 766.706 10.51803 1.35％

3:15 8619 151.08 149.501 777.11 766.5971 10.5129 1.35％

3:30 8618 151.08 149.501 777.107 766.5941 10.51292 1.35％

3:45 8649 151.08 149.501 777.213 766.6959 10.51711 1.35％

4:00 8600 151.08 149.501 777.047 766.5366 10.51038 1.35％

4:15 8584 151.08 149.501 776.994 766.4862 10.50778 1.35％

4:30 8634 151.08 149.501 777.161 766.6465 10.5145 1.35％

4:45 8613 151.08 149.501 777.09 766.578 10.51205 1.35％

5:00 8614 151.08 149.501 777.094 766.581 10.51302 1.35％

5:15 8606 153.08 151.501 783.448 770.8821 12.56589 1.60％

5:30 8584 153.08 151.501 783.38 770.8176 12.5624 1.60％

5:45 8619 153.83 152.251 786.228 772.5655 13.66253 1.74％

6:00 8786 155.827 154.251 793.239 777.5067 15.73231 1.98％

6:15 8820 155.827 154.251 793.358 777.6186 15.73942 1.98％

6:30 8855 156.052 154.451 793.824 778.181 15.64296 1.97％

6:45 9056 156.177 154.576 794.82 779.2152 15.60475 1.96％

7:00 9290 165.172 163.5757 819.931 800.7794 19.15161 2.34％

7:15 9628 175.297 173.7007 851.853 827.4984 24.35455 2.86％

7:30 10106 175.302 173.7007 853.874 829.1133 24.76074 2.90％

7:45 10649 175.302 173.7007 855.671 830.8632 24.80785 2.90％

8:00 11068 180.547 178.9495 872.018 845.6273 26.39067 3.03％

8:15 11760.9 196.6 195.0218 918.723 893.2584 25.46465 2.77％

8:30 12402.9 208.09 206.497 952.68 932.1823 20.49773 2.15％

8:45 12940.9 215.491 213.897 978.681 959.7067 18.97428 1.94％

9:00 13336.9 216.626 215.022 983.797 964.9564 18.84062 1.92％

9:15 13502.9 216.626 215.022 984.849 965.9886 18.86042 1.92％

9:30 13646.9 216.626 215.022 985.186 966.3081 18.87789 1.92％

9:45 13814.9 216.626 215.022 986.426 967.5278 18.89821 1.92％

10:00 13848.9 216.851 215.247 987.467 968.563 18.90399 1.91％

10:15 13848.9 216.851 215.247 987.467 968.563 18.90399 1.91％

10:30 13882.9 216.851 215.247 987.741 968.8332 18.90785 1.91％

10:45 13882.9 216.851 215.247 987.741 968.8332 18.90785 1.91％

11:00 13916.9 217.626 216.022 990.675 971.7644 18.91058 1.91％

11:15 13885.9 216.851 215.247 989.278 968.8573 20.42066 2.06％

11:30 13680.9 216.626 215.022 987.285 966.8857 20.3993 2.07％

11:45 13380.9 216.626 215.0223 985.36 964.9382 20.42176 2.07％

12:00 13118.9 216.626 215.022 983.874 963.4182 20.45582 2.08％

12:15 12980.9 216.626 215.0217 982.979 962.4807 20.49826 2.09％

12:30 12875.9 216.626 215.0217 982.285 961.763 20.52196 2.09％

12:45 12763.9 216.626 215.022 981.576 961.0252 20.55082 2.09％

13:00 12702.9 216.626 215.0215 981.15 960.5635 20.58648 2.10％

13:15 12868.9 216.626 215.022 981.851 961.2873 20.56374 2.09％

13:30 13000.9 216.626 215.022 982.044 961.4637 20.58034 2.10％

13:45 13000.9 216.626 215.022 982.044 961.4637 20.58034 2.10％

14:00 13000.9 216.626 215.022 982.044 961.4637 20.58034 2.10％

14:15 13046.9 216.626 215.022 982.28 961.6955 20.5845 2.10％

14:30 13076.9 216.626 215.022 982.44 961.8507 20.58926 2.10％

14:45 13074.9 216.626 215.022 982.429 961.8407 20.58834 2.10％

15:00 13074.9 216.626 215.022 982.429 961.8407 20.58834 2.10％

15:15 13000.9 216.626 215.0215 981.923 961.3004 20.62263 2.10％

15:30 12966.9 216.626 215.0215 981.748 961.129 20.61899 2.10％

15:45 12864.9 216.626 215.0217 981.146 960.5 20.64598 2.10％

16:00 12694.9 216.626 215.0225 980.196 959.484 20.71205 2.11％

16:15 12838.9 216.626 215.0225 980.621 959.8902 20.73082 2.11％

16:30 13079.9 216.626 215.0217 981.966 961.2923 20.6737 2.11％

16:45 13203.9 216.626 215.0215 982.726 962.0775 20.64846 2.10％

17:00 13340.9 216.626 215.022 983.373 962.766 20.607 2.10％

17:15 13476.9 216.626 215.022 984.224 963.6047 20.61934 2.10％

17:30 13544.9 216.626 215.0217 984.768 964.1822 20.58576 2.09％

17:45 13578.9 216.626 215.0217 984.246 963.6631 20.58288 2.09％

18:00 13612.9 216.626 215.0217 984.083 963.5008 20.58217 2.09％

18:15 13646.9 216.626 215.0217 984.293 963.7075 20.58553 2.09％

18:30 13716.9 216.626 215.1215 984.842 964.6248 20.21725 2.05％

18:45 13782.9 216.626 215.0217 985.286 964.7235 20.56246 2.09％

19:00 13918.9 216.626 215.0217 986.411 965.9352 20.47585 2.08％

19:15 14168.9 217.001 215.397 989.862 969.3755 20.48654 2.07％

19:30 14186.9 217.451 215.847 991.524 971.0639 20.46014 2.06％

19:45 14069.9 216.876 215.272 989.784 968.1951 21.58891 2.18％

20:00 13885.9 216.626 215.022 988.391 966.2426 22.14841 2.24％

20:15 13602.9 216.626 215.0228 986.915 964.7276 22.18743 2.25％

20:30 13505.9 216.626 215.0228 987.74 965.5703 22.16974 2.24％

20:45 13281.9 216.626 215.0228 986.082 963.8899 22.19208 2.25％

21:00 13082.9 216.626 215.0235 984.729 962.4626 22.26641 2.26％

21:15 12880.9 216.626 215.023 983.476 961.1834 22.29256 2.27％

21:30 12722.9 216.626 215.0233 982.576 960.2591 22.3169 2.27％

21:45 12644.9 216.626 215.0233 982.393 960.0706 22.32239 2.27％

22:00 12371.9 216.626 215.0235 981.115 958.7814 22.33362 2.28％

22:15 12038.9 216.626 215.023 980 957.6625 22.3375 2.28％

22:30 11714 211.248 209.6505 967.142 938.3361 28.80589 2.98％

22:45 11400 199.325 197.7255 925.482 898.5413 26.94072 2.91％

23:00 11008 194.077 192.4755 910.983 882.122 28.86103 3.17％

23:15 10572 194.077 192.4755 910.353 881.6099 28.7431 3.16％

23:30 10068 187.102 185.5005 883.294 860.1607 23.13333 2.62％

23:45 9528 172.227 170.6005 838.78 818.6855 20.0945 2.40％

Whole day amounts to 18,897 18743.79 88386.091 86499.59 1886.496 2.13%

Fig. 4 is the whole network load chart (transverse axis is the time among the figure, and the longitudinal axis is electric weight MW); (transverse axis is the time to Fig. 5 among the figure, and the longitudinal axis is the coal consumption amount ton for the total system total consumption of coal comparison diagram before and after optimizing; 1 line always consumes coal for the system before optimizing, and 2 lines are that the system after optimizing always consumes coal), as shown in the figure: the system after the optimization always consumes coal and always consumes the coal reduction obviously than the system before optimizing.System always consumes 88386.091 tons in coal before optimizing, and optimization back system always consumes 86499.59 tons in coal, optimizes the back and practices thrift 1886.496 tons of coals, and the coal consumption reduction is about 2.134%.

Explanation is at last; Above embodiment is only unrestricted in order to technical scheme of the present invention to be described; Although with reference to preferred embodiment the present invention is specified, those of ordinary skill in the art should be appreciated that and can make amendment or be equal to replacement technical scheme of the present invention; And not breaking away from the aim and the scope of technical scheme of the present invention, it all should be encompassed in the middle of the claim scope of the present invention.

Claims (3)

1. a method that realizes dynamic energy-saving distribution of electrical power system is characterized in that, comprises the steps:

1) the burden with power P of acquisition node in the EMS system of each generator node from certain unit electrical network _D, branch road ij constraints P _Ijmax, last one actual the exerting oneself of generator bus i constantly

The actual P that exerts oneself of this moment generator bus i _Gi, generator bus i the minimum P that exerts oneself _Gimin, generator bus i EIAJ P _Gimax, branch road quantity NT and generator quantity NG;

2) data of gathering are brought into the target function of multi-period dynamic energy-saving distribution, the calculating generator adjustment amount of exerting oneself:

min F(t)＝h ₁F ₁(t)+h ₂F ₂(t)

Constraints: $\underset{i\&Element;\mathrm{NG}}{\mathrm{\&Sigma;}}{P}_{\mathrm{Gi}}\left(t\right)=\underset{k\&Element;\mathrm{ND}}{\mathrm{\&Sigma;}}{P}_{\mathrm{Dk}}\left(t\right)+P_L\left(t\right)$

P _Gi(t)-P _Gi(t-1)≤ΔP _GRCimax i∈NG

P _Gi(t-1)-P _Gi(t)≤ΔP _GRCimax

|P _ij(t)|≤P _ijmax ij∈NT

Wherein, h ₁: the weight factor of energy consumption in the very first time section target function

h ₂: the weight factor of generator output adjustment in the second time period target function

P _Dk(t): load point k is at the load value of period t;

T: running time section;

P _Gi(t): generator i is in the meritorious output of section t running time;

P _Gi(t-1): generator i is in the meritorious output of section t-1 running time;

Δ P _GRCimax: take into account the meritorious adjusting constraint of the generator i of generator climbing speed;

F ₁(t) be the target function of minimum comprehensive coal consumption of first period:

$\min F_1=\underset{i=1}{\overset{\mathrm{NG}}{\mathrm{\&Sigma;}}}{F}_i\left(P_{\mathrm{Gi}}\right)+{\mathrm{\&gamma;P}}_L$

F ₂(t) be the minimum target function of the second period generator output adjustment amount:

$\min F_2=\underset{i=1}{\overset{\mathrm{NG}}{\mathrm{\&Sigma;}}}{(P_{\mathrm{Gi}}-P_{\mathrm{Gi}}^0)}^2$

The target function of minimum comprehensive coal consumption of first period and the minimum bound for objective function of the second period generator output adjustment amount comprise meritorious balance, generator output restriction and the constraint of circuit trend:

$\underset{i\&Element;\mathrm{NG}}{\mathrm{\&Sigma;}}{P}_{\mathrm{Gi}}=\underset{k\&Element;\mathrm{ND}}{\mathrm{\&Sigma;}}{P}_{\mathrm{Dk}}+P_L$

P _Gimin≤P _Gi≤P _Gimax

|P _ij|≤P _ijmax

Wherein, system's total load

P _Ij: the trend on the branch road ij, carry out trend through Niu Lafa and calculate

P _Ijmax: the heap(ed) capacity of circuit and transformer

The total consumption of coal of

all generating sets of system

P _L: via net loss

F _i: the fuel consumption of generator unit i

γ: active loss is converted to the conversion coefficient of coal consumption

3) application load distribution factor method is assigned to system's total load the load bus of each power plant; The said sharing of load factor is to be obtained through being connected to check with corrected Calculation with the SCADA database by the generator output adjustment amount that calculates; The method of said check and correction be with the same day in the real time data each concrete load point at the load value of each period as the sharing of load factor, the load that combines system's total load of prediction to carry out each concrete load point again calculates;

4) with other unit electrical network through carrying out equivalence, the generation schedule that must make new advances with the joining interconnection of this unit electrical network;

5) new generation schedule is fed back in the EMS system of each generator node in this unit electrical network, as next step implementation plan.

2. the method for realization dynamic energy-saving distribution of electrical power system according to claim 1 is characterized in that: said h ₁Value be 1, h ₂Value be 1.25.

3. the method for realization dynamic energy-saving distribution of electrical power system according to claim 1; It is characterized in that: in said step 4), said equivalence is to be duty value or equivalent power supplys such as this unit electrical network is borderline with the Power Exchange between this unit electrical network and other unit electrical network is equivalent.

Images