CN106204331B - Embedded electromechanical transient simulation multiport equivalent network component parameters quick calculation method - Google Patents
Embedded electromechanical transient simulation multiport equivalent network component parameters quick calculation method Download PDFInfo
- Publication number
- CN106204331B CN106204331B CN201610567304.5A CN201610567304A CN106204331B CN 106204331 B CN106204331 B CN 106204331B CN 201610567304 A CN201610567304 A CN 201610567304A CN 106204331 B CN106204331 B CN 106204331B
- Authority
- CN
- China
- Prior art keywords
- matrix
- formula
- interface
- impedance
- hybrid simulation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 47
- 230000001052 transient effect Effects 0.000 title claims abstract description 18
- 238000004364 calculation method Methods 0.000 title claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims abstract description 49
- 238000000034 method Methods 0.000 claims abstract description 16
- 230000006835 compression Effects 0.000 claims abstract description 7
- 238000007906 compression Methods 0.000 claims abstract description 7
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 4
- 230000007257 malfunction Effects 0.000 claims abstract description 4
- 238000012856 packing Methods 0.000 claims abstract description 4
- 230000008030 elimination Effects 0.000 claims description 3
- 238000003379 elimination reaction Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000000819 phase cycle Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The present invention relates to embedded electromechanical transient simulation multiport equivalent network component parameters quick calculation methods, belong to electric system electromechanics-electromagnetic transient hybrid simulation analysis technical field；This method comprises: determining power network topology to a certain normal operating condition of power grid or malfunction, network node voltage equation matrix --- admittance battle array Y is formed, Factorization is carried out to Y matrix, obtains the matrix of consequence A of packing factor decomposition；The upper and lower triangular portions of matrix A are solved respectively, obtain solution vector, then is reconstructed into compression Impedance Matrix M and obtains inverse matrix compression admittance battle array E；According to each element of matrix E, the equivalent multiport network of the correspondence hybrid simulation subnetting interface practical topology of containing parameter is formed；Practical topology terminates self-impedance, mutual impedance, calculating in electromechanics-electromagnetism hybrid simulation interface equivalence.For this method calculating process without simplification, calculated result is accurate, substantially increases the efficiency and calculating speed of calculation method.
Description
Technical field
The invention belongs to electric system electromechanics-electromagnetic transient hybrid simulation analysis technical fields, in particular to a kind of embedded
Electromechanical transient simulation multiport equivalent network component parameters quick calculation method.
Background technique
Electrical-Machanical Transient Simulation of Power System is mainly used for analyzing the stability of electric system, is also used to analyze electric system
After being interfered under one normal operating condition by certain, original operation transition or transition can be returned to after the regular hour
The problem of to a new steady operational status.At present towards bulk power grid electromechanical transient simulation software or program, mostly use greatly dilute
Thin admittance matrix stores the network information, based on this efficient information, quickly and accurately calculates multiport equivalent network element ginseng
Number, is electric system electromechanics-electromagnetic transient hybrid simulation interface equivalence necessary step and link, it can also be used to which electric system is imitative
Very middle Equivalent Network simplifies.
Solving equivalent impedance using interface node three short circuit current and steady-state load flow is also a kind of common method, but
Be in systems in practice its steady parameter when the generator equivalence potential that uses and system occur to take when three-phase shortcircuit
Generator equivalence transient potential is usually present difference, majority of case can not strict guarantee equivalence potential it is invariable,
There is certain deviation when solving equivalent impedance.Due to having used the symmetric fault information based on trend impedance, for solving negative phase-sequence
With zero sequence equivalent impedance, need just converse negative phase-sequence and Zero sequence parameter by a series of transformation, calculating speed is slower, and not
The a variety of load grid structures of actual electric network can be adapted to automatically.
Summary of the invention
The purpose of the present invention is the computational efficiency shortcomings to overcome prior art, propose a kind of embedded electromechanical transient
Multiport equivalent network component parameters quick calculation method is emulated, for this method calculating process without simplification, calculated result is accurate, significantly
Improve the efficiency and calculating speed of calculation method.
The invention proposes a kind of method that embedded electromechanical transient simulation multiport equivalent network component parameters quickly calculate,
Specifically includes the following steps:
1) note full electric network node number is N, and electric system electromechanics-electromagnetic transient hybrid simulation interface quantity is n, in full electricity
The serial number netted in N number of node is denoted as { i1,i2,…,in, n≤N；
2) it to a certain normal operating condition of full electric network or malfunction, determines power network topology, forms network node voltage side
Journey matrix --- admittance battle array is denoted as Y, as shown in formula (1):
In formula, Y11Y22Y33…YNNFor self-admittance, it is other be transadmittance；
3) Factorization is carried out to Y matrix, obtains the matrix of consequence of packing factor decomposition, is denoted as A, as shown in formula (2):
A=L+D+U (2)
In formula, A is the factor table matrix of matrix Y, and L is the lower triangular portions of factor table matrix A, and D is factor table matrix A
Main diagonal element part, U are the upper triangular portions of factor table matrix A；
4) L, D, the U obtained by step 3) concurrent n thread simultaneously respectively, respectively to the former generation based on matrix L, matrix U
Back substitution solve, as shown in formula (3):
Yxk=bk, k=1,2 ..., n (3)
In formula,DefinitionIt is i-thkA interface element is the N-dimensional column vector that 1.0, other elements are zero, ik
∈1,2,…,n；xkFor corresponding solution vector；
By step 1), 2), 3) Y, the b obtainedkIt substitutes into formula (3), solution obtains solution vector x1x2…xn；
5) solution vector is enabled to be reconstructed into compression Impedance MatrixUtilize Gaussian elimination
Method obtains the inverse matrix compression admittance battle array E=M of M-1；
6) according to each element of matrix E, the equivalence for forming the correspondence hybrid simulation subnetting interface practical topology of containing parameter is more
Port network；It specifically includes:
6-1) to self-impedance:
Electromechanics-electromagnetism hybrid simulation k-th of interface of subnetting is to self-impedance, k ∈ 1,2 ..., n, as shown in formula (4):
In formula (4), Ekk、EklFor element in matrix E；
6-2) mutual impedance:
If Ekl≠ 0, wherein k ≠ l and Ekl=Elk, wherein k ≠ l, then k-th of electromechanics-electromagnetism hybrid simulation subnetting interface and
Mutual impedance value between first of electromechanics-electromagnetism hybrid simulation subnetting interface, as shown in formula (5):
If 6-3) Ekl≠ 0, wherein k ≠ l and Ekl≠Elk, wherein k ≠ l, then k-th of electromechanics-electromagnetism hybrid simulation subnetting connects
It connects for an ideal transformer with one group of mutual impedance between mouth and first of electromechanics-electromagnetism hybrid simulation subnetting interface, mutual impedance
In k-th of electromechanics-electromagnetism hybrid simulation interface side, mutual impedance is calculated in l hybrid simulation interface side, impedance value by formula (5)
It arrives；
The side the ideal transformer k side ratio l no-load voltage ratio, as shown in formula (6):
If 6-4) Ekl=Elk=0, wherein k ≠ l, then do not exchange between k-th and first of hybrid simulation subnetting interface
Contact, no mutual impedance branch；
7) the multiport equivalence program of electromechanical side is embedded in practical topology in electromechanics-electromagnetism hybrid simulation interface equivalence etc.
It is worth network element to self-impedance, mutual impedance, calculating terminates.
The features of the present invention and the utility model has the advantages that
The present invention calculates each parameter value using the physical definition based on equivalent network element, and calculating process is calculated without simplification
As a result accurate, and measure is solved by the parallel network equation of preparatory Factorization and task, substantially increase calculation method
Efficiency and calculating speed.In addition, the present invention can adapt to actual electric network Various Complex grid structure automatically, as included nothing in system
The asynchronous subnet that alternating current circuit is directly connected to provides reasonable effective result.
Detailed description of the invention
Fig. 1 is the overall procedure block diagram of the method for the present invention.
Fig. 2 is the multiport equivalence schematic diagram of the method for the present invention.
Specific embodiment
A kind of embedded electromechanical transient simulation multiport equivalent network component parameters quick calculation method proposed by the present invention, under
Detailed description are as follows in conjunction with the drawings and the specific embodiments in face:
A kind of embedded electromechanical transient simulation multiport equivalent network component parameters quick calculation method proposed by the present invention,
Flow diagram as shown in Figure 1, specifically includes the following steps:
1) note the whole network node number is N, and electric system electromechanics-electromagnetic transient hybrid simulation interface quantity is n, in the whole network N
Serial number in a node is denoted as { i1,i2,…,in, n≤N；
2) it to a certain normal operating condition of power grid or malfunction, determines power network topology, forms network node voltage equation
Matrix --- admittance battle array is denoted as Y, as shown in formula (1):
In formula, Y11Y22Y33…YNNFor self-admittance, it is other be transadmittance；
3) Factorization is carried out to Y matrix, obtains the matrix of consequence of packing factor decomposition, is denoted as A, as shown in formula (2):
A=L+D+U (2)
In formula, A is the factor table matrix of matrix Y, and L is the lower triangular portions of factor table matrix A, and D is factor table matrix A
Main diagonal element part, U are the upper triangular portions of factor table matrix A；
4) L, D, the U obtained by step 3) concurrent n thread simultaneously respectively, respectively to the former generation based on matrix L, matrix U
Back substitution solve, as shown in formula (3):
Yxk=bk, k=1,2 ..., n (3)
In formula,DefinitionIt is i-thkA interface element is the N-dimensional column vector that 1.0, other elements are zero, ik
∈1,2,…,n；xkFor corresponding solution vector；
By step 1), 2), 3) Y, the b obtainedkIt substitutes into formula (3), solution obtains solution vector x1x2…xn；
5) solution vector is enabled to be reconstructed into compression Impedance MatrixUtilize Gaussian elimination
Method obtains the inverse matrix compression admittance battle array E=M of M-1；
6) according to each element of matrix E, the equivalence for forming the correspondence hybrid simulation subnetting interface practical topology of containing parameter is more
Port network；As shown in Fig. 2, wherein Z11、Z22、Z22…、ZkkFor self-impedance, other is mutual impedance；It specifically includes:
6-1) to self-impedance:
Electromechanics-electromagnetism hybrid simulation k-th of interface of subnetting is to self-impedance, k ∈ 1,2 ..., n, as shown in formula (4):
In formula (4), Ekk、EklFor element in matrix E；
6-2) mutual impedance:
If Ekl≠ 0 (k ≠ l) and Ekl=Elk(k ≠ l), then k-th of electromechanics-electromagnetism hybrid simulation subnetting interface and first
Mutual impedance value between electromechanics-electromagnetism hybrid simulation subnetting interface, as shown in formula (5):
If 6-3) Ekl≠ 0 (k ≠ l) and Ekl≠Elk(k ≠ l), then k-th of electromechanics-electromagnetism hybrid simulation subnetting interface and
It connects for an ideal transformer with one group of mutual impedance between l electromechanics-electromagnetism hybrid simulation subnetting interface, mutual impedance is in kth
A electromechanics-electromagnetism hybrid simulation interface side, mutual impedance are calculated in l hybrid simulation interface side, impedance value by formula (5)；
The side the ideal transformer k side ratio l no-load voltage ratio, as shown in formula (6):
If 6-4) Ekl=Elk=0 (k ≠ l), then without exchange connection between k-th and first of hybrid simulation subnetting interface
Network, no mutual impedance branch；
7) the multiport equivalence program of electromechanical side is embedded in practical topology in electromechanics-electromagnetism hybrid simulation interface equivalence etc.
It is worth network element to self-impedance, mutual impedance, calculating terminates.
Claims (1)
1. a kind of embedded electromechanical transient simulation multiport equivalent network component parameters quick calculation method, which is characterized in that the party
Method specifically includes the following steps:
1) note full electric network node number is N, and electric system electromechanics-electromagnetic transient hybrid simulation interface quantity is n, in full electric network N
Serial number in a node is denoted as { i1,i2,…,in, n≤N；
2) it to a certain normal operating condition of full electric network or malfunction, determines power network topology, forms network node voltage equation square
Battle array --- admittance battle array is denoted as Y, as shown in formula (1):
In formula, Y11 Y22 Y33 … YNNFor self-admittance, it is other be transadmittance；
3) Factorization is carried out to Y matrix, obtains the matrix of consequence of packing factor decomposition, is denoted as A, as shown in formula (2):
A=L+D+U (2)
In formula, A is the factor table matrix of matrix Y, and L is the lower triangular portions of factor table matrix A, and D is the master couple of factor table matrix A
Angle member part, U are the upper triangular portions of factor table matrix A；
4) L, D, the U obtained by step 3) concurrent n thread simultaneously respectively, respectively to time of former generation, matrix U based on matrix L
In generation, solves, as shown in formula (3):
Yxk=bk, k=1,2 ..., n (3)
In formula,DefinitionIt is i-thkA interface element is the N-dimensional column vector that 1.0, other elements are zero, ik∈1,
2,…,n；xkFor corresponding solution vector；
By step 1), 2), 3) Y, the b obtainedkIt substitutes into formula (3), solution obtains solution vector x1x2…xn；
5) solution vector is enabled to be reconstructed into compression Impedance MatrixUsing Gaussian elimination method, obtain
Inverse matrix to M compresses admittance battle array E=M-1；
6) according to each element of matrix E, the equivalent multiport of the correspondence hybrid simulation subnetting interface practical topology of containing parameter is formed
Network；It specifically includes:
6-1) to self-impedance:
Electromechanics-electromagnetism hybrid simulation k-th of interface of subnetting is to self-impedance, k ∈ 1,2 ..., n, as shown in formula (4):
In formula (4), Ekk、EklFor element in matrix E, wherein k ≠ l,；
6-2) mutual impedance:
If Ekl≠ 0, wherein k ≠ l, and Ekl=Elk, wherein k ≠ l, then k-th of electromechanics-electromagnetism hybrid simulation subnetting interface and l
Mutual impedance value between a electromechanics-electromagnetism hybrid simulation subnetting interface, as shown in formula (5):
If 6-3) Ekl≠ 0, wherein k ≠ l and Ekl≠Elk, wherein k ≠ l, then k-th of electromechanics-electromagnetism hybrid simulation subnetting interface and
It is connected with one group of mutual impedance between first of electromechanics-electromagnetism hybrid simulation subnetting interface for an ideal transformer, mutual impedance is the
K electromechanics-electromagnetism hybrid simulation interface side, mutual impedance are calculated in l hybrid simulation interface side, impedance value by formula (5)；
The side the ideal transformer k side ratio l no-load voltage ratio, as shown in formula (6):
If 6-4) Ekl=Elk=0, wherein k ≠ l, then get in touch between k-th and first of hybrid simulation subnetting interface without exchange,
Without mutual impedance branch；
7) the multiport equivalence program for being embedded in electromechanical side is calculated to practical topology in electromechanics-electromagnetism hybrid simulation interface equivalence etc.
It is worth network element to self-impedance, mutual impedance, calculating terminates.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610567304.5A CN106204331B (en) | 2016-07-19 | 2016-07-19 | Embedded electromechanical transient simulation multiport equivalent network component parameters quick calculation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610567304.5A CN106204331B (en) | 2016-07-19 | 2016-07-19 | Embedded electromechanical transient simulation multiport equivalent network component parameters quick calculation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106204331A CN106204331A (en) | 2016-12-07 |
| CN106204331B true CN106204331B (en) | 2019-07-12 |
Family
ID=57493029
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610567304.5A Active CN106204331B (en) | 2016-07-19 | 2016-07-19 | Embedded electromechanical transient simulation multiport equivalent network component parameters quick calculation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106204331B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111965478B (en) * | 2020-06-30 | 2023-10-20 | 中国电力科学研究院有限公司 | Method and system for determining short-circuit current by considering phase shifting effect of transformer |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102521452A (en) * | 2011-12-14 | 2012-06-27 | 中国电力科学研究院 | Computing system of large power grid closed loop |
| CN102592003A (en) * | 2011-12-15 | 2012-07-18 | 天津大学 | Data exchange method for electromechanical and electromagnetic transient hybrid simulation |
| CN103605829A (en) * | 2013-10-22 | 2014-02-26 | 华中科技大学 | Equivalent model method for performing electromagnetic transient simulation on alternating current/direct current hybrid power grid |
-
2016
- 2016-07-19 CN CN201610567304.5A patent/CN106204331B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102521452A (en) * | 2011-12-14 | 2012-06-27 | 中国电力科学研究院 | Computing system of large power grid closed loop |
| CN102592003A (en) * | 2011-12-15 | 2012-07-18 | 天津大学 | Data exchange method for electromechanical and electromagnetic transient hybrid simulation |
| CN103605829A (en) * | 2013-10-22 | 2014-02-26 | 华中科技大学 | Equivalent model method for performing electromagnetic transient simulation on alternating current/direct current hybrid power grid |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106204331A (en) | 2016-12-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101505061B (en) | Computation method capable of tracking Davinan equivalence parameter base on time domain simulation | |
| CN103559347A (en) | Method for establishing electromagnetic transient simulation model of large-scale AC-DC (Alternating Current - Direct Current) power system | |
| Yu et al. | An unscented particle filtering approach to decentralized dynamic state estimation for DFIG wind turbines in multi-area power systems | |
| CN101527455B (en) | Interconnected electric network distributed current calculating method on the basis of alternation and iteration of current module | |
| CN103500245B (en) | A kind of field road transient state-transient state coupling simulation method based on m ulti-loop m ethod | |
| CN103066593A (en) | Three-phase load flow calculation method of weak looped distribution network comprising multi-type distributed power sources | |
| CN105938578A (en) | Large-scale photovoltaic power station equivalent modeling method based on clustering analysis | |
| CN102708225B (en) | Fragmentation debugging method for electromagnetic transient simulation of large AC (alternating current)-DC (direct current) power grid | |
| CN105205244B (en) | Loop Closing Operation analogue system based on electromechanics-electromagnetism hybrid simulation technology | |
| CN109802392B (en) | Large-scale power distribution network load flow calculation method and device | |
| CN104578049B (en) | A kind of transient power quality analysis system of electromechanical electromagnetic transient hybrid simulation | |
| CN109066669A (en) | Simulating method for electric power system, electric system simulation device and electronic equipment | |
| CN109190320A (en) | Parallel Heterogeneous simulation method suitable for ADPSS double-fed blower closed-loop test | |
| CN104361796B (en) | Electric energy metering device simulation training platform | |
| CN103956735A (en) | Harmonic power flow analysis method of distributed power generation system | |
| Fang et al. | Realization of electromechanical transient and electromagnetic transient real time hybrid simulation in power system | |
| CN106204331B (en) | Embedded electromechanical transient simulation multiport equivalent network component parameters quick calculation method | |
| CN106410811B (en) | Iteration small impedance branches endpoint changes the tidal current computing method of Jacobian matrix for the first time | |
| CN108649597A (en) | The analysis method that a kind of failure moving method and failure influence HVDC commutation failures | |
| CN105375475A (en) | Power system fast transient stability simulation method based on real-time dynamic equivalence | |
| CN111339624B (en) | PSASP and EMTP/ATP-based short-circuit current direct-current component calculation method | |
| CN104750994B (en) | Method for judging transient instability dominance based on outer region power ratio | |
| CN108181835B (en) | A kind of electric system DC grid real-time emulation method and device | |
| Abooshahab et al. | Disturbance and state estimation in partially known power networks | |
| CN107064624A (en) | Alternating-current charging pile electric energy gauging method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20210706 Address after: 100176 11th floor, building 1, No.2 Yuncheng street, Beijing Economic and Technological Development Zone, Daxing District, Beijing Patentee after: BEIJING ZHIZHONG ENERGY INTERNET RESEARCH INSTITUTE Co.,Ltd. Patentee after: China Southern Power Grid Research Institute Co.,Ltd. Patentee after: TSINGHUA University Address before: 100176 11th floor, building 1, No.2 Yuncheng street, Beijing Economic and Technological Development Zone, Daxing District, Beijing Patentee before: BEIJING ZHIZHONG ENERGY INTERNET RESEARCH INSTITUTE Co.,Ltd. Patentee before: CSG POWER GRID TECHNOLOGY RESEARCH CENTER Patentee before: TSINGHUA University |