CN101710702A - 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 to improve the quality of economic growth and the important channel of benefit, and be the fundamental state policy of country, be 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 that the optimization generation schedule of the whole network total consumption of coal minimum of security constraint is satisfied in realization, promptly ensures the comprehensive total consumption of coal minimum of the equivalent coal consumption sum of gross coal consumption rate and via net loss after the via net loss conversion is standard coal consumption.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

At the prior art above shortcomings, 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 maximum output 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 that generator output is adjusted 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 target function of the second period generator output adjustment amount minimum:

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 by 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 maximum output 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 by 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, described h ₁Value be 1, h ₂Value be 1.25;

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

Further, in described step 4), described 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 by 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, promptly ensure the comprehensive total consumption of coal minimum of the equivalent coal consumption sum of gross coal consumption rate and via net loss.

Description of drawings

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

Fig. 2 is the topological structure schematic diagram 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 the drawings and specific embodiments 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 schematic diagram of power system network, and Fig. 3 is the external system isopleth map of electrical network, as shown in the figure.The concrete steps of the method for realization dynamic energy-saving distribution of electrical power system are as follows:

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 maximum output 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 is

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)

As 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)

As 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

As some moment t:

The generator node number	????a	????b	????c
				????85	????0.00065	????0.01390	????41.84502

The generator node number	????a	????b	????c
				????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 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 at 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)$

As some moment t:

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

As some moment t:

The generator node number	????P Gimin	????P Gi	????P Gimax
				????85	????150	????170	????300

The generator node number	????P Gimin	????P Gi	????P Gimax
				????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 by 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 maximum output 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 finished on the basis of following two periods generation schedule optimization:

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 that generator output is adjusted 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, carry out the load calculating of each concrete load bus again in conjunction with system's total load of prediction 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 calculate by following formula 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 by 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 carries out each concrete load point in conjunction with system's total load of prediction calculates again.

If system's total load that the t period is predicted is PD (t), the load of the load point k that obtains by real time data is PRTDk (t), and so, load point k can calculate by following formula 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) the external system equivalence in the dynamic energy-saving distribution

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 mutually, 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.

By external province electrical network equivalence, 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 as follows:

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 example 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

Illustrate:

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 maximum output of gaining merit

The meritorious minimum of PMIN-generator is exerted oneself

QMAX-generator reactive maximum output

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 maximum output of Qgcup-generator failure

The idle minimum of Qgclo-generator failure is exerted oneself

Ownregn-node location

The real system example 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

Illustrate:

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 situations to embody the network topology circuit and change (transformer that has four kinds of situations 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 represents that this branch road cut-offs (corresponding situation (3)) at this time point

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

BRSTAT=2 represents 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 example 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

Illustrate:

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 maximum output of gaining merit

AGCSTATUS-generator AGC state:

AGCSTATUS=1, expression AGC state, promptly adjustable unit

AGCSTATUS=0 represents 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 example 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

Illustrate:

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 represents 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 example 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

Illustrate:

The TIME-time period

TotalLoad-this time period system total load

The real system example 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 saves 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 the present invention is had been described in detail with reference to preferred embodiment, those of ordinary skill in the art is to be understood that, 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 solution of the present invention, it all should be encompassed in the middle of the claim scope of the present invention.

Claims (4)

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, 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 maximum output 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??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 that generator output is adjusted 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{\&Sigma;}}}{F}_i\left(P_{\mathrm{Gi}}\right)+\mathrm{\&gamma;}{P}_L$

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

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

Constraints comprises 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??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 by 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 maximum output 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 by 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: described 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 is characterized in that: in described step 3), the generator output adjustment amount that calculates goes out the sharing of load factor by being connected with the SCADA database to check with corrected Calculation.

4. the method for realization dynamic energy-saving distribution of electrical power system according to claim 1, it is characterized in that: in described step 4), described 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