CN107238779A - A kind of fault distance-finding method of active power distribution network - Google Patents
A kind of fault distance-finding method of active power distribution network Download PDFInfo
- Publication number
- CN107238779A CN107238779A CN201710353758.7A CN201710353758A CN107238779A CN 107238779 A CN107238779 A CN 107238779A CN 201710353758 A CN201710353758 A CN 201710353758A CN 107238779 A CN107238779 A CN 107238779A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msubsup
- fault
- section
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/08—Locating faults in cables, transmission lines, or networks
- G01R31/081—Locating faults in cables, transmission lines, or networks according to type of conductors
- G01R31/086—Locating faults in cables, transmission lines, or networks according to type of conductors in power transmission or distribution networks, i.e. with interconnected conductors
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
- Y04S10/52—Outage or fault management, e.g. fault detection or location
Abstract
The invention discloses a kind of fault distance-finding method of active power distribution network, it comprises the following steps：The data of active power distribution network are obtained first, secondly each power supply section is judged to confirm fault section, then Injection Current of the trouble point downstream to trouble point is calculated using forward-backward sweep method, furthermore calculate the voltage phasor and electric current phasor of fault section downstream, and the nominal position of trouble point is calculated, trouble point is finally calculated to the distance of transformer station.The present invention considers the access of distributed power source multiple spot, need only be by using the synchronizing current phasor of local voltage electric current phasor and distributed power source with regard to the fault localization of active power distribution network can be carried out, and achieve that positioning for fault section by comparing the electric current of two sections of section, by calculating fault current and realizing the positioning distance measuring for failure to the contribution degree of fault current from trouble point downstream electric current, it is adaptable to the fault localization in the case of the access of multiple spot distributed power source.
Description
Technical field
The present invention relates to a kind of fault distance-finding method of active power distribution network, belong to distribution network automated technical field.
Background technology
Various distributed power sources constantly access power distribution network, and the permeability of distributed power source is also constantly lifted, and what this all made has
The FLT of source power distribution network turns into one of guarantee power supply reliability and the key technology of power supply quality.In regular distribution net
In, it can't detect fault current at trouble point downstream switch using fault current is able to detect that at the upstream switch of trouble point
Phenomenon realizes the positioning of fault section.And in the active power distribution network of distributed power source hypersynchronous, what distributed power source was provided
Short circuit current flow may make the fault current for flowing through trouble point downstream switch be more than overcurrent setting valve, cause conventional fault location
Method fails.
Both at home and abroad for influence and corresponding solution of the distributed power source to active power distribution network fault location research into
Fruit mainly concentrates on using impedance method, synchronous machine approximate model method or carries out event using the method for synchronizing voltage electric current phasor
Barrier positioning, the synchronized phasor localization method that is based on above-mentioned at present at switch, it is necessary to install synchronous phasor measuring device and high speed number additional
According to acquisition system, realize that fault section is determined by being in communication with each other between the phasor difference and intelligent terminal of two sections of alternative route
Position, but the distance that failure occurs can not be determined, only access one these approaches increases amount of calculation, and only in consideration system
The situation of individual distributed power source, does not account for distributed power source multiple spot access situation.
The content of the invention
For the deficiency of above-mentioned multi-objective optimization algorithm, the present invention proposes a kind of distributed power source multiple spot that is applied to and accessed
In the case of active power distribution network fault distance-finding method, it can realize the fault localization of active power distribution network.
The present invention solves its technical problem and adopted the technical scheme that：
A kind of fault distance-finding method for active power distribution network that the present invention is provided, it is characterized in that, comprise the following steps：Obtain first
The data of active power distribution network are obtained, secondly each power supply section is judged to confirm fault section, generation is pushed back before then utilizing
Method calculates Injection Current of the trouble point downstream to trouble point, furthermore calculates the voltage phasor and electric current phase of fault section downstream
Amount, and the nominal position of trouble point is calculated, trouble point is finally calculated to the distance of transformer station.
Preferably, the data of the source power distribution network include the electric parameter of topology, load and circuit, and transformer station in event
The synchronizing current phasor of voltage x current phasor and distributed power source during barrier.
The fault distance-finding method for another active power distribution network that the present invention is provided, it is characterized in that, comprise the following steps：
Step 1, the data of active power distribution network are obtained；
The data of the active power distribution network at least include data below：Electric parameter, the transformer station of topology, load and circuit
The synchronizing current phasor of voltage x current phasor and distributed power source in failure；
If failure occurs between node k and k+1, then the voltage of fault point is expressed as：
Wherein, VfFor the voltage of fault point, IfFor the electric current of fault point, Φ is three rank unit matrixs, m be failure away from
From Zk,k+1For the series impedance matrix of circuit, YkFor the shunt admittance matrix of circuit,For fault section initiating terminal k voltage,For fault section initiating terminal k electric current；
The electric current of fault point is expressed as：
Wherein,The electric current of trouble point is injected for trouble point upstream,The electric current of trouble point is injected for trouble point downstream；
Trouble point upstream Injection Current is：
The resistance for obtaining following fault point by formula (1) and formula (2) is：
Wherein,WithThe voltage of the phase of fault point i-th and jth phase, λ are represented respectivelyiAnd λjFault point is represented respectively
The state of each phase, the i.e. phase fault are then 1, are otherwise 0；
Step 2, power supply section is carried out judging to confirm fault section, ifThen IkFor failure prior fault section upstream
Hold the electric current phasor at k；
Step 3, if the ratio between electric current and downstream electric current of the fault section upstream end of power supply section η are less than or equal to reference
Step 4 is then gone to during value ε, step 2 is otherwise gone to and carries out judging next power supply section；
Step 4, Injection Current of the trouble point downstream to trouble point is calculated using forward-backward sweep method；
According to pushing back for method before in Load flow calculation, calculated using Fibonacci search algorithm and obtain trouble point downstream pair
The Injection Current of trouble point
In backward steps, the voltage at the end of kth+1 is expressed as：
Wherein,WithThe three-phase voltage vector at the end of kth+1 and kth end is represented respectively,Represent the three-phase at the end of kth+1
Current phasor, A and B are impedance matrix；
In backward steps, kth end is expressed as with the voltage at the end of kth+1, the relation of electric current phasor：
Wherein,The three-phase current vector at kth end is represented, c and d are impedance matrix；
When fault section connects distributed power source, if the synchronization measured by intelligent electronic device at the distributed power source
Electric current phasor is IIED, then the current phasor at the end of kth+1 be expressed as：
Before push away backward steps go to reference voltage terminal error less than reference value ε when, obtain
Step 5, the voltage x current phasor of fault section downstream is calculated, and failure is calculated using Fibonacci search algorithm
Q sections of downstream Injection Currents in section；
In the case of many contacts of many segmentations, the state at the end of fault section kth+1 is calculated and updated by formula (8), i.e.,
Wherein,For the three-phase current phasor at the end of kth+1,Every other section of the end of kth+1 electric current is pooled to during for failure
Phasor sum；
When being distributed formula power supply, formula (9) is rewritten as：
In formula,The fault current provided when being failure by distributed power source；
Due toIt is the function on fault section internal fault apart from m, therefore the fault section of determination is equally divided into n
Section, then every section of downstream Injection Current is expressed asThe downstream injection of each fault section is calculated using Fibonacci search algorithm
Electric current；
Step 6, if the ratio between electric current and downstream electric current of the fault section upstream end of power supply section η are less than or equal to reference
Step 7 is then gone to during value ε, the downstream Injection Current that the 4th step calculates next fault section is otherwise gone to；
Step 7, the distance that failure occurs is calculated；
Work as satisfaction
And during q≤n (11),
Then nominal position of the failure in section is：
In formula, mntSection nominal length during make it that formula (4) value is minimum, i.e. 0≤mnt≤ 1, mqAnd mq-1It is expressed as
Q sections and q-1 sections of fault section arrives the nominal range of section upstream end points；
Trouble point is to the distance of transformer station：
Wherein, diFor the length of i-th section of feeder line, dkRepresent the length of fault section feeder line.
Preferably, sentencing in the ratio between electric current and the downstream electric current of fault section upstream end of power supply section η and reference value ε
According to as follows：
η≤ε (14)
Wherein, ε ∈ [0.3,0.5] and for definite value.
Beneficial effects of the present invention are as follows：
The present invention considers the access of distributed power source multiple spot, need only be by using local voltage electric current phasor and distributed electrical
The synchronizing current phasor in source, and can be real by the electric current for comparing two sections of section with regard to that can carry out the fault localization of active power distribution network
Referring now to the positioning of fault section；It is by installing intelligent terminal IED monitorings at each distributed power source on the line
Synchronizing current phasor, by calculate fault current and from trouble point downstream electric current the contribution degree of fault current is realized for
The positioning distance measuring of failure, it is adaptable to the fault localization in the case of the access of multiple spot distributed power source.
Compared with prior art, the present invention has the following advantages that：
1st, the situation of distributed power source multiple spot access is taken into full account, the feelings without being limited to the access of distributed power source single-point
Condition；
2nd, must be active with regard to that can carry out by using the synchronizing current phasor of local voltage, electric current phasor and distributed power source
Fault Location for Distribution Network；
3rd, only fault section need to just be can determine that by comparing section upstream and downstream electric current；
4th, main website is only involved in the localization of fault after fault section is determined, fault section location is not required to main website participation
5th, the investment of equipment need not additionally be increased.
Brief description of the drawings
Fig. 1 is the method flow diagram of the embodiment of the present invention 1；
Fig. 2 is the faulty line model schematic of the embodiment of the present invention 2；
Fig. 3 is the method flow diagram of the embodiment of the present invention 2.
Embodiment
For the technical characterstic for illustrating this programme can be understood, below by embodiment, and its accompanying drawing is combined, to this hair
It is bright to be described in detail.Following disclosure provides many different embodiments or example is used for realizing the different knots of the present invention
Structure.In order to simplify disclosure of the invention, hereinafter the part and setting of specific examples are described.In addition, the present invention can be with
Repeat reference numerals and/or letter in different examples.This repetition is that for purposes of simplicity and clarity, itself is not indicated
Relation between various embodiments are discussed and/or set.It should be noted that part illustrated in the accompanying drawings is not necessarily to scale
Draw.Present invention omits the description to known assemblies and treatment technology and process to avoid being unnecessarily limiting the present invention.
In order to overcome the deficiencies in the prior art, the present invention is by calculating contribution of the distribution power on circuit to fault current
Spend to realize the fault localization of active power distribution network, it is applied to the active power distribution network failure in the case of distributed power source multiple spot is accessed
Ranging.
Embodiment 1
As shown in figure 1, a kind of fault distance-finding method of active power distribution network of the present invention, it comprises the following steps：Obtain first
The data of active power distribution network are obtained, secondly each power supply section is judged to confirm fault section, generation is pushed back before then utilizing
Method calculates Injection Current of the trouble point downstream to trouble point, furthermore calculates the voltage phasor and electric current phase of fault section downstream
Amount, and the nominal position of trouble point is calculated, trouble point is finally calculated to the distance of transformer station.
Preferably, the data of the source power distribution network include the electric parameter of topology, load and circuit, and transformer station in event
The synchronizing current phasor of voltage x current phasor and distributed power source during barrier.
Embodiment 2
As shown in Figures 2 and 3, the fault distance-finding method of another active power distribution network of the invention, it comprises the following steps：
Step 1, the data of active power distribution network are obtained；
Step 2, power supply section is carried out judging to confirm fault section, ifThen IkFor failure prior fault section upstream
Hold the electric current phasor at k；
Step 3, if the ratio between electric current and downstream electric current of the fault section upstream end of power supply section η are less than or equal to reference
Step 4 is then gone to during value ε, step 2 is otherwise gone to and carries out judging next power supply section；
Step 4, Injection Current of the trouble point downstream to trouble point is calculated using forward-backward sweep method；
Step 5, the voltage x current phasor of fault section downstream is calculated, and failure is calculated using Fibonacci search algorithm
Q sections of downstream Injection Currents in section；
Step 6, if the ratio between electric current and downstream electric current of the fault section upstream end of power supply section η are less than or equal to reference
Step 7 is then gone to during value ε, the downstream Injection Current that the 4th step calculates next fault section is otherwise gone to；
Step 7, the distance that failure occurs is calculated.
The specific implementation of embodiment 2 is as follows：
If failure occurs between node k and k+1, then the voltage of fault point can be expressed as：
Wherein, VfFor the voltage of fault point, IfFor the electric current of fault point, Φ is three rank unit matrixs, m be failure away from
From Zk,k+1For the series impedance matrix of circuit, YkFor the shunt admittance matrix of circuit,For fault section initiating terminal k voltage,For fault section initiating terminal k electric current.And the electric current of fault point can be expressed as：
Wherein,The electric current of trouble point is injected for trouble point upstream,The electric current of trouble point is injected for trouble point downstream.
The circuit model of failure is as shown in Figure 2.
Trouble point upstream Injection Current is：
The resistance that can obtain following fault point by formula (1) and formula (2) is：
WhereinWithThe voltage of the phase of fault point i-th and jth phase, λ are represented respectivelyiAnd λjRepresent that fault point is each respectively
The state of phase, the i.e. phase fault are then 1, are otherwise 0.
According to pushing back for method before in Load flow calculation, it can be calculated and be obtained under trouble point using Fibonacci search algorithm
Swim the Injection Current to trouble pointWhen physical fault occurs, all loads can be equivalent to resistance, and all distributions
Formula power supply can be equivalent to current source, therefore the STU at distributed power source can obtain its synchronizing current phasor by measurement.
During forward-backward sweep method, the voltage at the end of kth+1 is represented by：
WhereinWithThe three-phase voltage vector at the end of kth+1 and kth end is represented respectively,Represent the three-phase at the end of kth+1
Current phasor, A and B are impedance matrix.In backward steps, the relation of kth end and the voltage x current phasor at the end of kth+1 can be with table
It is shown as：
WhereinThe three-phase current vector at kth end is represented, c and d are impedance matrix.If fault section connects distributed power source,
If the synchronizing current phasor measured by intelligent electronic device (IED) at distributed power source is IIED, then the end of kth+1
Current phasor is represented by
It is current when pushing away backward steps and going to reference voltage terminal error less than ε, so that it may to obtain
In the case of many contacts of many segmentations, the state at the end of fault section kth+1 can be calculated and updated by formula (8),
I.e.
Wherein,Only to consider three-phase current phasor at this section at the end of kth+1,The end of kth+1 is pooled to during for failure
Every other section of electric current phasor sum.Formula power supply is distributed, then formula (9) can be rewritten as
In formulaThe fault current provided when being failure by distributed power source.
Due toIt is the function on fault section internal fault apart from m, therefore the fault section of determination is equally divided into n
Section, then every section of downstream Injection Current is expressed asEvery section of downstream Injection Current is calculated using Fibonacci search algorithm, works as satisfaction
When, then nominal position of the failure in section is：
M in formulantSection nominal length during make it that formula (4) value is minimum, i.e. 0≤mnt≤1 mq, mq-1It is expressed as event
Hinder the nominal range that q sections and q-1 sections of section arrives section upstream end points.Therefore distance of the failure away from transformer station can be obtained
For：
Wherein, diFor the length of i-th section of feeder line, dkRepresent the length of fault section feeder line.
In the power system of prior art, the voltage variety of supply line should not caused by after distributed power source access
The ratio between system short-circuit electric current is less than when rated current and grid entry point more than the 5% of supply line, i.e. distributed power source are short-circuit
5%, while distributed power source short circuit current flow is not more than 6 times of its rated current, therefore it is concluded that：Actual active
In power distribution network, the short circuit current flow that distributed power source is provided is not over the 30% of the short circuit current flow that system is provided.Also imply that
The ratio between electric current and downstream electric current of fault section upstream end are not more than 0.3.Measured in view of load current and current transformer
The influence of the factors such as error, certain nargin is generally left to fault location criterion, so setting identification fault section criterion such as
Under：
η≤ε (14)
Wherein ε ∈ [0.3,0.5] and for definite value, the value also can determine whether interval where fault section internal fault point.
The present invention carries out the idiographic flow of fault location ranging as shown in figure 3, the active power distribution network fault localization step is：
1) data of active power distribution network are obtained, such as topology, the electric parameter of load and circuit, transformer station are in failure
The synchronizing current phasor of voltage x current phasor and distributed power source；
2) each power supply section is judged, to confirm fault section, wherein settingAnd IkFor failure prior fault
Electric current phasor at the k of section upstream end；
If 3) meet formula (14) formula, go to the 4th step and carry out position of failure point judgement in fault section, otherwise go to the 2nd step
Judge next section；
4) Injection Current of the trouble point downstream to trouble point is calculated using forward-backward sweep method；
5) the voltage x current phasor of fault section downstream is calculated, and fault section is calculated using Fibonacci search algorithm
In q sections of downstream Injection Currents；
If 6) meet formula (14), the nominal position that the 7th step calculates trouble point is gone to, the 4th step is otherwise gone to and calculates failure
The downstream Injection Current of next section of section；
7) distance that failure occurs is calculated according to formula (12) and (13), that is, obtains trouble point to the distance of transformer station.
Synchronizing current phasor, profit are monitored by installing intelligent terminal IED at each distributed power source on the line
With the distance-finding method of the present invention by calculating fault current and the contribution degree of fault current being realized from trouble point downstream electric current
For the positioning distance measuring of failure.
The present invention has taken into full account the situation of distributed power source multiple spot access, need only be by using local voltage electric current phasor
And method of the synchronizing current phasor of distributed power source with regard to active power distribution network fault localization can be carried out, and by comparing section
Two sections of electric current achieves that the positioning for fault section.
Simply the preferred embodiment of the present invention described above, for those skilled in the art,
Without departing from the principles of the invention, some improvements and modifications can also be made, these improvements and modifications are also regarded as this hair
Bright protection domain.
Claims (4)
1. a kind of fault distance-finding method of active power distribution network, it is characterized in that, comprise the following steps：Active power distribution network is obtained first
Data, are secondly judged to confirm fault section each power supply section, are then calculated using forward-backward sweep method under trouble point
The Injection Current to trouble point is swum, furthermore calculates the voltage phasor and electric current phasor of fault section downstream, and calculates trouble point
Nominal position, finally calculate trouble point to transformer station distance.
2. a kind of fault distance-finding method of active power distribution network as claimed in claim 1, it is characterized in that, the number of the source power distribution network
According to include topology, load and circuit electric parameter, and voltage x current phasor and distributed power source of the transformer station in failure
Synchronizing current phasor.
3. a kind of fault distance-finding method of active power distribution network, it is characterized in that, comprise the following steps：
Step 1, the data of active power distribution network are obtained；
The data of the active power distribution network at least include data below：Topology, the electric parameter of load and circuit, transformer station are in event
The synchronizing current phasor of voltage x current phasor and distributed power source during barrier；
If failure occurs between node k and k+1, then the voltage of fault point is expressed as：
<mrow>
<msub>
<mi>V</mi>
<mi>f</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>&Phi;</mi>
<mo>+</mo>
<mn>0.5</mn>
<msup>
<mi>m</mi>
<mn>2</mn>
</msup>
<msub>
<mi>Z</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>Y</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<msubsup>
<mi>V</mi>
<mi>f</mi>
<mi>k</mi>
</msubsup>
<mo>-</mo>
<msub>
<mi>mZ</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<msubsup>
<mi>I</mi>
<mi>f</mi>
<mi>k</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, VfFor the voltage of fault point, IfFor the electric current of fault point, Φ is three rank unit matrixs, and m is fault distance,
Zk,k+1For the series impedance matrix of circuit, YkFor the shunt admittance matrix of circuit,For fault section initiating terminal k voltage,
For fault section initiating terminal k electric current；
The electric current of fault point is expressed as：
<mrow>
<msub>
<mi>I</mi>
<mi>f</mi>
</msub>
<mo>=</mo>
<msubsup>
<mi>I</mi>
<mi>f</mi>
<mi>U</mi>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>I</mi>
<mi>f</mi>
<mi>D</mi>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein,The electric current of trouble point is injected for trouble point upstream,The electric current of trouble point is injected for trouble point downstream；
Trouble point upstream Injection Current is：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>I</mi>
<mi>f</mi>
<mi>U</mi>
</msubsup>
<mo>=</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>mY</mi>
<mi>k</mi>
</msub>
<mo>+</mo>
<mn>0.25</mn>
<msup>
<mi>m</mi>
<mn>3</mn>
</msup>
<msub>
<mi>Y</mi>
<mi>k</mi>
</msub>
<msub>
<mi>Z</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>Y</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<msubsup>
<mi>V</mi>
<mi>f</mi>
<mi>k</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mrow>
<mi>&Phi;</mi>
<mo>+</mo>
<mn>0.5</mn>
<msup>
<mi>m</mi>
<mn>2</mn>
</msup>
<msub>
<mi>Z</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<msub>
<mi>Y</mi>
<mi>k</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<msubsup>
<mi>I</mi>
<mi>f</mi>
<mi>k</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
The resistance for obtaining following fault point by formula (1) and formula (2) is：
<mrow>
<msub>
<mi>x</mi>
<mi>f</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
<msub>
<mi>V</mi>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
</msub>
<mo>-</mo>
<msub>
<mi>&lambda;</mi>
<mi>j</mi>
</msub>
<msub>
<mi>V</mi>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
</msub>
</mrow>
<mrow>
<msub>
<mi>&lambda;</mi>
<mi>i</mi>
</msub>
<msub>
<mi>I</mi>
<msub>
<mi>f</mi>
<mi>i</mi>
</msub>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein,WithThe voltage of the phase of fault point i-th and jth phase, λ are represented respectivelyiAnd λjEach phase of fault point is represented respectively
State, i.e., the phase fault then be 1, be otherwise 0；
Step 2, power supply section is carried out judging to confirm fault section, ifThen IkFor failure prior fault section upstream end k
The electric current phasor at place；
Step 3, if the ratio between electric current and downstream electric current of the fault section upstream end of power supply section η are less than or equal to reference value ε
Shi Ze goes to step 4, otherwise goes to step 2 and carries out judging next power supply section；
Step 4, Injection Current of the trouble point downstream to trouble point is calculated using forward-backward sweep method；
According to pushing back for method before in Load flow calculation, calculated using Fibonacci search algorithm and obtain trouble point downstream to failure
The Injection Current of point
In backward steps, the voltage at the end of kth+1 is expressed as：
<mrow>
<msubsup>
<mi>V</mi>
<mi>f</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>AV</mi>
<mi>f</mi>
<mi>k</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>BI</mi>
<mi>f</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein,WithThe three-phase voltage vector at the end of kth+1 and kth end is represented respectively,Represent the three-phase current at the end of kth+1
Vector, A and B are impedance matrix；
In backward steps, kth end is expressed as with the voltage at the end of kth+1, the relation of electric current phasor：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>V</mi>
<mi>f</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>AV</mi>
<mi>f</mi>
<mi>k</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>BI</mi>
<mi>f</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>I</mi>
<mi>f</mi>
<mi>k</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mi>cV</mi>
<mi>f</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>dI</mi>
<mi>f</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein,The three-phase current vector at kth end is represented, c and d are impedance matrix；
When fault section connects distributed power source, if the synchronizing current measured by intelligent electronic device at the distributed power source
Phasor is IIED, then the current phasor at the end of kth+1 be expressed as：
<mrow>
<msubsup>
<mi>I</mi>
<mi>f</mi>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<mo>-</mo>
<msubsup>
<mi>I</mi>
<mi>f</mi>
<mrow>
<mi>I</mi>
<mi>E</mi>
<mi>D</mi>
</mrow>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
Before push away backward steps go to reference voltage terminal error less than reference value ε when, obtain
Step 5, the voltage x current phasor of fault section downstream is calculated, and fault section is calculated using Fibonacci search algorithm
In q sections of downstream Injection Currents；
In the case of many contacts of many segmentations, the state at the end of fault section kth+1 is calculated and updated by formula (8), i.e.,
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>V</mi>
<mi>f</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mover>
<mi>I</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>d</mi>
<mi>l</mi>
</msub>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>b</mi>
<mi>l</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>c</mi>
<mi>l</mi>
</msub>
</mrow>
</mtd>
<mtd>
<msub>
<mi>a</mi>
<mi>l</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>V</mi>
<mi>f</mi>
<mi>k</mi>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<msubsup>
<mi>I</mi>
<mi>f</mi>
<mi>k</mi>
</msubsup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>I</mi>
<mi>f</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<msub>
<mover>
<mi>I</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>I</mi>
<mrow>
<mi>&Sigma;</mi>
<mo>,</mo>
<mi>f</mi>
</mrow>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein,For the three-phase current phasor at the end of kth+1,Every other section of the end of kth+1 electric current phasor is pooled to during for failure
Sum；
When being distributed formula power supply, formula (9) is rewritten as：
<mrow>
<msubsup>
<mi>I</mi>
<mi>f</mi>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>=</mo>
<msub>
<mover>
<mi>I</mi>
<mo>&OverBar;</mo>
</mover>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msubsup>
<mi>I</mi>
<mrow>
<mi>&Sigma;</mi>
<mo>,</mo>
<mi>f</mi>
</mrow>
<mrow>
<mi>k</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>I</mi>
<mi>f</mi>
<mrow>
<mi>D</mi>
<mi>G</mi>
</mrow>
</msubsup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula,The fault current provided when being failure by distributed power source；
Due toIt is the function on fault section internal fault apart from m, therefore the fault section of determination is equally divided into n sections, then often
Section downstream Injection Current is expressed asThe downstream Injection Current of each fault section is calculated using Fibonacci search algorithm；
Step 6, if the ratio between electric current and downstream electric current of the fault section upstream end of power supply section η are less than or equal to reference value ε
Shi Ze goes to step 7, otherwise goes to the downstream Injection Current that the 4th step calculates next fault section；
Step 7, the distance that failure occurs is calculated；
Work as satisfaction
And q≤n (11)
When, then nominal position of the failure in section is：
<mrow>
<msup>
<mi>m</mi>
<mrow>
<mi>n</mi>
<mi>t</mi>
</mrow>
</msup>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>m</mi>
<mrow>
<mi>q</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>m</mi>
<mi>q</mi>
</msub>
</mrow>
<mn>2</mn>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
In formula, mntSection nominal length during make it that formula (4) value is minimum, i.e. 0≤mnt≤ 1, mqAnd mq-1It is expressed as failure
Q sections and q-1 sections of section arrives the nominal range of section upstream end points；
Trouble point is to the distance of transformer station：
<mrow>
<mi>L</mi>
<mo>=</mo>
<munderover>
<mi>&Sigma;</mi>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<msub>
<mi>d</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<msup>
<mi>m</mi>
<mrow>
<mi>n</mi>
<mi>t</mi>
</mrow>
</msup>
<msub>
<mi>d</mi>
<mi>k</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein, diFor the length of i-th section of feeder line, dkRepresent the length of fault section feeder line.
4. a kind of fault distance-finding method of active power distribution network as claimed in claim 3, it is characterized in that, in the failure of power supply section
The ratio between the electric current of section upstream end and downstream electric current η and reference value ε criterion are as follows：
η≤ε (14)
Wherein, ε ∈ [0.3,0.5] and for definite value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710353758.7A CN107238779B (en) | 2017-05-18 | 2017-05-18 | A kind of fault distance-finding method of active power distribution network |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710353758.7A CN107238779B (en) | 2017-05-18 | 2017-05-18 | A kind of fault distance-finding method of active power distribution network |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107238779A true CN107238779A (en) | 2017-10-10 |
| CN107238779B CN107238779B (en) | 2019-10-11 |
Family
ID=59984515
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710353758.7A Active CN107238779B (en) | 2017-05-18 | 2017-05-18 | A kind of fault distance-finding method of active power distribution network |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107238779B (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107843810A (en) * | 2017-11-01 | 2018-03-27 | 东南大学 | A kind of active power distribution network fault section tuning on-line method based on state estimation |
| CN108387818A (en) * | 2018-01-23 | 2018-08-10 | 中国石油大学(华东) | A kind of fault distance-finding method suitable for the tree-shaped catalogue containing distributed generation resource |
| CN109375055A (en) * | 2018-09-28 | 2019-02-22 | 国网河北省电力有限公司电力科学研究院 | A kind of detection method of phasor measurement system, system, device and storage medium |
| CN111527661A (en) * | 2017-12-29 | 2020-08-11 | Abb电网瑞士股份公司 | Fault location in multi-terminal tapped lines |
| CN111525441A (en) * | 2020-04-30 | 2020-08-11 | 国网湖南省电力有限公司 | Uninterrupted operation method for operation point |
| CN113848431A (en) * | 2021-10-22 | 2021-12-28 | 广西电网有限责任公司电力科学研究院 | Power distribution network line fault positioning method and system |
| CN115308536A (en) * | 2022-09-29 | 2022-11-08 | 西华大学 | Mu PMU-based DG-containing power distribution network fault section identification method |
| CN116626540A (en) * | 2023-07-25 | 2023-08-22 | 国网山东省电力公司莱芜供电公司 | Method, system, terminal and storage medium for judging broken line fault interval |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102253307A (en) * | 2011-04-19 | 2011-11-23 | 中国电力科学研究院 | Method for positioning single-phase grounding failure of distribution network containing distributed power supply |
| CN103208818A (en) * | 2013-04-16 | 2013-07-17 | 山东科技大学 | Distribution network loop closing current calculating method based on distributed power supply |
| CN103344882A (en) * | 2013-07-02 | 2013-10-09 | 国家电网公司 | Distributed power source contained power distribution network fault location method and system |
| CN104459465A (en) * | 2014-11-26 | 2015-03-25 | 国家电网公司 | Fault section locating method for power distribution network containing distributed power supplies |
| US9389268B1 (en) * | 2013-01-28 | 2016-07-12 | The Florida State University Research Foundation, Inc. | System for locating a fault in a power distribution system comprising at least one DC voltage distributed resource having a controllable voltage source converter |
| CN105759173A (en) * | 2016-03-31 | 2016-07-13 | 中国科学院电工研究所 | Adaptive fault section positioning method of power distribution network containing distributed power supply |
| CN106383296A (en) * | 2016-11-02 | 2017-02-08 | 华北电力大学（保定） | Improved impedance type active distribution network fault location algorithm based on phasor analysis |
-
2017
- 2017-05-18 CN CN201710353758.7A patent/CN107238779B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102253307A (en) * | 2011-04-19 | 2011-11-23 | 中国电力科学研究院 | Method for positioning single-phase grounding failure of distribution network containing distributed power supply |
| US9389268B1 (en) * | 2013-01-28 | 2016-07-12 | The Florida State University Research Foundation, Inc. | System for locating a fault in a power distribution system comprising at least one DC voltage distributed resource having a controllable voltage source converter |
| CN103208818A (en) * | 2013-04-16 | 2013-07-17 | 山东科技大学 | Distribution network loop closing current calculating method based on distributed power supply |
| CN103344882A (en) * | 2013-07-02 | 2013-10-09 | 国家电网公司 | Distributed power source contained power distribution network fault location method and system |
| CN104459465A (en) * | 2014-11-26 | 2015-03-25 | 国家电网公司 | Fault section locating method for power distribution network containing distributed power supplies |
| CN105759173A (en) * | 2016-03-31 | 2016-07-13 | 中国科学院电工研究所 | Adaptive fault section positioning method of power distribution network containing distributed power supply |
| CN106383296A (en) * | 2016-11-02 | 2017-02-08 | 华北电力大学（保定） | Improved impedance type active distribution network fault location algorithm based on phasor analysis |
Non-Patent Citations (1)
| Title |
|---|
| 戴志辉 等: "基于改进阻抗法的有源配电网故障测距算法", 《POWER SYSTEM TECHNOLOGY》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107843810A (en) * | 2017-11-01 | 2018-03-27 | 东南大学 | A kind of active power distribution network fault section tuning on-line method based on state estimation |
| CN107843810B (en) * | 2017-11-01 | 2019-12-17 | 东南大学 | Active power distribution network fault section online positioning method based on state estimation |
| CN111527661A (en) * | 2017-12-29 | 2020-08-11 | Abb电网瑞士股份公司 | Fault location in multi-terminal tapped lines |
| CN108387818A (en) * | 2018-01-23 | 2018-08-10 | 中国石油大学(华东) | A kind of fault distance-finding method suitable for the tree-shaped catalogue containing distributed generation resource |
| CN109375055A (en) * | 2018-09-28 | 2019-02-22 | 国网河北省电力有限公司电力科学研究院 | A kind of detection method of phasor measurement system, system, device and storage medium |
| CN109375055B (en) * | 2018-09-28 | 2020-10-09 | 国网河北省电力有限公司电力科学研究院 | Detection method, system and device of phasor measurement system and storage medium |
| CN111525441A (en) * | 2020-04-30 | 2020-08-11 | 国网湖南省电力有限公司 | Uninterrupted operation method for operation point |
| CN113848431A (en) * | 2021-10-22 | 2021-12-28 | 广西电网有限责任公司电力科学研究院 | Power distribution network line fault positioning method and system |
| CN115308536A (en) * | 2022-09-29 | 2022-11-08 | 西华大学 | Mu PMU-based DG-containing power distribution network fault section identification method |
| CN116626540A (en) * | 2023-07-25 | 2023-08-22 | 国网山东省电力公司莱芜供电公司 | Method, system, terminal and storage medium for judging broken line fault interval |
| CN116626540B (en) * | 2023-07-25 | 2023-10-13 | 国网山东省电力公司莱芜供电公司 | Method, system, terminal and storage medium for judging broken line fault interval |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107238779B (en) | 2019-10-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107238779A (en) | A kind of fault distance-finding method of active power distribution network | |
| CN107340455B (en) | Single-phase broken line fault identification method and application of high-voltage circuit of power distribution network | |
| Liao et al. | Online optimal transmission line parameter estimation for relaying applications | |
| Liao et al. | Fault-location algorithms without utilizing line parameters based on the distributed parameter line model | |
| Liao | Fault location for single-circuit line based on bus-impedance matrix utilizing voltage measurements | |
| Du et al. | On-line estimation of transmission line parameters, temperature and sag using PMU measurements | |
| CN111796164B (en) | Fault location method for full-parallel AT traction network | |
| CN107064736B (en) | A kind of Fault Locating Method connecing inverse distributed power power distribution network containing more T | |
| CN109494696B (en) | Power distribution network asymmetric fault positioning and isolating method and system based on adaptive reclosing | |
| CN105762777B (en) | A kind of longitudinal protection method that inverse distributed power power distribution network is connect containing more T | |
| Cordova et al. | Fault location identification in smart distribution networks with distributed generation | |
| CN107271841B (en) | Same-tower double-circuit T-shaped power transmission line fault location method based on positive sequence current difference | |
| Majumdar et al. | A three-phase state estimation in unbalanced distribution networks with switch modelling | |
| CN111542980A (en) | Parametric-free identification of fault location in multi-terminal power transmission line | |
| CN105429137A (en) | Distribution network loop closing power supply switching method | |
| CN110927519A (en) | Active power distribution network fault positioning method based on mu PMU measurement value | |
| Duong et al. | Online voltage stability monitoring based on PMU measurements and system topology | |
| CN113972634B (en) | Micro-grid differential protection method for self-adaptive zero-sequence energy injection | |
| CN103412190B (en) | Switch-class device state evaluation method based on parameter on-line identification | |
| CN107046285A (en) | A kind of state of electric distribution network appraisal procedure based on hybrid measurement | |
| CN109586254B (en) | Voltage time type feeder line reverse closing locking method for half-cycle voltage transient quantity | |
| CN107565549B (en) | Power system network topology analysis method based on synchronous phasor measurement | |
| CN113659548B (en) | Power distribution network pilot protection method and system based on positive sequence fault component energy direction | |
| Katic et al. | Distribution state estimation field testing | |
| WO2021114804A1 (en) | Small hydropower station rapid grid connection method considering deviation and trend changes |
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 |