CN102903016A - Distributed power generation planning method - Google Patents

Abstract

The invention provides a distributed power generation planning method and belongs to the technical field of power energy saving. The method is mainly used for coordinating relation between power grid planning and distributed power supply configuration. The method includes first introducing a node sensitivity coefficient array and a unit loading increment cost array, building a power grid node/ branch margin capacity cost model and a calculation method thereof, quantizing effect of the distributed power supply grid tying on the power grid node/ branch margin capacity cost and finally conducting benefit calculation and evaluation on two distributed power supply property right models owned by a user and a power grid company based on a quantization model of effect of the distributed power supply grid tying to the power grid node/ branch limit capacity cost.

Description

A kind of distributed generation planning method

Technical field

The invention belongs to the electric power energy-saving technical field, particularly a kind of distributed generation planning method based on power distribution network node/branch road limit Capacity Cost.

Background technology

The present distributed power source object of planning is mainly take indicators of costs such as investment operation expenses as main, reliability, security and feature of environmental protection index are auxiliaryly (or to be used as the sub-goal of low weight coefficient, or as the constraint condition processing), seldom consider utilization factor and the abundant intensity of the network equipment, this will be unfavorable for utilizing distributed generation technology that system resource is distributed, improves rationally the utilization ratio of power supply capacity and the economical operation that realizes power distribution network, under existing situation, this problem is to coordinate problem demanding prompt solution in distribution network planning and the distributed power source layoutprocedure.Searching can reflect truly that the index of node capacity tensity or the line load order of severity is the primary problem that solves.At present for the research of power distribution network node/branch road limit Capacity Cost portrayed the node specific power increase/the equipment dilatation investment that reduces to cause in advance/postpone the cost/income that produces, equipment in the capacity expansion relates to is responsible for circuit, switch and the transformer etc. that this node injecting power transmits, yet it adopts branch road active power to describe line load level and residue of network organization transmission line capability, system losses have been ignored, can not take into account idle transmission and impact, and be only applicable to DC power flow calculating.

Summary of the invention

The objective of the invention is in order to overcome the deficiencies in the prior art, a kind of distributed generation planning method is proposed, by setting up a kind of power distribution network node/branch road limit Capacity Cost model, and be applied to coordinate in distribution network planning and the distributed power source allocation problem, be applicable to AC power flow and consider distributed generation planning method power-supply unit utilization factor and abundant intensity and the impact of reflection reactive power flow.

For achieving the above object, the present invention adopts following technical scheme:

A kind of distributed generation planning method may further comprise the steps:

1) introduces node sensitivity coefficient matrix and specific load increment cost matrix, set up power distribution network node/branch road limit Capacity Cost model;

2) having quantized distributed power source is incorporated into the power networks on the impact of power distribution network node/branch road limit Capacity Cost;

3) be incorporated into the power networks on the quantitative model of power distribution network node/branch road limit Capacity Cost impact based on distributed power source, the user provided for oneself carry out benefit with two kinds of distributed power source property rights of grid company pattern and find the solution and assess.

The introducing node sensitivity coefficient matrix of described step 1) and specific load increment cost matrix, the implementation of having set up power distribution network node/branch road limit Capacity Cost model is:

11) node sensitivity coefficient matrix

The support level that the load of node i increases the electric current increase of branch road k adopts sensitivity coefficient γ _IkExpression:

${\mathrm{\γ}}_{\mathrm{ik}}=\frac{{\∂I}_k}{{\∂P}_{\mathrm{di}}}$

I in the formula _kThe electric current that represents k bar branch road, P _DiThe load that represents i bar branch road;

Sensitivity matrix γ is a sparse matrix; Be expressed as: γ=γ _Ik(Pd); P in the formula _dCharacterize the network load distribution matrix;

12) branch current variable quantity

Load fluctuation on other node affects the curent change of branch road k, the electric current I of branch road k _kBe changed to: ${\mathrm{dI}}_k=\underset{i\∈N_k}{\mathrm{\Σ}}{\mathrm{\γ}}_{\mathrm{ik}}\left(P_d\right)d{P}_{\mathrm{di}};$

N in the formula _kFor utilizing branch road k, need carry out the load bus collection of power delivery, dP _Di=δ _i(t) dt, wherein δ _i(t) be the load growth rate of node i in the different periods;

