CN109829216A - A kind of spacer units model equivalent method for power system real-time simulation - Google Patents
A kind of spacer units model equivalent method for power system real-time simulation Download PDFInfo
- Publication number
- CN109829216A CN109829216A CN201910055552.5A CN201910055552A CN109829216A CN 109829216 A CN109829216 A CN 109829216A CN 201910055552 A CN201910055552 A CN 201910055552A CN 109829216 A CN109829216 A CN 109829216A
- Authority
- CN
- China
- Prior art keywords
- equivalent
- coefficient
- node
- spacer units
- multivalue
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 125000006850 spacer group Chemical group 0.000 title claims abstract description 42
- 238000004088 simulation Methods 0.000 title claims abstract description 30
- 239000011159 matrix material Substances 0.000 claims description 13
- 238000012545 processing Methods 0.000 claims description 8
- 230000005611 electricity Effects 0.000 claims description 3
- 238000002347 injection Methods 0.000 claims description 2
- 239000007924 injection Substances 0.000 claims description 2
- 238000013461 design Methods 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 239000000243 solution Substances 0.000 description 6
- 230000009466 transformation Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention discloses a kind of spacer units model equivalent method for power system real-time simulation.This interval model of element equivalent method is internal equivalent to spacer units progress common in electric system, is two-value conductance by the element equivalent in addition to current transformer and ground resistance；It is eliminated in advance using internal node of the variable star angle to spacer units, reduces the dimension of network equation；And technology is prestored using multivalue coefficient, give multivalue equivalent admittance addressing method, multivalue voltage coefficient addressing method and the addressing circuit multiplexing method of spacer units.The method that the present invention designs reduces the time needed for network equation solves, and accelerates simulation velocity, effectively increases power system real-time simulation scale.
Description
Technical field
The invention belongs to technical field of electric power automation, are related to a kind of spacer units model etc. of power system real-time simulation
Efficacious prescriptions method.
Background technique
Electric system simulation analysis is the basis of planning and design of power system and management and running, and wherein real-time simulation can be to two
Secondary device carries out hardware in loop experiment and carries out technical training to power system staff.
State variable technique and the nodal method of analysis are two kinds of classical emulation method for solving.State variable technique passes through prediction correction
Method implicitly solves each differential equation, has larger flexibility, and the solution to non-linear element in the selection of algorithm
Accuracy is higher, but computational efficiency is lower and cannot flexibly switch simulation step length.And the nodal method of analysis is according in electric system
The voltage-current relationship of primary element uses numerical integrating by differential equation equivalence for corresponding concomitant circuit model, and passes through
The nodal method of analysis solves linear circuit, finally obtains the instantaneous voltage value of each node.Its main feature is that simulation efficiency is high and accurate
Property it is high, while can guarantee good numerical stability, have become the mainstream algorithm of power system real-time simulation.
But there are a big defects for the existing nodal method of analysis for power system real-time simulation: existing in electric system big
Measure the spacer units being made of breaker, disconnecting switch, current transformer.The presence of spacer units can greatly increase network equation
Dimension increase simulation time to increase the calculation amount of solution, be unfavorable for realize real-time simulation.
Summary of the invention
Purpose of the invention is to overcome the shortcomings in the prior art, provide it is a kind of design rationally, net can be effectively reduced
Network equation solution calculation amount, to improve the spacer units model equivalent method of power system real-time simulation efficiency.
The purpose of the present invention is what is be achieved through the following technical solutions:
A kind of spacer units model equivalent method for power system real-time simulation, comprising the following steps:
(1) spacer units model is carried out internal equivalent；By the conductance in addition to current transformer and ground resistance with two
It is worth the description of conductance model, does variable star angle by center node of internal node, the internal node of equivalent model is eliminated, thus will
Spacer units model is equivalent to the pure resistance network of only boundary node；
(2) by the external characteristics equation of YU=I describe inside spacer units it is equivalent after simulation model, wherein Y is interval
The equivalent node admittance matrix of unit, U are node voltage vector, and I is node Injection Current vector, and each equivalent admittance in Y is used
Multivalue coefficient prestores method processing, wherein it is that possible take all of each coefficient of node admittance matrix that the multivalue coefficient, which prestores method,
Value prestores in advance, without operational components in simulation process, but the value of each coefficient is directly obtained by addressing circuit
Method；Addressing circuit obtains actually taking for multivalue coefficient by the on state of each conductance inside statistical interval unit equivalent model
Value；
(3) by superposition theorem and homogeneous theorem it is found that spacer units internal node voltages be equal to each boundary node voltage with
The sum of products of corresponding voltage coefficient equally prestores method processing voltage coefficient using multivalue coefficient；
(4) it is sought accordingly for the equivalent node admittance matrix of spacer units and the corresponding multivalue coefficient configuration of voltage coefficient
Location circuit；The equivalent circuit of each phase sub-network of spacer units has symmetry in topology, and identical addressing circuit can be not
It is used when with the addressing of multivalue coefficient；
(5) the corresponding multivalue coefficient of the equivalent node admittance matrix and voltage coefficient of spacer units is pre-stored in storage unit
In, operational components directly can obtain each coefficient value by addressing circuit during real-time simulation.
Compared with prior art, the beneficial effects brought by the technical solution of the present invention are as follows:
1, eliminate spacer units internal node using variable star angle, to spacer units carry out it is internal it is equivalent after only surplus boundary save
Point reduces the dimension of network equation, reduces serial degree when solving network equation.
2, spacer units external characteristics node admittance matrix and voltage coefficient description after will be equivalent, it is pre- using multivalue coefficient
Technical treatment node admittance matrix and voltage coefficient are deposited, operational components pressure is alleviated.
3, for different spacer units, addressing circuit is multiplexed using its symmetry, improves hardware resource
Utilization rate.
Detailed description of the invention
Fig. 1 is the inlet wire interval schematic diagram of the invention being connected with double-bus；
Fig. 2 is spacer units equivalent circuit schematic diagram of the invention；
Fig. 3 is the internal equivalent preceding schematic diagram of spacer units A phase sub-network (hereinafter referred to as sub-network) of the invention；
Fig. 4 is that sub-network of the invention eliminates the schematic diagram after node 4a；
Fig. 5 is the equivalent admittance Y of sub-network of the invention1a,2dAddressing circuit schematic diagram；
Specific embodiment
The present invention is described in further detail below in conjunction with the drawings and specific embodiments.It should be appreciated that described herein
Specific embodiment be only used to explain the present invention, be not intended to limit the present invention.
The present invention proposes a kind of spacer units model equivalent method for power system real-time simulation, including following step
It is rapid:
(1) the problem of a large amount of existing spacer units influence power system real-time simulation speed in electric system is proposed.For
This is internal to the progress of spacer units model equivalent, by the two-value conductance model of the conductance in addition to current transformer and ground resistance
Description.Variable star angle is done by center node of internal node later, the internal node of equivalent model is eliminated, thus will interval list
Meta-model is equivalent to the pure resistance network of only boundary node.
(2) it after equivalent inside the simulation model of spacer units, can be described with the external characteristics equation of type such as YU=I.Node is led
Each equivalent admittance in matrix Y of receiving can prestore method processing with multivalue coefficient.Wherein, it is to lead node that the multivalue coefficient, which prestores method,
All possible values of each coefficient of matrix of receiving prestore in advance, without operational components in simulation process, but pass through addressing electricity
The method that road directly obtains the value of each coefficient.The conducting that addressing circuit passes through conductance each inside statistical interval unit equivalent model
State obtains the practical value of multivalue coefficient.
(3) by superposition theorem and homogeneous theorem it is found that spacer units internal node voltages be equal to each boundary node voltage with
The sum of products of corresponding voltage coefficient.Therefore voltage coefficient also uses multivalue coefficient to prestore method processing.
(4) equivalent circuit of each phase sub-network of spacer units in topology have symmetry, it is identical seek a yard circuit can
To be used in the addressing of different multivalue coefficients.
Function of the invention and effect are described further separately below:
As shown in Figure 1, giving the connected inlet wire spacer units schematic diagram of double-bus, Fig. 2 is its equivalent circuit.Wherein,
G2g、G3g、G4gFor grounding conductance, G34a、G34b、G34cFor the equivalent conductance of current transformer, remaining conductance is all two-value conductance.For
The value of grounding conductance is set arbitrarily, it is equivalent to carry out network internal using node 1,5,6,2d, 3d, 4d as boundary.Before equivalent
Model have 3 independent sub-networks, each sub-network only includes 7 two-value conductances, to only surplus after equivalent inside each sub-network
Boundary node reduces the dimension of network equation, reduces serial degree when solving network equation.
The equivalent process regarding internal node as and eliminating in advance inside sub-network, by doing star by center node of internal node
The pre- cancellation process of angular transformation analysis internal node.Illustrate so that A phase sub-network eliminates the process of node 4a in advance as an example, Fig. 3 is A phase
Equivalent preceding equivalent circuit inside sub-network, Fig. 4 are the equivalent circuit eliminated after node 4a.After internal node 4a is eliminated, node
6a, 5a, 4d, 3a are Angle connection by star-like connection relation transformation before, increase 6 equivalent conductances, the number of each equivalent conductance newly
Value is equal to the sum of all conductances is divided by between the product and internal node and connected node of conductance between internal node and interdependent node, example
Such as G36aCalculation formula
Successively doing variable star angle as center node using node 4a, 3a, 2a can be completed mistake equivalent inside A phase sub-network
Journey.When the values dictate of two-value conductance in sub-network is 10-8/108When S, by the formula analysis of variable star angle it is found that after equivalent
The value of equivalent conductance can up to 128 in sub-network, therefore and the storage pressure of generation can not put up with.Equivalent conductance is such as
The reason of as many as this is that (each two-value conductance is equal to 10 on path when the path between two boundary nodes is not turned on8When S), they it
Between equivalent conductance be a minimum and nonzero value, and the minimum is mutual between any two because being influenced by remaining two-value conductance
Not identical, it is nonsensical for recording these minimums.If the path between regulation internal node and Arbitrary Boundaries node is not led
Its node voltage is zero when logical, and all two-value conductances in Fig. 2 are defined as 0/108S, then when the road between two boundary nodes
When diameter is not turned on, the equivalent conductance between them is only zero, and the value of equivalent conductance each in this way can greatly reduce.
The equivalent rear pure resistance network for only surplus boundary node inside the simulation model of spacer units, can be with type such as YU=I's
The description of external characteristics equation.By analyzing internal equivalent process it can be found that when a certain internal node to equivalent admittance interdependent node
Path conducting when, then will affect the numerical value of the equivalent admittance using the internal node as the variable star angle of center node, boundary is electric
The value for leading (each two-value conductance between internal node and the boundary node that is connected) determines this influence.Utilize this rule, node
Each equivalent admittance in admittance matrix Y can be handled with multivalue Y-factor method Y, with the equivalent admittance Y of A phase sub-network1a,2dFor explanation it is more
The treatment process of value coefficient method.It is 1. boundary node 3d, 4d, 5a, 6a to the path of 1. (node 2a) by 1a to 2d path reference number
2. 3. 4. 5. respectively marked as.The on-off of path 1. determines Y1a,2dIt whether is zero, remaining path represents boundary conductance in star
To the influence of equivalent conductance in angular transformation.Path in homogeneous variable star angle is divided into one kind, then is 2. one kind, is 3. 4. 5.
It is a kind of.According to the sequence that the number of all kinds of middle path conductings is arranged and combined from small to large, each case is calculated in advance
Lower Y1a,2dSpecific value and by multivalue coefficient store, Y1a,2dShared 1+2 × 4=9 value.All kinds of middle guiding paths need to only be counted
The number of diameter can calculate Y1a,2dOffset, addressing circuit is as shown in Figure 5.
All equivalent admittances can be completed to address in aforementioned manners, general addressing formula
Wherein, D is the actual address of the multivalue coefficient, D0For base address, n is the number of class behind route classification, siIt is i-th
The number that path is connected in a class, mjFor the total number in path in j-th of class.When handling equivalent self-admittance addressing, path is 1.
In only self-admittance interdependent node, s1Be constantly equal to is 1.
In order to realize, electric current calculates the function of judging with current over-zero in breaker in current transformer, needs each imitative
True step size computation goes out each internal node voltages of spacer units.By superposition theorem with homogeneous theorem it is found that internal node voltages are equal to
The sum of products of each boundary node voltage and corresponding voltage coefficient.Voltage coefficient also can use the processing of multivalue Y-factor method Y, in terms of
The voltage coefficient B of external node 1a when calculating the voltage of internal node 3a1a,3aFor illustrate treatment process, node will be removed in Fig. 3
All boundary nodes ground connection outside 1a, when node 1a order position ideal voltage source, the voltage and B of node 3a1a,3aNumerically
It is equal.Label is distinguished in the path of node 1a to other boundary nodes, as long as thering is one to lead in all paths comprising node 3a
Lead to, then the voltage non-zero of node 4a, B1a,4aNon-zero.B is calculated according to the step identical as processing equivalent admittance1a,3aAll possibility
Value and by multivalue coefficient store.B1a,3a2 × 2 × 4=16 value is shared, wherein zero stores twice.Zero is deposited
Although storage repeatedly wastes some data spaces, avoid judge voltage coefficient whether the trouble for being zero, make to decode electric
Road is more succinct, and can make the different internal nodes same decoding circuit corresponding with the voltage coefficient of the same external node.
The general accessing formula of multivalue voltage coefficient:
When carrying out power system real-time simulation calculating, operational components need to match for each phase sub-network of different interval unit
Set multiple addressing circuits.When the simulation model of electric system changes, different addressing circuits is consolidated according to conception
Change within hardware.The equivalent circuit of each phase sub-network of spacer units in topology have symmetry, it is identical seek a yard circuit can
To use in the addressing of different multivalue parameters, need to only adjust each time-varying conductance in time-varying influences the position in word.Use multivalue
When parametric method handles equivalent admittance all in sub-network shown in Fig. 3, voltage coefficient, 13 kinds of addressing circuits, table 1 need to be only configured
Give the corresponding relationship between each multivalue parameter and addressing circuit number.
The mapping table of 1 multivalue coefficient of table and addressing circuit number
The present invention is not limited to embodiments described above.Above the description of specific embodiment is intended to describe and say
Bright technical solution of the present invention, the above mentioned embodiment is only schematical, is not restrictive.This is not being departed from
In the case of invention objective and scope of the claimed protection, those skilled in the art may be used also under the inspiration of the present invention
The specific transformation of many forms is made, within these are all belonged to the scope of protection of the present invention.
Claims (1)
1. a kind of spacer units model equivalent method for power system real-time simulation, which comprises the following steps:
(1) spacer units model is carried out internal equivalent；By the two-value electricity of the conductance in addition to current transformer and ground resistance
The description of guided mode type, does variable star angle by center node of internal node, and the internal node of equivalent model is eliminated, thus will interval
Model of element is equivalent to the pure resistance network of only boundary node；
(2) by the external characteristics equation of YU=I describe inside spacer units it is equivalent after simulation model, wherein Y is spacer units
Equivalent node admittance matrix, U is node voltage vector, and I is node Injection Current vector, each equivalent admittance multivalue in Y
Coefficient prestores method processing, wherein it is to propose all possible values of each coefficient of node admittance matrix that the multivalue coefficient, which prestores method,
Before prestore, without operational components in simulation process, but the method for directly obtaining the value of each coefficient by addressing circuit；
Addressing circuit obtains the practical value of multivalue coefficient by the on state of each conductance inside statistical interval unit equivalent model；
(3) by superposition theorem with homogeneous theorem it is found that spacer units internal node voltages are equal to each boundary node voltage and corresponding
The sum of products of voltage coefficient equally prestores method processing voltage coefficient using multivalue coefficient；
(4) for the equivalent node admittance matrix of spacer units and the corresponding addressing electricity of the corresponding multivalue coefficient configuration of voltage coefficient
Road；The equivalent circuit of each phase sub-network of spacer units has symmetry in topology, and identical addressing circuit can be more in difference
Value coefficient uses when addressing；
(5) the corresponding multivalue coefficient of the equivalent node admittance matrix and voltage coefficient of spacer units is pre-stored in storage unit,
Operational components directly can obtain each coefficient value by addressing circuit during real-time simulation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910055552.5A CN109829216A (en) | 2019-01-22 | 2019-01-22 | A kind of spacer units model equivalent method for power system real-time simulation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910055552.5A CN109829216A (en) | 2019-01-22 | 2019-01-22 | A kind of spacer units model equivalent method for power system real-time simulation |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109829216A true CN109829216A (en) | 2019-05-31 |
Family
ID=66861113
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910055552.5A Pending CN109829216A (en) | 2019-01-22 | 2019-01-22 | A kind of spacer units model equivalent method for power system real-time simulation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109829216A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111581738A (en) * | 2020-04-16 | 2020-08-25 | 天津大学 | Component method for real-time simulation element of power system |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080250369A1 (en) * | 2007-04-09 | 2008-10-09 | Ming-Hong Lai | Method of estimating the signal delay in a VLSI circuit |
| CN103207941A (en) * | 2013-04-27 | 2013-07-17 | 清华大学 | Transient analysis method and transient analysis system under integrated circuit power supply network full-parameter model |
| CN104298822A (en) * | 2014-10-09 | 2015-01-21 | 天津大学 | Nonlinearity-considered electromagnetic transient real-time simulation method of transformer substation |
| CN104638771A (en) * | 2015-02-15 | 2015-05-20 | 国家电网公司 | Quantitative analysis method for short-term reliability of process-level network of intelligent substation |
| CN105224754A (en) * | 2015-10-14 | 2016-01-06 | 清华大学 | A kind of simulation of power electronic method based on Interpolation compensation current switch model |
-
2019
- 2019-01-22 CN CN201910055552.5A patent/CN109829216A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080250369A1 (en) * | 2007-04-09 | 2008-10-09 | Ming-Hong Lai | Method of estimating the signal delay in a VLSI circuit |
| CN103207941A (en) * | 2013-04-27 | 2013-07-17 | 清华大学 | Transient analysis method and transient analysis system under integrated circuit power supply network full-parameter model |
| CN104298822A (en) * | 2014-10-09 | 2015-01-21 | 天津大学 | Nonlinearity-considered electromagnetic transient real-time simulation method of transformer substation |
| CN104638771A (en) * | 2015-02-15 | 2015-05-20 | 国家电网公司 | Quantitative analysis method for short-term reliability of process-level network of intelligent substation |
| CN105224754A (en) * | 2015-10-14 | 2016-01-06 | 清华大学 | A kind of simulation of power electronic method based on Interpolation compensation current switch model |
Non-Patent Citations (2)
| Title |
|---|
| BINGDA ZHANG.ETC: ""Modeling of Power System Simulation Based on FRTDS"", 《ENERGIES》 * |
| 王潇;张炳达;陈雄;: "电力系统实时仿真中细粒度并行实现", 天津大学学报(自然科学与工程技术版) * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111581738A (en) * | 2020-04-16 | 2020-08-25 | 天津大学 | Component method for real-time simulation element of power system |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109145957A (en) | The identification and processing method and processing device of power distribution network abnormal index based on big data | |
| CN103928925B (en) | Power distribution network load flow calculation method based on forward-backward sweep | |
| CN103917992B (en) | For the method and apparatus for the synapse weight for storing neutral net using memory by probability mode | |
| CN108628947A (en) | A kind of business rule matched processing method, device and processing equipment | |
| Jubril et al. | Solving multi-objective economic dispatch problem via semidefinite programming | |
| CN106936127B (en) | A kind of line load regression analysis and prediction technique and system | |
| CN108737149A (en) | Network topology real machine simulation method and system | |
| CN106886616A (en) | A kind of automatic subnetting method of extensive electro-magnetic transient grid simulation | |
| CN103872681A (en) | Online real-time loop closing method based on integration of major network and distribution network | |
| CN105354656A (en) | Partition decoupling based distributed parallel computing method and system for distribution network state estimation | |
| CN104636801A (en) | Transmission line audible noise prediction method based on BP neural network optimization | |
| CN107464048A (en) | A kind of plan security check method a few days ago based on research state | |
| CN104318006B (en) | A kind of optimization method of virtual terminals of intelligent substation design | |
| CN104050536B (en) | Topology analysis method for power distribution network crossing voltage grades | |
| CN103793590A (en) | GPU-based computation method for quickly solving power flow in distribution networks | |
| CN106202662A (en) | A kind of distribution network structure figure draws mapping method automatically | |
| CN109474023A (en) | Intelligent distribution network section real time updating method, system, storage medium and terminal | |
| CN103346556B (en) | A kind of distribution network loop method for rapidly positioning | |
| CN110019420A (en) | A kind of data sequence prediction technique and calculate equipment | |
| CN103777092B (en) | A kind of relay protection test system based on cloud and method thereof | |
| CN108471353A (en) | A method of the NE capacity analysis based on deep neural network algorithm and prediction | |
| CN108460455A (en) | Model treatment method and device | |
| CN110082612A (en) | Ultra-high/extra-high voltage converter station DC fields Electric Field Simulation method, apparatus and storage medium | |
| CN112750051A (en) | Random forest algorithm-based phase sequence identification method and device for transformer area and terminal equipment | |
| CN110413722A (en) | Address choice method, apparatus and non-transient storage medium |
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 | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190531 |
|
| WD01 | Invention patent application deemed withdrawn after publication |