CN105225033A - A kind of substation energy efficiency evaluating method based on analytical hierarchy process - Google Patents

Abstract

The invention discloses a kind of substation energy efficiency evaluating method based on analytical hierarchy process, step comprises: obtain n desired value for m transformer station to be assessed; Original index data matrix is set up and standardization generates achievement data matrix to be assessed based on gray relative analysis method; Set up target layers structure, the relative weighting of each energy efficiency evaluation index relative the superiors index of agriculture products hierarchical structure bottom, obtains relative weighting matrix; For achievement data matrix to be assessed, the correlation coefficient between the reference value calculating energy efficiency evaluation index corresponding in the desired value of each energy efficiency evaluation index after each transformer station to be assessed standardization and reference line respectively, obtains incidence coefficient matrix; Calculate the integrated assessment result matrix of m transformer station to be assessed.The present invention can not only analyze transformer station efficiency rationality degree and also can amount of energy saving size before and after assay Transfomer Substation Reconstruction, have that data assessment accuracy is high, the advantage of assessment high-speed and high-efficiency.

Description

A kind of substation energy efficiency evaluating method based on analytical hierarchy process

Technical field

The present invention relates to intelligent grid field, be specifically related to a kind of substation energy efficiency evaluating method based on analytical hierarchy process.

Background technology

In actual power distribution network, transformer station, as the hinge link between generating, transmission of electricity, distribution, accounts for the major part of whole distribution energy consumption.For ensureing substation safety, economy, green operation, requiring that the apparatus factor to transformer station, operation factor, architectural factors and environmental factor carry out comprehensive energy efficiency evaluation, namely setting up energy efficiency evaluation index system.The conventional method of current substation energy efficiency assessment has levels analytic approach.Each index divides by analytical hierarchy process, and calculates the relative weighting between different layers index.Finally draw the weight of each solution layer index for destination layer index.But analytical hierarchy process by the efficiency rationality degree analyzing transformer station, and can only really cannot draw the size of amount of energy saving before and after Transfomer Substation Reconstruction.

Summary of the invention

The technical problem to be solved in the present invention: for the problems referred to above of prior art, there is provided a kind of can not only analyze transformer station efficiency rationality degree and also can amount of energy saving size before and after assay Transfomer Substation Reconstruction, data assessment accuracy is high, the substation energy efficiency evaluating method based on analytical hierarchy process of assessment high-speed and high-efficiency.

In order to solve the problems of the technologies described above, the technical solution used in the present invention is:

Based on a substation energy efficiency evaluating method for analytical hierarchy process, step comprises:

1) for m transformer station to be assessed, the desired value of n the energy efficiency evaluation index of specifying is obtained respectively;

2) for m transformer station to be assessed, the desired value of n energy efficiency evaluation index of each transformer station to be assessed is set up the original index data matrix R shown in formula (1) based on gray relative analysis method; Obtain the capable R of optimal value generating reference of each energy efficiency evaluation index respectively ₀, by original index data matrix R and reference line R ₀in each data carry out standardization and combination obtain the achievement data matrix S to be assessed shown in formula (2);

$R=\left[\begin{array}{cccc}{r}_{11}& r_{12}& ...& r_{1n}\\ r_{21}& r_{22}& ...& r_{2n}\\ \·& \·& & \·\\ \·& \·& & \·\\ \·& \·& & \·\\ r_{m1}& r_{m2}& ...& r_{mn}\end{array}\right]---\left(1\right)$

In formula (1), R is the original index data matrix set up, r _mnit is the desired value of m transformer station n-th to be assessed energy efficiency evaluation index;

$S=\left[\begin{array}{cccc}1& 1& ...& 1\\ s_{11}& s_{12}& ...& s_{1n}\\ s_{21}& s_{22}& ...& s_{2n}\\ \·& \·& & \·\\ \·& \·& & \·\\ \·& \·& & \·\\ s_{m1}& s_{m2}& ...& s_{mn}\end{array}\right]---\left(2\right)$

In formula (2), S is the achievement data matrix to be assessed obtained, the first behavior reference line R in achievement data matrix S to be assessed ₀result after standardization, s _mnbe m the desired value of transformer station n-th energy efficiency evaluation index to be assessed after standardization;

3) for achievement data matrix S to be assessed, desired value and the reference line R of each energy efficiency evaluation index after each transformer station to be assessed standardization is calculated respectively ₀correlation coefficient between the reference value of the energy efficiency evaluation index of middle correspondence, obtains the incidence coefficient matrix β of m transformer station to be assessed;

4) desired value for described n energy efficiency evaluation index sets up target layers structure based on the classification of energy efficiency evaluation index, based on the relative weighting of each energy efficiency evaluation index relative the superiors index of analytical hierarchy process agriculture products hierarchical structure bottom, obtain the relative weighting matrix ω be made up of the relative weighting of n energy efficiency evaluation index relative the superiors index;

5) calculate the integrated assessment result matrix of m transformer station to be assessed according to formula (3), and wherein the energy efficiency evaluation result of i-th transformer station to be assessed such as formula shown in (4);

X＝β·ω(3)

In formula (3), X is the integrated assessment result matrix of m transformer station to be assessed, and β is the incidence coefficient matrix of m transformer station to be assessed, and ω is the relative weighting matrix that the relative weighting of n energy efficiency evaluation index is formed;

$x_i=\[{\mathrm{\β}}_i\left(1\right),{\mathrm{\β}}_i\left(2\right),...,{\mathrm{\β}}_i\left(n\right)\]\·\left[\begin{array}{c}{\mathrm{\ω}}_1\\ {\mathrm{\ω}}_2\\ \·\\ \·\\ \·\\ {\mathrm{\ω}}_n\end{array}\right]---\left(4\right)$

In formula (4), x _ibe the energy efficiency evaluation result of i-th transformer station to be assessed, β _in () is desired value and the reference line R of the n-th energy efficiency evaluation index of i-th transformer station to be assessed ₀correlation coefficient between the reference value of the energy efficiency evaluation index of middle correspondence, ω _nit is the relative weighting of the target indicator layer of the n-th energy efficiency evaluation index relative indicatrix hierarchical structure.

Preferably, described step 1) in the desired value of n energy efficiency evaluation index of specifying comprise electrical equipment information altogether, system running state, building structure information, environmental impact factor four classification, energy efficiency evaluation index under described electrical equipment information classification comprises nature cooling main transformer number of units accounting A1, air-conditioning Energy Efficiency Ratio A2, energy-saving lamp general power and light fixture general power accounting A3, switch board temp and humidity regulator energy-saving index A4, air-conditioning Temperature and Humidity Control mode energy efficiency indexes A5, illumination Temperature and Humidity Control mode energy efficiency indexes A6, ventilation Temperature and Humidity Control mode energy efficiency indexes A7, switch board Temperature and Humidity Control mode energy efficiency indexes A8, stand with Variant number energy saving index A9, the horizontal index A10 of ageing equipment, energy efficiency evaluation index under described system running state classification comprising load shape coefficient B1, main transformer economical operation index B2, busbar voltage deviation B3, power factor B4, alternate laod unbalance coefficient B 5, standing with becoming economical operation index B6, described building structure information comprises building orientation index C1, shape coefficient of architecture C2, heat transfer across wall index C3, the heat insulation index C4 of window, window sunshade index C5, described environmental impact factor comprises heating degree number D1 and air-conditioning and to subsist several D2.

Preferably, described step 1) in classification n energy efficiency evaluation index desired value in, the computing function expression formula of the described main transformer of cooling naturally number of units accounting A1 is such as formula shown in (5);

A1＝n/N(5)

In formula (5), A1 is nature cooling main transformer number of units accounting, and n is the main transformer number of units of nature cooling; N is the total number of units of main transformer of transformer station;

The computing function expression formula of described air-conditioning Energy Efficiency Ratio A2 is such as formula shown in (6);

$A2=\underset{i=1}{\overset{n}{\Σ}}{\mathrm{EER}}_i\·\frac{P_{out.i}}{P_{out.total}}---\left(6\right)$

In formula (6), n is the total number of units of transformer station's air-conditioning; EER _ifor i-th air-conditioning Energy Efficiency Ratio of transformer station; P _out.totalfor substation has air-conditioning total value to customize cold; P _out.ibe i-th specified refrigerating capacity of air-conditioning; Wherein i-th air-conditioning trapped energy theory of transformer station _icomputing function expression formula such as formula shown in (7);

EER＝P _out/P _in(7)

In formula (7), EER is the air-conditioning Energy Efficiency Ratio of transformer station's air-conditioning, P _outfor the refrigerating capacity of air-conditioning under rated input power, P _infor air-conditioning rated input power;

The computing function expression formula of described energy-saving lamp general power and light fixture general power accounting A3 is such as formula shown in (8);

$A3=\frac{\underset{i=1}{\overset{m}{\Σ}}{P}_{es.i}}{\underset{j=1}{\overset{n}{\Σ}}{P}_j},(m\≤n)---\left(8\right)$

In formula (8), A3 is energy-saving lamp general power and light fixture general power accounting, and n is transformer station's light fixture sum, and m is transformer station's energy saving lamp sum, P _jfor the power of transformer station's jth light fixture; P _es.ifor the power of transformer station's i-th energy saving lamp;

The computing function expression formula of described switch board temp and humidity regulator energy-saving index A4 is such as formula shown in (9);

$A4=\frac{0.5n_2+n_3}{n_1+n_2+n_3}---\left(9\right)$

In formula (9), be A4 switch board temp and humidity regulator energy-saving index, n ₁for the number of heating plate in switch board temp and humidity regulator, n ₂for the number with fan heating plate in switch board temp and humidity regulator, n ₃for the number of heat exchanger in switch board temp and humidity regulator;

The computing function expression formula of described air-conditioning Temperature and Humidity Control mode energy efficiency indexes A5, illumination Temperature and Humidity Control mode energy efficiency indexes A6, ventilation Temperature and Humidity Control mode energy efficiency indexes A7, switch board Temperature and Humidity Control mode energy efficiency indexes A8 is such as formula shown in (10);

$Ai=\frac{P_{auto}}{P},(i=5,6,7,8)---\left(10\right)$

In formula (10), Ai (i=5,6,7,8) calculated air-conditioning Temperature and Humidity Control mode energy efficiency indexes A5 or illumination Temperature and Humidity Control mode energy efficiency indexes A6 or ventilation Temperature and Humidity Control mode energy efficiency indexes A7 or switch board Temperature and Humidity Control mode energy efficiency indexes A8 is represented; P _autorepresent that control mode is the equipment general power automatically controlled, P represents the equipment general power of all control modes;

The computing function expression formula of described station Variant number energy saving index A9 is such as formula shown in (11);

$A9=\frac{\underset{i=1}{\overset{n}{\Σ}}{\mathrm{SCORE}}_i}{n}---\left(11\right)$

In formula (11), A9 is station Variant number energy saving, and n is transformer station's distribution transforming number of units; SCORE _ithe energy saving score that i-th station determined for composing submeter according to the station Variant number preset becomes;

The computing function expression formula of the horizontal index A10 of described ageing equipment is such as formula shown in (12);

$A10=\frac{\underset{i=1}{\overset{n}{\Σ}}(1+N_i/80)}{n}---\left(12\right)$

In formula (12), A10 is the horizontal index of ageing equipment, and n is the quantity of transformer station's main electric power equipment, and described main electric power equipment comprises main transformer, stands with change, shnt capacitor, isolating switch and disconnecting link, N _ibe the operation year number of i-th main electric power equipment, and when the operation year number of main electric power equipment is greater than 20 years N _ivalue get 32;

The computing function expression formula of described load shape coefficient B1 is such as formula shown in (13);

$B1=\frac{I_{eff}}{I_{ar}}---\left(13\right)$

In formula (13), B1 is load shape coefficient, I _efffor the Payload of fixed time section, I _arfor the average load of fixed time section;

The computing function expression formula of described main transformer economical operation index B2 is such as formula shown in (14);

$B2=\left\{\begin{array}{c}\mathrm{\&beta;},\mathrm{\&beta;}<{\mathrm{\&beta;}}_J\\ 1.33{\mathrm{\&beta;}}_J^2+\frac{0.75-\mathrm{\&beta;}}{0.75-{\mathrm{\&beta;}}_J}({\mathrm{\&beta;}}_J-1.33{\mathrm{\&beta;}}_J^2),{\mathrm{\&beta;}}_J\&le;\mathrm{\&beta;}<0.75\\ {\mathrm{\&beta;}}_J^2+\frac{1-\mathrm{\&beta;}}{0.25}(1.33{\mathrm{\&beta;}}_J^2-{\mathrm{\&beta;}}_J^2),0.75\&le;\mathrm{\&beta;}<1\end{array}\right.---\left(14\right)$

In formula (14), B2 is main transformer economical operation index, and β is transformer Rate of average load; β _jfor transformer economic operation point load factor;

The computing function expression formula of described busbar voltage deviation B3 is such as formula shown in (15);

$B3=\underset{j=1}{\overset{m}{\&Sigma;}}{\mathrm{\&Delta;u}}_{Tj}/m---\left(15\right)$

In formula (15), B3 is busbar voltage deviation, and m represents transformer substation voltage number of degrees, Δ u _tjrepresent the average of jth electric pressure bus integral point voltage deviation absolute value within the statistics phase, Δ u _tjcomputing function expression formula such as formula shown in (16);

${\mathrm{\&Delta;u}}_T=\underset{i=1}{\overset{n}{\&Sigma;}}\frac{|u_i-u_N|}{u_N}/n---\left(16\right)$

In formula (16), Δ u _tfor the average of certain electric pressure bus integral point voltage deviation absolute value within the statistics phase, u _ifor the integral point magnitude of voltage of bus, n is integral point number in the statistics phase, u _nfor bus rated voltage;

The computing function expression formula of described power factor B4 is such as formula shown in (17);

$B4=\underset{i=1}{\overset{n}{\&Sigma;}}{\mathrm{\&eta;}}_i\&CenterDot;\frac{P_i}{P_{total}}---\left(17\right)$

In formula (17), B4 is power factor, and n is the transforming plant main transformer number of windings, η _ifor the power factor of main transformer i side, P _ifor the active power of main transformer i side, P _totalfor the total active power of transformer station, wherein the power factor η of main transformer i side _ifor the active-power P of main transformer i side _idivided by the applied power S of main transformer i side _iobtain;

The computing function expression formula of described alternate laod unbalance coefficient B 5 is such as formula shown in (18) or formula (19);

In formula (18), B5 is alternate laod unbalance coefficient, S _mfor the maximum load power in three-phase, S _xfor the minimum load power in three-phase, for three-phase average load power;

In formula (19), B5 is alternate laod unbalance coefficient, I _mfor the maximum load current in three-phase, I _xfor the minimum load current in three-phase, for three-phase average load current;

Described station becomes the computing function expression formula of economical operation index B6 such as formula shown in (20);

$B6=\left\{\begin{array}{c}\mathrm{\&beta;},\mathrm{\&beta;}<{\mathrm{\&beta;}}_J\\ 1.33{\mathrm{\&beta;}}_J^2+\frac{0.75-\mathrm{\&beta;}}{0.75-{\mathrm{\&beta;}}_J}({\mathrm{\&beta;}}_J-1.33{\mathrm{\&beta;}}_J^2),{\mathrm{\&beta;}}_J\&le;\mathrm{\&beta;}<0.75\\ {\mathrm{\&beta;}}_J^2+\frac{1-\mathrm{\&beta;}}{0.25}(1.33{\mathrm{\&beta;}}_J^2-{\mathrm{\&beta;}}_J^2),0.75\&le;\mathrm{\&beta;}<1\end{array}\right.---\left(20\right)$

In formula (20), B6 is station change economical operation index, and β is transformer Rate of average load; β _jfor transformer economic operation point load factor;

When the building orientation of transformer station is north-south, the value of building orientation index C1 is taken as 1, when the building orientation of transformer station is East and West direction, the value of building orientation index C1 is taken as 0, when the building orientation of transformer station is between north-south and East and West direction, according to towards angle being the value that building orientation index C1 composes between 0 and 1;

The computing function expression formula of described shape coefficient of architecture C2 is such as formula shown in (21);

$C2=\frac{F_0}{V_0}---\left(21\right)$

In formula (21), C2 is shape coefficient of architecture, F ₀for the external surface area of building; V ₀for architectural volume;

The computing function expression formula of described heat transfer across wall index C3 is such as formula shown in (22);

$C3=\frac{1}{R_i+(R_1+R_2+...+R_n)+R_e}---\left(22\right)$

In formula (22), C3 is heat transfer across wall index, R _ifor the heat exchange of building enclosure inside surface hinders, R _efor the heat exchange of building enclosure outside surface hinders, R _nfor the thermal resistance of the n-th layer material of building enclosure;

The heat insulation index C4 of described window specifically refers to the heat-proof quality score of the construction window glass determined according to the heat insulation index grade form of window preset;

Described window sunshade index C5 specifically pointer to the pre-set constant value without sunshade measure, indoor sunshade measure, outdoor sunshade measure three kinds of modes of transformer station;

The computing function expression formula of described heating degree number D1 is such as formula shown in (23);

$D1=\underset{i=1}{\overset{n}{\&Sigma;}}(t_R-t_{m,i})---\left(23\right)$

In formula (23), D1 is heating degree number, and n is days of heating period or calculates number of days, t _rfor indoor reference temperature, t _m,ifor the outdoor daily mean temperature of Heating Period;

Described air-conditioning subsists the computing function expression formula of several D2 such as formula shown in (24);

$D2=\underset{i=1}{\overset{n}{\&Sigma;}}(t_{m,i}-t_R)---\left(24\right)$

In formula (24), D2 is that air-conditioning is subsisted number, and n is days of heating period or calculates number of days, t _rfor indoor reference temperature, t _m,ifor the outdoor daily mean temperature of Heating Period.

Preferably, described step 1) Plays process specifically refers to: the type judging each energy efficiency evaluation index, if the type of current energy efficiency evaluation index is the direct index that value is the bigger the better, then formula (25) is adopted to calculate the standardization result of the desired value of current energy efficiency evaluation index; If the type of the current energy efficiency evaluation index negative index that to be value the smaller the better, then adopt formula (26) to calculate the standardization result of the desired value of current energy efficiency evaluation index; If the type of the current energy efficiency evaluation index moderate osculant index of being advisable that is value, then adopt formula (27) to calculate the standardization result of the desired value of current energy efficiency evaluation index;

$z=\frac{x_{ij}-x_{min}}{x_{max}-x_{min}}---\left(25\right)$

$z=\frac{x_{max}-x_{ij}}{x_{max}-x_{min}}---\left(26\right)$

$z=\left\{\begin{array}{cc}\frac{x_{ij}-x_{\min }}{U_1-x_{\min }}& x_{\min }<x_{ij}<U_1\\ 1& U_1\&le;x_{ij}\&le;U_2\\ \frac{x_{\max }-x_{ij}}{x_{\max }-U_2}& U_2<x_{ij}<x_{\max }\end{array}\right.---\left(27\right)$

In formula (25), (26) and (27), z is the standardization result of the desired value of current energy efficiency evaluation index, x _minfor the minimum value of the desired value of current energy efficiency evaluation index, x _maxfor the maximal value of the desired value of current energy efficiency evaluation index, x _ijfor the desired value of current energy efficiency evaluation index; [U ₁, U ₂] for preset current energy efficiency evaluation index optimal zone between.

Preferably, described step 3) in specifically refer to the desired value and the reference line R that calculate each energy efficiency evaluation index after each transformer station to be assessed standardization according to formula (30) ₀correlation coefficient between the reference value of the energy efficiency evaluation index of middle correspondence;

${\mathrm{\&beta;}}_i\left(j\right)=\frac{\underset{i}{min}\underset{j}{min}|s_{0j}-s_{ij}|+\mathrm{\&rho;}\underset{i}{\max }\underset{j}{\max }|s_{0j}-s_{ij}|}{|s_{0j}-s_{ij}|+\mathrm{\&rho;}\underset{i}{\max }\underset{j}{\max }|s_{0j}-s_{ij}|}---\left(30\right)$

In formula (30), β _ij () is the desired value s of the jth energy efficiency evaluation index after the standardization of i-th transformer station to be assessed _ijwith reference line R ₀the reference value s of the energy efficiency evaluation index of middle correspondence _0jbetween correlation coefficient, ρ is the resolution ratio of each energy efficiency evaluation index, and the interval of the resolution ratio ρ of each energy efficiency evaluation index is ρ ∈ [0,1].

Preferably, described step 4) detailed step comprise:

4.1) desired value for described n energy efficiency evaluation index sets up target layers structure based on the classification of energy efficiency evaluation index, the superiors of described mark hierarchical structure are the target indicator layer that substation energy efficiency assessment result is formed, orlop is the scheme indicator layer that n energy efficiency evaluation index is formed, and is the intermediary outcomes layer that the classification of energy efficiency evaluation index is formed between described scheme indicator layer and target indicator layer; Energy efficiency evaluation index under each being classified, as one group of energy efficiency evaluation index, selects one group of energy efficiency evaluation index as current group of energy efficiency evaluation index;

4.2) determine the important coefficient of each group energy efficiency evaluation index another group energy efficiency evaluation index relative in current group of energy efficiency evaluation index, obtain the judgment matrix of the current group of energy efficiency evaluation index be made up of all important coefficients;

4.3) the eigenvalue of maximum λ of the judgment matrix of current group of energy efficiency evaluation index is calculated _max, be 1 for benchmark with the weight sum of all energy efficiency evaluation indexs of current group of energy efficiency evaluation index, by described eigenvalue of maximum λ _maxcorresponding proper vector is normalized, and obtains the relative weighting of each energy efficiency evaluation index relative the superiors index in current group of energy efficiency evaluation index;

4.4) judge whether that all groups of energy efficiency evaluation indexs are disposed, if be not yet disposed, then select next group energy efficiency evaluation index as current group of energy efficiency evaluation index, and redirect performs step 4.2); Be disposed else if, then redirect performs step 4.5);

4.5) relative weighting of each energy efficiency evaluation index relative the superiors index in all groups of energy efficiency evaluation indexs is combined, obtain the relative weighting matrix ω be made up of the relative weighting of n energy efficiency evaluation index relative the superiors index, wherein n energy efficiency evaluation index relative the superiors index relative weighting and be 1.

Preferably, described step 4.2) in when to determine in current group of energy efficiency evaluation index the importance of each group energy efficiency evaluation index another group energy efficiency evaluation index relatively, specifically refer to the importance determining another group energy efficiency evaluation index relatively of each group energy efficiency evaluation index in current group of energy efficiency evaluation index according to nine grades of scaling laws.

Preferably, described when to determine in current group of energy efficiency evaluation index the important coefficient of each group energy efficiency evaluation index another group energy efficiency evaluation index relatively according to nine grades of scaling laws, each group energy efficiency evaluation index relatively another importance organizing energy efficiency evaluation index comprises nine grades of significance levels, described nine grades of significance levels comprise identical, important a little, obviously important, strongly important, intergrade significance level in extremely important Pyatyi significance level and described five kinds of Pyatyi significance levels in the middle of arbitrary neighborhood Pyatyi significance level, and important coefficient corresponding to described nine grades of significance levels is followed successively by 1 ~ 9 according to importance degree.

Preferably, described step 4.3) also comprise the step of the judgment matrix of current group of energy efficiency evaluation index being carried out to consistency desired result, and if verification pass through, then redirect performs step 4.4), otherwise redirect performs step 4.2 again).

Preferably, the detailed step that the described judgment matrix to current group of energy efficiency evaluation index carries out consistency desired result comprises:

4.3.1) judge the exponent number of the judgment matrix of current group of energy efficiency evaluation index, if the exponent number of the judgment matrix of current group of energy efficiency evaluation index is no more than 2, then directly redirect performs step 4.4), otherwise redirect performs next step;

4.3.2) the inconsistency degree quantitative index of the judgment matrix of current group of energy efficiency evaluation index is calculated according to formula (28);

$CI=\frac{{\mathrm{\&lambda;}}_{max}-n}{n-1}---\left(28\right)$

In formula (28), CI is the inconsistency degree quantitative index of the judgment matrix of current group of energy efficiency evaluation index, λ _maxfor the eigenvalue of maximum of the judgment matrix of current group of energy efficiency evaluation index, n is the quantity of all energy efficiency evaluation indexs;

4.3.3) search default Aver-age Random Consistency Index and energy efficiency evaluation index schedule of quantities according to the quantity n of energy efficiency evaluation index, obtain the Aver-age Random Consistency Index value RI of the judgment matrix of current group of energy efficiency evaluation index;

4.3.4) the inconsistency degree quantitative index modified value CR of the judgment matrix of current group of energy efficiency evaluation index is calculated according to formula (29);

$CR=\frac{CI}{RI}---\left(29\right)$

In formula (29), CR is the inconsistency degree quantitative index modified value of the judgment matrix of current group of energy efficiency evaluation index, CI is the inconsistency degree quantitative index of the judgment matrix of current group of energy efficiency evaluation index, and RI is the Aver-age Random Consistency Index value of the judgment matrix of current group of energy efficiency evaluation index;

4.3.5) judge that the inconsistency degree quantitative index modified value CR of the judgment matrix of current group of energy efficiency evaluation index is less than 0.1 and whether sets up, if set up, judge that carrying out consistency desired result to the judgment matrix of current group of energy efficiency evaluation index passes through, otherwise judge that carrying out consistency desired result to the judgment matrix of current group of energy efficiency evaluation index does not pass through.

The substation energy efficiency evaluating method that the present invention is based on analytical hierarchy process has following advantage:

The desired value of n energy efficiency evaluation index after 1, the present invention utilizes m transformer station to be assessed, the standardization of each transformer station to be assessed sets up achievement data matrix S to be assessed based on gray relative analysis method, calculates desired value and the reference line R of each energy efficiency evaluation index after each transformer station to be assessed standardization respectively ₀correlation coefficient between the reference value of the energy efficiency evaluation index of middle correspondence, obtains the incidence coefficient matrix β of m transformer station to be assessed, simultaneously, desired value for n energy efficiency evaluation index sets up target layers structure based on the classification of energy efficiency evaluation index, based on the relative weighting of each energy efficiency evaluation index relative the superiors index of analytical hierarchy process agriculture products hierarchical structure bottom, obtain the relative weighting matrix ω be made up of the relative weighting of n energy efficiency evaluation index relative the superiors index, the integrated assessment result matrix X of m transformer station to be assessed is calculated based on X=β ω, the correlativity of m transformer station to be assessed can be made full use of, the efficiency rationality degree of transformer station can not only be analyzed, and can amount of energy saving size before and after assay Transfomer Substation Reconstruction, there is data assessment accuracy high, the advantage of assessment high-speed and high-efficiency.

2, the disposable substation energy efficiency assessment that can complete m transformer station to be assessed of the present invention, has the advantage of assessment high-speed and high-efficiency, is particularly suitable for the energy efficiency evaluation of transformer station in enormous quantities.

Accompanying drawing explanation

Fig. 1 is the basic procedure schematic diagram of embodiment of the present invention method.

Fig. 2 is the schematic diagram of the target layers structure set up in the embodiment of the present invention.

Embodiment

As shown in Figure 1, the present embodiment comprises based on the step of the substation energy efficiency evaluating method of analytical hierarchy process:

1) for m transformer station to be assessed, the desired value of n the energy efficiency evaluation index of specifying is obtained respectively;

2) for m transformer station to be assessed, the desired value of n energy efficiency evaluation index of each transformer station to be assessed is set up the original index data matrix R shown in formula (1) based on gray relative analysis method; Obtain the capable R of optimal value generating reference of each energy efficiency evaluation index respectively ₀, by original index data matrix R and reference line R ₀in each data carry out standardization and combination obtain the achievement data matrix S to be assessed shown in formula (2);

$R=\left[\begin{array}{cccc}{r}_{11}& r_{12}& ...& r_{1n}\\ r_{21}& r_{22}& ...& r_{2n}\\ \&CenterDot;& \&CenterDot;& & \&CenterDot;\\ \&CenterDot;& \&CenterDot;& & \&CenterDot;\\ \&CenterDot;& \&CenterDot;& & \&CenterDot;\\ r_{m1}& r_{m2}& ...& r_{mn}\end{array}\right]---\left(1\right)$

In formula (1), R is the original index data matrix set up, r _mnit is the desired value of m transformer station n-th to be assessed energy efficiency evaluation index;

$S=\left[\begin{array}{cccc}1& 1& ...& 1\\ s_{11}& s_{12}& ...& s_{1n}\\ s_{21}& s_{22}& ...& s_{2n}\\ \&CenterDot;& \&CenterDot;& & \&CenterDot;\\ \&CenterDot;& \&CenterDot;& & \&CenterDot;\\ \&CenterDot;& \&CenterDot;& & \&CenterDot;\\ s_{m1}& s_{m2}& ...& s_{mn}\end{array}\right]---\left(2\right)$

In formula (2), S is the achievement data matrix to be assessed obtained, the first behavior reference line R in achievement data matrix S to be assessed ₀result after standardization, s _mnbe the desired value of the n-th energy efficiency evaluation index after standardization of m transformer station to be assessed;

3) for achievement data matrix S to be assessed, desired value and the reference line R of each energy efficiency evaluation index after each transformer station to be assessed standardization is calculated respectively ₀correlation coefficient between the reference value of the energy efficiency evaluation index of middle correspondence, obtains the incidence coefficient matrix β of m transformer station to be assessed;

4) desired value for n energy efficiency evaluation index sets up target layers structure based on the classification of energy efficiency evaluation index, based on the relative weighting of each energy efficiency evaluation index relative the superiors index of analytical hierarchy process agriculture products hierarchical structure bottom, obtain the relative weighting matrix ω be made up of the relative weighting of n energy efficiency evaluation index relative the superiors index;

5) calculate the integrated assessment result matrix of m transformer station to be assessed according to formula (3), and wherein the energy efficiency evaluation result of i-th transformer station to be assessed such as formula shown in (4);

X＝β·ω(3)

In formula (3), X is the integrated assessment result matrix of m transformer station to be assessed, and β is the incidence coefficient matrix of m transformer station to be assessed, and ω is the relative weighting matrix that the relative weighting of n energy efficiency evaluation index is formed;

$x_i=\&lsqb;{\mathrm{\&beta;}}_i\left(1\right),{\mathrm{\&beta;}}_i\left(2\right),...,{\mathrm{\&beta;}}_i\left(n\right)\&rsqb;\&CenterDot;\left[\begin{array}{c}{\mathrm{\&omega;}}_1\\ {\mathrm{\&omega;}}_2\\ \&CenterDot;\\ \&CenterDot;\\ \&CenterDot;\\ {\mathrm{\&omega;}}_n\end{array}\right]---\left(4\right)$

In formula (4), x _ibe the energy efficiency evaluation result of i-th transformer station to be assessed, β _in () is desired value and the reference line R of the n-th energy efficiency evaluation index of i-th transformer station to be assessed ₀correlation coefficient between the reference value of the energy efficiency evaluation index of middle correspondence, ω _nit is the relative weighting of the target indicator layer of the n-th energy efficiency evaluation index relative indicatrix hierarchical structure.The energy efficiency evaluation result x of i-th transformer station to be assessed _ilarger, then the efficiency of i-th transformer station to be assessed is close to the efficiency of desirable transformer station, illustrates that the efficiency of i-th transformer station to be assessed is higher, thus can discharge each quality participating in assessment substation energy efficiency, finally obtains the result assessed.

See Fig. 2, step 1) in the desired value of n energy efficiency evaluation index of specifying comprise electrical equipment information altogether, system running state, building structure information, environmental impact factor four classification, energy efficiency evaluation index under electrical equipment information classification comprises nature cooling main transformer number of units accounting A1, air-conditioning Energy Efficiency Ratio A2, energy-saving lamp general power and light fixture general power accounting A3, switch board temp and humidity regulator energy-saving index A4, air-conditioning Temperature and Humidity Control mode energy efficiency indexes A5, illumination Temperature and Humidity Control mode energy efficiency indexes A6, ventilation Temperature and Humidity Control mode energy efficiency indexes A7, switch board Temperature and Humidity Control mode energy efficiency indexes A8, stand with Variant number energy saving index A9, the horizontal index A10 of ageing equipment, energy efficiency evaluation index under system running state classification comprising load shape coefficient B1, main transformer economical operation index B2, busbar voltage deviation B3, power factor B4, alternate laod unbalance coefficient B 5, standing with becoming economical operation index B6, building structure information comprises building orientation index C1, shape coefficient of architecture C2, heat transfer across wall index C3, the heat insulation index C4 of window, window sunshade index C5, environmental impact factor comprises heating degree number D1 and air-conditioning and to subsist several D2.It should be noted that, energy efficiency evaluation choose targets is more, more comprehensive, then will be more accurate to the energy efficiency evaluation assessment of transformer station, certainly, above-mentioned energy efficiency evaluation index in the present embodiment is only exemplary illustration, the embodiment of the present embodiment technical scheme not merely depends on the concrete selection form of above-mentioned energy efficiency evaluation index, those skilled in the art also can according to above-mentioned enlightenment, carry out increasing new energy efficiency evaluation index to the selection of energy efficiency evaluation index, deletion energy efficiency evaluation index or part energy efficiency evaluation index is carried out changing etc., its principle is identical with the present embodiment, therefore no longer carry out expansion explanation at this.

In the present embodiment, step 1) in classification n energy efficiency evaluation index desired value in:

◇ cools the computing function expression formula of main transformer number of units accounting A1 naturally such as formula shown in (5);

A1＝n/N(5)

In formula (5), A1 is nature cooling main transformer number of units accounting, and n is the main transformer number of units of nature cooling; N is the total number of units of main transformer of transformer station; The main transformer type of cooling mainly contains nature cooling, forcing functions mode, naturally cools main transformer number of units accounting A1 can embody transformer station efficiency situation from the angle of the main transformer type of cooling.

The computing function expression formula of ◇ air-conditioning Energy Efficiency Ratio A2 is such as formula shown in (6);

$A2=\underset{i=1}{\overset{n}{\&Sigma;}}{\mathrm{EER}}_i\&CenterDot;\frac{P_{out.i}}{P_{out.total}}---\left(6\right)$

In formula (6), n is the total number of units of transformer station's air-conditioning; EER _ifor i-th air-conditioning Energy Efficiency Ratio of transformer station; P _out.totalfor substation has air-conditioning total value to customize cold; P _out.ibe i-th specified refrigerating capacity of air-conditioning; Wherein i-th air-conditioning trapped energy theory of transformer station _icomputing function expression formula such as formula shown in (7);

EER＝P _out/P _in(7)

In formula (7), EER is the air-conditioning Energy Efficiency Ratio of transformer station's air-conditioning, P _outfor the refrigerating capacity of air-conditioning under rated input power, P _infor air-conditioning rated input power; Air-conditioning Energy Efficiency Ratio is larger, represent that the efficiency of air-conditioning is higher, therefore after the trapped energy theory calculating the every platform air-conditioning of transformer station, adopt calculated with weighted average method to go out the Energy Efficiency Ratio of whole transformer station air-conditioning again, the ratio that weighting coefficient gets every platform air conditioner refrigerating amount and all air-conditioning overall refrigerating effects obtains air-conditioning Energy Efficiency Ratio A2.

The computing function expression formula of ◇ energy-saving lamp general power and light fixture general power accounting A3 is such as formula shown in (8);

$A3=\frac{\underset{i=1}{\overset{m}{\&Sigma;}}{P}_{es.i}}{\underset{j=1}{\overset{n}{\&Sigma;}}{P}_j},(m\&le;n)---\left(8\right)$

In formula (8), A3 is energy-saving lamp general power and light fixture general power accounting, and n is transformer station's light fixture sum, and m is transformer station's energy saving lamp sum, P _jfor the power of transformer station's jth light fixture; P _es.ifor the power of transformer station's i-th energy saving lamp; For substation illumination energy saving, require that the life-span of light fixture is long, luminescence efficiency is high, the utilance of light fixture is high, to reduce illuminator light fixture quantity, improves light fixture service efficiency, to reach energy-saving and cost-reducing object.Therefore energy saving lamp should be selected as far as possible, as Metal halogen lamp, high-pressure mercury lamp, electrodeless florescent lamp, LED etc., in this present embodiment, use the ratio of transformer station's energy-saving lamp general power and light fixture general power to weigh the energy saving of lighting.

The computing function expression formula of ◇ switch board temp and humidity regulator energy-saving index A4 is such as formula shown in (9);

$A4=\frac{0.5n_2+n_3}{n_1+n_2+n_3}---\left(9\right)$

In formula (9), be A4 switch board temp and humidity regulator energy-saving index, n ₁for the number of heating plate in switch board temp and humidity regulator, n ₂for the number with fan heating plate in switch board temp and humidity regulator, n ₃for the number of heat exchanger in switch board temp and humidity regulator; Generally the device regulated humiture is all housed in the switch board of transformer station, for equipment component in bin provides stable suitable running temperature, thus ensures its serviceable life and electric property, avoid the impact being subject to climatic factor.Conventional temp and humidity regulator has heating plate, band fan heating plate, heat exchanger.In general, the efficiency of heating plate is minimum, and band fan heating plate Energy Efficiency Ratio heating plate is slightly high, and heat exchanger efficiency is the highest among three.Adopt the mode of marking accordingly, the energy saving of heating plate is 0 point, and the energy saving of band fan heating plate is 0.5 point, and the energy saving of heat exchanger is 1 point, and the score of all switch board temp and humidity regulators of transformer station is averaged, obtain switch board temp and humidity regulator energy-saving index A4.

The control mode of transformer station's air-conditioning, illumination, ventilation and switch board moisture control system on substation energy efficiency impact closely.Compare with manual control mode, air-conditioning and ventilating system adopt humiture automatically to control, illuminator adopts acoustic control, the automatic control mode such as light-operated, switch board Temperature and Humidity Control is controlled automatically by testing environment temperature and condensation situation of change, adopt above automatic control mode all can reduce operation hours, improve the efficiency of equipment.When control mode is Non-follow control, desired value is taken as 0, represents that efficiency is minimum; When control mode for desired value when automatically controlling gets 1, represent that efficiency is best.But, in actual transformer station, equipment component Non-follow control may be there is, situation that equipment component controls automatically.Such as, exterior lighting adopts light-operated automatic control, and room lighting adopts Non-follow control; Because switch board manufacturer is different, also there is some switch board Non-follow control temp and humidity regulator, some switch board adopts automatic control mode, therefore expresses for this kind of desired value automatic control equipment power accounting in the present embodiment.

The computing function expression formula of ◇ air-conditioning Temperature and Humidity Control mode energy efficiency indexes A5, illumination Temperature and Humidity Control mode energy efficiency indexes A6, ventilation Temperature and Humidity Control mode energy efficiency indexes A7, switch board Temperature and Humidity Control mode energy efficiency indexes A8 is such as formula shown in (10);

$Ai=\frac{P_{auto}}{P},(i=5,6,7,8)---\left(10\right)$

In formula (10), Ai (i=5,6,7,8) calculated air-conditioning Temperature and Humidity Control mode energy efficiency indexes A5 or illumination Temperature and Humidity Control mode energy efficiency indexes A6 or ventilation Temperature and Humidity Control mode energy efficiency indexes A7 or switch board Temperature and Humidity Control mode energy efficiency indexes A8 is represented; P _autorepresent that control mode is the equipment general power automatically controlled, P represents the equipment general power of all control modes.

◇ station by the computing function expression formula of Variant number energy saving index A9 such as formula shown in (11);

$A9=\frac{\underset{i=1}{\overset{n}{\&Sigma;}}{\mathrm{SCORE}}_i}{n}---\left(11\right)$

In formula (11), A9 is station Variant number energy saving, and n is transformer station's distribution transforming number of units; SCORE _ithe energy saving score that i-th station determined for composing submeter according to the station Variant number preset becomes; Whether stand by change energy saving mainly from being energy-saving substation transformer angle consideration, adopt scoring to main distribution transforming range of models energy saving marking assignment, in the present embodiment, it is specifically as shown in table 1 that the station Variant number preset composes submeter.

Table 1: the station Variant number preset composes submeter.

Distribution transforming model	64 standards	73 standards	86 standards	S7	S9	S11	S13	S15
									Score	1	2	3	4	5	6	7	8

The computing function expression formula of the horizontal index A10 of ◇ ageing equipment is such as formula shown in (12);

$A10=\frac{\underset{i=1}{\overset{n}{\&Sigma;}}(1+N_i/80)}{n}---\left(12\right)$

In formula (12), A10 is the horizontal index of ageing equipment, and n is the quantity of transformer station's main electric power equipment, and main electric power equipment comprises main transformer, stands with change, shnt capacitor, isolating switch and disconnecting link, N _ibe the operation year number of i-th main electric power equipment, and when the operation year number of main electric power equipment is greater than 20 years N _ivalue get 32; Converting station electric power ageing equipment not only increases the probability broken down when equipment runs, and equipment energy consumption also can be made to increase, and reduces equipment work efficiency.Use the mean value of transformer station's main electric power ageing equipment loss factor as the horizontal index A10 of ageing equipment in the present embodiment.

The computing function expression formula of ◇ load shape coefficient B1 is such as formula shown in (13);

$B1=\frac{I_{eff}}{I_{ar}}---\left(13\right)$

In formula (13), B1 is load shape coefficient, I _efffor the Payload of fixed time section, I _arfor the average load of fixed time section; Load shape coefficient B1 is defined as Payload I in a period of time _effwith average load I _arratio.The size of B1 value is relevant with the lasting load curve that straight line changes, B1 value is larger, the amplitude of variation of load curve is larger, peak and the low ebb difference of curve are also larger, load current is larger by the loss produced during power equipment, when B1 value close to 1 time, load curve is approximately a flat curve, minimum on the loss of power equipment impact.Wherein, Payload I _effcomputing function expression formula for such as formula shown in (13-1); Average load I _arcomputing function expression formula for such as formula shown in (13-2);

$I_{eff}=\sqrt{\underset{0}{\overset{T}{\&Integral;}}{I}^2\left(t\right)dt/T}---(13-1)$

$I_{ar}=\underset{0}{\overset{T}{\&Integral;}}I\left(t\right)dt/T---(13-2)$

In formula (13-1) and formula (13-2), I _efffor the Payload of fixed time section, I _arfor the average load of fixed time section, T represents fixed time section, and I (t) represents the load current of t, therefore Payload I in fixed time section T _effwith average load I _arin the general available statistical phase, integral point load current carries out simplification calculating.

The computing function expression formula of ◇ main transformer economical operation index B2 is such as formula shown in (14);

$B2=\left\{\begin{array}{c}\mathrm{\&beta;},\mathrm{\&beta;}<{\mathrm{\&beta;}}_J\\ 1.33{\mathrm{\&beta;}}_J^2+\frac{0.75-\mathrm{\&beta;}}{0.75-{\mathrm{\&beta;}}_J}({\mathrm{\&beta;}}_J-1.33{\mathrm{\&beta;}}_J^2),{\mathrm{\&beta;}}_J\&le;\mathrm{\&beta;}<0.75\\ {\mathrm{\&beta;}}_J^2+\frac{1-\mathrm{\&beta;}}{0.25}(1.33{\mathrm{\&beta;}}_J^2-{\mathrm{\&beta;}}_J^2),0.75\&le;\mathrm{\&beta;}<1\end{array}\right.---\left(14\right)$

In formula (14), B2 is main transformer economical operation index, and β is transformer Rate of average load; β _jfor transformer economic operation point load factor; From formula (14), the value of main transformer economical operation index B2 is larger, represents that main transformer runs more economical.For two-winding transformer in the present embodiment, provide the calculating formula of main transformer economical operation index, three winding main transformer economical operation index calculating method is similar with it.Division according to transformer economic operation interval: when transformer Rate of average load β is less than economical operating point load factor β _jtime, get B2=β; When transformer Rate of average load β is positioned between the right half-court of optimal economic operational area, i.e. β _jduring≤β <0.75, between left half-court Rate of average load being mapped to optimal economic operational area, as economical operation desired value, namely such as formula shown in (14-1); When 0.75≤β≤1, when being namely positioned at the right interval of Economic moving region, Rate of average load is mapped to the left interval of Economic moving region, as economical operation desired value, namely such as formula shown in (14-2);

$B2=1.33{\mathrm{\&beta;}}_J^2+\frac{0.75-\mathrm{\&beta;}}{0.75-{\mathrm{\&beta;}}_J}({\mathrm{\&beta;}}_J-1.33{\mathrm{\&beta;}}_J^2)---(14-1)$

$B2={\mathrm{\&beta;}}_J^2+\frac{1-\mathrm{\&beta;}}{0.25}(1.33{\mathrm{\&beta;}}_J^2-{\mathrm{\&beta;}}_J^2)---(14-2)$

Formula (14-1) and each variable in (14-2) are identical with in formula (14), do not repeat them here.

The computing function expression formula of ◇ busbar voltage deviation B3 is such as formula shown in (15);

$B3=\underset{j=1}{\overset{m}{\&Sigma;}}{\mathrm{\&Delta;u}}_{Tj}/m---\left(15\right)$

In formula (15), B3 is busbar voltage deviation, and m represents transformer substation voltage number of degrees, Δ u _tjrepresent the average of jth electric pressure bus integral point voltage deviation absolute value within the statistics phase, Δ u _tjcomputing function expression formula such as formula shown in (16);

${\mathrm{\&Delta;u}}_T=\underset{i=1}{\overset{n}{\&Sigma;}}\frac{|u_i-u_N|}{u_N}/n---\left(16\right)$

In formula (16), Δ u _tfor the average of certain electric pressure bus integral point voltage deviation absolute value within the statistics phase, u _ifor the integral point magnitude of voltage of bus, n is integral point number in the statistics phase, u _nfor bus rated voltage;

In transformer station, the virtual voltage of a certain bus and the number percent of difference to bus rated voltage of this bus rated voltage are called the voltage deviation of this bus, are shown below.Busbar voltage deviation is one of important indicator investigating electric energy quality for substation.Cause the factor of substation bus bar voltage deviation to have idle underpower, reactive-load compensation are excessive, line transmission apart from long, electric load is overweight and kick the beam, wherein reactive power deficiency is the main cause causing voltage deviation.Therefore, adopt the voltage deviation Δ u of each electric pressure bus of transformer station to weigh transformer substation voltage quality and reactive-load compensation situation, the computing function expression formula of voltage deviation Δ u is such as formula shown in (16-1); When substation bus bar voltage data is each integral point voltage u in statistics phase T _itime, its voltage deviation adopts formula (16) to calculate.

$\mathrm{\&Delta;}u=\frac{u-u_N}{u_N}---(16-1)$

In formula (16-1), Δ u is busbar voltage deviation, and u is bus virtual voltage, u _nfor bus rated voltage.When considering each electric pressure bus of transformer station, busbar voltage Deviation Indices is each electric pressure busbar voltage deviation average, namely such as formula shown in (15).

The computing function expression formula of ◇ power factor B4 is such as formula shown in (17);

$B4=\underset{i=1}{\overset{n}{\&Sigma;}}{\mathrm{\&eta;}}_i\&CenterDot;\frac{P_i}{P_{total}}---\left(17\right)$

In formula (17), B4 is power factor, and n is the transforming plant main transformer number of windings, η _ifor the power factor of main transformer i side, P _ifor the active power of main transformer i side, P _totalfor the total active power of transformer station, wherein the power factor η of main transformer i side _ifor the active-power P of main transformer i side _idivided by the applied power S of main transformer i side _iobtain;

Improve transformer station power factor for reduction electric energy loss, increasing economic efficiency has a very important role, the power factor of power supply department to transforming plant main transformer has certain standard-required.The power factor of main transformer i side is such as formula shown in (17-1);

In formula (17-1), η _iwith be the power factor of main transformer i side, P _ifor the active power of main transformer i side; S _ifor the applied power of main transformer i side.After the power factor calculating each electric pressure side of main transformer, then adopt calculated with weighted average method to go out the power factor of whole transformer station, weighting coefficient gets the ratio that the burden with power of each side accounts for total burden with power, shown in (17).

The computing function expression formula of the alternate laod unbalance coefficient B 5 of ◇ is such as formula shown in (18) or formula (19);

In formula (18), B5 is alternate laod unbalance coefficient, S _mfor the maximum load power in three-phase, S _xfor the minimum load power in three-phase, for three-phase average load power;

In formula (19), B5 is alternate laod unbalance coefficient, I _mfor the maximum load current in three-phase, I _xfor the minimum load current in three-phase, for three-phase average load current;

To the loss of the power equipment such as main transformer and circuit be caused to increase by the threephase load imbalance of transformer station, work efficiency declines.With the calculated value of alternate laod unbalance coefficient B 5 as this index in the present embodiment, the alternate laod unbalance of the larger explanation of desired value of alternate laod unbalance coefficient B 5 is more serious.

◇ station is with becoming the computing function expression formula of economical operation index B6 such as formula shown in (20);

$B6=\left\{\begin{array}{c}\mathrm{\&beta;},\mathrm{\&beta;}<{\mathrm{\&beta;}}_J\\ 1.33{\mathrm{\&beta;}}_J^2+\frac{0.75-\mathrm{\&beta;}}{0.75-{\mathrm{\&beta;}}_J}({\mathrm{\&beta;}}_J-1.33{\mathrm{\&beta;}}_J^2),{\mathrm{\&beta;}}_J\&le;\mathrm{\&beta;}<0.75\\ {\mathrm{\&beta;}}_J^2+\frac{1-\mathrm{\&beta;}}{0.25}(1.33{\mathrm{\&beta;}}_J^2-{\mathrm{\&beta;}}_J^2),0.75\&le;\mathrm{\&beta;}<1\end{array}\right.---\left(20\right)$

In formula (20), B6 is station change economical operation index, and β is transformer Rate of average load; β _jfor transformer economic operation point load factor; Stand and to provide safeguard with becoming productive life electricity consumption in station, standing belongs to substation transformer with change, and general capacity is less, be in light running and hot stand-by duty for a long time, and therefore its loss is comparatively large, and efficiency is lower.Aforesaid main transformer economical operation index calculating method is also applicable to station completely and becomes economical operation.

◇ is when the building orientation of transformer station is north-south, the value of building orientation index C1 is taken as 1, when the building orientation of transformer station is East and West direction, the value of building orientation index C1 is taken as 0, when the building orientation of transformer station is between north-south and East and West direction, according to towards angle being the value that building orientation index C1 composes between 0 and 1;

The computing function expression formula of ◇ shape coefficient of architecture C2 is such as formula shown in (21);

$C2=\frac{F_0}{V_0}---\left(21\right)$

In formula (21), C2 is shape coefficient of architecture, F ₀for the external surface area of building; V ₀for architectural volume; Shape coefficient of building is the external surface area of building and the ratio of its volume.Shape coefficient of architecture index reflects the size of transformer station's architectural exterior-protecting construction area of dissipation, and Shape Coefficient desired value is larger, and illustrate that architectural exterior-protecting construction area of dissipation is larger, heat consump-tion in transfer is corresponding also larger.When adding up outer surface of building and amassing, ground and staircase partition wall etc. should be got rid of without the area of cooling and heating load.

The computing function expression formula of ◇ heat transfer across wall index C3 is such as formula shown in (22);

$C3=\frac{1}{R_i+(R_1+R_2+...+R_n)+R_e}---\left(22\right)$

In formula (22), C3 is heat transfer across wall index, R _ifor the heat exchange of building enclosure inside surface hinders, R _efor the heat exchange of building enclosure outside surface hinders, R _nfor the thermal resistance of the n-th layer material of building enclosure.

The energy resource consumption of thermal and insulating performance to building of transformer station's architectural exterior-protecting construction is in close relations.Good building enclosure can the outdoor temperature difference in holding chamber, can the outer heat of high temperature of air lock in summer, and winter can reduce indoor thermal loss, maintenance indoor temperature.The major parameter weighing heat transfer across wall performance is heat transfer across wall resistance and heat transfer coefficient.From the building enclosure side resistance be delivered to suffered by opposite side, heat is called that heat transfer across wall hinders, homogenous material heat transfer across wall resistance calculating formula is such as formula shown in (22-1);

R＝δ/λ _c(22-1)

In formula (22-1), R is the resistance of homogenous material heat transfer across wall, and unit is (m ²k)/W; δ is layer thickness, and unit is m; λ _cfor the coefficient of heat conductivity of material, unit is W/ (mK).When the material of building enclosure has multilayer, its resistance of heat transfer is each layer resistance of heat transfer sum, and the resistance of heat transfer that therefore building enclosure is total is such as formula shown in (22-2);

R ₀＝R _i+(R ₁+R ₂++R _n)+R _e(22-2)

In formula (22-2), R ₀for total resistance of heat transfer; R _i, R _ebe respectively the heat exchange resistance of inside surface and outside surface; R ₁, R ₂, R _nfor the thermal resistance of layers of material.

Inverse is got to the total resistance of heat transfer of building enclosure, the Coefficient K of building enclosure can be obtained such as formula shown in (22-3);

$K=\frac{1}{R_0},W/(m^2\&CenterDot;K)---(22-3)$

In formula (22-3), K is that (unit is W/ (m for the heat transfer coefficient of building enclosure ²k)), R ₀for total resistance of heat transfer; On this basis, get enclosure structure heat transfer coefficient (heat transfer across wall index C3) in the present embodiment and go as index, and the computing function expression formula of heat transfer across wall index C3 is such as formula shown in (22); From the index calculate formula shown in formula (22), desired value is less, and represent that heat transfer coefficient is less, the thermal and insulating performance of transformer station's architectural exterior-protecting construction is better, and building efficiency is higher.

The heat insulation index C4 of ◇ window specifically refers to the heat-proof quality score of the construction window glass determined according to the heat insulation index grade form of window preset; The efficiency impact of heat-proof quality on air-conditioning and heating in transformer station's buildings of construction window glass is very large.Carry out efficiency scoring according to the window-glass energy saving of current main Types herein, refer to target value as this, marking the results are shown in following table 2.Score value higher expression window heat-proof quality is better.

Table 2: the heat insulation index grade form of window.

Type of glass	Low-e glass	Coated glass	Antisolar glass	Double glazing	Simple glass
						Score	5	4	3	2	1

◇ window sunshade index C5 specifically pointer to the pre-set constant value without sunshade measure, indoor sunshade measure, outdoor sunshade measure three kinds of modes of transformer station; Window sunshade is to preventing transformer station's Summer Indoor the temperature rises effect obviously.Can be divided into indoor sunshade and outdoor sunshade according to Types ofshading, indoor sunshade generally adopts the mode such as curtain, blind; Outdoor sunshade generally adopts the mode such as sunshade, roller shutter.Indoor sunshade is owing to being enter indoor laggard row relax to sunlight, and shaded effects generally do not have outdoor shaded effects good.When transformer station's building window is without sunshade measure, this desired value gets 0; Have family internal sunshade measure time, get 0.5; 1 is got during outdoor sunshade, specifically also can see shown in formula (22-4).

In formula (22-4), C5 is window sunshade index.

The computing function expression formula of ◇ heating degree number D1 is such as formula shown in (23);

$D1=\underset{i=1}{\overset{n}{\&Sigma;}}(t_R-t_{m,i})---\left(23\right)$

In formula (23), D1 is heating degree number, and unit is DEG C id or dd, n is days of heating period or calculating number of days, and unit is d, t _rfor indoor reference temperature, unit is DEG C (China's most area generally gets 19 DEG C), t _m,ifor the outdoor daily mean temperature of Heating Period, unit is DEG C; Heating degree number D1 in the present embodiment gets heating degree number HDD value, is and deducts outdoor daily mean temperature with the number of degrees of indoor basic air temperature, and cumulative result of suing for peace day by day.

◇ air-conditioning subsists the computing function expression formula of several D2 such as formula shown in (24);

$D2=\underset{i=1}{\overset{n}{\&Sigma;}}(t_{m,i}-t_R)---\left(24\right)$

In formula (24), D2 is that air-conditioning is subsisted number, and n is days of heating period or calculates number of days, t _rfor indoor reference temperature (China t _rgenerally get 25 DEG C), t _m,ifor the outdoor daily mean temperature of Heating Period, its unit definition is identical with formula (23).Air-conditioning in the present embodiment several D2 that subsists gets air-conditioning and to subsist several CDD value, deducts indoor basic air temperature with outdoor daily mean temperature, and cumulative result of suing for peace day by day.

In the present embodiment, step 2) Plays process specifically refers to: the type judging each energy efficiency evaluation index, if the type of current energy efficiency evaluation index is the direct index that value is the bigger the better, then formula (25) is adopted to calculate the standardization result of the desired value of current energy efficiency evaluation index; If the type of the current energy efficiency evaluation index negative index that to be value the smaller the better, then adopt formula (26) to calculate the standardization result of the desired value of current energy efficiency evaluation index; If the type of the current energy efficiency evaluation index moderate osculant index of being advisable that is value, then adopt formula (27) to calculate the standardization result of the desired value of current energy efficiency evaluation index;

$z=\frac{x_{ij}-x_{min}}{x_{max}-x_{min}}---\left(25\right)$

$z=\frac{x_{max}-x_{ij}}{x_{max}-x_{min}}---\left(26\right)$

$z=\left\{\begin{array}{cc}\frac{x_{ij}-x_{\min }}{U_1-x_{\min }}& x_{\min }<x_{ij}<U_1\\ 1& U_1\&le;x_{ij}\&le;U_2\\ \frac{x_{\max }-x_{ij}}{x_{\max }-U_2}& U_2<x_{ij}<x_{\max }\end{array}\right.---\left(27\right)$

In formula (25), (26) and (27), z is the standardization result of the desired value of current energy efficiency evaluation index, x _minfor the minimum value of the desired value of current energy efficiency evaluation index, x _maxfor the maximal value of the desired value of current energy efficiency evaluation index, x _ijfor the desired value of current energy efficiency evaluation index; [U ₁, U ₂] for preset current energy efficiency evaluation index optimal zone between.If [x _min, x _max] be the constant interval of certain evaluation index value, x _minfor the minimum value that this index is desirable, x _maxfor the maximal value that this index is desirable, then formula (25) ~ formula (27) is adopted evaluation index to be standardized as dimensionless number between [0,1].

The basis of gray relative analysis method builds gray relative factor space, supposes that X is gray relative factor space, its sequence comprised is x _i, x _i=(x _i(1), x _i(2), x _i(n)), i ∈ I={1,2, m}, meet k ∈ K={1,2, n}.Then definable sequence x _iwith reference sequences x ₀correlation degree such as formula shown in (27-1);

$\mathrm{\&gamma;}(x_0,x_i)=\frac{1}{n}\underset{k=1}{\overset{n}{\&Sigma;}}\mathrm{\&gamma;}\left(x_0\right(k),x_i(k\left)\right)---(27-1)$

In formula (27-1), γ (x ₀(k), x _i(k)) be x ₀and x _iat the correlation coefficient of k point, n is sequence x _iin the element number that comprises.If γ is (x ₀, x _i) > γ (x ₀, x _j), then x is described _icompare x _jby force, that is: with reference line x ₀relatively, x _iimpact than x _jgreatly.The original index data matrix R shown in formula (1) is set up for the desired value of n energy efficiency evaluation index after m transformer station to be assessed, the standardization of each transformer station to be assessed based on gray relative analysis method in the present embodiment.Reference line R ₀corresponding ideal object is R ₀'=(r ₀₁, r ₀₂..., r _0n), wherein each element r _0jbe the optimal value of an energy efficiency evaluation index, and reference line R ₀be ideal object R ₀' standardization after result, therefore reference line R ₀in result after standardization, each element is 1, therefore and by original index data matrix R and reference line R ₀combine after standardization, can obtain the achievement data matrix S to be assessed shown in formula (2).

Make S _i=(s _i1, s _i2..., s _in), i=0,1 ..., m, then S ₀as with reference to row, calculate S successively _ieach index and reference line S ₀the correlation coefficient β of corresponding index _i(j), i=1,2 ..., m; J=1,2 ..., n.In the present embodiment, step 3) in specifically refer to the desired value and the reference line R that calculate each energy efficiency evaluation index after each transformer station to be assessed standardization according to formula (30) ₀correlation coefficient between the reference value of the energy efficiency evaluation index of middle correspondence;

${\mathrm{\&beta;}}_i\left(j\right)=\frac{\underset{i}{min}\underset{j}{min}|s_{0j}-s_{ij}|+\mathrm{\&rho;}\underset{i}{\max }\underset{j}{\max }|s_{0j}-s_{ij}|}{|s_{0j}-s_{ij}|+\mathrm{\&rho;}\underset{i}{\max }\underset{j}{\max }|s_{0j}-s_{ij}|}---\left(30\right)$

In formula (30), β _ij () is the desired value s of the jth energy efficiency evaluation index after the standardization of i-th transformer station to be assessed _ijwith reference line R ₀the reference value s of the energy efficiency evaluation index of middle correspondence _0jbetween correlation coefficient, ρ is the resolution ratio of each energy efficiency evaluation index, and the interval of the resolution ratio ρ of each energy efficiency evaluation index is ρ ∈ [0,1], and general value is ρ=0.5.

Finally, step 3) obtain the incidence coefficient matrix β of m transformer station to be assessed such as formula shown in (30-1);

$\mathrm{\&beta;}=\left[\begin{array}{cccc}{\mathrm{\&beta;}}_1\left(1\right)& {\mathrm{\&beta;}}_1\left(2\right)& ...& {\mathrm{\&beta;}}_1\left(n\right)\\ {\mathrm{\&beta;}}_2\left(1\right)& {\mathrm{\&beta;}}_2\left(2\right)& ...& {\mathrm{\&beta;}}_2\left(n\right)\\ \&CenterDot;& \&CenterDot;& & \&CenterDot;\\ \&CenterDot;& \&CenterDot;& & \&CenterDot;\\ \&CenterDot;& \&CenterDot;& & \&CenterDot;\\ {\mathrm{\&beta;}}_m\left(1\right)& {\mathrm{\&beta;}}_m\left(2\right)& ...& {\mathrm{\&beta;}}_m\left(n\right)\end{array}\right]---(30-1)$

In formula (30-1), β is the incidence coefficient matrix of m transformer station to be assessed, n energy efficiency evaluation index of each behavior transformer station to be assessed with reference line R ₀correlation coefficient between the reference value of the energy efficiency evaluation index of middle correspondence, such as β _m(n) represent the n-th energy efficiency evaluation index of m transformer station to be assessed with reference line R ₀correlation coefficient between the reference value of the energy efficiency evaluation index of middle correspondence.Sum-average arithmetic is carried out to each index degree of association of each transformer station to be assessed, namely tries to achieve the correlation degree size between each transformer station to be assessed and desirable transformer station, shown in (30-2);

$x_i=\frac{1}{n}\underset{j=1}{\overset{n}{\&Sigma;}}{\mathrm{\&beta;}}_i\left(j\right)---(30-2)$

In formula (30-2), x _irepresent the correlation degree size between i-th transformer station to be assessed and desirable transformer station, n is the total quantity of energy efficiency evaluation index, β _ij () is the desired value s of the jth energy efficiency evaluation index after the standardization of i-th transformer station to be assessed _ijwith reference line R ₀the reference value s of the energy efficiency evaluation index of middle correspondence _0jbetween correlation coefficient.On this basis, according to x _ithe value good and bad degree that can treat evaluation object evaluate.The x of certain evaluation object _ibe worth larger, illustrate that this evaluation object is more close to ideal object, thus can illustrate that this object is more outstanding in all objects to be assessed.And analytical hierarchy process, for determining each index weights, when calculating resolution ratio ρ and each evaluation object degree of association, consider the impact of index weights, thus improve the confidence level of gray relative analysis method.

In the present embodiment, step 4) detailed step comprise:

4.1) desired value for n energy efficiency evaluation index sets up target layers structure based on the classification of energy efficiency evaluation index, as shown in Figure 2, the superiors of mark hierarchical structure are the target indicator layer that substation energy efficiency assessment result is formed, orlop is the scheme indicator layer that n energy efficiency evaluation index is formed, and is the intermediary outcomes layer that the classification of energy efficiency evaluation index is formed between scheme indicator layer and target indicator layer; Energy efficiency evaluation index under each being classified, as one group of energy efficiency evaluation index, selects one group of energy efficiency evaluation index as current group of energy efficiency evaluation index; According to the target layers structure set up, the membership of last layer and lower one deck can be determined.

4.2) determine the important coefficient of each group energy efficiency evaluation index another group energy efficiency evaluation index relative in current group of energy efficiency evaluation index, obtain the judgment matrix of the current group of energy efficiency evaluation index be made up of all important coefficients; Be supplied in the present embodiment and weight is set to distinguish the significance level of different index, great majority are comprised to the challenge of multi objective, distinguish the significance level of index and be not easy, analytical hierarchy process adopts the method compared between two to determine the important coefficient of each group energy efficiency evaluation index another group energy efficiency evaluation index relative in current group of energy efficiency evaluation index.Known see Fig. 2, the desired value of n the energy efficiency evaluation index that the present embodiment is specified comprises electrical equipment information, system running state, building structure information, environmental impact factor four classification altogether, therefore four corresponding four groups of energy efficiency evaluation indexs of classification, obtain judgment matrix A, B, C, D of the current group of energy efficiency evaluation index be made up of all important coefficients altogether, for judgment matrix A, its form is such as formula shown in (30-3);

$A=\left[\begin{array}{ccccc}{a}_{11}& a_{12}& a_{13}& ...& a_{1n}\\ a_{21}& a_{22}& a_{23}& ...& a_{2n}\\ a_{31}& a_{32}& a_{33}& ...& a_{3n}\\ \&CenterDot;& \&CenterDot;& \&CenterDot;& \&CenterDot;& \&CenterDot;\\ \&CenterDot;& \&CenterDot;& \&CenterDot;& \&CenterDot;& \&CenterDot;\\ \&CenterDot;& \&CenterDot;& \&CenterDot;& \&CenterDot;& \&CenterDot;\\ a_{n1}& a_{n2}& a_{n3}& ...& a_{nn}\end{array}\right]---(30-3)$

In formula (30-3), a _ijrepresent index x _irelative to x _jimportance, and a _ji=1/a _ij.

In the present embodiment, the concrete value of judgment matrix A is as shown in table 3; The concrete value of judgment matrix B is as shown in table 4; The concrete value of judgment matrix C is as shown in table 5; The concrete value of judgment matrix D is as shown in table 6; The judgment matrix value of the current group of energy efficiency evaluation index be made up of all important coefficients is as shown in table 7.

Table 3: the value table of judgment matrix A;

A	A1	A2	A3	A4	A5	A6	A7	A8	A9	A10
											A1	1	1	3	9	5	5	6	9	3	6
A2	1	1	4	7	5	6	8	9	2	4

A3	1/3	1/4	1	4	1	3	4	6	1/5	3
											A4	1/9	1/7	1/4	1	1/5	1/3	1	2	1/7	1
A5	1/5	1/5	1	5	1	4	6	6	1/2	5
											A6	1/5	1/6	1/3	3	1/4	1	4	6	1/3	3
A7	1/6	1/8	1/4	1	1/6	1/4	1	1/2	1/7	1
											A8	1/9	1/9	1/6	1/2	1/6	1/6	2	1	1/7	2
A9	1/3	1/2	5	7	2	3	7	7	1	6
											A10	1/6	1/4	1/3	1	1/5	1/3	1	1/2	1/6	1

Table 4: the value table of judgment matrix B;

B	B1	B2	B3	B4	B5	B6
							B1	1	1	4	4	8	8
B2	1	1	4	4	8	8
							B3	1/4	1/4	1	1	4	4
B4	1/4	1/4	1	1	4	4
							B5	1/8	1/8	1/4	1/4	1	1
B6	1/8	1/8	1/4	1/4	1	1

Table 5: the value table of judgment matrix C;

C	C1	C2	C3	C4	C5
						C1	1	1/7	1/9	1/5	1/3
C2	7	1	1/3	3	5
						C3	9	3	1	5	7
C4	5	1/3	1/5	1	3
						C5	3	1/5	1/7	1/3	1

Table 6: the value table of judgment matrix D;

D	D1	D2
			D1	1	1
D2	1	1

Table 7: the value table of the judgment matrix of the current group of energy efficiency evaluation index be made up of all important coefficients;

Indicator layer	A	B	C	D
					A	1	3	5	9
B	1/3	1	3	7
					C	1/5	1/3	1	5
D	1/9	1/7	1/5	1

4.3) the eigenvalue of maximum λ of the judgment matrix of current group of energy efficiency evaluation index is calculated _max, be 1 for benchmark with the weight sum of all energy efficiency evaluation indexs of current group of energy efficiency evaluation index, by eigenvalue of maximum λ _maxcorresponding proper vector is normalized, and obtains the relative weighting of each energy efficiency evaluation index relative the superiors index in current group of energy efficiency evaluation index; For judgment matrix A, because judgment matrix A describes the relative importance between each index, from matrix angle analysis, in the eigenwert of judgment matrix A and proper vector, namely contain the relative importance information between each index.Therefore, the eigenvalue of maximum λ of desirable judgment matrix A _maxcharacteristic of correspondence vector is as the relative weighting of each level index.Meanwhile, consider that each layer index weights sum should be 1, therefore need by eigenvalue of maximum λ _maxcorresponding proper vector is normalized.For judgment matrix A, there is formula (30-4);

A·β＝λ·β(30-4)

In formula (30-4), A is judgment matrix; λ is characteristic root; β is proper vector.Therefore, eigenvalue of maximum λ is got _maxand characteristic of correspondence vector, and be normalized, what obtain is exactly the weighted value of lower floor's index for upper strata index.

4.4) judge whether that all groups of energy efficiency evaluation indexs are disposed, if be not yet disposed, then select next group energy efficiency evaluation index as current group of energy efficiency evaluation index, and redirect performs step 4.2); Be disposed else if, then redirect performs step 4.5);

4.5) relative weighting of each energy efficiency evaluation index relative the superiors index in all groups of energy efficiency evaluation indexs is combined, obtain the relative weighting matrix ω be made up of the relative weighting of n energy efficiency evaluation index relative the superiors index, wherein n energy efficiency evaluation index relative the superiors index relative weighting and be 1.

In the present embodiment, step 4.2) in when to determine in current group of energy efficiency evaluation index the importance of each group energy efficiency evaluation index another group energy efficiency evaluation index relatively, specifically refer to the importance determining another group energy efficiency evaluation index relatively of each group energy efficiency evaluation index in current group of energy efficiency evaluation index according to nine grades of scaling laws.

In the present embodiment, when determining the important coefficient of each group energy efficiency evaluation index another group energy efficiency evaluation index relative in current group of energy efficiency evaluation index according to nine grades of scaling laws, each group energy efficiency evaluation index relatively another importance organizing energy efficiency evaluation index comprises nine grades of significance levels, nine grades of significance levels comprise identical, important a little, obviously in important, strongly important, extremely important Pyatyi significance level and five kinds of Pyatyi significance levels in the middle of arbitrary neighborhood Pyatyi significance level intergrade significance level, and the important coefficient a that nine grades of significance levels are corresponding _ij1 ~ 9 is followed successively by according to importance degree, specifically as shown in table 8.

Table 8: nine grades of scaling law reference tables.

a ij	Explanation	a ij	Explanation
				1	Index x iWith x jImportance is identical	7	Index x iThan index x jStrongly important
3	Index x iThan index x jImportant a little	9	Index x iThan index x jExtremely important
				5	Index x iThan index x jObviously important	2、4、6、8	The intermediate state of corresponding above adjacent judgement

For judgment matrix A, for the judgment matrix A drawn with nine grades of scaling laws, due to the impact of evaluator individual preference, may cause occurring inconsistent situation between each index importance, such as index A1 is more important than index A2, index A2 is more important than index A3, and index A3 is more important than index A1.When index is less, this inconsistency can by checking and avoiding one by one; When index is more, then need by carrying out consistency desired result to find to judgment matrix and getting rid of problems.Therefore in the present embodiment, step 4.3) also comprise the step of the judgment matrix of current group of energy efficiency evaluation index being carried out to consistency desired result, and if verification pass through, then redirect performs step 4.4), otherwise redirect performs step 4.2 again).

In the present embodiment, the detailed step judgment matrix of current group of energy efficiency evaluation index being carried out to consistency desired result comprises:

4.3.1) judge the exponent number of the judgment matrix of current group of energy efficiency evaluation index, if the exponent number of the judgment matrix of current group of energy efficiency evaluation index is no more than 2, then directly redirect performs step 4.4), otherwise redirect performs next step;

4.3.2) the inconsistency degree quantitative index of the judgment matrix of current group of energy efficiency evaluation index is calculated according to formula (28);

$CI=\frac{{\mathrm{\&lambda;}}_{max}-n}{n-1}---\left(28\right)$

In formula (28), CI is the inconsistency degree quantitative index of the judgment matrix of current group of energy efficiency evaluation index, λ _maxfor the eigenvalue of maximum of the judgment matrix of current group of energy efficiency evaluation index, n is the quantity of all energy efficiency evaluation indexs;

4.3.3) search default Aver-age Random Consistency Index and energy efficiency evaluation index schedule of quantities according to the quantity n of energy efficiency evaluation index, obtain the Aver-age Random Consistency Index value RI of the judgment matrix of current group of energy efficiency evaluation index;

4.3.4) the inconsistency degree quantitative index modified value CR of the judgment matrix of current group of energy efficiency evaluation index is calculated according to formula (29);

$CR=\frac{CI}{RI}---\left(29\right)$

In formula (29), CR is the inconsistency degree quantitative index modified value of the judgment matrix of current group of energy efficiency evaluation index, CI is the inconsistency degree quantitative index of the judgment matrix of current group of energy efficiency evaluation index, and RI is the Aver-age Random Consistency Index value of the judgment matrix of current group of energy efficiency evaluation index.

In the present embodiment, the corresponding relation of the Aver-age Random Consistency Index value RI of the judgment matrix of current group of energy efficiency evaluation index and the quantity n of energy efficiency evaluation index is as shown in table 9;

Table 9: the corresponding relation value table of the quantity n of Aver-age Random Consistency Index RI and energy efficiency evaluation index;

n	1	2	3	4	5	6	7	8	9	10	11
												RI	0	0	0.58	0.9	1.12	1.24	1.32	1.41	1.45	1.49	1.51

4.3.5) judge that the inconsistency degree quantitative index modified value CR of the judgment matrix of current group of energy efficiency evaluation index is less than 0.1 and whether sets up, if set up, judge that carrying out consistency desired result to the judgment matrix of current group of energy efficiency evaluation index passes through, otherwise judge that carrying out consistency desired result to the judgment matrix of current group of energy efficiency evaluation index does not pass through.

The above is only the preferred embodiment of the present invention, protection scope of the present invention be not only confined to above-described embodiment, and all technical schemes belonged under thinking of the present invention all belong to protection scope of the present invention.It should be pointed out that for those skilled in the art, some improvements and modifications without departing from the principles of the present invention, these improvements and modifications also should be considered as protection scope of the present invention.

Claims (10)

1., based on a substation energy efficiency evaluating method for analytical hierarchy process, it is characterized in that step comprises:

1) for m transformer station to be assessed, the desired value of n the energy efficiency evaluation index of specifying is obtained respectively;

2) for m transformer station to be assessed, the desired value of n energy efficiency evaluation index of each transformer station to be assessed is set up the original index data matrix R shown in formula (1) based on gray relative analysis method; Obtain the capable R of optimal value generating reference of each energy efficiency evaluation index respectively ₀, by original index data matrix R and reference line R ₀in each data carry out standardization and combination obtain the achievement data matrix S to be assessed shown in formula (2);

$R=\left[\begin{array}{cccc}{r}_{11}& r_{12}& ...& r_{1n}\\ r_{21}& r_{22}& ...& r_{2n}\\ \&CenterDot;& \&CenterDot;& & \&CenterDot;\\ \&CenterDot;& \&CenterDot;& & \&CenterDot;\\ \&CenterDot;& \&CenterDot;& & \&CenterDot;\\ r_{m1}& r_{m2}& ...& r_{mn}\end{array}\right]---\left(1\right)$

In formula (1), R is the original index data matrix set up, r _mnit is the desired value of m transformer station n-th to be assessed energy efficiency evaluation index;

$S=\left[\begin{array}{cccc}1& 1& ...& 1\\ s_{11}& s_{12}& ...& s_{1n}\\ s_{21}& s_{22}& ...& s_{2n}\\ \&CenterDot;& \&CenterDot;& & \&CenterDot;\\ \&CenterDot;& \&CenterDot;& & \&CenterDot;\\ \&CenterDot;& \&CenterDot;& & \&CenterDot;\\ s_{m1}& s_{m2}& ...& s_{mn}\end{array}\right]---\left(2\right)$

In formula (2), S is the achievement data matrix to be assessed obtained, the first behavior reference line R in achievement data matrix S to be assessed ₀result after standardization, s _mnbe m the desired value of transformer station n-th energy efficiency evaluation index to be assessed after standardization;

3) for achievement data matrix S to be assessed, desired value and the reference line R of each energy efficiency evaluation index after each transformer station to be assessed standardization is calculated respectively ₀correlation coefficient between the reference value of the energy efficiency evaluation index of middle correspondence, obtains the incidence coefficient matrix β of m transformer station to be assessed;

4) desired value for described n energy efficiency evaluation index sets up target layers structure based on the classification of energy efficiency evaluation index, based on the relative weighting of each energy efficiency evaluation index relative the superiors index of analytical hierarchy process agriculture products hierarchical structure bottom, obtain the relative weighting matrix ω be made up of the relative weighting of n energy efficiency evaluation index relative the superiors index;

5) calculate the integrated assessment result matrix of m transformer station to be assessed according to formula (3), and wherein the energy efficiency evaluation result of i-th transformer station to be assessed such as formula shown in (4);

X＝β·ω(3)

In formula (3), X is the integrated assessment result matrix of m transformer station to be assessed, and β is the incidence coefficient matrix of m transformer station to be assessed, and ω is the relative weighting matrix that the relative weighting of n energy efficiency evaluation index is formed;

$x_i=\&lsqb;{\mathrm{\&beta;}}_i\left(1\right),{\mathrm{\&beta;}}_i\left(2\right),...,{\mathrm{\&beta;}}_i\left(n\right)\&rsqb;\&CenterDot;\left[\begin{array}{c}{\mathrm{\&omega;}}_1\\ {\mathrm{\&omega;}}_2\\ \&CenterDot;\\ \&CenterDot;\\ \&CenterDot;\\ {\mathrm{\&omega;}}_n\end{array}\right]---\left(4\right)$

In formula (4), x _ibe the energy efficiency evaluation result of i-th transformer station to be assessed, β _in () is desired value and the reference line R of the n-th energy efficiency evaluation index of i-th transformer station to be assessed ₀correlation coefficient between the reference value of the energy efficiency evaluation index of middle correspondence, ω _nit is the relative weighting of the target indicator layer of the n-th energy efficiency evaluation index relative indicatrix hierarchical structure.

2. the substation energy efficiency evaluating method based on analytical hierarchy process according to claim 1, it is characterized in that, described step 1) in the desired value of n energy efficiency evaluation index of specifying comprise electrical equipment information altogether, system running state, building structure information, environmental impact factor four classification, energy efficiency evaluation index under described electrical equipment information classification comprises nature cooling main transformer number of units accounting A1, air-conditioning Energy Efficiency Ratio A2, energy-saving lamp general power and light fixture general power accounting A3, switch board temp and humidity regulator energy-saving index A4, air-conditioning Temperature and Humidity Control mode energy efficiency indexes A5, illumination Temperature and Humidity Control mode energy efficiency indexes A6, ventilation Temperature and Humidity Control mode energy efficiency indexes A7, switch board Temperature and Humidity Control mode energy efficiency indexes A8, stand with Variant number energy saving index A9, the horizontal index A10 of ageing equipment, energy efficiency evaluation index under described system running state classification comprising load shape coefficient B1, main transformer economical operation index B2, busbar voltage deviation B3, power factor B4, alternate laod unbalance coefficient B 5, standing with becoming economical operation index B6, described building structure information comprises building orientation index C1, shape coefficient of architecture C2, heat transfer across wall index C3, the heat insulation index C4 of window, window sunshade index C5, described environmental impact factor comprises heating degree number D1 and air-conditioning and to subsist several D2.

3. the substation energy efficiency evaluating method based on analytical hierarchy process according to claim 2, it is characterized in that, described step 1) in classification n energy efficiency evaluation index desired value in, the computing function expression formula of the described main transformer number of units of cooling naturally accounting A1 is such as formula shown in (5);

A1＝n/N(5)

In formula (5), A1 is nature cooling main transformer number of units accounting, and n is the main transformer number of units of nature cooling; N is the total number of units of main transformer of transformer station;

The computing function expression formula of described air-conditioning Energy Efficiency Ratio A2 is such as formula shown in (6);

$A2=\underset{i=1}{\overset{n}{\&Sigma;}}{\mathrm{EER}}_i\&CenterDot;\frac{P_{out.i}}{P_{out.total}}---\left(6\right)$

In formula (6), n is the total number of units of transformer station's air-conditioning; EER _ifor i-th air-conditioning Energy Efficiency Ratio of transformer station; P _out.totalfor substation has air-conditioning total value to customize cold; P _out.ibe i-th specified refrigerating capacity of air-conditioning; Wherein i-th air-conditioning trapped energy theory of transformer station _icomputing function expression formula such as formula shown in (7);

EER＝P _out/P _in(7)

In formula (7), EER is the air-conditioning Energy Efficiency Ratio of transformer station's air-conditioning, P _outfor the refrigerating capacity of air-conditioning under rated input power, P _infor air-conditioning rated input power;

The computing function expression formula of described energy-saving lamp general power and light fixture general power accounting A3 is such as formula shown in (8);

$A3=\frac{\underset{i=1}{\overset{m}{\&Sigma;}}{P}_{es.i}}{\underset{j=1}{\overset{n}{\&Sigma;}}{P}_j},(m\&le;n)---\left(8\right)$

In formula (8), A3 is energy-saving lamp general power and light fixture general power accounting, and n is transformer station's light fixture sum, and m is transformer station's energy saving lamp sum, P _jfor the power of transformer station's jth light fixture; P _es.ifor the power of transformer station's i-th energy saving lamp;

The computing function expression formula of described switch board temp and humidity regulator energy-saving index A4 is such as formula shown in (9);

$A4=\frac{0.5n_2+n_3}{n_1+n_2+n_3}---\left(9\right)$

In formula (9), be A4 switch board temp and humidity regulator energy-saving index, n ₁for the number of heating plate in switch board temp and humidity regulator, n ₂for the number with fan heating plate in switch board temp and humidity regulator, n ₃for the number of heat exchanger in switch board temp and humidity regulator;

The computing function expression formula of described air-conditioning Temperature and Humidity Control mode energy efficiency indexes A5, illumination Temperature and Humidity Control mode energy efficiency indexes A6, ventilation Temperature and Humidity Control mode energy efficiency indexes A7, switch board Temperature and Humidity Control mode energy efficiency indexes A8 is such as formula shown in (10);

$Ai=\frac{P_{auto}}{P},(i=5,6,7,8)---\left(10\right)$

In formula (10), Ai (i=5,6,7,8) calculated air-conditioning Temperature and Humidity Control mode energy efficiency indexes A5 or illumination Temperature and Humidity Control mode energy efficiency indexes A6 or ventilation Temperature and Humidity Control mode energy efficiency indexes A7 or switch board Temperature and Humidity Control mode energy efficiency indexes A8 is represented; P _autorepresent that control mode is the equipment general power automatically controlled, P represents the equipment general power of all control modes;

The computing function expression formula of described station Variant number energy saving index A9 is such as formula shown in (11);

$A9=\frac{\underset{i=1}{\overset{n}{\&Sigma;}}{\mathrm{SCORE}}_i}{n}---\left(11\right)$

In formula (11), A9 is station Variant number energy saving, and n is transformer station's distribution transforming number of units; SCORE _ithe energy saving score that i-th station determined for composing submeter according to the station Variant number preset becomes;

The computing function expression formula of the horizontal index A10 of described ageing equipment is such as formula shown in (12);

$A10=\frac{\underset{i=1}{\overset{n}{\&Sigma;}}(1+N_i/80)}{n}---\left(12\right)$

In formula (12), A10 is the horizontal index of ageing equipment, and n is the quantity of transformer station's main electric power equipment, and described main electric power equipment comprises main transformer, stands with change, shnt capacitor, isolating switch and disconnecting link, N _ibe the operation year number of i-th main electric power equipment, and when the operation year number of main electric power equipment is greater than 20 years N _ivalue get 32;

The computing function expression formula of described load shape coefficient B1 is such as formula shown in (13);

$B1=\frac{I_{eff}}{I_{ar}}---\left(13\right)$

In formula (13), B1 is load shape coefficient, I _efffor the Payload of fixed time section, I _arfor the average load of fixed time section;

The computing function expression formula of described main transformer economical operation index B2 is such as formula shown in (14);

