CN103529275B - Area power grid Analysis of Short-Circuit Current method and apparatus - Google Patents
Area power grid Analysis of Short-Circuit Current method and apparatus Download PDFInfo
- Publication number
- CN103529275B CN103529275B CN201310511797.7A CN201310511797A CN103529275B CN 103529275 B CN103529275 B CN 103529275B CN 201310511797 A CN201310511797 A CN 201310511797A CN 103529275 B CN103529275 B CN 103529275B
- Authority
- CN
- China
- Prior art keywords
- region
- matrix
- borderline region
- circuit current
- branch road
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 54
- 238000004458 analytical method Methods 0.000 title claims abstract description 43
- 238000004364 calculation method Methods 0.000 claims abstract description 45
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims description 114
- 230000008030 elimination Effects 0.000 claims description 10
- 238000003379 elimination reaction Methods 0.000 claims description 10
- 230000002093 peripheral effect Effects 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 5
- 238000010586 diagram Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 8
- 230000008859 change Effects 0.000 description 5
- 239000000203 mixture Substances 0.000 description 3
- 101100390985 Schizosaccharomyces pombe (strain 972 / ATCC 24843) fmc1 gene Proteins 0.000 description 2
- 101100460197 Schizosaccharomyces pombe (strain 972 / ATCC 24843) new13 gene Proteins 0.000 description 2
- 101100460198 Schizosaccharomyces pombe (strain 972 / ATCC 24843) new14 gene Proteins 0.000 description 2
- 101100459356 Schizosaccharomyces pombe (strain 972 / ATCC 24843) new18 gene Proteins 0.000 description 2
- 101100460203 Schizosaccharomyces pombe (strain 972 / ATCC 24843) new2 gene Proteins 0.000 description 2
- 101100460202 Schizosaccharomyces pombe (strain 972 / ATCC 24843) new25 gene Proteins 0.000 description 2
- 101100460204 Schizosaccharomyces pombe (strain 972 / ATCC 24843) new4 gene Proteins 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 101100167780 Schizosaccharomyces pombe (strain 972 / ATCC 24843) coa2 gene Proteins 0.000 description 1
- 101100172625 Schizosaccharomyces pombe (strain 972 / ATCC 24843) erh1 gene Proteins 0.000 description 1
- 101100228799 Schizosaccharomyces pombe (strain 972 / ATCC 24843) gem6 gene Proteins 0.000 description 1
- 101100322386 Schizosaccharomyces pombe (strain 972 / ATCC 24843) new8 gene Proteins 0.000 description 1
- 101100460205 Schizosaccharomyces pombe (strain 972 / ATCC 24843) new9 gene Proteins 0.000 description 1
- 101100255212 Schizosaccharomyces pombe (strain 972 / ATCC 24843) rsa3 gene Proteins 0.000 description 1
- 101100210181 Schizosaccharomyces pombe (strain 972 / ATCC 24843) vts1 gene Proteins 0.000 description 1
- 101100273916 Schizosaccharomyces pombe (strain 972 / ATCC 24843) wip1 gene Proteins 0.000 description 1
- 241000396836 Trinodes Species 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 239000004020 conductor Substances 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000002787 reinforcement Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Abstract
The invention provides a kind of area power grid Analysis of Short-Circuit Current method and apparatus, comprise: Region dividing is carried out to network system: regional classification Water demand being calculated short circuit current level is survey region, being borderline region by there is the regional classification of directly getting in touch with survey region, is perimeter by the regional classification except survey region and borderline region; Borderline region and perimeter are carried out equivalent simplification, forms the borderline region after equivalence, obtain the network system after the equivalence that the borderline region after by survey region and equivalence forms; According to the network system after equivalence, the short circuit current level of analytical calculation survey region.The invention solves in prior art and be difficult to while guarantee calculation of short-circuit current precision, the technical matters of data volume during minimizing Analysis of Short-Circuit Current, reach on the basis ensureing calculation of short-circuit current precision, data volume during minimizing Analysis of Short-Circuit Current and the technique effect of calculated amount.
Description
Technical field
The present invention relates to technical field of power systems, particularly a kind of area power grid Analysis of Short-Circuit Current method and apparatus.
Background technology
Calculation of short-circuit current; be a requisite basic calculating in the design and running of electric system, the result of calculation of short-circuit current is generally applied to the determination etc. of the selection of electrical equipment and current－carrying conductor, the determination of protective relaying device, Main Wire Connection Scheme and the method for operation.In addition, carry out electric system transient stability analysis, determine that transmission line of electricity also all must carry out calculation of short-circuit current in interference of communication etc.
In recent years along with the continuous expansion of electric system scale; contact between electrical network strengthens gradually; the development of extensive interconnected network can realize the reasonable development of resource and distribute, improves the reliability of power supply and the economy of operation rationally on the one hand; on the other hand because the reinforcement of electric network composition makes electrical distance shorten; the short-circuit current of part plant stand increases more; even exceeded the rupturing capacity of switch, short-circuit current problem has more and more caused the attention of power planning, management and running department.
At present, time calculation of short-circuit current is carried out to certain electrical network, often only pay close attention to the short circuit current level of this electrical network location (survey region can be called), but, because mutually can have an impact and disturb between electrical network, how to ensureing the accuracy of survey region being carried out to calculation of short-circuit current, with regard to needing, computation and analysis being carried out to the region in electric system except survey region.If all carry out analytical calculation to the data of whole electric power networks, so computing velocity will be very slow, and also need to safeguard the electric network data in other area except survey region, this will greatly increase data volume and the calculated amount of Analysis of Short-Circuit Current simultaneously.
Simultaneously, along with the continuous expansion of electrical network scale, in order to ensure the accuracy of calculation of short-circuit current, the data volume of operation of power networks personnel when carrying out Analysis of Short-Circuit Current and calculating becomes more and more huger, the data volume of carrying out data analysis also gets more and more, and so also just must result in and become more and more higher to the requirement of simulation software and allocation of computer.Meanwhile, owing to being subject to the restriction of the algorithm of program own, also there is certain requirement to the scale of electrical network in some business softwares.
How while guarantee calculation of short-circuit current precision, reduce data volume during Analysis of Short-Circuit Current, not yet propose effective solution at present.
Summary of the invention
The invention provides a kind of area power grid Analysis of Short-Circuit Current method and apparatus, be difficult to while guarantee calculation of short-circuit current precision at least to solve in prior art, the technical matters of data volume during minimizing Analysis of Short-Circuit Current.
Embodiments provide a kind of area power grid Analysis of Short-Circuit Current method, comprise: Region dividing is carried out to network system: regional classification Water demand being calculated short circuit current level is survey region, being borderline region by there is the regional classification of directly getting in touch with survey region, is perimeter by the regional classification except survey region and borderline region; Borderline region and perimeter are carried out equivalent simplification, forms the borderline region after equivalence, obtain the network system after the equivalence that the borderline region after by survey region and equivalence forms; According to the network system after equivalence, the short circuit current level of analytical calculation survey region.
In one embodiment, borderline region and perimeter are carried out equivalent simplification, form the borderline region after equivalence, obtain the network system after the equivalence that the borderline region after by survey region and equivalence forms, comprising: according to the bus admittance matrix of the sequential build network system of perimeter, borderline region and survey region; Utilize the part in Gaussian elimination method cancellation bus admittance matrix peripheral region, form the equivalent posterior nodal point admittance matrix only retaining borderline region and survey region; The variable quantity matrix of the matrix element that borderline region is corresponding in certainty equivalents posterior nodal point admittance matrix; According to the connection branch road of the generator increased in the borderline region after described variable quantity matrix certainty equivalents, increase, and the reactance parameter of the reactance parameter of generator increased and the connection branch road of increase; According to the generator of the increase determined, the connection branch road of increase, the reactance parameter of the generator of increase and the reactance parameter of connection branch road increased carry out network reconfiguration form equivalence after network system.
In one embodiment, according to the reactance parameter of the connection branch road increased in the borderline region after described variable quantity matrix certainty equivalents, comprising: the node number determining the connection branch road two ends increased; Get the negative inverse of element corresponding with the node number at two ends in variable quantity matrix, as the reactance parameter of the connection branch road increased.
In one embodiment, according to the reactance parameter of the generator increased in the borderline region after described variable quantity matrix certainty equivalents, comprising: the node number determining the node be connected with the generator increased; By this node number in the admittance of the connection branch road be directly connected with the node of this node number in perimeter, variable quantity matrix be expert on each element obtain resultant admittance as sum operation; Get the reactance parameter of inverse as the generator increased of resultant admittance.
In one embodiment, the variable quantity matrix of the matrix element that borderline region is corresponding in certainty equivalents posterior nodal point admittance matrix, comprising: the number determining borderline region interior joint; When the number of borderline region interior joint is less than predetermined threshold, the method adopting linear equation to solve calculates variable quantity matrix.
In one embodiment, described predetermined threshold is 6.
The embodiment of the present invention additionally provides a kind of area power grid Analysis of Short-Circuit Current device, comprise: Region dividing module, for carrying out Region dividing to network system: regional classification Water demand being calculated short circuit current level is survey region, being borderline region by there is the regional classification of directly getting in touch with survey region, is perimeter by the regional classification except survey region and borderline region; Equivalent Simplification module, for borderline region and perimeter are carried out equivalent simplification, forms the borderline region after equivalence, obtains the network system after the equivalence that the borderline region after by survey region and equivalence forms; Analysis module, for according to the network system after equivalence, the short circuit current level of analytical calculation survey region.
In one embodiment, described Equivalent Simplification module comprises: construction unit, for the bus admittance matrix of the sequential build network system according to perimeter, borderline region and survey region; Simplified element, for utilizing the part in Gaussian elimination method cancellation bus admittance matrix peripheral region, forms the equivalent posterior nodal point admittance matrix only retaining borderline region and survey region; Matrix determining unit, for the variable quantity matrix of matrix element corresponding to borderline region in certainty equivalents posterior nodal point admittance matrix; Parameter determination unit, for the connection branch road according to the generator increased in the borderline region after described variable quantity matrix certainty equivalents, increase, and the reactance parameter of the reactance parameter of generator increased and the connection branch road of increase; Reconfiguration unit, for the connection branch road according to the generator increased, increase, the reactance parameter of the generator of increase and the reactance parameter of connection branch road increased carry out network reconfiguration form equivalence after network system.
In one embodiment, described parameter determination unit comprises: first segment period determination subelement, for determining the node number at the connection branch road two ends increased; First parameter determination subelement, for getting the negative inverse of element corresponding with the node number at two ends in variable quantity matrix, as the reactance parameter of the connection branch road increased.
In one embodiment, described parameter determination unit comprises: second section period determination subelement, for the node number of the node that the generator determined with increase is connected; Summing elements, for by this node number in the admittance of the connection branch road be directly connected with the node of this node number in perimeter, variable quantity matrix be expert on each element obtain resultant admittance as sum operation; Second parameter determination subelement, for getting the reactance parameter of inverse as the generator increased of resultant admittance.
In one embodiment, described matrix determining unit comprises: nodes determination subelement, for determining the number of borderline region interior joint; Computation subunit, for when the number of borderline region interior joint is less than predetermined threshold, the method adopting linear equation to solve calculates variable quantity matrix.
In one embodiment, described predetermined threshold is 6.
In embodiments of the present invention, it is survey region that Water demand in network system is calculated the regional classification of short circuit current level, then be borderline region by there is the regional classification of directly getting in touch with survey region, be perimeter by the regional classification except survey region and borderline region, then be the network system of the borderline region after only comprising survey region and equivalence by this network system Equivalent Simplification, finally by the short circuit current level of the network system analytical calculation survey region after equivalence, because no longer there is perimeter in the network system after Equivalent Simplification, therefore required when analytical calculation short circuit current level data volume will obviously reduce.Efficiently solve by the way in prior art and be difficult to while guarantee calculation of short-circuit current precision, the technical matters of data volume during minimizing Analysis of Short-Circuit Current, reach on the basis ensureing calculation of short-circuit current precision, data volume during minimizing Analysis of Short-Circuit Current and the technique effect of calculated amount.
Accompanying drawing explanation
Accompanying drawing described herein is used to provide a further understanding of the present invention, forms a application's part, does not form limitation of the invention.In the accompanying drawings:
Fig. 1 is the process flow diagram of the area power grid Analysis of Short-Circuit Current method of the embodiment of the present invention;
Fig. 2 is the process flow diagram carrying out Equivalent Network simplification of the embodiment of the present invention;
Fig. 3 is the process flow diagram of the area power grid Analysis of Short-Circuit Current method specific embodiment of the embodiment of the present invention;
Fig. 4 is the network structure Region dividing schematic diagram of the embodiment of the present invention;
Fig. 5 be the embodiment of the present invention each region between the schematic diagram of line admittance of interconnection;
Fig. 6 is the network connection figure of certain node in borderline region after the Equivalent Simplification of the embodiment of the present invention;
Fig. 7 is that the network structure of the IEEE39 node of the embodiment of the present invention divides schematic diagram;
Fig. 8 is the structured flowchart of the area power grid Analysis of Short-Circuit Current device of the embodiment of the present invention.
Embodiment
For making the object, technical solutions and advantages of the present invention clearly understand, below in conjunction with embodiment and accompanying drawing, the present invention is described in further details.At this, exemplary embodiment of the present invention and illustrating for explaining the present invention, but not as a limitation of the invention.
Provide a kind of area power grid Analysis of Short-Circuit Current in embodiments of the present invention, as shown in Figure 1, the method comprises the following steps:
Step 101: Region dividing is carried out to network system: regional classification Water demand being calculated short circuit current level is survey region, being borderline region by there is the regional classification of directly getting in touch with survey region, is perimeter by the regional classification except survey region and borderline region;
Step 102: borderline region and perimeter are carried out equivalent simplification, forms the borderline region after equivalence, obtains the network system after the equivalence that the borderline region after by survey region and equivalence forms;
Step 103: according to the network system after equivalence, the short circuit current level of analytical calculation survey region.
In the above-described embodiments, it is survey region that Water demand in network system is calculated the regional classification of short circuit current level, then be borderline region by there is the regional classification of directly getting in touch with survey region, be perimeter by the regional classification except survey region and borderline region, then be the network system of the borderline region after only comprising survey region and equivalence by this network system Equivalent Simplification, finally by the short circuit current level of the network system analytical calculation survey region after equivalence, because no longer there is perimeter in the network system after Equivalent Simplification, therefore required when analytical calculation short circuit current level data volume will obviously reduce.Efficiently solve by the way in prior art and be difficult to while guarantee calculation of short-circuit current precision, the technical matters of data volume during minimizing Analysis of Short-Circuit Current, reach on the basis ensureing calculation of short-circuit current precision, data volume during minimizing Analysis of Short-Circuit Current and the technique effect of calculated amount.
In one embodiment, the Equivalent Simplification of network system can be carried out by the bus admittance matrix of network system, in one embodiment, borderline region and perimeter are carried out equivalent simplification by above-mentioned steps 102, form the borderline region after equivalence, obtain the network system after the equivalence that the borderline region after by survey region and equivalence forms, can as shown in Figure 2, comprise the following steps:
Step 201: according to the bus admittance matrix of the sequential build network system of perimeter, borderline region and survey region;
Step 202: the part utilizing Gaussian elimination method cancellation bus admittance matrix peripheral region, forms the equivalent posterior nodal point admittance matrix only retaining borderline region and survey region;
Step 203: the variable quantity matrix of the matrix element that borderline region is corresponding in certainty equivalents posterior nodal point admittance matrix;
Step 204: according to the connection branch road of the generator increased in the borderline region after described variable quantity matrix certainty equivalents, increase, and the reactance parameter of the reactance parameter of generator increased and the connection branch road of increase;
Step 205: according to the generator of the increase determined, the connection branch road of increase, the reactance parameter of the reactance parameter of the generator of increase, the connection branch road of increase and generator carry out network reconfiguration with the connected mode be connected between branch road and form the network system after equivalence.
Namely, first build the bus admittance matrix of network system, then matrix element corresponding to bus admittance matrix peripheral region is eliminated by Gaussian elimination method, after Gaussian elimination, matrix element corresponding to perimeter will by equivalence in matrix element corresponding to borderline region, the matrix element that in bus admittance matrix after finally can finding Equivalent Simplification, only borderline region is corresponding there occurs change, can directly according to the generator increased in borderline region in the electric network composition after the change certainty equivalents of matrix element corresponding to borderline region be connected branch road, and increase generator and the reactance parameter being connected branch road.When carrying out Equivalent Simplification and reconstruct network system, original node can be kept in borderline region and to be connected branch road all constant, but by newly-increased generator be connected branch road and add relevant position in borderline region to.
Concrete, the reactance parameter of generator increased according to the borderline region after described variable quantity matrix certainty equivalents and the reactance parameter of the connection branch road of increase can comprise:
1) node number at the connection branch road two ends increased is determined; Get the negative inverse of element corresponding with the node number at two ends in variable quantity matrix, as the reactance parameter of the connection branch road increased;
2) node number of the node be connected with the generator increased is determined; By this node number in the admittance of the connection branch road be directly connected with the node of this node number in perimeter, variable quantity matrix be expert on each element obtain resultant admittance as sum operation; Get the reactance parameter of inverse as the generator increased of resultant admittance.
In order to be specifically described the principle of above-mentioned steps 201 to 205, below in conjunction with a concrete example, the area power grid Analysis of Short-Circuit Current method that the embodiment of the present invention provides is specifically described.
In the present embodiment, area power grid Analysis of Short-Circuit Current method as shown in Figure 3, comprises the following steps:
Step 301: Region dividing
Concrete, choose Water demand and calculate the area of short circuit current level as survey region, there is the area directly got in touch with as borderline region with survey region, the area except survey region and borderline region is as perimeter.After network divides, survey region is connected by interconnection with borderline region, and borderline region is connected by interconnection with perimeter, and survey region and perimeter then do not exist to be contacted directly.
Step 302: networked examination
Concrete, according to the sequential build total system bus admittance matrix Y of perimeter, borderline region and survey region.Such as: Gaussian elimination method can be utilized to carry out the Equivalent Simplification of network, with the part in cancellation Y peripheral region, form the bus admittance matrix after the equivalence only retaining borderline region and survey region
Step 303: borderline region network reconfiguration
The equivalent posterior nodal point admittance matrix generated in interconnection parameter between research on utilization region and borderline region and above-mentioned steps 302
in Partial Elements, borderline region network is reconstructed.Wherein, above-mentioned border networks reconstruct comprises: interconnection reconstruct and generator reconstruct, also needs to calculate the reactance parameter of newly-increased element (comprise newly-increased generator and be connected branch road with newly-increased) in network after reconstructing simultaneously.
Step 304: carry out Analysis of Short-Circuit Current according to the network after reconstruct.
Below above-mentioned steps 302 and step 303 are specifically described:
In above-mentioned steps 302 networked examination, when carrying out analytical calculation to electrical network, in order to reduce the workload of calculator memory and reduction personnel, the part carrying out analyzing in network, can not needed to carry out waiting the simplify processes of value.Be the schematic diagram of network structure Region dividing as shown in Figure 4, for system as shown in Figure 4, the set of node of survey region represents with A, and the set of node of borderline region represents with B, and the set of node of perimeter represents with C.
According to the order of perimeter, borderline region and survey region, define system-wide bus admittance matrix Y, according to bus admittance matrix and the reciprocal feature of nodal impedance matrix, shown in the equation formulations 1 of partitioned matrix form:
Wherein, Y represents bus admittance matrix, and Z represents nodal impedance matrix, E representation unit matrix.
At present, what power industry mainly adopted when calculating short-circuit current is equivalent voltage source method, and the short-circuit impedance namely by solving trouble spot calculates short-circuit current value.According to the character of nodal impedance matrix, element Z on its diagonal line
kkfor the resulting impedance over the ground of entering viewed from whole network from node k, also claim self-impedance or input impedance, the Z when electrical network generation metallic short circuit
kkbe short-circuit impedance.
Due to when carrying out area power grid calculation of short-circuit current, only needing the short circuit current level being concerned about each plant stand in survey region, namely only needing to calculate partitioned matrix Z
aAelement on middle diagonal line, therefore, the computing formula of its correspondence as shown in Equation 2.
Utilize Gaussian elimination method cancellation perimeter C, can obtain:
Contrast above-mentioned formula 2 and above-mentioned formula 3 can be found out, the matrix element that after equivalent, in equation, only borderline region is corresponding there occurs change, namely only has Y
bBthere occurs change, Y before and after equivalent
bBvariable quantity be:
By above-mentioned derivation, Y before and after Equivalent Simplification
bBvariation delta Y
bBcharacterize the impact of external node on internal node, the character according to matrix transpose can draw Δ Y
bBand Y
bBfor symmetric matrix, therefore borderline region can be reconstructed into the electric power networks that with the addition of newly-increased branch road.
In the process of above-mentioned steps 303 borderline region network reconfiguration, according to the situation of change of the equivalence course interior joint admittance matrix derived above, input data due to calculation of short-circuit current software are model parameters of power system component, therefore the matrix form after by equivalence is needed to convert actual power network wiring structure to, i.e. network reconfiguration problem.Network reconfiguration in calculation of short-circuit current process comprises interconnection reconstruct and generator reconstruct, also needs the reactance parameter calculating newly-increased element simultaneously.
For the system shown in above-mentioned Fig. 4, specifically describe the concrete methods of realizing of borderline region network reconfiguration in equivalence course.Called after b1, b2, the b3 respectively of 3 nodes in borderline region, between each region, the line admittance of interconnection as shown in Figure 5.
For the node b1 in borderline region, the method that the interconnection reconstruct after introducing Equivalent Network in borderline region network reconfiguration and generator reconstruct, solution procedure comprises the following steps:
S1: before equivalence, Y
bBthe self-admittance of middle b1 node is:
Y
bB(b1, b1)=y1+y4+y5+y6(formula 5)
The transadmittance of b1 node is:
S2: after equivalence, by the known matrix Y of formula 4
bBthe row element at middle b1 node place exists newly-increased element, suppose the Δ Y calculated by formula 4
bBthe element that middle b1 is expert at is:
Δ Y
bB(b1)=[Δ y1 Δ y2 Δ y3] (formula 7)
Can be drawn by formula 5 to 7, matrix after equivalent
the row element at middle b1 node place is:
According to requirement and the Δ Y of networked examination
bB,
feature, can in the hope of the network connection of Equivalent Simplification posterior nodal point b1 as shown in Figure 6, as can be seen from the transadmittance part in formula 8, between node b1 and node b2, newly-increased admittance is the-circuit of Δ y2, between node b1 and node b3, newly-increased admittance is the-circuit of Δ y3, as can be seen from the self-admittance part of formula 8, after equivalent there is new branch road over the ground in b1 node, from the short-circuit impedance concept in calculation of short-circuit current, branch road equivalence can become generator branch road over the ground, and the subtranient reactance of newly-increased generator is X ' '
d(being converted to admittance is y ' ').
As seen from Figure 6, for b1 node, the newly-increased line parameter circuit value to b2, b3 node can be determined by the transadmittance part in formula 7, then can not directly obtain for equivalent generator reactance parameter, can be drawn by the network structure in Fig. 6, after equivalent, the self-admittance of b1 node is:
In contrast equation 8 and formula 9
expression formula, can obtain:
Y ' '=(y5+y6)+(Δ y1+ Δ y2+ Δ y3) (formula 10)
Be that node b1 to be connected with perimeter the admittance sum (the admittance sum being eliminated circuit be namely connected with node b1) of circuit as can be seen from the Section 1 in formula 10, y ' ' expression formula, Section 2 is matrix Δ Y
bBeach element sum that interior joint b1 is expert at.
As can be seen from above-mentioned analysis, only need the variable quantity of the matrix element that borderline region is corresponding in the admittance matrix after calculating Equivalent Simplification, just the reactance parameter of the generator be connected by interconnection with each node in borderline region can be determined simply, and the reactance parameter of the connection branch road increased between each node in borderline region, thus realize network reconfiguration.
Can draw newly-increased branch road relevant to node b1 in borderline region and the reactance parameter of newly-increased generator according to above-mentioned Fig. 6, formula 8 and formula 10, concrete parameter is as shown in table 1:
Table 1
The processing mode of borderline region interior joint b2, b3 is identical with the processing mode of node b1, does not repeat them here.
Can be found out by above-mentioned Equivalent Simplification step, only need calculate Δ Y
bB, just simply can obtain the reactance parameter that borderline region increases element newly.As can be seen from above-mentioned formula 4, Δ Y
bBsolve in need partitioned matrix Y
cCcarry out inversion operation.In the electrical network of reality, consider that the scale of borderline region is general less, be generally 4 to 6 nodes, therefore can be carried out the inversion operation of alternative matrix by the method for solving of linear equation.In one embodiment, in certainty equivalents posterior nodal point admittance matrix, the variable quantity matrix of the matrix element that borderline region is corresponding can comprise: the number determining borderline region interior joint; When the number of borderline region interior joint is less than predetermined threshold, the method adopting linear equation to solve calculates variable quantity matrix.Wherein, this predetermined threshold can value be 4 to 6, and concrete chooses according to actual needs, and the present invention is not construed as limiting this.That is, the mode solved by linear equation substitutes the mode of matrix inversion, effectively can reduce calculated amount like this when nodes is less.
Below for IEEE39 node, analyzing calculation of short-circuit current external network Equivalent Simplification method, as shown in Figure 7, is that the network structure of this IEEE39 node divides schematic diagram.
As shown in Figure 7, borderline region comprises: new14, new18 and new25 tri-nodes, wherein, and reference capacity SB=100MVA, reference voltage V B=Vav=345kV, according to: the order of new14, new18, new25 calculates Δ Y
bB, result of calculation as shown in Equation 11:
By the computing formula in above-mentioned steps S3, according to Δ Y in formula 11
bBexploitation to obtain the reactance value of element and each newly-increased element that each boundary node newly increases as shown in table 2.
Table 2
PSD-SCCP software is adopted to carry out three short circuit current scanning to the network before and after equivalent, do not consider load, line resistance, line charging power and transformer non-standard no-load voltage ratio during calculating, before and after equivalent, in survey region, each node short-circuit current scanning result is as shown in table 3.
Table 3
| Plant stand | Before equivalence/kA | After equivalent/kA | Difference/kA |
| new1 | 9.783 | 9.783 | 0 |
| new2 | 15.820 | 15.819 | -0.001 |
| new3 | 12.428 | 12.426 | -0.002 |
| new4 | 11.331 | 11.330 | -0.001 |
| new5 | 11.599 | 11.599 | 0 |
| new6 | 11.950 | 11.950 | 0 |
| new7 | 9.511 | 9.511 | 0 |
| new8 | 9.893 | 9.893 | 0 |
| new9 | 9.472 | 9.472 | 0 |
| new10 | 11.047 | 11.047 | 0 |
| new11 | 10.770 | 10.770 | 0 |
| new12 | 4.657 | 4.656 | 0 |
| new13 | 10.396 | 10.395 | -0.001 |
As can be seen from above-mentioned table 3, before and after equivalent, only the short-circuit current value of node new2, new3, new4, new13 has small difference, further analysis can be found out that these four nodes and borderline region exist and directly get in touch with, its short-circuit current error can be summed up as round-off error, when namely adding the reactance of newly-increased branch road and newly-increased equivalent generator in PSD-SCCP software, calculating data are rounded up, thus causes the calculation of short-circuit current value of neighbouring node to occur small difference.The method being demonstrated the embodiment of the present invention provided by the analysis result on IEEE39 node system can ensure the computational accuracy of short-circuit current.
State on the invention in embodiment and give in area power grid calculation of short-circuit current process, Equivalent Simplification process is carried out to external electrical network, and after Equivalent Simplification, reconstruct the power network wiring form of borderline region, the form of newly-increased circuit and newly-increased generator is become to be incorporated in borderline region by external network equivalence, thus form the electric network data only retaining survey region and borderline region, thus significantly reduce the electrical network scale of calculation of short-circuit current, reduce the requirement to calculator memory and simulation software, improve the work efficiency of operation of power networks analyst.
Based on same inventive concept, additionally provide a kind of area power grid Analysis of Short-Circuit Current device in the embodiment of the present invention, as described in the following examples.The principle of dealing with problems due to area power grid Analysis of Short-Circuit Current device is similar to area power grid Analysis of Short-Circuit Current method, therefore the enforcement of area power grid Analysis of Short-Circuit Current device see the enforcement of area power grid Analysis of Short-Circuit Current method, can repeat part and repeats no more.Following used, term " unit " or " module " can realize the software of predetermined function and/or the combination of hardware.Although the device described by following examples preferably realizes with software, hardware, or the realization of the combination of software and hardware also may and conceived.Fig. 8 is a kind of structured flowchart of the area power grid Analysis of Short-Circuit Current device of the embodiment of the present invention, as shown in Figure 8, comprising: Region dividing module 801, Equivalent Simplification module 802 and analysis module 803, be described this structure below.
Region dividing module 801, for carrying out Region dividing to network system: regional classification Water demand being calculated short circuit current level is survey region, being borderline region by there is the regional classification of directly getting in touch with survey region, is perimeter by the regional classification except survey region and borderline region;
Equivalent Simplification module 802, for borderline region and perimeter are carried out equivalent simplification, forms the borderline region after equivalence, obtains the network system after the equivalence that the borderline region after by survey region and equivalence forms;
Analysis module 803, for according to the network system after equivalence, the short circuit current level of analytical calculation survey region.
In one embodiment, above-mentioned Equivalent Simplification module 802 comprises: construction unit, for the bus admittance matrix of the sequential build network system according to perimeter, borderline region and survey region; Simplified element, for utilizing the part in Gaussian elimination method cancellation bus admittance matrix peripheral region, forms the equivalent posterior nodal point admittance matrix only retaining borderline region and survey region; Matrix determining unit, for the variable quantity matrix of matrix element corresponding to borderline region in certainty equivalents posterior nodal point admittance matrix; Parameter determination unit, for the connection branch road according to the generator increased in the borderline region after described variable quantity matrix certainty equivalents, increase, and the reactance parameter of the reactance parameter of generator increased and the connection branch road of increase; Reconfiguration unit, for the connection branch road according to the generator increased, increase, the reactance parameter of the reactance parameter of the generator of increase, the connection branch road of increase and generator carry out network reconfiguration with the connected mode be connected between branch road and form the network system after equivalence.
In one embodiment, described parameter determination unit comprises: first segment period determination subelement, for determining the node number at the connection branch road two ends increased; First parameter determination subelement, for getting the negative inverse of element corresponding with the node number at two ends in variable quantity matrix, as the reactance parameter of the connection branch road increased.
In one embodiment, described parameter determination unit comprises: second section period determination subelement, for the node number of the node that the generator determined with increase is connected; Summing elements, for by this node number in the admittance of the connection branch road be directly connected with the node of this node number in perimeter, variable quantity matrix be expert on each element obtain resultant admittance as sum operation; Second parameter determination subelement, for getting the reactance parameter of inverse as the generator increased of resultant admittance.
In one embodiment, matrix determining unit comprises: nodes determination subelement, for determining the number of borderline region interior joint; Computation subunit, for when the number of borderline region interior joint is less than predetermined threshold, the method adopting linear equation to solve calculates variable quantity matrix.In specifically implementing at one, above-mentioned predetermined threshold can value be 6.
In another embodiment, additionally provide a kind of software, this software is for performing the technical scheme described in above-described embodiment and preferred implementation.
In another embodiment, additionally provide a kind of storage medium, store above-mentioned software in this storage medium, this storage medium includes but not limited to: CD, floppy disk, hard disk, scratch pad memory etc.
From above description, can find out, the embodiment of the present invention achieves following technique effect: it is survey region that Water demand in network system is calculated the regional classification of short circuit current level, then be borderline region by there is the regional classification of directly getting in touch with survey region, be perimeter by the regional classification except survey region and borderline region, then be the network system of the borderline region after only comprising survey region and equivalence by this network system Equivalent Simplification, finally by the short circuit current level of the network system analytical calculation survey region after equivalence, because no longer there is perimeter in the network system after Equivalent Simplification, therefore required when analytical calculation short circuit current level data volume will obviously reduce.Efficiently solve by the way in prior art and be difficult to while guarantee calculation of short-circuit current precision, the technical matters of data volume during minimizing Analysis of Short-Circuit Current, reach on the basis ensureing calculation of short-circuit current precision, data volume during minimizing Analysis of Short-Circuit Current and the technique effect of calculated amount.
Obviously, those skilled in the art should be understood that, each module of the above-mentioned embodiment of the present invention 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, and in some cases, step shown or described by can performing with the order be different from herein, 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 embodiment of 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 embodiment of 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 (10)
1. an area power grid Analysis of Short-Circuit Current method, is characterized in that, comprising:
Region dividing is carried out to network system: regional classification Water demand being calculated short circuit current level is survey region, being borderline region by there is the regional classification of directly getting in touch with survey region, is perimeter by the regional classification except survey region and borderline region;
According to the bus admittance matrix of the sequential build network system of perimeter, borderline region and survey region;
Utilize the part in Gaussian elimination method cancellation bus admittance matrix peripheral region, form the equivalent posterior nodal point admittance matrix only retaining borderline region and survey region;
The variable quantity matrix of the matrix element that borderline region is corresponding in certainty equivalents posterior nodal point admittance matrix;
According to the connection branch road of the generator increased in the borderline region after described variable quantity matrix certainty equivalents, increase, and the reactance parameter of the reactance parameter of generator increased and the connection branch road of increase;
According to the generator of the increase determined, the connection branch road of increase, the reactance parameter of the generator of increase and the reactance parameter of connection branch road increased carry out network reconfiguration form equivalence after network system;
According to the network system after equivalence, the short circuit current level of analytical calculation survey region.
2. the method for claim 1, is characterized in that, according to the reactance parameter of the connection branch road increased in the borderline region after described variable quantity matrix certainty equivalents, comprising:
Determine the node number at the connection branch road two ends increased;
Get the negative inverse of element corresponding with the node number at two ends in variable quantity matrix, as the reactance parameter of the connection branch road increased.
3. the method for claim 1, is characterized in that, according to the reactance parameter of the generator increased in the borderline region after described variable quantity matrix certainty equivalents, comprising:
Determine the node number of the node be connected with the generator increased;
By this node number in the admittance of the connection branch road be directly connected with the node of this node number in perimeter, variable quantity matrix be expert on each element obtain resultant admittance as sum operation;
Get the reactance parameter of inverse as the generator increased of resultant admittance.
4. the method for claim 1, is characterized in that, the variable quantity matrix of the matrix element that borderline region is corresponding in certainty equivalents posterior nodal point admittance matrix, comprising:
Determine the number of borderline region interior joint;
When the number of borderline region interior joint is less than predetermined threshold, the method adopting linear equation to solve calculates variable quantity matrix.
5. method as claimed in claim 4, it is characterized in that, described predetermined threshold is 6.
6. an area power grid Analysis of Short-Circuit Current device, is characterized in that, comprising:
Region dividing module, for carrying out Region dividing to network system: regional classification Water demand being calculated short circuit current level is survey region, being borderline region by there is the regional classification of directly getting in touch with survey region, is perimeter by the regional classification except survey region and borderline region;
Construction unit, for the bus admittance matrix of the sequential build network system according to perimeter, borderline region and survey region;
Simplified element, for utilizing the part in Gaussian elimination method cancellation bus admittance matrix peripheral region, forms the equivalent posterior nodal point admittance matrix only retaining borderline region and survey region;
Matrix determining unit, for the variable quantity matrix of matrix element corresponding to borderline region in certainty equivalents posterior nodal point admittance matrix;
Parameter determination unit, for the connection branch road according to the generator increased in the borderline region after described variable quantity matrix certainty equivalents, increase, and the reactance parameter of the reactance parameter of generator increased and the connection branch road of increase;
Reconfiguration unit, for the connection branch road according to the generator increased, increase, the reactance parameter of the generator of increase and the reactance parameter of connection branch road increased carry out network reconfiguration form equivalence after network system;
Analysis module, for according to the network system after equivalence, the short circuit current level of analytical calculation survey region.
7. device as claimed in claim 6, it is characterized in that, described parameter determination unit comprises:
First segment period determination subelement, for determining the node number at the connection branch road two ends increased;
First parameter determination subelement, for getting the negative inverse of element corresponding with the node number at two ends in variable quantity matrix, as the reactance parameter of the connection branch road increased.
8. device as claimed in claim 6, it is characterized in that, described parameter determination unit comprises:
Second section period determination subelement, for the node number of the node that the generator determined with increase is connected;
Summing elements, for by this node number in the admittance of the connection branch road be directly connected with the node of this node number in perimeter, variable quantity matrix be expert on each element obtain resultant admittance as sum operation;
Second parameter determination subelement, for getting the reactance parameter of inverse as the generator increased of resultant admittance.
9. device as claimed in claim 6, it is characterized in that, described matrix determining unit comprises:
Nodes determination subelement, for determining the number of borderline region interior joint;
Computation subunit, for when the number of borderline region interior joint is less than predetermined threshold, the method adopting linear equation to solve calculates variable quantity matrix.
10. device as claimed in claim 9, it is characterized in that, described predetermined threshold is 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310511797.7A CN103529275B (en) | 2013-10-25 | 2013-10-25 | Area power grid Analysis of Short-Circuit Current method and apparatus |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310511797.7A CN103529275B (en) | 2013-10-25 | 2013-10-25 | Area power grid Analysis of Short-Circuit Current method and apparatus |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103529275A CN103529275A (en) | 2014-01-22 |
| CN103529275B true CN103529275B (en) | 2016-01-20 |
Family
ID=49931446
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310511797.7A Active CN103529275B (en) | 2013-10-25 | 2013-10-25 | Area power grid Analysis of Short-Circuit Current method and apparatus |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103529275B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104391825A (en) * | 2014-11-10 | 2015-03-04 | 南昌大学 | Method of quickly solving node impedance matrix of electric power system based on Gaussian elimination method |
| CN105260949B (en) * | 2015-10-22 | 2019-07-09 | 江苏省电力公司南京供电公司 | A kind of distribution network short circuit calculation method based on matrixing |
| CN108306270B (en) * | 2018-02-27 | 2019-08-06 | 华翔翔能电气股份有限公司 | Photo-voltaic power supply power distribution network short circuit current control method, intelligent terminal and storage medium |
| CN109390943B (en) * | 2018-12-27 | 2021-09-28 | 国电南瑞科技股份有限公司 | Whole-network short-circuit current analysis method and system based on external network online equivalence |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3748521B2 (en) * | 2001-05-30 | 2006-02-22 | 三菱電機株式会社 | Load current detecting method and apparatus and electric power steering apparatus |
| CN100541212C (en) * | 2007-11-08 | 2009-09-16 | 国网南京自动化研究院 | Grid disturbance on-line automatic identification method based on phasor measuring technique |
| CN102937675B (en) * | 2012-10-22 | 2015-12-02 | 江苏省电力公司常州供电公司 | The method of work of Real-time Electrical Distribution Network Data collection and fault location system |
| CN203054144U (en) * | 2012-12-03 | 2013-07-10 | 长安大学 | Range unit for measuring short circuit of transmission line |
| CN103076533B (en) * | 2012-12-28 | 2015-08-12 | 华北电力大学（保定） | The analytical approach of line parameter circuit value and fault disturbance in electrical network |
-
2013
- 2013-10-25 CN CN201310511797.7A patent/CN103529275B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN103529275A (en) | 2014-01-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Milano et al. | Dynamic REI equivalents for short circuit and transient stability analyses | |
| Medjroubi et al. | Open data in power grid modelling: new approaches towards transparent grid models | |
| Trias | The holomorphic embedding load flow method | |
| Lee et al. | Transmission expansion by branch-and-bound integer programming with optimal cost-capacity curves | |
| Shu et al. | A parallel transient stability simulation for power systems | |
| Almeida et al. | Harmonic state estimation through optimal monitoring systems | |
| Saleh et al. | Optimal microgrids placement in electric distribution systems using complex network framework | |
| Lin et al. | Division algorithm and interconnection strategy of restoration subsystems based on complex network theory | |
| CN103529275B (en) | Area power grid Analysis of Short-Circuit Current method and apparatus | |
| CN101719182A (en) | Parallel partition electromagnetic transient digital simulation method of AC and DC power system | |
| Yuan et al. | Graph computing-based wls fast decoupled state estimation | |
| Zhao et al. | Graph-based Preconditioning Conjugate Gradient Algorithm for" N-1" Contingency Analysis | |
| CN109066650B (en) | Power system vulnerability assessment method and terminal equipment | |
| Deakin et al. | Stochastic hosting capacity in LV distribution networks | |
| Jayawardene et al. | A graph theory-based clustering method for power system networks | |
| CN101404040B (en) | Computation resource partition method for power system real-time simulation based on subgraph isomorphism | |
| Ding et al. | Graph theory based splitting strategies for power system islanding operation | |
| Barros et al. | A fog model for dynamic load flow analysis in smart grids | |
| Alayande et al. | An alternative algorithm for solving generation‐to‐load matching and loss allocation problems | |
| Zhang et al. | An efficient approach for robust SCOPF considering load and renewable power uncertainties | |
| Fu et al. | Efficient Topology Assessment for Integrated Transmission and Distribution Network with 10,000+ Inverter-based Resources | |
| Al‐Gabalawy | Deep analysis of the influence of the different power system structures on the performance of the energy storage systems | |
| Alghamdi et al. | PV hosting capacity calculation using particle swarm optimization | |
| Parihar et al. | Power flow analysis of balanced radial distribution system with composite load model | |
| US20210392050A1 (en) | Methods and Systems for Evaluating Data Transportability in Distribution Grids |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |