CN106204331A - Embedded electromechanical transient simulation multiport equivalent network component parameters quick calculation method - Google Patents

Abstract

The present invention relates to embedded electromechanical transient simulation multiport equivalent network component parameters quick calculation method, belong to power system electromechanics electromagnetic transient hybrid simulation analysis technical field；The method includes: normal operating condition a certain to electrical network or malfunction, determines power network topology, forms network node voltage equation matrix admittance battle array Y, Y matrix is carried out factorisation, obtains the matrix of consequence A that packing factor is decomposed；Respectively the upper and lower triangular portions of matrix A is solved, obtain solution vector, then reconstruct compression Impedance Matrix M and obtain inverse matrix compression admittance battle array E；According to each element of matrix E, form the equivalent multiport network of the corresponding hybrid simulation subnetting interface practical topology of containing parameter；In dynamo-electric electromagnetism hybrid simulation interface equivalence, self-impedance, mutual impedance, calculating are terminated by practical topology.This method calculating process is without simplifying, and result of calculation is accurate, substantially increases the efficiency of computational methods and calculates speed.

Description

Embedded electromechanical transient simulation multiport equivalent network component parameters quick calculation method

Technical field

The invention belongs to power system electromechanics-electromagnetic transient hybrid simulation analysis technical field, embedded particularly to one Electromechanical transient simulation multiport equivalent network component parameters quick calculation method.

Background technology

Electrical-Machanical Transient Simulation of Power System is mainly used in analyzing the stability of power system, is also used for analyzing power system After being disturbed by certain under one normal operating condition, original operation transition or transition can be returned to after the regular hour Problem to a new steady operational status.At present towards bulk power grid electromechanical transient simulation software or program, mostly use dilute Dredge admittance matrix storage the network information, based on this efficient information, quickly and accurately calculate multiport equivalent network element ginseng Number, is that power system electromechanics-electromagnetic transient hybrid simulation interface equivalence necessary step and link are it can also be used to power system is imitated Middle Equivalent Network simplifies.

Utilizing seam node three short circuit current and steady-state load flow to solve equivalent impedance is also a kind of common method, but Be in systems in practice its steady parameter time use electromotor equivalence electromotive force with system generation three-phase shortcircuit time take Electromotor equivalence transient potential is usually present difference, majority of case can not strict guarantee equivalence electromotive force invariable, Certain deviation is had when solving equivalent impedance.Owing to employing symmetric fault information based on trend impedance, for solving negative phase-sequence With zero sequence equivalent impedance, need just can be conversed by a series of conversion negative phase-sequence and Zero sequence parameter, calculate speed relatively slow, and not Can automatically adapt to actual electric network multiple loaded webs shelf structure.

Summary of the invention

It is an object of the invention to the computational efficiency weak point for overcoming prior art, it is proposed that a kind of embedded electromechanical transient Emulation multiport equivalent network component parameters quick calculation method, this method calculating process is without simplifying, and result of calculation is accurate, significantly Improve the efficiency of computational methods and calculate speed.

The present invention proposes a kind of method that embedded electromechanical transient simulation multiport equivalent network component parameters quickly calculates, Specifically include following steps:

1) note full electric network node number is N, and power system electromechanics-electromagnetic transient hybrid simulation interface quantity is n, at full electricity The sequence number netted in N number of node is designated as { i₁,i₂,…,i_n, n≤N；

2) normal operating condition a certain to full electric network or malfunction, determines power network topology, forms network node voltage side Journey matrix admittance battle array, is designated as Y, as shown in formula (1):

$Y=\left[\begin{array}{ccccc}{Y}_{11}& Y_{12}& Y_{13}& \mathrm{...}& Y_{1N}\\ Y_{21}& Y_{22}& Y_{23}& \mathrm{...}& Y_{2N}\\ Y_{31}& Y_{32}& Y_{33}& \mathrm{...}& Y_{3N}\\ .& .& .& .& .\\ .& .& .& .& .\\ .& .& .& .& .\\ Y_{N1}& Y_{N2}& Y_{N3}& \mathrm{...}& Y_{NN}\end{array}\right]---\left(1\right)$

In formula, Y₁₁Y₂₂Y₃₃…Y_NNFor self-admittance, other be transadmittance；

3) Y matrix is carried out factorisation, obtains the matrix of consequence that packing factor is decomposed, be designated 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 is the upper triangular portions of factor table matrix A；

4) by step 3) L, D, U of obtaining the most simultaneously and concurrently n thread, respectively to former generation based on matrix L, matrix U Back substitution solve, as shown in formula (3):

Yx_k=b_k, k=1,2 ..., n (3)

In formula,DefinitionIt is i-th_kIndividual seam element is 1.0, other element be zero N-dimensional column vector, i_k ∈1,2,…,n；x_kFor corresponding solution vector；

By step 1), 2), 3) Y, b of obtaining_kSubstitute into formula (3), solve and obtain solution vector x₁x₂…x_n；

5) solution vector is made to reconstruct compression Impedance MatrixUtilize Gaussian elimination Method, obtains inverse matrix compression admittance battle array E=M of M^-1；

6) according to each element of matrix E, the equivalence of the corresponding hybrid simulation subnetting interface practical topology forming containing parameter is many Port network；Specifically include:

6-1) to self-impedance:

Electromechanics-electromagnetism hybrid simulation subnetting kth interface to self-impedance, k ∈ 1,2 ..., n, as shown in formula (4):

$Z_{kk}=\frac{1}{E_{kk}+\underset{l=1,l\≠k}{\overset{n}{\Σ}}{E}_{kl}}---\left(4\right)$

In formula (4), E_kk、E_klFor element in matrix E；

6-2) mutual impedance:

If E_kl≠ 0, wherein k ≠ l and E_kl=E_lk, wherein k ≠ l, then kth electromechanics-electromagnetism hybrid simulation subnetting interface and Mutual impedance value between the l electromechanics-electromagnetism hybrid simulation subnetting interface, as shown in formula (5):

$Z_{kl}=-\frac{1}{E_{kl}}---\left(5\right)$

If 6-3) E_kl≠ 0, wherein k ≠ l and E_kl≠E_lk, wherein k ≠ l, then kth electromechanics-electromagnetism hybrid simulation subnetting connects It is an ideal transformer and one group of mutual impedance series connection, mutual impedance between mouth and the l electromechanics-electromagnetism hybrid simulation subnetting interface At kth electromechanics-electromagnetism hybrid simulation interface side, mutual impedance is calculated by formula (5) at l hybrid simulation interface side, resistance value Arrive；

Ideal transformer k side is than l side no-load voltage ratio, as shown in formula (6):

${\mathrm{\ρ}}_{kl}=\sqrt{\frac{E_{kl}}{E_{lk}}}---\left(6\right)$

If 6-4) E_kl=E_lk=0, wherein k ≠ l, then not exchange between kth and the l hybrid simulation subnetting interface Contact, without mutual impedance branch road；

7) the multiport equivalence program of dynamo-electric side it is embedded in practical topology in electromechanics-electromagnetism hybrid simulation interface equivalence etc. Self-impedance, mutual impedance, calculating are terminated by value network element.

The feature of the present invention and beneficial effect:

The present invention uses physical definition based on equivalent network element to calculate each parameter value, and the process that calculates, without simplifying, calculates Result is accurate, and solves measure by the network equation of factorisation in advance and tasks in parallel, substantially increases computational methods Efficiency and calculating speed.Additionally, the present invention can adapt to actual electric network Various Complex grid structure automatically, comprise nothing as in system The asynchronous subnet that alternating current circuit is directly connected to, provides the most effective result.

Accompanying drawing explanation

Fig. 1 is the overall procedure block diagram of the inventive method.

Fig. 2 is the multiport equivalence schematic diagram of the inventive method.

Detailed description of the invention

The one embedded electromechanical transient simulation multiport equivalent network component parameters quick calculation method that the present invention proposes, under Face describes in detail as follows in conjunction with the drawings and the specific embodiments:

The one embedded electromechanical transient simulation multiport equivalent network component parameters quick calculation method that the present invention proposes, its FB(flow block) is as it is shown in figure 1, specifically include following steps:

1) note the whole network node number is N, and power system electromechanics-electromagnetic transient hybrid simulation interface quantity is n, at the whole network N Sequence number in individual node is designated as { i₁,i₂,…,i_n, n≤N；

2) normal operating condition a certain to electrical network or malfunction, determines power network topology, forms network node voltage equation Matrix admittance battle array, is designated as Y, as shown in formula (1):

$Y=\left[\begin{array}{ccccc}{Y}_{11}& Y_{12}& Y_{13}& \mathrm{...}& Y_{1N}\\ Y_{21}& Y_{22}& Y_{23}& \mathrm{...}& Y_{2N}\\ Y_{31}& Y_{32}& Y_{33}& \mathrm{...}& Y_{3N}\\ .& .& .& .& .\\ .& .& .& .& .\\ .& .& .& .& .\\ Y_{N1}& Y_{N2}& Y_{N3}& \mathrm{...}& Y_{NN}\end{array}\right]---\left(1\right)$

In formula, Y₁₁Y₂₂Y₃₃…Y_NNFor self-admittance, other be transadmittance；

3) Y matrix is carried out factorisation, obtains the matrix of consequence that packing factor is decomposed, be designated 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 is the upper triangular portions of factor table matrix A；

4) by step 3) L, D, U of obtaining the most simultaneously and concurrently n thread, respectively to former generation based on matrix L, matrix U Back substitution solve, as shown in formula (3):

Yx_k=b_k, k=1,2 ..., n (3)

In formula,DefinitionIt is i-th_kIndividual seam element is 1.0, other element be zero N-dimensional column vector, i_k ∈1,2,…,n；x_kFor corresponding solution vector；

By step 1), 2), 3) Y, b of obtaining_kSubstitute into formula (3), solve and obtain solution vector x₁x₂…x_n；

5) solution vector is made to reconstruct compression Impedance MatrixUtilize Gaussian elimination Method, obtains inverse matrix compression admittance battle array E=M of M^-1；

6) according to each element of matrix E, the equivalence of the corresponding hybrid simulation subnetting interface practical topology forming containing parameter is many Port network；As in figure 2 it is shown, wherein Z₁₁、Z₂₂、Z₂₂…、Z_kkFor self-impedance, other is mutual impedance；Specifically include:

6-1) to self-impedance:

Electromechanics-electromagnetism hybrid simulation subnetting kth interface to self-impedance, k ∈ 1,2 ..., n, as shown in formula (4):

$Z_{kk}=\frac{1}{E_{kk}+\underset{l=1,l\≠k}{\overset{n}{\Σ}}{E}_{kl}}---\left(4\right)$

In formula (4), E_kk、E_klFor element in matrix E；

6-2) mutual impedance:

If E_kl≠ 0 (k ≠ l) and E_kl=E_lk(k ≠ l), then kth electromechanics-electromagnetism hybrid simulation subnetting interface and l are individual Mutual impedance value between electromechanics-electromagnetism hybrid simulation subnetting interface, as shown in formula (5):

$Z_{kl}=-\frac{1}{E_{kl}}---\left(5\right)$

If 6-3) E_kl≠ 0 (k ≠ l) and E_kl≠E_lk(k ≠ l), then kth electromechanics-electromagnetism hybrid simulation subnetting interface and Being an ideal transformer and one group of mutual impedance series connection between l electromechanics-electromagnetism hybrid simulation subnetting interface, mutual impedance is in kth Individual electromechanics-electromagnetism hybrid simulation interface side, mutual impedance is calculated by formula (5) at l hybrid simulation interface side, resistance value；

Ideal transformer k side is than l side no-load voltage ratio, as shown in formula (6):

${\mathrm{\ρ}}_{kl}=\sqrt{\frac{E_{kl}}{E_{lk}}}---\left(6\right)$

If 6-4) E_kl=E_lk=0 (k ≠ l), then do not exchange connection between kth and the l hybrid simulation subnetting interface Network, without mutual impedance branch road；

7) the multiport equivalence program of dynamo-electric side it is embedded in practical topology in electromechanics-electromagnetism hybrid simulation interface equivalence etc. Self-impedance, mutual impedance, calculating are terminated by value network element.

Claims (1)

1. an embedded electromechanical transient simulation multiport equivalent network component parameters quick calculation method, it is characterised in that the party Method specifically includes following steps:

1) note full electric network node number is N, and power system electromechanics-electromagnetic transient hybrid simulation interface quantity is n, at full electric network N Sequence number in individual node is designated as { i₁,i₂,…,i_n, n≤N；

2) normal operating condition a certain to full electric network or malfunction, determines power network topology, forms network node voltage equation square Battle array admittance battle array, is designated as Y, as shown in formula (1):

$Y=\left[\begin{array}{ccccc}{Y}_{11}& Y_{12}& Y_{13}& \mathrm{...}& Y_{1N}\\ Y_{21}& Y_{22}& Y_{23}& \mathrm{...}& Y_{2N}\\ Y_{31}& Y_{32}& Y_{33}& \mathrm{...}& Y_{3N}\\ .& .& .& .& .\\ .& .& .& .& .\\ .& .& .& .& .\\ Y_{N1}& Y_{N2}& Y_{N3}& \mathrm{...}& Y_{NN}\end{array}\right]---\left(1\right)$

In formula, Y₁₁ Y₂₂ Y₃₃ … Y_NNFor self-admittance, other be transadmittance；

3) Y matrix is carried out factorisation, obtains the matrix of consequence that packing factor is decomposed, be designated 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 that the master of factor table matrix A is right Angle unit part, U is the upper triangular portions of factor table matrix A；

4) by step 3) L, D, U of obtaining the most simultaneously and concurrently n thread, respectively to former generation based on matrix L, matrix U time In generation, solves, as shown in formula (3):

Yx_k=b_k, k=1,2 ..., n (3)

In formula,DefinitionIt is i-th_kIndividual seam element is 1.0, other element be zero N-dimensional column vector, i_k∈1, 2,…,n；x_kFor corresponding solution vector；

By step 1), 2), 3) Y, b of obtaining_kSubstitute into formula (3), solve and obtain solution vector x₁ x₂ … x_n；

5) solution vector is made to reconstruct compression Impedance MatrixUtilize Gaussian elimination method, Inverse matrix to M compresses admittance battle array E=M^-1；

6) according to each element of matrix E, the equivalent multiport of the corresponding hybrid simulation subnetting interface practical topology of containing parameter is formed Network；Specifically include:

6-1) to self-impedance:

Electromechanics-electromagnetism hybrid simulation subnetting kth interface to self-impedance, k ∈ 1,2 ..., n, as shown in formula (4):

$Z_{kk}=\frac{1}{E_{kk}+\underset{1=1,l\≠k}{\overset{n}{\Σ}}{E}_{ki}}---\left(4\right)$

In formula (4), E_kk、E_klFor element, wherein k ≠ l in matrix E,；

6-2) mutual impedance:

If E_kl≠ 0, wherein k ≠ l, and E_kl=E_lk, wherein k ≠ l, then kth electromechanics-electromagnetism hybrid simulation subnetting interface and l Mutual impedance value between individual electromechanics-electromagnetism hybrid simulation subnetting interface, as shown in formula (5):

$Z_{kl}=-\frac{1}{E_{kl}}---\left(5\right)$

If 6-3) E_kl≠ 0, wherein k ≠ l and E_kl≠E_lk, wherein k ≠ l, then kth electromechanics-electromagnetism hybrid simulation subnetting interface and Being an ideal transformer and one group of mutual impedance series connection between l electromechanics-electromagnetism hybrid simulation subnetting interface, mutual impedance is the K electromechanics-electromagnetism hybrid simulation interface side, mutual impedance is calculated by formula (5) at l hybrid simulation interface side, resistance value；

Ideal transformer k side is than l side no-load voltage ratio, as shown in formula (6):

${\mathrm{\ρ}}_{kl}=\sqrt{\frac{E_{kl}}{E_{lk}}}---\left(6\right)$

If 6-4) E_kl=E_lk=0, wherein k ≠ l, then do not exchange contact between kth and the l hybrid simulation subnetting interface, Branch road without mutual impedance；

7) the multiport equivalence program of dynamo-electric side it is embedded in practical topology equivalence net in electromechanics-electromagnetism hybrid simulation interface equivalence Self-impedance, mutual impedance, calculating are terminated by network element.