$B2=\left\{\begin{array}{c}\mathrm{\&beta;},\mathrm{\&beta;}<{\mathrm{\&beta;}}_J\\ 1.33{\mathrm{\&beta;}}_J^2+\frac{0.75-\mathrm{\&beta;}}{0.75-{\mathrm{\&beta;}}_J}({\mathrm{\&beta;}}_J-1.33{\mathrm{\&beta;}}_J^2),{\mathrm{\&beta;}}_J\&le;\mathrm{\&beta;}<0.75\\ {\mathrm{\&beta;}}_J^2+\frac{1-\mathrm{\&beta;}}{0.25}(1.33{\mathrm{\&beta;}}_J^2-{\mathrm{\&beta;}}_J^2),0.75\&le;\mathrm{\&beta;}<1\end{array}\right.---\left(14\right)$

In formula (14), B2 is main transformer economical operation index, and β is transformer Rate of average load; β _jfor transformer economic operation point load factor;

The computing function expression formula of described busbar voltage deviation B3 is such as formula shown in (15);

$B3=\underset{j=1}{\overset{m}{\&Sigma;}}{\mathrm{\&Delta;u}}_{Tj}/m---\left(15\right)$

In formula (15), B3 is busbar voltage deviation, and m represents transformer substation voltage number of degrees, △ u _tjrepresent the average of jth electric pressure bus integral point voltage deviation absolute value within the statistics phase, △ u _tjcomputing function expression formula such as formula shown in (16);

${\mathrm{\&Delta;u}}_T=\underset{i=1}{\overset{n}{\&Sigma;}}\frac{|u_i-u_N|}{u_N}/n---\left(16\right)$

In formula (16), △ u _tfor the average of certain electric pressure bus integral point voltage deviation absolute value within the statistics phase, u _ifor the integral point magnitude of voltage of bus, n is integral point number in the statistics phase, u _nfor bus rated voltage;

The computing function expression formula of described power factor B4 is such as formula shown in (17);

$B4=\underset{i=1}{\overset{n}{\&Sigma;}}{\mathrm{\&eta;}}_i\&CenterDot;\frac{P_i}{P_{total}}---\left(17\right)$

In formula (17), B4 is power factor, and n is the transforming plant main transformer number of windings, η _ifor the power factor of main transformer i side, P _ifor the active power of main transformer i side, P _totalfor the total active power of transformer station, wherein the power factor η of main transformer i side _ifor the active-power P of main transformer i side _idivided by the applied power S of main transformer i side _iobtain;

The computing function expression formula of described alternate laod unbalance coefficient B 5 is such as formula shown in (18) or formula (19);

In formula (18), B5 is alternate laod unbalance coefficient, S _mfor the maximum load power in three-phase, S _xfor the minimum load power in three-phase, for three-phase average load power;

In formula (19), B5 is alternate laod unbalance coefficient, I _mfor the maximum load current in three-phase, I _xfor the minimum load current in three-phase, for three-phase average load current;

Described station becomes the computing function expression formula of economical operation index B6 such as formula shown in (20);

$B6=\left\{\begin{array}{c}\mathrm{\&beta;},\mathrm{\&beta;}<{\mathrm{\&beta;}}_J\\ 1.33{\mathrm{\&beta;}}_J^2+\frac{0.75-\mathrm{\&beta;}}{0.75-{\mathrm{\&beta;}}_J}({\mathrm{\&beta;}}_J-1.33{\mathrm{\&beta;}}_J^2),{\mathrm{\&beta;}}_J\&le;\mathrm{\&beta;}<0.75\\ {\mathrm{\&beta;}}_J^2+\frac{1-\mathrm{\&beta;}}{0.25}(1.33{\mathrm{\&beta;}}_J^2-{\mathrm{\&beta;}}_J^2),0.75\&le;\mathrm{\&beta;}<1\end{array}\right.---\left(20\right)$

In formula (20), B6 is station change economical operation index, and β is transformer Rate of average load; β _jfor transformer economic operation point load factor;

When the building orientation of transformer station is north-south, the value of building orientation index C1 is taken as 1, when the building orientation of transformer station is East and West direction, the value of building orientation index C1 is taken as 0, when the building orientation of transformer station is between north-south and East and West direction, according to towards angle being the value that building orientation index C1 composes between 0 and 1;

The computing function expression formula of described shape coefficient of architecture C2 is such as formula shown in (21);

$C2=\frac{F_0}{V_0}---\left(21\right)$

In formula (21), C2 is shape coefficient of architecture, F ₀for the external surface area of building; V ₀for architectural volume;

The computing function expression formula of described heat transfer across wall index C3 is such as formula shown in (22);

$C3=\frac{1}{R_i+(R_1+R_2+...+R_n)+R_e}---\left(22\right)$

In formula (22), C3 is heat transfer across wall index, R _ifor the heat exchange of building enclosure inside surface hinders, R _efor the heat exchange of building enclosure outside surface hinders, R _nfor the thermal resistance of the n-th layer material of building enclosure;

The heat insulation index C4 of described window specifically refers to the heat-proof quality score of the construction window glass determined according to the heat insulation index grade form of window preset;

Described window sunshade index C5 specifically pointer to the pre-set constant value without sunshade measure, indoor sunshade measure, outdoor sunshade measure three kinds of modes of transformer station;

The computing function expression formula of described heating degree number D1 is such as formula shown in (23);

$D1=\underset{i=1}{\overset{n}{\&Sigma;}}(t_R-t_{m,i})---\left(23\right)$

In formula (23), D1 is heating degree number, and n is days of heating period or calculates number of days, t _rfor indoor reference temperature, t _m,ifor the outdoor daily mean temperature of Heating Period;

Described air-conditioning subsists the computing function expression formula of several D2 such as formula shown in (24);

$D2=\underset{i=1}{\overset{n}{\&Sigma;}}(t_{m,i}-t_R)---\left(24\right)$

In formula (24), D2 is that air-conditioning is subsisted number, and n is days of heating period or calculates number of days, t _rfor indoor reference temperature, t _m,ifor the outdoor daily mean temperature of Heating Period.

4. the substation energy efficiency evaluating method based on analytical hierarchy process according to claim 3, it is characterized in that, described step 2) Plays process specifically refers to: the type judging each energy efficiency evaluation index, if the type of current energy efficiency evaluation index is the direct index that value is the bigger the better, then formula (25) is adopted to calculate the standardization result of the desired value of current energy efficiency evaluation index; If the type of the current energy efficiency evaluation index negative index that to be value the smaller the better, then adopt formula (26) to calculate the standardization result of the desired value of current energy efficiency evaluation index; If the type of the current energy efficiency evaluation index moderate osculant index of being advisable that is value, then adopt formula (27) to calculate the standardization result of the desired value of current energy efficiency evaluation index;

$z=\frac{x_{ij}-x_{min}}{x_{max}-x_{min}}---\left(25\right)$

$z=\frac{x_{max}-x_{ij}}{x_{max}-x_{min}}---\left(26\right)$

$z=\left\{\begin{array}{cc}\frac{x_{ij}-x_{\min }}{U_1-x_{\min }}& x_{\min }<x_{ij}<U_1\\ 1& U_1<x_{ij}<U_2\\ \frac{x_{\max }-x_{ij}}{x_{\max }-U_2}& U_2<x_{ij}<x_{\max }\end{array}\right.---\left(27\right)$

In formula (25), (26) and (27), z is the standardization result of the desired value of current energy efficiency evaluation index, x _minfor the minimum value of the desired value of current energy efficiency evaluation index, x _maxfor the maximal value of the desired value of current energy efficiency evaluation index, x _ijfor the desired value of current energy efficiency evaluation index; [U ₁, U ₂] for preset current energy efficiency evaluation index optimal zone between.

5. the substation energy efficiency evaluating method based on analytical hierarchy process according to claim 4, it is characterized in that, described step 3) in specifically refer to the desired value and the reference line R that calculate each energy efficiency evaluation index after each transformer station to be assessed standardization according to formula (30) ₀correlation coefficient between the reference value of the energy efficiency evaluation index of middle correspondence;

${\mathrm{\&beta;}}_i\left(j\right)=\frac{\underset{i}{min}\underset{j}{min}|s_{0j}-s_{ij}|+\mathrm{\&rho;}\underset{i}{\max }\underset{j}{\max }|s_{0j}-s_{ij}|}{|s_{0j}-s_{ij}|+\mathrm{\&rho;}\underset{i}{\max }\underset{j}{\max }|s_{0j}-s_{ij}|}---\left(30\right)$

In formula (30), β _ij () is the desired value s of the jth energy efficiency evaluation index after the standardization of i-th transformer station to be assessed _ijwith reference line R ₀the reference value s of the energy efficiency evaluation index of middle correspondence _0jbetween correlation coefficient, ρ is the resolution ratio of each energy efficiency evaluation index, and the interval of the resolution ratio ρ of each energy efficiency evaluation index is ρ ∈ [0,1].

6. according to the substation energy efficiency evaluating method based on analytical hierarchy process in Claims 1 to 5 described in any one, it is characterized in that, described step 4) detailed step comprise:

4.1) desired value for described n energy efficiency evaluation index sets up target layers structure based on the classification of energy efficiency evaluation index, the superiors of described mark hierarchical structure are the target indicator layer that substation energy efficiency assessment result is formed, orlop is the scheme indicator layer that n energy efficiency evaluation index is formed, and is the intermediary outcomes layer that the classification of energy efficiency evaluation index is formed between described scheme indicator layer and target indicator layer; Energy efficiency evaluation index under each being classified, as one group of energy efficiency evaluation index, selects one group of energy efficiency evaluation index as current group of energy efficiency evaluation index;

4.2) determine the important coefficient of each group energy efficiency evaluation index another group energy efficiency evaluation index relative in current group of energy efficiency evaluation index, obtain the judgment matrix of the current group of energy efficiency evaluation index be made up of all important coefficients;

4.3) the eigenvalue of maximum λ of the judgment matrix of current group of energy efficiency evaluation index is calculated _max, be 1 for benchmark with the weight sum of all energy efficiency evaluation indexs of current group of energy efficiency evaluation index, by described eigenvalue of maximum λ _maxcorresponding proper vector is normalized, and obtains the relative weighting of each energy efficiency evaluation index relative the superiors index in current group of energy efficiency evaluation index;

4.4) judge whether that all groups of energy efficiency evaluation indexs are disposed, if be not yet disposed, then select next group energy efficiency evaluation index as current group of energy efficiency evaluation index, and redirect performs step 4.2); Be disposed else if, then redirect performs step 4.5);

4.5) relative weighting of each energy efficiency evaluation index relative the superiors index in all groups of energy efficiency evaluation indexs is combined, obtain the relative weighting matrix ω be made up of the relative weighting of n energy efficiency evaluation index relative the superiors index, wherein n energy efficiency evaluation index relative the superiors index relative weighting and be 1.

7. the substation energy efficiency evaluating method based on analytical hierarchy process according to claim 6, it is characterized in that, described step 4.2) in when to determine in current group of energy efficiency evaluation index the importance of each group energy efficiency evaluation index another group energy efficiency evaluation index relatively, specifically refer to the importance determining another group energy efficiency evaluation index relatively of each group energy efficiency evaluation index in current group of energy efficiency evaluation index according to nine grades of scaling laws.

8. the substation energy efficiency evaluating method based on analytical hierarchy process according to claim 7, it is characterized in that, described when to determine in current group of energy efficiency evaluation index the important coefficient of each group energy efficiency evaluation index another group energy efficiency evaluation index relatively according to nine grades of scaling laws, each group energy efficiency evaluation index relatively another importance organizing energy efficiency evaluation index comprises nine grades of significance levels, described nine grades of significance levels comprise identical, important a little, obviously important, strongly important, intergrade significance level in extremely important Pyatyi significance level and described five kinds of Pyatyi significance levels in the middle of arbitrary neighborhood Pyatyi significance level, and important coefficient corresponding to described nine grades of significance levels is followed successively by 1 ~ 9 according to importance degree.

9. the substation energy efficiency evaluating method based on analytical hierarchy process according to claim 8, it is characterized in that, described step 4.3) also comprise the step of the judgment matrix of current group of energy efficiency evaluation index being carried out to consistency desired result, if and verification is passed through, then redirect performs step 4.4), otherwise redirect performs step 4.2 again).

10. the substation energy efficiency evaluating method based on analytical hierarchy process according to claim 9, is characterized in that, the detailed step that the described judgment matrix to current group of energy efficiency evaluation index carries out consistency desired result comprises:

4.3.1) judge the exponent number of the judgment matrix of current group of energy efficiency evaluation index, if the exponent number of the judgment matrix of current group of energy efficiency evaluation index is no more than 2, then directly redirect performs step 4.4), otherwise redirect performs next step;

4.3.2) the inconsistency degree quantitative index of the judgment matrix of current group of energy efficiency evaluation index is calculated according to formula (28);

$CI=\frac{{\mathrm{\&lambda;}}_{max}-n}{n-1}---\left(28\right)$

In formula (28), CI is the inconsistency degree quantitative index of the judgment matrix of current group of energy efficiency evaluation index, λ _maxfor the eigenvalue of maximum of the judgment matrix of current group of energy efficiency evaluation index, n is the quantity of all energy efficiency evaluation indexs;

4.3.3) search default Aver-age Random Consistency Index and energy efficiency evaluation index schedule of quantities according to the quantity n of energy efficiency evaluation index, obtain the Aver-age Random Consistency Index value RI of the judgment matrix of current group of energy efficiency evaluation index;

4.3.4) the inconsistency degree quantitative index modified value CR of the judgment matrix of current group of energy efficiency evaluation index is calculated according to formula (29);

$CR=\frac{CI}{RI}---\left(29\right)$

In formula (29), CR is the inconsistency degree quantitative index modified value of the judgment matrix of current group of energy efficiency evaluation index, CI is the inconsistency degree quantitative index of the judgment matrix of current group of energy efficiency evaluation index, and RI is the Aver-age Random Consistency Index value of the judgment matrix of current group of energy efficiency evaluation index;

4.3.5) judge that the inconsistency degree quantitative index modified value CR of the judgment matrix of current group of energy efficiency evaluation index is less than 0.1 and whether sets up, if set up, judge that carrying out consistency desired result to the judgment matrix of current group of energy efficiency evaluation index passes through, otherwise judge that carrying out consistency desired result to the judgment matrix of current group of energy efficiency evaluation index does not pass through.