CN106208082B - A kind of bilingual coupling tidal current computing method of power distribution network - Google Patents
A kind of bilingual coupling tidal current computing method of power distribution network Download PDFInfo
- Publication number
- CN106208082B CN106208082B CN201610839531.9A CN201610839531A CN106208082B CN 106208082 B CN106208082 B CN 106208082B CN 201610839531 A CN201610839531 A CN 201610839531A CN 106208082 B CN106208082 B CN 106208082B
- Authority
- CN
- China
- Prior art keywords
- phase
- power
- formula
- decoupling
- distribution network
- 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
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 bilingual coupling tidal current computing methods of power distribution network.Efficient computational problem of the present invention to solve three-phase distribution net trend is combined to distributed generation resource, line parameter circuit value asymmetry, the adaptability of high R/X ratio and weakly loops.The present invention directly carries out alternate decoupling to distribution line three-phase using alternate decoupling compensation model first；Then reduce the influence of the R/X ratio of each phase using complex field mark change technology, then introduce fast decoupled power flow algorithm and carry out decomposition operation in phase.Two steps decoupling in this method reduces Jacobian matrix dimension and remains unchanged it in calculating process, enormously simplifies computational complexity.Unsymmetrical distribution net decoupling Load flow calculation problem can be effectively treated in the present invention, have the characteristics that calculating speed is fast and adaptable, there is good engineering application value.
Description
Technical field
The present invention relates to a kind of bilingual coupling tidal current computing methods of power distribution network, belong to Power System Analysis field.
Technical background
Distribution power system load flow calculation is to carry out the basis of distribution system analysis.In recent years, with a large amount of accesses of distributed generation resource
And the fast development of active distribution network technology, conventional electrical distribution net from " Dan Yuan " mode to " multi-source " Mode change, system
Asymmetry is increasingly apparent, and network structure is more complicated, thus traditional distribution power flow analysis method be difficult to adapt to it is current
Power distribution network.In order to further increase the computational efficiency of current power distribution network trend, enhance the adaptability of algorithm, it is necessary to find one kind
Efficient distribution power flow fast algorithm.
As the innovatory algorithm of classical Newton method, fast decoupled power flow algorithm is fast by its convergence, saves the advantage of memory
It is applied widely in power transmission network Load flow calculation, but since the characteristics of route high R/X ratio, fast decoupled is calculated in power distribution network
Method not can be used directly in power distribution network.The complex field mark change technology risen in recent years efficiently solves high R/ in distribution line
The problem of X ratio, and then fast decoupled power flow algorithm is successfully introduced into distribution system Load flow calculation, but this method is not
Consider the three-phase imbalance problem of power distribution network.Therefore, for power supply, route, load three-phase imbalance in modern power distribution systems with
And looped network phenomenon, it is necessary to study a kind of quick, efficient, adaptable Power Flow Calculation Methods For Distribution Network.
Summary of the invention
Goal of the invention: the present invention provides a kind of bilingual coupling trend of power distribution network for the technical issues of solution needed for the prior art
Calculation method.
Technical solution: the present invention is a kind of bilingual coupling tidal current computing method of power distribution network, comprising the following steps:
Step 1: power distribution network three-phase data is carried out data separating, is separated into A phase dependency number by input power distribution network three-phase data
According to, B phase related data and C phase related data；
Step 2: each phase data of the ABC isolated being subjected to complex field mark respectively and is changed, and establishes respective PQ points respectively
Jacobian matrix needed for solution；
Step 201: power reference value is selected in complex fieldVoltage reference valueAre as follows:
In formula, SbaseFor apparent energy, UbaseFor benchmark voltage magnitude, e is the bottom of natural logrithm, and j is imaginary unit, φs
For corresponding power reference angle, φvFor corresponding voltage reference angle；
Then acquire corresponding impedance reference value are as follows:
In formula, ZbaseFor impedance reference value modulus value, subscript * indicates conjugation；
So, the mark of power changes form in complex field are as follows:
In formula, m=1,2,3..., t, t are injecting power number；Pm, QmIndicate m-th injection active power, reactive power；
Subscript cpu indicates that the mark under complex field is changed；Subscript pu indicates that the mark under traditional real number field is changed；
The mark of line impedance changes form in complex field are as follows:
αi=arctan (Xi/Ri) (6)
In formula, i=1,2,3 ..., l, l are circuitry number；Ri,XiRespectively indicate i-th of resistance and reactance；αiIt indicates i-th
The tangent angle of reactance and resistance；Subscript cpu indicates that the mark under complex field is changed；Subscript pu indicates the mark under traditional real number field
Change, ZiFor the impedance of branch i；
Find out from formula (5), the impedance value after the change of complex field subscript increases compared with the impedance value after the change of real number field subscript
A power reference angle φ is addeds, this be complex field fast decoupled trend used in the power distribution network one it is vital because
Element；By the experience of traditional power transmission network Load flow calculation, line parameter circuit value R < < X is the necessary condition of fast decoupled Load flow calculation,
It is that traditional power transmission network fast decoupled power flow algorithm is difficult to adapt to the major reason of power distribution network, how is described in detail below with multiple
The method that number field mark is changed, solves the problems, such as high R/X ratio in distribution line；
Step 202: by formula (5) it is found that complex field mark corresponding resistance and reactance after changing are as follows:
To make Ri,cpu<<Xi,cpu, then
That is,
By formula (9) it is found that αiIt is to be determined by the intrinsic parameter of route, therefore select a suitable power reference angle φsMake
?So that Ri,cpu<<Xi,cpu, to enhance fast decoupled power flow algorithm in the adaptation of distribution network
Property；Voltage reference angle φ is noticed by formula (1)~(9)vSelection, have no effect on, therefore take φv=0 to mitigate the complexity of problem
Degree calculates to simplify；
Step 203: complex field power reference angle φsSelection to fast decoupled power flow algorithm distribution network application have
Decisive role, from the property of tangent function:
By formula (10) it is found that as (αi+φs) level off to pi/2 left side when, tan (αi+φs) level off to+∞；As (αi+
φs) level off to pi/2 right side when, tan (αi+φs) leveling off to-∞, analysis can obtain, both approach modes can meet quickly
The requirement of decomposition computation, only as selected φsSo that tan (the α of part routei+φs) level off to-∞ when, these routes R/X
Ratio will appear the case where negative value, this does not meet the normal physical phenomenon that line parameter circuit value is shown, therefore select in complex field
Power reference angle φsWhen, it is intended merely to reduce the influence of R/X ratio as much as possible, matched with enhancing fast decoupled Load flow calculation
The applicability of power grid, there is no actual physical significances for each parameter of route after the change of complex field mark；
Power reference angle φ is chosen using the method for optimizations, initially set up following objective function:
In formula, l is circuitry number；
According to formula (11), function f (φs) indicate the α of all routesiPhase shift φsLater, to the left or to the right close to the journey of pi/2
Degree；Institute is in the hope of making objective function f (φs) the smallest φ of valuesAs our optimal solutions to be looked for, enable f (φs) to φsLocal derviation
It is zero, φ can be acquiredsOptimal solution are as follows:
In formula, l is circuitry number；
Step 3: the voltage magnitude of each phase and phase angle being set to initial value, and iteration count is set to 1；
Step 4: current compensation is carried out to A phase using decoupling compensation model, a fast decoupled is then carried out and calculates, amendment
The voltage magnitude and phase angle of A phase part；Similarly, the calculating similar with A is carried out to B phase, C phase, and corrects B phase C phase portion respectively
The voltage magnitude and phase angle divided；
Step 5: judging whether the maximum value of ABC phase correction amount absolute value reaches convergent requirement, if reaching convergent requirement
Then stop calculating, if not up to convergent requirement, continues cycling through, until reaching convergent requirement；Finally output is as a result, knot
Beam calculates.
The utility model has the advantages that the present invention is compared with prior art: alternate decouple of one kind proposed by the present invention is quickly divided with complex field
The bilingual coupling of three-phase distribution net trend may be implemented in the bilingual coupling three-phase rapid tidal current algorithm of the active distribution network that resolving Algorithm combines
It efficiently calculates, combines to distributed generation resource, line parameter circuit value asymmetry, the adaptability of high R/X ratio and weakly loops.This method
Alternate decoupling is directly carried out to distribution line three-phase using alternate decoupling compensation model first；Then change skill using complex field mark
Art reduces the influence of the R/X ratio of each phase, then introduces fast decoupled power flow algorithm and carries out decomposition operation in phase.In this method
The decoupling of two steps reduces Jacobian matrix dimension and remains unchanged it in calculating process, and it is complicated to enormously simplify calculating
Property.
Detailed description of the invention
Fig. 1 is the bilingual coupling Load flow calculation flow chart of power distribution network of the present invention；
Fig. 2 is distribution line three-phase π type Type Equivalent Circuit Model schematic diagram of the present invention；
Fig. 3 is A phase decoupling compensation equivalent-circuit model schematic diagram of the present invention；
Fig. 4 is the voltage magnitude schematic diagram of IEEE13 node system of the present invention；
Fig. 5 is the relational graph of IEEE34 node system four kinds of method the number of iterations and convergence precision；
Fig. 6 is the average calculation times schematic diagram of each method.
Specific embodiment
Implementation of the invention is described further below in conjunction with attached drawing and example, but implementation of the invention and comprising unlimited
In this.
A kind of bilingual coupling tidal current computing method of power distribution network, comprising the following steps:
Step 1: power distribution network three-phase data is carried out data separating, is separated into A phase dependency number by input power distribution network three-phase data
According to, B phase related data and C phase related data；
Step 2: each phase data of the ABC isolated being subjected to complex field mark respectively and is changed, and establishes respective PQ points respectively
Constant Jacobian matrix matrix needed for solution；
Further, step 2 the following steps are included:
Step 201: power reference value is selected in complex fieldVoltage reference valueAre as follows:
In formula, SbaseFor apparent energy, UbaseFor benchmark voltage magnitude, e is the bottom of natural logrithm, and j is imaginary unit, φs
For corresponding power reference angle, φvFor corresponding voltage reference angle.
Then acquire corresponding impedance reference value are as follows:
In formula, ZbaseFor impedance reference value modulus value, subscript * indicates conjugation.
So, the mark of power changes form in complex field are as follows:
In formula, m=1,2,3..., t, t are injecting power number；Pm, QmIndicate m-th injection active power, reactive power；
Subscript cpu indicates that the mark under complex field is changed；Subscript pu indicates that the mark under traditional real number field is changed.
The mark of line impedance changes form in complex field are as follows:
αi=arctan (Xi/Ri) (6)
In formula, i=1,2,3 ..., l, l are circuitry number；Ri,XiRespectively indicate i-th of resistance and reactance；αiIt indicates i-th
The tangent angle of reactance and resistance；Subscript cpu indicates that the mark under complex field is changed；Subscript pu indicates the mark under traditional real number field
Change.
It can be seen that the impedance value after complex field subscript is changed and the impedance value phase after the change of real number field subscript from formula (5)
Than increasing a power reference angle φs, this be complex field fast decoupled trend can used in the power distribution network one to closing
Important factor.By the experience of traditional power transmission network Load flow calculation, line parameter circuit value R < < X be fast decoupled Load flow calculation must
It wants condition, and traditional power transmission network fast decoupled power flow algorithm to be difficult to adapt to the major reason of power distribution network, is described in detail below
The method how to change with complex field mark, solves the problems, such as high R/X ratio in distribution line.
Step 202: by formula (5) it is found that complex field mark corresponding resistance and reactance after changing are as follows:
To make Ri,cpu<<Xi,cpu, then
That is,
By formula (9) it is found that αiIt is to be determined by the intrinsic parameter of route, therefore a suitable power reference angle can be selected
φsSo thatSo that Ri,cpu<<Xi,cpu, to enhance fast decoupled power flow algorithm in distribution network
Adaptability.By formula (1)~(9) it may be noted that voltage reference angle φvSelection, have no effect on, therefore take φv=0 is asked with mitigating
The complexity of topic calculates to simplify.
Step 203: complex field power reference angle φsSelection to fast decoupled power flow algorithm distribution network application have
Decisive role.From the property of tangent function:
By formula (10) it is found that as (αi+φs) level off to pi/2 left side when, tan (αi+φs) level off to+∞；As (αi+
φs) level off to pi/2 right side when, tan (αi+φs) leveling off to-∞, analysis can obtain, both approach modes can meet quickly
The requirement of decomposition computation, only as selected φsSo that tan (the α of part routei+φs) level off to-∞ when, these routes R/X
Ratio will appear the case where negative value, this does not meet the normal physical phenomenon that line parameter circuit value is shown.Therefore it is selected in complex field
Power reference angle φsWhen, it is intended merely to reduce the influence of R/X ratio as much as possible, matched with enhancing fast decoupled Load flow calculation
The applicability of power grid, there is no actual physical significances for each parameter of route after the change of complex field mark.
The present invention chooses power reference angle φ using the method for optimizations, initially set up following objective function:
In formula, l is circuitry number.
According to formula (11), function f (φs) indicate the α of all routesiPhase shift φsLater, to the left or to the right close to the journey of pi/2
Degree.Institute is in the hope of making objective function f (φs) the smallest φ of valuesAs our optimal solutions to be looked for.Enable f (φs) to φsLocal derviation
It is zero, φ can be acquiredsOptimal solution are as follows:
In formula, l is circuitry number.
Step 3: the voltage magnitude of each phase and phase angle being set to initial value, and iteration count is set to 1；
Step 4: current compensation is carried out to A phase using decoupling compensation model, a fast decoupled is then carried out and calculates, amendment
The voltage magnitude and phase angle of A phase part.Similarly, the calculating similar with A is carried out to B phase, C phase, and corrects B phase C phase portion respectively
The voltage magnitude and phase angle divided.
Further, step 4 the following steps are included:
Step 401: shown in the accurate three-phase π type Type Equivalent Circuit Model attached drawing 2 of distribution network line.In figure,For route
Direct-to-ground capacitance part；For the admittance part of route, then the impedance Z of routeabcIt can indicate are as follows:
Z in formulaaa,Zbb,ZccRespectively A, B, C phase self-impedance part, Zab=ZbaFor AB mutual impedance, Zac=ZcaFor AC
Mutual impedance, Zbc=ZcbFor BC mutual impedance；
It can derive the admittance equations of distribution line before decoupling are as follows:
Wherein,
In formulaRespectively A, B, C phase self-admittance part,For the alternate transadmittance of AB,For
The alternate transadmittance of AC,For the alternate transadmittance of BC；In formulaRespectively A, B, C phase self-capacitance part,For the alternate mutual capacitance of AB,For the alternate mutual capacitance of AC,For the alternate mutual capacitance of BC；
In view of the asymmetry of power distribution network three-phase line parameter, the i.e. inequality of three-phase mutual impedance:
Decoupling based on order components can not eliminate the coupling between each sequence, thus the present invention using alternate decoupling compensation model into
The alternate decoupling of row three-phase.Formula (14) is directly unfolded, and the element on off-diagonal is moved on to the right of equation, through whole
After reason, further derive:
In formula, k is the number of iterations,Indicate the compensation electric current of A phase branch ij,Indicate the compensation of branch i-node
Electric current,Indicate the compensation electric current of branch j node.3 institute of A phase decoupling compensation equivalent-circuit model attached drawing can be obtained by formula (17)
Show.
Step 402: according to A phase decoupling compensation current equivalence circuit model, the section at A phase line both ends after decoupling can be found out
Point injecting compensating electric current:
In formula, k is the number of iterations,Indicate the equivalent injecting compensating electric current of node i,Indicate that node j's is equivalent
Injecting compensating electric current.
Finally by the node injecting compensating electric current at A phase line both ends after decoupling, it is converted to corresponding compensation injecting power and adds
Enter into two end node injecting power of route, complete the alternate decoupling work of the first step:
In formula, k is the number of iterations, Δ Sa,iIndicate the equivalent injecting compensating power of node i, Δ Sa,jExpression node j etc.
Injecting compensating power is imitated, subscript * indicates conjugation.
Similarly, the decoupling compensation model of A phase can be extrapolated to B phase and C phase, and details are not described herein.Adopt alternate Decoupled Model
Later, the coupling between ABC three-phase is then correspondingly converted to the compensation electric current of each phase, to realize asymmetric three-phase line
Alternate decoupling, can then be converted to the calculating of each phase accordingly for the calculating of power distribution network Three-phase Power Flow, accelerate computational efficiency.
Step 5: judging whether the maximum value of ABC phase correction amount absolute value reaches convergent requirement, if reaching convergent requirement
Then stop calculating, if not up to convergent requirement, continues cycling through, until reaching convergent requirement.Finally output is as a result, knot
Beam calculates.
Embodiment
The example test macro used has selected 3 for IEEE13, IEEE34, IEEE123 three-phase imbalance distribution system
Kind method is compared with method (hereinafter referred to as this paper algorithm) proposed by the invention.It for convenience of elaboration, now provides: herein
Algorithm, that is, alternate decoupling (pu)+PQ decomposition method (cpu)；Method 1 is conventional three-phase newton power flow algorithm (pu)；Method 2 is alternate
It decouples (pu)+single-phase Newton method (pu)；Method 3 is alternate decoupling (pu)+PQ decomposition method (pu).The area of this paper algorithm and method 3
It is not that this paper algorithm uses complex field PQ decomposition method in second step phase in Uncoupled procedure, and is solved in the second step phase of method 3
Coupling process then uses real number field PQ decomposition method.
Attached drawing 4 is the voltage magnitude calculated result of IEEE13 three-phase unbalance system, selects conventional three-phase newton tide herein
Benchmark of the flow algorithm as context of methods calculated result, convergence precision 10-5.In figure " x ,+, * " respectively indicates A phase, B phase
With the calculated result of C phase conventional three-phase Newton method, " ▽, zero, " respectively indicates the calculating knot of A phase, B phase and C phase context of methods
Fruit.It can be seen that the calculated result of the calculated result and conventional three-phase Newton method of context of methods is almost the same from attached drawing 4, verifying
The validity of context of methods.
Attached drawing 5 is the relational graph of IEEE34 node system four kinds of method the number of iterations and convergence precision, wherein abscissa table
What is shown is the number of iterations, and that ordinate indicates is the negative logarithm-log of convergence precision10(·).It can be seen that from attached drawing 5
In IEEE34 bus test system, as classical trend solution, method 1 has good convergence property, iterates to calculate by 5 times
Afterwards, precision can reach 10-5.Due to the influence of R/X ratio route high in route, method 3 cannot be effective in the test macro
Convergence demonstrates traditional fast decoupled power flow algorithm in distribution power system load flow calculation the problem of bad adaptability.Context of methods and side
Method 2 has similar convergence property, mostly changes in approximately linear, wherein method 2 can reach 10 in 10 times or so computational accuracies-5, and context of methods can reach 10 in 12 times or so computational accuracies-5.As it can be seen that method proposed in this paper has certain engineering real
Use prospect.
Attached drawing 6 is the comparison of three kinds of method average calculation times in test macro IEEE34, IEEE123, convergence precision
It is set to 10-5, average calculation times are the average results that every kind of method runs 100 times.This paper algorithm leads to it can be seen from attached drawing 6
The alternate decoupling of ABC, PQ decoupling in each phase phase of ABC are crossed, the dimension of Jacobian matrix is greatly reduced, reduces Jacobian matrix
Memory space.In addition have benefited from the constant of Jacobian matrix, mentioned algorithm, which effectively saves, herein calculates the time, improves meter
Calculate efficiency.From attached drawing 5 and attached drawing 6 as can be seen that although this paper algorithm is on constringency performance and method 2 is very close, herein
Algorithm is substantially better than method 2 in computational efficiency.Therefore in general, this paper algorithm has great advantage in calculated performance,
It is able to achieve three-phase imbalance power distribution network fast and efficiently Load flow calculation.
Table 1 illustrates the situation of change that complex field mark in IEEE34 test macro changes front and back B phase part line parameter circuit value.
1 IEEE34 test macro complex field mark of table changes the variation of front and back B phase part line parameter circuit value
As can be seen from Table 1, route after complex field mark is changed | R/X | value has obtained a degree of reduction.Partial line
There is the phenomenon that R value is negative in road parameter, and this phenomenon is just consistent with the analysis result of the second trifle herein, belongs to and just often occurs
As complex field mark changes the influence for being intended merely to reduce route R/X ratio as much as possible, is existed with enhancing fast decoupled Load flow calculation
The applicability of power distribution network, there is no actual physical significances for each parameter of route after the change of complex field mark.
It can be seen that the comprehensive of the bilingual coupling power flow algorithm of power distribution network proposed in the present invention from the result of sample calculation analysis
Other algorithms can be substantially better than, have the characteristics that calculating is efficient, convergence is linear, adaptable.Therefore in three-phase imbalance distribution
Within the scope of the Practical of net Load flow calculation, the bilingual coupling power flow algorithm of power distribution network proposed in the present invention has good application
Value.
Claims (1)
1. a kind of bilingual coupling tidal current computing method of power distribution network, which comprises the following steps:
Step 1: power distribution network three-phase data is carried out data separating, is separated into A phase related data, B by input power distribution network three-phase data
Phase related data and C phase related data；
Step 2: each phase data of the ABC isolated being subjected to complex field mark respectively and is changed, and establishes respective PQ respectively and decomposes institute
The Jacobian matrix needed；
Step 201: power reference value is selected in complex fieldVoltage reference valueAre as follows:
In formula, SbaseFor apparent energy, UbaseFor benchmark voltage magnitude, e is the bottom of natural logrithm, and j is imaginary unit, φsFor phase
The power reference angle answered, φvFor corresponding voltage reference angle；
Then acquire corresponding impedance reference value are as follows:
In formula, ZbaseFor impedance reference value modulus value, subscript * indicates conjugation；
So, the mark of power changes form in complex field are as follows:
In formula, m=1,2,3..., t, t are injecting power number；Pm, QmIndicate m-th injection active power, reactive power；Subscript
Cpu indicates that the mark under complex field is changed；Subscript pu indicates that the mark under traditional real number field is changed；
The mark of line impedance changes form in complex field are as follows:
αi=arctan (Xi/Ri) (6)
In formula, i=1,2,3 ..., l, l are circuitry number；Ri,XiRespectively indicate i-th of resistance and reactance；αiIndicate i-th of reactance
With the tangent angle of resistance；Subscript cpu indicates that the mark under complex field is changed；Subscript pu indicates that the mark under traditional real number field is changed, Zi
For the impedance of branch i；
Find out from formula (5), the impedance value after the change of complex field subscript increases compared with the impedance value after the change of real number field subscript
One power reference angle φs, this is complex field fast decoupled trend vital factor used in power distribution network；By
The experience of traditional power transmission network Load flow calculation, line parameter circuit value R < < X are the necessary conditions of fast decoupled Load flow calculation, and are passed
System power transmission network fast decoupled power flow algorithm is difficult to adapt to the major reason of power distribution network, how is described in detail below with complex field
The method changed is marked, solves the problems, such as high R/X ratio in distribution line；
Step 202: by formula (5) it is found that complex field mark corresponding resistance and reactance after changing are as follows:
To make Ri,cpu<<Xi,cpu, then
That is,
By formula (9) it is found that αiIt is to be determined by the intrinsic parameter of route, therefore select a suitable power reference angle φsSo thatSo that Ri,cpu<<Xi,cpu, to enhance fast decoupled power flow algorithm in the adaptability of distribution network；
Voltage reference angle φ is noticed by formula (1)~(9)vSelection, have no effect on, therefore take φv=0 to mitigate the complexity of problem,
It is calculated to simplify；
Step 203: complex field power reference angle φsSelection to fast decoupled power flow algorithm distribution network application have certainly
Qualitative effect, from the property of tangent function:
By formula (10) it is found that as (αi+φs) level off to pi/2 left side when, tan (αi+φs) level off to+∞；As (αi+φs) become
When being bordering on the right side of pi/2, tan (αi+φs) leveling off to-∞, analysis can obtain, both approach modes can meet fast decoupled meter
The requirement of calculation, only as selected φsSo that tan (the α of part routei+φs) level off to-∞ when, the R/X ratio meeting of these routes
There is the case where negative value, this does not meet the normal physical phenomenon that line parameter circuit value is shown, therefore power base is selected in complex field
Quasi- angle φsWhen, it is intended merely to reduce the influence of R/X ratio as much as possible, to enhance fast decoupled Load flow calculation in power distribution network
Applicability, there is no actual physical significances for each parameter of route after the change of complex field mark；
Power reference angle φ is chosen using the method for optimizations, initially set up following objective function:
In formula, l is circuitry number；
According to formula (11), function f (φs) indicate the α of all routesiPhase shift φsLater, to the left or to the right close to the degree of pi/2；
Institute is in the hope of making objective function f (φs) the smallest φ of valuesAs our optimal solutions to be looked for, enable f (φs) to φsLocal derviation be
Zero, φ can be acquiredsOptimal solution are as follows:
In formula, l is circuitry number；
Step 3: the voltage magnitude of each phase and phase angle being set to initial value, and iteration count is set to 1；
Step 4: current compensation being carried out to A phase using decoupling compensation model, a fast decoupled is then carried out and calculates, correct A phase
Partial voltage magnitude and phase angle；Similarly, the calculating similar with A is carried out to B phase, C phase, and corrects B phase C phase part respectively
Voltage magnitude and phase angle；
The decoupling compensation model are as follows:
Step 401: in the accurate three-phase π type Type Equivalent Circuit Model of distribution network line,For the direct-to-ground capacitance part of route；
For the admittance part of route, then the impedance Z of routeabcIt indicates are as follows:
Z in formulaaa,Zbb,ZccRespectively A, B, C phase self-impedance part, Zab=ZbaFor AB mutual impedance, Zac=ZcaIt is alternate for AC
Mutual impedance, Zbc=ZcbFor BC mutual impedance；
It can derive the admittance equations of distribution line before decoupling are as follows:
Wherein,
In formulaRespectively A, B, C phase self-admittance part,For the alternate transadmittance of AB,For AC phase
Between transadmittance,For the alternate transadmittance of BC；In formulaRespectively A, B, C phase self-capacitance part,For
The alternate mutual capacitance of AB,For the alternate mutual capacitance of AC,For the alternate mutual capacitance of BC；
In view of the asymmetry of power distribution network three-phase line parameter, the i.e. inequality of three-phase mutual impedance:
Decoupling based on order components can not eliminate the coupling between each sequence, therefore the phase of three-phase is carried out using alternate decoupling compensation model
Between decouple, formula (14) is directly unfolded, and the element on off-diagonal is moved on to the right of equation, after collated, into one
Step derives:
In formula, k is the number of iterations,Indicate the compensation electric current of A phase branch ij,Indicate the compensation electric current of branch i-node,The compensation electric current for indicating branch j node, can be obtained shown in A phase decoupling compensation equivalent-circuit model as formula (17)；
Step 402: according to A phase decoupling compensation current equivalence circuit model, finding out the node injection benefit at A phase line both ends after decoupling
Repay electric current:
In formula, k is the number of iterations,Indicate the equivalent injecting compensating electric current of node i,Indicate that the equivalent injection of node j is mended
Repay electric current；
Finally by the node injecting compensating electric current at A phase line both ends after decoupling, it is converted to corresponding compensation injecting power and is added to
In two end node injecting power of route, the alternate decoupling work of the first step is completed:
In formula, k is the number of iterations, Δ Sa,iIndicate the equivalent injecting compensating power of node i, Δ Sa,jIndicate the equivalent note of node j
Enter and compensate power, subscript * indicates conjugation；
Similarly, the decoupling compensation model of A phase is extrapolated to B phase and C phase, the coupling after alternate Decoupled Model, between ABC three-phase
Conjunction is then correspondingly converted to the compensation electric current of each phase, so that the alternate decoupling of asymmetric three-phase line is realized, for distribution
The calculating of net Three-phase Power Flow is then converted to the calculating of each phase accordingly, accelerates computational efficiency；
Step 5: judging whether the maximum value of ABC phase correction amount absolute value reaches convergent requirement, stop if reaching convergent requirement
It only calculates, if not up to convergent requirement, continues cycling through, until reaching convergent requirement；Finally output is as a result, terminate meter
It calculates.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610839531.9A CN106208082B (en) | 2016-09-21 | 2016-09-21 | A kind of bilingual coupling tidal current computing method of power distribution network |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610839531.9A CN106208082B (en) | 2016-09-21 | 2016-09-21 | A kind of bilingual coupling tidal current computing method of power distribution network |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106208082A CN106208082A (en) | 2016-12-07 |
| CN106208082B true CN106208082B (en) | 2019-03-12 |
Family
ID=58067227
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610839531.9A Active CN106208082B (en) | 2016-09-21 | 2016-09-21 | A kind of bilingual coupling tidal current computing method of power distribution network |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106208082B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107104430A (en) * | 2017-05-05 | 2017-08-29 | 河海大学 | A kind of bilingual coupling method for estimating state of power distribution network |
| CN108964047A (en) * | 2018-08-07 | 2018-12-07 | 河海大学 | A kind of State Estimation for Distribution Network |
| CN109245107A (en) * | 2018-10-26 | 2019-01-18 | 贵州电网有限责任公司 | A kind of uneven distribution power system load flow calculation considering distributed generation resource access |
| CN110728032B (en) * | 2019-09-20 | 2021-04-27 | 河海大学 | Quick power flow calculation method for electricity-heat interconnection comprehensive energy system considering ring network |
| CN113655262B (en) * | 2021-04-12 | 2024-01-16 | 重庆大学 | Voltage measurement self-decoupling method based on multidimensional equivalent capacitance calculation |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103065059A (en) * | 2013-01-29 | 2013-04-24 | 河海大学 | Method for calculating power flow of radial power distribution network based on variable substitution |
| CN103353918A (en) * | 2013-04-24 | 2013-10-16 | 国家电网公司 | Method for selecting initial values of three-phase Newton-Raphson method power flow calculation of power distribution network |
| CN105071387A (en) * | 2015-08-13 | 2015-11-18 | 河海大学 | Power distribution network rapid decomposition state estimation method based on complex number domain per unit theory |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| TWI264864B (en) * | 2005-04-08 | 2006-10-21 | Univ Chang Gung | Power flow calculation method of power grid with unified power flow controller |
-
2016
- 2016-09-21 CN CN201610839531.9A patent/CN106208082B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103065059A (en) * | 2013-01-29 | 2013-04-24 | 河海大学 | Method for calculating power flow of radial power distribution network based on variable substitution |
| CN103353918A (en) * | 2013-04-24 | 2013-10-16 | 国家电网公司 | Method for selecting initial values of three-phase Newton-Raphson method power flow calculation of power distribution network |
| CN105071387A (en) * | 2015-08-13 | 2015-11-18 | 河海大学 | Power distribution network rapid decomposition state estimation method based on complex number domain per unit theory |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106208082A (en) | 2016-12-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106208082B (en) | A kind of bilingual coupling tidal current computing method of power distribution network | |
| Hameed et al. | A modified backward/forward sweep load flow method for islanded radial microgrids | |
| Ju et al. | An extension of FBS three-phase power flow for handling PV nodes in active distribution networks | |
| Penido et al. | Three-phase power flow based on four-conductor current injection method for unbalanced distribution networks | |
| Balamurugan et al. | Review of power flow studies on distribution network with distributed generation | |
| JP5815382B2 (en) | Power flow analysis for balanced distribution systems | |
| Mohn et al. | Tracing PV and QV curves with the help of a CRIC continuation method | |
| CN103094900B (en) | Distributed generation power distribution network three-phase load flow calculation method taking phase sequence mixing method into consideration | |
| CN111625914A (en) | Short-circuit current calculation method and system | |
| CN103971026A (en) | General method for calculating tide of positive power distribution networks | |
| CN104794531A (en) | MDCOPF (modified direct current optimal power flow) method based on grid loss equivalent load model | |
| CN105162099A (en) | Operational curved surface method for determining asymmetric short-circuit current of distributed generation access power grid | |
| Abaali et al. | Comparison of newton Raphson and gauss seidel methods for power flow analysis | |
| Sun et al. | A fast and robust load flow method for distribution systems with distributed generations | |
| CN109149583A (en) | Active power distribution network is succinctly pushed forward back substitution tidal current computing method | |
| Saleh | The Formulation of a Power Flow Using $ d\text {--} q $ Reference Frame Components—Part II: Unbalanced $3\phi $ Systems | |
| CN111049144A (en) | Active power distribution network load flow calculation PV node processing method based on constant change rate matrix | |
| Yang et al. | Three-phase power flow calculations using initial voltage estimation method for unbalanced distribution networks | |
| CN105896604A (en) | Calculation method of ultimate capacity of distributed photovoltaic power generation system connected to power network | |
| CN112688332B (en) | Hybrid alternating current-direct current microgrid power flow analysis method and system based on sequence component coordinates | |
| Obrenić et al. | A novel intervals-based algorithm for the distribution short-circuit calculation | |
| Zhu et al. | Integrated power flow solution for interconnected transmission and distribution network | |
| CN111900738B (en) | Three-phase unbalanced load flow calculation method based on compensation algorithm | |
| Presada et al. | State estimation in power systems with FACTS devices and PMU measurements | |
| Nurdin et al. | Backward Forward Sweep Algorithm for Unbalanced Three-Phase Power Flow Analysis in Distribution Systems Containing Voltage Regulator |
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 |