CN109342885A - A kind of localization method and system of DC distribution net line fault - Google Patents
A kind of localization method and system of DC distribution net line fault Download PDFInfo
- Publication number
- CN109342885A CN109342885A CN201811353661.7A CN201811353661A CN109342885A CN 109342885 A CN109342885 A CN 109342885A CN 201811353661 A CN201811353661 A CN 201811353661A CN 109342885 A CN109342885 A CN 109342885A
- Authority
- CN
- China
- Prior art keywords
- mathematical model
- line
- faulty line
- value
- fitness function
- 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
-
- 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/088—Aspects of digital computing
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Locating Faults (AREA)
Abstract
The invention discloses a kind of localization method of DC distribution net line fault, the localization method includes: firstly, establishing the mathematical model of faulty line；Then, fitness function is constructed using the mathematical model of the faulty line；Finally, identifying based on genetic algorithm to the fitness function, the abort situation of the faulty line is determined.The present invention eliminates the influence of transition resistance by the mathematical model of built faulty line, improve the accuracy of positioning, and construct fitness function, parameter identification problem is converted by fault-location problem, and then parameter identification is carried out using genetic algorithm, it determines abort situation, avoids as the wrong interference caused by positioning of data acquisition, further improve the accuracy of positioning.
Description
Technical field
The present invention relates to DC distribution net field, in particular to the localization method of a kind of DC distribution net line fault and it is
System.
Background technique
With a large amount of renewable energy access power grid, based on voltage-source type (Voltage Sourced Converters,
VSC) the flexible direct current power distribution network of inverter becomes research hotspot, and when abnormal and failure occurs in direct current system, failure process is rapid,
It is very harmful, higher requirement is proposed to fault location.It is positioned for DC system fault, can be divided mainly into traveling wave method and non-
Traveling wave method two major classes.Carrying out fault location principle using traveling wave is by detecting transient state travelling wave between fault point and measurement point
Transmission time difference realizes fault location, and this method positioning accuracy is high, strong robustness.But this method manually needs in implementation process
The amount of asking is larger, it is difficult to automation is realized, and it is more demanding to sample frequency, when wavefront amplitude is restricted, it will
Positioning is caused to fail.And then there is scholar to propose that non-traveling wave method carries out fault location, the non-traveling wave method master of extensive discussions is obtained at present
Have: a method of fault location is realized using single-end electrical quantity, this method positioning accuracy is higher, and calculation amount is small, but vulnerable to
The influence of opposite end inverter adjustment effect, the accuracy applied in practical projects are to be improved.Therefore, direct current system how is improved
The anti-transition resistance ability for fault location of uniting is strong, and then improves the accuracy of fault location, becomes a technology urgently to be resolved
Problem.
Summary of the invention
The object of the present invention is to provide the localization methods and system of a kind of DC distribution net line fault, to improve direct current system
The anti-transition resistance ability for fault location of uniting is strong, and then improves the accuracy of fault location.
To achieve the above object, the present invention provides following schemes:
A kind of localization method of DC distribution net line fault, the localization method include the following steps:
Establish the mathematical model of faulty line；
Fitness function is constructed using the mathematical model of the faulty line；
The fitness function is identified based on genetic algorithm, determines the abort situation of the faulty line.
Optionally, the mathematical model for establishing faulty line, specifically includes:
Obtain the first equivalent circuit of the faulty line；
According to first equivalent circuit, the first starting point mathematical model of the faulty line, the first are established
Beginning mathematical model is shown below:
According to first equivalent circuit, the first end mathematical model of the faulty line, the first end are established
Mathematical model is shown below:
By first starting point mathematical model and the first end mathematical model simultaneous, the faulty line is obtained
First mathematical model, first mathematical model are shown below:
Wherein, I1, R1, L1, C1, Vdc1Respectively close to the equivalent current values of the starting point of faulty line, resistance value, inductance
The voltage value of value, capacitance and capacitor both ends；I2, R2, L2, C2, Vdc2Respectively close to the equivalent current of the end of faulty line
Value, resistance value, inductance value, capacitance capacitor both ends voltage value, RfFor transition resistance value；If unit length line resistance is r,
Then R1=r × D1, R2=r × D2If unit length line inductance is l, then L1=l × D1, L2=l × D2, D1For close to starting point
Route length, D1For by the length of the route of proximal end, D=D1+D2, D is total line length.
Optionally, the mathematical model using the faulty line constructs fitness function, specifically includes:
By the first mathematical model discretization of the faulty line, the first discrete model, first discrete model are obtained
It is shown below:
Wherein,K indicates the number of iterations；
First discrete model is optimized using GA algorithm, obtains the first fitness function, described first adapts to
Degree function is shown below:
In formula, S (R1,R2,L1,L2) it is the first fitness function, fk(R1,R2,L1,L2) indicate kth time the first mathematics of iteration
The functional value of model.
Optionally, the mathematical model for establishing faulty line, specifically includes:
Obtain the second equivalent circuit of the faulty line；
According to second equivalent circuit, the second starting point mathematical model of the faulty line is established, described second
Beginning mathematical model is shown below:
According to second equivalent circuit, the second end mathematical model of the faulty line, the second end are established
Mathematical model is shown below:
By second starting point mathematical model and the second end mathematical model simultaneous, the faulty line is obtained
Second mathematical model, second mathematical model are shown below:
Wherein, I1, R1, L1, C1, Vdc1Respectively close to the equivalent current values of the starting point of faulty line, resistance value, inductance
The voltage value of value, capacitance and capacitor both ends；I2, R2, L2, C2, Vdc2Respectively close to the equivalent current of the end of faulty line
Value, resistance value, inductance value, capacitance capacitor both ends voltage value, RfFor transition resistance value；If unit length line resistance is r,
Then R1=r × D1, R2=r × D2If unit length line inductance is l, then L1=l × D1, L2=l × D2, D1For close to starting point
Route length, D1For by the length of the route of proximal end, D=D1+D2, D is total line length.
Optionally, the mathematical model using the faulty line constructs fitness function, specifically includes:
By the second mathematical model discretization of the faulty line, the second discrete model, second discrete model are obtained
It is shown below:
Wherein,K indicates the number of iterations；
Second discrete model is optimized using GA algorithm, obtains the second fitness function, described second adapts to
Degree function is shown below:
In formula, S'(R1,R2,L1,L2) it is the second fitness function, fk(R1,R2,L1,L2) indicate the kth time number of iteration second
Learn the functional value of model.
Optionally, described that the fitness function is identified based on genetic algorithm, determine the event of the faulty line
Hinder position, specifically include:
Calculate the fitness function value of each discrete location；
The discrete location that fitness function value is greater than preset first threshold value is chosen, as candidate fault locations；
Intersection and mutation operation are carried out to the candidate fault locations, the abort situation after being optimized；
The fitness function value of abort situation after calculating each optimization, the number of iterations increase by 1,；
The maximum position of fitness function value and last iteration obtain suitable in abort situation after judging the optimization
Whether the difference of the maximum position of response functional value is less than second threshold, obtains the first judging result；
If first judging result be it is yes, using the maximum position of the fitness function value as abort situation；
If first judging result be it is no, judge whether the number of iterations is less than maximum number of iterations, obtain the
Two judging results；
If second judging result be it is yes, return step " chooses fitness function value and is greater than preset first threshold value
Discrete location, as candidate fault locations "；
If second judging result be it is no, using the maximum position of the fitness function value as abort situation.
A kind of positioning system of DC distribution net line fault, the positioning system include:
Mathematical model establishes module, for establishing the mathematical model of faulty line；
Fitness function constructing module, for constructing fitness function using the mathematical model of the faulty line；
Abort situation determining module determines the event for identifying based on genetic algorithm to the fitness function
Hinder the abort situation of route.
Optionally, the mathematical model establishes module, specifically includes:
First equivalent circuit acquisition submodule, for obtaining the first equivalent circuit of the faulty line；
First starting point mathematical model setting up submodule, for establishing the fault wire according to first equivalent circuit
The first starting point mathematical model on road, first starting point mathematical model are shown below:
First end mathematical model setting up submodule, for establishing the faulty line according to first equivalent circuit
First end mathematical model, the first end mathematical model is shown below:
First mathematical model setting up submodule is used for first starting point mathematical model and the first end mathematics
Model simultaneous obtains the first mathematical model of the faulty line, and first mathematical model is shown below:
Wherein, I1, R1, L1, C1, Vdc1Respectively close to the equivalent current values of the starting point of faulty line, resistance value, inductance
The voltage value of value, capacitance and capacitor both ends；I2, R2, L2, C2, Vdc2Respectively close to the equivalent current of the end of faulty line
Value, resistance value, inductance value, capacitance capacitor both ends voltage value, RfFor transition resistance value；If unit length line resistance is r,
Then R1=r × D1, R2=r × D2If unit length line inductance is l, then L1=l × D1, L2=l × D2, D1For close to starting point
Route length, D1For by the length of the route of proximal end, D=D1+D2, D is total line length.
Optionally, the fitness function constructing module, specifically includes:
First mathematical model discretization submodule, for obtaining the first mathematical model discretization of the faulty line
First discrete model, first discrete model are shown below:
Wherein,K indicates the number of iterations；
First discrete model optimizes submodule, for being optimized using GA algorithm to first discrete model, obtains
First fitness function, first fitness function are shown below；
In formula, S (R1,R2,L1,L2) it is the first fitness function, fk(R1,R2,L1,L2) indicate kth time the first mathematics of iteration
The functional value of model.
Optionally, the mathematical model establishes module, specifically includes:
Second equivalent circuit acquisition submodule, for obtaining the second equivalent circuit of the faulty line；
Second starting point mathematical model setting up submodule, for establishing the fault wire according to second equivalent circuit
The second starting point mathematical model on road, second starting point mathematical model are shown below:
Second end mathematical model setting up submodule, for establishing the faulty line according to second equivalent circuit
Second end mathematical model, the second end mathematical model is shown below:
Second mathematical model setting up submodule is used for second starting point mathematical model and the second end mathematics
Model simultaneous obtains the second mathematical model of the faulty line, and second mathematical model is shown below:
Wherein, I1, R1, L1, C1, Vdc1Respectively close to the equivalent current values of the starting point of faulty line, resistance value, inductance
The voltage value of value, capacitance and capacitor both ends；I2, R2, L2, C2, Vdc2Respectively close to the equivalent current of the end of faulty line
Value, resistance value, inductance value, capacitance capacitor both ends voltage value, RfFor transition resistance value；If unit length line resistance is r,
Then R1=r × D1, R2=r × D2If unit length line inductance is l, then L1=l × D1, L2=l × D2, D1For close to starting point
Route length, D1For by the length of the route of proximal end, D=D1+D2, D is total line length.
The specific embodiment provided according to the present invention, the invention discloses following technical effects:
The invention discloses a kind of localization method of DC distribution net line fault, the localization method includes: firstly, building
The mathematical model of vertical faulty line；Then, fitness function is constructed using the mathematical model of the faulty line；Finally, being based on
Genetic algorithm identifies the fitness function, determines the abort situation of the faulty line.The present invention, which passes through, to be built
The mathematical model of faulty line eliminates the influence of transition resistance, improves the accuracy of positioning, and constructs fitness function, will
Fault-location problem is converted into parameter identification problem, and then carries out parameter identification using genetic algorithm, determines abort situation, avoids
It is interfered caused by positioning since data acquisition is wrong, further improves the accuracy of positioning.
Detailed description of the invention
It in order to more clearly explain the embodiment of the invention or the technical proposal in the existing technology, below will be to institute in embodiment
Attached drawing to be used is needed to be briefly described, it should be apparent that, the accompanying drawings in the following description is only some implementations of the invention
Example, for those of ordinary skill in the art, without any creative labor, can also be according to these attached drawings
Obtain other attached drawings.
Fig. 1 is a kind of flow chart of the localization method of DC distribution net line fault provided by the invention；
Fig. 2 is the circuit diagram of the first equivalent circuit provided by the invention；
Fig. 3 is the circuit diagram of the second equivalent circuit provided by the invention；
Fig. 4 is the flow chart provided by the invention identified based on genetic algorithm to the fitness function；
Fig. 5 is a kind of structure chart of the positioning system of DC distribution net line fault provided by the invention；
Fig. 6 is the experimental system figure of the present invention provided by the invention localization method built and system；
Fig. 7 is the simulated line figure of the localization method of the invention built and system；
Fig. 8 is the intereelectrode short-circuit fault location provided of the invention and the convergence rate experiment knot that monopolar grounding fault positions
Fruit figure；
Fig. 9 is intereelectrode short-circuit fault location provided by the invention and the accuracy experimental result that monopolar grounding fault positions
Figure；
Figure 10 is the experimental result that localization method and system of the invention provided by the invention position 18 kinds of line faults
Figure.
Specific embodiment
The object of the present invention is to provide the localization methods and system of a kind of DC distribution net line fault, to improve direct current system
The anti-transition resistance ability for fault location of uniting is strong, and then improves the accuracy of fault location.
In order to make the foregoing objectives, features and advantages of the present invention clearer and more comprehensible, with reference to the accompanying drawing and specific real
Mode is applied to be described in further detail invention.
As shown in Figure 1, a kind of localization method of DC distribution net line fault, the localization method include the following steps:
Step 101, the mathematical model of faulty line is established.
When the failure of the faulty circuit is intereelectrode short-circuit failure, the mathematical modulo of faulty line is established described in step 101
Type specifically includes:
When intereelectrode short-circuit failure, due to the self-protection function of IGBT, the locking of instantaneous current conversion device occurs for failure, and exchange can be ignored
The electric current of side influences, and the electric current on route all derives from the electric discharge of capacitor.The first equivalent circuit of the faulty line is obtained,
As shown in Figure 2.
According to first equivalent circuit, the first starting point mathematical model of the faulty line, the first are established
Beginning mathematical model is shown below:
According to first equivalent circuit, the first end mathematical model of the faulty line, the first end are established
Mathematical model is shown below:
By first starting point mathematical model and the first end mathematical model simultaneous, the faulty line is obtained
First mathematical model, first mathematical model are shown below:
Wherein, I1, R1, L1, C1, Vdc1Respectively close to the equivalent current values of the starting point of faulty line, resistance value, inductance
The voltage value of value, capacitance and capacitor both ends；I2, R2, L2, C2, Vdc2Respectively close to the equivalent current of the end of faulty line
Value, resistance value, inductance value, capacitance capacitor both ends voltage value, RfFor transition resistance value；If unit length line resistance is r,
Then R1=r × D1, R2=r × D2If unit length line inductance is l, then L1=l × D1, L2=l × D2, D1For close to starting point
Route length, D1For by the length of the route of proximal end, D=D1+D2, D is total line length.
When the failure of the faulty line is monopolar grounding fault, the mathematical modulo of faulty line is established described in step 101
Type specifically includes:
Obtain the second equivalent circuit of the faulty line；At this point, second equivalent circuit is as shown in Figure 3.
According to second equivalent circuit, the second starting point mathematical model of the faulty line is established, described second
Beginning mathematical model is shown below:
According to second equivalent circuit, the second end mathematical model of the faulty line, the second end are established
Mathematical model is shown below:
By second starting point mathematical model and the second end mathematical model simultaneous, the faulty line is obtained
Second mathematical model, second mathematical model are shown below:
Wherein, I1, R1, L1, C1, Vdc1Respectively close to the equivalent current values of the starting point of faulty line, resistance value, inductance
The voltage value of value, capacitance and capacitor both ends；I2, R2, L2, C2, Vdc2Respectively close to the equivalent current of the end of faulty line
Value, resistance value, inductance value, capacitance capacitor both ends voltage value, RfFor transition resistance value；If unit length line resistance is r,
Then R1=r × D1, R2=r × D2If unit length line inductance is l, then L1=l × D1, L2=l × D2, D1For close to starting point
Route length, D1For by the length of the route of proximal end, D=D1+D2, D is total line length.
Step 102, fitness function is constructed using the mathematical model of the faulty line.
When the failure of the faulty circuit is intereelectrode short-circuit failure, step 102, the number using the faulty line
Construction of A Model fitness function is learned, is specifically included:
By the first mathematical model discretization of the faulty line, the first discrete model, first discrete model are obtained
It is shown below:
Wherein,K indicates the number of iterations；
First discrete model is optimized using GA algorithm, obtains the first fitness function, described first adapts to
Degree function is shown below:
In formula, S (R1,R2,L1,L2) it is the first fitness function, fk(R1,R2,L1,L2) indicate kth time the first mathematics of iteration
The functional value of model.
When the failure of the faulty line is monopolar grounding fault, the number of the faulty line is utilized described in step 102
Construction of A Model fitness function is learned, is specifically included:
By the second mathematical model discretization of the faulty line, the second discrete model, second discrete model are obtained
It is shown below:
Wherein,K indicates the number of iterations；
Second discrete model is optimized using GA algorithm, obtains the second fitness function, described second adapts to
Degree function is shown below:
In formula, S'(R1,R2,L1,L2) it is the second fitness function, fk(R1,R2,L1,L2) indicate the kth time number of iteration second
Learn the functional value of model.
Step 103, the fitness function is identified based on genetic algorithm, determines the fault bit of the faulty line
It sets.
The fitness function is identified based on genetic algorithm described in step 103, determines the event of the faulty line
Hinder position, as shown in figure 4, specifically including:
Two lateral capacitance both end voltage of faulty line is sampled, parameters are initialized, is counted according to formula (10) or (7)
Calculate the fitness function value of each discrete location；Wherein when line fault is intereelectrode short-circuit failure, public affairs are sought using formula (8)
Parameter in formula (10) seeks the parameter in formula (7) using formula (9) when line fault is monopolar grounding fault.
The discrete location that fitness function value is greater than preset first threshold value is chosen, as candidate fault locations；It specifically includes
Discrete location is ranked up and is selected according to the size of fitness function value.
Intersection and mutation operation are carried out to the candidate fault locations, the abort situation after being optimized；
The fitness function value of abort situation after calculating each optimization, the number of iterations increase by 1,；
The maximum position of fitness function value and last iteration obtain suitable in abort situation after judging the optimization
Whether the difference of the maximum position of response functional value is less than second threshold, obtains the first judging result；
If first judging result be it is yes, using the maximum position of the fitness function value as abort situation；
If first judging result be it is no, judge whether the number of iterations is less than maximum number of iterations, obtain the
Two judging results；
If second judging result be it is yes, return step " chooses fitness function value and is greater than preset first threshold value
Discrete location, as candidate fault locations "；
If second judging result be it is no, using the maximum position of the fitness function value as abort situation.
As shown in figure 5, the present invention also provides a kind of positioning system of DC distribution net line fault, the positioning system
Include:
Mathematical model establishes module 501, for establishing the mathematical model of faulty line.
When the failure of the faulty circuit is intereelectrode short-circuit failure, the mathematical model establishes module 501, specific to wrap
It includes:
First equivalent circuit acquisition submodule, for obtaining the first equivalent circuit of the faulty line；
First starting point mathematical model setting up submodule, for establishing the fault wire according to first equivalent circuit
The first starting point mathematical model on road, first starting point mathematical model are shown below:
First end mathematical model setting up submodule, for establishing the faulty line according to first equivalent circuit
First end mathematical model, the first end mathematical model is shown below:
First mathematical model setting up submodule is used for first starting point mathematical model and the first end mathematics
Model simultaneous obtains the first mathematical model of the faulty line, and first mathematical model is shown below:
Wherein, I1, R1, L1, C1, Vdc1Respectively close to the equivalent current values of the starting point of faulty line, resistance value, inductance
The voltage value of value, capacitance and capacitor both ends；I2, R2, L2, C2, Vdc2Respectively close to the equivalent current of the end of faulty line
Value, resistance value, inductance value, capacitance capacitor both ends voltage value, RfFor transition resistance value；If unit length line resistance is r,
Then R1=r × D1, R2=r × D2If unit length line inductance is l, then L1=l × D1, L2=l × D2, D1For close to starting point
Route length, D1For by the length of the route of proximal end, D=D1+D2, D is total line length.
When the failure of the faulty line is monopolar grounding fault, the mathematical model establishes module 401, specific to wrap
It includes:
Second equivalent circuit acquisition submodule, for obtaining the second equivalent circuit of the faulty line；
Second starting point mathematical model setting up submodule, for establishing the fault wire according to second equivalent circuit
The second starting point mathematical model on road, second starting point mathematical model are shown below:
Second end mathematical model setting up submodule, for establishing the faulty line according to second equivalent circuit
Second end mathematical model, the second end mathematical model is shown below:
Second mathematical model setting up submodule is used for second starting point mathematical model and the second end mathematics
Model simultaneous obtains the second mathematical model of the faulty line, and second mathematical model is shown below:
Wherein, I1, R1, L1, C1, Vdc1Respectively close to the equivalent current values of the starting point of faulty line, resistance value, inductance
The voltage value of value, capacitance and capacitor both ends；I2, R2, L2, C2, Vdc2Respectively close to the equivalent current of the end of faulty line
Value, resistance value, inductance value, capacitance capacitor both ends voltage value, RfFor transition resistance value；If unit length line resistance is r,
Then R1=r × D1, R2=r × D2If unit length line inductance is l, then L1=l × D1, L2=l × D2, D1For close to starting point
Route length, D1For by the length of the route of proximal end, D=D1+D2, D is total line length.
Fitness function constructing module 502, for constructing fitness function using the mathematical model of the faulty line；When
When the failure of the faulty circuit is intereelectrode short-circuit failure, the fitness function constructing module 502 is specifically included:
First mathematical model discretization submodule, for obtaining the first mathematical model discretization of the faulty line
First discrete model, first discrete model are shown below:
Wherein,K indicates the number of iterations；
First discrete model optimizes submodule, for being optimized using GA algorithm to first discrete model, obtains
First fitness function, first fitness function are shown below；
In formula, S (R1,R2,L1,L2) it is the first fitness function, fk(R1,R2,L1,L2) indicate kth time the first mathematics of iteration
The functional value of model.
Abort situation determining module 503, for being identified based on genetic algorithm to the fitness function, described in determination
The abort situation of faulty line.
Embodiment 1:
The present invention has built 6 end ring net HILS (HWIL simulation) experimental system as shown in FIG. 6.The system is by model
The RT-LAB real-time simulator of OP5600, the DSP operation device of model TMS320F28335 and host computer etc. are constituted, and are built
Simulated line is as shown in fig. 7, fault setting on the route between G-VSC and W-VSC.
G-VSC uses double circle controling mode, and outer ring uses the sagging control of voltage, and inner ring uses constant DC current control side
Formula.W-VSC is operated under maximum power control mode, in certain situations it is desirable to downrating.Energy-storage module is in charging
Or discharge condition.Meanwhile power-balance and stable operation in order to ensure system, energy-storage units play the part of entire DC distribution net
Balance nodes role, when necessary be in isolated operation mode.Solar panel is incorporated to DC distribution by DC-DC converter
Net, outer ring use maximal power tracing control mode, and inner ring uses constant voltage control mode.DC bus-bar voltage 500V is obtained,
Current transformer DC bus capacitor 2mF, unit length line resistance value are 0.0139 Ω/km, and unit length line electricity inductance value is 0.159
Ω/km, line length 10km, sample frequency 25us.
The parameter and population of genetic algorithm initializes: parameter to be identified has 4, therefore searching for dimension is 4, parameter value range
For the product of the long route road value of unit and route whole audience length；Population Size is usually 5-10 times of dimension, due to fault location pair
Algorithm speed proposes high requirement, therefore Population Size takes 20,22 after needing that failure is taken to occur when sampling voltage
Point.Crossover probability takes 0.6, and mutation probability 0.01, generation gap takes 0.95, and for the quick-action for guaranteeing protection, maximum genetic algebra is taken
30 times.Initial population determination process is as follows: firstly, sampling to 22 groups of voltage values after breaking down.Secondly, according to equation
(9) and (10) are it is found that the quantity that 22 groups of voltage values are corresponding with efficacious prescriptions journey is 20 groups.4 unknown numbers of demand solution, therefore can calculate
4 initial populations of C20.Finally, 20 groups of parameter value composition initial populations are therefrom randomly selected.
Embodiment 2:
Carry out parameter identification using genetic algorithm, for monopolar grounding fault situation and bipolar short circuit fault condition respectively into
Row emulation.With the increase of the number of iterations,Convergent as shown in figure 8, with mistake in figure
Crossing resistance is 10 Ω, in case where actual location distance is 5km.It knows at iteration the 20th time, for intereelectrode short-circuit failure
When with monopolar grounding fault, the value of 1/S restrains, and can guarantee convergence rate while convergence is good.
Embodiment 3:
Judge the accuracy of parameter identification, need to carry out curve fitting to formula (8) and formula (9), be write as function shape
Formula such as formula (11)-(12).Obtained result is finally recognized to be brought into formula, when with bipolar short-circuit conditions, transition resistance 10
Ω, in case where actual location distance is 5km.Fitting result is obtained as shown in figure 9, when abscissa represents using formula (11)
Between, the value of ordinate representative function y.
Emulation data are essentially coincided with matched curve, illustrate that identification result is good.
Embodiment 4:
Failure is respectively set for 18 kinds of situations, is carried out curve fitting according to (11)-(12), obtains Figure 10.Abscissa generation
18 kinds of situations of table, average per moment that ordinate represents 20 moment in each case correspond to the value of y.The formula of y indicates such as
(13).The value of j depends on fault type, when fault type is intereelectrode short-circuit failure, j=1；Conversely, j=2.
As shown in Figure 10, every kind of failure corresponding per moment be averaged y value between [- 0.2-0.2], error is very
It is small.
The specific embodiment provided according to the present invention, the invention discloses following technical effects:
The invention discloses a kind of localization method of DC distribution net line fault, the localization method includes: firstly, building
The mathematical model of vertical faulty line；Then, fitness function is constructed using the mathematical model of the faulty line；Finally, being based on
Genetic algorithm identifies the fitness function, determines the abort situation of the faulty line.The present invention, which passes through, to be built
The mathematical model of faulty line eliminates the influence of transition resistance, improves the accuracy of positioning, and constructs fitness function, will
Fault-location problem is converted into parameter identification problem, and then carries out parameter identification using genetic algorithm, determines abort situation, avoids
It is interfered caused by positioning since data acquisition is wrong, further improves the accuracy of positioning.It provides through the invention
Embodiment is it is found that positioning accuracy provided by the invention is very high, and error is 1% hereinafter, eliminating sample information in sampling process has
Interference accidentally has very strong robustness.
The invention, which can reduce production structure change bring, to be influenced, and existing product in the market is preferably tracked.Pass through
This method, index calculate in each year enter index computer capacity in agricultural product all by agricultural product macroscopic view situation upper one year carry out
Newly-increased or superseded, this method can preferably solve that each agricultural products produces every year, consumption structure changes bring otherness and asks
Topic.
Each embodiment in this specification is described in a progressive manner, the highlights of each of the examples are with other
The difference of embodiment, the same or similar parts in each embodiment may refer to each other.For system disclosed in embodiment
For, since it is corresponded to the methods disclosed in the examples, so being described relatively simple, related place is said referring to method part
It is bright.
Specific examples are used herein to describe the principles and implementation manners of the present invention, the explanation of above embodiments
Method and its core concept of the invention are merely used to help understand, described embodiment is only that a part of the invention is real
Example is applied, instead of all the embodiments, based on the embodiments of the present invention, those of ordinary skill in the art are not making creation
Property labour under the premise of every other embodiment obtained, shall fall within the protection scope of the present invention.
Claims (10)
1. a kind of localization method of DC distribution net line fault, which is characterized in that the localization method includes the following steps:
Establish the mathematical model of faulty line；
Fitness function is constructed using the mathematical model of the faulty line；
The fitness function is identified based on genetic algorithm, determines the abort situation of the faulty line.
2. a kind of localization method of DC distribution net line fault according to claim 1, which is characterized in that the foundation
The mathematical model of faulty line, specifically includes:
Obtain the first equivalent circuit of the faulty line；
According to first equivalent circuit, the first starting point mathematical model of the faulty line, first starting point are established
Mathematical model is shown below:
According to first equivalent circuit, the first end mathematical model of the faulty line, the first end mathematics are established
Model is shown below:
By first starting point mathematical model and the first end mathematical model simultaneous, the first of the faulty line is obtained
Mathematical model, first mathematical model are shown below:
Wherein, I1, R1, L1, C1, Vdc1Respectively close to the equivalent current values of the starting point of faulty line, resistance value, inductance value, electricity
The voltage value of capacitance and capacitor both ends；I2, R2, L2, C2, Vdc2Equivalent current values, electricity respectively close to the end of faulty line
Resistance value, inductance value, capacitance capacitor both ends voltage value, RfFor transition resistance value；If unit length line resistance is r, then R1=
r×D1, R2=r × D2If unit length line inductance is l, then L1=l × D1, L2=l × D2, D1For close to the route of starting point
Length, D1For by the length of the route of proximal end, D=D1+D2, D is total line length.
3. a kind of localization method of DC distribution net line fault according to claim 2, which is characterized in that the utilization
The mathematical model of the faulty line constructs fitness function, specifically includes:
By the first mathematical model discretization of the faulty line, the first discrete model is obtained, first discrete model is as follows
Shown in formula:
Wherein,K indicates the number of iterations；
First discrete model is optimized using GA algorithm, obtains the first fitness function, the first fitness letter
Number is shown below:
In formula, S (R1,R2,L1,L2) it is the first fitness function, fk(R1,R2,L1,L2) indicate kth time the first mathematical model of iteration
Functional value.
4. a kind of localization method of DC distribution net line fault according to claim 1, which is characterized in that the foundation
The mathematical model of faulty line, specifically includes:
Obtain the second equivalent circuit of the faulty line；
According to second equivalent circuit, the second starting point mathematical model of the faulty line, second starting point are established
Mathematical model is shown below:
According to second equivalent circuit, the second end mathematical model of the faulty line, the second end mathematics are established
Model is shown below:
By second starting point mathematical model and the second end mathematical model simultaneous, the second of the faulty line is obtained
Mathematical model, second mathematical model are shown below:
Wherein, I1, R1, L1, C1, Vdc1Respectively close to the equivalent current values of the starting point of faulty line, resistance value, inductance value, electricity
The voltage value of capacitance and capacitor both ends；I2, R2, L2, C2, Vdc2Equivalent current values, electricity respectively close to the end of faulty line
Resistance value, inductance value, capacitance capacitor both ends voltage value, RfFor transition resistance value；If unit length line resistance is r, then R1=
r×D1, R2=r × D2If unit length line inductance is l, then L1=l × D1, L2=l × D2, D1For close to the route of starting point
Length, D1For by the length of the route of proximal end, D=D1+D2, D is total line length.
5. a kind of localization method of DC distribution net line fault according to claim 4, which is characterized in that the utilization
The mathematical model of the faulty line constructs fitness function, specifically includes:
By the second mathematical model discretization of the faulty line, the second discrete model is obtained, second discrete model is as follows
Shown in formula:
Wherein,K indicates the number of iterations；
Second discrete model is optimized using GA algorithm, obtains the second fitness function, the second fitness letter
Number is shown below:
In formula, S'(R1,R2,L1,L2) it is the second fitness function, fk(R1,R2,L1,L2) indicate kth time the second mathematical modulo of iteration
The functional value of type.
6. a kind of localization method of DC distribution net line fault according to claim 3 or 5, which is characterized in that described
The fitness function is identified based on genetic algorithm, the abort situation of the faulty line is determined, specifically includes:
Calculate the fitness function value of each discrete location；
The discrete location that fitness function value is greater than preset first threshold value is chosen, as candidate fault locations；
Intersection and mutation operation are carried out to the candidate fault locations, the abort situation after being optimized；
The fitness function value of abort situation after calculating each optimization, the number of iterations increase by 1,；
The fitness that the maximum position of fitness function value and last iteration obtain in abort situation after judging the optimization
Whether the difference of the maximum position of functional value is less than second threshold, obtains the first judging result；
If first judging result be it is yes, using the maximum position of the fitness function value as abort situation；
If first judging result be it is no, judge whether the number of iterations is less than maximum number of iterations, obtain second and sentence
Disconnected result；
If second judging result be it is yes, return step " choose fitness function value be greater than preset first threshold value it is discrete
Position, as candidate fault locations "；
If second judging result be it is no, using the maximum position of the fitness function value as abort situation.
7. a kind of positioning system of DC distribution net line fault, which is characterized in that the positioning system includes:
Mathematical model establishes module, for establishing the mathematical model of faulty line；
Fitness function constructing module, for constructing fitness function using the mathematical model of the faulty line；
Abort situation determining module determines the fault wire for identifying based on genetic algorithm to the fitness function
The abort situation on road.
8. a kind of positioning system of DC distribution net line fault according to claim 7, which is characterized in that the mathematics
Model building module specifically includes:
First equivalent circuit acquisition submodule, for obtaining the first equivalent circuit of the faulty line；
First starting point mathematical model setting up submodule, for establishing the faulty line according to first equivalent circuit
First starting point mathematical model, first starting point mathematical model are shown below:
First end mathematical model setting up submodule, for according to first equivalent circuit, establishing the of the faulty line
One end mathematical model, the first end mathematical model are shown below:
First mathematical model setting up submodule is used for first starting point mathematical model and the first end mathematical model
Simultaneous obtains the first mathematical model of the faulty line, and first mathematical model is shown below:
Wherein, I1, R1, L1, C1, Vdc1Respectively close to the equivalent current values of the starting point of faulty line, resistance value, inductance value, electricity
The voltage value of capacitance and capacitor both ends；I2, R2, L2, C2, Vdc2Equivalent current values, electricity respectively close to the end of faulty line
Resistance value, inductance value, capacitance capacitor both ends voltage value, RfFor transition resistance value；If unit length line resistance is r, then R1=
r×D1, R2=r × D2If unit length line inductance is l, then L1=l × D1, L2=l × D2, D1For close to the route of starting point
Length, D1For by the length of the route of proximal end, D=D1+D2, D is total line length.
9. a kind of positioning system of DC distribution net line fault according to claim 8, which is characterized in that the adaptation
Function construction module is spent, is specifically included:
First mathematical model discretization submodule, for obtaining first for the first mathematical model discretization of the faulty line
Discrete model, first discrete model are shown below:
Wherein,K indicates the number of iterations；
First discrete model optimizes submodule, for optimizing using GA algorithm to first discrete model, obtains first
Fitness function, first fitness function are shown below；
In formula, S (R1,R2,L1,L2) it is the first fitness function, fk(R1,R2,L1,L2) indicate kth time the first mathematical model of iteration
Functional value.
10. a kind of positioning system of DC distribution net line fault according to claim 7, which is characterized in that the number
Model building module is learned, is specifically included:
Second equivalent circuit acquisition submodule, for obtaining the second equivalent circuit of the faulty line；
Second starting point mathematical model setting up submodule, for establishing the faulty line according to second equivalent circuit
Second starting point mathematical model, second starting point mathematical model are shown below:
Second end mathematical model setting up submodule, for according to second equivalent circuit, establishing the of the faulty line
Two end mathematical models, the second end mathematical model are shown below:
Second mathematical model setting up submodule is used for second starting point mathematical model and the second end mathematical model
Simultaneous obtains the second mathematical model of the faulty line, and second mathematical model is shown below:
Wherein, I1, R1, L1, C1, Vdc1Respectively close to the equivalent current values of the starting point of faulty line, resistance value, inductance value, electricity
The voltage value of capacitance and capacitor both ends；I2, R2, L2, C2, Vdc2Equivalent current values, electricity respectively close to the end of faulty line
Resistance value, inductance value, capacitance capacitor both ends voltage value, RfFor transition resistance value；If unit length line resistance is r, then R1=
r×D1, R2=r × D2If unit length line inductance is l, then L1=l × D1, L2=l × D2, D1For close to the route of starting point
Length, D1For by the length of the route of proximal end, D=D1+D2, D is total line length.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811353661.7A CN109342885B (en) | 2018-11-14 | 2018-11-14 | Method and system for positioning line fault of direct-current power distribution network |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811353661.7A CN109342885B (en) | 2018-11-14 | 2018-11-14 | Method and system for positioning line fault of direct-current power distribution network |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109342885A true CN109342885A (en) | 2019-02-15 |
| CN109342885B CN109342885B (en) | 2020-10-30 |
Family
ID=65315195
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811353661.7A Expired - Fee Related CN109342885B (en) | 2018-11-14 | 2018-11-14 | Method and system for positioning line fault of direct-current power distribution network |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109342885B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110556803A (en) * | 2019-10-08 | 2019-12-10 | 上海科技大学 | direct current transmission and distribution line relay protection method based on dynamic state estimation |
| CN111124884A (en) * | 2019-11-20 | 2020-05-08 | 北京航空航天大学 | Multi-fault decoupling and fault positioning method based on genetic algorithm |
| CN110323726B (en) * | 2019-07-17 | 2021-05-14 | 国网江苏省电力有限公司 | Self-adaptive line protection method and device for direct-current power distribution network |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102981099A (en) * | 2012-12-10 | 2013-03-20 | 辽宁省电力有限公司沈阳供电公司 | Location method for single-phase earth fault of power distribution network based on genetic algorithm and location device |
| CN103267926A (en) * | 2013-03-19 | 2013-08-28 | 中国石油大学(华东) | Data-gram (DG)-containing power distribution network fault distance measurement for fault feature matching based on differential evolution algorithm |
| TW201537186A (en) * | 2014-03-19 | 2015-10-01 | Univ Nat Kaohsiung Normal | Method for power line outage identification |
| CN106405319A (en) * | 2015-07-30 | 2017-02-15 | 南京理工大学 | Rough set electric power system fault diagnosis method based on heuristic information |
| CN107091972A (en) * | 2017-07-05 | 2017-08-25 | 东南大学 | A kind of active power distribution network Fault Locating Method based on improvement population |
| CN108037414A (en) * | 2017-12-11 | 2018-05-15 | 福州大学 | A kind of electrical power distribution network fault location method based on hierarchical mode and intelligent checking algorithm |
| CN108206541A (en) * | 2018-01-30 | 2018-06-26 | 国网上海市电力公司 | A kind of distribution network electric energy quality disturbance source locating method containing distributed generation resource |
-
2018
- 2018-11-14 CN CN201811353661.7A patent/CN109342885B/en not_active Expired - Fee Related
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102981099A (en) * | 2012-12-10 | 2013-03-20 | 辽宁省电力有限公司沈阳供电公司 | Location method for single-phase earth fault of power distribution network based on genetic algorithm and location device |
| CN103267926A (en) * | 2013-03-19 | 2013-08-28 | 中国石油大学(华东) | Data-gram (DG)-containing power distribution network fault distance measurement for fault feature matching based on differential evolution algorithm |
| TW201537186A (en) * | 2014-03-19 | 2015-10-01 | Univ Nat Kaohsiung Normal | Method for power line outage identification |
| CN106405319A (en) * | 2015-07-30 | 2017-02-15 | 南京理工大学 | Rough set electric power system fault diagnosis method based on heuristic information |
| CN107091972A (en) * | 2017-07-05 | 2017-08-25 | 东南大学 | A kind of active power distribution network Fault Locating Method based on improvement population |
| CN108037414A (en) * | 2017-12-11 | 2018-05-15 | 福州大学 | A kind of electrical power distribution network fault location method based on hierarchical mode and intelligent checking algorithm |
| CN108206541A (en) * | 2018-01-30 | 2018-06-26 | 国网上海市电力公司 | A kind of distribution network electric energy quality disturbance source locating method containing distributed generation resource |
Non-Patent Citations (2)
| Title |
|---|
| 殷豪 等: "纵横交叉算法在配电网故障定位中的应用", 《电力***保护与控制》 * |
| 郭壮志 等: "基于遗传算法的配电网故障定位", 《电网技术》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110323726B (en) * | 2019-07-17 | 2021-05-14 | 国网江苏省电力有限公司 | Self-adaptive line protection method and device for direct-current power distribution network |
| CN110556803A (en) * | 2019-10-08 | 2019-12-10 | 上海科技大学 | direct current transmission and distribution line relay protection method based on dynamic state estimation |
| CN110556803B (en) * | 2019-10-08 | 2021-11-26 | 上海科技大学 | Direct current transmission and distribution line relay protection method based on dynamic state estimation |
| CN111124884A (en) * | 2019-11-20 | 2020-05-08 | 北京航空航天大学 | Multi-fault decoupling and fault positioning method based on genetic algorithm |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109342885B (en) | 2020-10-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108923398B (en) | A kind of DC distribution network protection method based on voltage characteristic traveling wave Similar measure | |
| CN108469576B (en) | Direct-current fault detection method for multi-terminal alternating-current and direct-current hybrid power distribution network | |
| CN103728535A (en) | Extra-high-voltage direct-current transmission line fault location method based on wavelet transformation transient state energy spectrum | |
| CN109342885A (en) | A kind of localization method and system of DC distribution net line fault | |
| CN104597378A (en) | Fault line-selection method of power distribution network containing DG based on transient state non-power-frequency zero-sequence currents | |
| CN109256970B (en) | MMC-MTDC transmission system monopolar grounding fault current calculation method | |
| CN104537580A (en) | Method for constructing power distribution network frame based on K-means cluster | |
| CN111856322B (en) | Direct-current distribution network bipolar short-circuit fault accurate positioning method and device based on MMC | |
| CN109598377B (en) | AC/DC hybrid power distribution network robust planning method based on fault constraint | |
| CN103605829A (en) | Equivalent model method for performing electromagnetic transient simulation on alternating current/direct current hybrid power grid | |
| CN103178508B (en) | Pilot protection method of VSC-HVDC (Voltage Source Converter-High Voltage Direct Current) power transmission circuit based on shunt capacitance parameter identification | |
| CN109918762B (en) | AC/DC power grid hybrid simulation comparison model construction method | |
| CN104393590B (en) | Electrical network Transient Instability pattern INTELLIGENT IDENTIFICATION method | |
| CN108196162A (en) | The DC distribution net monopolar grounding fault wireline inspection method of model and parameters identification | |
| CN107271853A (en) | Electrical power distribution automatization system distribution low current grounding localization method and system | |
| CN104599193A (en) | Rule base based single-phase ground fault positioning method of distribution network | |
| CN110350515A (en) | A kind of flexible direct current power network modeling method suitable for analysis of the fault current | |
| CN110489806A (en) | Electromagnetic transient modeling and calculation method comprising multivoltage source type current transformer power grid | |
| CN108802570B (en) | Fault detection system and method for alternating current-direct current series-parallel micro-grid | |
| CN106779269A (en) | Reliability assessment method for secondary system of intelligent substation based on VC Method | |
| CN106033894B (en) | Method for judging stability of extra-high voltage direct current multi-drop point power grid | |
| CN104901328B (en) | Based on the soft straight control model automatic identifying method of the multiterminal under complex control network | |
| CN102496075B (en) | A kind of online data integration method based on internal memory | |
| CN108957234A (en) | Power distribution network earthing wire-selecting method and system based on volume measuring point and failure longitudinal direction feature | |
| CN113447758A (en) | Single-phase earth fault distance measuring method for multi-branch current collection circuit of wind power plant |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201030 Termination date: 20211114 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |