CN103904696A - Capacity allocation optimization method for island/offshore platform independent new energy resource micro-grid - Google Patents

Abstract

The invention belongs to the field of electric system micro-grids and relates to a capacity allocation optimization method for an island/offshore platform independent new energy resource micro-grid. The method includes the steps that the annual load hour average power of the island/offshore platform is acquired through actual measurement or estimation; local annual natural condition hour average data are acquired from a meteorological department; the annual predicated generation power of each power generation device per unit capacity is calculated; the number of all power generation devices and the number of all energy storage devices are initialized; the gross generation of all the power generation devices and the total power consumption of loads in the time period are calculated; the annual load power shortage ratio and energy surplus ratio are calculated. According to the method, correction factors for influence of special natural conditions on system parameters and correction factors for influence of special geographical conditions on cost are adopted, the ocean power generation devices and diesel engines are additionally arranged, and therefore the micro-grid can provide electric power better through renewable energy sources and the power supply reliability of the system is improved as well.

Description

A kind of capacity configuration optimization method for the micro-electrical network of island/offshore platform independence new forms of energy

Technical field

The invention belongs to the micro-electrical network of electric power system field, relate to the most preferably access method of capacity configuration of island/ocean platform independent micro-grid system.

Background technology

Island independent micro-grid system is a kind of being applicable to away from the island on land or the electric power system of offshore platform, uses multiple distributed new generating equipment and traditional generating equipment that electric energy is provided.In system, the capacity of every distributed power generation equipment need to be selected according to local Natural resources condition and loading condiction with suiting measures to local conditions.But due to fluctuation and the uncertainty of natural energy resources, must add energy storage device with the guarantee system ability of operation continuously.Generally speaking in system, can comprise photovoltaic power generation equipment, wind-driven generator, power generation with marine energy equipment (comprising marine tidal-current energy, wave energy etc.), diesel engine, miniature gas turbine, energy storage device (as storage battery, flywheel energy storage, super capacitor energy-storage etc.), and relevant power-converting device.

The quality that so numerous its capacity of electrification energy storage equipment is chosen directly has influence on economy and the reliability of whole system, so a kind of outstanding Optimal Capacity choosing method is most important.The present invention, just for this factor, takes into full account the characteristic of ocean, island, has designed targetedly a kind of capacity configuration optimization method for the micro-electrical network of island/offshore platform independence new forms of energy.The result that method in application the present invention draws can be used as the foundation of the micro-electrical network master-plan of island/offshore platform.

In the domestic existing research of optimizing about independent micro-electrical network distributed power source capacity, taking economy, reliability, the quality of power supply etc. as target, choose single index or many indexs wind light mutual complementing, wind-light storage complementary power generation system are optimized to calculating.Patent " the stored energy capacitance choosing method of wind, light, storage micro-grid system " (application number 201010601680.4) has proposed a kind of method of optimizing specially stored energy capacitance in the situation that generating equipment capacity is definite.Patent " a kind of isolated microgrid distributed power source capacity and the optimization method of layouting " (application number 20120371966.7) has proposed a kind of use genetic algorithm, solves the optimized method of independent micro-net capacity.

But the research object of above research be all the micro-electrical network of the conventional independence in land, does not optimize targetedly for ocean, island micro-grid system, its method and result are also and be not suitable for the capacity optimization computation of the micro-electrical network in island.In addition due to the particularity of natural resources in ocean, the combination of conventional distributed power generation can not be given full play to the traditional characteristic of its complementarity and mutual benefit, Novel generating equipment add the selection that is inevitable.

Summary of the invention

The object of the present invention is to provide a kind of calculating more accurately for the capacity configuration optimization method of the micro-electrical network of island/offshore platform independence new forms of energy.

The object of the present invention is achieved like this:

(1) obtain island/offshore platform load hourly average power of the whole year, root is determined the maximum annual short of electricity rate DPSP of permission _maxallow annual energy surplus ratio REPG with maximum _max, the generating equipment type that selected island/offshore platform autonomous power supply system comprises, at least comprises diesel engine, photovoltaic power generation equipment, wind power plant, ocean current generating device, the one in batteries to store energy;

(2) obtain local annual natural conditions hourly average data from meteorological department, data comprise temperature, intensity of illumination, wind speed, ocean current speed;

(3), with reference to the model of generating equipment, calculate each generating equipment unit capacity prediction generated output of the whole year:

The output power model of photovoltaic power generation equipment:

$I_{\mathrm{pv}}=I_{\mathrm{sc}}\·\{1-C_1[\exp \left(\frac{V_{\mathrm{pv}}-\mathrm{\ΔV}}{C_2\·V_{\mathrm{oc}}}\right)-1\left]\right\}+\mathrm{\ΔI}$

P＝V _pv·I _pv

In formula: I _pVfor photovoltage model best operating point electric current under arbitrary condition, I _scfor photovoltage model short circuit current, V _pvfor photovoltage model optimum operation point voltage under arbitrary condition, V _ocfor photovoltaic battery panel open circuit voltage,

The output power model of wind power plant:

$P_{\mathrm{wt}}\left(v\right)=\left\{\begin{array}{c}a\·v^3-b\·P_R;v_{\mathrm{ci}}\≤v\≤v_r\\ P_R;v_r\≤v\≤v_{\mathrm{co}}\\ 0;\mathrm{otherwise}\end{array}\right.$

$a=\frac{P_R}{{v_r}^3-{v_{c_i}}^3},b=\frac{{v_{c_i}}^3}{{v_r}^3-{v_{c_i}}^3}$

In formula: P _rfor the rated power of blower fan, ν _rfor the rated wind speed of blower fan,

for the incision wind speed of blower fan, v _cofor the cut-out wind speed of blower fan;

The output power model of ocean current machine generating equipment:

$P_{\mathrm{oc}}\left(v\right)=\left\{\begin{array}{c}a\·v^3-b\·P_R;v_{\mathrm{ci}}\≤v\≤v_r\\ P_R;v_r\≤v\≤v_{\mathrm{co}}\\ 0;\mathrm{otherwise}\end{array}\right.$

$a=\frac{P_R}{{v_r}^3-{v_{c_i}}^3},b=\frac{{v_{c_i}}^3}{{v_r}^3-{v_{c_i}}^3}$

In formula: P _rfor the rated power of ocean current machine, ν _rfor the nominal flow rate of ocean current machine, for the incision flow velocity of ocean current machine, v _cofor the flow velocity that cuts out of ocean current machine;

(4) quantity of initial all generating equipments and energy storage device, initial number can be decided to be the minimum value of reality in can Application Range;

(5) establish time t=1, characterize section sometime;

(6) calculate interior all gross generations of generating equipments of this time period, the total electricity consumption of load, the relatively size of the two;

(7) if generating total amount is greater than electricity consumption total amount in this time period, calculates energy storage device charging and whether reach the energy storage upper limit, if reach the upper limit, calculate the energy residue in this time period, if do not reach the upper limit, calculate the SOC after energy storage device charging;

The computational methods of EPG and charging SOC:

EPG(t)＝E _Gen(t)-[E _L(t) _ηinv+(SOC _max-SOC(t-1))/η _B]

In formula: E _gen(t) gross energy producing for renewable energy power generation, E _l(t) gross energy consuming for load, η _invfor the efficiency of inverter, SOC _maxfor the energy storage device SOC upper limit, SOC (t-1) is the SOC value of previous moment energy storage device, η _bfor energy storage device charge efficiency

SOC(t)＝SOC(t-1)·(1-σ)+[E _Gen(t)-E _L(t)/η _inv]·η _B

In formula: SOC (t) is the state-of-charge of storage battery in the time of time t, and SOC (t-1) is the state-of-charge of storage battery in the time of time t-1, and σ is self-discharge rate;

(8) if generating total amount is less than electricity consumption total amount in this time period, calculates energy storage device electric discharge and whether reach energy storage lower limit, if reach lower limit, calculate the energy disappearance DPS in this time period, if do not reach lower limit, calculate the SOC after energy storage device electric discharge,

The computational methods of DPS and electric discharge SOC:

DPS(t)＝E _L(t)-[E _Gen(t)+SOC(t-1)-SOC _min]η _inv

In formula: SOC _minfor energy storage device SOC lower limit:

SOC(t)＝SOC(t-1)·(1-σ)-[E _L(t)/η _inv-E _Gen(t)]；

(9) in the time of t≤8760, return to step 6, t is from increasing 1 simultaneously; Otherwise continue to carry out following steps;

(10) calculate annual load short of electricity rate (DPSP) and energy surplus ratio (REPG):

$\mathrm{DPSP}=\underset{t=1}{\overset{T}{\mathrm{\Σ}}}\mathrm{DPS}\left(t\right)/\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{E}_L\left(t\right)$

$\mathrm{REPG}=\underset{t=1}{\overset{T}{\mathrm{\Σ}}}\mathrm{EPG}\left(t\right)/\underset{t=1}{\overset{T}{\mathrm{\Σ}}}{E}_L\left(t\right);$

(11) judge whether DPSP and REPG meet constraints, if met simultaneously, the configuration combination that record is now selected, photovoltaic capacity, fan capacity, trend machine capacity, energy storage device capacity; Otherwise step after continuing to carry out, constraints is:

DPSP≤DPSP _max

REPG≤REPG _max；

(12) if the quantity of generating equipment of now selecting and the capacity of energy storage device do not reach the maximum of reality in can Application Range, quantity is corresponding to increasing 1, returns to step 5 and carries out; Otherwise step after continuing to carry out;

(13) calculate and contrast overall life cycle cost LCC and the levelized cost of energy LUEC of the combination that all meets constraints, select as required the combination of economy the best, the final result for optimizing:

LCC＝C _initial+C _o&m+E _energy+R _replace+S _salvage

$\mathrm{LUEC}=\frac{(C_{\mathrm{initial}}+C_{\mathrm{replace}}+C_{O\&M})\×\mathrm{CRF}}{\underset{t=1}{\overset{8760}{\mathrm{\Σ}}}{E}_{\mathrm{gen}}\left(t\right)}$

In formula: C _initialfor the initial cost of system, C _{o & m}for the cost of equipment operation and maintenance, E _energyfor energy cost, annual all fuel costs, R _replacefor the cost of alternative electric weight in the lifetime of system phase, S _salvagefor residual value, the one's last year net value in the end of the year.

Beneficial effect of the present invention is:

The present invention is according on the basis of existing similar invention, for the special circumstances on island, ocean, has increased the modifying factor on system parameters impact on special natural condition; Increase the modifying factor of special geological surrounding to cost impact; Increase power generation with marine energy equipment (tidal current generator, wave energy generator etc.) and diesel engine, made it both can better utilize regenerative resource that electric power is provided, strengthened again the power supply reliability of system.More than improve and innovate the micro-net capacity of mixing that can make for ocean, island and plan that calculating is more accurate, more practical.

Brief description of the drawings

Fig. 1 is the structure chart of a kind of typical island independent distribution formula electric power system.

Fig. 2 is the flow chart of realizing a kind of mode that capacity distributes rationally.

Embodiment

Below in conjunction with accompanying drawing, the invention will be further described.

For above problem, the present invention proposes a kind of micro-net capacity optimization method for the unique geographical environment in ocean.Propose, for the improved distributed power generation model of marine environment, in the situation that meeting constraints, to select the result of economy optimum.Its result more meets ocean feature, more approaches actual conditions.The micro-grid generation system in island involved in the present invention, comprises photovoltaic cell, wind-driven generator, energy by ocean current generator, wave energy generator, diesel power generation equipment, batteries to store energy equipment, with and relevant power-converting device.Target of the present invention is to provide a kind of capacity configuration optimization method for island/offshore platform, under the condition of known natural data situation, taking economy as target, many kinds of parameters index is constraint, obtain the optimum capacity configuration of above various Blast Furnace Top Gas Recovery Turbine Unit (TRT) and energy storage device, and can obtain the indices such as economy, the fail safe data under different configuring conditions.

A kind of embodiment that stand alone generating system capacity optimization method is mixed on the island that the present invention proposes is as follows:

In step (1), establish DPSP _max=0.1, REPG _max=0.5.The generating equipment of selecting has photovoltaic power generation equipment, wind-driven generator, ocean current generator, lead-acid battery energy storage.

In step (2), obtain the data of the somewhere intensity of illumination of the whole year, wind speed, temperature, ocean current speed

In step (3), the particularity of considering the natural conditions of ocean, has proposed improved illumination, wind speed model, and added additional factor in the device model relevant to cost calculation, model is revised, more the cost of installation and use of equipment in ocean in closing to reality.

1, the correction of intensity of illumination

Globalradiation on photovoltaic battery panel inclined plane be the intensity of illumination sum of direct projection, sky diffuse reflection and the ground return of this period on inclined plane:

G _G,tilt＝G _dir,tilt+G _diff,tilt+G _refl,tilt

Above three parameters can be by the intensity of illumination G of this period on horizontal plane _{g, hor}calculate.Wherein, on inclined plane, the account form of direct projection and diffuse reflection intensity of illumination is identical with conventional method:

$G_{\mathrm{diff},\mathrm{tilt}}=\left[G_{\mathrm{diff},\mathrm{<figure-callout\; id="1"\; label="hor"\; filenames="HDA0000474494680000021.png"\; state="\{\{state\}\}">hor</figure-callout>}}\frac{1+\cos }{2}\right]\&CenterDot;[1+F\&times;{\sin }^3(\mathrm{\&beta;}/2)]\&CenterDot;[1+F\&times;{\cos }^2\mathrm{\&theta;}\&CenterDot;{\cos }^3{\mathrm{\&theta;}}_z]$

$G_{\mathrm{dir},\mathrm{tilt}}=G_{\mathrm{dir},\mathrm{hor}}\frac{\cos \mathrm{\&theta;}}{\cos {\mathrm{\&theta;}}_z}$

Ground return intensity of illumination G on inclined plane _{refl, tilt}with overall intensity of illumination G on horizontal plane _{g, hor}relevant.Due to the particularity of ocean condition, marine illumination reflection is stronger, when calculating, need in traditional calculating formula, add modifying factor α:

$G_{\mathrm{refl},\mathrm{tilt}}=\mathrm{\&rho;}\frac{G_{G,\mathrm{hor}(1-\cos )}}{2}\&CenterDot;\mathrm{\&alpha;}$

Its value is for being greater than 1, and relevant with photovoltaic apparatus installation ground periphery ocean surface area, sea is more, and its value is larger.

2, the correction of wind speed

The correction of wind speed is that a wind speed size conversion of measuring on anemometer height for handle is the air speed data at the firm setting height(from bottom) of wind power generation wheel place.Like this, the energy output that calculates wind energy conversion system output is just more accurate, the variation of wind speed and setting height(from bottom) have certain relation, generally highly larger, air speed value is also just larger, and the size of the size of wind speed and the roughness on earth's surface and temperature also has very large relation, in 500 meters of altitude ranges overhead, the variation of wind speed in vertical height, can press facial index formula and calculate:

$\frac{v}{v_0}={\left(\frac{H}{H_0}\right)}^{\mathrm{\&alpha;}}$

α is earth's surface coefficient of friction, is decided by the stability of atmosphere and the roughness on earth's surface, and under ocean condition, its value is generally 0.1～0.13.

3, the correction of atmospheric density

The size of atmospheric density is along with the temperature of the temperature of air, steam and the variation of height above sea level and different, this just need to revise accordingly to the atmospheric density of concrete infield, due to Island topographic structure more complicated, air humidity is also larger, therefore the atmospheric density on island is influential to the energy output of wind-driven generator, so be necessary atmospheric density to revise, to being modified to of atmospheric density:

$\mathrm{\&rho;}={\mathrm{\&rho;}}_0\&times;\frac{1}{1+0.00366t}\&times;\left(\frac{p-0.378e}{1000}\right)$

In the time there is no moisture measurement, use the Ideal-Gas Equation to solve density, computing formula is as follows:

$\mathrm{\&rho;}=\frac{p}{\mathrm{RT}}$

4, the correction of equipment cost

Cost correction herein comprises installation cost and working service cost.In this natural conditions situation in ocean, while calculating the equipment cost of electricity generation system, not only to consider purchase cost and the installation cost of equipment, also need to, for the special natural environment in ocean, consider the maintenance cost that it is annual.

In the present invention, for installation cost, because the place to use of equipment is away from land, its cost of transportation is relevant with the value of distance apart from land and equipment itself.In installation, consider to need the feature of the anti-high wind of the anti-salt sand prevention of waterproof, equipment all to need through special processing (as three anti-paints etc.).

In the present invention, in the economic model for photovoltaic power generation equipment, wind power plant, add system maintenance cost element, and set up its Mathematical Modeling:

C _(O&M)0＝ζC _i0

In step (4), determine that the reality of equipment can select the scope of quantity, photovoltaic power generation equipment is 0-20, and wind power plant is 0-10, and ocean current generating device is 0-10.

In step (5), initial time is made as t=1.

In step (6), calculate whole gross generations of generating equipment and the total electricity consumption of load in this time period.The relatively size of the two.

In step (7), the energy storage upper limit of energy storage device is to be determined by the total capacity of energy storage device.

In step (8), if generating total amount is less than electricity consumption total amount in this time period, calculates energy storage device electric discharge and whether reach energy storage lower limit.If reach lower limit, calculate the energy disappearance (DPS) in this time period, if do not reach lower limit, calculate the SOC after energy storage device electric discharge.

In step (9), in the time of t > 8760, represent that annual data are complete as calculated.

(10) calculate annual load short of electricity rate (DPSP) and energy surplus ratio (REPG):

$\mathrm{DPSP}=\underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}\mathrm{DPS}\left(t\right)/\underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}{E}_L\left(t\right)$

$\mathrm{REPG}=\underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}\mathrm{EPG}\left(t\right)/\underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}{E}_L\left(t\right);$

(11) judge whether DPSP and REPG meet constraints, if met simultaneously, the configuration combination (photovoltaic capacity, fan capacity, trend machine capacity, energy storage device capacity) that record is now selected; Otherwise step after continuing to carry out.Constraints is:

DPSP≤DPSP _max

REPG≤REPG _max

(12) if the quantity of generating equipment of now selecting and the capacity of energy storage device do not reach the maximum of reality in can Application Range, quantity is corresponding to increasing 1, returns to step 5 and carries out; Otherwise step after continuing to carry out.

(13) calculating and contrast overall life cycle cost (LCC) and the levelized cost of energy (LUEC) of the combination that all meets constraints, select as required the combination of economy the best, is the final result of optimizing.

LCC＝C _initial+C _o&m+E _energy+R _replace+S _salvage

$\mathrm{LUEC}=\frac{(C_{\mathrm{initial}}+C_{\mathrm{replace}}+C_{O\&M})\&times;\mathrm{CRF}}{\underset{t=1}{\overset{8760}{\mathrm{\&Sigma;}}}{E}_{\mathrm{gen}}\left(t\right)}.$

Claims (1)

1. for a capacity configuration optimization method for the micro-electrical network of island/offshore platform independence new forms of energy, it is characterized in that:

(1) obtain island/offshore platform load hourly average power of the whole year, root is determined the maximum annual short of electricity rate DPSP of permission _maxallow annual energy surplus ratio REPG with maximum _max, the generating equipment type that selected island/offshore platform autonomous power supply system comprises, at least comprises diesel engine, photovoltaic power generation equipment, wind power plant, ocean current generating device, the one in batteries to store energy;

(2) obtain local annual natural conditions hourly average data from meteorological department, data comprise temperature, intensity of illumination, wind speed, ocean current speed;

(3), with reference to the model of generating equipment, calculate each generating equipment unit capacity prediction generated output of the whole year:

The output power model of photovoltaic power generation equipment:

$I_{\mathrm{pv}}=I_{\mathrm{sc}}\&CenterDot;\{1-C_1[\exp \left(\frac{V_{\mathrm{pv}}-\mathrm{\&Delta;V}}{C_2\&CenterDot;V_{\mathrm{oc}}}\right)-1\left]\right\}+\mathrm{\&Delta;I}$

P＝V _pv·I _pv

In formula: I _pVfor photovoltage model best operating point electric current under arbitrary condition, I _scfor photovoltage model short circuit current, V _pvfor photovoltage model optimum operation point voltage under arbitrary condition, V _ocfor photovoltaic battery panel open circuit voltage,

