CN105160450A - Energy-efficiency-included power plant power supply program method in consideration of environment factors - Google Patents

Abstract

The invention relates to an energy-efficiency-included power plant power supply program method in consideration of environment factors. The method comprises steps that, 1), an energy-efficiency-included power plant power supply program target function in consideration of environment factors is established, the target function aims to minimize investment cost, and the program investment cost comprises energy efficiency power plant cost, routine power plant set power generation cost, CO2 emission cost and SO2 emission cost; 2), constrained conditions of energy-efficiency-included power plant power supply program in consideration of the environment factors area acquired; 3), an energy-efficiency-included power plant power supply program model in consideration of the environment factors is established according to the target function and the constrained conditions; and 4), an optimal solution of the energy-efficiency-included power plant power supply program model in consideration of the environment factors is acquired through an interior point method, the program investment cost is minimized, and an establishment type and capacity of an energy-efficiency power plant are determined. Compared with the prior art, the method has advantages of advanced algorithm and comprehensive consideration.

Description

A kind of consider environmental factor containing energy efficiency power plant power source planning method

Technical field

The present invention relates to power source planning field, especially relate to a kind of consider environmental factor containing energy efficiency power plant power source planning method.

Background technology

Along with going from bad to worse of environmental pollution and being becoming tight day of energy supply, energy-saving and emission-reduction are worldwide subject to people's attention, and walking low-carbon (LC) path to economic development becomes the inevitable choice that the world realizes sustainable development.China is the country based on thermal power generation, maintaining sustained and rapid growth in recent years due to China's economy, some areas are made to have occurred the situation of a large amount of short of electricity, simple dependence enlarger pool-size and enlarging generating plant are difficult to the need for electricity meeting China, but also may bring the problems such as unit operation efficiency is low, primary energy consumption is large, environmental pollution is serious.In order to maintain the Sustainable Healthy Development of China's economy, clearly propose in " the comprehensive programme of work of " 12 " energy-saving and emission-reduction " that State Council issues to control energy resource consumption total amount and total emissions of major pollutants, and propose to strengthen demand Side Management, widely popularize energy efficiency power plant, as can be seen here, energy efficiency power plant construction is the important content of current China power construction.

Energy efficiency power plant refers to by adopting energy-efficient equipment and the approach such as product, optimization power mode, form the total action scheme of certain area, industry or enterprise's electrification rebuild plan, reach the object identical with Jian Xin power plant, the demand of minimizing is treated as the quantity of electricity that " virtual plant " provides, realize the object of energy conservation and pollution reduction.Compared with actual new power plant construction, energy efficiency power plant is carried out electric energy optimizing and obtains in former electric power system, and the no longer extra resource such as land occupation, consumption of fossil fuels, has huge Social benefit and economic benefit.

Traditional power source planning method can not consider the impact of environmental factor on programmed cost, and one of purport of energy efficiency power plant is exactly the discharge of decreasing pollution thing, along with environmental protection ideas affects all trades and professions gradually, therefore, in planning, consider that necessary Environmental costs become new requirement and emphasis.

Summary of the invention

Object of the present invention be exactly provide to overcome defect that above-mentioned prior art exists a kind of algorithm advanced, consider comprehensively to consider environmental factor containing energy efficiency power plant power source planning method.

Object of the present invention can be achieved through the following technical solutions:

That considers environmental factor contains an energy efficiency power plant power source planning method, comprises the following steps:

1) set up the objective function containing energy efficiency power plant power source planning considering environmental factor, this objective function minimizes for making planning cost of investment, and described planning cost of investment comprises energy efficiency power plant cost, conventional power plant unit generation cost, CO ₂discharge costs and SO ₂discharge costs;

2) constraint condition containing energy efficiency power plant power source planning considering environmental factor is obtained;

3) according to objective function and constraint condition set up consider environmental factor containing energy efficiency power plant Generation Expansion Planning Model;

4) adopt interior point method to obtain the optimum solution containing energy efficiency power plant Generation Expansion Planning Model of consideration environmental factor, make planning cost of investment minimum, thus determine the type that energy efficiency power plant is built and capacity.

Described step 1) in consideration environmental factor containing the objective function of energy efficiency power plant Generation Expansion Planning Model be:

$\min f=C_{\mathrm{CPP}}+C_{\mathrm{EPP}}+C_{{\mathrm{CO}}_2}+C_{{\mathrm{SO}}_2}$

$C_{CPP}=\underset{i=1}{\overset{N_C}{\Σ}}\frac{c_{{\mathrm{CPP}}_i}}{{\mathrm{\η}}_i}{P}_{{\mathrm{CPP}}_i}T$

$C_{EPP}=\underset{j=1}{\overset{N_E}{\Σ}}{c}_{{\mathrm{EPP}}_j}{P}_{{\mathrm{EPP}}_j}$

$C_{{\mathrm{CO}}_2}=\underset{i=1}{\overset{N_C}{\Σ}}{\mathrm{\α}}_C{\mathrm{\β}}_C{\mathrm{\λ}}_C{\mathrm{bP}}_{{\mathrm{CPP}}_i}{\mathrm{TC}}_C$

$C_{{\mathrm{SO}}_2}=\underset{i=1}{\overset{N_C}{\Σ}}{\mathrm{\α}}_S{\mathrm{\β}}_S{\mathrm{\λ}}_{S}b(1-{\mathrm{\ϵ}}_S)P_{{\mathrm{CPP}}_i}{\mathrm{TC}}_S$

Wherein, C _cPPfor the conventional power plant unit operation fuel cost of T hour, for the unit of fuel cost of conventional power plant unit i, for the capacity of conventional power plant unit i, N _cfor the quantity of unit in conventional power plant, η _ifor the av eff that conventional power plant unit i is annual, T is conventional power plant unit annual utilization hours, C _ePPfor the construction total cost of energy efficiency power plant, for the unit construction cost of energy efficiency power plant j, for the capacity of energy efficiency power plant j, N _efor the quantity of energy efficiency power plant, with be respectively the CO that conventional power plant runs T hour ₂and SO ₂discharge costs, α _cfor fossil fuel carbon content rate, α _sfor fossil fuel sulfur-bearing rate, β _cfor C is converted into CO ₂conversion coefficient, β _sfor S is converted into SO ₂conversion coefficient, λ _cfor the release rate of C in fossil fuel, λ _sfor the release rate of S in fossil fuel, C _cfor CO ₂discharge unit price, C _sfor SO ₂discharge unit price, b is gross coal consumption rate, ε _sfor the synthesis desulfurating efficiency of relevant device.

Described step 2) in constraint condition comprise equality constraint and inequality constrain condition.

Described equality constraint is Line Flow constraint, and Line Flow is constrained to active power flow and calculates, and calculating formula is:

P _C-P _L+P _E-Bθ＝0

Wherein, P _cfor conventional power plant capacity, P _lfor load, P _efor energy efficiency power plant capacity, B is system admittance matrix, and θ is node voltage phase angle.

Described inequality constrain condition comprises node voltage constraint, conventional power plant capacity-constrained, energy efficiency power plant capacity-constrained, CO ₂exhaust emission constraint, SO ₂exhaust emission constraint and energy efficiency power plant cost constraint,

Node voltage constraint inequality is:

0.9≤V _k≤1.1

Conventional power plant i capacity-constrained inequality is:

$\underset{\‾}{P_{{\mathrm{CPP}}_i}}\≤P_{{\mathrm{CPP}}_i}\≤\stackrel{\‾}{P_{{\mathrm{CPP}}_i}}$

Energy efficiency power plant j capacity-constrained inequality is:

$0\≤P_{{\mathrm{EPP}}_j}\≤\stackrel{\‾}{P_{{\mathrm{EPP}}_j}}$

CO ₂exhaust emission constraint constraint inequality is:

$\underset{i=1}{\overset{N_C}{\Σ}}{\mathrm{\α}}_C{\mathrm{\β}}_C{\mathrm{\λ}}_C{\mathrm{bP}}_{{\mathrm{CPP}}_i}T\≤\stackrel{\‾}{A_C}$

SO ₂exhaust emission constraint constraint inequality is:

$\underset{i=1}{\overset{N_C}{\Σ}}{\mathrm{\α}}_S{\mathrm{\β}}_S{\mathrm{\λ}}_{S}b(1-{\mathrm{\ϵ}}_S)P_{{\mathrm{CPP}}_i}T\≤\stackrel{\‾}{A_S}$

Energy efficiency power plant cost constraint inequality is:

$\underset{j=1}{\overset{N_E}{\Σ}}{C}_{{\mathrm{EPP}}_j}{P}_{{\mathrm{EPP}}_j}\≤\stackrel{\‾}{C_E}$

Wherein, V _kfor the voltage perunit value of a kth node, for conventional power plant lower bound of capacity, for conventional power plant maximum size, for the maximum size of energy efficiency power plant, for CO ₂the discharge upper limit, for SO ₂the discharge upper limit, for the energy efficiency power plant cost upper limit.

Compared with prior art, the present invention has the following advantages:

One, algorithm is advanced: the Generation Expansion Planning Model containing energy efficiency power plant of the present invention is a linear model, Primal-dual Interior Point is selected in model solution algorithm, the feature of primal-dual interior method is new for constrained objective function by structure---penalty is defined in feasible zone, and in feasible zone, ask the extreme point of penalty, exploration point when namely solving unconstrained problem is always inner at feasible zone, like this, in the process of sequence unconstrained optimization problem solving lnner guide, the solution always feasible solution of the serial unconstrained optimization problem tried to achieve, thus in the optimum solution of the former constrained optimization problem of the inner Step wise approximation of feasible zone, it is insensitive to problem scale that interior point method has iterations, convergence rapidly, strong robustness, computational accuracy is high, result and arranging of initial value contact the advantages such as not tight.

Two, consider comprehensively: the present invention set up consider environmental factor containing during energy efficiency power plant Generation Expansion Planning Model by CO ₂discharge costs and SO ₂discharge costs is taken into account and is comprehensively analyzed, and the result obtained can reflect the impact of environmental factor, has made positive effect to containing energy efficiency power plant power source planning.

Accompanying drawing explanation

Fig. 1 is invention process flow diagram of the present invention.

Fig. 2 is the process flow diagram of interior point method of the present invention.

Fig. 3 is final system grid structure figure.

Embodiment

Below in conjunction with the drawings and specific embodiments, the present invention is described in detail.

Embodiment:

As shown in Figure 1, a kind of consider environmental factor containing energy efficiency power plant power source planning method, comprise the following steps:

S1 sets up the objective function containing energy efficiency power plant power source planning considering environmental factor, and this objective function minimizes for making planning cost of investment, and planning cost of investment comprises energy efficiency power plant cost, conventional power plant unit generation cost, CO ₂discharge costs and SO ₂discharge costs;

S2 obtains the constraint condition containing energy efficiency power plant power source planning considering environmental factor;

S3 according to objective function and constraint condition set up consider environmental factor containing energy efficiency power plant Generation Expansion Planning Model;

S4 adopts interior point method to obtain the optimum solution containing energy efficiency power plant Generation Expansion Planning Model of consideration environmental factor, and make planning cost of investment minimum, the parameter of this optimum solution comprises conventional power plant capacity P _g, energy efficiency power plant capacity P _e, node voltage V and node voltage phase angle theta, thus determine the type that energy efficiency power plant is built and capacity.

Step 1) in consideration environmental factor containing the objective function of energy efficiency power plant Generation Expansion Planning Model be:

$\min f=C_{\mathrm{CPP}}+C_{\mathrm{EPP}}+C_{{\mathrm{CO}}_2}+C_{{\mathrm{SO}}_2}$

$C_{CPP}=\underset{i=1}{\overset{N_C}{\Σ}}\frac{c_{{\mathrm{CPP}}_i}}{{\mathrm{\η}}_i}{P}_{{\mathrm{CPP}}_i}T$

$C_{EPP}=\underset{j=1}{\overset{N_E}{\Σ}}{c}_{{\mathrm{EPP}}_j}{P}_{{\mathrm{EPP}}_j}$

$C_{{\mathrm{CO}}_2}=\underset{i=1}{\overset{N_C}{\Σ}}{\mathrm{\α}}_C{\mathrm{\β}}_C{\mathrm{\λ}}_C{\mathrm{bP}}_{{\mathrm{CPP}}_i}{\mathrm{TC}}_C$

$C_{{\mathrm{SO}}_2}=\underset{i=1}{\overset{N_C}{\Σ}}{\mathrm{\α}}_S{\mathrm{\β}}_S{\mathrm{\λ}}_{S}b(1-{\mathrm{\ϵ}}_S)P_{{\mathrm{CPP}}_i}{\mathrm{TC}}_S$

Wherein, C _cPPfor the conventional power plant unit operation fuel cost of T hour, for the unit of fuel cost of conventional power plant unit i, for the capacity of conventional power plant unit i, N _cfor the quantity of unit in conventional power plant, η _ifor the av eff that conventional power plant unit i is annual, T is conventional power plant unit annual utilization hours, C _ePPfor the construction total cost of energy efficiency power plant, for the unit construction cost of energy efficiency power plant j, for the capacity of energy efficiency power plant j, N _efor the quantity of energy efficiency power plant, with be respectively the CO that conventional power plant runs T hour ₂and SO ₂discharge costs, α _cfor fossil fuel carbon content rate, α _sfor fossil fuel sulfur-bearing rate, β _cfor C is converted into CO ₂conversion coefficient, β _sfor S is converted into SO ₂conversion coefficient, λ _cfor the release rate of C in fossil fuel, λ _sfor the release rate of S in fossil fuel, C _cfor CO ₂discharge unit price, C _sfor SO ₂discharge unit price, b is gross coal consumption rate, ε _sfor the synthesis desulfurating efficiency of relevant device.

Step 2) in constraint condition comprise equality constraint and inequality constrain condition.

Equality constraint is Line Flow constraint, and Line Flow is constrained to active power flow and calculates, and calculating formula is:

P _G-P _L+P _E-Bθ＝0

Wherein, P _gfor conventional power plant capacity, P _lfor load, P _efor energy efficiency power plant capacity, B is system admittance matrix, and θ is node voltage phase angle.

Described inequality constrain condition comprises node voltage constraint, conventional power plant capacity-constrained, energy efficiency power plant capacity-constrained, CO ₂exhaust emission constraint, SO ₂exhaust emission constraint and energy efficiency power plant cost constraint,

Node voltage constraint inequality is:

0.9≤V _k≤1.1

Conventional power plant i capacity-constrained inequality is:

$\underset{\‾}{P_{{\mathrm{CPP}}_i}}\≤P_{{\mathrm{CPP}}_i}\≤\stackrel{\‾}{P_{{\mathrm{CPP}}_i}}$

Energy efficiency power plant j capacity-constrained inequality is:

$0\≤P_{{\mathrm{EPP}}_j}\≤\stackrel{\‾}{P_{{\mathrm{EPP}}_j}}$

CO ₂exhaust emission constraint constraint inequality is:

$\underset{i=1}{\overset{N_C}{\Σ}}{\mathrm{\α}}_C{\mathrm{\β}}_C{\mathrm{\λ}}_C{\mathrm{bP}}_{{\mathrm{CPP}}_i}T\≤\stackrel{\‾}{A_C}$

SO ₂exhaust emission constraint constraint inequality is:

$\underset{i=1}{\overset{N_C}{\Σ}}{\mathrm{\α}}_S{\mathrm{\β}}_S{\mathrm{\λ}}_{S}b(1-{\mathrm{\ϵ}}_S)P_{{\mathrm{CPP}}_i}T\≤\stackrel{\‾}{A_S}$

Energy efficiency power plant cost constraint inequality is:

$\underset{j=1}{\overset{N_E}{\Σ}}{C}_{{\mathrm{EPP}}_j}{P}_{{\mathrm{EPP}}_j}\≤\stackrel{\‾}{C_E}$

Wherein, V _kfor the voltage perunit value of a kth node, for conventional power plant lower bound of capacity, for conventional power plant maximum size, for the maximum size of energy efficiency power plant, for CO ₂the discharge upper limit, for SO ₂the discharge upper limit, for the energy efficiency power plant cost upper limit.

Interior point method comprises the following steps:

1, obtain and meet the initial solution space of constraint condition, and following form can be reduced to what consider environmental factor containing energy efficiency power plant Generation Expansion Planning Model:

obj.min.f(x)

s.t.h(x)＝0

$\underset{\‾}{g}\≤g\left(x\right)\≤\stackrel{\‾}{g}$

Wherein, x=[P _c, P _e, V, θ], P _cfor conventional power plant capacity matrix, P _efor energy efficiency power plant capacity matrix, V is node voltage matrix, and θ is node phase angle matrix.

Meet u>0, l>0 by introducing slack variable u and l, u and l, inequality constrain be converted into equality constraint:

$g\left(x\right)+u=\stackrel{\‾}{g}$

g(x)-l＝ g

2, being Lagrangian function by model conversation is:

$L=f\left(x\right)-y^{T}h\left(x\right)-z^T\[g\left(x\right)-l-\underset{\‾}{g}\]-w^T\[g\left(x\right)+u-\stackrel{\‾}{g}\]-\mathrm{\μ}\underset{i=1}{\overset{r}{\Σ}}ln\left(l_i\right)-\mathrm{\μ}\underset{i=1}{\overset{r}{\Σ}}ln\left(u_i\right)$

In formula, y, z, w are Lagrange multiplier and meet y ≠ 0, z>0, w<0; μ is Discontinuous Factors, and definition duality gap Gap=l ^tz-u ^tw, then

3, solve above-mentioned disturbance KKT equation with Newton-Raphson method, obtain following system of equations:

Wherein, $L_x^{\&prime;}=L_x+{\&dtri;}_{x}g\left(x\right)\&lsqb;L^{-1}(L_l^{\mathrm{\&mu;}}+{\mathrm{ZL}}_z)+U^{-1}(L_u^{\mathrm{\&mu;}}-{\mathrm{WL}}_w)\&rsqb;$

$H^{\&prime;}=-\&lsqb;{\&dtri;}_x^{2}f\left(x\right)-{\&dtri;}_x^{2}h\left(x\right)y-{\&dtri;}_x^{2}g\left(x\right)(z+w)\&rsqb;-{\&dtri;}_{x}g\left(x\right)\&lsqb;L^{-1}Z-U^{-1}W\&rsqb;{\&dtri;}_x^{T}g\left(x\right)$

4, the correction that above-mentioned system of equations obtains kth time iteration is solved:

$\left[\begin{array}{c}{x}^{(k+1)}\\ l^{(k+1)}\\ u^{(k+1)}\end{array}\right]=\left[\begin{array}{c}{x}^{\left(k\right)}\\ l^{\left(k\right)}\\ u^{\left(k\right)}\end{array}\right]+{\mathrm{\&alpha;}}_p\left[\begin{array}{c}\mathrm{\&Delta;}x\\ \mathrm{\&Delta;}l\\ \mathrm{\&Delta;}u\end{array}\right]$

$\left[\begin{array}{c}{y}^{(k+1)}\\ z^{(k+1)}\\ w^{(k+1)}\end{array}\right]=\left[\begin{array}{c}{y}^{\left(k\right)}\\ z^{\left(k\right)}\\ w^{\left(k\right)}\end{array}\right]+{\mathrm{\&alpha;}}_d\left[\begin{array}{c}\mathrm{\&Delta;}y\\ \mathrm{\&Delta;}z\\ \mathrm{\&Delta;}w\end{array}\right]$

Wherein, α _pand α _dbe respectively original steps and antithesis step-length, and

${\mathrm{\&alpha;}}_p=0.9995min\{min(\frac{-l_i}{{\mathrm{\&Delta;l}}_i},{\mathrm{\&Delta;l}}_i<0;\frac{-u_i}{{\mathrm{\&Delta;u}}_i},{\mathrm{\&Delta;u}}_i<0),1\}(i=1,2...,r)$

Through successive ignition, when duality gap, Gap meets Gap<10 ^-6time, if calculate convergence, then the value of conventional power plant capacity variable, energy efficiency power plant capacity variable, node voltage variable and each variable of node phase angle variable in the optimal value that can obtain making cost of investment minimum and the corresponding initial solution space of this minimum value.The process flow diagram of interior point method as shown in Figure 2.

Example of the present invention adopts 18 node systems, carries out instance analysis.Final system rack as shown in Figure 3.This embodiment utilizes Primal-dual Interior Point, carries out iterative: first, inequality constrain is converted into equality constraint, and set up Lagrangian function to proposed Generation Expansion Planning Model; Then local derviation is asked to above-mentioned Lagrangian function, obtain optimal condition; Newton Algorithm is adopted to optimal condition, obtains update equation group; Solve above-mentioned update equation group, and solve step-length, upgrade original variable and Lagrange multiplier, after k iteration, obtain optimum solution.This inventive embodiment obtains optimum solution after 19 iteration, and result shows, after adding energy efficiency power plant, the cost of investment that power source planning is total in forcasted years is 4376137.8 ten thousand yuan, and wherein, energy efficiency power plant cost of investment is 373747.1 ten thousand yuan, conventional power plant operating cost is 2704521.8 ten thousand yuan, CO ₂discharge costs is 1558769.2 ten thousand yuan, SO ₂discharge costs is 239099.7 ten thousand yuan; By contrast, if do not add energy efficiency power plant, the cost of investment that power source planning is total in forcasted years is 4977769.7 ten thousand yuan, and wherein, conventional power plant operating cost is 2990075.3 ten thousand yuan, CO ₂discharge costs is 1723349.8 ten thousand yuan, SO ₂discharge costs is 264344.6 ten thousand yuan.

Claims (5)

1. that considers environmental factor contains an energy efficiency power plant power source planning method, it is characterized in that, comprises the following steps:

1) set up the objective function containing energy efficiency power plant power source planning considering environmental factor, this objective function minimizes for making planning cost of investment, and described planning cost of investment comprises energy efficiency power plant cost, conventional power plant unit generation cost, CO ₂discharge costs and SO ₂discharge costs;

2) constraint condition containing energy efficiency power plant power source planning considering environmental factor is obtained;

3) according to objective function and constraint condition set up consider environmental factor containing energy efficiency power plant Generation Expansion Planning Model;

4) adopt interior point method to obtain the optimum solution containing energy efficiency power plant Generation Expansion Planning Model of consideration environmental factor, make planning cost of investment minimum, thus determine the type that energy efficiency power plant is built and capacity.

2. according to claim 1 a kind of consider environmental factor containing energy efficiency power plant power source planning method, it is characterized in that, described step 1) in the objective function containing energy efficiency power plant Generation Expansion Planning Model of consideration environmental factor be:

$\min f=C_{CPP}+C_{EPP}+C_{{\mathrm{CO}}_2}+C_{{\mathrm{SO}}_2}$

$C_{CPP}=\underset{i=1}{\overset{N_C}{\&Sigma;}}\frac{c_{{\mathrm{CPP}}_i}}{{\mathrm{\&eta;}}_i}{P}_{{\mathrm{CPP}}_i}T$

$C_{EPP}=\underset{j=1}{\overset{N_E}{\&Sigma;}}{c}_{{\mathrm{EPP}}_j}{P}_{{\mathrm{EPP}}_j}$

$C_{{\mathrm{CO}}_2}=\underset{i=1}{\overset{N_C}{\&Sigma;}}{\mathrm{\&alpha;}}_C{\mathrm{\&beta;}}_C{\mathrm{\&lambda;}}_C{\mathrm{bP}}_{{\mathrm{CPP}}_i}{\mathrm{TC}}_C$

$C_{{\mathrm{SO}}_2}=\underset{i=1}{\overset{N_C}{\&Sigma;}}{\mathrm{\&alpha;}}_S{\mathrm{\&beta;}}_S{\mathrm{\&lambda;}}_{S}b(1-{\mathrm{\&epsiv;}}_S)P_{{\mathrm{CPP}}_i}{\mathrm{TC}}_S$

