Summary of the invention
The present invention mainly solves the technical problem existing for prior art; A kind of micro-capacitance sensor power source planning Scheme Optimum Seeking Methods taking into account dynamic characters is proposed first.First set up random process model to be described in the dynamic characters of micro-capacitance sensor power source planning, in economic index, add time-varying parameter based on Real Option Theory, propose energy prices relative risk index.And introduce a class power benefit index, the complementary benefit of multiple intermittent power supply cogeneration is described.Finally, use the weight of the economic index of analytic hierarchy process (AHP) clear and definite power source planning scheme, reliability index, complementary performance indicator, environmental protection index, overall merit is made with preferred to the power source planning scheme of the micro-capacitance sensor of different space-time characteristic.
Above-mentioned technical problem of the present invention is mainly solved by following technical proposals:
Take into account a micro-capacitance sensor power source planning Scheme Optimum Seeking Methods for dynamic characters, it is characterized in that,
Step 1: the programme treating selection sets up characterization technique progress and the probabilistic random process model of fuel price, if the starting year of certain power source planning scheme is t, by adjustment random process model, calculating simulation goes out the desired value I of t Construction of Unit cost
t, and programme starts the fuel cost desired value P of latter 1 year
t+n.
Step 2: use Monte Carlo Method of Stochastic, by the equilibrium data of this micro-capacitance sensor hour to calculate, then calculate year short of electricity probability and expect a LOLP, lack amount of power supply year and expect that LOLE, year short of electricity frequency expect that LOLF, mean hours short of electricity amount expect that EENS is as reliability evaluation index; Calculate the Financial cost C of the Life cycle of micro-capacitance sensor
cF, the operation of each power supply and fuel cost O
ti; Calculate energy prices risk cost C
risk.
Step 3: for the micro-grid system containing multiple intermittent energy electricity generation system, according to the result of production simulation, calculate the complementary performance indicator of multiple intermittent energy, as complementary gain capacity C
m, complementary gain degree λ
un, generation of electricity by new energy accounting K
n, abandon resource accounting K
a, energy storage device utilance K
s, and K is compared in the generating of environmental protection index high-carbon
carbonwith year pollutant discharge amount E
p.
Step 4: the power source planning scheme treating selection based on analytic hierarchy process (AHP) carries out overall merit, need to calculate its each refer to target value, then the index of different schemes is compared, is normalized, form an index coefficient matrix μ (x
w).
Step 5: policymaker, by after carrying out com-parison and analysis to the importance degree between two between index set, builds suitable index weights matrix and judgment matrix, calculates its eigenvalue of maximum and characteristic vector, determine rational weight matrix W
k.
Step 6: last, is calculated the comprehensive evaluation value T of micro-capacitance sensor power source planning scheme to be selected, selects cost performance preferred plan from all optional programs by linear weighted function summation.
In above-mentioned a kind of micro-capacitance sensor power source planning Scheme Optimum Seeking Methods taking into account dynamic characters, in described step 1, the desired value E [I of unit cost
t] along with the expression formula of the change of t be in time:
E[I
t]=I
0e
-λt(1-φ)=I
0e
-γt
In formula, I
0represent the Construction of Unit cost of t=0, if λ t recording technique innovation number of times, φ ∈ [0,1) be the constant that table technology levies innovation degree.Parameter lambda is technological innovation rate.In formula, γ=λ t (1-φ).
Then the fuel price desired value in certain moment following is:
E[P
t]=P
0e
μt
In formula, μ represents the drift rate of Brownian motion process, P
0represent the fuel cost of t=0.
In above-mentioned a kind of micro-capacitance sensor power source planning Scheme Optimum Seeking Methods taking into account dynamic characters, in described step 2, the Financial cost mathematic(al) representation of Life cycle is:
In formula, N represents power supply type number, x
ibe the number of i-th kind of power supply, it is initial from J that J represents the project life cycle,
be initial construction cost project period of i-th kind of power supply, C
ptibe the carbon emission punishment cost of i-th kind of power supply, T is the life cycle time limit, O
tibe operating cost during i-th kind of power supply T,
represent maintenance cost, r is discount rate.In formula, unit is at the construction cost I of t
tsee the Construction of Unit cost E [I in step 1
t].
Operation and the fuel cost expression formula of each power supply are:
In formula,
the fuel cost proportionality coefficient of each high-carbon energy, E
tibe the energy output of i-th kind of power supply at t.P
tfor the price of t in international crude petroleum, see the desired value E [P of the fuel price in step 1 at t
t].
Energy prices obey random process at the change procedure in future, and be the variable of a change at random, if do not take into account this change at the power source planning initial stage, then program results has deviation, and definition energy prices risk cost characterizes this deviation:
In formula, O
varduring for taking into account following annual energy price volatility, the fuel cost in the Life cycle of planned project, O
confor the fuel cost when supposition future source of energy price is constant.P
conassuming that constant energy prices, the general fuel cost by the planning starting year substitutes into and calculates.
In above-mentioned a kind of micro-capacitance sensor power source planning Scheme Optimum Seeking Methods taking into account dynamic characters, in described step 3, definition wind-light combined power generation system complementary gain capacity is
C
m=C
un-C
wind-C
PV
In formula, C
unthe credible capacity of this wind-light combined power generation system, C
windthe credible capacity of wind-powered electricity generation in system, C
pVit is the credible capacity of photovoltaic generation in system.
Definition complementary gain degree is complementary gain capacity C
mcapacity C credible with wind-light combined power generation system
unratio:
The ratio that definition generation of electricity by new energy amount accounts for year total capacity requirement is generation of electricity by new energy accounting:
In formula, ∑ E
nfor the year effective output of all new forms of energy units of micro-capacitance sensor, ∑ E
loadfor this total capacity requirement electricity in micro-capacitance sensor year.
It is that calculating formula is as follows by the ratio abandoning energy output and new forms of energy gross generation that resource accounting index is abandoned in definition:
In formula, ∑ E
nAfor meet situation due to load and energy storage does not remain active volume time, energy output that new forms of energy are abandoned, ∑ E
nAfor the year effective output of all new forms of energy units of micro-capacitance sensor.
The utilance of definition energy storage device is as follows:
In formula, ∑ E
exfor energy storage device be used for the accumulative annual energy output sent of support load time, ∑ E
easeyear for energy storage device adds up available energy output.
Therefore, tool of the present invention has the following advantages: propose a kind of micro-capacitance sensor power source planning Scheme Optimum Seeking Methods taking into account dynamic characters first.First set up random process model to be described in the dynamic characters of micro-capacitance sensor power source planning, in economic index, add time-varying parameter based on Real Option Theory, propose energy prices relative risk index.And introduce a class power benefit index, the complementary benefit of multiple intermittent power supply cogeneration is described.Finally, use the weight of the economic index of analytic hierarchy process (AHP) clear and definite power source planning scheme, reliability index, complementary performance indicator, environmental protection index, overall merit can be made with preferred to the power source planning scheme of the micro-capacitance sensor of different space-time characteristic.
Embodiment:
Below in conjunction with accompanying drawing, technical scheme of the present invention is described in further detail:
The present invention considers economic index, reliability index, performance indicator, environmental protection index carry out preferably (accompanying drawing 1) micro-capacitance sensor capital project.First the random process model of technological progress and fuel price in future time is established, in order to describe the dynamic space-time characteristic of micro-capacitance sensor investment.In the economic index of micro-capacitance sensor, add relevant time-varying parameter, and introduce a class power benefit index, describe the complementary benefit of multiple intermittent power supply cogeneration.Then use analytic hierarchy process (AHP) to carry out overall merit to above index, in the power source planning scheme never with the micro-capacitance sensor of space-time characteristic, get optimal case.
1. the economic index computational process taking into account dynamic characters is as follows:
Assuming that the starting year of certain power source planning project is t, improve in t production technology, this project need bear Construction of Unit expense is then I
t, then this part cost is the cost of investment of the whole project of technology then.Technological progress is the external cause that cost reduces, and the progress of technology over time process is random, assuming that during t > 0, the construction cost of unit is:
In formula, I
0represent the Construction of Unit cost of t=0, N
tpoisson stochastic variable, if λ t recording technique innovation number of times, φ ∈ [0,1) be the constant that table technology levies innovation degree.Parameter lambda is the rate of technological innovation.Obviously, the desired value I of unit cost
texponential relationship is had along with time t:
E[I
t]=I
0e
-λt(1-φ)=I
0e
-γt
In formula, γ=λ t (1-φ).Therefore, if adjustment parameter lambda or φ, reduce technology while the arrival rate λ as increase technological innovation and levy innovation degree φ, then variable λ (1-φ) may be constant.In addition, adjustment parameter lambda or φ, can change the probability distribution of following Construction of Unit cost, I
tvariance as follows:
As can be seen from the above equation, when λ levels off to infinity, φ level off to 1 time, variance Var [I
t] close to zero, namely increase λ and φ and the uncertainty of technological progress can be made to reduce.On the other hand, adjustment φ level off to 0 time, the uncertainty of technological progress is maximum.That is, can the degree of uncertainty of adjustment technology progress by adjustment λ and φ.
In micro-capacitance sensor, traditional high-carbon unit such as gas turbine, diesel engine need to rely on fossil fuel to generate electricity, and therefore the price of fuel directly has influence on the operating cost of conventional rack.In this article, the fuel cost P of project t
tbe assumed to meet geometric Brownian motion process in the change process in future:
In formula, μ and σ represents drift rate and the fluctuation parameters of this Brownian motion process respectively, and dz is the often step variation of standard Brownian motion process.
Assuming that 0≤μ < r, wherein r is discount rate.According to the mathematical characteristic of geometric Brownian motion process, then the fuel price desired value in certain moment following is:
E[P
t]=P
0e
μt
Then the uncertainty of fuel price can be adjusted by parameter σ, can be learnt by above formula, and the adjustment of parameter σ can not change the motion path of Brownian motion process.
The economic mathematical model of the Life cycle of micro-capacitance sensor is described below:
In formula, N represents power supply type number, x
ibe the number of i-th kind of power supply, it is initial from J that J represents the project life cycle,
be initial construction cost project period of i-th kind of power supply, C
ptibe the carbon emission punishment cost of i-th kind of power supply, T is the life cycle time limit, O
tibe operating cost during i-th kind of power supply T,
represent maintenance cost, r is discount rate.In formula, unit is at the construction cost I of t
tsee the Construction of Unit cost E [I in random process model used above
t].
To high-carbon energy, annual operation and fuel cost are proportional to International Crude Oil then.
In formula,
the fuel cost proportionality coefficient of each high-carbon energy, E
tibe the energy output of i-th kind of power supply at t.P
tfor the price of t in international crude petroleum, see the desired value E [P at t in fuel price random process model used above
t].
Because wind-resources and light resources are in time and complementarity geographically, make combining wind and light to generate electricity than being used alone wind power generation or photovoltaic generation more can make up due to the loss that brings of energy intermittence.Definition complementary gain capacity and complementary degree weigh this benefit herein.
2. performance indicator computational process is as follows:
If wind-light combined power generation system complementary gain capacity is:
C
m=C
un-C
wind-C
PV
In formula, C
unthe credible capacity of this wind-light combined power generation system, C
windthe credible capacity of wind-powered electricity generation in system, C
pVit is the credible capacity of photovoltaic generation in system.From definition, if complementary gain capacity C
mbe greater than 0, then system benefits is in complementary characteristic, C
mlarger, complementary characteristic is better.
Definition complementary gain degree is complementary gain capacity C
mcapacity C credible with wind-light combined power generation system
unratio:
λ in formula
unfor the complementary gain degree of wind-light combined power generation system.From definition, if complementary gain degree λ
unlarger, system complementary characteristic is better.
The ratio that definition generation of electricity by new energy amount accounts for year total capacity requirement is generation of electricity by new energy accounting:
In formula, ∑ E
nfor the year effective output of all new forms of energy units of micro-capacitance sensor, ∑ E
loadfor this total capacity requirement electricity in micro-capacitance sensor year.
If what the active volume of energy storage was less than this moment new energy source machine group is imbued with when exerting oneself, the situation of " abandoning wind " and " abandoning light " will be there will be.Resource accounting index is abandoned in definition is herein that calculating formula is as follows by the ratio abandoning energy output and new forms of energy gross generation:
In formula, ∑ E
nAfor meet situation due to load and energy storage does not remain active volume time, energy output that new forms of energy are abandoned, ∑ E
nAfor the year effective output of all new forms of energy units of micro-capacitance sensor.
The residue active volume of energy storage device along with distributed power source different with the match condition of load and fluctuate, the power that the charge-discharge electric power of energy storage device and energy storage and power supply or load exchange.The utilance of definition energy storage device is as follows:
In formula, ∑ E
exfor energy storage device be used for the accumulative annual energy output sent of support load time, ∑ E
easeyear for energy storage device adds up available energy output.
Energy prices obey random process at the change procedure in future, and be the variable of a change at random, if do not take into account this change at the power source planning initial stage, then program results has deviation, define energy prices risk cost herein and characterize this deviation:
In formula, O
varduring for taking into account following annual energy price volatility, the fuel cost in the Life cycle of planned project, O
confor the fuel cost when supposition future source of energy price is constant.P
conassuming that constant energy prices, the general fuel cost by the planning starting year substitutes into and calculates.
After having calculated the indices of all planning capital projects to be selected, Integrated comparative need be done based on the index of analytic hierarchy process (AHP) to different schemes.
3. the preferred computational process of scheme based on analytic hierarchy process (AHP) is as follows:
Sort to index according to importance degree, structure judges weight ratio matrix, and then Judgement Matricies B, and calculates its eigenvalue of maximum λ
max, and corresponding characteristic vector as:
x=[x
1,x
2,x
3,…x
n]
T
The consistency of test and judge matrix, the degree of consistency of metric matrix B is expressed as
as C (B)≤0.1, think that the compatibility of judgment matrix B is better, the eigenvalue of maximum λ of matrix B
maxcharacteristic of correspondence vector is x=[x
1, x
2, x
3... x
n]
tbe exactly weight vectors W=[w
1, w
2, w
3... w
n]
t.
Because characteristic value characteristic of correspondence vector is general not unique, so carry out the normalization of characteristic vector.Normalization formula is as follows, the factor of evaluation for being the bigger the better:
Factor of evaluation for the smaller the better:
Bring gained Combining weights into formula, calculated the comprehensive evaluation value of micro-capacitance sensor power source planning scheme by linear weighted function summation:
In formula: μ (x
w) be the normalized value of actual index, W
kfor waiting the weighted value asking each index;
4., in order to verify the beneficial effect of the inventive method, carried out following emulation experiment:
Certain the actual Island electrical network adding wind-light storage electricity generation system is adopted to be example.The distributed electrical Source Type studied herein has blower fan, photovoltaic cell, miniature gas turbine 3 kinds.
To with reference in document [20] the crew qiting scheme that 4 kinds meet this micro-capacitance sensor power supply reliability, in table 1.These 4 kinds of schemes have all made the year short of electricity of present case micro-capacitance sensor hour in required scope.
Table 1
The allocation plan (unit: MW) of table 1 distributed power source unit
New forms of energy installation wherein in scheme 1 is minimum, and high-carbon machine kludge is maximum, does not have storage battery; Scheme 2 is contrary, is an installation scheme being partial to new forms of energy; Scheme 3 and scheme 4 are two kinds of half-way houses in new forms of energy and traditional high-carbon energy.
Suppose that the planning starting year of this micro-capacitance sensor was from 2014, as shown in table 2 to the attribute decision table of 4 kinds of schemes:
Table 2
Table 2 scheme attribute decision table
Weight vector W corresponding to its index obtains by analytic hierarchy process (AHP), namely
W=(W
1,W
2,…w
20)=(0.0822,0.0822,0.1643,
0.0548,0.1643,0.0822,0.0235,0.0183,
0.0235,0.0235,0.0274,0.0183,0.0205,
0.0411,0.0235,0.0041,0.0088,0.0044,
0.0044,0.0044)
According to Weight summation, the index comprehensive assessed value finally obtaining 4 kinds of schemes is
T
1=0.451,T
2=0.520,T
3=0.789,T
4=0.467
Carry out sequence according to assessed value size to such scheme can obtain:
T
3>T
2>T
4>T
1
Result shows, in this island microgrid planning example, after having taken into account the dynamic characters of evaluation index, considers reliability, economy and environmental protection index, adopts scheme 3 best.
In order to the dynamic characters in micro-capacitance sensor power source planning scheme is described, suppose that micro-capacitance sensor planning is initial from the different times, the scheme 1 in his-and-hers watches 1 does not have storage battery relatively more extreme, and hereafter Main Analysis scheme 2 is to scheme 4.
Accompanying drawing 3 is the Life cycle costs from initial micro-capacitance sensor programme in the different following not the same year:
Simulation result shows, for the micro-capacitance sensor in this example, from 2014 during planning construction, adopts the Life cycle cost of power configuration scheme 2 the highest; But from 2022 when this locality planning micro-capacitance sensor, adopt the Life cycle cost of power configuration scheme 2 to economize most.
According to the prediction to WeiLai Technology Progress Rate and the price of international energy of american energy Information Management Bureau and Mo Te MacDonald company, space-time characteristic can be divided into low-carbon (LC) scene, balance scene and high-carbon scene.Carried out the degree of uncertainty of adjustment technology progress by adjustment λ and φ, adjustment parameter σ adjusts the uncertainty of fuel price, simulates respectively to balance scene and high-carbon scene.
Based on these two kinds of scenes, respectively in emulation table 1 the execution starting year of rear 3 kinds of programmes be respectively 2011,2012, until the situation of 2039.
Its balance scene emulation obtain future 3 kinds of programmes overall target changing trend diagram as shown in Figure 4:
High-carbon scene simulation obtains the overall target changing trend diagram of following 3 kinds of programmes as accompanying drawing 5.
Can be seen by above emulation, after taking into account multi-space characteristic, the evaluation index that the comprehensive evaluation index of scheme 2 exceedes scheme 3 is following trend.In balance scene, scheme 2 exceeded scheme 3 in 2021, and in high-carbon scene, scheme 2 is deferred to 2032 and exceedes scheme 3.This is because the investment construction cost of unit in high-carbon scene and fuel price have less uncertainty, and in scheme 2, generation of electricity by new energy accounting is higher, will benefit from this uncertainty.Scientific and technological progress and support on policy can accelerate high-carbon scene to balance scene conversion, make generation of electricity by new energy have more cost performance.
Specific embodiment described in the present invention is only to the explanation for example of the present invention's spirit.Those skilled in the art can make various amendment or supplement or adopt similar mode to substitute to described specific embodiment, but can't depart from spirit of the present invention or surmount the scope that appended claims defines.