CN104700317A - Assessment method and device of grid partition power supply mode - Google Patents

Assessment method and device of grid partition power supply mode Download PDF

Info

Publication number
CN104700317A
CN104700317A CN201310654861.7A CN201310654861A CN104700317A CN 104700317 A CN104700317 A CN 104700317A CN 201310654861 A CN201310654861 A CN 201310654861A CN 104700317 A CN104700317 A CN 104700317A
Authority
CN
China
Prior art keywords
matrix
area division
sub area
index
powering mode
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201310654861.7A
Other languages
Chinese (zh)
Inventor
杨楠
刘景延
袁清芳
马恒瑞
孙健
陈璟昊
王海云
李亚东
高明伟
李亮玉
周宝玉
刘文颖
张东英
林俐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
State Grid Beijing Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Beijing Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, State Grid Beijing Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201310654861.7A priority Critical patent/CN104700317A/en
Publication of CN104700317A publication Critical patent/CN104700317A/en
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply

Landscapes

  • Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses an assessment method and device of a grid partition power supply mode. The assessment method includes calculating grades of indexes of the grid partition power supply mode to obtain a standard assessment matrix; calculating weights of the indexes to obtain a weight matrix; according to the weight matrix and the standard assessment matrix, calculating a comprehensive assessment matrix of the grid partition power supply mode. By the use of the assessment method and device, the problem that the grid partition power supply plan fails to be comprehensively assessed is solved, and accordingly an optimal grid partition plan is reasonably selected.

Description

The appraisal procedure of sub area division powering mode and device
Technical field
The present invention relates to Electric Power Network Planning field, in particular to a kind of appraisal procedure and device of sub area division powering mode.
Background technology
It is the safe operation and the reliable power supply major action that ensure important big city electrical network that sub area division runs, and is also the developing direction of the following method of operation.Sub area division runs this decision process to be needed to consider the factor such as short circuit current level, power supply safety reliability before and after power transmitting capability, subregion, should, from each side comprehensive evaluation pros and cons of operation of power networks, select optimally partitioned mutual for mode.This decision process is multiple goal, a unstructuredness optimization problem, be difficult to be described as complete, unified mathematical model with mathematical linguistics, the uncoupling loop scheme that present stage is commonly used is generally rely on engineering experience to judge, or be optimal objective with single factors, determine final uncoupling loop scheme, and then manage to allow other factor meet the condition of this uncoupling loop scheme.
For cannot the problem of comprehensive assessment sub area division power supply plan in correlation technique, at present effective solution be not yet proposed.
Summary of the invention
Fundamental purpose of the present invention is the appraisal procedure and the device that provide a kind of sub area division powering mode, cannot the problem of comprehensive assessment sub area division power supply plan to solve in prior art.
To achieve these goals, according to an aspect of the present invention, provide a kind of appraisal procedure of sub area division powering mode, comprising: the scoring calculating each index of sub area division powering mode, obtains standard evaluation matrix; Calculate the weights of each index, obtain weight matrix; And according to the comprehensive grading matrix of weight matrix and standard evaluation matrix computations sub area division powering mode.
Further, calculate the scoring of each index of sub area division powering mode, obtain standard evaluation matrix and comprise: according to the scoring of each index of predefine formulae discovery sub area division powering mode, obtain Evaluations matrix and standardization Evaluations matrix obtains standard evaluation matrix wherein, r represents the quantity of sub area division powering mode, and p represents the quantity of index, c ' ijrepresent the scoring of a jth index of i-th sub area division powering mode, c ijthe scoring of a jth index of i-th sub area division powering mode after representing standardization.
Further, standardization Evaluations matrix obtains standard evaluation matrix and comprises: according to formula or calculate each element in Evaluations matrix, obtain standard evaluation matrix, wherein i gets 1 to r successively, and j gets 1 to p successively, c ' ijrepresent the maximal value of the i-th row of Evaluations matrix, c ' ijrepresent the minimum value of Evaluations matrix i-th row.
Further, calculate the weights of each index, obtain weight matrix and comprise: the relative priority relation value receiving index, obtains precedence relation matrix wherein, b xyrepresent the preferred value of an xth index relative to y index, wherein, x and y gets 1 to p all successively; And the weights of each index are calculated according to precedence relation matrix, obtain weight matrix.
Further, calculate the weights of each index according to precedence relation matrix, obtain weight matrix and comprise: the eigenvalue of maximum λ calculating precedence relation matrix max; Calculate the eigenvalue of maximum λ of precedence relation matrix maxcharacteristic of correspondence vector T=[t 1t 2t q]; And according to formula calculate the weight w of each index i, obtain weight matrix W=[w 1w 2w q].
Further, comprise according to the comprehensive grading matrix of weight matrix and standard evaluation matrix computations sub area division powering mode: each element corresponding membership function set V=(v calculating standard evaluation Matrix C 1, v 2..., v m) in the relevant of each membership function be subordinate to angle value, obtain relevant Subject Matrix R=[R 1r 2r m], wherein, membership function set is predefined function set, for comprising the relevant matrix being subordinate to angle value of p element of the i-th row of standard evaluation Matrix C, r ijkrepresent c ija corresponding kth membership function v kbe relevantly subordinate to angle value, k gets 1 to m successively; And according to formula S i=WR ia calculates the comprehensive grading of sub area division powering mode, obtains comprehensive grading matrix S=[S 1s 2s m], wherein, S ibe the comprehensive grading of i-th sub area division powering mode, A=[a 1a 2a m], a kfor corresponding v kscore value.
Further, calculating the weights of each index according to precedence relation matrix, after obtaining weight matrix, the appraisal procedure of sub area division powering mode also comprises: the consistance of inspection precedence relation matrix; And if check out the consistance of precedence relation matrix nonconforming, again obtain precedence relation matrix.
Further, the consistance of precedence relation matrix is checked to comprise: according to formula calculate the consistance score value CI of precedence relation matrix; According to formula calculate the random consistance score value CR of precedence relation matrix, wherein the value of RI is determined by p; And the size of more random consistance score value CR and preset value, wherein, if random consistance score value CR is less than or equal to preset value, then determine that the consistance of precedence relation matrix meets the requirements.
To achieve these goals, according to a further aspect in the invention, a kind of apparatus for evaluating of sub area division powering mode is provided, the appraisal procedure of any one sub area division powering mode of this apparatus for evaluating for performing foregoing of the present invention and providing.
To achieve these goals, according to a further aspect in the invention, providing a kind of apparatus for evaluating of sub area division powering mode, comprising: the first computing unit, for calculating the scoring of each index of sub area division powering mode, obtaining standard evaluation matrix; Second computing unit, for calculating the weights of each index, obtains weight matrix; And the 3rd computing unit, for the comprehensive grading matrix according to weight matrix and standard evaluation matrix computations sub area division powering mode.
Further, the first computing unit comprises: the first computation subunit, for the scoring of each index according to predefine formulae discovery sub area division powering mode, obtains Evaluations matrix and normalizer unit, obtain standard evaluation matrix for standardization Evaluations matrix wherein, r represents the quantity of sub area division powering mode, and p represents the quantity of index, c ' ijrepresent the scoring of a jth index of i-th sub area division powering mode, c ijthe scoring of a jth index of i-th sub area division powering mode after representing standardization.
Further, normalizer unit comprises: the first computing module, for according to formula or calculate each element in Evaluations matrix, obtain standard evaluation matrix, wherein i gets 1 to r successively, and j gets 1 to p successively, c ' ijrepresent the maximal value of the i-th row of Evaluations matrix, c ' ijrepresent the minimum value of Evaluations matrix i-th row.
Further, the second computing unit comprises: receiving subelement, for receiving the relative priority relation value of index, obtaining precedence relation matrix wherein, b xyrepresent the preferred value of an xth index relative to y index, wherein x and y gets 1 to p all successively; And second computation subunit, for calculating the weights of each index according to precedence relation matrix, obtain weight matrix.
Further, the second computation subunit comprises: the second computing module, for calculating the eigenvalue of maximum λ of precedence relation matrix max; 3rd computing module, for calculating the eigenvalue of maximum λ of precedence relation matrix maxcharacteristic of correspondence vector T=[t 1t 2t q]; And the 4th computing module, for according to formula calculate the weight w of each index i, obtain weight matrix W=[w 1w 2w q].
Further, the 3rd computing unit comprises: the 3rd computation subunit, for calculating each element corresponding membership function set V=(v of standard evaluation Matrix C 1, v 2..., v m) in the relevant of each membership function be subordinate to angle value, obtain relevant Subject Matrix R=[R 1r 2r m], wherein, membership function set is predefined function set, for comprising the relevant matrix being subordinate to angle value of p element of the i-th row of standard evaluation Matrix C, r ijkrepresent c ija corresponding kth membership function v kbe relevantly subordinate to angle value, k gets 1 to m successively; And the 4th computation subunit, for according to formula S i=WR ia calculates the comprehensive grading of sub area division powering mode, obtains comprehensive grading matrix S=[S 1s 2s m], wherein, S ibe the comprehensive grading of i-th sub area division powering mode, A=[a 1a 2a m], a kfor corresponding v kscore value.
Further, the apparatus for evaluating of sub area division powering mode also comprises: verification unit, for checking the consistance of precedence relation matrix, wherein, if checking out the consistance of precedence relation matrix nonconforming, receiving subelement and again obtaining precedence relation matrix.
Further, verification unit comprises: the 5th computation subunit, for according to formula calculate the consistance score value CI of precedence relation matrix; 6th computation subunit, for according to formula calculate the random consistance score value CR of precedence relation matrix, wherein the value of RI is determined by p; And compare subelement, for the size of more random consistance score value CR and preset value, wherein, if random consistance score value CR is less than or equal to preset value, then determine that the consistance of precedence relation matrix meets the requirements.
The present invention adopts the scoring of each index calculating sub area division powering mode, obtains standard evaluation matrix; Calculate the weights of each index, obtain weight matrix; And according to the comprehensive grading matrix of weight matrix and standard evaluation matrix computations sub area division powering mode.By calculating scoring and the weights of each index, obtain weight matrix and standard evaluation matrix, recycling weight matrix and the comprehensive grading matrix of standard evaluation matrix to sub area division powering mode calculate, solving cannot the problem of comprehensive assessment sub area division power supply plan, and then reaches the effect of sub area division scheme of choose reasonable optimum.
Accompanying drawing explanation
The accompanying drawing forming a application's part is used to provide a further understanding of the present invention, and schematic description and description of the present invention, for explaining the present invention, does not form inappropriate limitation of the present invention.In the accompanying drawings:
Fig. 1 is the process flow diagram of the appraisal procedure of sub area division powering mode according to the embodiment of the present invention; And
Fig. 2 is the schematic diagram of the apparatus for evaluating of sub area division powering mode according to the embodiment of the present invention.
Embodiment
It should be noted that, when not conflicting, the embodiment in the application and the feature in embodiment can combine mutually.Below with reference to the accompanying drawings and describe the present invention in detail in conjunction with the embodiments.
The embodiment of the present invention provides a kind of appraisal procedure of sub area division powering mode, is specifically introduced below to the appraisal procedure of the sub area division powering mode of the embodiment of the present invention:
Fig. 1 is the process flow diagram of the appraisal procedure of the sub area division powering mode of the embodiment of the present invention.As shown in Figure 1, the method comprises following step S102 to step S106:
Step S102, calculates the scoring of each index of sub area division powering mode, obtains standard evaluation matrix.Particularly, before calculating, first formulate every evaluation index, evaluation index can comprise line load level, node voltage level, main-transformer Rate of average load, overall short circuit current level, the out-of-limit level of short-circuit current and apparatus overload level.Formulate a computing formula for each evaluation index, such as, following index calculate formula can be formulated:
(1) the horizontal index λ of line load l:
wherein, in formula, α is transmission line of electricity set of fingers; S lfor the power of branch road l; for S lthe upper limit.
(2) the horizontal λ of node voltage u:
λ U = Σ β ( U i - U i SP Δ U i lim ) 2 , Wherein β is bus set U i sp = ( U i H + U i L ) / 2 , Δ U i lim = ( U i H - U i L ) / 2 , U ibe i-th busbar voltage amplitude, with be respectively its upper and lower bound.
(3) main-transformer Rate of average load λ s:
wherein, z is the main-transformer number that in system, load factor exceedes a certain value, L ifor exceeding the load factor of each main-transformer of this threshold value.
(4) overall short circuit current level index λ sC:
wherein, in formula, γ is the set of bus isolating switch; I sCfor the short-circuit current value after subregion, I 0for the specified breaking current of isolating switch.
(5) the out-of-limit horizontal λ of short-circuit current z:
wherein, D ibe the short-circuit current level of i-th circuit, d is the number of lines.
(6) the horizontal λ of apparatus overload p:
wherein, e is the main-transformer number transshipped in system, l ifor the load factor of main-transformer.
In different sub area division powering mode, independent variable above in (1) to (6) item is different, thus the desired value obtained is also different, such as there is the index that p different, p computing formula is formulated for p index, have the powering mode that r different, each powering mode all has p desired value, is designated as by the Evaluations matrix calculated standardization is carried out to C ' and obtains standard evaluation matrix in standardisation process, be that profit evaluation model or cost type have the different computing formula of employing to calculate each element in Evaluations matrix according to index, obtain standard evaluation matrix, the computing formula of profit evaluation model is c ij = c ij ′ - min i c ij ′ max i c ij ′ - min i c ij ′ , The computing formula of cost type is c ij = c ij ′ - min i c ij ′ max i c ij ′ - min i c ij ′ . Profit evaluation model index refers to that this refers to that target value other scheme relative is the bigger the better, and cost type index is that this refers to that target value is the smaller the better relative to other scheme.First agriculture products is profit evaluation model or cost type, then selects corresponding formulae discovery, obtains standard evaluation matrix.Wherein the span of i is the integer from 1 to r, that is, i gets 1 to r successively, and the span of j is the integer from 1 to p, that is, j gets 1 to p successively, c ' ijrepresent the maximal value of the i-th row of Evaluations matrix, c ' ijrepresent the minimum value of Evaluations matrix i-th row.
Step S104, calculates the weights of each index, obtains weight matrix.Particularly, first by assessment experts, score value is provided to the relative priority relation between each index, and relative priority relation value is received, obtain precedence relation matrix wherein, p is the quantity of index, b xyrepresent the preferred value of an xth index relative to y index, the span of x and y is the integer of 1 to p, that is, x and y gets 1 to p all successively.And then calculate the eigenvalue of maximum λ of precedence relation matrix max, and calculate eigenvalue of maximum λ maxcharacteristic of correspondence vector T=[t 1t 2t q], according to formula calculate the weight w of index described in each i, namely obtain weight matrix W=[w 1w 2w q].
Particularly, electrical network assessment experts provides preferred value to each index, and can adopt 1 ~ 9 value method of Saaty, when adopting in this way, precedence relation matrix B has following characteristics: b xy>0, b xx=1, b xy=1/b yx.Such matrix is called just reciprocal type matrix.
Step S106, according to the comprehensive grading matrix of weight matrix and standard evaluation matrix computations sub area division powering mode.
Particularly, first, each element corresponding membership function set V=(v of standard evaluation Matrix C is calculated 1, v 2..., v m) in the relevant of each membership function be subordinate to angle value, obtain relevant Subject Matrix R=[R 1r 2r m], wherein, membership function set is predefined function set, for comprising the relevant matrix being subordinate to angle value of an i-th row p element of standard evaluation Matrix C, R ifor i-th row of Evaluations matrix C each element be subordinate to angle value relative to the relevant of each membership function, that is, R ibe the relevant subordinated-degree matrix of i-th sub area division powering mode, wherein, r ijkrepresent c ijto a kth membership function v kbe relevantly subordinate to angle value, the span of k is the integer of 1 to m, and k gets 1 to m successively.
Particularly, illustrate further for 5 membership functions below:
v 2 = c ij + 0.05 0.8 , 0 ≤ c ij ≤ 0.75 1.55 - c ij 0.8 , 0.75 ≤ c ij ≤ 1 ,
v 3 = c ij + 0 . 3 0.8 , 0 ≤ c ij ≤ 0.5 1.3 - c ij 0.8 , 0.5 ≤ c ij ≤ 1 ,
v 4 = c ij + 0.55 0.8 , 0 ≤ c ij ≤ 0.25 1.05 - c ij 0.8 , 0.25 ≤ c ij ≤ 1 ,
Corresponding V=(v 1, v 2, v 3, v 4, v 5) formulate score matrix be A, A=(9,7,5,3,1),
Then, according to formula S i=WR ia calculates the comprehensive grading of each sub area division powering mode, obtains comprehensive grading matrix S=[S 1s 2s m], wherein, S ibe the comprehensive grading of i-th sub area division powering mode, A=[a 1a 2a m], a kfor corresponding v kscore value.For stating conveniently, the comprehensive grading of a sub area division powering mode once can be calculated.
The appraisal procedure of the sub area division powering mode of the embodiment of the present invention, by calculating scoring and the weights of each index, obtain weight matrix and standard evaluation matrix, recycling weight matrix and the comprehensive grading matrix of standard evaluation matrix to sub area division powering mode calculate, solving cannot the problem of comprehensive assessment sub area division powering mode, and then reaches the effect of sub area division powering mode of choose reasonable optimum.
Further, the weights of each index are being calculated according to precedence relation matrix, after obtaining weight matrix, the appraisal procedure of the sub area division powering mode that the embodiment of the present invention provides also comprises: the consistance of inspection precedence relation matrix, if check out the consistance of precedence relation matrix nonconforming, again obtain precedence relation matrix.Rationally effective to ensure the matrix formed.
Particularly, first according to formula calculate the consistance score value CI of precedence relation matrix, then according to formula calculate the random consistance score value CR of precedence relation matrix, wherein, the value of RI is determined by p.Be the precedence relation matrix of 1 to 10 for index quantity p, RI value is as shown in table 1 below:
Table 1
P 1 2 3 4 5 6 7 8 9 10
RI 0.00 0.00 0.58 0.90 1.12 1.24 1.36 1.41 1.45 1.49
More than the situation of 10, RI value is also have p well-determined, can be obtained by the mode of tabling look-up.Therefore obvious index quantity p of the present invention is not defined as within 10.
Wherein, after calculating consistance score value CI, if CI is 0, then can determine that the consistance of precedence relation matrix meets the requirements, CI value is larger, illustrates that the consistance of precedence relation matrix is poorer.Further introducing random consistance score value CR, the relatively size of CR and preset value, wherein, if random consistance score value CR is less than or equal to preset value, then determine that the consistance of precedence relation matrix meets the requirements, in embodiments of the present invention, preset value can be set to 0.1, namely, if CR≤0.1, determine that the consistance of precedence relation matrix meets the requirements, otherwise, again obtain precedence relation matrix.
The embodiment of the present invention additionally provides a kind of apparatus for evaluating of sub area division powering mode, this apparatus for evaluating is mainly used in the appraisal procedure performing any one sub area division powering mode that embodiment of the present invention foregoing provides, and is specifically introduced below to the apparatus for evaluating of the sub area division powering mode that the embodiment of the present invention provides:
Fig. 2 is the schematic diagram of the apparatus for evaluating of sub area division powering mode according to the embodiment of the present invention, and as shown in Figure 2, this apparatus for evaluating mainly comprises the first computing unit 10, second computing unit 20 and the 3rd computing unit 30, wherein:
First computing unit 10 is for calculating the scoring of each index of sub area division powering mode, obtain standard evaluation matrix, particularly, before calculating, first formulate every evaluation index, evaluation index can comprise line load level, node voltage level, main-transformer Rate of average load, overall short circuit current level, the out-of-limit level of short-circuit current and apparatus overload level.Formulate a computing formula for each evaluation index, such as, following index calculate formula can be formulated:
(1) the horizontal index λ of line load l:
wherein, in formula, α is transmission line of electricity set of fingers; S lfor the power of branch road l; for S lthe upper limit.
(2) the horizontal λ of node voltage u:
λ U = Σ β ( U i - U i SP Δ U i lim ) 2 , Wherein β is bus set U i sp = ( U i H + U i L ) / 2 , Δ U i lim = ( U i H - U i L ) / 2 , U ibe i-th busbar voltage amplitude, with be respectively its upper and lower bound.
(3) main-transformer Rate of average load λ s:
wherein, z is the main-transformer number that in system, load factor exceedes a certain value, L ifor exceeding the load factor of each main-transformer of this threshold value.
(4) overall short circuit current level index λ sC:
wherein, in formula, γ is the set of bus isolating switch; I sCfor the short-circuit current value after subregion, I 0for the specified breaking current of isolating switch.
(5) the out-of-limit horizontal λ of short-circuit current z:
wherein, D ibe the short-circuit current level of i-th circuit, d is the number of lines.
(6) the horizontal λ of apparatus overload p:
wherein, e is the main-transformer number transshipped in system, l ifor the load factor of main-transformer.
Wherein, first computing unit 10 mainly comprises the first computation subunit and normalizer unit, in different sub area division powering mode, independent variable above in (1) to (6) item is different, thus the desired value obtained is also different, such as there is the index that p different, p computing formula is formulated for p index, there is the powering mode that r different, each powering mode all has p desired value, first computation subunit is used for the scoring according to above-mentioned each index of predefine formulae discovery, and the Evaluations matrix obtained is designated as normalizer unit is used for carrying out standardization to C ' and obtains standard evaluation matrix particularly, normalizer unit comprises the first computing module, it is that profit evaluation model or cost type have the different computing formula of employing to calculate each element in Evaluations matrix that this first computing module is used for according to index, and obtain standard evaluation matrix, the computing formula of profit evaluation model is c ij = c ij ′ - min i c ij ′ max i c ij ′ - min i c ij ′ , The computing formula of cost type is c ij = c ij ′ - min i c ij ′ max i c ij ′ - min i c ij ′ . Profit evaluation model index refers to that this refers to that target value other scheme relative is the bigger the better, and cost type index is that this refers to that target value is the smaller the better relative to other scheme.First agriculture products is profit evaluation model or cost type, then selects corresponding formulae discovery, obtains standard evaluation matrix.Wherein the span of i is the integer from 1 to r, that is, i gets 1 to r successively, and the span of j is the integer from 1 to p, that is, j gets 1 to p successively, c ' ijrepresent the maximal value of the i-th row of Evaluations matrix, c ' ijrepresent the minimum value of Evaluations matrix i-th row.
Second computing unit 20 is for calculating the weights of each index, obtain weight matrix, particularly, second computing unit 20 mainly comprises reception subelement and the second computation subunit, for the calculating of each index weight value, first by assessment experts, score value is provided to the relative priority relation between each index, and by reception subelement, relative priority relation value is received, obtain precedence relation matrix wherein, p is the quantity of index, b xyrepresent the preferred value of an xth index relative to y index, the span of x and y is the integer of 1 to p, that is, x and y gets 1 to p all successively.Then calculated the weights of each index by the second computation subunit according to precedence relation matrix, obtain weight matrix.
Wherein, the second computation subunit mainly comprises the second computing module, the 3rd computing module and the 4th computing module, and wherein, the second computing module is for calculating the eigenvalue of maximum λ of precedence relation matrix max, the 3rd computing module is for calculating eigenvalue of maximum λ maxcharacteristic of correspondence vector T=[t 1t 2t q], the 4th computing module is used for according to formula calculate the weight w of index described in each i, namely obtain weight matrix W=[w 1w 2w q].
Particularly, electrical network assessment experts provides preferred value to each index, and can adopt 1 ~ 9 value method of Saaty, when adopting in this way, precedence relation matrix B has following characteristics: b xy>0, b xx=1, b xy=1/b yx.Such matrix is called just reciprocal type matrix.
3rd computing unit 30 is for the comprehensive grading matrix according to weight matrix and standard evaluation matrix computations sub area division powering mode, and wherein, the 3rd computing unit 30 mainly comprises the 3rd computation subunit and the 4th computation subunit.
3rd computation subunit is for calculating each element corresponding membership function set V=(v of standard evaluation Matrix C 1, v 2..., v m) in the relevant of each membership function be subordinate to angle value, obtain relevant Subject Matrix R=[R 1r 2r m], wherein, membership function set is predefined function set, for comprising the relevant matrix being subordinate to angle value of an i-th row p element of standard evaluation Matrix C, R ifor i-th row of Evaluations matrix C each element be subordinate to angle value relative to the relevant of each membership function, that is, R ibe the relevant subordinated-degree matrix of i-th sub area division powering mode, wherein, r ijkrepresent c ijto a kth membership function v kbe relevantly subordinate to angle value, the span of k is the integer of 1 to m, and k gets 1 to m successively.
Particularly, illustrate further for 5 membership functions below:
v 2 = c ij + 0.05 0.8 , 0 ≤ c ij ≤ 0.75 1.55 - c ij 0.8 , 0.75 ≤ c ij ≤ 1 ,
v 3 = c ij + 0 . 3 0.8 , 0 ≤ c ij ≤ 0.5 1.3 - c ij 0.8 , 0.5 ≤ c ij ≤ 1 ,
v 4 = c ij + 0.55 0.8 , 0 ≤ c ij ≤ 0.25 1.05 - c ij 0.8 , 0.25 ≤ c ij ≤ 1 ,
Corresponding V=(v 1, v 2, v 3, v 4, v 5) formulate score matrix be A, A=(9,7,5,3,1),
4th computation subunit is used for according to formula S i=WR ia calculates the comprehensive grading of each sub area division powering mode, obtains comprehensive grading matrix S=[S 1s 2s m], wherein, S ibe the comprehensive grading of i-th sub area division powering mode, A=[a 1a 2a m], a kfor corresponding v kscore value.For stating conveniently, the comprehensive grading of a sub area division powering mode once can be calculated.
The apparatus for evaluating of the sub area division powering mode of the embodiment of the present invention, by calculating scoring and the weights of each index, obtain weight matrix and standard evaluation matrix, recycling weight matrix and the comprehensive grading matrix of standard evaluation matrix to sub area division powering mode calculate, solving cannot the problem of comprehensive assessment sub area division powering mode, and then reaches the effect of sub area division powering mode of choose reasonable optimum.
Further, the apparatus for evaluating of the sub area division powering mode that the embodiment of the present invention provides also comprises verification unit, the weights of each index are being calculated according to precedence relation matrix, after obtaining weight matrix, verification unit is for checking the consistance of precedence relation matrix, if check out the consistance of precedence relation matrix nonconforming, then again obtain precedence relation matrix by reception subelement.Rationally effective to ensure the matrix formed.
Particularly, verification unit mainly comprises the 5th computation subunit, the 6th computation subunit and compares subelement, wherein:
5th computation subunit is used for according to formula calculate the consistance score value CI of precedence relation matrix.
6th computation subunit is used for according to formula calculate the random consistance score value CR of precedence relation matrix, wherein the value of RI is determined by p.Be the precedence relation matrix of 1 to 10 for index quantity p, RI value as shown in Table 1.
Relatively subelement is used for the size of more random consistance score value CR and preset value, wherein, if random consistance score value CR is less than or equal to preset value, then determines that the consistance of precedence relation matrix meets the requirements.
Wherein, after calculating consistance score value CI, if CI is 0, then can determine that the consistance of precedence relation matrix meets the requirements, CI value is larger, illustrates that the consistance of precedence relation matrix is poorer.Further introducing random consistance score value CR, the relatively size of CR and preset value, wherein, if random consistance score value CR is less than or equal to preset value, then determine that the consistance of precedence relation matrix meets the requirements, in embodiments of the present invention, preset value can be set to 0.1, namely, if CR≤0.1, determine that the consistance of precedence relation matrix meets the requirements, otherwise, again obtain precedence relation matrix.
As can be seen from the above description, present invention achieves comprehensive assessment sub area division powering mode, reach the effect of the sub area division powering mode of choose reasonable optimum.
It should be noted that, can perform in the computer system of such as one group of computer executable instructions in the step shown in the process flow diagram of accompanying drawing, and, although show logical order in flow charts, but in some cases, can be different from the step shown or described by order execution herein.
Obviously, those skilled in the art should be understood that, above-mentioned of the present invention each module or each step can realize with general calculation element, they can concentrate on single calculation element, or be distributed on network that multiple calculation element forms, alternatively, they can realize with the executable program code of calculation element, thus, they can be stored and be performed by calculation element in the storage device, or they are made into each integrated circuit modules respectively, or the multiple module in them or step are made into single integrated circuit module to realize.Like this, the present invention is not restricted to any specific hardware and software combination.
The foregoing is only the preferred embodiments of the present invention, be not limited to the present invention, for a person skilled in the art, the present invention can have various modifications and variations.Within the spirit and principles in the present invention all, any amendment done, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (16)

1. an appraisal procedure for sub area division powering mode, is characterized in that, comprising:
Calculate the scoring of each index of described sub area division powering mode, obtain standard evaluation matrix;
Calculate the weights of each described index, obtain weight matrix; And
The comprehensive grading matrix of sub area division powering mode according to described weight matrix and described standard evaluation matrix computations.
2. the appraisal procedure of sub area division powering mode according to claim 1, is characterized in that, calculates the scoring of each index of described sub area division powering mode, obtains standard evaluation matrix and comprises:
According to the scoring of each described index of sub area division powering mode described in predefine formulae discovery, obtain Evaluations matrix and
Evaluations matrix described in standardization obtains described standard evaluation matrix wherein, r represents the quantity of described sub area division powering mode, and p represents the quantity of described index, c ' ijrepresent the scoring of a jth index of i-th described sub area division powering mode, c ijthe scoring of a jth index of i-th described sub area division powering mode after representing standardization.
3. the appraisal procedure of sub area division powering mode according to claim 2, is characterized in that, Evaluations matrix described in standardization obtains described standard evaluation matrix and comprises:
According to formula c ij = c ij ′ - min i c ij ′ max i c ij ′ - min i c ij ′ Or c ij = max i c ij ′ - c ij ′ max i c ij ′ - min i c ij ′ Calculate each element in described Evaluations matrix, obtain described standard evaluation matrix, wherein i gets 1 to r successively, and j gets 1 to p successively, c ' ijrepresent the maximal value of the i-th row of described Evaluations matrix, c ' ijrepresent the minimum value of described Evaluations matrix i-th row.
4. the appraisal procedure of sub area division powering mode according to claim 1, is characterized in that, calculates the weights of each described index, obtains weight matrix and comprises:
Receive the relative priority relation value of described index, obtain precedence relation matrix wherein, b xyrepresent the preferred value of an xth described index relative to the individual described index of y, wherein, x and y gets 1 to p all successively; And
Calculate the weights of each described index according to described precedence relation matrix, obtain described weight matrix.
5. the appraisal procedure of sub area division powering mode according to claim 4, is characterized in that, calculates the weights of each described index according to described precedence relation matrix, obtains described weight matrix and comprises:
Calculate the eigenvalue of maximum λ of described precedence relation matrix max;
Calculate the described eigenvalue of maximum λ of described precedence relation matrix maxcharacteristic of correspondence vector T=[t 1t 2t q]; And
According to formula calculate the weight w of index described in each i, obtain described weight matrix W=[w 1w 2w q].
6. the appraisal procedure of sub area division powering mode according to claim 2, is characterized in that, according to described weight matrix and described standard evaluation matrix computations, the comprehensive grading matrix of sub area division powering mode comprises:
Calculate each element corresponding membership function set V=(v of described standard evaluation Matrix C 1, v 2..., v m) in the relevant of each membership function be subordinate to angle value, obtain relevant Subject Matrix R=[R 1r 2r m], wherein, described membership function set is predefined function set, for comprising the relevant matrix being subordinate to angle value of p element of the i-th row of described standard evaluation Matrix C, r ijkrepresent c ija corresponding kth membership function v kbe relevantly subordinate to angle value, k gets 1 to m successively; And
According to formula S i=WR ia calculates the comprehensive grading of described sub area division powering mode, obtains comprehensive grading matrix S=[S 1s 2s m], wherein, S ibe the comprehensive grading of i-th described sub area division powering mode, A=[a 1a 2a m], a kfor corresponding v kscore value.
7. the appraisal procedure of sub area division powering mode according to claim 5, it is characterized in that, calculating the weights of each described index according to described precedence relation matrix, after obtaining described weight matrix, the appraisal procedure of described sub area division powering mode also comprises:
Check the consistance of described precedence relation matrix; And
If check out the consistance of described precedence relation matrix nonconforming, again obtain described precedence relation matrix.
8. the appraisal procedure of sub area division powering mode according to claim 7, is characterized in that, checks the consistance of described precedence relation matrix to comprise:
According to formula calculate the consistance score value CI of described precedence relation matrix;
According to formula calculate the random consistance score value CR of described precedence relation matrix, wherein the value of RI is determined by p; And
The size of more described random consistance score value CR and preset value, wherein, if described random consistance score value CR is less than or equal to described preset value, then determines that the consistance of described precedence relation matrix meets the requirements.
9. an apparatus for evaluating for sub area division powering mode, is characterized in that, comprising:
First computing unit, for calculating the scoring of each index of described sub area division powering mode, obtains standard evaluation matrix;
Second computing unit, for calculating the weights of each described index, obtains weight matrix; And
3rd computing unit, for the comprehensive grading matrix of sub area division powering mode according to described weight matrix and described standard evaluation matrix computations.
10. the apparatus for evaluating of sub area division powering mode according to claim 9, is characterized in that, described first computing unit comprises:
First computation subunit, for the scoring of each the described index according to sub area division powering mode described in predefine formulae discovery, obtains Evaluations matrix and
Normalizer unit, obtains described standard evaluation matrix for Evaluations matrix described in standardization wherein, r represents the quantity of described sub area division powering mode, and p represents the quantity of described index, c ' ijrepresent the scoring of a jth index of i-th described sub area division powering mode, c ijthe scoring of a jth index of i-th described sub area division powering mode after representing standardization.
The apparatus for evaluating of 11. sub area division powering modes according to claim 10, is characterized in that, described normalizer unit comprises:
First computing module, for according to formula c ij = c ij ′ - min i c ij ′ max i c ij ′ - min i c ij ′ Or c ij = max i c ij ′ - c ij ′ max i c ij ′ - min i c ij ′ Calculate each element in described Evaluations matrix, obtain described standard evaluation matrix, wherein i gets 1 to r successively, and j gets 1 to p successively, c ' ijrepresent the maximal value of the i-th row of described Evaluations matrix, c ' ijrepresent the minimum value of described Evaluations matrix i-th row.
The apparatus for evaluating of 12. sub area division powering modes according to claim 9, is characterized in that, described second computing unit comprises:
Receiving subelement, for receiving the relative priority relation value of described index, obtaining precedence relation matrix wherein, b xyrepresent the preferred value of an xth described index relative to the individual described index of y, wherein x and y gets 1 to p all successively; And
Second computation subunit, for calculating the weights of each described index according to described precedence relation matrix, obtains described weight matrix.
The apparatus for evaluating of 13. sub area division powering modes according to claim 12, is characterized in that, described second computation subunit comprises:
Second computing module, for calculating the eigenvalue of maximum λ of described precedence relation matrix max;
3rd computing module, for calculating the described eigenvalue of maximum λ of described precedence relation matrix maxcharacteristic of correspondence vector T=[t 1t 2t q]; And
4th computing module, for according to formula calculate the weight w of index described in each i, obtain described weight matrix W=[w 1w 2w q].
The apparatus for evaluating of 14. sub area division powering modes according to claim 10, is characterized in that, described 3rd computing unit comprises:
3rd computation subunit, for calculating each element corresponding membership function set V=(v of described standard evaluation Matrix C 1, v 2..., v m) in the relevant of each membership function be subordinate to angle value, obtain relevant Subject Matrix R=[R 1r 2r m], wherein, described membership function set is predefined function set, for comprising the relevant matrix being subordinate to angle value of p element of the i-th row of described standard evaluation Matrix C, r ijkrepresent c ija corresponding kth membership function v kbe relevantly subordinate to angle value, k gets 1 to m successively; And
4th computation subunit, for according to formula S i=WR ia calculates the comprehensive grading of described sub area division powering mode, obtains comprehensive grading matrix S=[S 1s 2s m], wherein, S ibe the comprehensive grading of i-th described sub area division powering mode, A=[a 1a 2a m], a kfor corresponding v kscore value.
The apparatus for evaluating of 15. sub area division powering modes according to claim 13, is characterized in that, the apparatus for evaluating of described sub area division powering mode also comprises:
Verification unit, for checking the consistance of described precedence relation matrix,
Wherein, if check out the consistance of described precedence relation matrix nonconforming, described reception subelement obtains described precedence relation matrix again.
The apparatus for evaluating of 16. sub area division powering modes according to claim 15, it is characterized in that, described verification unit comprises:
5th computation subunit, for according to formula calculate the consistance score value CI of described precedence relation matrix;
6th computation subunit, for according to formula calculate the random consistance score value CR of described precedence relation matrix, wherein the value of RI is determined by p; And
Relatively subelement, for the size of more described random consistance score value CR and preset value, wherein, if described random consistance score value CR is less than or equal to described preset value, then determines that the consistance of described precedence relation matrix meets the requirements.
CN201310654861.7A 2013-12-05 2013-12-05 Assessment method and device of grid partition power supply mode Pending CN104700317A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310654861.7A CN104700317A (en) 2013-12-05 2013-12-05 Assessment method and device of grid partition power supply mode

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310654861.7A CN104700317A (en) 2013-12-05 2013-12-05 Assessment method and device of grid partition power supply mode

Publications (1)

Publication Number Publication Date
CN104700317A true CN104700317A (en) 2015-06-10

Family

ID=53347407

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310654861.7A Pending CN104700317A (en) 2013-12-05 2013-12-05 Assessment method and device of grid partition power supply mode

Country Status (1)

Country Link
CN (1) CN104700317A (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040167731A1 (en) * 2002-12-23 2004-08-26 Abb, Inc. Failure rate adjustment for electric power network reliability analysis
CN101930491A (en) * 2010-08-18 2010-12-29 华北电力大学 Method for comprehensively evaluating measures for limiting short-circuit current

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040167731A1 (en) * 2002-12-23 2004-08-26 Abb, Inc. Failure rate adjustment for electric power network reliability analysis
CN101930491A (en) * 2010-08-18 2010-12-29 华北电力大学 Method for comprehensively evaluating measures for limiting short-circuit current

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
于会泉: "电网短路电流限制措施的优化配置研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 *
吴耀文等: "大规模风电特高压专用通道落点优选方法", 《中国电机工程学报》 *

Similar Documents

Publication Publication Date Title
CN100541965C (en) Transmission Expansion Planning in Electric method based on the certainty bilevel programming model
CN103258299B (en) A kind of many direct currents concentrate the receiving end Net Frame of Electric Network optimization method of feed-in
AlHajri et al. Exploiting the radial distribution structure in developing a fast and flexible radial power flow for unbalanced three-phase networks
CN103311926B (en) Power system cascading failure simulation method based on unified power flow controller
CN104992266A (en) Method of determining power grid node importance degree and system thereof
CN102623987B (en) Multiple-DC (direct current)-droppoint selection method based on multiple feed-in short circuit ratios
CN104598725A (en) Evaluation method for power distribution network adaptability of distributed power supply after grid connection
CN106600139A (en) Reliability evaluation method of power distribution network
CN104901308B (en) Power system critical circuits discrimination method
CN105870935B (en) Radial distribution networks idle work optimization method based on clustering algorithm
CN106960394A (en) A kind of AC-DC hybrid power grid ability to transmit electricity appraisal procedure based on Monte Carlo
CN106532717A (en) Comprehensive sensitivity analysis-based circuit overload load-shedding coordination optimization method
CN105303460A (en) Identification method of key nodes and key branches in power grid
CN102609792A (en) Extra-high-voltage alternating current and direct current transmission mode applicability selecting method and device using same
CN106156885A (en) A kind of electrical network weaker overall integrated evaluating method
CN105389670A (en) Power grid node importance degree determination method
CN109615245A (en) A kind of method of pair of Electric Power Network Planning total evaluation
CN104376506A (en) Electric power system risk assessment method based on deterministic analysis
CN107679289A (en) A kind of dynamic passive compensation collocation method for reducing multi-infeed HVDC commutation failure risk
CN106548413A (en) A kind of power system energy storage fitness-for-service assessment method and system
CN107622360A (en) A kind of critical circuits recognition methods for considering subjective and objective factor
CN106655201A (en) Security domain-based safe optimization and control method for electric power thermal stability
CN107425527A (en) A kind of THE UPFC static security prevention and control method
CN107749618A (en) A kind of short circuit current current limliting optimization method in power network
CN104156879A (en) Electric power system protection and control method based on risk assessment

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20150610