Wherein, C _cPPfor the conventional power plant unit operation fuel cost of T hour, for the unit of fuel cost of conventional power plant unit i, for the capacity of conventional power plant unit i, N _cfor the quantity of unit in conventional power plant, η _ifor the av eff that conventional power plant unit i is annual, T is conventional power plant unit annual utilization hours, C _ePPfor the construction total cost of energy efficiency power plant, for the unit construction cost of energy efficiency power plant j, for the capacity of energy efficiency power plant j, N _efor the quantity of energy efficiency power plant, with be respectively the CO that conventional power plant runs T hour ₂and SO ₂discharge costs, α _cfor fossil fuel carbon content rate, α _sfor fossil fuel sulfur-bearing rate, β _cfor C is converted into CO ₂conversion coefficient, β _sfor S is converted into SO ₂conversion coefficient, λ _cfor the release rate of C in fossil fuel, λ _sfor the release rate of S in fossil fuel, C _cfor CO ₂discharge unit price, C _sfor SO ₂discharge unit price, b is gross coal consumption rate, ε _sfor the synthesis desulfurating efficiency of relevant device.

3. according to claim 1 a kind of consider environmental factor containing energy efficiency power plant power source planning method, it is characterized in that, described step 2) in constraint condition comprise equality constraint and inequality constrain condition.

4. according to claim 3 a kind of consider environmental factor containing energy efficiency power plant power source planning method, it is characterized in that, described equality constraint is Line Flow constraint, and Line Flow is constrained to active power flow and calculates, and calculating formula is:

P _C-P _L+P _E-Bθ＝0

Wherein, P _cfor conventional power plant capacity, P _lfor load, P _efor energy efficiency power plant capacity, B is system admittance matrix, and θ is node voltage phase angle.

5. according to claim 3 a kind of consider environmental factor containing energy efficiency power plant power source planning method, it is characterized in that, described inequality constrain condition comprises node voltage constraint, conventional power plant capacity-constrained, energy efficiency power plant capacity-constrained, CO ₂exhaust emission constraint, SO ₂exhaust emission constraint and energy efficiency power plant cost constraint,

Node voltage constraint inequality is:

0.9≤V _k≤1.1

Conventional power plant i capacity-constrained inequality is:

$\underset{\&OverBar;}{P_{{\mathrm{CPP}}_i}}\&le;{P_{CPP}}_{{}_i}\&le;\stackrel{\&OverBar;}{P_{{\mathrm{CPP}}_i}}$

Energy efficiency power plant j capacity-constrained inequality is:

$0\&le;P_{{\mathrm{EPP}}_j}\&le;\stackrel{\&OverBar;}{P_{{\mathrm{EPP}}_j}}$

CO ₂exhaust emission constraint constraint inequality is:

$\underset{i=1}{\overset{N_C}{\&Sigma;}}{\mathrm{\&alpha;}}_C{\mathrm{\&beta;}}_C{\mathrm{\&lambda;}}_C{\mathrm{bP}}_{{\mathrm{CPP}}_i}T\&le;\stackrel{\&OverBar;}{A_C}$

SO ₂exhaust emission constraint constraint inequality is:

$\underset{i=1}{\overset{N_C}{\&Sigma;}}{\mathrm{\&alpha;}}_S{\mathrm{\&beta;}}_S{\mathrm{\&lambda;}}_{S}b(1-{\mathrm{\&epsiv;}}_S)P_{{\mathrm{CPP}}_i}T\&le;\stackrel{\&OverBar;}{A_S}$

Energy efficiency power plant cost constraint inequality is:

$\underset{j=1}{\overset{N_E}{\&Sigma;}}{C}_{{\mathrm{EPP}}_j}{P}_{{\mathrm{EPP}}_j}\&le;\stackrel{\&OverBar;}{C_E}$

Wherein, V _kfor the voltage perunit value of a kth node, for conventional power plant lower bound of capacity, for conventional power plant maximum size, for the maximum size of energy efficiency power plant, for CO ₂the discharge upper limit, for SO ₂the discharge upper limit, for the energy efficiency power plant cost upper limit.