CN104281749A - Innovation graph method based DG (distributed generation) power distribution network included topology identification method - Google Patents
Innovation graph method based DG (distributed generation) power distribution network included topology identification method Download PDFInfo
- Publication number
- CN104281749A CN104281749A CN201410530548.7A CN201410530548A CN104281749A CN 104281749 A CN104281749 A CN 104281749A CN 201410530548 A CN201410530548 A CN 201410530548A CN 104281749 A CN104281749 A CN 104281749A
- Authority
- CN
- China
- Prior art keywords
- value
- branch road
- topology
- calculating
- new breath
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses an innovation graph method based DG (distributed generation) power distribution network included topology identification method. The method specifically includes the steps: modeling a power distribution network and determining links and branches; calculating state estimation results of k time; calculating innovation values with measuring branch circuits; calculating the calculated innovation values of the whole network; calculating difference vectors of the branches equipped with measurement; judging whether part of the branch circuits are great in difference vectors and exist independently or not; judging whether the branch circuits with large difference vectors form bad circuits or not; regaining the calculated innovation values of the whole network after eliminating impact of bad data; calculating estimation trend; judging whether the DG is achieved in small-power situation or not; judging whether innovation values of links of the DG are not zero; judging whether branch circuits with calculated innovation absolute values great and estimation tread values small and close to zero exist or not; recognizing topology error branch circuits according to the judgment results. The method is applicable to the DG power distribution network included topology identification, low in measurement configuration requirement and easy for engineering practice.
Description
Technical field
The present invention relates to a kind of topology Identification method of power distribution network, particularly a kind of topology Identification method of the power distribution network to the shortage measurement information containing DG.
Background technology
The interconnected mode of current electric system centralized power generation, long distance power transmission and bulk power grid shows great superiority in resource allocation and system stability, but there is Local Accidents and easily spread drawbacks such as causing large area blackout.Along with the exhaustion of conventional energy resources and environmental pollution increasingly sharpen, environmental protection, efficiently and flexibly distributed power generation (Distributed Generation, DG) obtain concern and the attention of countries in the world, and development in recent years is rapid.
Distributed power generation generally refer to generated output in the thousands of watts of miniaturizations to 50 megawatts, scale, be dispersed in the compact electrical generating systems of the access distribution system near user.Islet operation is the new method of operation occurred in power distribution network along with distributed generation technology development, refers to by the independent running status of powering to sub-load of DG.After distribution is broken down, under the prerequisite ensureing power system safety and stability, maintain DG normal power supply as much as possible under island operation state, to reduce power failure area, improve power supply reliability.But consider voltage and the frequency instability of islet operation, therefore the electric energy requiring powerful application scenario and DG to send in standard " IEEE Std 1547-2003 " also needs to deliver to electrical network except satisfied regional load supply, there is isolated island in that case, require that DG works in island operation state; And the electric energy that low power application scenario and DG send only meets some areas workload demand, if there is island state, start island protect function, DG stops load supplying towards periphery.
In the actual telecontrol system of power distribution network, once switch remote signalling does not conform to switch physical location, be then called network structure mistake, that is network topology mistake.The prerequisite that correct power distribution network topology is state estimation, localization of fault and fault recovery etc. are applied, the topological structure of mistake can cause state estimation result incorrect, even do not restrain, make yardman cannot obtain the running status of electrical network, thus the operations such as load transfer cannot be carried out, and Topology Error also possibility causing trouble positioning result mistake, cause the consequences such as power supply quality is poor, power supply reliability reduction, therefore, significant to carrying out topology Identification containing the power distribution network of DG.
At present, the method for topology Identification mainly contains Sudden Changing Rate startup method, network structure search identification method and Parameter Estimation Method.The remote measurement that Sudden Changing Rate startup method Main Basis distribution automation terminal is uploaded and remote signalling amount are carried out analysis and are judged Topology Error section; Network structure search identification method is a kind of method of state estimation, and possible distribution network textural association is enumerated and carried out state estimation calculating by one by one, until search is met the network structure of criterion; Parameter Estimation Method obtains the parameter of branch road by parameter estimation and then judges this membership, reaches the object of topology Identification.It is high to distribution automation hardware requirement that Sudden Changing Rate starts method, and current most of distribution automation fails to reach this requirement; Network structure search identification method need carry out repeatedly state estimation and calculate, and inapplicable to multiple Topology Error situation; Parameter Estimation Method is very high to the requirement of measure configuration, makes it inapplicable to measuring relatively deficient power distribution network.And said method does not all consider the impact that distributed power source accesses, and inapplicable for the power distribution network of metric data scarcity.
Summary of the invention
The technical issues that need to address of the present invention are to provide a kind of lower to metric data requirement, and the method for topology Identification exactly can be carried out to the power distribution network containing D G, occur to avoid the fault handling time situation that is long, power supply reliability reduction caused because of Topology Error problem.
For solving the problems of the technologies described above, the technical solution used in the present invention is as follows.
Based on the topology Identification method containing DG power distribution network of innovation graph approach, comprise the following steps:
Step one, to Modeling of Distribution Network, determines chord and tree;
Step 2, calculates the state estimation result in k moment;
Step 3, calculates the new breath value of all measurement branch roads;
Step 4, the new breath value according to chord calculates the calculating newly breath value obtaining the whole network;
Step 5, props up the difference vector of the difference calculating its true new breath value and calculate new breath value to the tree being equipped with measurement;
Step 6, judges whether to exist partial branch difference vector very large and independently to exist, carries out step 8 if exist; Step 9 is carried out if do not exist;
Step 7, judges whether that there is the very large branch road of difference vector forms bad loop, carries out step 8 if exist; Step 9 is carried out if do not exist;
Step 8, obtains the calculating newly breath value of the whole network again according to judged result after getting rid of the impact of bad data;
Step 9, the state estimation result according to the k moment tries to achieve estimation trend with the new breath value of calculating;
Step 10, judges whether contained DG power distribution network is miniwatt occasion; If large-power occasions carries out step 11; If miniwatt occasion carries out step 12;
Step 11, the calculating judging whether to exist certain branch road newly breath value is very large but estimate that trend value is very little of zero; If the calculating that there is certain branch road newly breath value is very large, but estimates that trend value is very little of zero, so this branch road exists and disconnects Topology Error suddenly, and this branch road is actual is off state, and topology Identification terminates; If do not exist, show that network is without Topology Error, topology Identification terminates;
Step 12, judges whether that the new breath value that there is DG place chord is non-vanishing; If the new breath value of DG place chord is non-vanishing carry out step 13, if the new breath value of DG place chord is zero, carry out step 14;
Step 13, if the new breath value of all DG places chord is zero, judges whether to exist that to calculate new breath absolute value very large, and estimate trend value very little close to zero branch road; If there is branch road, to meet new breath value absolute value comparatively large, estimate trend value close to zero condition, so Topology Error branch road is the most downstream branch road meeting above-mentioned condition, and this branch road virtual condition is off, and topology Identification terminates; If do not exist, show that, without Topology Error in network, topology Identification terminates;
Step 14, if the new breath value that there is DG place chord is non-vanishing, be then the situation of DG at Topology Error downstream branch, should the new breath value of DG place chord be changed into its opposite number, again calculate that the whole network calculates the estimation trend value that new breath value and step 9 calculate each branch road, if newly breath value is very large in the calculating of certain branch road according to step 4, but estimate that trend value is very little of zero, then the existence of this branch road disconnects Topology Error suddenly, and this branch road is actual is off state, and topology Identification terminates.
The state estimation result calculating the k moment described in step 2 of the present invention carries out state estimation calculating with the measurement amount in k moment and network topology structure, and state estimation calculates according to formula one:
formula one
In formula,
for state variable, branch current magnitudes and phase angle are chosen as its state variable for power distribution network,
for measurement amount,
for measuring function,
for measuring weight matrix;
Separate the nonlinear equation obtained with Niu Lafa, kth time iteration amount is:
Wherein
, for measuring Jacobian matrix,
for information matrix, iterative computation until
be less than setting threshold value (
), quantity of state now
required by being, and then voltage, the Power Flow Information of whole network can be obtained.
New breath value described in step 3 of the present invention obtains according to formula two:
formula two
In formula
for k+1 moment measurement amount,
for the k moment measures estimated value,
be the new breath value in k+1 moment.
Difference vector described in step 5 of the present invention obtains according to formula three:
formula three
In formula,
for newly being ceased the calculating newly breath value calculating gained by chord,
for the true new breath value of the branch road that step 3 obtains.
Estimation trend described in step 9 of the present invention obtains according to formula four:
formula four
In formula
for the k moment measures estimated value and state estimation result.
Owing to have employed above technical scheme, the invention technological progress is as follows.
The present invention is directed to power distribution network vast in territory, measure deficient feature, the innovation graph approach with forecast function is utilized to detect, identify the abnormal event such as bad data and Topology Error, and the access for DG is carried out on innovation graph approach the impact that power distribution network topology Identification causes and is proposed corresponding solution, thus guarantee that distribution network topological structure is accurate, and then reduce because topology is forbidden caused fault handling time long, the problem that power supply reliability reduces, compensate for existing discrimination method requires high to metric data redundance, the deficiencies such as calculation of complex, be specially adapted to the topology Identification containing DG power distribution network.The present invention is directed to the criterion that the capacity of DG and on-position propose the impact of innovation graph approach topology Identification to adapt with it, all can accurate recognition to the various Topology Errors containing DG power distribution network, adapt to the trend of increasing distributed power source access in power distribution network; Principle is simple, less demanding to measure configuration, is easy to engineering practice.
Accompanying drawing explanation
Fig. 1 is process flow diagram of the present invention;
Fig. 2 is the circuit diagram of the present embodiment.
Embodiment
Below in conjunction with the drawings and specific embodiments, the present invention is further elaborated.
In power distribution network topology Identification process, Topology Error can be defined as switch reality and not conjugate and the change of terminal information upload reflection on off state, namely report in Topology Error and real system and change in topology occurs, but because the causes such as communication disruption cause main station system not make respective reaction, namely do not report Topology Error.Report Topology Error situation easily processes, do not report in Topology Error situation again to disconnect suddenly and the most difficult identification of non-reporting errors, therefore only topology Identification method is set forth for this situation, and from succinct convenient consideration, will disconnect suddenly and not report that Topology Error is referred to as disconnecting Topology Error suddenly.
Based on the topology Identification method containing DG power distribution network of innovation graph approach, its process flow diagram as shown in Figure 1, specifically comprises the following steps:
Step one, to Modeling of Distribution Network, determine chord and tree:
Electricity distribution network model need set up ground node, and all constructs a branch road over the ground to each node.Measure because substation transformer each in distribution is all equipped with and can load bus be equivalent to, therefore choose all nodes except feeder line root node over the ground branch road over the ground branch road be chord, tree then elected as by all the other branch roads of this feeder line.
Step 2, calculates the state estimation result in k moment:
Carry out state estimation calculating with the measurement amount in k moment and network topology structure, the raw data mart modeling of distribution is become ripe data by state estimation, determines accuracy and the degree of functioning of other analytical applications to a certain extent.
Distributed power source can be considered as PQ node, PV node, PI node, PQ(V respectively according to the difference of its control characteristic in state estimation calculates) node and I node.
Use at present application maximum be that basic weighted least-squares method carries out state estimation, state estimation calculates according to formula one:
formula one
In formula,
for state variable, branch current magnitudes and phase angle are chosen as its state variable for power distribution network,
for measurement amount,
for measuring function,
for measuring weight matrix.
Objective function gets minimum value
on its extreme point, therefore have
.
Separate the nonlinear equation obtained with Niu Lafa, kth time iteration amount is:
Wherein
, for measuring Jacobian matrix,
for information matrix, iterative computation until
be less than setting threshold value (
), quantity of state now
required by being, and then voltage, the Power Flow Information of whole network can be obtained.
Step 3, calculates the new breath value of all measurement branch roads:
By the measurement amount in k+1 moment
deduct the state estimation result in k moment
namely measure estimated value, obtain the new breath value in k+1 moment
.
New breath value obtains according to formula two:
formula two
In formula
for k+1 moment measurement amount,
for the k moment measures estimated value,
be the new breath value in k+1 moment.
To be made up of new breath value and the figure meeting circuit theory is innovation graph.Follow first, second law of Kirchhoff between new breath value in innovation graph, when some branch currents represented by new breath value of known complete network structure and some, branch current every other in innovation graph can calculate.Select a tree, form independent loop calculate branch current, distribution network model as shown in Figure 1 according to the structural system of tree, dotted line represents chord, and solid line represents tree.Chord electric current is made to be known loop current vector
, just can calculate branch currents all in innovation graph
for:
In formula, C is branch road-path incidence matrix.
In electric system, custom adopts branch power, when employing branch road active power replaces branch current, and can obtain following formula when adopting DC Model to calculate:
In formula
for chord truly newly ceases,
for branch road calculates new breath.
Step 4, the new breath value according to chord calculates the calculating newly breath value obtaining the whole network, obtains the calculating newly breath value that all trees prop up
.
Step 5, props up for being equipped with the tree measured, and its measured value and the difference measuring estimated value are true new breath value
, simultaneously using the difference of the calculating newly breath value and the true new breath value that are calculated gained by chord as difference vector.
Bad data in power distribution network does not follow circuital law, thus by loop method from chord newly breath value calculate the calculating that obtains newly breath value be consistent data acquisition, when they with those purely according to measuring value deduct measure true new breath value that estimated value obtains compare time, namely difference can reveal.It is difference vector that definition calculates new breath value with the difference of true new breath value.
Difference vector obtains according to formula three:
formula three
In formula,
for newly being ceased the calculating newly breath value calculating gained by chord,
for the true new breath value of the branch road that step 3 obtains.
When larger difference vector element is in stochastic distribution state, can judges that the tree that bad data is positioned at difference vector larger props up, this measuring value correction be become the difference of its value and difference vector; Larger difference vector is all there is on the branch road of all measured values in a loop except chord, when namely forming bad loop, then illustrate that the measurement of this chord is bad data, this measuring value correction is become its value deducts difference vector equally, can get rid of the harmful effect that may cause topology Identification of bad data.
Step 6, judges whether to exist partial branch difference vector very large and independently to exist, if having partial branch difference vector very large and independently to exist, then can be determined as a tree bad data, carry out step 8; Step 9 is carried out if do not exist.
Step 7, judges whether that there is the very large branch road of difference vector forms bad loop, if there is the very large branch road of difference vector to form bad loop, then can think and form bad loop, and quarter, dish was not for existing chord bad data, carried out step 8; Step 9 is carried out if do not exist.
Step 8, if there is bad data, is modified to measuring value and is deducted difference vector, namely returned step 3.After the impact getting rid of bad data, recalculate the new breath value obtaining each chord, and calculate the calculating newly breath value obtaining the whole network
.
Step 9, the state estimation result according to the k moment tries to achieve estimation trend with the new breath value of calculating; The state estimation result in the k moment of gained in new breath value and step 2 is calculated by gained in step 4
be added, its result is as estimation trend
.
Estimation trend obtains according to formula four:
formula four
In formula
for the k moment measures estimated value and state estimation result.
Step 10, judges whether the DG power distribution network that feeder line is got involved is miniwatt occasion; If large-power occasions carries out step 11; If miniwatt occasion carries out step 12.
Step 11, if high-power applications occasion, the calculating judging whether to exist certain branch road newly breath value is very large but estimate that trend value is very little of zero.If the calculating that there is certain branch road newly breath value is very large, but estimates that trend value is very little of zero, so this branch road exists and disconnects Topology Error suddenly, and this branch road is actual is off state, and topology Identification terminates; If do not exist, show that network is without Topology Error, topology Identification terminates.
Step 12, if low-power applications occasion, then needs the relative position judging DG and Topology Error branch road, and the new breath value namely by judging whether to exist DG place chord is non-vanishing to be determined.If the new breath value of DG place chord is non-vanishing carry out step 13, if the new breath value of DG place chord is zero, carry out step 14.
Step 13, if the new breath value of all DG places chord is zero, judge whether to exist calculate new breath absolute value very large and estimate trend value very little close to zero branch road.If exist branch road meet new breath value absolute value comparatively large and estimate trend value close to zero condition, so Topology Error branch road is the most downstream branch road meeting above-mentioned condition, and this branch road virtual condition is off, and topology Identification terminates; If do not exist, show that, without Topology Error in network, topology Identification terminates.
Step 14, if the new breath value that there is DG place chord is non-vanishing, be then the situation of DG at Topology Error downstream branch, should the new breath value of DG place chord be changed into its opposite number, again calculate that the whole network calculates the estimation trend value that new breath value and step 9 calculate each branch road, if newly breath value is very large in the calculating of certain branch road according to step 4, but estimate that trend value is very little of zero, then the existence of this branch road disconnects Topology Error suddenly, and this branch road is actual is off state, and topology Identification terminates.
The principle of the method for the invention is as follows:
From circuit theory, any branch road
all available additional electrical potential source
and line impedance
represent, if this branch road is closed, then
; If disconnect, according to principle,displacement,
size equals this leg open voltage, and this branch current is zero.According to superposition principle, the running status in AC Ioad flow model k+1 moment equals the running status in k moment and the new breath network diagram sum in k+1 moment.
In innovation graph, the Injection Current of all nodes is zero, and the potential source in network is only the additional potential of circuit.When there is not Topology Error, the additional potential of all circuits is all zero, and all new breath amounts are also zero; And when there is unexpected disconnection Topology Error, the potential source of Topology Error branch road equals its open-circuit voltage, and new breath value is non-vanishing, and amplitude is larger near Topology Error branch road.The running status in AC Ioad flow model k+1 moment, namely estimates trend.When occurring to disconnect suddenly Topology Error, the estimation trend value of Topology Error branch road is zero.
Consider the impact of error, judge that disconnecting suddenly Topology Error should meet two conditions: one is that to calculate the absolute value of new breath amount very large; Two is estimate that the absolute value of trend is less.If certain branch road meets above-mentioned two conditions simultaneously, then this branch road is the branch road of Topology Error, and switch virtual condition should be off.In like manner known when closing suddenly and do not report Topology Error, criterion one with to disconnect suddenly and non-reporting errors is consistent, criterion two for the absolute value of this branch road estimation trend larger.Meet above-mentioned two conditions and then illustrate that the virtual condition of this switch is closed.
The DG that feeder line accesses is divided into powerful application scenario and low-power applications occasion two kinds of situations.After circuit disconnects, require that DG works in island operation state under powerful application scenario; If there is island state during low power application scenario, start island protect function, DG stops load supplying towards periphery.Thus, innovation graph approach not only can be relevant with the capacity of distributed power source to carrying out topology Identification containing the power distribution network of DG, and the relative position also by itself and Topology Error branch road affects.
Illustrate that it affects and proposes corresponding solution with 5 node 4 branch road feeder line schematic diagram shown in Fig. 2 below.Figure interior joint 4 accesses distributed power source, its exert oneself into
, negative sign represents that it sends electric energy.
Situation one: under low-power applications occasion, DG is positioned at Topology Error downstream branch
Suppose branch road
occur to disconnect Topology Error suddenly, do not consider power attenuation and DG stoppage in transit, result of calculation is as shown in table 1.
Table 1
The estimation trend of branch road 2,3,4 is 0 as seen from the table, calculates new breath value and is then respectively
,
with
if,
then there will be
with
calculating newly breath value be respectively
with
, according to two criterions of innovation graph approach identification Topology Error,
with
estimation trend be zero, therefore calculate the larger branch road of new breath value then for Topology Error branch road, now innovation graph approach identification result mistake.For solving this problem, the new breath value of DG chord being over the ground modified to its opposite number and DG by-DG, calculating new breath after obtaining revising in table, can see
the absolute value calculating new breath after revising is very large, and estimates that trend value is 0, and therefore Topology Error branch road is
.
Situation two: under high-power applications occasion, DG is positioned at Topology Error downstream branch
Same hypothesis branch road
occur to disconnect Topology Error suddenly, do not consider power attenuation and DG islet operation, result of calculation is as shown in table 2.
Table 2
Due to branch road
dG islet operation after disconnecting, for
with
continue power supply,
then powered by major network, therefore over the ground branch road 0-2,0-3,0-5 newly breath value be zero, only have the new breath value of DG place chord 0-4 to be over the ground
.Can obtaining branch road 2 from table 2, to calculate new breath absolute value very large and estimate that trend value is zero, and therefore innovation graph approach correctly picks out branch road 2 Topology Error.In this situation, the access of DG does not cause harmful effect to innovation graph approach topology Identification.
Situation three: under low-power applications occasion, DG is positioned at Topology Error branch road upstream
Suppose that branch road 4 occurs to disconnect Topology Error suddenly, think that DG does not consider under PQ control mode that DG exerts oneself change, do not consider power attenuation, result of calculation is as shown in table 3 simultaneously.
Table 3
The absolute value that known each branch road calculates newly breath from upper table is all equal is
.If
then the estimation trend of branch road 2 is 0, the display of innovation graph approach topology Identification result branch road 2, branch road 4 Topology Error; If
then the estimation trend of branch road 1 is 0, and innovation graph approach topology Identification result then shows branch road 1, branch road 4 Topology Error.Innovation graph approach topology Identification result is all inaccurate in above-mentioned two situations.For reaching the object of correct identification, when occurring that on certain supply path, the new breath order of magnitude of a more than branch road is close, and when estimating that trend value is the situation of zero, the most downstream branch road simultaneously meeting above-mentioned two conditions is Topology Error branch road.If in upper table
or
, because branch road 4 is in the most downstream of supply path and the absolute value meeting new breath value is comparatively large, estimate the criterion that trend is less, therefore innovation graph correctly judges that branch road 4 is as Topology Error branch road.
Situation four: under high-power applications occasion, DG is positioned at Topology Error branch road upstream
In this situation of same hypothesis, branch road 4 occurs to disconnect Topology Error suddenly, and the result of calculation obtained is consistent with table 3, but under high-power applications occasion, due to
, thus in upper table, the absolute value of the estimation trend of each branch road is all greater than 0, by innovation graph approach known branch road 4 Topology Error.So time DG access on innovation graph approach topology Identification result without impact.
In sum, DG accesses the impact of power distribution network on innovation graph approach topology Identification and is mainly manifested in: 1. under miniwatt occasion, Topology Error downstream branch has DG, should change the new breath value of DG branch road over the ground into its opposite number, utilize innovation graph approach identification just can obtain correct result; 2. under miniwatt occasion, there is DG Topology Error branch road upstream, criterion should be changed into meet new breath value absolute value comparatively large simultaneously, estimate trend be 0 most downstream branch road be Topology Error branch road.In other situations, DG accesses innovation graph approach topology Identification criterion without impact.
With 5 node 4 branch road feeder line schematic diagram shown in Fig. 2, the impact of bad data on Topology Error criterion is described equally, to know large-power occasions that DG accesses the criterion of distribution to innovation graph approach topology Identification and do not cause harmful effect from above-mentioned analysis, be only described for low-power applications occasion when therefore considering the affecting of bad data.
Suppose that branch road 2 disconnects suddenly, node 3 measures as bad data, from
become
.Do not consider power attenuation, result of calculation is as shown in table 4:
Table 4
Can see from upper table, branch road 1, branch road 2 difference vector are larger and equal, with node 3 over the ground chord form bad loop, therefore node 3 measures as bad data, and should revise becomes its measuring value and deduct difference vector.After getting rid of the impact of bad data, gained table content is consistent with table 1.
More than analyze known, based on innovation graph approach containing DG power distribution network topology Identification method to identification bad data and Topology Error very effective.
Claims (5)
1., based on the topology Identification method containing DG power distribution network of innovation graph approach, it is characterized in that comprising the following steps:
Step one, to Modeling of Distribution Network, determines chord and tree;
Step 2, calculates the state estimation result in k moment;
Step 3, calculates the new breath value of all measurement branch roads;
Step 4, the new breath value according to chord calculates the calculating newly breath value obtaining the whole network;
Step 5, props up the difference vector of the difference calculating its true new breath value and calculate new breath value to the tree being equipped with measurement;
Step 6, judges whether to exist partial branch difference vector very large and independently to exist, carries out step 8 if exist; Step 9 is carried out if do not exist;
Step 7, judges whether that there is the very large branch road of difference vector forms bad loop, carries out step 8 if exist; Step 9 is carried out if do not exist;
Step 8, obtains the calculating newly breath value of the whole network again according to judged result after getting rid of the impact of bad data;
Step 9, the state estimation result according to the k moment tries to achieve estimation trend with the new breath value of calculating;
Step 10, judges whether contained DG power distribution network is miniwatt occasion; If large-power occasions carries out step 11; If miniwatt occasion carries out step 12;
Step 11, the calculating judging whether to exist certain branch road newly breath value is very large but estimate that trend value is very little of zero; If the calculating that there is certain branch road newly breath value is very large, but estimates that trend value is very little of zero, so this branch road exists and disconnects Topology Error suddenly, and this branch road is actual is off state, and topology Identification terminates; If do not exist, show that network is without Topology Error, topology Identification terminates;
Step 12, judges whether that the new breath value that there is DG place chord is non-vanishing; If the new breath value of DG place chord is non-vanishing carry out step 13, if the new breath value of DG place chord is zero, carry out step 14;
Step 13, if the new breath value of all DG places chord is zero, judges whether to exist that to calculate new breath absolute value very large, and estimate trend value very little close to zero branch road; If there is branch road, to meet new breath value absolute value comparatively large, estimate trend value close to zero condition, so Topology Error branch road is the most downstream branch road meeting above-mentioned condition, and this branch road virtual condition is off, and topology Identification terminates; If do not exist, show that, without Topology Error in network, topology Identification terminates;
Step 14, if the new breath value that there is DG place chord is non-vanishing, be then the situation of DG at Topology Error downstream branch, should the new breath value of DG place chord be changed into its opposite number, again calculate that the whole network calculates the estimation trend value that new breath value and step 9 calculate each branch road, if newly breath value is very large in the calculating of certain branch road according to step 4, but estimate that trend value is very little of zero, then the existence of this branch road disconnects Topology Error suddenly, and this branch road is actual is off state, and topology Identification terminates.
2. the topology Identification method containing DG power distribution network based on innovation graph approach according to claim 1, it is characterized in that: the state estimation result calculating the k moment described in step 2 carries out state estimation calculating with the measurement amount in k moment and network topology structure, and state estimation calculates according to formula one:
formula one
In formula,
for state variable, branch current magnitudes and phase angle are chosen as its state variable for power distribution network,
for measurement amount,
for measuring function,
for measuring weight matrix;
Separate the nonlinear equation obtained with Niu Lafa, kth time iteration amount is:
Wherein
, for measuring Jacobian matrix,
for information matrix, iterative computation until
be less than setting threshold value (
), quantity of state now
required by being, and then voltage, the Power Flow Information of whole network can be obtained.
3. the topology Identification method containing DG power distribution network based on innovation graph approach according to claim 1, is characterized in that: the new breath value described in step 3 obtains according to formula two:
formula two
In formula
for k+1 moment measurement amount,
for the k moment measures estimated value,
be the new breath value in k+1 moment.
4. the topology Identification method containing DG power distribution network based on innovation graph approach according to claim 1, is characterized in that: the difference vector described in step 5 obtains according to formula three:
formula three
In formula,
for newly being ceased the calculating newly breath value calculating gained by chord,
for the true new breath value of the branch road that step 3 obtains.
5. the topology Identification method containing DG power distribution network based on innovation graph approach according to claim 1, is characterized in that: the estimation trend described in step 9 obtains according to formula four:
formula four
In formula
for the k moment measures estimated value and state estimation result.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410530548.7A CN104281749B (en) | 2014-10-10 | 2014-10-10 | The topology Identification method of the power distribution network containing DG based on innovation graph approach |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410530548.7A CN104281749B (en) | 2014-10-10 | 2014-10-10 | The topology Identification method of the power distribution network containing DG based on innovation graph approach |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104281749A true CN104281749A (en) | 2015-01-14 |
| CN104281749B CN104281749B (en) | 2017-12-08 |
Family
ID=52256619
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410530548.7A Active CN104281749B (en) | 2014-10-10 | 2014-10-10 | The topology Identification method of the power distribution network containing DG based on innovation graph approach |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104281749B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107909197A (en) * | 2017-11-13 | 2018-04-13 | 国网福建省电力有限公司 | A kind of big branch method of operation statistical analysis technique based on feeder line tree |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101593977A (en) * | 2009-04-17 | 2009-12-02 | 哈尔滨工业大学 | Method for estimating branch parameters of power network |
| CN102034026A (en) * | 2010-12-16 | 2011-04-27 | 哈尔滨工业大学 | Innovation graph approach based identification method of sudden load change of electric system |
| CN102136105A (en) * | 2011-03-22 | 2011-07-27 | 哈尔滨工业大学 | Phase measurement unit-based power grid information graph parameter estimation method |
| CN102968691A (en) * | 2012-11-20 | 2013-03-13 | 广东电网公司 | Power-cut scheme decision making system and application method thereof |
-
2014
- 2014-10-10 CN CN201410530548.7A patent/CN104281749B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101593977A (en) * | 2009-04-17 | 2009-12-02 | 哈尔滨工业大学 | Method for estimating branch parameters of power network |
| CN102034026A (en) * | 2010-12-16 | 2011-04-27 | 哈尔滨工业大学 | Innovation graph approach based identification method of sudden load change of electric system |
| CN102136105A (en) * | 2011-03-22 | 2011-07-27 | 哈尔滨工业大学 | Phase measurement unit-based power grid information graph parameter estimation method |
| CN102968691A (en) * | 2012-11-20 | 2013-03-13 | 广东电网公司 | Power-cut scheme decision making system and application method thereof |
Non-Patent Citations (3)
| Title |
|---|
| RENATO RIZZO ETAL.: "An Innovative Reconfigurable Integrated Converter Topology Suitable for Distributed Generation", 《ENERGIES》 * |
| 张艳军 等: "新息图状态估计分块算法", 《哈尔滨工程大学学报》 * |
| 林瑞星 等: "改进新息图法在地区电网AVC系统中的应用", 《电力系统自动化》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107909197A (en) * | 2017-11-13 | 2018-04-13 | 国网福建省电力有限公司 | A kind of big branch method of operation statistical analysis technique based on feeder line tree |
| CN107909197B (en) * | 2017-11-13 | 2021-08-31 | 国网福建省电力有限公司 | Statistical analysis method for operation mode of large branch based on feeder tree |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104281749B (en) | 2017-12-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Della Giustina et al. | Electrical distribution system state estimation: measurement issues and challenges | |
| Schenato et al. | Bayesian linear state estimation using smart meters and pmus measurements in distribution grids | |
| Yuan et al. | A comparative study of measurement-based Thevenin equivalents identification methods | |
| Wang et al. | A new node voltage stability index based on local voltage phasors | |
| Alzate et al. | Distribution system state estimation to support coordinated voltage-control strategies by using smart meters | |
| CN108448568A (en) | Power distribution network admixture method of estimation based on a variety of time cycle measurement data | |
| Rahnamay-Naeini et al. | On the role of power-grid and communication-system interdependencies on cascading failures | |
| CN103413044A (en) | Substation measurement information based method for estimating local topology of electric power system | |
| CN113078630A (en) | Low-voltage distribution network topology identification method based on real-time measurement data | |
| CN107679745A (en) | A kind of typical rack feeder line power supply reliability evaluation method | |
| Yuan et al. | Real-time wide-area loading margin sensitivity (WALMS) in power systems | |
| Kim et al. | LSTM-based fault direction estimation and protection coordination for networked distribution system | |
| CN105610168A (en) | Application method of node voltage stability index of multi-feed-in alternating current and direct current system | |
| CN105305507A (en) | PMU (phasor measurement Unit) distribution optimization method and system for AC-DC interconnected power grid | |
| Eghbali et al. | Multi-area state estimation based on PMU measurements in distribution networks | |
| Duong et al. | Online voltage stability monitoring based on PMU measurements and system topology | |
| Su et al. | Novel voltage-to-power sensitivity estimation for phasor measurement unit-unobservable distribution networks based on network equivalent | |
| Goh et al. | Dynamic estimation of power system stability in different Kalman filter implementations | |
| CN109858061A (en) | Distribution network equivalence and simplifying method towards voltage power sensitivity estimation | |
| CN103647284B (en) | A kind of voltage stabilization Forecasting Methodology of discontinuity surface problem when solving list | |
| CN104281749A (en) | Innovation graph method based DG (distributed generation) power distribution network included topology identification method | |
| Duong et al. | A method for real time voltage stability monitoring in sub-transmission networks | |
| CN104821577A (en) | Three-phase four-wire system distribution network robust estimation method based on intelligent electric meter measurement | |
| Polster et al. | Voltage stability monitoring using a modified thevenin impedance | |
| Gomathi et al. | Optimal location of PMUs for complete observability of power system network |
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 | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 510620 No.2, Tianhe south 2nd Road, Guangzhou, Guangdong Patentee after: GUANGZHOU POWER SUPPLY Co.,Ltd. Patentee after: China Southern Power Grid Digital Grid Technology (Guangdong) Co.,Ltd. Patentee after: GUANGZHOU JIAYUAN ELECTRIC POWER TECHNOLOGY Co.,Ltd. Address before: 510620 No.2, Tianhe south 2nd Road, Guangzhou, Guangdong Patentee before: GUANGZHOU POWER SUPPLY Co.,Ltd. Patentee before: GUANGZHOU KE TENG INFORMATION TECHNOLOGY Co.,Ltd. Patentee before: GUANGZHOU JIAYUAN ELECTRIC POWER TECHNOLOGY Co.,Ltd. |
|
| CP01 | Change in the name or title of a patent holder |