The current increases amount of branch road k: ${\mathrm{dI}}_k=\underset{i\∈N_k}{\mathrm{\Σ}}{\mathrm{\γ}}_{\mathrm{ik}}\left(P_d\right){\mathrm{\δ}}_i\left(t\right)\mathrm{dt};$

13) the branch road increase-volume time

By the curent change trend of branch road k, obtain branch road k increase-volume time point τ according to following formula _k

${\∫}_0^{{\mathrm{\τ}}_k}\underset{i\∈N_k}{\mathrm{\Σ}}{\mathrm{\γ}}_{\mathrm{ik}}\left(P_d\right){\mathrm{\δ}}_i\left(t\right)\mathrm{dt}=I_k^{\mathrm{cap}}-I_k$

In the formula,

Be the rated capacity of branch road k, I _kBe the electric current under the time coordinate initial point;

14) branch road increase-volume investment

The increase-volume cost of investment need be converted to present worth, the employing formula

Or

The calculating of discounting; In the formula: C _kBe the increase-volume cost of branch road k, r is discount rate, PV _kDiscount value for the increase-volume cost;

15) specific load increment cost matrix

151) node load changes

The load increment of supposing node i is △ P _Di, corresponding new network load distribution matrix is The increase-volume time point that then branch road k is new By formula

Determine, in the formula, For branch road k at the network load distribution matrix

Under current value;

152) specific load increment cost coefficient

Node i is to the specific load increment cost IC of branch road k _KiFor:

In the formula, pf _DiBe the load power factor of node i, R _APBe year value coefficients such as fund, relevant with discount rate with the increase-volume cost of investment receipts time limit; As △ P _Di→ 0 o'clock, IC _KiEquivalence is:

153) specific load increment cost matrix

The arbitrary load bus of system is to the specific load increment cost IC of all branch roads _Ki, obtain specific load increment cost matrix IC and be: IC=(IC _Ki) _{B * N}In the formula, B is a way, and N is the load bus number;

16) node limit Capacity Cost

The marginal Capacity Cost of node i is that node i is supported the increment cost sum LMCC of unit of the branch road that this node power is carried to all _i:

17) branch road limit Capacity Cost

The marginal Capacity Cost of branch road k has utilized branch road k for all and has carried out the node of load power transmission to the increment cost sum BMCC of unit of this branch road _k:

Described step 2) quantification distributed power source be incorporated into the power networks on the impact of power distribution network node/branch road limit Capacity Cost, its quantification manner is:

Behind the distributed power source access power distribution network, to its access node and on every side node load all play the power supporting role, distributed power source is regarded as a negative load, its access exerts an influence to the network load distribution situation, has changed load distribution matrix P _dWith specific load increment cost matrix IC, thus so that the marginal Capacity Cost LMCC of each node _iMarginal Capacity Cost BMCC with each branch road _kVariation has also occured, node limit Capacity Cost variable quantity C _LMCCFor: $C_{\mathrm{LMCC}}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}({\mathrm{LMCC}}_{\mathrm{DG},i}-\mathrm{LM}{\mathrm{CC}}_i);$

LMCC in the formula _{DG, i}The be incorporated into the power networks marginal Capacity Cost of posterior nodal point i of expression distributed power source;

Branch road limit Capacity Cost variable quantity C _BMCCFor:

BMCC in the formula _{DG, k}The be incorporated into the power networks marginal Capacity Cost of rear branch road k of expression distributed power source; The marginal Capacity Cost LMCC of each node _iMarginal Capacity Cost BMCC with each branch road _kVariation be equivalent.

Being incorporated into the power networks on the quantitative model of power distribution network node/branch road limit Capacity Cost impact based on distributed power source of described step 3), the user provided for oneself carry out benefit with two kinds of distributed power source property rights of grid company pattern and find the solution with the implementation of assessing and be:

31) user provides the Capacity Benefit of distributed power source for oneself

Grid company is according to C _LMCCOr C _BMCCQuantize to assess the Capacity Benefit B that the user provides distributed power source for oneself _DG, grid company need compensate user distribution formula power construction expense B _DG:

B _DG=-S _DG×C _LMCC=-S _DG×C _BMCC

Wherein, S _DGRated capacity for distributed power source;

32) the power distribution network increase-volume based on distributed power source delays planning

Consider distributed power source to the increase-volume retarding action of power distribution network, the adoption rate factor is determined its income of bringing to the contribution rate of unit cost, and mathematical model is:

$\min \underset{g=1}{\overset{n_{\mathrm{DG}}}{\mathrm{\Σ}}}{S}_{\mathrm{DG},g}$

In the formula, n _DGBe distributed power source configuration sum; Satisfy trend constraint condition:

$S.T.P_{\mathrm{Gi}}-P_{\mathrm{Di}}=V_i\underset{j=1}{\overset{N}{\mathrm{\Σ}}}(G_{\mathrm{ij}}{V}_j\cos {\mathrm{\θ}}_{\mathrm{ij}}+B_{\mathrm{ij}}{V}_j\sin {\mathrm{\θ}}_{\mathrm{ij}})$

$Q_{\mathrm{Gi}}-Q_{\mathrm{Di}}=V_i\underset{j=1}{\overset{N}{\mathrm{\Σ}}}(G_{\mathrm{ij}}{V}_j\cos {\mathrm{\θ}}_{\mathrm{ij}}-B_{\mathrm{ij}}{V}_j\cos {\mathrm{\θ}}_{\mathrm{ij}})$

Meritorious and the idle units limits condition of power supply:

P _Gi,min≤P _Gi≤P _Gi,max

Q _Gi,min≤Q _Gi≤Q _Gi,max

Node voltage bound constraint condition:

V _imin≤V _i≤V _imax

Branch road rated current constraint condition, considered the inverse current that the distributed power source access causes:

$\left|I_k\right|\≤I_k^{\mathrm{cap}}$

Distributed power source planning total volume retrains, and is no more than the alpha proportion of system's total load:

$\underset{g=1}{\overset{n_{\mathrm{DG}}}{\mathrm{\Σ}}}{P}_{\mathrm{DG},g}\≤\mathrm{\α}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{P}_{\mathrm{Di}}$

Distributed power source delays income to the capacity of network and accounts for the big or small as follows of its total cost:

C _L≥βC _DG

C _DGBe the unit integrated cost of distributed power source, β is that income is to the contribution rate of cost.

Beneficial effect of the present invention has the following advantages: the present invention is based on power distribution network node limit Capacity Cost (LMCC-AC) computation model and the method for AC power flow, be not only applicable to distribution-free formula power supply situation, also be applicable to single distributed power source and many distributed power sources are incorporated into the power networks scene and need change power supply and electrical network parameter, and it is also conceivable that the reactive power effect of distributed power source.

Description of drawings

Fig. 1 is radial distribution feeder synoptic diagram (not containing DG);

Fig. 2 is the increase-volume time point τ that asks for a k road _k

Fig. 3 is radial distribution feeder synoptic diagram (containing DG);

Fig. 4 delays planning for the power distribution network increase-volume based on DG.

Embodiment

The present invention is described in further detail below in conjunction with embodiment and accompanying drawing, but embodiments of the present invention are not limited to this.

(1) power distribution network node/branch road limit Capacity Cost modeling

1) node sensitivity coefficient matrix

Radial distribution feeder synoptic diagram as shown in Figure 1 supposes that the growth pattern of load is the constant power factor increase, and the support level that the load of node i increases the electric current increase of branch road k adopts sensitivity coefficient γ _IkExpression:

${\mathrm{\γ}}_{\mathrm{ik}}=\frac{{\∂I}_k}{{\∂P}_{\mathrm{di}}}$

For radial power distribution network circuit, the load variations of node i is only influential to branch road (comprising circuit and the transformer branch road) electric current of supporting this load power supply, and therefore, sensitivity matrix γ is a sparse matrix.γ _IkNot only characterize the meritorious variation of node i, and included the reactive power fluctuation impact.Apparently, γ _IkRelevant with the network load distribution, γ is not normal matrix, can be expressed as: γ=γ _Ik(P _d); P in the formula _dCharacterize the network load distribution matrix.

2) branch current variable quantity

Load fluctuation on other node affects the electric current I of branch road k _kBe changed to:

N in the formula _kFor utilizing branch road k, need carry out the load bus collection of power delivery, dP _Di=δ _i(t) dt, wherein δ _i(t) be the load growth rate of node i in the different periods; Obtain the current increases amount of branch road k:

3) the branch road increase-volume time

By the curent change trend of branch road k, obtain branch road k increase-volume time point τ according to formula _k

${\∫}_0^{{\mathrm{\τ}}_k}\underset{i\∈N_k}{\mathrm{\Σ}}{\mathrm{\γ}}_{\mathrm{ik}}\left(P_d\right){\mathrm{\δ}}_i\left(t\right)\mathrm{dt}=I_k^{\mathrm{cap}}-I_k$

In the formula,

Be the rated capacity of branch road k, I _kBe the electric current under the time coordinate initial point (reference time).Calculate increase-volume time point τ based on AC power flow _kCan pass through the dichotomy method rapid solving, concrete calculation process as shown in Figure 2.

4) branch road increase-volume investment

The increase-volume cost of investment need be converted to present worth, and the employing formula is discounted calculating.

${\mathrm{PV}}_k=\frac{C_k}{{(1+r)}^{{\mathrm{\τ}}_k}}$

In the formula, C _kBe the increase-volume cost of branch road k, r is discount rate, PV _kDiscount value for the increase-volume cost.Also can adopt the discount computing formula of continuous compound rate form, namely

5) specific load increment cost matrix

51) node load changes

The load increment of supposing node i is △ P _Di, corresponding new network load distribution matrix is

The increase-volume time point that then branch road k is new By formula

Determine, in the formula,

For branch road k at the network load distribution matrix Under current value;

52) specific load increment cost coefficient

The increase-volume time point of the network equipment is not quite similar, and its specific load increment cost is also different, for ease of relatively and calculate, can adopt and wait the year value representation.Therefore, node i is to the specific load increment cost IC of branch road k _KiFor: ${\mathrm{IC}}_{\mathrm{ki}}=\frac{{\mathrm{\ΔPV}}_k}{\mathrm{\Δ}{S}_{\mathrm{di}}}\×P_{\mathrm{AP}}=\frac{{\mathrm{PV}}_k^{*}-{\mathrm{PV}}_k}{{\mathrm{\ΔP}}_{\mathrm{di}}/{\mathrm{pf}}_{\mathrm{di}}}\×R_{\mathrm{AP}}$

In the formula, pf _DiBe the load power factor of node i, R _APBe year value coefficients such as fund, relevant with discount rate with the increase-volume cost of investment receipts time limit; As △ P _Di→ 0 o'clock, IC _KiEquivalence is:

53) specific load increment cost matrix

Flow process in the same way, the arbitrary load bus of computing system is to the specific load increment cost IC of all branch roads _Ki, obtain specific load increment cost matrix IC and be: IC=(IC _Ki) _{B * N}In the formula, B is a way, and N is the load bus number.

6) node limit Capacity Cost

The marginal Capacity Cost of node i is that node i is supported the increment cost sum LMCC of unit of the branch road that this node power is carried to all _i:

B wherein _iFor by the set of fingers of load bus i utilization with through-put power.By aforementioned analysis as can be known, node i is zero to the specific load increment cost of not supporting the branch road that this node load is powered, so LMCC _iEquivalence is: ${\mathrm{LMGG}}_i=\underset{k\∈B}{\mathrm{\Σ}}{\mathrm{IC}}_{\mathrm{ki}}.$

7) branch road limit Capacity Cost

The marginal Capacity Cost of branch road k has utilized branch road k for all and has carried out the node of load power transmission to the increment cost sum BMCC of unit of this branch road _k,

With LMCC _iDiscussion similar, so BMCC _kCan equivalence be:

(2) distributed power source is incorporated into the power networks on the quantification of power distribution network node/branch road limit Capacity Cost impact

Behind the distributed power source access power distribution network, to its access node and on every side node load all play the power supporting role, distributed power source is regarded as a negative load, its access exerts an influence to the network load distribution situation, has changed load distribution matrix P _dWith specific load increment cost matrix IC, thus so that the marginal Capacity Cost LMCC of each node _iMarginal Capacity Cost BMCC with each branch road _kVariation has also occured, node limit Capacity Cost variable quantity C _LMCCFor: $C_{\mathrm{LMCC}}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}({\mathrm{LMCC}}_{\mathrm{DG},i}-{\mathrm{LMCC}}_i);$

LMCC in the formula _{DG, i}The be incorporated into the power networks marginal Capacity Cost of posterior nodal point i of expression distributed power source;

Branch road limit Capacity Cost variable quantity C _BMCCFor:

BMCC in the formula _{DG, k}The be incorporated into the power networks marginal Capacity Cost of rear branch road k of expression distributed power source; The marginal Capacity Cost LMCC of each node _iMarginal Capacity Cost BMCC with each branch road _kVariation be equivalent.Prove as follows:

$C_L=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}({\mathrm{LMCC}}_{\mathrm{DG},i}-{\mathrm{LMCC}}_i)=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}(\underset{k=1}{\overset{B}{\mathrm{\Σ}}}{\mathrm{IC}}_{\mathrm{DG},\mathrm{ki}}-\underset{k=1}{\overset{B}{\mathrm{\Σ}}}{\mathrm{IC}}_{\mathrm{ki}})$

$=\underset{k=1}{\overset{B}{\mathrm{\Σ}}}(\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{\mathrm{IC}}_{\mathrm{DG},\mathrm{ki}}-\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{\mathrm{IC}}_{\mathrm{ki}})=\underset{k=1}{\overset{B}{\mathrm{\Σ}}}({\mathrm{BMCC}}_{\mathrm{DG},k}-{\mathrm{BMCC}}_k)=C_B$

(3) performance evaluation of distributed power source configuration

Be incorporated into the power networks on the quantitative model of power distribution network node/branch road limit Capacity Cost impact based on distributed power source, the user provided for oneself carry out benefit with two kinds of all distributed power source property right patterns of grid company and find the solution and assess:

1) user provides the Capacity Benefit of distributed power source for oneself

Grid company is according to C _LMCCOr C _BMCCQuantize to assess the Capacity Benefit B that the user provides distributed power source for oneself _DG, grid company need compensate user distribution formula power construction expense B _DG:

B _DG=-S _DG×C _LMCC=-S _DG×C _BMCC

Wherein, S _DGRated capacity for distributed power source;

2) the power distribution network increase-volume based on distributed power source delays planning

Consider distributed power source to the increase-volume retarding action of power distribution network, the adoption rate factor is determined its income of bringing (system node limit Capacity Cost variable quantity C _LMCC) to the contribution rate of unit cost, mathematical model is:

$\min \underset{g=1}{\overset{n_{\mathrm{DG}}}{\mathrm{\Σ}}}{S}_{\mathrm{DG},g}$

In the formula, n _DGBe distributed power source configuration sum; Satisfy trend constraint condition:

$S.T.P_{\mathrm{Gi}}-P_{\mathrm{Di}}=V_i\underset{j=1}{\overset{N}{\mathrm{\Σ}}}(G_{\mathrm{ij}}{V}_j\cos {\mathrm{\θ}}_{\mathrm{ij}}+B_{\mathrm{ij}}{V}_j\sin {\mathrm{\θ}}_{\mathrm{ij}})$

$Q_{\mathrm{Gi}}-Q_{\mathrm{Di}}=V_i\underset{j=1}{\overset{N}{\mathrm{\Σ}}}(G_{\mathrm{ij}}{V}_j\cos {\mathrm{\θ}}_{\mathrm{ij}}-B_{\mathrm{ij}}{V}_j\cos {\mathrm{\θ}}_{\mathrm{ij}})$

Meritorious and the idle units limits condition of power supply:

P _Gi,min≤P _Gi≤P _Gi,max

Q _Gi,min≤Q _Gi≤Q _Gi,max

Node voltage bound constraint condition:

V _imin≤V _i≤V _imax

Branch road rated current constraint condition, considered the inverse current that the distributed power source access causes:

$\left|I_k\right|\≤I_k^{\mathrm{cap}}$

Distributed power source planning total volume retrains, and is no more than the alpha proportion of system's total load:

$\underset{g=1}{\overset{n_{\mathrm{DG}}}{\mathrm{\Σ}}}{P}_{\mathrm{DG},g}\≤\mathrm{\α}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{P}_{\mathrm{Di}}$

Distributed power source delays income to the capacity of network and accounts for the big or small as follows of its total cost:

C _L≥βC _DG

C _DGBe the unit integrated cost of distributed power source, β is that income is to the contribution rate of cost.

Based on heuristic thought, distributed power source priority access LMCC _iMaximum node adopts algorithm flow as shown in Figure 4 to carry out model solution.

