CN110086209A - The determination method and device of operation states of electric power system - Google Patents
The determination method and device of operation states of electric power system Download PDFInfo
- Publication number
- CN110086209A CN110086209A CN201910373763.3A CN201910373763A CN110086209A CN 110086209 A CN110086209 A CN 110086209A CN 201910373763 A CN201910373763 A CN 201910373763A CN 110086209 A CN110086209 A CN 110086209A
- Authority
- CN
- China
- Prior art keywords
- node
- admittance matrix
- electric system
- matrix
- shunt capacitance
- 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.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/26—Arrangements for eliminating or reducing asymmetry in polyphase networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/50—Arrangements for eliminating or reducing asymmetry in polyphase networks
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The embodiment of the present invention provides a kind of determination method and device of operation states of electric power system, firstly, the number to generator node, shunt capacitance node and balance nodes carries out static optimization, the first admittance matrix is generated, and half dynamic optimization is carried out, generate the second admittance matrix；Then, the size for judging the nonzero element number generated in the first admittance matrix and the second admittance matrix, the admittance matrix for selecting the nonzero element number generated less carry out Load flow calculation to electric system；Finally, determining the operating status of electric system according to calculation of tidal current.It can be seen that, by comparing using the nonzero element number generated in the admittance matrix after static optimization and half dynamic optimization, the admittance matrix for selecting the nonzero element number generated less carries out Load flow calculation to electric system, it can be improved Load flow calculation efficiency, and then guarantee to quickly determine out the operation states of electric power system of less node.
Description
Technical field
The present invention relates to Power System Analysis technical field more particularly to a kind of determination methods of operation states of electric power system
And device.
Background technique
When handling electric system problem, such as: when planning power grid, overhauling, require to determine electric system
Operating status, to guarantee that electric system can operate normally.When electric system when something goes wrong, it is necessary to quickly determine electric power
The operating status of system, finds out rapidly problem, to solve the problems, such as electric system, keeps power system recovery normal.
And in practical applications, the electric system for less node, the side of the operating status due to determining electric system
The limitation of method can quickly determine out the operating status of the electric system sometimes, and can not just quickly determine out sometimes should
The operating status of electric system.As it can be seen that the electric system for less node, is that not can guarantee to quickly determine out electric system
Operating status.Therefore, the operating status for the electric system that one can guarantee to quickly determine out less node is found
The problem of method is current urgent need to resolve.
Summary of the invention
In view of the above problems, the purpose of the embodiment of the present invention is that provide a kind of operation states of electric power system determination method and
Device can guarantee the operating status for quickly determining out the electric system of less node.
In a first aspect, the embodiment of the present invention provides a kind of determination method of operation states of electric power system, which comprises
Determine generators in power systems node location, shunt capacitance node location, balance nodes position, incoming line parameter and transformation
Device parameter；Based on the generator node location, the shunt capacitance node location and the balance nodes position, to the hair
The number of motor node, the shunt capacitance node and the balance nodes carries out static optimization, and according to the incoming line
Parameter and the transformer parameter generate the first admittance matrix of the electric system；Based on the generator node location, institute
Shunt capacitance node location and the balance nodes position are stated, to the generator node, the shunt capacitance node and described
The numbers of balance nodes carries out half dynamic optimization, and according to the incoming line parameter and the transformer parameter, described in generation
Second admittance matrix of electric system；Judge whether the nonzero element number generated in first admittance matrix is greater than described
The nonzero element number generated in two admittance matrixs；If so, being carried out according to second admittance matrix to the electric system
Load flow calculation；If it is not, carrying out Load flow calculation to the electric system then according to first admittance matrix；According to Load flow calculation
As a result, determining the operating status of the electric system.
Second aspect, the embodiment of the present invention provide a kind of determining device of operation states of electric power system, and described device includes:
Determining module, be configured to determine that generators in power systems node location, shunt capacitance node location, balance nodes position,
Incoming line parameter and transformer parameter；First optimization module is configured as based on the generator node location, the parallel connection
The generator node, the shunt capacitance node and the balance are saved in capacitive node position and the balance nodes position
The number of point carries out static optimization, and according to the incoming line parameter and the transformer parameter, generates the electric system
The first admittance matrix；Second optimization module is configured as based on the generator node location, shunt capacitance node position
Set with the balance nodes position, to the number of the generator node, the shunt capacitance node and the balance nodes into
Half dynamic optimization of row, and according to the incoming line parameter and the transformer parameter, generate the electric system second is led
Receive matrix；It is described to be configured as judging whether the nonzero element number generated in first admittance matrix is greater than for judgment module
The nonzero element number generated in second admittance matrix；First computing module is configured as if so, according to second admittance
Matrix carries out Load flow calculation to the electric system；Second computing module is configured as if it is not, then according to first admittance
Matrix carries out Load flow calculation to the electric system；Execution module is configured as determining the electricity according to calculation of tidal current
The operating status of Force system.
The third aspect, the embodiment of the present invention provide a kind of electronic equipment, and the electronic equipment includes: at least one processing
Device；And at least one processor, the bus being connected to the processor；Wherein, the processor, memory pass through described total
Line completes mutual communication；The processor is used to call program instruction in the memory, with execute said one or
Method in multiple technical solutions.
Fourth aspect, the embodiment of the present invention provide a kind of computer readable storage medium, and the storage medium includes storage
Program, wherein equipment where controlling the storage medium in described program operation executes said one or multiple technical sides
Method in case.
The determination method and device of operation states of electric power system provided in an embodiment of the present invention, sends out in determining electric system
After motor node location, shunt capacitance node location, balance nodes position, incoming line parameter and transformer parameter, firstly, base
In generator node location, shunt capacitance node location and balance node location, to generator node, shunt capacitance node peace
The number of weighing apparatus node carries out static optimization, and according to incoming line parameter and transformer parameter, generate electric system first is led
Receive matrix, and, based on generator node location, shunt capacitance node location and balance node location, to generator node, simultaneously
The number for joining capacitive node and balance nodes carries out half dynamic optimization, and according to incoming line parameter and transformer parameter, generates
Second admittance matrix of electric system；Then, judge whether the nonzero element number generated in the first admittance matrix is greater than second
The nonzero element number generated in admittance matrix；If so, carrying out Load flow calculation to electric system according to the second admittance matrix；
If it is not, carrying out Load flow calculation to electric system then according to the first admittance matrix；Finally, determining electric power according to calculation of tidal current
The operating status of system.As it can be seen that by comparing using the non-zero generated in the admittance matrix after static optimization and half dynamic optimization
Element number, the admittance matrix for selecting the nonzero element number generated less carry out Load flow calculation to electric system.Due to non-zero
The computational efficiency of the less admittance matrix of element is high, therefore the efficiency that the electric power system tide that can be improved less node calculates,
And then can guarantee to quickly determine out the operation states of electric power system of less node, so that the electric system of less node is asked
Topic can be handled timely and effectively.
Detailed description of the invention
By reading the following detailed description of the preferred embodiment, various other advantages and benefits are common for this field
Technical staff will become clear.The drawings are only for the purpose of illustrating a preferred embodiment, and is not considered as to the present invention
Limitation.And throughout the drawings, the same reference numbers will be used to refer to the same parts.In the accompanying drawings:
Fig. 1 is the flow diagram one of the determination method of the operation states of electric power system in the embodiment of the present invention；
Fig. 2 is the flow diagram two of the determination method of the operation states of electric power system in the embodiment of the present invention；
Fig. 3 is the schematic diagram of the π type equivalent circuit in the embodiment of the present invention；
Fig. 4 is the network diagram of the electric system in the embodiment of the present invention；
Fig. 5 is the network of 5 node power systems in the embodiment of the present invention；
Fig. 6 is the structural schematic diagram of the determining device of the operation states of electric power system in the embodiment of the present invention；
Fig. 7 is the structural schematic diagram of the electronic equipment in the embodiment of the present invention.
Specific embodiment
The exemplary embodiment that the present invention will be described in more detail below with reference to accompanying drawings.Although showing the present invention in attached drawing
Exemplary embodiment, it being understood, however, that may be realized in various forms the present invention without should be by embodiments set forth here
It is limited.It is to be able to thoroughly understand the present invention on the contrary, providing these embodiments, and can be by the scope of the present invention
It is fully disclosed to those skilled in the art.
The embodiment of the invention provides a kind of determination methods of operation states of electric power system, in practical applications, above-mentioned true
Determining method is mainly to be executed by sensor.Fig. 1 is the determination method of the operation states of electric power system in the embodiment of the present invention
Flow diagram one, shown in Figure 1, this method may include:
S101: generators in power systems node location, shunt capacitance node location, balance nodes position, input are determined
Line parameter circuit value and transformer parameter.
Herein, electric system by power plant, send power transformation route, power supply and distribution and the electric energy that forms of the links such as electricity consumption it is raw
Production and consumption system.
In order to determine generators in power systems node location, shunt capacitance node location, balance nodes position, input
Line parameter circuit value and transformer parameter can obtain the initial data of electric system from the department of keeping electric system initial data,
These initial data can be the parameters of electric power of each link and the topological structure of electric system etc. in electric system.Obtaining electric power
After the initial data of system, generators in power systems node location, simultaneously can be got by the initial data of electric system
Join capacitive node position, balance nodes position, incoming line parameter and transformer parameter, so as to the admittance square of following needs system
The formation of battle array.
S102: based on generator node location, shunt capacitance node location and balance node location, to generator node,
The number of shunt capacitance node and balance nodes carries out static optimization, and according to incoming line parameter and transformer parameter, generates
First admittance matrix of electric system.
It, can by the corresponding admittance matrix of electric system since each electric system is corresponding with an admittance matrix
Load flow calculation is carried out to the electric system, and then determines the operating status of the electric system.And admittance matrix is that height is sparse
Matrix, the sparse degree of the admittance matrix of electric system will have a direct impact on the computational efficiency of admittance matrix, and then influence electric power
The Load flow calculation efficiency of system, and then influence the determination efficiency of operation states of electric power system.Therefore, the sparse of admittance matrix is improved
Degree is highly important for the determination efficiency for improving operation states of electric power system.
And one of improving the core key point of the sparse degree of admittance matrix is exactly to save to node each in electric system
Point number optimization.Incoming line parameter and transformer ginseng based on the node after progress node numbering optimization and the electric system
Number, the admittance matrix after generating optimization.The sparse degree of admittance matrix can be improved to the greatest extent, and then improves electric system
Load flow calculation efficiency, and then improve operation states of electric power system determination efficiency.
It in practical applications, can be using static optimization to electric system in the case of less for electric system interior joint
In each node number carry out node numbering optimization.So-called static optimization counts electric power exactly before node is numbered
Each node goes out line number in system.So-called line number out, is exactly the quantity for the branch that present node is connect with adjacent node.Then it presses
Node each in electric system is numbered in the sequence of line number from less to more out.The number of each node carries out in electric system
After static optimization, in conjunction with the incoming line parameter and transformer parameter of electric system, it will be able to obtain after static optimization
The first admittance matrix.
It should be understood that so-called electric system interior joint is less, can be when the number of nodes in electric system is less than certain
When one preset value, so that it may assert that the node in the electric system is less.The preset value can be 10,20 or 50 etc..For upper
The size of preset value is stated, can be not specifically limited herein depending on the actual complex situation of electric system.
Although static optimization is for promoting the sparse degree of admittance matrix however, the logic of static optimization is simple
Effect is poor.So other than considering the mode of static optimization, it is also necessary to consider other optimal ways to each in electric system
The number of node optimizes.
S103: based on generator node location, shunt capacitance node location and balance node location, to generator node,
The number of shunt capacitance node and balance nodes carries out half dynamic optimization, and according to incoming line parameter and transformer parameter, raw
At the second admittance matrix of electric system.
Herein, other than carrying out node numbering optimization to node each in electric system by the way of static optimization, also
Node numbering optimization can be carried out to node each in electric system by the way of half dynamic optimization.So-called half dynamic optimization, just
It is first then to be eliminated, then correct still unnumbered node only to the least node of line number is numbered out in electric system
Go out line number, at least go out the corresponding node of line number by revised and be numbered again, and eliminated, and so on, Zhi Dao electricity
Until all nodes are all numbered in Force system.The number of each node carries out half dynamic optimization in electric system
Afterwards, in conjunction with the incoming line parameter and transformer parameter of electric system, it will be able to obtain second after more than half dynamic optimizations
Admittance matrix.
It should be understood that step S102 and step S103 have no sequencing in execution, it can be and first carry out step
S102 executes step S103 again, is also possible to first carry out step S103 and executes step S102 again, can also be and is performed simultaneously step
S102 and step S103, is not specifically limited herein.
Although half dynamic optimization is complicated compared to the logic of static optimization, half dynamic optimization is to each in electric system
The effect of optimization of node serial number is preferable, can be improved the sparse degree of admittance matrix, and then improve the computational efficiency of admittance matrix,
The final determination efficiency for improving operation states of electric power system.
However, in some cases, using half dynamic optimization when carrying out node numbering optimization to node each in electric system
Effect instead than using static optimization effect it is poor.
Such as: for wherein there is certain electric system A of 10 nodes, node volume is carried out to node each in electric system A
Number optimization, the effect using static optimization is just better than using the effect of half dynamic optimization.And for wherein exist 5 nodes certain
Electric system B is carried out node numbering optimization to node each in electric system B, is not just used using the effect of half dynamic optimization
The effect of static optimization is good.
Therefore, for wherein existing for the electric system of less node, a kind of optimal node numbering optimization side is selected
Formula for improving the computational efficiency of admittance matrix, and then improves electric power system tide computational efficiency, and then improves and determine power train
The efficiency for operating status of uniting, is very necessary.
S104: judge whether the nonzero element number that generates in the first admittance matrix is greater than and generated in the second admittance matrix
Nonzero element number；If so, executing S105 and S107；If it is not, then executing S106 and S107.
S105: according to the second admittance matrix, Load flow calculation is carried out to electric system.
S106: according to the first admittance matrix, Load flow calculation is carried out to electric system.
S107: according to calculation of tidal current, the operating status of electric system is determined.
Herein, since the calculation of tidal current according to electric system determines that the operating status of electric system is existing skill
Art, therefore for how to determine according to the calculation of tidal current of electric system the operating status of electric system, details are not described herein.
Due to the nonzero element generated in optimized admittance matrix number number, influence whether generate lead
Receive the sparse degree of matrix.Specifically, the number of the nonzero element generated in admittance matrix is fewer, the sparse degree of admittance matrix
It is higher, cause the computational efficiency of admittance matrix higher.Therefore, the admittance matrix for selecting the nonzero element number generated few, into
The Load flow calculation of row electric system can be improved the Load flow calculation efficiency of electric system, and then improve the fortune for determining electric system
The efficiency of row state, so that the problem of enabling the electric system of less node is timely and effectively handled.
As shown in the above, the determination method of operation states of electric power system provided in an embodiment of the present invention is determining electricity
Generator node location, shunt capacitance node location, balance nodes position, incoming line parameter and transformer parameter in Force system
Afterwards, firstly, based on generator node location, shunt capacitance node location and balance node location, to generator node, electricity in parallel
The number for holding node and balance nodes carries out static optimization, and according to incoming line parameter and transformer parameter, generates power train
First admittance matrix of system, and, based on generator node location, shunt capacitance node location and balance node location, to hair
The number of motor node, shunt capacitance node and balance nodes carries out half dynamic optimization, and according to incoming line parameter and transformation
Device parameter generates the second admittance matrix of electric system；Then, judging the nonzero element number generated in the first admittance matrix is
It is no to be greater than the nonzero element number generated in the second admittance matrix；If so, being carried out according to the second admittance matrix to electric system
Load flow calculation；If it is not, carrying out Load flow calculation to electric system then according to the first admittance matrix；Finally, according to Load flow calculation knot
Fruit determines the operating status of electric system.As it can be seen that by comparing using in the admittance matrix after static optimization and half dynamic optimization
The nonzero element number of generation, the admittance matrix for selecting the nonzero element number generated less carry out trend meter to electric system
It calculates.Since the computational efficiency of the less admittance matrix of nonzero element is high, the electric power system tide of less node can be improved
The efficiency of calculating, and then can guarantee to quickly determine out the operation states of electric power system of less node, so that less node
The problem of electric system, can timely and effectively be handled.
Further, as the refinement and extension of method shown in Fig. 1, the embodiment of the invention also provides a kind of electric system
Method for determining running state.Fig. 2 is the process signal of the determination method of the operation states of electric power system in the embodiment of the present invention
Figure two, shown in Figure 2, this method may include:
S201: generators in power systems node location, shunt capacitance node location, balance nodes position, input are determined
Line parameter circuit value and transformer parameter.
S202: based on generator node location, shunt capacitance node location and balance node location, generator section is counted
The line number out of point, shunt capacitance node and balance nodes.
Herein, after generator node location, shunt capacitance node location and balance node location has been determined, just quite
In being aware of the network of electric system, according to the network of electric system, generator node, shunt capacitance node just can know that
With the line number out of balance nodes.
However, in order to improve out the statistical efficiency of line number, and then improve the determination efficiency of the operating status of electric system.?
In practical application, generator node and balance nodes can be considered as ground nodes, each node goes out in statistics electric system
When line number, can not count ground nodes goes out line number, and then on the basis of not influencing node optimization effect, additionally it is possible to improve
The statistical efficiency of line number out.
At the same time, go out line number for the ease of each node of statistics, shunt capacitance node can be considered as to the equivalent electricity of π type
Road.Fig. 3 is the schematic diagram of the π type equivalent circuit in the embodiment of the present invention, shown in Figure 3, wherein I, J are node, and R is electricity
Resistance, X is reactance.So, it is possible simplified balance nodes goes out the statistics difficulty of line number, and then improves out the statistical efficiency of line number.
S203: according to generator node, shunt capacitance node and balance nodes go out line number and incoming line parameter and
Transformer parameter generates original admittance matrix.
Wherein, original admittance matrix is the admittance matrix generated without node numbering optimization.
S204: according to generator node, shunt capacitance node and balance nodes go out line number and incoming line parameter and
Transformer parameter generates the first admittance matrix.
Wherein, the first admittance matrix is the node serial number of generator node, shunt capacitance node and balance nodes by quiet
Admittance matrix after state optimization.
Specifically, step S204 may include:
Step 1: according to the sequence of line number out from less to more, to generator node, shunt capacitance node and balance nodes into
Row number obtains the first number group.
It should be understood that these nodes can be numbered in any order for going out the identical node of line number.
And in order to improve the computational efficiency of admittance matrix, it is saved for going out the identical generator node of line number, shunt capacitance node and balance
Point can preferentially be numbered shunt capacitance node, then generator node and balance nodes are numbered.
Illustrative: assuming that the line number that goes out of generator node A is 2, the line number that goes out of generator node B is 3, shunt capacitance section
The line number that goes out of point C is 2, and the line number that goes out of shunt capacitance node D is 4, and the line number that goes out of balance nodes E is 3.According to line number out by less to
More sequences is ranked up each node, obtains: shunt capacitance node C, generator node A, generator node B, balance nodes
E, shunt capacitance node D.It is 1 by shunt capacitance node C number, is 2 by generator node A number, generator node B is numbered
It is 3, is 4 by balance nodes E number, is 5 by shunt capacitance node D number, these numbers are formed the first number group.
Step 2: according to the first number group, incoming line parameter and transformer parameter, the first admittance matrix is generated.
However, for the number for the nonzero element being further reduced in the first admittance matrix, can by generator node and
Balance nodes are all considered as ground nodes, when the number of each node optimizes in electric system, exclude ground nodes, such as
This, can be further reduced the number of the nonzero element in the first admittance matrix.
Specifically, continuing to execute following step three to step 6 after having executed above-mentioned steps one.
Step 3: being based on generator node location, and statistics generator node goes out line number.
Step 4: according to the sequence of line number out from less to more, being numbered generator node, obtains the second number group.
Step 5: the generator in the first number group is numbered, and replaces with corresponding generator number in the second number group.
Step 6: according to replacement result, incoming line parameter and transformer parameter, the first admittance square of electric system is generated
Battle array.
Then a upper example, firstly, the line number that goes out based on generator node A is 2, the line number that goes out of generator node B is 3,
Each node is ranked up according to the sequence of line number out from less to more, is obtained: generator node A, generator node B.Then, will
Generator node A number is 1, is 2 by generator node B number, these numbers are formed the second number group.Then, it uses
Number in second number group is replaced the number of the corresponding node in the first number group, is obtained: the number of generator node A
It is 1, the number of generator node B is 2, and the number of shunt capacitance node C is 1, and the number of shunt capacitance node D is 5, balance section
The number of point E is 4.Finally, generating first according to the node serial number and incoming line parameter and transformer parameter that newly obtain and leading
Receive matrix.
Herein, incoming line parameter includes: the first node I, tail node J, resistance R, reactance of each branch in electric system
X.Transformer parameter includes: no-load voltage ratio K.It is to need to input first node I, tail node J, resistance during forming admittance matrix
R, reactance X and no-load voltage ratio K, can form admittance matrix.Wherein, first node I needs to be greater than tail node J, for ground branch, tail
Node J=0, for non-transformer branch, no-load voltage ratio K=1.
S205: according to generator node, shunt capacitance node and balance nodes go out line number and incoming line parameter and
Transformer parameter generates the second admittance matrix.
Wherein, the second admittance matrix is the node serial number of generator node, shunt capacitance node and balance nodes through more than half
Admittance matrix after dynamic optimization.
Specifically, step S205 may include:
Step 1: storing and recording the least node of line number out, and delete out the least node of line number and with go out
The branch of the least node connection of line number.
Step 2: that continues to count other nodes in electric system in addition to deletion of node goes out line number, most to line number out
Few node is stored and recorded, until node all in electric system has been recorded.
Step 3: according to the sequencing of record, generator node, shunt capacitance node and balance nodes are compiled
Number.
Step 4: according to number result, incoming line parameter and transformer parameter, the second admittance square of electric system is generated
Battle array.
Illustratively, Fig. 4 is the network diagram of the electric system in the embodiment of the present invention, shown in Figure 4, in the electricity
In Force system, there are node 401, node 402, node 403 and nodes 404, wherein the line number that goes out of node 401 is 3, node 402
Go out line number be 2, node 403 go out line number be 2, node 404 go out line number be 1.Firstly, will currently go out the least node of line number
404 are stored and recorded, and deletion of node 404 and branch adjacent thereto, at this point, the outlet of electric system interior joint 401
Number is 2, and the line number that goes out of node 402 is 2, and the line number that goes out of node 403 is 2；Then, selection currently goes out the least node of line number, if
Current at least line number is corresponding with multiple nodes out, then arbitrarily selects a node, i.e., store and record node 401, and
Deletion of node 401 and branch adjacent thereto；And so on, until node 401, node 402, node 403 and node 404 all by
It has been recorded that, be 1 by the number of node 404 according to the sequence (node 404, node 401, node 402, node 403) of record then,
It is 2 by the number of node 401, is 3 by the number of node 402, is 4 by the number of node 403；Finally, according to number result, incoming line
Parameter and transformer parameter generate the second admittance matrix of electric system.
Herein it should be understood that step S203, step S204 and step S205 have no it is successive execute sequence, can be
It is performed simultaneously, is also possible to not be performed simultaneously, be not specifically limited herein.
In addition, the specific generation type of the second admittance matrix is identical as the specific generation type of the first admittance matrix, therefore
This is repeated no more.
After foring original admittance matrix, the first admittance matrix and the second admittance matrix, it is necessary to by the first admittance
The nonzero element number generated in matrix and the second admittance matrix compares, to determine the higher admittance square of sparse degree
Battle array, and then carry out Load flow calculation.
S206: judge whether the nonzero element number that generates in the first admittance matrix is greater than and generated in the second admittance matrix
Nonzero element number；If so, executing S207 and S209；If it is not, then executing S208 and S209.
S207: according to the second admittance matrix, Load flow calculation is carried out to electric system.
S208: according to the first admittance matrix, Load flow calculation is carried out to electric system.
S209: according to calculation of tidal current, the operating status of electric system is determined.
Specifically, S206 may include:
Step 1: the original nonzero element number in original admittance matrix, the first non-zero in the first admittance matrix are counted
The second nonzero element number in element number and the second admittance matrix.
Step 2: subtraction process is carried out to the first nonzero element number and original nonzero element number, obtains first number.
Step 3: subtraction process is carried out to the second nonzero element number and original nonzero element number, obtains second number.
Step 4: judge whether first number is greater than second number；If so, executing S207 and S209；If it is not, then executing
S208 and S209.
Illustratively, it is assumed that original nonzero element number is 4 in original admittance matrix, the first non-zero in the first admittance matrix
Element number is 6, and the second nonzero element number is 8 in the second admittance matrix.By the first nonzero element number 6 and original non-zero entry
Plain number 4 is subtracted each other, and first number 2 is obtained, and the second nonzero element number 8 and original nonzero element number 4 are subtracted each other, obtain second
Number 4, it can be determined that go out less than second number 4 of first number 2, therefore, it is determined the non-zero generated in the second admittance matrix
Element number is less, i.e. the sparse degree of the second admittance matrix is higher, and then carries out electric system according to the second admittance matrix
Load flow calculation determines the operating status of electric system finally according to calculation of tidal current.
And counting the original nonzero element number in original admittance matrix, the first nonzero element in the first admittance matrix
When the second nonzero element number in number and the second admittance matrix, in order to improve statistical efficiency, above-mentioned steps one specifically may be used
To include:
Step A: triangle decomposition is carried out to original admittance matrix, by the upper triangular matrix after original admittance matrix triangle decomposition
In nonzero element number as original nonzero element number.
Step B: triangle decomposition is carried out to the first admittance matrix, by the upper triangular matrix after the first admittance matrix triangle decomposition
In nonzero element number as the first nonzero element number.
Step C: triangle decomposition is carried out to the second admittance matrix, by the upper triangular matrix after the second admittance matrix triangle decomposition
In nonzero element number as the second nonzero element number.
Wherein, step A, step B and step C sequentially has no sequencing executing, and can be and is performed simultaneously, can also be with
It is not to be performed simultaneously, is not specifically limited herein.
Since admittance matrix is symmetrical matrix, pass through the nonzero element in upper triangular matrix in statistics admittance matrix
Number, it will be able to know nonzero element number in admittance matrix, the standard that nonzero element number counts in it ensure that admittance matrix
While true property, additionally it is possible to improve the efficiency of statistics.
In addition, needing after generating the first admittance matrix and the second admittance matrix by the first admittance matrix and second
Admittance matrix is stored, in case subsequent use.And in order to save memory space, it can only store the first admittance matrix triangle point
The upper triangular matrix after upper triangular matrix and the second admittance matrix triangle decomposition after solution.
Specifically, may include: for the first admittance matrix
Step 1: triangle decomposition is carried out to the first admittance matrix, the upper triangular matrix after obtaining triangle decomposition.
Step 2: upper triangular matrix is decomposed into the first array, the second array and third array.
Wherein, the first array is used to store the diagonal element in upper triangular matrix, and the second array is for storing three angular moments
Off-diagonal element in battle array, third array are used to store the number of nonzero element in each row of upper triangular matrix.
Step 3: the first array of storage, the second array and third array.
In practical applications, the first array can be expressed as YD, and the second array can be expressed as Y, and third array can be with table
It is shown as YDS.When needing to calculate the first admittance matrix, it is only necessary to according to YD, Y, YDS, so that it may quickly know the first admittance
The number of each row nonzero element and position in matrix.
It should be understood that the specific storage mode phase of the specific storage mode of the second admittance matrix and the first admittance matrix
Together, therefore details are not described herein.
So far, the process of the entire operating status for determining electric system is just completed.
Finally, illustrating the process for determining the operating status of electric system again with the example of a 5 node power systems.
Fig. 5 is the network of 5 node power systems in the embodiment of the present invention, shown in Figure 5, is deposited in the electric system
There are 5 nodes, and can determine the position of this 5 nodes, and by the initial data of the electric system, can determine
Incoming line parameter and transformer parameter, referring specifically to table 1, table 1 is the incoming line parameter and transformer of 5 node power systems
Parameter.
Table 1
| Branch first node I | Branch tail node J | Resistance R | Reactance X | No-load voltage ratio K |
| 1 | 0 | 0 | -4 | 1 |
| 2 | 0 | 0 | -2 | 1 |
| 2 | 1 | 0.04 | 0.25 | 1 |
| 3 | 0 | 0 | -4 | 1 |
| 3 | 1 | 0.1 | 0.35 | 1 |
| 3 | 2 | 0.08 | 0.3 | 1 |
| 4 | 2 | 0 | 0.015 | 1.05 |
| 5 | 3 | 0 | 0.03 | 1.05 |
Next, according to the position of this 5 nodes, incoming line parameter and transformer parameter, respectively to this 5 nodes
Number carry out static optimization and half dynamic optimization, the after the first admittance matrix and half dynamic optimization after generating static optimization
Two admittance matrixs；First admittance matrix is subjected to triangle decomposition, obtains the upper triangular matrix of the first admittance matrix, and three from this
Array YD, the array Y of off-diagonal element and the number YDS of each row nonzero element that diagonal element is obtained in angular moment battle array, join respectively
It is shown in Table 2, table 3 and table 4, table 2 is the array YD of diagonal element in upper triangular matrix, and table 3 is off-diagonal element in upper triangular matrix
Array Y, table 4 is the number YDS of each row nonzero element in upper triangular matrix, wherein G is the real part of admittance, i.e. conductance,
B is the imaginary part of admittance, i.e. susceptance, so-called admittance, is exactly the general designation of conductance and susceptance, and admittance defines in power electronics
For the inverse of impedance (impedance), symbol Y, unit is Siemens, referred to as western (S).Admittance is also a plural number, by reality
Number part (conductance G) and imaginary part (susceptance B) composition, i.e. Y=G+jB.
Table 2
| G11 | B11 |
| G22 | B22 |
| G33 | B33 |
| G44 | B44 |
| G55 | B55 |
Table 3
| G21 | B21 |
| G31 | B31 |
| G32 | B32 |
| G42 | B42 |
| G53 | B53 |
Table 4
| 0 | 1 | 2 | 1 | 1 |
Likewise, treatment process of second admittance matrix after carrying out triangle decomposition is identical as the first admittance matrix, therefore
This is repeated no more.
Then, pass through YD, Y, YDS of YD, Y, YDS of the first admittance matrix and the second admittance matrix, it will be able to judge
The size of the nonzero element number generated in first admittance matrix and the nonzero element number generated in the second admittance matrix, in turn
According to judging result, the admittance matrix for selecting the nonzero element number generated few carries out Load flow calculation, and then according to Load flow calculation
As a result, determining the operating status of the 5 node power system.
Based on the same inventive concept, as an implementation of the above method, the embodiment of the invention also provides a kind of power trains
The determining device for operating status of uniting.Fig. 6 is that the structure of the determining device of the operation states of electric power system in the embodiment of the present invention is shown
It is intended to, shown in Figure 6, which may include: determining module 601, be configured to determine that generators in power systems section
Point position, shunt capacitance node location, balance nodes position, incoming line parameter and transformer parameter；First optimization module
602, it is configured as based on generator node location, shunt capacitance node location and balance node location, to generator node, simultaneously
The number for joining capacitive node and balance nodes carries out static optimization, and according to incoming line parameter and transformer parameter, generates electricity
First admittance matrix of Force system；Second optimization module 603 is configured as based on generator node location, shunt capacitance node
Position and balance node location carry out half dynamic optimization to the number of generator node, shunt capacitance node and balance nodes, and
According to incoming line parameter and transformer parameter, the second admittance matrix of electric system is generated；Judgment module 604, is configured as
Judge whether the nonzero element number generated in the first admittance matrix is greater than the nonzero element number generated in the second admittance matrix；
First computing module 605 is configured as if so, carrying out Load flow calculation to electric system according to the second admittance matrix；Second meter
Module 606 is calculated, is configured as if it is not, carrying out Load flow calculation to electric system then according to the first admittance matrix；Execution module 607,
It is configured as determining the operating status of electric system according to calculation of tidal current.
Based on previous embodiment, judgment module, be configured as based on generator node location, shunt capacitance node location and
Balance nodes position, and according to incoming line parameter and transformer parameter, the original admittance matrix of electric system is generated, it is original to lead
Matrix of receiving is the admittance matrix without optimization；Count original nonzero element number in original admittance matrix, the first admittance square
The first nonzero element number in battle array and the second nonzero element number in the second admittance matrix；To the first nonzero element
It is several to carry out subtraction process with original nonzero element number, obtain first number；To the second nonzero element number and original non-zero entry
Plain number carries out subtraction process, obtains second number；Judge whether first number is greater than second number.
Based on previous embodiment, judgment module is configured as carrying out triangle decomposition to original admittance matrix, by original admittance
Nonzero element number in upper triangular matrix after matrix triangle decomposition is as original nonzero element number；To the first admittance matrix
Triangle decomposition is carried out, using the nonzero element number in the upper triangular matrix after the first admittance matrix triangle decomposition as the first non-zero
Element number；Triangle decomposition is carried out to the second admittance matrix, it will be in the upper triangular matrix after the second admittance matrix triangle decomposition
Nonzero element number is as the second nonzero element number.
Based on previous embodiment, the first optimization module is configured as based on generator node location, shunt capacitance node position
Set and balance node location, the line number out of statistics generator node, shunt capacitance node and balance nodes；According to line number out by few
To more sequences, generator node, shunt capacitance node and balance nodes are numbered, obtain the first number group；Based on hair
Motor node location, statistics generator node go out line number；According to the sequence of line number out from less to more, generator node is carried out
Number, obtains the second number group；By the generator number in the first number group, corresponding generator in the second number group is replaced with
Number；According to replacement result, incoming line parameter and transformer parameter, the first admittance matrix of electric system is generated.
Based on previous embodiment, incoming line parameter include: the first node, tail node of each branch in electric system, resistance,
Reactance；Transformer parameter includes: no-load voltage ratio.
Based on previous embodiment, which can also include: the first decomposing module, be configured as to the first admittance matrix into
Row triangle decomposition, the upper triangular matrix after obtaining triangle decomposition；Upper triangular matrix is decomposed into the first array, the second array and
Three arrays, wherein the first array is used to store the diagonal element in upper triangular matrix, and the second array is for storing upper triangular matrix
In off-diagonal element, third array is used to store the number of nonzero element in each row of upper triangular matrix；Store the first array,
Two arrays and third array.
Based on previous embodiment, the second optimization module is configured as based on generator node location, shunt capacitance node position
Set and balance node location, the line number out of statistics generator node, shunt capacitance node and balance nodes；It is least to line number out
Node is stored and recorded, and the branch deleting out the least node of line number and connecting with the least node of line number out；Continue
Count other nodes in electric system in addition to deletion of node goes out line number, and the least node of line number out is stored and remembered
Record, until node all in electric system has been recorded；According to the sequencing of record, to generator node, shunt capacitance
Node and balance nodes are numbered；According to number result, incoming line parameter and transformer parameter, the of electric system is generated
Two admittance matrixs.
Based on previous embodiment, the second optimization module is configured as based on generator node location, shunt capacitance node position
Node location is set and balanced, generator node and balance nodes are considered as ground nodes, it is equivalent that shunt capacitance node is considered as π type
Circuit；Statistics is considered as the line number out of the shunt capacitance node after π type equivalent circuit.
Based on previous embodiment, which can also include: the second decomposing module, be configured as to the second admittance matrix into
Row triangle decomposition, the upper triangular matrix after obtaining triangle decomposition；Upper triangular matrix is decomposed into the first array, the second array and
Three arrays, wherein the first array is used to store the diagonal element in upper triangular matrix, and the second array is for storing upper triangular matrix
In off-diagonal element, third array is used to store the number of nonzero element in each row of upper triangular matrix；Store the first array,
Two arrays and third array.
It need to be noted that: the description of apparatus above embodiment, be with the description of above method embodiment it is similar,
With the similar beneficial effect of same embodiment of the method.For undisclosed technical detail in apparatus of the present invention embodiment, please refer to
The description of embodiment of the present invention method and understand.
Based on the same inventive concept, the embodiment of the invention also provides a kind of electronic equipment.Fig. 7 is in the embodiment of the present invention
Electronic equipment structural schematic diagram, shown in Figure 7, which may include: at least one processor 701；With
And at least one processor 702, the bus 703 being connect with processor 701；Wherein, processor 701, memory 702 pass through bus
703 complete mutual communication；Processor 701 is used to call program instruction in memory 702, to execute said one or more
Method in a embodiment.
It need to be noted that: the description of the above electronic equipment embodiment, the description with above method embodiment are classes
As, there is with embodiment of the method similar beneficial effect.For not draped over one's shoulders in the embodiment of the electronic equipment of the embodiment of the present invention
The technical detail of dew please refers to the description of embodiment of the present invention method and understands.
Based on the same inventive concept, the embodiment of the invention also provides a kind of computer readable storage medium, above-mentioned calculating
Machine readable storage medium storing program for executing includes the program of storage, wherein in program operation, equipment where control storage medium executes above-mentioned one
Method in a or multiple embodiments.
It need to be noted that: the description of the above computer readable storage medium embodiment, with above method embodiment
Description be it is similar, have with embodiment of the method similar beneficial effect.Computer-readable for the embodiment of the present invention is deposited
Undisclosed technical detail in the embodiment of storage media, please refers to the description of embodiment of the present invention method and understands.
It should be understood by those skilled in the art that, embodiments herein can provide as method, system or computer program
Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the application
Apply the form of example.Moreover, it wherein includes the computer of computer usable program code that the application, which can be used in one or more,
The computer program implemented in usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.) produces
The form of product.
The application is referring to method, the process of equipment (system) and computer program product according to the embodiment of the present application
Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions
The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs
Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real
The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates,
Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or
The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one
The step of function of being specified in a box or multiple boxes.
In a typical configuration, calculating equipment includes one or more processors (CPU), input/output interface, net
Network interface and memory.
Memory may include the non-volatile memory in computer-readable medium, random access memory (RAM) and/
Or the forms such as Nonvolatile memory, such as read-only memory (ROM) or flash memory (flash RAM).Memory is computer-readable Jie
The example of matter.
Computer-readable medium includes permanent and non-permanent, removable and non-removable media can be by any method
Or technology come realize information store.Information can be computer readable instructions, data structure, the module of program or other data.
The example of the storage medium of computer includes, but are not limited to phase change memory (PRAM), static random access memory (SRAM), moves
State random access memory (DRAM), other kinds of random access memory (RAM), read-only memory (ROM), electric erasable
Programmable read only memory (EEPROM), flash memory or other memory techniques, read-only disc read only memory (CD-ROM) (CD-ROM),
Digital versatile disc (DVD) or other optical storage, magnetic cassettes, tape magnetic disk storage or other magnetic storage devices
Or any other non-transmission medium, can be used for storage can be accessed by a computing device information.As defined in this article, it calculates
Machine readable medium does not include temporary computer readable media (transitory media), such as the data-signal and carrier wave of modulation.
It should also be noted that, the terms "include", "comprise" or its any other variant are intended to nonexcludability
It include so that the process, method, commodity or the equipment that include a series of elements not only include those elements, but also to wrap
Include other elements that are not explicitly listed, or further include for this process, method, commodity or equipment intrinsic want
Element.In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that including element
There is also other identical elements in process, method, commodity or equipment.
It will be understood by those skilled in the art that embodiments herein can provide as method, system or computer program product.
Therefore, complete hardware embodiment, complete software embodiment or embodiment combining software and hardware aspects can be used in the application
Form.It is deposited moreover, the application can be used to can be used in the computer that one or more wherein includes computer usable program code
The shape for the computer program product implemented on storage media (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.)
Formula.
The above is only embodiments herein, are not intended to limit this application.To those skilled in the art,
Various changes and changes are possible in this application.It is all within the spirit and principles of the present application made by any modification, equivalent replacement,
Improve etc., it should be included within the scope of the claims of this application.
Claims (10)
1. a kind of determination method of operation states of electric power system, which is characterized in that the described method includes:
Determine generators in power systems node location, shunt capacitance node location, balance nodes position, incoming line parameter and
Transformer parameter；
Based on the generator node location, the shunt capacitance node location and the balance nodes position, to the power generation
The number of machine node, the shunt capacitance node and the balance nodes carries out static optimization, and is joined according to the incoming line
The several and transformer parameter, generates the first admittance matrix of the electric system；
Based on the generator node location, the shunt capacitance node location and the balance nodes position, to the power generation
The number of machine node, the shunt capacitance node and the balance nodes carries out half dynamic optimization, and according to the incoming line
Parameter and the transformer parameter generate the second admittance matrix of the electric system；
Judge whether the nonzero element number generated in first admittance matrix is greater than to generate in second admittance matrix
Nonzero element number；
If so, carrying out Load flow calculation to the electric system according to second admittance matrix；
If it is not, carrying out Load flow calculation to the electric system then according to first admittance matrix；
According to calculation of tidal current, the operating status of the electric system is determined.
2. the method according to claim 1, wherein the non-zero generated in judgement first admittance matrix
Whether element number is greater than the nonzero element number generated in second admittance matrix, comprising:
Based on the generator node location, the shunt capacitance node location and the balance nodes position, and according to described
Incoming line parameter and the transformer parameter generate the original admittance matrix of the electric system, the original admittance matrix
For without the admittance matrix of optimization；
Count the original nonzero element number in the original admittance matrix, the first nonzero element in first admittance matrix
The second nonzero element number in number and second admittance matrix；
Subtraction process is carried out to the first nonzero element number and the original nonzero element number, obtains first number；
Subtraction process is carried out to the second nonzero element number and the original nonzero element number, obtains second number；
Judge whether first number is greater than second number.
3. according to the method described in claim 2, it is characterized in that, the original non-zero counted in the original admittance matrix
Second in the first nonzero element number and second admittance matrix in element number, first admittance matrix is non-
Neutral element number, comprising:
Triangle decomposition is carried out to the original admittance matrix, it will be in the upper triangular matrix after the original admittance matrix triangle decomposition
Nonzero element number as original nonzero element number；
Triangle decomposition is carried out to first admittance matrix, it will be in the upper triangular matrix after the first admittance matrix triangle decomposition
Nonzero element number as the first nonzero element number；
Triangle decomposition is carried out to second admittance matrix, it will be in the upper triangular matrix after the second admittance matrix triangle decomposition
Nonzero element number as the second nonzero element number.
4. the method according to claim 1, wherein described be based on the generator node location, the parallel connection
The generator node, the shunt capacitance node and the balance are saved in capacitive node position and the balance nodes position
The number of point carries out static optimization, and according to the incoming line parameter and the transformer parameter, generates the electric system
The first admittance matrix, comprising:
Based on the generator node location, the shunt capacitance node location and the balance nodes position, the hair is counted
The line number out of motor node, the shunt capacitance node and the balance nodes；
According to the sequence of line number out from less to more, to the generator node, the shunt capacitance node and the balance nodes
It is numbered, obtains the first number group；
Based on the generator node location, count the generator node goes out line number；
According to the sequence of line number out from less to more, the generator node is numbered, obtains the second number group；
By the generator number in the first number group, corresponding generator number in the second number group is replaced with；
According to replacement result, the incoming line parameter and the transformer parameter, the first admittance of the electric system is generated
Matrix.
5. according to the method described in claim 4, it is characterized in that, the incoming line parameter includes: in the electric system
The first node of each branch, tail node, resistance, reactance；The transformer parameter includes: no-load voltage ratio.
6. according to the method described in claim 4, it is characterized in that, it is described according to replacement result, the incoming line parameter and
The transformer parameter, after the first admittance matrix for generating the electric system, the method also includes:
Triangle decomposition is carried out to first admittance matrix, the upper triangular matrix after obtaining triangle decomposition；
The upper triangular matrix is decomposed into the first array, the second array and third array, wherein first array is for depositing
The diagonal element in the upper triangular matrix is stored up, second array is used to store the nondiagonal element in the upper triangular matrix
Element, the third array are used to store the number of nonzero element in each row of the upper triangular matrix；
Store first array, second array and the third array.
7. the method according to claim 1, wherein described be based on the generator node location, the parallel connection
The generator node, the shunt capacitance node and the balance are saved in capacitive node position and the balance nodes position
The number of point carries out half dynamic optimization, and according to the incoming line parameter and the transformer parameter, generates the power train
Second admittance matrix of system, comprising:
Based on the generator node location, the shunt capacitance node location and the balance nodes position, the hair is counted
The line number out of motor node, the shunt capacitance node and the balance nodes；
The least node of line number out is stored and recorded, and delete it is described go out the least node of line number and with it is described go out line number
The branch of least node connection；
That continues to count other nodes in the electric system in addition to deletion of node goes out line number, to the least node of line number out
It is stored and recorded, until node all in the electric system has been recorded；
According to the sequencing of record, the generator node, the shunt capacitance node and the balance nodes are compiled
Number；
According to number result, the incoming line parameter and the transformer parameter, the second admittance of the electric system is generated
Matrix.
8. the method according to the description of claim 7 is characterized in that described be based on the generator node location, the parallel connection
Capacitive node position and the balance nodes position count the generator node, the shunt capacitance node and the balance
Node goes out line number, comprising:
Based on the generator node location, the shunt capacitance node location and the balance nodes position, by the power generation
Machine node and the balance nodes are considered as ground nodes, and the shunt capacitance node is considered as π type equivalent circuit；
Statistics is considered as the line number out of the shunt capacitance node after π type equivalent circuit.
9. the method according to the description of claim 7 is characterized in that it is described according to number result, the incoming line parameter and
The transformer parameter, after the second admittance matrix for generating the electric system, the method also includes:
Triangle decomposition is carried out to second admittance matrix, the upper triangular matrix after obtaining triangle decomposition；
The upper triangular matrix is decomposed into the first array, the second array and third array, wherein first array is for depositing
The diagonal element in the upper triangular matrix is stored up, second array is used to store the nondiagonal element in the upper triangular matrix
Element, the third array are used to store the number of nonzero element in each row of the upper triangular matrix；
Store first array, second array and the third array.
10. a kind of determining device of operation states of electric power system, which is characterized in that described device includes:
Determining module is configured to determine that generators in power systems node location, shunt capacitance node location, balance nodes position
It sets, incoming line parameter and transformer parameter；
First optimization module is configured as based on the generator node location, the shunt capacitance node location and described flat
Weigh node location, carries out static optimization to the number of the generator node, the shunt capacitance node and the balance nodes,
And according to the incoming line parameter and the transformer parameter, the first admittance matrix of the electric system is generated；
Second optimization module is configured as based on the generator node location, the shunt capacitance node location and described flat
Weigh node location, and it is excellent to carry out half dynamic to the number of the generator node, the shunt capacitance node and the balance nodes
Change, and according to the incoming line parameter and the transformer parameter, generates the second admittance matrix of the electric system；
Judgment module is configured as judging whether the nonzero element number generated in first admittance matrix is greater than described second
The nonzero element number generated in admittance matrix；
First computing module is configured as if so, carrying out trend meter to the electric system according to second admittance matrix
It calculates；
Second computing module is configured as if it is not, carrying out trend meter to the electric system then according to first admittance matrix
It calculates；
Execution module is configured as determining the operating status of the electric system according to calculation of tidal current.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910373763.3A CN110086209B (en) | 2019-05-07 | 2019-05-07 | Method and device for determining running state of power system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910373763.3A CN110086209B (en) | 2019-05-07 | 2019-05-07 | Method and device for determining running state of power system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110086209A true CN110086209A (en) | 2019-08-02 |
| CN110086209B CN110086209B (en) | 2020-12-22 |
Family
ID=67418919
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910373763.3A Active CN110086209B (en) | 2019-05-07 | 2019-05-07 | Method and device for determining running state of power system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110086209B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102231521A (en) * | 2011-06-24 | 2011-11-02 | 中国电力科学研究院 | Power grid operation state identification method in distribution network self-healing control |
| CN102832616A (en) * | 2012-08-29 | 2012-12-19 | 国核电力规划设计研究院 | Node voltage calculation method and device for power system |
| CN104901308A (en) * | 2015-06-26 | 2015-09-09 | 太原理工大学 | Power system critical circuit identification method |
-
2019
- 2019-05-07 CN CN201910373763.3A patent/CN110086209B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102231521A (en) * | 2011-06-24 | 2011-11-02 | 中国电力科学研究院 | Power grid operation state identification method in distribution network self-healing control |
| CN102832616A (en) * | 2012-08-29 | 2012-12-19 | 国核电力规划设计研究院 | Node voltage calculation method and device for power system |
| CN104901308A (en) * | 2015-06-26 | 2015-09-09 | 太原理工大学 | Power system critical circuit identification method |
Non-Patent Citations (1)
| Title |
|---|
| 邢洁等: "《使潮流方程直接可解的PMU配置方法的改进》", 《电网技术》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110086209B (en) | 2020-12-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110909920B (en) | Power transmission network capacity planning optimization method and system considering multiple fault scenes | |
| Gen et al. | Network design techniques using adapted genetic algorithms | |
| CN104835018B (en) | The reliability estimation method and device of medium voltage distribution network rack | |
| CN106503844B (en) | A kind of power circuit path optimization method using genetic algorithm | |
| Zhao et al. | [Retracted] An Improved SPEA2 Algorithm with Adaptive Selection of Evolutionary Operators Scheme for Multiobjective Optimization Problems | |
| CN106326987A (en) | Multi-objective optimization method and multi-objective optimization device | |
| CN110445167A (en) | A kind of optimization method and system of photovoltaic access distribution | |
| Shang et al. | Spatial scheduling optimization algorithm for block assembly in shipbuilding | |
| CN108764645A (en) | A kind of photo-voltaic power generation station low voltage crossing performance estimating method and system | |
| CN104951442B (en) | A kind of method and apparatus of definitive result vector | |
| Rahimipour Behbahani et al. | Reconfiguration of harmonic polluted distribution network using modified discrete particle swarm optimization equipped with smart radial method | |
| CN105096218B (en) | A kind of topological structure of electric analysis method that distribution archipelago concurrent collaborative calculates | |
| CN110086209A (en) | The determination method and device of operation states of electric power system | |
| CN110516408B (en) | Four-line one-library design management system based on Demo3D simulation | |
| Wang et al. | Pareto optimization of power system reconstruction using NSGA-II algorithm | |
| CN107370149A (en) | Determine the method and system of power system load reduction | |
| CN112183933A (en) | Verification assembly line scheduling method and system based on hybrid variable neighborhood evolution algorithm | |
| CN107437804A (en) | Relate to the check method and device of net definite value | |
| CN110021936A (en) | The method for determining running state of power distribution network containing ring and device | |
| CN114358447B (en) | Flexible interconnection device configuration method and device, electronic equipment and storage medium | |
| CN110428164A (en) | A kind of operation of power networks state monitoring method, device, equipment and readable storage medium storing program for executing | |
| Loukil et al. | Parallel hybrid genetic algorithms for solving Q3AP on computational grid | |
| Saldarriaga-Cortés et al. | Integrated planning of electricity and natural gas systems under uncertain hydro inflows: A multi-objetive approach | |
| Taghavi-farda et al. | Multi-objective group scheduling with learning effect in the cellular manufacturing system | |
| CN109933858B (en) | Core division parallel simulation method for power distribution network |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |