CN108896853A - Network voltage temporary frequency reducing time appraisal procedure and device with distributed generation resource - Google Patents
Network voltage temporary frequency reducing time appraisal procedure and device with distributed generation resource Download PDFInfo
- Publication number
- CN108896853A CN108896853A CN201810901265.7A CN201810901265A CN108896853A CN 108896853 A CN108896853 A CN 108896853A CN 201810901265 A CN201810901265 A CN 201810901265A CN 108896853 A CN108896853 A CN 108896853A
- Authority
- CN
- China
- Prior art keywords
- bus
- fault
- impedance
- voltage dip
- evaluated
- 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
Abstract
This application involves a kind of temporary frequency reducing of the network voltage with distributed generation resource time appraisal procedure and devices.This method includes the node admittance matrix calculate node impedance matrix according to system;The self-impedance of fault point and mutual impedance after short circuit occurs in computing system;Calculate evaluated voltage dip amplitude of the bus under different faults type;Determine the voltage dip duration of the evaluated bus;The temporary frequency reducing time of bus is evaluated according to the voltage dip amplitude of the evaluated bus and voltage dip duration statistics.The present invention improves the accuracy to the voltage dip frequency assessment containing distributed power supply system, and institute can allow failure that position distribution rule occurs according to historical failure information and more tally with the actual situation using kernel density estimation method, keep the assessment result of voltage dip more accurate and meet reality.
Description
Technical field
This application involves automation fields, more particularly to a kind of temporary frequency reducing of the network voltage with distributed generation resource time
Appraisal procedure and device.
Background technique
With the development of infant industry and the raising of industrial automatization, more novel electric power electric equipment are put into
It uses.And sensitive equipment is extremely sensitive to voltage dip, single voltage dip event all can bring extreme loss to user.Cause
This analyze to the short trouble of electric system and accurately the assessment voltage dip generation frequency is of great significance.
There are two main classes for voltage dip frequency appraisal procedure:Survey statistic law and stochastic appraisal method.Actual measurement statistic law passes through
It measures to count the feature and the frequency that determine and temporarily drop.The time and economic cost that measurement method is spent are higher, and assessment result is difficult
To be generalized to other systems.And stochastic appraisal method assesses voltage dip by modeling and simulating, has generalization and predictability.
The regularity of distribution of different location probability of malfunction has larger impact to the assessment of the voltage dip frequency on route.Due to disliking
Transmission line short-circuit fault caused by all meetings of the factors such as bad weather, branch or animal occurs, it is difficult to predict that position occurs in short trouble
It sets, position, which occurs, for line fault has uncertainty.And the big multipair abort situation distribution of stochastic appraisal hair that existing literature proposes
Rule carries out simple hypothesis, this can cause the estimation of the voltage dip frequency deviation occur, influences the accuracy of temporary frequency reducing time assessment.
Summary of the invention
Based on this, it is necessary in view of the above technical problems, provide a kind of accuracy that can be improved temporary frequency reducing time assessment
Network voltage temporary frequency reducing time appraisal procedure and device with distributed generation resource.
First aspect present invention provides a kind of temporary frequency reducing of the network voltage with distributed generation resource time appraisal procedure, the side
Method includes:
According to the node admittance matrix calculate node impedance matrix of system;
The self-impedance of fault point and mutual impedance after being occurred according to short circuit in the nodal impedance matrix computing system;
Evaluated voltage dip width of the bus under different faults type is calculated according to the self-impedance and the mutual impedance
Value;
Determine the voltage dip duration of the evaluated bus;
Mother is evaluated according to the voltage dip amplitude of the evaluated bus and voltage dip duration statistics
The temporary frequency reducing time of line.
Optionally, the self-impedance of fault point and mutual impedance by calculating as follows after short circuit occurs in the system:
Wherein, i indicates to need to assess the bus number of the voltage dip frequency, and short trouble occurs between bus m and bus n
The position f on route, x indicate the position f and bus m normalized cumulant, (1-x) expression the position f at a distance from bus n, r 0,
1 or 2, respectively indicate zero sequence, positive sequence and negative phase-sequence; For the self-impedance of nodal impedance matrix node m, n, i;For the mutual impedance between each node;For the impedance of route m-n.
Optionally,
The fault type includes single-phase grounding fault, two-phase short-circuit fault, line to line fault ground fault and three-phase
At least one of short trouble;
When the single-phase grounding fault, the voltage dip amplitude of the evaluated bus is calculated by following formula:
When the two-phase short-circuit fault, the voltage dip amplitude of the evaluated bus is calculated by following formula:
When the line to line fault ground fault, the voltage dip amplitude of the evaluated bus is calculated by following formula:
When the three phase short circuit fault, the voltage dip amplitude of the evaluated bus is calculated by following formula:
Wherein, a is twiddle factor, value ej120°;The voltage magnitude of i point A phase before occurring for short trouble;The voltage magnitude of f point before occurring for short trouble.
Optionally, when main protection works normally, the voltage dip duration of the evaluated bus passes through following public
Formula calculates:
Wherein, [L1,L2] be the main protection of route both ends common protection scope section;
Or,
When main protection failure tripping, the voltage dip duration of the evaluated bus is calculated by following formula:
t2(x, i)=T2; (13)
Wherein, T1Fault time is removed for main protection;T2Fault time is removed for back-up protection.
Optionally, the temporary frequency reducing time of the evaluated bus is calculated by following formula:
Wherein, M and N is respectively bus sum and route sum, NB1And NL1When respectively bus and power line main protection act
The voltage dip frequency;NB2And NL2Voltage dip frequency when protection act is respectively backlogged on bus and route;K is event
Hinder type;δB,KWith δL,KIt is the probability that bus and four class short trouble of route occur respectively;L is line length;P1Main protection event
The probability of barrier;P2It is in the probability in railroad section representated by fault point for Location, G is to meet amplitude to be in
[Vmin,Vmax] interior, duration T1Or T2, position be in [L1,L2] within fault point number;F is to meet amplitude to be in
[Vmin,Vmax] interior, duration T2, position be in [L1,L2] except fault point number.
Optionally, the method also includes:
According to the history short trouble data of transmission line of electricity, short trouble when counting each short trouble by following formula
Point arrives the distance of route head end:
Wherein, LlineFor the total track length, LfaultFor short circuit trouble point to route head end distance, LnorAfter normalization
Short circuit trouble point is to route head end distance;
Determine the kernel function of corresponding kernel density estimation method;
The window width of the kernel density estimation method is calculated by following formula:
Wherein, σ is the n sample data X of stochastic variable x1, X2, X3..., XnStandard deviation;
Location of short circuit probability density is calculated by following formula:
Wherein, K () is kernel function, and h is window width;
The probability that a certain section of route is broken down is calculated by following formula:
Wherein, failure is in route section [l1,l2] in probability by integrating to obtain to probability density function.
Second aspect of the present invention also provides a kind of temporary frequency reducing of the network voltage with distributed generation resource time assessment device, packet
It includes:
First computing module, for the node admittance matrix calculate node impedance matrix according to system;
Second computing module, the resistance certainly for fault point after being occurred according to short circuit in the nodal impedance matrix computing system
Anti- and mutual impedance;
Third computing module calculates evaluated bus under different faults type according to the self-impedance and the mutual impedance
Voltage dip amplitude;
Determining module, for determining the voltage dip duration of the evaluated bus;
Statistical module, for the voltage dip amplitude and voltage dip duration statistics according to the evaluated bus
The temporary frequency reducing time of the evaluated bus.
Optionally, the self-impedance of fault point and mutual impedance by calculating as follows after short circuit occurs in the system:
Wherein, i indicates to need to assess the bus number of the voltage dip frequency, and short trouble occurs between bus m and bus n
The position f on route, x indicate the position f and bus m normalized cumulant, (1-x) expression the position f at a distance from bus n, r 0,
1 or 2, respectively indicate zero sequence, positive sequence and negative phase-sequence; For the self-impedance of nodal impedance matrix node m, n, i;For the mutual impedance between each node;For the impedance of route m-n.
Optionally,
The fault type includes single-phase grounding fault, two-phase short-circuit fault, line to line fault ground fault and three-phase
At least one of short trouble;
When the single-phase grounding fault, the voltage dip amplitude of the evaluated bus is calculated by following formula:
When the two-phase short-circuit fault, the voltage dip amplitude of the evaluated bus is calculated by following formula:
When the line to line fault ground fault, the voltage dip amplitude of the evaluated bus is calculated by following formula:
When the three phase short circuit fault, the voltage dip amplitude of the evaluated bus is calculated by following formula:
Wherein, a is twiddle factor, value ej120°;The voltage magnitude of i point A phase before occurring for short trouble;The voltage magnitude of f point before occurring for short trouble.
Optionally, when main protection works normally, the voltage dip duration of the evaluated bus passes through following public
Formula calculates:
Wherein, [L1,L2] be the main protection of route both ends common protection scope section;
Or,
When main protection failure tripping, the voltage dip duration of the evaluated bus is calculated by following formula:
t2(x, i)=T2; (13)
Wherein, T1Fault time is removed for main protection;T2Fault time is removed for back-up protection.
Optionally, the temporary frequency reducing time of the evaluated bus is calculated by following formula:
Wherein, M and N is respectively bus sum and route sum, NB1And NL1When respectively bus and power line main protection act
The voltage dip frequency;NB2And NL2Voltage dip frequency when protection act is respectively backlogged on bus and route;K is event
Hinder type;δB,KWith δL,KIt is the probability that bus and four class short trouble of route occur respectively;L is line length;P1Main protection event
The probability of barrier;P2It is in the probability in railroad section representated by fault point for Location, G is to meet amplitude to be in
[Vmin,Vmax] interior, duration T1Or T2, position be in [L1,L2] within fault point number;F is to meet amplitude to be in
[Vmin,Vmax] interior, duration T2, position be in [L1,L2] except fault point number.
Optionally, the method also includes:
According to the history short trouble data of transmission line of electricity, short trouble when counting each short trouble by following formula
Point arrives the distance of route head end:
Wherein, LlineFor the total track length, LfaultFor short circuit trouble point to route head end distance, LnorAfter normalization
Short circuit trouble point is to route head end distance;
Determine the kernel function of corresponding kernel density estimation method;
The window width of the kernel density estimation method is calculated by following formula:
Wherein, σ is the n sample data X of stochastic variable x1, X2, X3..., XnStandard deviation;
Location of short circuit probability density is calculated by following formula:
Wherein, K () is kernel function, and h is window width;
The probability that a certain section of route is broken down is calculated by following formula:
Wherein, failure is in route section [l1,l2] in probability by integrating to obtain to probability density function.
Third aspect present invention also provides a kind of computer equipment, including processor and memory, deposits on the memory
Computer program instructions are contained, realize any one of first aspect or first aspect institute when described program instruction is executed by processor
The temporary frequency reducing time appraisal procedure of the network voltage with distributed generation resource stated.
Fourth aspect present invention one kind is computer-readable to be stored in medium, is stored with computer program on the medium and refers to
It enables, realizes that first aspect or first aspect are described in any item with distributed electrical when described program instruction is executed by processor
The network voltage in source temporary frequency reducing time appraisal procedure.
In network voltage provided by the invention with distributed generation resource temporary frequency reducing time appraisal procedure, the present invention is according to history
Fault message estimates the distribution of Location using kernel density estimation method, avoids caused by subjective hypothesis partially
Difference improves the accuracy to the voltage dip frequency assessment containing distributed power supply system, and used kernel density estimation method
It can be allowed failure that position distribution rule occurs according to historical failure information more to tally with the actual situation, make the assessment result of voltage dip more
Add accurate and meets reality.
Detailed description of the invention
Fig. 1 is the schematic diagram for the network voltage temporary frequency reducing time appraisal procedure that the present invention has distributed generation resource;
Fig. 2 a is one embodiment figure of distributed generation resource equivalent circuit of the invention;
Fig. 2 b is one embodiment figure of distributed generation resource equivalent circuit of the invention;
Fig. 3 is a schematic diagram of voltage dip amplitude assessment models in the present invention;
Fig. 4 is the schematic diagram for the network voltage temporary frequency reducing time assessment device that the present invention has distributed generation resource;
Fig. 5 is the internal structure chart of computer equipment in one embodiment.
Specific embodiment
It is with reference to the accompanying drawings and embodiments, right in order to which the objects, technical solutions and advantages of the application are more clearly understood
The application is further elaborated.It should be appreciated that specific embodiment described herein is only used to explain the application, not
For limiting the application.
Current traditional stochastic appraisal method only has studied the amplitude Characteristics of voltage dip, only assumes that the duration is one fixed
Value obeys random distribution, while some voltage dip appraisal procedures do not consider main protection tripping situation, this will lead to assessment
Result information deletes mistake.Simultaneously in some voltage dip frequency assessments, line failure probability distribution subjectivity is assumed to be
Distribution situation, such as be uniformly distributed, normal distribution, exponential distribution etc..And the actual regularity of distribution may be with point of simple hypothesis
Cloth rule is different, this may result in temporarily drop assessment result and biggish error occurs.Invention new method is according to history event
The probability that barrier information is broken down using different location on kernel density estimation method estimation route, so that assessment result is more acurrate and objective
It sees.And the case where protection fault condition causes voltage dip duration extension is considered in frequency assessment, in frequency result
Contain voltage dip Duration Information.
In one embodiment, referring to Fig. 1, Fig. 1 is the network voltage temporary frequency reducing time that the present invention has distributed generation resource
The schematic diagram of appraisal procedure, as shown in Figure 1, providing a kind of temporary frequency reducing of the network voltage with distributed generation resource time assessment side
Method includes the following steps:
S101, the node admittance matrix calculate node impedance matrix according to system.
Wherein, the position accessed according to the node admittance matrix of a certain network system and distributed generation resource, to the system
Node admittance matrix is modified.Distributed generation resource can be equivalent to the series model of voltage source and impedance, as shown in Figure 2 a, Fig. 2 a
It is one embodiment figure of distributed generation resource equivalent circuit of the invention.Due to using node matrix equation to be calculated, by mould of connecting
Type promise is transformed to parallel model to keep system node admittance matrix order constant, and parallel model is as shown in Figure 2 b, and Fig. 2 b is
One embodiment figure of distributed generation resource equivalent circuit of the invention.Access distributed generation resource is equivalent to be increased in access node
One admittance only need to increase admittance value in the corresponding access node of node admittance matrix.After correcting node admittance matrix, to square
Battle array inverts and can obtain nodal impedance matrix.
S102, the self-impedance according to fault point after short circuit generation in the nodal impedance matrix computing system and mutual impedance.
The temporary frequency reducing time that need to assess bus is assessed using fault position method.Between being waited on the line according to fault position method
Every taking as fault point, self-impedance and the mutual impedance of each fault point are calculated.As shown in figure 3, Fig. 3 is that voltage is temporary in the present invention
One schematic diagram of range of decrease value assessment models, i are the buses number for needing to assess the voltage dip frequency, and short trouble occurs in bus
The position f on route between m and bus n.It is x that point f and the normalized cumulant of bus m, which occur, for short circuit, is at a distance from bus n
(1-x).Then the self-impedance of fault point f and f are to the mutual impedance of i
In formula:R=0,1,2 respectively indicate zero sequence, positive sequence and negative phase-sequence;For nodal impedance matrix node
M, the self-impedance of n, i;For the mutual impedance between each node;For the impedance of route m-n.
S103, to calculate evaluated voltage of the bus under different faults type according to the self-impedance and the mutual impedance temporary
Range of decrease value.
Wherein, consider four kinds of short trouble types, then single-phase grounding fault, two-phase short-circuit fault, line to line fault connect
When earth fault, three phase short circuit fault, the voltage dip amplitude of bus i is followed successively by
In formula, a is twiddle factor, value ej120°;The voltage magnitude of i point A phase before occurring for short trouble;The voltage magnitude of f point before occurring for short trouble.WithPerunit value can calculate when simplification be taken as 1.
S104, the voltage dip duration for determining the evaluated bus.
Wherein, the duration of voltage dip is determined by protection excision fault time, by main protection system failure possibility
It influences, when failure is removed by bus backup protection system for damping, the duration of voltage dip then can be longer.
When main protection works normally
When main protection failure tripping
t2(x, i)=T2 (8)
In formula:[L1,L2] be the main protection of route both ends common protection scope section;T1Fault time is removed for main protection;
T2Fault time is removed for back-up protection.
S105, it is commented according to the voltage dip amplitude of the evaluated bus and voltage dip duration statistics
Estimate the temporary frequency reducing time of bus.
Short trouble causes voltage dip amplitude V and continues on all short circuit trouble points or certain bus on determining certain route
After time t, the frequency that the assessment busbar voltage of statistics temporarily drops is
In formula:NB1And NL1Voltage dip frequency when respectively bus and power line main protection act;NB2And NL2It is respectively female
Voltage dip frequency when protection act is backlogged on line and route;K is fault type;δB,KWith δL,KIt is bus and route respectively
The probability that four class short troubles occur;L is line length;P1The probability of main protection failure;P2Failure is in for Location
Probability in the representative railroad section of point, is found out by (18) formula;G is to meet amplitude to be in [Vmin,Vmax] in, the duration be
T1Or T2, position be in [L1,L2] within fault point number;F is to meet amplitude to be in [Vmin,Vmax] in, the duration be
T2, position be in [L1,L2] except fault point number.
Then consider all position failures of the whole network cause assess bus the voltage dip frequency be
In formula:M and N is respectively bus sum and route sum.
Optionally, Failure probability distribution on kernel density estimation method estimation route can be used in the present invention, specifically, can be used
Such as under type:
1) statistical history fault message
According to the history short trouble data of transmission line of electricity, short circuit trouble point is to route head end when counting each short trouble
Distance.The distance of statistics is normalized, as shown in (1) formula, wherein LlineFor the total track length, LfaultIt is short
Road fault point is to route head end distance, LnorFor the short circuit trouble point after normalization to route head end distance.
Short circuit trouble point is set as stochastic variable x to route head end distance after normalization, wherein 0≤x≤1.Normalized
Short circuit trouble point afterwards is then the n sample data X of stochastic variable x to route head end distance1, X2, X3..., Xn。
2) kernel function used in kernel density estimation method is chosen
It is 1 that kernel function, which has integral, and mean value is the properties such as 0.Common kernel function has homogeneous nucleus, nucleus vestibularis triangularis, secondary core and height
This core etc..By taking Gaussian kernel as an example, function is
Influence of the selection of kernel function to estimation is little, general to choose Gaussian kernel.
3) window width of kernel density estimation method is chosen
If window width selection is excessive, estimated probability density function f (x) can be more smooth, make its sample data compared with
It is influenced by occasional case smaller in the case where few, but partial data loss in detail can be made.If selection is too small, can reflect more
More raw sample data minutias, but curve f (x) can be not smooth enough.Window width can carry out primary Calculation according to the following formula.
Wherein, σ is the n sample data X of stochastic variable x1, X2, X3..., XnStandard deviation.
Through (3) formula calculating after, window width can be adjusted according to last estimated result, keep probability density function more smooth or
Person reflects more data details features.
4) location of short circuit probability density function is calculated
The probability density function of stochastic variable x can be found out by following formula
In formula, K () is kernel function, and h is window width.
5) probability that a certain section of route is broken down is calculated
As shown in formula (5), failure is in route section [l1,l2] in probability can by probability density function product
Get.
There is the network voltage temporary frequency reducing time appraisal procedure of distributed generation resource to be illustrated the present invention above, below it is right
There is the present invention network voltage temporary frequency reducing time assessment device of distributed generation resource to be illustrated, referring to Fig. 4, Fig. 4 Fig. 4 is this hair
The schematic diagram of the bright network voltage with distributed generation resource temporary frequency reducing time assessment device, the device include:
First computing module 401, for the node admittance matrix calculate node impedance matrix according to system;
Second computing module 402, for fault point after being occurred according to short circuit in the nodal impedance matrix computing system
Self-impedance and mutual impedance;
Third computing module 403, for calculating evaluated bus in different event according to the self-impedance and the mutual impedance
Hinder the voltage dip amplitude under type;
Determining module 404, for determining the voltage dip duration of the evaluated bus;
Statistical module 405, for according to the evaluated bus voltage dip amplitude and the voltage dip duration
Count the temporary frequency reducing time of the evaluated bus.
Optionally, the self-impedance of fault point and mutual impedance by calculating as follows after short circuit occurs in the system:
Wherein, i indicates to need to assess the bus number of the voltage dip frequency, and short trouble occurs between bus m and bus n
The position f on route, x indicate the position f and bus m normalized cumulant, (1-x) expression the position f at a distance from bus n, r 0,
1 or 2, respectively indicate zero sequence, positive sequence and negative phase-sequence; For the self-impedance of nodal impedance matrix node m, n, i;For the mutual impedance between each node;For the impedance of route m-n.
Optionally,
The fault type includes single-phase grounding fault, two-phase short-circuit fault, line to line fault ground fault and three-phase
At least one of short trouble;
When the single-phase grounding fault, the voltage dip amplitude of the evaluated bus is calculated by following formula:
When the two-phase short-circuit fault, the voltage dip amplitude of the evaluated bus is calculated by following formula:
When the line to line fault ground fault, the voltage dip amplitude of the evaluated bus is calculated by following formula:
When the three phase short circuit fault, the voltage dip amplitude of the evaluated bus is calculated by following formula:
Wherein, a is twiddle factor, value ej120°;The voltage magnitude of i point A phase before occurring for short trouble;The voltage magnitude of f point before occurring for short trouble.
Optionally, when main protection works normally, the voltage dip duration of the evaluated bus passes through following public
Formula calculates:
Wherein, [L1,L2] be the main protection of route both ends common protection scope section;
Or,
When main protection failure tripping, the voltage dip duration of the evaluated bus is calculated by following formula:
t2(x, i)=T2; (13)
Wherein, T1Fault time is removed for main protection;T2Fault time is removed for back-up protection.
Optionally, the temporary frequency reducing time of the evaluated bus is calculated by following formula:
Wherein, M and N is respectively bus sum and route sum, NB1And NL1When respectively bus and power line main protection act
The voltage dip frequency;NB2And NL2Voltage dip frequency when protection act is respectively backlogged on bus and route;K is event
Hinder type;δB,KWith δL,KIt is the probability that bus and four class short trouble of route occur respectively;L is line length;P1Main protection event
The probability of barrier;P2It is in the probability in railroad section representated by fault point for Location, G is to meet amplitude to be in
[Vmin,Vmax] interior, duration T1Or T2, position be in [L1,L2] within fault point number;F is to meet amplitude to be in
[Vmin,Vmax] interior, duration T2, position be in [L1,L2] except fault point number.
Optionally, the method also includes:
According to the history short trouble data of transmission line of electricity, short trouble when counting each short trouble by following formula
Point arrives the distance of route head end:
Wherein, LlineFor the total track length, LfaultFor short circuit trouble point to route head end distance, LnorAfter normalization
Short circuit trouble point is to route head end distance;
Determine the kernel function of corresponding kernel density estimation method;
The window width of the kernel density estimation method is calculated by following formula:
Wherein, σ is the n sample data X of stochastic variable x1, X2, X3..., XnStandard deviation;
Location of short circuit probability density is calculated by following formula:
Wherein, K () is kernel function, and h is window width;
The probability that a certain section of route is broken down is calculated by following formula:
Wherein, failure is in route section [l1,l2] in probability by integrating to obtain to probability density function.
In one embodiment, a kind of computer equipment is provided, internal structure chart can be as shown in Figure 3.The calculating
Machine equipment includes processor, memory, network interface and the database connected by system bus.Wherein, the computer equipment
Processor for provide calculate and control ability.The memory of the computer equipment includes non-volatile memory medium, memory
Reservoir.The non-volatile memory medium is stored with operating system, computer program and database.The built-in storage is non-volatile
The operation of operating system and computer program in storage medium provides environment.The database of the computer equipment is believed for user
Breath, business information etc..The network interface of the computer equipment is used to communicate with external terminal by network connection.The computer
To realize the network voltage temporary frequency reducing time appraisal procedure with distributed generation resource when program is executed by processor.
It will be understood by those skilled in the art that structure shown in Fig. 3, only part relevant to application scheme is tied
The block diagram of structure does not constitute the restriction for the computer equipment being applied thereon to application scheme, specific computer equipment
It may include perhaps combining certain components or with different component layouts than more or fewer components as shown in the figure.
In one embodiment, a kind of computer equipment, including memory and processor are provided, which is stored with
Computer program, the processor realize the step in above-mentioned each embodiment of the method when executing computer program.
In one embodiment, a kind of computer readable storage medium is provided, computer program is stored thereon with, is calculated
Machine program realizes the step in above-mentioned each embodiment of the method when being executed by processor.
Those of ordinary skill in the art will appreciate that realizing all or part of the process in above-described embodiment method, being can be with
Relevant hardware is instructed to complete by computer program, the computer program can be stored in a non-volatile computer
In read/write memory medium, the computer program is when being executed, it may include such as the process of the embodiment of above-mentioned each method.Wherein,
To any reference of memory, storage, database or other media used in each embodiment provided herein,
Including non-volatile and/or volatile memory.Nonvolatile memory may include read-only memory (ROM), programming ROM
(PROM), electrically programmable ROM (EPROM), electrically erasable ROM (EEPROM) or flash memory.Volatile memory may include
Random access memory (RAM) or external cache.By way of illustration and not limitation, RAM is available in many forms,
Such as static state RAM (SRAM), dynamic ram (DRAM), synchronous dram (SDRAM), double data rate sdram (DDRSDRAM), enhancing
Type SDRAM (ESDRAM), synchronization link (Synchlink) DRAM (SLDRAM), memory bus (Rambus) direct RAM
(RDRAM), direct memory bus dynamic ram (DRDRAM) and memory bus dynamic ram (RDRAM) etc..
Each technical characteristic of above embodiments can be combined arbitrarily, for simplicity of description, not to above-described embodiment
In each technical characteristic it is all possible combination be all described, as long as however, the combination of these technical characteristics be not present lance
Shield all should be considered as described in this specification.
The several embodiments of the application above described embodiment only expresses, the description thereof is more specific and detailed, but simultaneously
It cannot therefore be construed as limiting the scope of the patent.It should be pointed out that coming for those of ordinary skill in the art
It says, without departing from the concept of this application, various modifications and improvements can be made, these belong to the protection of the application
Range.Therefore, the scope of protection shall be subject to the appended claims for the application patent.
Claims (10)
1. a kind of network voltage with distributed generation resource temporary frequency reducing time appraisal procedure, which is characterized in that the method includes:
According to the node admittance matrix calculate node impedance matrix of system;
The self-impedance of fault point and mutual impedance after being occurred according to short circuit in the nodal impedance matrix computing system;
Evaluated voltage dip amplitude of the bus under different faults type is calculated according to the self-impedance and the mutual impedance;
Determine the voltage dip duration of the evaluated bus;
Bus is evaluated according to the voltage dip amplitude of the evaluated bus and voltage dip duration statistics
Temporary frequency reducing time.
2. the network voltage according to claim 1 with distributed generation resource temporary frequency reducing time appraisal procedure, which is characterized in that
The self-impedance of fault point f and f pass through the mutual impedance of i after short circuit occurs in the system calculates as follows:
Wherein, i indicates to need to assess the bus number of the voltage dip frequency, and the route between bus m and bus n occurs for short trouble
On the position f, x indicate the position f and bus m normalized cumulant, (1-x) expression the position f at a distance from bus n, r 0,1 or
2, respectively indicate zero sequence, positive sequence and negative phase-sequence; For the self-impedance of nodal impedance matrix node m, n, i;For the mutual impedance between each node;For the impedance of route m-n.
3. the network voltage according to claim 2 with distributed generation resource temporary frequency reducing time appraisal procedure, which is characterized in that
The fault type includes single-phase grounding fault, two-phase short-circuit fault, line to line fault ground fault and three phase short circuit fault
At least one of;
When the single-phase grounding fault, the voltage dip amplitude of the evaluated bus is calculated by following formula:
When the two-phase short-circuit fault, the voltage dip amplitude of the evaluated bus is calculated by following formula:
When the line to line fault ground fault, the voltage dip amplitude of the evaluated bus is calculated by following formula:
When the three phase short circuit fault, the voltage dip amplitude of the evaluated bus is calculated by following formula:
Wherein, a is twiddle factor, value ej120;The voltage magnitude of i point A phase before occurring for short trouble;It is short
The voltage magnitude of f point before road failure occurs.
4. the network voltage according to claim 1 with distributed generation resource temporary frequency reducing time appraisal procedure, which is characterized in that
When main protection works normally, the voltage dip duration of the evaluated bus is calculated by following formula:
Wherein, [L1,L2] be the main protection of route both ends common protection scope section;
Or,
When main protection failure tripping, the voltage dip duration of the evaluated bus is calculated by following formula:
t2(x, i)=T2; (13)
Wherein, T1Fault time is removed for main protection;T2Fault time is removed for back-up protection.
5. the network voltage according to claim 1 with distributed generation resource temporary frequency reducing time appraisal procedure, which is characterized in that
The temporary frequency reducing time of the evaluated bus is calculated by following formula:
Wherein, M and N is respectively bus sum and route sum, NB1And NL1Electricity when respectively bus and power line main protection act
Press temporary frequency reducing time;NB2And NL2Voltage dip frequency when protection act is respectively backlogged on bus and route;K is failure classes
Type;δB,KWith δL,KIt is the probability that bus and four class short trouble of route occur respectively;L is line length;P1Main protection failure
Probability;P2It is in the probability in railroad section representated by fault point for Location, G is to meet amplitude to be in [Vmin,
Vmax] interior, duration T1Or T2, position be in [L1,L2] within fault point number;F is to meet amplitude to be in [Vmin,
Vmax] interior, duration T2, position be in [L1,L2] except fault point number.
6. the network voltage according to claim 5 with distributed generation resource temporary frequency reducing time appraisal procedure, which is characterized in that
The method also includes:
According to the history short trouble data of transmission line of electricity, short circuit trouble point is arrived when counting each short trouble by following formula
The distance of route head end:
Wherein, LlineFor the total track length, LfaultFor short circuit trouble point to route head end distance, LnorFor the short circuit after normalization
Fault point is to route head end distance;
Determine the kernel function of corresponding kernel density estimation method;
The window width of the kernel density estimation method is calculated by following formula:
Wherein, σ is the n sample data X of stochastic variable x1, X2, X3..., XnStandard deviation;
Location of short circuit probability density is calculated by following formula:
Wherein, K () is kernel function, and h is window width;
The probability that a certain section of route is broken down is calculated by following formula:
Wherein, failure is in route section [l1,l2] in probability by integrating to obtain to probability density function.
7. a kind of network voltage with distributed generation resource temporary frequency reducing time assessment device, which is characterized in that including:
First computing module, for the node admittance matrix calculate node impedance matrix according to system;
Second computing module, the self-impedance of fault point and mutual impedance after occurring for short circuit in computing system;
Third computing module calculates evaluated voltage dip amplitude of the bus under different faults type;
Determining module, for determining the voltage dip duration of the evaluated bus;
Statistical module, for according to the voltage dip amplitude of the evaluated bus and voltage dip duration statistics
The temporary frequency reducing time of evaluated bus.
8. the network voltage according to claim 7 with distributed generation resource temporary frequency reducing time assessment device, which is characterized in that
The self-impedance of fault point and mutual impedance by calculating as follows after short circuit occurs in the system:
Wherein, i indicates to need to assess the bus number of the voltage dip frequency, and the route between bus m and bus n occurs for short trouble
On the position f, x indicate the position f and bus m normalized cumulant, (1-x) expression the position f at a distance from bus n, r 0,1 or
2, respectively indicate zero sequence, positive sequence and negative phase-sequence; For the self-impedance of nodal impedance matrix node m, n, i;For the mutual impedance between each node;For the impedance of route m-n.
9. a kind of computer equipment, which is characterized in that including processor and memory, computer journey is stored on the memory
The tool as described in any one of claim 1 to claim 6 is realized in sequence instruction when described program instruction is executed by processor
It is distributed the network voltage temporary frequency reducing time appraisal procedure of formula power supply.
10. one kind is computer-readable to be stored in medium, which is characterized in that computer program instructions are stored on the medium, when
Described program instruction, which is realized when being executed by processor as described in any one of claim 1 to claim 6, has distribution
The network voltage of power supply temporary frequency reducing time appraisal procedure.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810901265.7A CN108896853B (en) | 2018-08-09 | 2018-08-09 | Power grid voltage sag frequency evaluation method and device with distributed power supply |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810901265.7A CN108896853B (en) | 2018-08-09 | 2018-08-09 | Power grid voltage sag frequency evaluation method and device with distributed power supply |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108896853A true CN108896853A (en) | 2018-11-27 |
CN108896853B CN108896853B (en) | 2021-06-22 |
Family
ID=64354620
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810901265.7A Active CN108896853B (en) | 2018-08-09 | 2018-08-09 | Power grid voltage sag frequency evaluation method and device with distributed power supply |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108896853B (en) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110456219A (en) * | 2019-08-12 | 2019-11-15 | 国网江西省电力有限公司电力科学研究院 | A kind of appraisal procedure of distribution line short voltage dip degree |
CN110909493A (en) * | 2019-12-30 | 2020-03-24 | 国网陕西省电力公司 | Voltage sag evaluation method based on sag domain analysis |
CN110907754A (en) * | 2019-12-02 | 2020-03-24 | 国网福建省电力有限公司莆田供电公司 | Fault line severity evaluation method based on PSD-BPA |
CN111260188A (en) * | 2020-01-08 | 2020-06-09 | 国网福建省电力有限公司 | Voltage sag treatment device configuration method |
CN111473962A (en) * | 2020-03-20 | 2020-07-31 | 广西电网有限责任公司电力科学研究院 | High-voltage circuit breaker operating spring performance evaluation method and system |
CN111983515A (en) * | 2020-08-28 | 2020-11-24 | 广东电网有限责任公司广州供电局 | Voltage sag frequency random evaluation method and device |
CN112234577A (en) * | 2020-10-21 | 2021-01-15 | 南方电网科学研究院有限责任公司 | Single-phase line break fault processing method, device and medium |
CN111830441B (en) * | 2020-07-17 | 2021-01-29 | 四川大学 | Method for estimating frequency of voltage sag caused by commissioning of transformer |
CN112345853A (en) * | 2020-10-30 | 2021-02-09 | 广东电网有限责任公司广州供电局 | Method for evaluating improvement degree of voltage sag characteristic of transformer substation bus by fault current limiter |
CN112731150A (en) * | 2020-12-08 | 2021-04-30 | 深圳供电局有限公司 | Voltage sag state estimation method and device, computer equipment and storage medium |
CN112765815A (en) * | 2021-01-21 | 2021-05-07 | 广西电网有限责任公司电力科学研究院 | Power grid lightning protection management optimization method based on artificial bee colony algorithm |
CN113281584A (en) * | 2021-04-21 | 2021-08-20 | 东南大学溧阳研究院 | Voltage sag source positioning method based on correction parameters |
CN113671290A (en) * | 2021-08-21 | 2021-11-19 | 四川大学 | Voltage sag trip rate evaluation method for non-rectangular VTC sensitive equipment |
CN114461982A (en) * | 2022-01-13 | 2022-05-10 | 四川大学 | Power transmission line protection characteristic identification and voltage sag duration estimation method |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104300532A (en) * | 2014-09-30 | 2015-01-21 | 国网山西省电力公司电力科学研究院 | Voltage sag evaluation process based on matrix factor |
CN105699849A (en) * | 2016-03-18 | 2016-06-22 | 国网上海市电力公司 | Voltage sag estimation method based on quantum-behaved particle swarm optimization algorithm |
CN105761161A (en) * | 2015-12-21 | 2016-07-13 | 国家电网公司 | AC and DC power distribution network power supply mode evaluation method |
CN105930976A (en) * | 2016-04-28 | 2016-09-07 | 华北电力大学 | Node voltage sag severity comprehensive assessment method based on weighted ideal point method |
CN106786572A (en) * | 2017-03-25 | 2017-05-31 | 西安科技大学 | The power system operating mode that configuration shnt capacitor alleviates voltage dip determines method |
-
2018
- 2018-08-09 CN CN201810901265.7A patent/CN108896853B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104300532A (en) * | 2014-09-30 | 2015-01-21 | 国网山西省电力公司电力科学研究院 | Voltage sag evaluation process based on matrix factor |
CN105761161A (en) * | 2015-12-21 | 2016-07-13 | 国家电网公司 | AC and DC power distribution network power supply mode evaluation method |
CN105699849A (en) * | 2016-03-18 | 2016-06-22 | 国网上海市电力公司 | Voltage sag estimation method based on quantum-behaved particle swarm optimization algorithm |
CN105930976A (en) * | 2016-04-28 | 2016-09-07 | 华北电力大学 | Node voltage sag severity comprehensive assessment method based on weighted ideal point method |
CN106786572A (en) * | 2017-03-25 | 2017-05-31 | 西安科技大学 | The power system operating mode that configuration shnt capacitor alleviates voltage dip determines method |
Non-Patent Citations (3)
Title |
---|
PIRJO HEINE等: "A Method for Estimating the Frequency and Cost of Voltage Sags", 《IEEE TRANSACTIONS ON POWER SYSTEMS》 * |
叶曦等: "不确定条件下计及线路保护动作特性的电压暂降频次评估", 《电力自动化设备》 * |
陈礼频等: "环网中计计距离保护时限特性的电压暂降随机评估", 《四川大学学报(工程科学版)》 * |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110456219A (en) * | 2019-08-12 | 2019-11-15 | 国网江西省电力有限公司电力科学研究院 | A kind of appraisal procedure of distribution line short voltage dip degree |
CN110907754B (en) * | 2019-12-02 | 2021-10-01 | 国网福建省电力有限公司莆田供电公司 | Fault line severity evaluation method based on PSD-BPA |
CN110907754A (en) * | 2019-12-02 | 2020-03-24 | 国网福建省电力有限公司莆田供电公司 | Fault line severity evaluation method based on PSD-BPA |
CN110909493A (en) * | 2019-12-30 | 2020-03-24 | 国网陕西省电力公司 | Voltage sag evaluation method based on sag domain analysis |
CN111260188A (en) * | 2020-01-08 | 2020-06-09 | 国网福建省电力有限公司 | Voltage sag treatment device configuration method |
CN111260188B (en) * | 2020-01-08 | 2022-06-07 | 国网福建省电力有限公司 | Voltage sag treatment device configuration method |
CN111473962A (en) * | 2020-03-20 | 2020-07-31 | 广西电网有限责任公司电力科学研究院 | High-voltage circuit breaker operating spring performance evaluation method and system |
CN111473962B (en) * | 2020-03-20 | 2022-08-16 | 广西电网有限责任公司电力科学研究院 | High-voltage circuit breaker operating spring performance evaluation method and system |
CN111830441B (en) * | 2020-07-17 | 2021-01-29 | 四川大学 | Method for estimating frequency of voltage sag caused by commissioning of transformer |
CN111983515A (en) * | 2020-08-28 | 2020-11-24 | 广东电网有限责任公司广州供电局 | Voltage sag frequency random evaluation method and device |
CN112234577A (en) * | 2020-10-21 | 2021-01-15 | 南方电网科学研究院有限责任公司 | Single-phase line break fault processing method, device and medium |
CN112234577B (en) * | 2020-10-21 | 2023-04-28 | 南方电网科学研究院有限责任公司 | Single-phase disconnection fault processing method, device and medium |
CN112345853A (en) * | 2020-10-30 | 2021-02-09 | 广东电网有限责任公司广州供电局 | Method for evaluating improvement degree of voltage sag characteristic of transformer substation bus by fault current limiter |
CN112345853B (en) * | 2020-10-30 | 2024-04-02 | 广东电网有限责任公司广州供电局 | Evaluation method for improvement degree of voltage sag characteristics of bus of transformer substation by fault current limiter |
CN112731150A (en) * | 2020-12-08 | 2021-04-30 | 深圳供电局有限公司 | Voltage sag state estimation method and device, computer equipment and storage medium |
CN112765815A (en) * | 2021-01-21 | 2021-05-07 | 广西电网有限责任公司电力科学研究院 | Power grid lightning protection management optimization method based on artificial bee colony algorithm |
CN113281584A (en) * | 2021-04-21 | 2021-08-20 | 东南大学溧阳研究院 | Voltage sag source positioning method based on correction parameters |
CN113671290A (en) * | 2021-08-21 | 2021-11-19 | 四川大学 | Voltage sag trip rate evaluation method for non-rectangular VTC sensitive equipment |
CN114461982A (en) * | 2022-01-13 | 2022-05-10 | 四川大学 | Power transmission line protection characteristic identification and voltage sag duration estimation method |
CN114461982B (en) * | 2022-01-13 | 2023-04-07 | 四川大学 | Power transmission line protection characteristic identification and voltage sag duration estimation method |
Also Published As
Publication number | Publication date |
---|---|
CN108896853B (en) | 2021-06-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108896853A (en) | Network voltage temporary frequency reducing time appraisal procedure and device with distributed generation resource | |
US20180375373A1 (en) | Impact increments-based state enumeration reliability assessment approach and device thereof | |
US7483826B2 (en) | Group-based BCU methods for on-line dynamical security assessments and energy margin calculations of practical power systems | |
CN109188138A (en) | A kind of voltage dip frequency appraisal procedure and device | |
Leite et al. | Distribution system state estimation using the Hamiltonian cycle theory | |
Timko et al. | Monte Carlo simulation of power system stability | |
Carpaneto et al. | Evaluation of the probability density functions of distribution system reliability indices with a characteristic functions-based approach | |
WO2023221587A1 (en) | Method for determining state of health of power battery of electric vehicle, and server | |
CN104463707B (en) | A kind of protection act information analysis method based on failure wave-recording | |
CN109946561A (en) | A kind of current transformer polarity testing method, apparatus, equipment and storage medium | |
Mahmood et al. | Formal reliability analysis of protective systems in smart grids | |
El-Kady | Probabilistic short-circuit analysis by Monte Carlo simulations | |
Qin et al. | A multiple model filtering approach to transmission line fault diagnosis | |
Barzegkar-Ntovom et al. | Methodology for evaluating equivalent models for the dynamic analysis of power systems | |
WO2020041302A1 (en) | A computer-aided technique for assessing infrastructure reliability and resilience and related systems, methods, and devices | |
CN109375145B (en) | Voltage loss fault monitoring method and device based on three-phase three-wire meter equipment | |
Ma et al. | Probabilistic vulnerability assessment based on power flow and voltage distribution | |
Baldwin et al. | Fault locating in distribution networks with the aid of advanced metering infrastructure | |
CN111313399B (en) | Method, device and system for confirming topological relation of regional power utilization network | |
CN106875026B (en) | Hybrid planning method for medium-and-long-term power transmission network expansion planning in power market environment | |
Pal et al. | Analysis and modeling of protection system hidden failures and its impact on power system cascading events | |
CN112633618A (en) | Reliability evaluation method and device for medium-voltage distribution network | |
JP5111258B2 (en) | FAILURE LOCATION DEVICE, FAILURE LOCATION METHOD, AND FAILURE LOCATION PROGRAM | |
EP4160867A1 (en) | Method for determining an optimal placement of measurement units for estimating the state of a physical power distribution grid | |
CN109558956A (en) | The evaluation method that goes offline, device, computer equipment and the storage medium of measuring equipment |
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 | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20200918 Address after: 510620 Tianhe District, Guangzhou, Tianhe South Road, No. two, No. 2, No. Applicant after: Guangzhou Power Supply Bureau of Guangdong Power Grid Co.,Ltd. Address before: 510620 Tianhe District, Guangzhou, Tianhe South Road, No. two, No. 2, No. Applicant before: GUANGZHOU POWER SUPPLY Co.,Ltd. |
|
TA01 | Transfer of patent application right | ||
GR01 | Patent grant | ||
GR01 | Patent grant |