Above-described embodiment is the better embodiment of the present invention; but embodiments of the present invention are not restricted to the described embodiments; other any do not deviate from change, the modification done under Spirit Essence of the present invention and the principle, substitutes, combination, simplify; all should be the substitute mode of equivalence, be included within protection scope of the present invention.

Claims (4)

1. distributed generation planning method is characterized in that may further comprise the steps:

1) introduces node sensitivity coefficient matrix and specific load increment cost matrix, set up power distribution network node/branch road limit Capacity Cost model;

2) having quantized distributed power source is incorporated into the power networks on the impact of power distribution network node/branch road limit Capacity Cost;

3) be incorporated into the power networks on the quantitative model of power distribution network node/branch road limit Capacity Cost impact based on distributed power source, the user provided for oneself carry out benefit with two kinds of distributed power source property rights of grid company pattern and find the solution and assess.

2. described distributed generation planning method according to claim 1 is characterized in that introducing node sensitivity coefficient matrix and the specific load increment cost matrix of described step 1), and the implementation of having set up power distribution network node/branch road limit Capacity Cost model is:

11) node sensitivity coefficient matrix

The support level that the load of node i increases the electric current increase of branch road k adopts sensitivity coefficient γ _IkExpression:

${\mathrm{\γ}}_{\mathrm{ik}}=\frac{{\∂I}_k}{{\∂P}_{\mathrm{di}}}$

I in the formula _kThe electric current that represents k bar branch road, P _DiThe load that represents i bar branch road;

Sensitivity matrix γ is a sparse matrix; Be expressed as: γ=γ _Ik(P _d); P in the formula _dCharacterize the network load distribution matrix;

12) branch current variable quantity

Load fluctuation on other node affects the curent change of branch road k, the electric current I of branch road k _kBe changed to: ${\mathrm{dI}}_k=\underset{i\∈N_k}{\mathrm{\Σ}}{\mathrm{\γ}}_{\mathrm{ik}}\left(P_d\right)d{P}_{\mathrm{di}};$

N in the formula _kFor utilizing branch road k, need carry out the load bus collection of power delivery, dP _Di=δ _i(t) dt, wherein δ _i(t) be the load growth rate of node i in the different periods;

The current increases amount of branch road k: ${\mathrm{dI}}_k=\underset{i\∈N_k}{\mathrm{\Σ}}{\mathrm{\γ}}_{\mathrm{ik}}\left(P_d\right){\mathrm{\δ}}_i\left(t\right)\mathrm{dt};$

13) the branch road increase-volume time

By the curent change trend of branch road k, obtain branch road k increase-volume time point τ according to following formula _k

${\∫}_0^{{\mathrm{\τ}}_k}\underset{i\∈N_k}{\mathrm{\Σ}}{\mathrm{\γ}}_{\mathrm{ik}}\left(P_d\right){\mathrm{\δ}}_i\left(t\right)\mathrm{dt}=I_k^{\mathrm{cap}}-I_k$

In the formula,

Be the rated capacity of branch road k, I _kBe the electric current under the time coordinate initial point;

14) branch road increase-volume investment

The increase-volume cost of investment need be converted to present worth, the employing formula

Or

The calculating of discounting; In the formula: C _kBe the increase-volume cost of branch road k, r is discount rate, PV _kDiscount value for the increase-volume cost;

15) specific load increment cost matrix

151) node load changes

The load increment of supposing node i is △ P _Di, corresponding new network load distribution matrix is

The increase-volume time point that then branch road k is new

By formula

Determine, in the formula,

For branch road k at the network load distribution matrix Under current value;

152) specific load increment cost coefficient

Node i is to the specific load increment cost IC of branch road k _KiFor:

In the formula, pf _DiBe the load power factor of node i, R _APBe year value coefficients such as fund, relevant with discount rate with the increase-volume cost of investment receipts time limit; As △ P _Di→ 0 o'clock, IC _KiEquivalence is:

153) specific load increment cost matrix

The arbitrary load bus of system is to the specific load increment cost IC of all branch roads _Ki, obtain specific load increment cost matrix IC and be: IC=(IC _Ki) _{B * N}In the formula, B is a way, and N is the load bus number;

16) node limit Capacity Cost

The marginal Capacity Cost of node i is that node i is supported the increment cost sum LMCC of unit of the branch road that this node power is carried to all _i:

17) branch road limit Capacity Cost

The marginal Capacity Cost of branch road k has utilized branch road k for all and has carried out the node of load power transmission to the increment cost sum BMCC of unit of this branch road _k:

3. described distributed generation planning method according to claim 1 is characterized in that described step 2) quantification distributed power source be incorporated into the power networks on the impact of power distribution network node/branch road limit Capacity Cost, its quantification manner is:

Behind the distributed power source access power distribution network, to its access node and on every side node load all play the power supporting role, distributed power source is regarded as a negative load, its access exerts an influence to the network load distribution situation, changed load distribution matrix Pd and specific load increment cost matrix IC, thus so that the marginal Capacity Cost LMCC of each node _iMarginal Capacity Cost BMCC with each branch road _kVariation has also occured, node limit Capacity Cost variable quantity C _LMCCFor: $C_{\mathrm{LMCC}}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}({\mathrm{LMCC}}_{\mathrm{DG},i}-\mathrm{LM}{\mathrm{CC}}_i);$

LMCC in the formula _{DG, i}The be incorporated into the power networks marginal Capacity Cost of posterior nodal point i of expression distributed power source;

Branch road limit Capacity Cost variable quantity C _BMCCFor:

BMCC in the formula _{DG, k}The be incorporated into the power networks marginal Capacity Cost of rear branch road k of expression distributed power source; The marginal Capacity Cost LMCC of each node _iMarginal Capacity Cost BMCC with each branch road _kVariation be equivalent.

4. described distributed generation planning method according to claim 1, it is characterized in that being incorporated into the power networks on the quantitative model of power distribution network node/branch road limit Capacity Cost impact based on distributed power source of described step 3), the user provided for oneself carry out benefit with two kinds of distributed power source property rights of grid company pattern and find the solution with the implementation of assessing and be:

31) user provides the Capacity Benefit of distributed power source for oneself

Grid company is according to C _LMCCOr C _BMCCQuantize to assess the Capacity Benefit B that the user provides distributed power source for oneself _DG, grid company need compensate user distribution formula power construction expense B _DG:

B _DG=-S _DG×C _LMCC=-S _DG×C _BMCC

Wherein, S _DGRated capacity for distributed power source;

32) the power distribution network increase-volume based on distributed power source delays planning

Consider distributed power source to the increase-volume retarding action of power distribution network, the adoption rate factor is determined its income of bringing to the contribution rate of unit cost, and mathematical model is:

$\min \underset{g=1}{\overset{n_{\mathrm{DG}}}{\mathrm{\Σ}}}{S}_{\mathrm{DG},g}$

In the formula, n _DGBe distributed power source configuration sum; Satisfy trend constraint condition:

$S.T.P_{\mathrm{Gi}}-P_{\mathrm{Di}}=V_i\underset{j=1}{\overset{N}{\mathrm{\Σ}}}(G_{\mathrm{ij}}{V}_j\cos {\mathrm{\θ}}_{\mathrm{ij}}+B_{\mathrm{ij}}{V}_j\sin {\mathrm{\θ}}_{\mathrm{ij}})$

$Q_{\mathrm{Gi}}-Q_{\mathrm{Di}}=V_i\underset{j=1}{\overset{N}{\mathrm{\Σ}}}(G_{\mathrm{ij}}{V}_j\cos {\mathrm{\θ}}_{\mathrm{ij}}-B_{\mathrm{ij}}{V}_j\cos {\mathrm{\θ}}_{\mathrm{ij}})$

Meritorious and the idle units limits condition of power supply:

P _Gi,min≤P _Gi≤P _Gi,max

Q _Gi,min≤Q _Gi≤Q _Gi,max

Node voltage bound constraint condition:

V _imin≤V _i≤V _imax

Branch road rated current constraint condition, considered the inverse current that the distributed power source access causes:

$\left|I_k\right|\≤I_k^{\mathrm{cap}}$

Distributed power source planning total volume retrains, and is no more than the alpha proportion of system's total load:

$\underset{g=1}{\overset{n_{\mathrm{DG}}}{\mathrm{\Σ}}}{P}_{\mathrm{DG},g}\≤\mathrm{\α}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{P}_{\mathrm{Di}}$

Distributed power source delays income to the capacity of network and accounts for the big or small as follows of its total cost:

C _L≥βC _DG

C _DGBe the unit integrated cost of distributed power source, β is that income is to the contribution rate of cost.

Images