CN109980650A - A kind of radial pattern distribution system tidal current computing method - Google Patents
A kind of radial pattern distribution system tidal current computing method Download PDFInfo
- Publication number
- CN109980650A CN109980650A CN201910424449.3A CN201910424449A CN109980650A CN 109980650 A CN109980650 A CN 109980650A CN 201910424449 A CN201910424449 A CN 201910424449A CN 109980650 A CN109980650 A CN 109980650A
- Authority
- CN
- China
- Prior art keywords
- matrix
- node
- branch
- distribution system
- row
- 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
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Abstract
The invention discloses a kind of tidal current computing methods of radial pattern distribution system, Load flow calculation is carried out for the distribution system in the direction of distribution net platform region one, first according to radial pattern distribution system the characteristics of, establish the circuit model of radial pattern distribution system, its topological structure is reasonably numbered, and using matrix form come the primary data of indication circuit, the especially digitization description of circuit topology relationship, make the logical relation clear and definite that figure connects in circuit model, the introducing of incidence matrix realizes the Programmed Design to path search process, and primary data matrix is used to input as unique data, keep entire programming directly simple；Each load bus voltage and each Branch Power Flow of power distribution network are solved in calculating process using the algebraic operation of matrix element, thus the program calculating used time is less, convergence rate is very fast.
Description
Technical field
The present invention relates to low-voltage distribution technical field more particularly to a kind of radial pattern distribution system tidal current computing methods.
Background technique
Currently, distribution power system load flow calculation is the important foundation of distribution network system analysis, economical operation.The network weight of power distribution network
Structure, troubleshooting, idle work optimization and state estimation, line loss analyzing etc. require to use the result of distribution power system load flow calculation.It is a kind of
The tidal current computing method of function admirable is the key that distribution system management.As power department is to the attention degree of power distribution network management
Constantly deepen, the Load flow calculation research specifically for power distribution network has developed extensively.Since power distribution network is that power transmission network and electric energy are used
The intermediate link at family, voltage class is low compared with power transmission network, and network structure is in tree-like, multiple-limb single radial knot when steady-state operation
Structure, route R/X value is higher, and most cases are greater than 1；Low voltage distribution transformer multidigit is in load center, from distribution transformer low pressure
Lateral multiple directions power supply, can be considered as infinitely great power supply, the power supply of multiple directions for distribution system higher level's substation exit bus
It can separate and be calculated, the trend of low-voltage network flows to load from distribution transformer low-pressure side, and the trend on route flows tool
There is one-way.
Domestic and foreign scholars propose various Power Flow Calculation Methods For Distribution Networks according to the characteristics of power distribution network, and such as Newton method improves PQ solution
Coupling method, circuit impedance method and forward-backward sweep method etc..But Newton method needs to form admittance matrix, and Jacobian matrix is diagonal excellent
Gesture has not existed, and is difficult to restrain；It improves PQ decoupling method and compensation technique is introduced to the biggish route of R/X value, this algorithm is complicated, funeral
The advantage that the original calculation amount of quick decoupling is small, convergence is reliable is lost；Circuit impedance method needs complicated node and branch number, than
It is relatively time-consuming.Comparatively, being pushed forward back substitution power flow algorithm takes full advantage of the design feature of network radially, this method physics is general
Read apparent, but due to distribution system huge structure, branch is more, to the search of hierarchical structure data in calculating process, affects
The calculating speed of trend all calculates the power loss of branch, is related to the calculating directly to data matrix, accounts for when data are pushed forward every time
It is bigger and relatively time-consuming with space.
Summary of the invention
It, can be fast and efficiently to matching the object of the present invention is to provide a kind of radial pattern distribution system tidal current computing method
Power grid progress Load flow calculation, fast convergence rate, used time are few.
The technical solution adopted by the present invention are as follows:
A kind of radial pattern distribution system tidal current computing method, specifically includes following steps:
Step 1 the characteristics of according to radial pattern distribution system, establishes the circuit model of distribution system, and distribution transformer is low
It presses side as potential node, and is used as balance nodes in Load flow calculation, equivalence is voltage magnitude and phase angle is constant known quantity
And assume that three-phase voltage is symmetrical；All loads are equivalent to branch endpoint node concentrated load on feeder line, and negative for invariable power, that is, PQ
Lotus, referred to as node load；Feeder line branch uses lumped parameter model, the characteristics of according to radial pattern distribution system, therefore power distribution system
Only one potential node in system, is balance nodes, remaining node is all load bus, and the trend in distribution system is from potential
Node, which sets out, flows to load bus.Each branch is connected with two nodes in system in distribution system, and trend flows out node
Referred to as branch beginning node, trend flow into node and are known as branch endpoint node, each node is at most connected with two branches.
Step 2, according to the circuit model of distribution system, in circuit model topological structure branch and node compile
Number, wherein for branch from the branch that is connected with potential node, number consecutively 1,2,3,4 ... b, b are the branch of distribution system
Number；Potential node number is 0, and the load bus number consecutively of corresponding branch end is 1,2,3,4 ..., and n-1, n are distribution system
Number of nodes；Then have b=n-1 according to the characteristics of radial pattern distribution system, in order to down-stream it is easy for the sake of, when number, makes branch
It numbers identical as connected endpoint node number；
And the topological structure and design parameter of circuit model are stated by primary data matrix D S, primary data matrix D S
For b row, the matrixes of 5 column, wherein the i-th behavior: DS (i :)=[x, NS (x), NR (x), Z (x), S (NR (x))]；I=1,
2,......,b；X is i-th branch number, and NS (x) is the number of i-th branch beginning node, and NR (x) is i-th branch end
The number of end node, then NR (x)=i；Z (x) is i-th branch impedance, and S (NR (x)) is the section that branch endpoint node is NR (x)
The complex power of point load；
Step 3 constructs the node branch incidence matrix NB of distribution system, node branch by primary data matrix D S
Incidence matrix NB is n row, the matrix of b column, jth row the i-th column element NB (j, i) are as follows:
Wherein, j=1,2 ..., n, i=1,2 ..., b；
Step 4 constructs the path flowed through by the trend of potential node to load bus according to node branch incidence matrix NB
Matrix P, path matrix P are (n-1) row, the matrix of b column, jth row the i-th column element P (j, i) are as follows:
Wherein, j=1,2 ..., n-1, i=1,2 ..., b；
It includes the following steps that trend from potential node 0 to load bus j, which flows through the judgement of branch:
1) (j+1) row element NB (j :) in NB matrix is traversed, if NB (j, y)=1, it may be determined that branch y endpoint node is
j；
2) y column element NB (:, y) in NB matrix is traversed, if NB (z, y)=- 1, it may be determined that branch y beginning node is (z-
1)；
3) z row element NB (z :) in NB matrix is traversed, if NB (z, t)=1, it may be determined that branch t endpoint node is (z-
1)；
4) above step 2,3 is repeated, until branch beginning node is potential node 0；
It is then potential node 0 ... branch t, node (z-1), branch through path from potential node 0 to the trend of load bus j
Road y, load bus j；
Step 5, according to path matrix P and primary data matrix D S, the trend for constructing potential node to load node is flowed through
Branch impedance matrix ZP and node load complex power matrix SB；Branch impedance matrix ZP is that (n-1) goes, the matrix of b column, the
J row the i-th column element ZP (j, i) are as follows:
ZP (j, i)=P (j, i) * Z (i)；The complex power matrix SB of node load is b row, the matrix of (n-1) column, the i-th row
Jth column element SB (i, j) are as follows:
SB (i, j)=P (i, j) * S (j)；
Wherein, Z (i)=DS (i, 4) is the impedance of i-th branch；S (j)=DS (i, 5) is that the node load of j node is multiple
Power, wherein branch i endpoint node is j；J=1,2 ..., n-1, i=1,2 ..., b；
Step 6 calculates the trend of each the load bus voltage and each branch of radial pattern distribution system, specifically includes:
1) Initial Voltage Value for setting each load bus is followed successively by V1, V2... Vj... Vn-1, wherein VjIndicate load section
The Initial Voltage Value of point j, structure node voltage matrix V, node voltage matrix V are 1 row, the matrix of (n-1) column, jth column element:
V (1, j)=Vj
Wherein, j=1,2 ..., n-1；
2) node voltage conjugate matrices VC may further be obtained, node voltage conjugate matrices are 1 row, the matrix of n-1 column,
Jth column element:
VC (1, j)=V*(1,j)
The complex power conjugate matrices of node load are SC, and the complex power conjugate matrices SC of node load is b row, (n-1) column
Matrix, the i-th row jth column element are as follows:
SC (i, j)=SB*(i,j)
Wherein, i=1,2 ..., b；J=1 ..., n-1；
3) it can be calculated current matrix LC, current matrix LC is b row, the matrix of (n-1) column, the i-th row jth column element are as follows:
LC (i, j)=SC (i, j)/VC (1, j)
Branch current matrix F C, branch current matrix F C are b row, the matrix of 1 column, the i-th row element are as follows:
Wherein, i=1,2 ..., b；J=1 ..., n-1；
4) pressure drop matrix D, overall presure drop matrix M and voltage matrix T be can be calculated；Pressure drop matrix D, overall presure drop matrix M and electricity
Pressure matrix T is n-1 row, the matrix of b column；
D (j, i)=ZP (j, i) * FC (i, 1)
T (j, i)=V0-M(j,i)
It can be calculated new node voltage matrix Vnew, new node voltage matrix Vnew is 1 row, the matrix of (n-1) column, jth
Column element: Vnew (1, j)=min (T (j, 1), T (j, 2) ... T (j, (n-1)))
Wherein, V0For voltage measuring value at radial pattern distribution system low-pressure side transformer outlet, i=1,2 ..., b；j
=1 ..., n-1；L=1,2 ..., i；
5) determine the condition of convergence, select the error precision of system for δ, determine | Vnew (1, j)-V (1, j) |≤δ whether at
Vertical, if invalid, V (1, j)=Vnew (1, j) simultaneously returns to step 2)；
Otherwise, the trend of radial pattern distribution system each load bus voltage and each branch can be obtained:
The voltage of node j is Uj=V (1, j)
I-th branch current is Ili=FC (i, 1)
I-th Branch Power Flow be
Wherein, i=1,2 ..., b；J=1 ..., n-1.
The present invention establishes the circuit model of distribution system according to the characteristics of distribution system, in the circuit model of distribution system
Node and branch be numbered, and using matrix form come the primary data of indication circuit, especially circuit topological structure
Digitization description, makes the logical relation clear and definite that figure connects in circuit model, and the introducing of incidence matrix is realized to path
The Programmed Design of search process, and primary data matrix is used to input as unique data, keep entire programming directly simple；
Each Branch Power Flow and each load bus voltage of power distribution network, thus journey are solved in calculating process using the algebraic operation of matrix element
The sequence calculating used time is less, and convergence rate is very fast.
Detailed description of the invention
In order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, to embodiment or will show below
There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this
Some embodiments of invention for those of ordinary skill in the art without creative efforts, can be with
It obtains other drawings based on these drawings.
Fig. 1 is flow chart of the invention；
Fig. 2 is the circuit model schematic of radial pattern distribution system described in the embodiment of the present invention；
Fig. 3 is the circuit model branch node number schematic diagram of radial pattern distribution system described in the embodiment of the present invention；
Fig. 4 is flowing through branch to the trend of load bus 9 slave potential node and sentencing for distribution system described in the embodiment of the present invention
Determine schematic diagram.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, those of ordinary skill in the art are obtained every other under that premise of not paying creative labor
Embodiment shall fall within the protection scope of the present invention.
It is as shown in Figure 1,2 and 3, of the invention including the following steps:
Step 1 the characteristics of according to radial pattern distribution system, establishes the circuit model of distribution system, and distribution transformer is low
It presses side as potential node, and is used as balance nodes in Load flow calculation, equivalence is voltage magnitude and phase angle is constant known quantity
And assume that three-phase voltage is symmetrical；All loads are equivalent to branch endpoint node concentrated load on feeder line, and negative for invariable power, that is, PQ
Lotus, referred to as node load；Feeder line branch uses lumped parameter model, the characteristics of according to radial pattern distribution system, therefore power distribution system
Only one potential node in system, is balance nodes, remaining node is all load bus, and the trend in distribution system is from potential
Node, which sets out, flows to load bus.Each branch is connected with two nodes in system in distribution system, and trend flows out node
Referred to as branch beginning node, trend flow into node and are known as branch endpoint node, each node is at most connected with two branches.
Step 2, according to the circuit model of distribution system, in circuit model topological structure branch and node compile
Number, wherein for branch from the branch that is connected with potential node, number consecutively 1,2,3,4 ... b, b are the branch of distribution system
Number；Potential node number is 0, and the load bus number consecutively of corresponding branch end is 1,2,3,4 ..., and n-1, n are distribution system
Number of nodes；In order to distinguish the two in figure, the differentiation for adding bracket is carried out, but the two is identical correspondence in actual use
Number.Then have b=n-1 according to the characteristics of radial pattern distribution system, in order to down-stream it is easy for the sake of, when number, makes branch number
It is identical as connected endpoint node number；In order to both distinguish in figure, carry out plus the differentiation of bracket, but in actual use
The two is identical reference numeral.
And the topological structure and design parameter of circuit model are stated by primary data matrix D S, primary data matrix D S
For b row, the matrixes of 5 column, wherein the i-th behavior: DS (i :)=[x, NS (x), NR (x), Z (x), S (NR (x))]；I=1,
2,......,b；X is i-th branch number, and NS (x) is the number of i-th branch beginning node, and NR (x) is i-th branch end
The number of end node, then NR (x)=i；Z (x) is i-th branch impedance, and S (NR (x)) is the section that branch endpoint node is NR (x)
The complex power of point load；
Step 3 constructs the node branch incidence matrix NB of distribution system, node branch by primary data matrix D S
Incidence matrix NB is n row, the matrix of b column, jth row the i-th column element NB (j, i) are as follows:
Wherein, j=1,2 ..., n, i=1,2 ..., b；
Step 4 constructs the path flowed through by the trend of potential node to load bus according to node branch incidence matrix NB
Matrix P, path matrix P are (n-1) row, the matrix of b column, jth row the i-th column element P (j, i) are as follows:
Wherein, j=1,2 ..., n-1, i=1,2 ..., b；
It includes the following steps that trend from potential node 0 to load bus j, which flows through the judgement of branch:
5) (j+1) row element NB (j :) in NB matrix is traversed, if NB (j, y)=1, it may be determined that branch y endpoint node is
j；
6) y column element NB (:, y) in NB matrix is traversed, if NB (z, y)=- 1, it may be determined that branch y beginning node is (z-
1)；
7) z row element NB (z :) in NB matrix is traversed, if NB (z, t)=1, it may be determined that branch t endpoint node is (z-
1)；
8) above step 2,3 is repeated, until branch beginning node is potential node 0；
It is then potential node 0 ... branch t, node (z-1), branch through path from potential node 0 to the trend of load bus j
Road y, load bus j；
Step 5, according to path matrix P and primary data matrix D S, the trend for constructing potential node to load node is flowed through
Branch impedance matrix ZP and node load complex power matrix SB；Branch impedance matrix ZP is that (n-1) goes, the matrix of b column, the
J row the i-th column element ZP (j, i) are as follows:
ZP (j, i)=P (j, i) * Z (i)；The complex power matrix SB of node load is b row, the matrix of (n-1) column, the i-th row
Jth column element SB (i, j) are as follows:
SB (i, j)=P (i, j) * S (j)；
Wherein, Z (i)=DS (i, 4) is the impedance of i-th branch；S (j)=DS (i, 5) is that the node load of j node is multiple
Power, wherein branch i endpoint node is j；J=1,2 ..., n-1, i=1,2 ..., b；
Step 6 calculates the trend of each the load bus voltage and each branch of radial pattern distribution system, specifically includes:
1) Initial Voltage Value for setting each load bus is followed successively by V1, V2... Vj... Vn-1, wherein VjIndicate load section
The Initial Voltage Value of point j, structure node voltage matrix V, node voltage matrix V are 1 row, the matrix of (n-1) column, jth column element:
V (1, j)=Vj
Wherein, j=1,2 ..., n-1；
2) node voltage conjugate matrices VC may further be obtained, node voltage conjugate matrices are 1 row, the matrix of n-1 column,
Jth column element:
VC (1, j)=V* (1, j)
The complex power conjugate matrices of node load are SC, and the complex power conjugate matrices SC of node load is b row, (n-1) column
Matrix, the i-th row jth column element are as follows:
SC (i, j)=SB* (i, j)
Wherein, i=1,2 ..., b；J=1 ..., n-1；
3) it can be calculated current matrix LC, current matrix LC is b row, the matrix of (n-1) column, the i-th row jth column element are as follows:
LC (i, j)=SC (i, j)/VC (1, j)
Branch current matrix F C, branch current matrix F C are b row, the matrix of 1 column, the i-th row element are as follows:
Wherein, i=1,2 ..., b；J=1 ..., n-1；
4) pressure drop matrix D, overall presure drop matrix M and voltage matrix T be can be calculated；Pressure drop matrix D, overall presure drop matrix M and electricity
Pressure matrix T is n-1 row, the matrix of b column；
D (j, i)=ZP (j, i) * FC (i, 1)
T (j, i)=V0-M(j,i)
It can be calculated new node voltage matrix Vnew, new node voltage matrix Vnew is 1 row, the matrix of (n-1) column, jth
Column element: Vnew (1, j)=min (T (j, 1), T (j, 2) ... T (j, (n-1)))
Wherein, V0For voltage measuring value at radial pattern distribution system low-pressure side transformer outlet, i=1,2 ..., b；j
=1 ..., n-1；L=1,2 ..., i；
5) determine the condition of convergence, select the error precision of system for δ, determine | Vnew (1, j)-V (1, j) |≤δ whether at
Vertical, if invalid, V (1, j)=Vnew (1, j) simultaneously returns to step 2)；
Otherwise, the trend of radial pattern distribution system each load bus voltage and each branch can be obtained:
The voltage of node j is Uj=V (1, j)
I-th branch current is Ili=FC (i, 1)
I-th Branch Power Flow be
Wherein, i=1,2 ..., b；J=1 ..., n-1.
In order to further illustrate the present invention, it is illustrated below with specific embodiment, please refers to attached drawing.
Referring to Fig. 1, a kind of radial pattern distribution system tidal current computing method specifically includes following steps:
Step 1, the circuit model for providing the radial pattern distribution system that one simplifies here can be found in Fig. 2.
Using distribution transformer low-pressure side as potential node 1, equivalence is voltage magnitude and phase angle is constant known quantity and vacation
It is symmetrical to determine three-phase voltage；All loads are equivalent to node load 2 on feeder line, and are invariable power, that is, PQ load, and feeder line branch uses
Lumped parameter model, therefore only one potential node 1 in power distribution network, remaining node are all load bus 3；It is each in distribution system
Branch 4 is connected with two nodes in system, and trend outflow node is known as branch beginning node 5, and trend flows into node and is known as
Branch endpoint node 6；
Step 2, in the topological structure of radial pattern distribution system circuit model as shown in Figure 2 branch and node into
Row number is as shown in Figure 3, wherein branch is from the branch that is connected with potential node, number consecutively 1,2 ... 9, distribution system
Circuitry number b be 9；Potential node number is 0, and corresponding branch endpoint node number consecutively is 1,3,4 ... 9, distribution system
Number of nodes n is 10, including a potential node, 9 load bus；For the sake of down-stream simplicity, make to prop up here when number
Road number is identical as connected endpoint node number；
And stated by primary data matrix D S circuit model topological relation and design parameter it is as follows:
Wherein, primary data matrix D S is 9 rows, and the matrix of 5 column, Z (i) is i-th branch impedance, and S (j) is endpoint node
For the complex power of the node load of j, wherein the endpoint node of branch i is j, then i=j；I=1,2 ..., 9, j=1,
2,......,9；
Step 3, the node branch incidence matrix NB that can construct distribution system by primary data matrix D S are as follows:
Step 4, the road flowed through according to node branch incidence matrix NB and construction by the trend of potential node to load node
Drive matrix P is as follows:
Referring to fig. 4, for flowing through the decision process of branch to 9 trend of node from potential node 0, specific decision process institute
Show, including the following steps:
1) the 10th row element NB (10 :) NB (10,9)=1 in NB matrix is traversed, it may be determined that 9 endpoint node of branch is 9；
2) the 9th column element NB (:, 9), NB (4,9)=- 1 in NB matrix are traversed, it may be determined that 9 beginning node of branch is 3；
3) the 4th row element NB (4 :), NB (4,3)=1 in NB matrix are traversed, it may be determined that 3 endpoint node of branch is 3；
4) the 3rd column element NB (:, 3), NB (3,3)=- 1 in NB matrix are traversed, it may be determined that 3 beginning node of branch is 2；
5) the 2nd row element NB (3 :), NB (3,2)=1 in NB matrix are traversed, it may be determined that 2 endpoint node of branch is 2；
6) the 2nd column element NB (:, 2), NB (2,2)=- 1 in NB matrix are traversed, it may be determined that 2 beginning node of branch is 1；
7) the 2nd row element NB (2 :), NB (2,1)=1 in NB matrix are traversed, it may be determined that 1 endpoint node of branch is 1；
8) the 1st column element NB (:, 1), NB (1,1)=- 1 in NB matrix are traversed, it may be determined that 1 beginning node of branch is 0；
Then from potential node 0 to the trend of node j through path be potential node 0, branch 1, node 1, branch 2, node 2,
Branch 3, node 3, branch 9, node 9.
Step 5, according to path matrix P and primary data matrix D S, the trend electric current of construction potential node to load node
The load complex power matrix SB of the branch impedance matrix ZP and node flowed through is as follows:
Step 6 calculates each branch current and each load bus voltage of power distribution network, specifically includes
1) other each load bus Initial Voltage Values are arbitrarily set and are followed successively by V1, V2…Vj…Vn-1, wherein VjIndicate the
J node voltage initial value, structure node voltage matrix V.
2) node voltage conjugate matrices VC may further be obtained, node voltage conjugate matrices are 1 row, the matrixes of 9 column,
Middle jth column element:
VC (1, j)=V*(1,j)
Node load complex power conjugate matrices be SC, node load complex power conjugate matrices SC be 9 rows, 9 column matrixes,
In the i-th row jth column element are as follows:
SC (i, j)=SB*(i,j)
Wherein, i=1,2 ..., 9；J=1 ..., 9；
3) current matrix LC, branch current matrix F C be can be calculated,
LC (i, j)=SC (i, j)/VC (1, j)
Wherein, i=1,2 ..., 9；J=1 ..., 9；
It is convenient to write, it might as well set here:
4) pressure drop matrix D, overall presure drop matrix M and voltage matrix T be can be calculated；Pressure drop matrix D, overall presure drop matrix M and electricity
Pressure matrix T is 9 rows, the matrix of 9 column；
D (j, i)=ZP (j, i) * FC (i, 1)
T (j, i)=V0-M(j,i)
It can be calculated new node voltage matrix Vnew, new node voltage matrix Vnew is 1 row, the matrix of 9 column, jth column member
Element: Vnew (1, j)=min (T (j, 1), T (j, 2) ... T (j, (n-1)))
Wherein,
M (1)=Z (1) * Il (1)
M (2)=Z (1) * Il (1)+Z (2) * Il (2)
M (3)=Z (1) * Il (1)+Z (2) * Il (2)+Z (3) * Il (3)
M (4)=Z (1) * Il (1)+Z (2) * Il (2)+Z (3) * Il (3)+Z (4) * Il (4)
M (5)=Z (1) * Il (1)+Z (5) * Il (5)
M (6)=Z (1) * Il (1)+Z (5) * Il (5)+Z (6) * Il (6)
M (7)=Z (1) * Il (1)+Z (7) * Il (7)
M (8)=Z (1) * Il (1)+Z (2) * Il (2)+Z (8) * Il (8)
M (9)=Z (1) * Il (1)+Z (2) * Il (2)+Z (3) * Il (3)+Z (9) * Il (9)
Wherein, V0For voltage measuring value at distribution system low-pressure side transformer outlet, i=1,2 ..., 9；J=
1,......,9；
L=1,2 ..., i；
Vnew(1,j)
=[V0-M(1) V0-M(2) V0-M(3) V0-M(4) V0-M(5) V0-M(6) V0-M(7) V0-M(8) V0-M
(9)]
5) determine the condition of convergence, select the error precision of system for δ, determine | Vnew (1, j)-V (1, j) |≤δ whether at
It is vertical,
If invalid, V (1, j)=Vnew (1, j) simultaneously returns to step 2)；
Otherwise, the trend of each the load bus voltage and each branch of radial pattern distribution system can be obtained:
The voltage of node j is Uj=V (1, j)
I-th branch current is Ili=FC (i, 1)
I-th branch complex power be
Wherein, i=1,2 ..., b；J=1 ..., n-1.
The present invention provides a kind of tidal current computing method of radial pattern distribution system, for matching for the direction of distribution net platform region one
Electric system carry out Load flow calculation, first according to radial pattern distribution system the characteristics of, establish the circuit model of radial pattern distribution system,
Its topological structure is reasonably numbered, and utilizes the initial data of matrix form storage distribution system, including topological structure
Relationship and electrical network parameter, and inputted as unique data, establish the node branch incidence matrix of distribution system and by electricity
The path matrix that the trend of gesture node to load node flows through further constructs the complex power of branch impedance matrix and node load
Matrix；It is last directly to solve each load bus voltage of distribution system and answering for each branch using the algebraic operation of matrix element
Power.
The foregoing is merely illustrative of the preferred embodiments of the present invention, is not intended to limit the invention, all in essence of the invention
Within mind and principle, any modification, equivalent replacement, improvement and so on be should all be included in the protection scope of the present invention.
Claims (6)
1. a kind of radial pattern distribution system tidal current computing method, it is characterised in that: specifically include following steps:
Step 1 the characteristics of according to radial pattern distribution system, establishes the circuit model of radial pattern distribution system: by distribution transformer
Low-pressure side is used as balance nodes as potential node in Load flow calculation, and equivalence is voltage magnitude and phase angle is constant known
It measures and assumes that three-phase voltage is symmetrical；All loads are equivalent to branch endpoint node concentrated load, referred to as node load on feeder line,
Feeder line branch uses lumped parameter model；Only one potential node in radial pattern distribution system, remaining node are all load sections
Point, the trend in distribution system flow to load bus from potential node；In distribution system in each branch and system
Two nodes be connected, trend outflow node is known as branch beginning node, and trend flows into node and is known as branch endpoint node, each
A node is at most connected with two branches；
Step 2, according to the circuit model of radial pattern distribution system, branch and node to topological structure in circuit model are carried out
Number, for branch from the branch that is connected with potential node, number consecutively 1,2,3 ... b, b are the circuitry number of distribution system；Electricity
Gesture node serial number is 0, and corresponding branch endpoint node number consecutively is 1,3,4 ..., and n-1, n are the number of nodes of distribution system；According to
, then there is b=n-1 in the characteristics of radial pattern distribution system, and it is identical as the endpoint node number that branch is connected to number seasonal branch number；
The topological structure and design parameter of circuit model are stated by primary data matrix D S；
Step 3 constructs the node branch incidence matrix NB of radial pattern distribution system according to primary data matrix D S；
Step 4 indicates the path flowed through by the trend of potential node to load bus according to node branch incidence matrix NB construction
Matrix P；
Step 5 constructs the branch that the trend of potential node to load bus flows through according to path matrix P and primary data matrix D S
The complex power matrix SB of road impedance matrix ZP and node load；
Step 6 calculates each load bus voltage and each Branch Power Flow in radial pattern distribution system.
2. radial pattern distribution system tidal current computing method according to claim 1, it is characterised in that: the step one
In, primary data matrix D S is b row, the matrixes of 5 column, wherein the i-th behavior: DS (i :)=[x, NS (x), NR (x), Z (x), S
(NR(x；))] i=1,2 ..., b；X is i-th branch number, and NS (x) is the number of i-th branch beginning node, NR
It (x) is the number of i-th branch endpoint node, therefore NR (x)=i, Z (x) they are i-th branch impedance, S (NR (x)) is end segment
Point is the node load complex power of NR (x).
3. radial pattern distribution system tidal current computing method according to claim 1, it is characterised in that: in the step three
Node branch incidence matrix NB is n row, the matrix of b column, jth row the i-th column element NB (j, i) are as follows:Wherein, j=1,2 ..., n, i=1,2 ..., b.
4. radial pattern distribution system tidal current computing method according to claim 1, it is characterised in that: in the step four
Path matrix P is (n-1) row, the matrix of b column, jth row the i-th column element P (j, i) are as follows:Wherein, j=1,2 ..., n-1, i=1,
2,......,b；Wherein, flowing through the judgement in path to the trend of load bus j from potential node 0 includes the following steps:
(j+1) row element NB (j :) is traversed in NB matrix if NB (j, y)=1, it may be determined that branch y endpoint node is j；
Y column element NB (:, y) in NB matrix is traversed, if NB (z, y)=- 1, it may be determined that branch y beginning node is (z-1)；
Z row element NB (z :) in NB matrix is traversed, if NB (z, t)=1, it may be determined that branch t endpoint node is (z-1)；
Above step is repeated, until determining that branch beginning node is potential node 0.
5. radial pattern distribution system tidal current computing method according to claim 1, it is characterised in that: the step five
In, branch impedance matrix ZP is (n-1) row, the matrix of b column, jth row the i-th column element ZP (j, i) are as follows:
ZP (j, i)=P (j, i) * Z (i)；The complex power matrix SB of node load is b row, the matrix of (n-1) column, the i-th row jth column
Element S B (i, j)=P (i, j) * S (j)；Wherein, Z (i)=DS (i, 4) is the impedance of i-th branch；S (j)=DS (i, 5) is j
The node load complex power of node, branch i endpoint node are j；J=1,2 ..., n-1, i=1,2 ..., b.
6. radial pattern distribution system tidal current computing method according to claim 1, it is characterised in that: step six tool
Body includes the following steps: to specifically include:
6-1, the Initial Voltage Value for setting each load bus are followed successively by V1, V2... Vj... Vn- 1, wherein VjIndicate load bus
The Initial Voltage Value of j, structure node voltage matrix V, node voltage matrix V are 1 row, and the matrix of (n-1) column, wherein jth column are first
Element:
V (1, j)=Vj
Wherein, j=1,2 ..., n-1；
6-2, may further obtain node voltage conjugate matrices VC, and node voltage conjugate matrices are 1 row, the matrix of n-1 column,
Middle jth column element:
VC (1, j)=V*(1,j)
The complex power conjugate matrices of node load are SC, and the complex power conjugate matrices SC of node load is b row, the square of (n-1) column
Battle array, wherein the i-th row jth column element are as follows:
SC (i, j)=SB*(i,j)
Wherein, i=1,2 ..., b；J=1 ..., n-1；
6-3, can be calculated load current matrix LC, and branch current matrix F C, load current matrix LC are b row, the square of (n-1) column
Battle array, wherein the i-th row jth column element are as follows:
LC (i, j)=SC (i, j)/VC (1, j)
Branch current matrix F C is b row, the matrix of 1 column, wherein the i-th row element are as follows:
Wherein, i=1,2 ..., b；J=1 ..., n-1；
6-4 can be calculated pressure drop matrix D, overall presure drop matrix M and voltage matrix T；Pressure drop matrix D, overall presure drop matrix M and voltage
Matrix T is n-1 row, the matrix of b column；
D (j, i)=ZP (j, i) * FC (i, 1)
T (j, i)=V0-M(j,i)
It can be calculated new node voltage matrix Vnew, new node voltage matrix Vnew is 1 row, the matrix of (n-1) column, jth column member
Element: Vnew (1, j)=min (T (j, 1), T (j, 2) ... T (j, (n-1)))
Wherein, V0For voltage measuring value at radial pattern distribution system low-pressure side transformer outlet, i=1,2 ..., b；J=
1,......,n-1；L=1,2 ..., i；
6-5 determines the condition of convergence, selects the error precision of system for δ, determines | Vnew (1, j)-V (1, j) | whether≤δ is true,
If invalid, V (1, j)=Vnew (1, j) simultaneously returns to step 6-2)；
Otherwise, each node voltage of radial pattern distribution system and each Branch Power Flow can be obtained:
The voltage of node j is Uj=V (1, j)
I-th branch current is Ili=FC (i, 1)
I-th Branch Power Flow be
Wherein, Z (i)=DS (i, 4) is the impedance of i-th branch；I=1,2 ..., b；J=1 ..., n-1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910424449.3A CN109980650B (en) | 2019-05-21 | 2019-05-21 | Load flow calculation method of radiation type power distribution system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910424449.3A CN109980650B (en) | 2019-05-21 | 2019-05-21 | Load flow calculation method of radiation type power distribution system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109980650A true CN109980650A (en) | 2019-07-05 |
| CN109980650B CN109980650B (en) | 2022-11-15 |
Family
ID=67073880
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910424449.3A Active CN109980650B (en) | 2019-05-21 | 2019-05-21 | Load flow calculation method of radiation type power distribution system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109980650B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116029251A (en) * | 2023-03-23 | 2023-04-28 | 青岛青软晶尊微电子科技有限公司 | Circuit wiring optimization method and device based on circuit performance equalization |
| CN116581754A (en) * | 2023-07-12 | 2023-08-11 | 国网江苏省电力有限公司电力科学研究院 | Method, device, equipment and medium for calculating real-time power flow of power distribution network |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH08214458A (en) * | 1995-02-02 | 1996-08-20 | Fuji Electric Co Ltd | Power flow calculation for radial distribution system |
| CN103150606A (en) * | 2013-01-22 | 2013-06-12 | 中国电力科学研究院 | Optimal power flow optimization method of distributed power supplies |
| CN103795057A (en) * | 2014-01-23 | 2014-05-14 | 华北电力大学（保定） | Power distribution network topology serial number generating method based on searching in power distribution network load flow calculation |
| CN103928925A (en) * | 2014-04-17 | 2014-07-16 | 国家电网公司 | Power distribution network load flow calculation method based on forward-backward sweep |
| CN106786602A (en) * | 2017-02-04 | 2017-05-31 | 中南大学 | A kind of distribution power system load flow calculation method |
| CN108599167A (en) * | 2018-01-15 | 2018-09-28 | 国网吉林省电力有限公司电力科学研究院 | A kind of linearisation tidal current computing method of radial distribution networks |
-
2019
- 2019-05-21 CN CN201910424449.3A patent/CN109980650B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH08214458A (en) * | 1995-02-02 | 1996-08-20 | Fuji Electric Co Ltd | Power flow calculation for radial distribution system |
| CN103150606A (en) * | 2013-01-22 | 2013-06-12 | 中国电力科学研究院 | Optimal power flow optimization method of distributed power supplies |
| CN103795057A (en) * | 2014-01-23 | 2014-05-14 | 华北电力大学（保定） | Power distribution network topology serial number generating method based on searching in power distribution network load flow calculation |
| CN103928925A (en) * | 2014-04-17 | 2014-07-16 | 国家电网公司 | Power distribution network load flow calculation method based on forward-backward sweep |
| CN106786602A (en) * | 2017-02-04 | 2017-05-31 | 中南大学 | A kind of distribution power system load flow calculation method |
| CN108599167A (en) * | 2018-01-15 | 2018-09-28 | 国网吉林省电力有限公司电力科学研究院 | A kind of linearisation tidal current computing method of radial distribution networks |
Non-Patent Citations (2)
| Title |
|---|
| 付敏 等: "辐射状配电网潮流的改进前推回代法", 《哈尔滨理工大学学报》 * |
| 高军彦 等: "辐射型配电网潮流算法综述", 《中国电力教育》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116029251A (en) * | 2023-03-23 | 2023-04-28 | 青岛青软晶尊微电子科技有限公司 | Circuit wiring optimization method and device based on circuit performance equalization |
| CN116581754A (en) * | 2023-07-12 | 2023-08-11 | 国网江苏省电力有限公司电力科学研究院 | Method, device, equipment and medium for calculating real-time power flow of power distribution network |
| CN116581754B (en) * | 2023-07-12 | 2023-09-26 | 国网江苏省电力有限公司电力科学研究院 | Method, device, equipment and medium for calculating real-time power flow of power distribution network |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109980650B (en) | 2022-11-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103310065B (en) | Meter and the intelligent network distribution reconstructing method of distributed power generation and energy-storage units | |
| CN110707706B (en) | Power transmission network planning method and system based on line power flow distribution | |
| CN107196315A (en) | The extendable power-less optimized controlling method of the power distribution network containing light-preserved system | |
| CN106532710A (en) | Microgrid power flow optimization method considering voltage stability constraint | |
| CN108631278B (en) | The Optimal Configuration Method of breaker and fault current limiter in a kind of looped network formula direct-current micro-grid | |
| CN109980651A (en) | A kind of looped network type distribution system tidal current computing method | |
| CN109980650A (en) | A kind of radial pattern distribution system tidal current computing method | |
| CN109494719A (en) | A kind of mesolow mixing power distribution network stratification impedance analysis method | |
| CN107546773A (en) | A kind of more micro-capacitance sensor dynamic networking methods in region based on graph theory | |
| CN108462194A (en) | A kind of wide area optimization method for low-voltage network three-phase load unbalance | |
| CN109995059B (en) | Low-voltage power distribution system three-phase load optimization balance configuration method based on genetic algorithm | |
| CN111079972A (en) | Method, device and medium for planning reliability of active power distribution network | |
| CN108830451A (en) | A kind of the convergence potential evaluation method and system of user side distributed energy storage | |
| CN105914752A (en) | Pilot node selection method based on clustering by fast search and density peaks | |
| CN110943465B (en) | Energy storage system site selection and volume fixing optimization method | |
| CN109829599A (en) | The assemblage classification method and device of power distribution network based on high proportion renewable energy | |
| CN103515964A (en) | Reactive compensation control method and reactive compensation control device | |
| CN109995039A (en) | A kind of three-phase imbalance low-voltage distribution system tidal current computing method | |
| CN105896547A (en) | Large power network graded voltage control method under wind power access | |
| CN104009469B (en) | A kind of low-voltage distribution network platform district network topology isopleth map collecting method | |
| CN105449672A (en) | Method for estimating total supply capability of 220KV loop ring network | |
| CN113300477B (en) | Optimization method for energy storage configuration of central urban power grid | |
| Talwariya et al. | Optimization of distribution networks with integration of distributed generators using cooperative Game Theory | |
| CN105048475B (en) | The points distributing method of multi-infeed DC receiving-end system dynamic reactive compensation device | |
| CN106055730A (en) | Method for electric power system carbon-flow simulation analysis adapted to large-scale power grid |
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 |