Summary of the invention
The invention provides a kind of demand response effect evaluation method based on Fuzzy Analysis, base oneself upon electrical network, social demand, using reliable, economy, environmental protection as the three norms of weighing demand response effect, and guarantee comparability and the cumulative property between index by unified computing method, finally with Field Using Fuzzy Comprehensive Assessment, solve.
The technical solution adopted in the present invention is as follows:
A demand response effect evaluation method based on Fuzzy Analysis, comprises the following steps:
1) determine demand response dispatching effect evaluation index, described index comprises reliability index, economic index and feature of environmental protection index;
2) determine demand response effect assessment grade and pass judgment on collection;
Demand response effect assessment grade is divided into 6 grades from high to low, is defined as and passes judgment on collection, by described step 1) each index set respectively 6 quantized intervals, 6 quantized intervals of described 6 grades and each index are corresponding one by one;
3) determine demand response effect assessment set of factors;
4) determine the judgment matrix R of single index set of factors
k(k=1,2,3);
5) determine the weight coefficient matrix of factor of evaluation collection;
6) the fuzzy evaluation result of computation requirement response scheduling effect.
Aforesaid step 1), in, described reliability index comprises response mode reliability index R
type, response time reliability index R
timewith response Capacity Reliability index R
cap, wherein:
Response mode reliability index R
typerefer to that the number of users by the scheduling of automatic response mode participation demand response accounts for the ratio of all participating user numbers, computing formula is suc as formula shown in (1):
Wherein, n is for participating in the total number of users of demand response scheduling; R
i(con) being i user's control sign, is 0-1 decision variable, when control mode is automatic response mode, and R
i(con)=1, otherwise R
i(con)=0;
Response time reliability index R
timerefer to that actual schedule capacity surpasses the time scale of desirable scheduling capacity 80%;
Response Capacity Reliability index R
caprefer at scheduling capacity and be greater than within the scope of 80% desirable scheduling capacity, the product of weight coefficient and scheduling capacity ratio.
Aforesaid step 1), in, described economic index comprises unit response Capacity Cost E
unitwith the electrical network peak-valley difference E reducing
gap, wherein:
Unit response Capacity Cost E
unitrefer to that under certain event, grid company is the cost that Demand-side unit response capacity is paid, computing formula is as follows,
Wherein, n is for participating in the total number of users of demand response scheduling, C
i(t) be the response capacity of i user t period;
response total amount for i user in whole dispatching cycle; p
ithe economic compensation giving for electrical network under unit quantity of electricity;
response total volume for Demand-side generation;
The electrical network peak-valley difference E reducing
gaprefer within the demand response event phase, because of the ratio that demand response event causes electrical network peak-valley difference to reduce, computing formula is as follows,
Wherein, P
for(max), P
for(min) electrical network maximum, the minimum load of prediction while not occurring for event, P (max), P (min) is that the electrical network of event in the phase is maximum, minimum load.
Aforesaid step 1), in, described feature of environmental protection index comprises the environmental pollutant discharge B of minimizing
c,Nwith user's energy consumption grade B
level, wherein:
The environmental pollutant discharge B reducing
c,Nrefer to the carbon monoxide, carbon dioxide, nitrogen carbonoxide, sulphuric dioxide, particle and the hydrocarbons discharge capacity that because carrying out demand response event, reduce, computing formula is as follows,
Wherein, n is for participating in the total number of users of demand response scheduling, α
j(j=1,2 ..., 6) and be pollutant discharge amount design factor, α
1for CO emission design factor, α
2cO2 emissions design factor, α
3nitrogen carbonoxide Emission amount calculation coefficient, α
4sO2 emissions design factor, α
5discharge quantity of particulate matter design factor, α
6hydrocarbons Emission amount calculation coefficient,
for response total volume;
User's energy consumption grade B
levelwith reference to country, provinces and cities, electrical network, put into effect energy consumption index appraisal standards, be divided into 1 grade of energy consumption, 2 grades of energy consumptions and 3 grades of energy consumptions, energy consumption level raises successively.
Aforesaid step 3), in, described response effectiveness evaluation set of factors is divided into two-layer:
The comprehensive evaluation factor collection on upper strata is: U={u
1, u
2, u
3}={ reliability index, economic index, feature of environmental protection index },
The single index set of factors of bottom is respectively: U
1={ u
11, u
12, u
13}={ response mode reliability index, response time reliability index, response Capacity Reliability index }, U
2={ u
21, u
22}={ unit response Capacity Cost, the electrical network peak-valley difference of minimizing }, U
3={ u
31, u
32the environmental pollutant discharge of }={ minimizing, user's energy consumption grade }.
Aforesaid step 4), in, judgment matrix is respectively:
The judgment matrix R of reliability index
1for:
The judgment matrix R of economic index
2for:
The judgment matrix R of feature of environmental protection index
3for:
Wherein, certain single factor of line display single index set of factors, the share that the opinion rating of this monofactorial evaluation collection is shared is shown in list.
Aforesaid step 5), in, the weight coefficient matrix of factor of evaluation collection is respectively:
The heavy matrix of coefficients of the described comprehensive evaluation factor centralization of state power is: A=|a
1a
2a
3|, and a
1+ a
2+ a
3=1;
Described single index set of factors weight coefficient matrix is: A
1=| a
11a
12a
13|, and a
11+ a
12+ a
13=1, A
2=| a
21a
22|, and a
21+ a
22=1, A
3=| a
31a
32|, and a
31+ a
32=1,
Wherein, a
1, a
2, a
3for each factor in comprehensive evaluation factor collection U is passed judgment on the weight coefficient of collection, a relatively
11, a
12, a
13for single index set of factors U
1in each factor relatively pass judgment on the weight coefficient of collection, a
21, a
22for single index set of factors U
2in each factor relatively pass judgment on the weight coefficient of collection, a
31, a
32for single index set of factors U
3in each factor relatively pass judgment on the weight coefficient of collection.
Aforesaid step 6), the fuzzy evaluation result of computation requirement response scheduling effect, comprises the following steps:
6-1) with compound operation, obtain single index set of factors U
kthe comprehensive evaluation result B of (k=1,2,3)
k:
B
k=A
k*R
k=[b
k1b
k2…b
k6];
6-2) first solve the judgment matrix R of comprehensive evaluation factor collection U:
Then obtain the final comprehensive evaluation result B:B=A*R=[b of comprehensive evaluation factor collection U
1b
2b
6]
Wherein, b
1, b
2, b
3, b
4, b
5and b
6represent 6 shares that evaluation rank is shared;
6-3) choose b
1, b
2, b
3, b
4, b
5and b
6middle greatest measure is as this demand response dispatching effect evaluation result.
Compared with prior art, tool of the present invention has the following advantages:
The present invention considers the multiclass factors such as safety, reliable, economy, environmental protection, concrete index account form has been proposed, in relative value mode, guarantee comparability and the cumulative property of index, desired data convenient sources, easy to understand, and solve by Field Using Fuzzy Comprehensive Assessment, can take into account different electrical network targets, there is application wide, the feature that applicability is strong.
Embodiment
Below in conjunction with the drawings and specific embodiments, the present invention is described in further detail.
As shown in Figure 5, the demand response effect evaluation method based on Fuzzy Analysis of the present invention, comprises the following steps:
1, determine demand response dispatching effect evaluation index, base oneself upon electrical network demand, gap between evaluation requirement response events and the re-set target of this event, and the actual of this event tells on, determine that demand response dispatching effect evaluation index comprises reliability index, economic index and feature of environmental protection index three major types index altogether, as shown in Figure 1, wherein
1.1, reliability index comprises response time reliability index R
time, response mode reliability index R
typewith response Capacity Reliability index R
cap, particularly,
(1) response mode reliability index
Response mode reliability index R
typerefer to that the number of users by the scheduling of automatic response mode participation demand response accounts for the ratio of all participating user numbers, computing formula is suc as formula shown in (1):
Wherein, n is for participating in the total number of users of demand response scheduling; R
i(con) being i user's control sign, is 0-1 decision variable, when control mode is automatic response mode, and R
i(con)=1, R under other control mode
i(con)=0.
(2) response time reliability index
Response time reliability index R
timerefer to that actual schedule capacity surpasses the time scale of desirable scheduling capacity 80%;
Consider that Demand-side resource controllability is poor, the less practical factor that waits of regulation and control gear, demand response ideal curve is generally piecewise function, and every section of function is constant.Fig. 2 is certain demand response event design sketch, carries out altogether demand response scheduling events twice.
In Fig. 2, Δ
1represent the 1. desirable scheduling capacity of time demand response scheduling events, Δ
2represent the 2. desirable scheduling capacity of time demand response scheduling events, straight line represents desirable dispatch curve, and dotted line represents 80% desirable dispatch curve, and curve represents actual schedule curve.The 1. in time event, demand response scheduling time is t
1~t
5if, assert that scheduling capacity surpasses 80% of desirable scheduling capacity and dispatches successfully, the actual demand response scheduling time is t so
2~t
3and t
4~t
5.The, 2. in time event, the actual demand response scheduling time is t
6~t
7.So, in Fig. 2, the 1. response time reliability index R of time event
time, 1with the 2. response time reliability index R of time event
time, 2be respectively:
In Fig. 3, curve 1 is almost a straight line, the most approaching with ideal curve, and response time reliability index is 1; Curve 2 has certain fluctuation, but minimum scheduling capacity is still greater than 80% desirable scheduling capacity, and response time reliability index is still 1; Larger, maximum, the minimum scheduling capacity difference ratio of curve 3 fluctuation reaches 50%, and scheduling capacity is extremely unstable, and now reliability index is less than 1 the response time.For electrical network, for promoting the accuracy of response results, avoid unnecessary economic compensation or impact, response curve preferred sequence is: curve 1 > curve 2 > curves 3.
(3) response Capacity Reliability index
Response Capacity Reliability index R
caprefer at scheduling capacity and be greater than within the scope of 80% desirable scheduling capacity, the product of weight coefficient and scheduling capacity ratio,
In Fig. 4, a, b, c, d, e, f, g, h all represents time period, Δ
1represent desirable scheduling capacity, dotted line represents 80% desirable dispatch curve, and the dotted line of two comparatively denses represents 90% and 110% desirable dispatch curve.Suppose: when scheduling capacity and desirable scheduling capacity differ ratio in 10%, weight coefficient is 1; When scheduling capacity and desirable scheduling capacity differ ratio 10%-20% or scheduling capacity be greater than desirable scheduling capacity 20% time, weight coefficient is 0.8; In other situation, weight coefficient is 0.In Fig. 4, response Capacity Reliability index is so:
1.2, economic index comprises unit response Capacity Cost E
unitwith the electrical network peak-valley difference E reducing
gap, wherein:
(1) unit response Capacity Cost E
unitrefer to that under certain event, grid company is the cost that Demand-side unit response capacity is paid, computing formula is suc as formula shown in (5),
Wherein, n is for participating in the total number of users of demand response scheduling; C
i(t) be the response capacity of i user t period;
response total amount for i user in whole dispatching cycle; p
ithe economic compensation giving for electrical network under unit quantity of electricity;
response total volume for Demand-side generation.
(2) the electrical network peak-valley difference E reducing
gaprefer within the demand response event phase, because of the ratio that demand response event causes electrical network peak-valley difference to reduce, computing formula is suc as formula shown in (6),
Wherein, P
for(max), P
for(min) electrical network maximum, the minimum load of prediction while not occurring for event, P (max), P (min) is that the electrical network of event in the phase is maximum, minimum load.
1.3, feature of environmental protection index comprises the environmental pollutant discharge B of minimizing
c,Nwith user's energy consumption grade B
level, wherein:
(1) the environmental pollutant discharge B reducing
c,Nrefer to the carbon monoxide, carbon dioxide, nitrogen carbonoxide, sulphuric dioxide, particle and the hydrocarbons discharge capacity that because carrying out demand response event, reduce, computing formula is suc as formula shown in (7),
Wherein, n is for participating in the total number of users of demand response scheduling; α
j(j=1,2 ..., 6) and be pollutant discharge amount design factor, α
1for CO emission design factor, α
2cO2 emissions design factor, α
3nitrogen carbonoxide Emission amount calculation coefficient, α
4sO2 emissions design factor, α
5discharge quantity of particulate matter design factor, α
6hydrocarbons Emission amount calculation coefficient,
for response total volume.
(2) user's energy consumption grade B
levelgenerally with reference to country, provinces and cities, electrical network etc., put into effect energy consumption index appraisal standards, be divided into 1 grade of energy consumption, 2 grades of energy consumptions and 3 grades of energy consumptions, energy consumption level raises successively.
2, determine demand response effect assessment grade and pass judgment on collection
Demand response effect assessment grade is divided into 6 grades from high to low, is defined as and passes judgment on collection, the present invention represents with black, red, orange, yellow, blue, green, and the demand response effect of each grade representative is as follows:
1. black: represent demand response poor effect, there is a big difference with electrical network Expected Results;
2. red: represent demand response weak effect, the overwhelming majority response period fails to meet electrical network Expected Results;
3. orange: represent that demand response effect is poor, the partial response period fails to meet electrical network Expected Results;
4. yellow: to represent that demand response effect is general, substantially meet electrical network Expected Results;
5. blue: to represent that demand response is respond well, meet electrical network Expected Results;
6. green: represent that demand response effect is fine, user responds the best, finely meet electrical network Expected Results.
Therefore passing judgment on collection can be expressed as: V={v
1, v
2, v
3, v
4, v
5, v
6}={ is black, red, orange, and Huang, indigo plant, green.
In example of calculation shows, determined 7 indexs of the first step are set respectively to 6 quantized intervals, corresponding one of each quantized interval is passed judgment on collection.As response time reliability index, can think that this index belongs to " indigo plant " efficiency in 90%-95% interval, be " orange " efficiency in 70%-90% interval, and other by that analogy.The efficiency level judgement of all single indexs is all to advise that by data with existing interval and expert assessment and evaluation determines jointly, as concerning " user's energy consumption grade " this index, difference because of electrical network characteristic, in eastern region, more than 40% one-level efficiency can judge that electrical network is " indigo plant " efficiency, and in west area, can only judge that electrical network is " green " efficiency.That is to say, the quantized interval that evaluation rank is corresponding allows expert to set out from different perspectives, according to data with existing suggestion, certain index is belonged to which class efficiency level and make self evaluation, finally the form of expression with degree of membership reflects the preference degree of expert to the whole efficiency of certain electrical network.
3, determine demand response effect assessment set of factors
The demand response Indexes of Evaluation Effect of the first step is divided into 2 layers:
The comprehensive evaluation factor collection on upper strata is: U={u
1, u
2, u
3}={ reliability index, economic index, feature of environmental protection index },
The single index set of factors of bottom is respectively: U
1={ u
11, u
12, u
13}={ response mode reliability index, response time reliability index, response Capacity Reliability index }, U
2={ u
21, u
22}={ unit response Capacity Cost, the electrical network peak-valley difference of minimizing }, U
3={ u
31, u
32the environmental pollutant discharge of }={ minimizing, user's energy consumption grade }.
4, determine the judgment matrix R of single index set of factors
k(k=1,2,3)
The judgment matrix R of reliability index
1for:
The judgment matrix R of economic index
2for:
The judgment matrix R of feature of environmental protection index
3for:
Wherein, certain single factor of line display single index set of factors, the share that the opinion rating of this monofactorial evaluation collection is shared is shown in list; As r
115represent that expert is to u
11evaluation in, evaluation rank is " indigo plant " portion, with number percent, represents; r
123represent that expert is to u
12evaluation in, evaluation rank is " orange " portion, with number percent, represents; r
131represent that expert is to u
13evaluation in, evaluation rank is " black " portion, represents, by that analogy with number percent.
5, determine the weight coefficient matrix of factor of evaluation collection
The heavy matrix of coefficients of the comprehensive evaluation factor centralization of state power is: A=|a
1a
2a
3|, and a
1+ a
2+ a
3=1;
Single index set of factors weight coefficient matrix is respectively: A
1=| a
11a
12a
13|, and a
11+ a
12+ a
13=1, A
2=| a
21a
22|, and a
21+ a
22=1, A
3=| a
31a
32|, and a
31+ a
32=1,
Wherein, a
1, a
2, a
3for each factor in comprehensive evaluation factor collection U is passed judgment on the weight coefficient of collection, a relatively
11, a
12, a
13for single index set of factors U
1in each factor relatively pass judgment on the weight coefficient of collection, a
21, a
22for single index set of factors U
2in each factor relatively pass judgment on the weight coefficient of collection, a
31, a
32for single index set of factors U
3in each factor relatively pass judgment on the weight coefficient of collection.
6, the fuzzy evaluation result of computation requirement response scheduling effect,
First with compound operation, obtain single index set of factors U
kthe comprehensive evaluation result B of (k=1,2,3)
k:
B
k=A
k*R
k=[b
k1b
k2…b
k6];
Then solve the judgment matrix R of comprehensive evaluation factor collection U:
Then obtain the final comprehensive evaluation result B:B=A*R=[b of comprehensive evaluation factor collection U
1b
2b
6]
Wherein, b
1, b
2, b
3, b
4, b
5and b
6be that 6 evaluation ranks are " black ", " red ", " orange ", " Huang ", " indigo plant " and " green " portion, represents with number percent.The general b that selects
1, b
2, b
3, b
4, b
5and b
6middle greatest measure is as this demand response dispatching effect evaluation result.
For example, suppose A=[0.60.30.1], A
1=[0.20.40.4], A
2=[0.50.5], A
3=[0.60.4]
So,
B
1=[000.060.50.40.04],B
2=[000.20.650.150],B
3=[000.20.760.040],
B=A*R=[000.1160.5710.2890.024]
In matrix B, b
max=0.571, through expert's comprehensive evaluation, think that this demand response dispatching effect evaluation result belongs to " Huang " look, i.e. " demand response effect is general, substantially meets electrical network Expected Results ".
Concerning electrical network, the assessment result that the identical index in different regions causes may be not identical, and this need take into full account the difference characteristics such as area, user.And Field Using Fuzzy Comprehensive Assessment is converted into quantitative evaluation according to the degree of membership theory of fuzzy mathematics qualitative evaluation, by fuzzy mathematics, to being subject to things or the object of many factors restriction, make an overall evaluation, there is result clear, the feature that systematicness is strong, can solve preferably problem fuzzy, that be difficult to quantification, be applicable to demand response different stages of development.
Finally should be noted that: above embodiment is only in order to illustrate that technical scheme of the present invention is not intended to limit, although the present invention is had been described in detail with reference to above-described embodiment, those of ordinary skill in the field are to be understood that: still can modify or be equal to replacement the specific embodiment of the present invention, and do not depart from any modification of spirit and scope of the invention or be equal to replacement, it all should be encompassed in the middle of claim scope of the present invention.