The output power model of wind power plant:

$P_{\mathrm{wt}}\left(v\right)=\left\{\begin{array}{c}a\&CenterDot;v^3-b\&CenterDot;P_R;v_{\mathrm{ci}}\&le;v\&le;v_r\\ P_R;v_r\&le;v\&le;v_{\mathrm{co}}\\ 0;\mathrm{otherwise}\end{array}\right.$

$a=\frac{P_R}{{v_r}^3-{v_{c_i}}^3},b=\frac{{v_{c_i}}^3}{{v_r}^3-{v_{c_i}}^3}$

In formula: P _rfor the rated power of blower fan, ν _rfor the rated wind speed of blower fan,

for the incision wind speed of blower fan, v _cofor the cut-out wind speed of blower fan;

The output power model of ocean current machine generating equipment:

$P_{\mathrm{oc}}\left(v\right)=\left\{\begin{array}{c}a\&CenterDot;v^3-b\&CenterDot;P_R;v_{\mathrm{ci}}\&le;v\&le;v_r\\ P_R;v_r\&le;v\&le;v_{\mathrm{co}}\\ 0;\mathrm{otherwise}\end{array}\right.$

$a=\frac{P_R}{{v_r}^3-{v_{c_i}}^3},b=\frac{{v_{c_i}}^3}{{v_r}^3-{v_{c_i}}^3}$

In formula: P _rfor the rated power of ocean current machine, ν _rfor the nominal flow rate of ocean current machine,

for the incision flow velocity of ocean current machine, v _cofor the flow velocity that cuts out of ocean current machine;

(4) quantity of initial all generating equipments and energy storage device, initial number can be decided to be the minimum value of reality in can Application Range;

(5) establish time t=1, characterize section sometime;

(6) calculate interior all gross generations of generating equipments of this time period, the total electricity consumption of load, the relatively size of the two;

(7) if generating total amount is greater than electricity consumption total amount in this time period, calculates energy storage device charging and whether reach the energy storage upper limit, if reach the upper limit, calculate the energy residue in this time period, if do not reach the upper limit, calculate the SOC after energy storage device charging;

The computational methods of EPG and charging SOC:

EPG(t)＝E _Gen(t)-[E _L(t)η _inv+(SOC _max-SOC(t-1))/η _B]

In formula: E _gen(t) gross energy producing for renewable energy power generation, E _l(t) gross energy consuming for load, η _invfor the efficiency of inverter, SOC _maxfor the energy storage device SOC upper limit, SOC (t-1) is the SOC value of previous moment energy storage device, η _bfor energy storage device charge efficiency

SOC(t)＝SOC(t-1)·(1-σ)+[E _Gen(t)-E _L(t)/η _inv]·η _B

In formula: SOC (t) is the state-of-charge of storage battery in the time of time t, and SOC (t-1) is the state-of-charge of storage battery in the time of time t-1, and σ is self-discharge rate;

(8) if generating total amount is less than electricity consumption total amount in this time period, calculates energy storage device electric discharge and whether reach energy storage lower limit, if reach lower limit, calculate the energy disappearance DPS in this time period, if do not reach lower limit, calculate the SOC after energy storage device electric discharge,

The computational methods of DPS and electric discharge SOC:

DPS(t)＝E _L(t)-[E _Gen(t)+SOC(t-1)-SOC _min]η _inv

In formula: SOC _minfor energy storage device SOC lower limit:

SOC(t)＝SOC(t-1)·(1-σ)-[E _L(t)/η _inv-E _Gen(t)]；

(9) in the time of t≤8760, return to step 6, t is from increasing 1 simultaneously; Otherwise continue to carry out following steps;

(10) calculate annual load short of electricity rate (DPSP) and energy surplus ratio (REPG):

$\mathrm{DPSP}=\underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}\mathrm{DPS}\left(t\right)/\underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}{E}_L\left(t\right)$

$\mathrm{REPG}=\underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}\mathrm{EPG}\left(t\right)/\underset{t=1}{\overset{T}{\mathrm{\&Sigma;}}}{E}_L\left(t\right);$

(11) judge whether DPSP and REPG meet constraints, if met simultaneously, the configuration combination (photovoltaic capacity, fan capacity, trend machine capacity, energy storage device capacity) that record is now selected; Otherwise step after continuing to carry out, constraints is:

DPSP≤DPSP _max

REPG≤REPG _max；

(12) if the quantity of generating equipment of now selecting and the capacity of energy storage device do not reach the maximum of reality in can Application Range, quantity is corresponding to increasing 1, returns to step 5 and carries out; Otherwise step after continuing to carry out;

(13) calculate and contrast overall life cycle cost LCC and the levelized cost of energy LUEC of the combination that all meets constraints, select as required the combination of economy the best, the final result for optimizing:

LCC＝C _initial+C _o&m+E _energy+R _replace+S _salvage

$\mathrm{LUEC}=\frac{(C_{\mathrm{initial}}+C_{\mathrm{replace}}+C_{O\&M})\&times;\mathrm{CRF}}{\underset{t=1}{\overset{8760}{\mathrm{\&Sigma;}}}{E}_{\mathrm{gen}}\left(t\right)}$

In formula: C _initialfor the initial cost of system, C _{o & m}for the cost of equipment operation and maintenance, E _energyfor energy cost, annual all fuel costs, R _replacefor the cost of alternative electric weight in the lifetime of system phase, S _salvagefor residual value, the one's last year net value in the end of the year.

Images