CN104376504B - A kind of distribution system probabilistic reliability appraisal procedure based on analytic method - Google Patents
A kind of distribution system probabilistic reliability appraisal procedure based on analytic method Download PDFInfo
- Publication number
- CN104376504B CN104376504B CN201410645375.3A CN201410645375A CN104376504B CN 104376504 B CN104376504 B CN 104376504B CN 201410645375 A CN201410645375 A CN 201410645375A CN 104376504 B CN104376504 B CN 104376504B
- Authority
- CN
- China
- Prior art keywords
- power
- load
- isolated island
- time
- distribution
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 133
- 238000004458 analytical method Methods 0.000 title claims abstract description 41
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000004146 energy storage Methods 0.000 claims abstract description 58
- 238000002955 isolation Methods 0.000 claims abstract description 35
- 238000005070 sampling Methods 0.000 claims abstract description 27
- 238000000205 computational method Methods 0.000 claims abstract description 11
- 238000004088 simulation Methods 0.000 claims abstract description 11
- 238000011282 treatment Methods 0.000 claims abstract description 7
- 238000011144 upstream manufacturing Methods 0.000 claims description 25
- 239000011159 matrix material Substances 0.000 claims description 20
- 230000008439 repair process Effects 0.000 claims description 10
- 238000005286 illumination Methods 0.000 claims description 9
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 230000005611 electricity Effects 0.000 claims description 8
- 238000010248 power generation Methods 0.000 claims description 7
- 238000011160 research Methods 0.000 claims description 5
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 238000012937 correction Methods 0.000 claims description 4
- 230000001419 dependent effect Effects 0.000 claims description 4
- 230000009466 transformation Effects 0.000 claims description 4
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 3
- 230000009471 action Effects 0.000 claims description 3
- 238000004891 communication Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 abstract description 9
- 230000006870 function Effects 0.000 description 14
- 238000011156 evaluation Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 240000002853 Nelumbo nucifera Species 0.000 description 4
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 4
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 239000004744 fabric Substances 0.000 description 3
- 230000004888 barrier function Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- 108091064702 1 family Proteins 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000007257 malfunction Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000005693 optoelectronics Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
Abstract
A kind of distribution system probabilistic reliability appraisal procedure based on analytic method:Multidomain treat-ment is carried out to given power distribution network, is divided into faulty section, isolated area, seamless isolated island area and without the zone of influence, sets up Failure Mode Effective Analysis Table storehouse, initialize each parameter;Initialization simulation clock, produces random number, the crash rate parameter in each element state model obtains minimum value non-failure operation time, calculates the Fault Isolation time and load turns the band time, promote simulation clock;Failure Mode Effective Analysis Table storehouse is inquired about, classification belonging to each cell is determined, judges whether to form isolated island, adopt and be treated variously for non-isolated island area for isolated island area;The energy storage device state-of-charge probability distribution obtained according to the computational methods of probabilistic reliability, sets up the state sampling of energy storage device;The probability distribution probability of each load point single failure index is overlapped, the probabilistic reliability index of each load point and system is calculated.The present invention improves analog rate while certain calculation accuracy is ensured, comprehensively reflects grid condition.
Description
Technical field
The present invention relates to a kind of distribution system probabilistic reliability appraisal procedure.It is more particularly to a kind of to be applied to containing distribution
The distribution system probabilistic reliability appraisal procedure based on analytic method of power supply.
Background technology
In face of the opportunity to develop that energy science and technology is innovated and energy system makes the transition, modern society's energy security consciousness is progressively carried
Height, distribution Power System Reliability and then as power consumer focal point;Meanwhile, between lifting distribution Power System Reliability and investment
Marginal cost relation, directly affects the economic benefit of electric power enterprise and society, in addition, the renewable distributed electrical such as blower fan, photovoltaic
The exert oneself influence of fluctuation and energy storage device operation characteristic of source more exacerbates the complex nature of the problem.Therefore, comment overall scientific
Estimate distribution Power System Reliability to have important practical significance.
After distributed power source (DG) access power distribution network, power network becomes the network that power supply is connected with load point more than one, distribution
Deep change will all occur for the structure and the method for operation of system.Therefore, how to consider this in reliability evaluation
The influence of new system architecture and the method for operation, and how to consider that distributed power source itself is exerted oneself the influence of fluctuation, it is to work as
Research focus in preceding Model in Reliability Evaluation of Power Systems field.
By analysis it can be found that traditional analysis method considers failure generation, Fault Isolation reparation and load power
The randomness of fluctuation, influence of the randomness to distribution network reliability while consideration DG exerts oneself, using the desired value of reliability index
The reliability of sign system, can not fully reflect the fluctuation and uncertainty of active power distribution network only by expectation description;Together
When, existing method, in order to embody randomness diversity, is taken out in power supply reliability after calculating failure using to DG and load sequential
The method processing of sample, with it is same when discontinuity surface data from the sample survey carries out balancing the load, multiple sampling is tried to achieve reliability and expected, so
Computational accuracy will be influenceed by the sampling interval.In the case of at this stage, how randomness and fluctuation containing distributed power source are entered
Row analysis, sets up suitable probabilistic reliability evaluation system, to evaluating reliability of distribution network important.
In traditional reliability evaluation, according to the characteristics of power distribution network " closed loop design, open loop operation ", normally
Power network is only powered by single power supply to load point during operation.When system interior element breaks down, in the load of fault feeder section
Point is because line powering is interrupted and causes to have a power failure, and the load point after fault feeder section is needed by determining whether there is connection
Can whether network or contact spare capacity be sufficient, and then draw and restore electricity.But after DG access distribution systems, network then becomes one
The structure that individual many power supplys are connected with load point, the basic characteristic of power distribution network is changed, and this is commented to the reliability of distribution system
Estimate process and bring many new influences and problem.
Most significantly affecting for distributed power source is not only only that deep become is occurred for the method for operation for causing distribution system by it
Change, but the uncertainties of DG in itself also have a major impact to power distribution network.Also DG outage model is included in fail-safe analysis,
Consider influence of the DG failures to system.As source element, DG outage model is more complicated than two state models, so DG
Access can be significantly increased state of electric distribution network scale, and the complexity and amount of calculation of fail-safe analysis increase therewith.
In addition, in order to reduce negative effect of the distributed power source to power network, different classes of distributed power source, energy storage are filled
Put, load and corresponding control device are linked into power distribution network, be to play distributed electrical source efficiency in terms of Perspective of Energy angle
Most effective mode.So, it is actually that the reliability evaluation containing distributed power source is asked to the reliability assessment of microgrid
The extension of topic, needs also exist for paying attention to.Thus more need new processing method and analysis means.
Therefore, the reliability estimation method that probability problem is considered based on analytic method is built, is practical problem urgently to be resolved hurrily,
With good application value and construction value.
The content of the invention
The technical problems to be solved by the invention are to provide a kind of to calculate the distribution class probability decision related to user
The distribution system probabilistic reliability appraisal procedure based on analytic method of property index.
The technical solution adopted in the present invention is:A kind of distribution system probabilistic reliability appraisal procedure based on analytic method,
Comprise the following steps:
1) distribution net work structure is inputted, includes each component information, line length, each load position and each load peak, break
Road device and disconnecting switch position, power distribution network communication relationship, carry out multidomain treat-ment to each several part in power distribution network, set up fault mode shadow
Ring analysis Table storehouse;Simulation step length is set as 1h, total simulated time N, N is positive integer;Set up blower fan power output, photovoltaic generation
The time series data sequence of system output power and each load power, it is 100% charged shape to initialize each energy storage device original state
The reliability index probability description matrix W of state, initialization power off time and scarce deliverykFor 1;
2) initialization simulation clock is 0, the random number between m 0-1 is randomly generated, according in each element state model
Crash rate parameter lambda try to achieve m non-failure operation time TTF, use TTFiRepresent the non-failure operation time of i-th of element;Look for
Go out the non-failure operation time TTF of minimumi, a random number is produced to i-th of element, according to i-th of element repair rate parameter μ
Try to achieve fault correction time TTRi;At the same time, Fault Isolation time ST and Fault Isolation are produced and turns band time SRT with load,
And simulation clock is advanced to the non-failure operation time TTF of i-th of elementi;
3) Failure Mode Effective Analysis Table storehouse is inquired about, classification belonging to each cell is determined, judges whether cell forms isolated island, it is right
The cell for forming isolated island calculates the probability point of power off time and the scarce delivery of load using the computational methods of isolated island probabilistic reliability
Cloth;Cell for not forming isolated island, determines that power off time and load lack power supply using the time series data sequence of load power
The probability distribution of amount.
4) the energy storage device state-of-charge probability distribution obtained according to the computational methods of isolated island probabilistic reliability, sets up distribution
The state sampling of energy storage device, return to step 2 when net is normally run) until reaching total simulated time N of setting;
5) the probability distribution probability of the reliability index of each load point single failure in distribution net work structure is overlapped,
Calculate the probabilistic reliability index of each load point and system.
Step 1) described in in power distribution network each several part carry out multidomain treat-ment, setting up Failure Mode Effective Analysis Table storehouse is,
Power distribution network after a certain section failure of establishing power network first, failure is divided into following 6 regions:
(1) faulty section refers to all load power failures, power off time in fault element place feeder line area, the faulty section
For the repair time of fault element;
(2) upstream isolated area is positioned at the upstream in fault feeder area, the feeder line area being connected with feeder line area by non-breaker, on
Trip isolated area has a power failure after a failure, is restored electricity after Fault Isolation, and power off time is the Fault Isolation time;
(3) refer to without the zone of influence and be connected with faulty section or upstream isolated area by breaker, and with the main electricity of active distribution system
Load point in the connected feeder line area in source, no influence is not influenceed by element fault, therefore is not had a power failure;
(4) the seamless isolated island area in upstream refers to by breaker be connected with faulty section or upstream isolated area, and includes distribution
The feeder line area of power supply;Timely action of the formation of active isolated island dependent on the seamless isolated island area outlet breaker in upstream, islet operation
Time is the Fault Isolation time, and whether load point has a power failure and power off time power balance situation in island is determined, after Fault Isolation,
Upstream is seamless, and isolated island area accesses distribution system again;
(5) isolation isolated island area in downstream refers to positioned at faulty section downstream, the feeder line area being connected with faulty section by non-breaker,
Before Fault Isolation, the downstream isolates all load power failures in isolated island area, until a wide range of with other downstream areas composition
Isolated island is formed, and the islet operation time is the repair time of fault element and the difference of Fault Isolation time;
(6) the seamless isolated island area in downstream refers to be located at faulty section downstream, isolates isolated island area with faulty section or downstream and passes through open circuit
The connected feeder line area of device, downstream is seamless, and isolated island area can immediately enter isolated island mode when failure occurs runs, during islet operation
Between be the Fault Isolation time, on the basis of subregion, one by one travel through power distribution network in each component information, construct correspondence specific fault
Power distribution network other region formational situations under element, set up Failure Mode Effective Analysis Table storehouse.
Step 1) described in foundation generation blower fan power output sequence be by wind speed v approximately using Weibull distribution describe,
It is located in view of most of the time wind speed between incision wind speed and rated wind speed, then blower fan power output PwProbability density function
It is expressed as with following formula:
In formula:K and C are the shape and scale parameter of wind speed Weibull distribution respectively;Coefficient a=Prvci/(vci-vr), b=
Pr/(vr-vci), PrIt is the rated output power of blower fan, vci,vrIt is the incision wind speed and rated wind speed of blower fan respectively.
Step 1) described in photovoltaic power generation system output power sequence of setting up be to use following manner:
According to the power output P of photovoltaic generating systemsImitated with intensity of illumination I, photovoltaic battery array area S and opto-electronic conversion
Relation between rate η, i.e. photovoltaic power generation system output power Ps=IS η, if intensity of illumination I obeys beta in certain period of time
It is distributed, the probability density for obtaining photovoltaic generating system is:
In formula:Ps,maxFor the peak power output of photovoltaic array;α, β are the form parameter of intensity of illumination beta distribution.
Step 1) described in each load power sequence of setting up be to use following manner:
It is superimposed according to different time dimension typical load curve, obtains the time series data sequence L of load powert=Lp×Pw
×Pd×Ph(t), LpBy research load point annual maximum, PwFor in year corresponding with t-th hour-week load curve
Value, PdFor the value in week-daily load curve corresponding with t-th hour, Ph(t) for it is corresponding with t-th hour -when it is negative
Value in lotus curve.
Step 3) described in isolated island probabilistic reliability computational methods, comprise the following steps:
(1) the isolated island cell formed is determined, distribution net work structure and relevant parameter in isolated island cell is read;
(2) according to the power sequence of blower fan power output, photovoltaic system power output and load point, the isolated island formation phase is set up
The probability Distribution Model of interior stochastic variable;According to the charge and discharge cycles of energy storage before failure, determine that the isolated island moment forms energy storage lotus
The probability distribution of electricity condition;
(3) each rank central moment of each stochastic variable is asked for by the Monte Carlo methods of sampling;
(4) 3 estimation points are chosen in each stochastic variable, position and the weight of estimation point is calculated, the estimation of the rank of n × 3 is formed
Point group matrix and weight coefficient matrix;
(5) reliability index under each estimation point group value is calculated one by one;
(6) form the valuation matrix of each load reliability index and energy storage charge state, calculated load reliability index and
Each rank moment of the orign of energy storage charge state;
(7) each rank central moment of target stochastic variable is calculated, Gram-Charlier series expansions is used, obtains islet operation
At the end of each load point reliability index and the probability distribution of energy storage charge state.
Step 4) described in energy storage device state sampling include two parts:One is taking out for the running status of energy storage device
Sample, two be the sampling of the state-of-charge of energy storage device;Wherein, the methods of sampling of energy storage device running status is batteries string
Parallel combination into various power combinations be ranked up, put in order after normalized in [0,1] is interval, and then produce one
[0,1] interval random number, determines energy storage device running status, the sampling of state-of-charge is by producing a random number, root
According to the distribution situation of different energy storage device state-of-charges, random number is converted into the value of corresponding state-of-charge.
Step 5) described in each load point of calculating and the probabilistic reliability of system refer to calibration method and be:Obtaining each load point
After reliability index probability distribution after primary fault, if X and Y represent same load point twice in the case of different faults can
By property index probability distribution, then the probability density function for having Z=X+Y is written as
The probabilistic reliability index of each load point and system is obtained according to the method for convolution algorithm.
A kind of distribution system probabilistic reliability appraisal procedure based on analytic method of the present invention, have studied distributed power source and connects
Enter the Failure Mode Effective Analysis process of rear distribution system, by classifying to network system, and using point estimations analysis
Power off time and scarce delivery in isolated island.To blower fan, photovoltaic, batteries decile on the premise of higher level's power supply capacity is sufficient
The state of cloth power supply carries out point estimations sampling, can reflect the randomness and fluctuation of distributed power source, and ensuring
Analog rate is improved while certain calculation accuracy, comprehensively reflects grid condition.
Brief description of the drawings
Fig. 1 is two state models of non-source element;
Fig. 2 is three condition element sampling schematic diagram;
Fig. 3 is cycle discharge and recharge strategy storage battery charge state situation of change;
Fig. 4 is the IEEE-RBTS distribution network systems of transformation;
Fig. 5 is active power distribution network probabilistic reliability estimation flow figure;
Fig. 6 is active power distribution network probabilistic reliability estimation flow figure isolated island evaluation part;
Fig. 7 a are No. 13 load point APOD probability density function figures in Fig. 4;
Fig. 7 b are No. 13 load point AENS probability density function figures in Fig. 4.
Embodiment
A kind of distribution system probabilistic reliability based on analytic method of the present invention is assessed with reference to embodiment and accompanying drawing
Method is described in detail.
A kind of distribution system probabilistic reliability appraisal procedure based on analytic method of the present invention is related to user for calculating
Distribution class probabilistic reliability index.And the concept in feeder line area is combined, distribution system after distributed power source access is have studied
Failure Mode Effective Analysis process, network system is classified, and remain lonely using didactic load cutting method
Power balance in island.On the premise of higher level's power supply capacity is sufficient to blower fan, photovoltaic, batteries distributed power supply shape
State carries out point estimations sampling, can improve analog rate while certain calculation accuracy is ensured.
The fault impact classification chart of system is established, after have found the feeder line area included by every class fault zone, system
The power-off condition of internal loading point just can according to the affiliated fault zone in feeder line area where each load point it is different respectively minute
Analysis.For faulty section, without the load point in the zone of influence and the class region of upstream isolated area three, its power-off condition can be determined directly:
The power off time of load point is fault correction time in faulty section;Do not have a power failure without the load point in the zone of influence;In the isolated area of upstream
The power off time of load point is the Fault Isolation time.And it is seamless for the seamless isolated island area in upstream, downstream isolation isolated island area and downstream
Isolated island area, the power balance that its power-off condition is depended in isolated island, it is impossible to directly determine.If distributed power source is total in isolated island
Exert oneself more than total load, then the load point in island will not have a power failure;And if in isolated island distributed power source undercapacity, just need
Carry out load reduction.
Distributed power source is handled using point estimations, while utilizing cumulative distribution function and probability density function energy
The probabilistic statistical characteristicses of complete description stochastic variable.But in some practical problems, it is difficult to determine or is not required to determine random change
The definite probability distribution (PDF or CDF) of amount, only can make analysis with its some numerical characteristic.The basic thought of point estimations
It is to select discrete estimation point in known continuous random variable according to certain principle, the estimation point formation of multiple random variables is estimated
Enumeration group, the numerical characteristic of output variable to be asked is estimated using these discrete estimation point groups, such as expectation, variance and each rank square
Deng.
By function y=h (x) in μxNearby carry out Taylor expansion
Note y's is desired for μy, can be obtained by above formula:
Further can in the hope of point estimations selected point specific weight.Same provable, y high-order moment of the orign also can use
Similar approach is calculated, and formula is as follows
Further genralrlization is into the function of many variables, the calculation formula of each rank moment of the orign
On this basis, each stochastic variable is handled.Load point Calculation of Reliability can pass through following non-linear letter
Number is represented:
Y=h (x)=h (x1,x2,…,xn) (5)
Wherein y is that m ties up target stochastic variable, ykRepresent the POD or ENS of each load point, or isolated island finish time energy storage
The SOC of device;X is n dimension input stochastic variables, xkRepresent that wind speed, light intensity, the isolated island of certain node form moment energy storage charge state
Or payload.For there is k1Individual load, k2Fans, k3Group photovoltaic and k4The isolated island of group energy storage device, has
If xkExpectation, standard deviation, v ranks moment of the orign and central moment be respectively μk、σk、αk,vAnd βk,v.These inputs are random to be become
The numerical characteristic of amount can be asked for by the probabilistic model to DG and load and energy storage device model sampling.
In each stochastic variable xk(k=1,2 ..., n) on choose 3 estimation points, be denoted as xk,1、xk,2And xk,3, expression formula
For
xk,i=μk+ξk,iσk, i=1,2,3 (7)
ξ is measured in position in formulak,iExpression formula be
Wherein λk,3And λk,4Respectively stochastic variable xkThe coefficient of skewness and coefficient of kurtosis, can be by σkAnd βk,vCalculate.
It can see by formula (8), 3 estimation points that each stochastic variable is taken are respectively mean variable value and its left and right neighborhood
Interior two points for determining distance, this 3 specified points can be used as to the stochastic variable probability distribution situation such as wind speed, light intensity, load
Estimation.
N dimension input stochastic variables have n × 3 estimation point.By the estimation point of each stochastic variable respectively with its dependent variable
Expect to be combined into estimation point group, can form the rank estimation point group matrix X of n × 3, matrix element is
Xk,i=[μ1,μ2,...,μk-1,xk,i,μk+1,...,μn], i=1,2,3 (9)
Meanwhile, with each estimation point (group) to that should have certain weight coefficient, weight coefficient represents stochastic variable probability point
The probability size of the point is taken as in cloth, ξ can be measured by the position of each estimation pointk,iCalculate.Correspondence estimation point group matrix X, generates n
× 3 rank weight coefficient matrix P, matrix element is
Obtain after estimation point group, it is believed that each stochastic variable is taken as determination value by this estimation point group in during isolated island, according to
Blower fan photovoltaic is exerted oneself, payload and energy storage charge status, and strategy is cut down with reference to load, carries out balancing the load, and then calculate
Load point reliability under the determination value.The probabilistic reliability of stochastic variable will thus be calculated and be converted into multiple estimation points
The determination Calculation of Reliability of group.
Use ylRepresent the POD or ENS, h of certain load pointlBalancing the load and the Calculation of Reliability process of this load point are represented, then
Certain estimation point group Xk,iLoad point reliability valuation under value is yl=hl(Xk,i), to each element in estimation point group matrix X one by one
Calculated, form the rank valuation matrix of n × 3.Then ylV ranks moment of the orign can be expressed as
The general of reliability index or energy storage SOC can be calculated by each rank moment of the orign by deploying series using Gram-Charlier
Rate density function and cumulative distribution function.ylStandardized random variableProbability density function be
Deploy the coefficient C of series in formulaiIt is on βk,vMultinomial, can be by formula (3-38).It is rightCarry out corresponding reactionary slogan, anti-communist poster
Standardization is converted, you can obtain reliability index or energy storage SOC probability density function f (yl), complete this isolated island probability decision
Property calculate.
As shown in figure 5, a kind of distribution system probabilistic reliability appraisal procedure based on analytic method of the present invention, is specifically included
Following steps:
1) distribution net work structure is inputted, includes each component information, line length, each load position and each load peak, break
Road device and disconnecting switch position, power distribution network communication relationship, carry out multidomain treat-ment to each several part in power distribution network, set up fault mode shadow
Ring analysis (FMEA) Table storehouse;Simulation step length is set as 1h, total simulated time N, N is positive integer;Set up blower fan power output, light
The time series data sequence of photovoltaic generating system power output and each load power, it is 100% to initialize each energy storage device original state
The reliability index probability description matrix W of state-of-charge, initialization power off time and scarce deliverykFor 1;Wherein,
Described carries out multidomain treat-ment to each several part in power distribution network, and setting up Failure Mode Effective Analysis (FMEA) Table storehouse is,
Power distribution network after a certain section failure of establishing power network first, failure can be divided into following 6 regions:
(1) faulty section refers to all load power failures, power off time in fault element place feeder line area, the faulty section
For the repair time of fault element;
(2) upstream isolated area is positioned at the upstream in fault feeder area, the feeder line area being connected with feeder line area by non-breaker, on
Trip isolated area has a power failure after a failure, is restored electricity after Fault Isolation, and power off time is the Fault Isolation time;
(3) refer to without the zone of influence and be connected with faulty section or upstream isolated area by breaker, and with the main electricity of active distribution system
Load point in the connected feeder line area in source (bus), no influence is not influenceed by element fault, therefore is not had a power failure;
(4) the seamless isolated island area in upstream refers to by breaker be connected with faulty section or upstream isolated area, and includes distribution
The feeder line area of power supply;Timely action of the formation of active isolated island dependent on the seamless isolated island area outlet breaker in upstream, islet operation
Time is the Fault Isolation time, and whether load point has a power failure and power off time power balance situation in island is determined, after Fault Isolation,
Upstream is seamless, and isolated island area can access distribution system again;
(5) isolation isolated island area in downstream refers to positioned at faulty section downstream, the feeder line area being connected with faulty section by non-breaker,
Before Fault Isolation, the downstream isolates all load power failures in isolated island area, until a wide range of with other downstream areas composition
Isolated island is formed, and the islet operation time is the repair time of fault element and the difference of Fault Isolation time;
(6) the seamless isolated island area in downstream refers to be located at faulty section downstream, isolates isolated island area with faulty section or downstream and passes through open circuit
The connected feeder line area of device, downstream is seamless, and isolated island area can immediately enter isolated island mode when failure occurs runs, during islet operation
Between be the Fault Isolation time, on the basis of subregion, one by one travel through power distribution network in each component information, construct correspondence specific fault
Power distribution network other region formational situations under element, set up Failure Mode Effective Analysis Table storehouse.
Existing research thinks that wind speed has statistical property, and generally positive skewness is distributed, the function for describing wind speed profile
Or curve is more, wherein Weibull and normal distribution are widely used.Described foundation generation blower fan power output sequence is to adopt
Use following manner:Wind speed v is approximately described using Weibull distribution, it is contemplated that most of the time wind speed is located at incision wind speed and volume
Determine between wind speed, then blower fan power output PwProbability density function approximately can be expressed as with following formula:
In formula:K and C are the shape and scale parameter of wind speed Weibull distribution respectively;Coefficient a=Prvci/(vci-vr), b=
Pr/(vr-vci), PrIt is the rated output power of blower fan, vci,vrIt is the incision wind speed and rated wind speed of blower fan respectively.
Described photovoltaic power generation system output power sequence of setting up is:According to the power output P of photovoltaic generating systemsWith light
According to the relation between intensity I, photovoltaic battery array area S and photoelectric transformation efficiency η, i.e. photovoltaic power generation system output power Ps=
IS η, if intensity of illumination I is approximate in certain period of time to obey beta distribution, the probability density for obtaining photovoltaic generating system is:
In formula:Ps,maxFor the peak power output of photovoltaic array;α, β are the form parameter of intensity of illumination beta distribution.
Described each load power sequence of setting up is to use following manner:It is folded according to different time dimension typical load curve
Plus, obtain the time series data sequence L of load powert=Lp×Pw×Pd×Ph(t), LpBy research load point whole year it is maximum
Value, PwFor the value in year corresponding with t-th hour-week load curve, PdFor week-daily load curve corresponding with t-th hour
In value, Ph(t) for it is corresponding with t-th hour -when load curve in value.
2) initialization simulation clock is 0, the random number between m 0-1 is randomly generated, according in each element state model
Crash rate parameter lambda try to achieve m non-failure operation time TTF, use TTFiRepresent the non-failure operation time of i-th of element;Look for
Go out the non-failure operation time TTF of minimumi, a random number is produced to i-th of element, according to i-th of element repair rate parameter μ
Try to achieve fault correction time TTRi;At the same time, Fault Isolation time ST and Fault Isolation are produced and turns band time SRT with load,
And simulation clock is advanced to the non-failure operation time TTF of i-th of elementi;
3) Failure Mode Effective Analysis (FMEA) Table storehouse is inquired about, classification belonging to each cell is determined, judges whether cell forms
Isolated island, calculates power off time using the computational methods of isolated island probabilistic reliability to the cell for forming isolated island and load lacks delivery
Probability distribution;Cell for not forming isolated island, power off time and load are determined using the time series data sequence of load power
Lack the probability distribution of delivery.
The computational methods of described isolated island probabilistic reliability, flow as shown in Figure 6 comprises the following steps:
(1) the isolated island cell formed is determined, distribution net work structure and relevant parameter in isolated island cell is read;
(2) according to the power sequence of blower fan power output, photovoltaic system power output and load point, the isolated island formation phase is set up
The probability Distribution Model of interior stochastic variable;According to the charge and discharge cycles of energy storage before failure, determine that the isolated island moment forms energy storage lotus
The probability distribution of electricity condition;
(3) each rank central moment of each stochastic variable is asked for by the Monte Carlo methods of sampling;
(4) 3 estimation points are chosen in each stochastic variable, position and the weight of estimation point is calculated, the estimation of the rank of n × 3 is formed
Point group matrix and weight coefficient matrix;
(5) reliability index under each estimation point group value is calculated one by one;
(6) form the valuation matrix of each load reliability index and energy storage charge state, calculated load reliability index and
Each rank moment of the orign of energy storage charge state;
(7) each rank central moment of target stochastic variable is calculated, Gram-Charlier series expansions is used, obtains islet operation
At the end of each load point reliability index and the probability distribution of energy storage charge state.
4) the energy storage device state-of-charge probability distribution obtained according to the computational methods of isolated island probabilistic reliability, sets up distribution
The state sampling of energy storage device, return to step 2 when net is normally run) until reaching total simulated time N of setting;
The state sampling of described energy storage device includes two parts:One is the sampling of the running status of energy storage device, and two are
The sampling of the state-of-charge of energy storage device;Because the energy storage device in each distributed power source is also by multiple battery connection in series-parallel
Form, wherein, the methods of sampling of energy storage device running status is that the various power combinations that batteries connection in series-parallel is combined into are entered
Put in order after row sequence, normalized in [0,1] is interval, and then produce one [0,1] interval random number, determine that energy storage is filled
Running status is put, as shown in Fig. 2 drop volume running status represents the failure of part battery.The sampling of state-of-charge, is by production
A raw random number, according to the distribution situation of different energy storage device state-of-charges, corresponding state-of-charge is converted into by random number
Value.
5) the probability distribution probability of the reliability index of each load point single failure in distribution net work structure is overlapped,
Calculate the probabilistic reliability index of each load point and system.
Described each load point of calculating and the probabilistic reliability of system, which refer to calibration method, is:Obtaining the once event of each load point
After reliability index probability distribution after barrier, if X and Y represent same load point, the reliability twice in the case of different faults refers to
Mark probability distribution, then the probability density function for having Z=X+Y can be written as
The probabilistic reliability index of each load point and system is obtained according to the method for convolution algorithm.
The present invention a kind of distribution system probabilistic reliability appraisal procedure based on analytic method in, probabilistic reliability calculate with
The stochastic variables such as wind speed, illumination, load and energy storage SOC are as input condition, with the reliability index of isolated island load and isolated island knot
Energy storage SOC etc., as target variable, is stochastic variable nonlinear conversion processes more than one during beam.Probability decision based on point estimation
Property computational methods, estimation point is chosen to input stochastic variable, estimation point group and corresponding valuation weight matrix is formed;Pass through Gram-
Charlier expansion series obtains the probabilistic reliability index of isolated island internal loading.
Under the premise of herein, the heuristic load that can be formulated in isolated island cuts down strategy:
1) assume first that all loads can be supplied, calculate each total delivery of load point during isolated island;
2) if DG exerts oneself can not meet existing workload demand with current energy storage EIAJ sum, current power supply is cut down negative
The minimum load of total delivery during isolated island in lotus, will be cut in load and cuts down the time that moment to isolated island terminates being designated as Tcut, and
Into step 3), otherwise into step 4);
3) power balance calculating, repeat step (2) are carried out again, until meeting formula (5-9);
4) changed according to energy storage SOC, constantly carry out judging in step (2), until isolated island terminates.
Model is cut down according to heuristic load, during isolated island can be calculated in be cut in successively with battery discharging
Load and its reduction moment, so as to calculate the POD and ENS of each load point.
MATLAB is the business mathematics software that MathWorks companies of the U.S. produce, and is one and can be used for algorithm development, data
The advanced techniques computational language and interactive environment of visualization, data analysis and numerical computations.The present invention is using MATLAB as base
Plinth, realizes the Distribution Network Reliability assessment models accessed containing distributed power source, and the present invention is applied wherein, and to scheme
Checking is tested to application effect based on IEEE RBTS power distribution network standard examples shown in 4.
Analyzed using the multiple-limb feeder line in the IEEE RBTS Bus6 of transformation as example.Improved system is as schemed
Shown in 4, including 1 section of bus, 30 feeder line sections, 26 nodes, 23 distribution transformings, 23 load points (LP1 to LP23), 5 distributions
Formula power supply, some breakers and disconnecting switch, no fuse.Each distributed power source includes that same model blower fan is some, 1 photovoltaic
Array and 1 batteries, design parameter are as follows.
1) blower fan:Separate unit blower fan rated power 335kW;Cut wind speed 2.5m/s;Rated wind speed 12.5m/s;Cut-out wind speed
25m/s;Fitting parameter A, B, C are respectively -39.58,6.37,2.02;Mean wind speed is 19.56m/s, and wind speed profile standard deviation is
10.06m/s。
2) photovoltaic array:Parameter RcAnd GstdRespectively 0.15kW/m2And 1kW/m2。
3) battery:Every piece of battery rating 3000Ah, rated voltage 2V (6kWh);Parameter c=0.317, α=1,
K=1.22, ηc=ηd=0.927, Imax=610A.
4) element failure rate:Feeder fault rate is 0.065 times/year × km, and distribution transforming fault rate is 0.015 times/year, switch event
Barrier rate is 0.006 times/year, and mean repair time is 5h, obeys exponential distribution.Fault Isolation takes constant with the load turn band time
Value 1h.As shown in Table 1 and Table 2, the number of users of each load point is 1 family for each feeder line segment length and each load point peak value.Separate unit
Fan trouble state probability Pd=7.3%;Photovoltaic array is identical with the state model parameter of batteries, malfunction probability Pd
=3.2%, drop volume state probability Pe=5%.
The feeder line segment length of table 1
Length (km) | Feeder line section sequence number |
0.6 | 7,13 |
0.75 | 9,27 |
0.8 | 21 |
0.9 | 4,10 |
1.6 | 3,5,8,15,20,28 |
2.5 | 2,6,18,23,26 |
2.8 | 1,12,16,22,25,30 |
3.2 | 11,17,19,24,29 |
3.5 | 14 |
The load data of table 2
Load point sequence number | Load peak (kW) | Load point sequence number | Load peak (kW) |
1,6 | 360.1 | 7,23 | 796.2 |
2 | 380.6 | 8,11,14,19 | 337.6 |
3,13,17 | 653.4 | 9,21 | 737.4 |
4,18 | 686.4 | 10,12,16,22 | 340.9 |
5 | 434.7 | 15,20 | 501.8 |
Result of calculation shows that, for the load point of formation isolated island cell, its APOD and AENS show certain fluctuation
Property, some load point is for example chosen, its APOD and AENS probability density function image are as shown in Figure 7;For orphan can not be formed
The load point of island cell, its APOD is determination value, and AENS probability distribution situation is similar with Fig. 7.It can be seen that, single load point
Reliability index the approximate Normal Distribution of probability distribution.Because the electric estimation technique chooses 3 estimations to every n-dimensional random variable n
Point, only meter and the rank of each stochastic variable 5 and following central moment, therefore when using Gram-Charlier series expansions in calculating
The distribution of reliability index true probability can not be simulated completely, occur in that the situation that local probability is negative value.
Claims (4)
1. a kind of distribution system probabilistic reliability appraisal procedure based on analytic method, it is characterised in that comprise the following steps:
1) distribution net work structure is inputted, includes each component information, line length, each load position and each load peak, breaker
And disconnecting switch position, power distribution network communication relationship, by 6 regions of power distribution network, multidomain treat-ment is carried out to each region in power distribution network, built
Vertical Failure Mode Effective Analysis Table storehouse;Simulation step length is set as 1h, total simulated time N, N is positive integer;Set up blower fan output
The time series data sequence of power, photovoltaic power generation system output power and each load power, initializes each energy storage device original state
For 100% state-of-charge, the reliability index probability description matrix W of initialization power off time and scarce deliverykFor 1;
Described foundation generation blower fan power output sequence is approximately to describe wind speed v using Weibull distribution, it is contemplated that big portion
Wind speed is located between incision wind speed and rated wind speed between timesharing, then blower fan power output PwProbability density function represented with following formula
For:
In formula:K and C are the shape and scale parameter of wind speed Weibull distribution respectively;Coefficient a=Prvci/(vci-vr), b=Pr/
(vr-vci), PrIt is the rated output power of blower fan, vci,vrIt is the incision wind speed and rated wind speed of blower fan respectively;
Described photovoltaic power generation system output power sequence of setting up is to use following manner:
According to the power output P of photovoltaic generating systemsWith intensity of illumination I, photovoltaic battery array area S and photoelectric transformation efficiency η it
Between relation, i.e. photovoltaic power generation system output power Ps=IS η, if intensity of illumination I obeys beta distribution in certain period of time,
The probability density for obtaining photovoltaic generating system is:
In formula:Ps,maxFor the peak power output of photovoltaic array;α, β are the form parameter of intensity of illumination beta distribution;
2) initialization simulation clock is 0, the random number between m 0-1 is randomly generated, according to the mistake in each element state model
Efficiency parameters λ tries to achieve m non-failure operation time TTF, uses TTFiRepresent the non-failure operation time of i-th of element;Find out most
Small non-failure operation time TTFi, a random number is produced to i-th of element, tried to achieve according to i-th of element repair rate parameter μ
Fault correction time TTRi;At the same time, Fault Isolation time ST and Fault Isolation are produced and turns band time SRT with load, and will
Simulation clock is advanced to the non-failure operation time TTF of i-th of elementi;
3) Failure Mode Effective Analysis Table storehouse is inquired about, classification belonging to each cell is determined, judges whether cell forms isolated island, to being formed
The cell of isolated island calculates power off time using the computational methods of isolated island probabilistic reliability and load lacks the probability distribution of delivery;It is right
In the cell for not forming isolated island, determine that power off time and load lack the general of delivery using the time series data sequence of load power
Rate is distributed;
4) the energy storage device state-of-charge probability distribution obtained according to the computational methods of isolated island probabilistic reliability, is setting up power distribution network just
The state sampling of energy storage device, return to step 2 when often running) until reaching total simulated time N of setting;
The state sampling of described energy storage device includes two parts:One is the sampling of the running status of energy storage device, and two be energy storage
The sampling of the state-of-charge of device;Wherein, the methods of sampling of energy storage device running status is that batteries connection in series-parallel is combined into
Various power combinations be ranked up, put in order after normalized in [0,1] is interval, and then it is interval to produce one [0,1]
Random number, determines energy storage device running status, the sampling of state-of-charge is by producing a random number, according to different energy storage
The distribution situation of device state-of-charge, random number is converted into the value of corresponding state-of-charge;
5) the probability distribution probability of the reliability index of each load point single failure in distribution net work structure is overlapped, calculated
The probabilistic reliability index of each load point and system;
Described each load point of calculating and the probabilistic reliability of system, which refer to calibration method, is:After each load point primary fault is obtained
Reliability index probability distribution after, if X and Y represent same load point, the reliability index twice in the case of different faults is general
Rate is distributed, then the probability density function for having Z=X+Y is written as
The probabilistic reliability index of each load point and system is obtained according to the method for convolution algorithm.
2. a kind of distribution system probabilistic reliability appraisal procedure based on analytic method according to claim 1, its feature exists
In step 1) described in power distribution network in each several part carry out multidomain treat-ment, it is to set first to set up Failure Mode Effective Analysis Table storehouse
Power distribution network after a certain section failure of power distribution network, failure is divided into following 6 regions:
(1) all load power failures in feeder line area, the faulty section where faulty section refers to fault element, power off time is event
Hinder the repair time of element;
(2) upstream isolated area be located at fault feeder area upstream, the feeder line area being connected with feeder line area by non-breaker, upstream every
Have a power failure after a failure from area, restored electricity after Fault Isolation, power off time is the Fault Isolation time;
(3) refer to without the zone of influence and be connected with faulty section or upstream isolated area by breaker, and with active distribution system main power source phase
Load point in feeder line area even, no influence is not influenceed by element fault, therefore is not had a power failure;
(4) the seamless isolated island area in upstream refers to by breaker be connected with faulty section or upstream isolated area, and includes distributed power source
Feeder line area;Timely action of the formation of active isolated island dependent on the seamless isolated island area outlet breaker in upstream, islet operation time
For the Fault Isolation time, whether load point has a power failure and power off time power balance situation in island is determined, after Fault Isolation, upstream
Seamless isolated island area accesses distribution system again;
(5) isolation isolated island area in downstream refers to be located at faulty section downstream, the feeder line area by non-breaker being connected with faulty section, failure
Before isolation, the downstream isolates all load power failures in isolated island area, until a wide range of isolated island constituted with other downstream areas
Formed, the islet operation time is the repair time of fault element and the difference of Fault Isolation time;
(6) the seamless isolated island area in downstream refers to be located at faulty section downstream, isolates isolated island area with faulty section or downstream and passes through breaker phase
Feeder line area even, downstream is seamless, and isolated island area can immediately enter isolated island mode when failure occurs runs, and the islet operation time is
The Fault Isolation time, on the basis of subregion, each component information in power distribution network is traveled through one by one, construct correspondence specific fault element
Under power distribution network other region formational situations, set up Failure Mode Effective Analysis Table storehouse.
3. a kind of distribution system probabilistic reliability appraisal procedure based on analytic method according to claim 1, its feature exists
In step 1) described in each load power sequence of setting up be to use following manner:
It is superimposed according to different time dimension typical load curve, obtains the time series data sequence L of load powert=Lp×Pw×Pd
×Ph(t), LpBy research load point annual maximum, PwFor the value in year corresponding with t-th hour-week load curve,
PdFor the value in week-daily load curve corresponding with t-th hour, Ph(t) for it is corresponding with t-th hour -when load curve
In value.
4. a kind of distribution system probabilistic reliability appraisal procedure based on analytic method according to claim 1, its feature exists
In step 3) described in isolated island probabilistic reliability computational methods, comprise the following steps:
(1) the isolated island cell formed is determined, distribution net work structure and relevant parameter in isolated island cell is read;
(2) according to the power sequence of blower fan power output, photovoltaic system power output and load point, set up in during isolated island formation
The probability Distribution Model of stochastic variable;According to the charge and discharge cycles of energy storage before failure, determine that the isolated island moment forms the charged shape of energy storage
Probability of state is distributed;
(3) each rank central moment of each stochastic variable is asked for by the Monte Carlo methods of sampling;
(4) 3 estimation points are chosen in each stochastic variable, position and the weight of estimation point is calculated, forms the rank estimation point group of n × 3
Matrix and weight coefficient matrix;
(5) reliability index under each estimation point group value is calculated one by one;
(6) the valuation matrix of each load reliability index and energy storage charge state, calculated load reliability index and energy storage are formed
Each rank moment of the orign of state-of-charge;
(7) each rank central moment of target stochastic variable is calculated, Gram-Charlier series expansions is used, obtains islet operation and terminate
When each load point reliability index and the probability distribution of energy storage charge state.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410645375.3A CN104376504B (en) | 2014-11-06 | 2014-11-06 | A kind of distribution system probabilistic reliability appraisal procedure based on analytic method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410645375.3A CN104376504B (en) | 2014-11-06 | 2014-11-06 | A kind of distribution system probabilistic reliability appraisal procedure based on analytic method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104376504A CN104376504A (en) | 2015-02-25 |
CN104376504B true CN104376504B (en) | 2017-10-27 |
Family
ID=52555395
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410645375.3A Active CN104376504B (en) | 2014-11-06 | 2014-11-06 | A kind of distribution system probabilistic reliability appraisal procedure based on analytic method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104376504B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2723575C1 (en) * | 2019-07-01 | 2020-06-16 | Федеральное государственное бюджетное учреждение "4 Центральный научно-исследовательский институт" Министерства обороны Российской Федерации | Method of estimating probabilities of accidents in items of rocket-and-space equipment using stochastic network models of occurrence and development of emergency situations |
Families Citing this family (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105203152B (en) * | 2014-06-27 | 2018-06-19 | 国家电网公司 | A kind of photovoltaic power generation equipment failure risk exponential forecasting device and Forecasting Methodology |
CN104636993B (en) * | 2015-03-05 | 2019-12-10 | 国网山东省电力公司日照供电公司 | Power distribution system reliability algorithm |
CN106780129B (en) * | 2015-05-29 | 2021-03-30 | 江苏省电力公司常州供电公司 | Reliability evaluation method for distribution network containing distributed photovoltaic power |
CN105071381B (en) * | 2015-07-28 | 2017-04-12 | 天津大学 | State enumeration reliability evaluation method and device based on influence increment |
CN105447618B (en) * | 2015-11-06 | 2018-08-10 | 清华大学 | A kind of electric system subregion reliability estimation method |
CN105787812B (en) * | 2016-03-17 | 2022-02-25 | 清华大学 | State analysis method of MMC flexible ring network control device |
CN107294121B (en) * | 2016-04-11 | 2021-04-06 | 中国电力科学研究院 | Method and system for acquiring typical working condition curve of energy storage system |
CN106354985B (en) * | 2016-10-26 | 2020-04-10 | 华中科技大学 | Power distribution system reliability assessment method considering distributed power supply |
CN108932573B (en) * | 2017-05-25 | 2022-04-08 | 国家能源投资集团有限责任公司 | Distributed energy system stability assessment method and device |
CN107545365B (en) * | 2017-08-25 | 2020-03-17 | 合肥工业大学 | Reliability evaluation method suitable for power distribution network containing high-permeability distributed power supply |
CN107506331B (en) * | 2017-08-25 | 2023-01-10 | 国网新疆电力公司经济技术研究院 | Micro-grid reliability calculation method based on time correlation and element running time |
CN108711886B (en) * | 2018-06-08 | 2021-08-31 | 国网福建省电力有限公司 | Method for generating garden distribution network time sequence operation sample |
CN109101390B (en) * | 2018-06-29 | 2021-08-24 | 平安科技(深圳)有限公司 | Timed task abnormity monitoring method based on Gaussian distribution, electronic device and medium |
CN110188998B (en) * | 2019-05-13 | 2021-03-16 | 湖北省电力勘测设计院有限公司 | Evaluation method for time sequence construction reliability of wind turbine generator and energy storage of power distribution network |
CN110365010A (en) * | 2019-06-24 | 2019-10-22 | 国网上海市电力公司 | The index selection method of evaluating reliability of distribution network containing DG based on 2m point estimations |
CN110610303B (en) * | 2019-08-23 | 2022-05-31 | 太原理工大学 | Direct-current power distribution network reliability assessment method considering source-load uncertainty |
CN110675043B (en) * | 2019-09-17 | 2023-03-24 | 深圳供电局有限公司 | Method and system for determining power grid power failure key line based on cascading failure model |
CN110749443B (en) * | 2019-11-27 | 2021-05-18 | 济南大学 | Rolling bearing fault diagnosis method and system based on high-order origin moment |
CN111340349B (en) * | 2020-02-21 | 2023-09-05 | 中国电子产品可靠性与环境试验研究所((工业和信息化部电子第五研究所)(中国赛宝实验室)) | Product reliability evaluation method, device, computer equipment and storage medium |
CN111429008B (en) * | 2020-03-25 | 2022-05-27 | 广东电网有限责任公司 | Reliability evaluation method, device and equipment of secondary system and storage medium |
CN111475953B (en) * | 2020-04-10 | 2023-05-05 | 广东电网有限责任公司 | Energy supply reliability influence analysis method, device equipment and storage medium |
CN111860959B (en) * | 2020-06-19 | 2024-03-29 | 齐丰科技股份有限公司 | Power system cascading failure prediction method |
CN111709587B (en) * | 2020-06-22 | 2022-05-24 | 国网山西省电力公司电力科学研究院 | Power distribution system state probability evaluation method based on probability-time sequence uncertainty |
CN111861029A (en) * | 2020-07-29 | 2020-10-30 | 国网上海市电力公司 | Power supply reliability assessment method considering island division and network reconstruction |
CN111882228A (en) * | 2020-07-31 | 2020-11-03 | 国网重庆市电力公司电力科学研究院 | Reliability evaluation method for power distribution network containing distributed power supply |
CN112165120B (en) * | 2020-10-12 | 2022-04-26 | 国网山东省电力公司潍坊供电公司 | Reliability evaluation method and system for active power distribution network containing distributed power supply |
CN112199864A (en) * | 2020-11-09 | 2021-01-08 | 江苏南通发电有限公司 | Reliability optimization method for industrial user type optical storage micro-grid |
CN112668173B (en) * | 2020-12-24 | 2022-06-10 | 国网江西省电力有限公司电力科学研究院 | Method for calculating 10kV line topological relation threshold based on skewed distribution |
CN114529070B (en) * | 2022-02-08 | 2024-04-19 | 天津大学 | Comprehensive energy microgrid optimal control method considering random power failure energy supply reliability |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US228154A (en) * | 1880-05-25 | Harvester |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7519506B2 (en) * | 2002-11-06 | 2009-04-14 | Antonio Trias | System and method for monitoring and managing electrical power transmission and distribution networks |
US8200372B2 (en) * | 2008-03-31 | 2012-06-12 | The Royal Institution For The Advancement Of Learning/Mcgill University | Methods and processes for managing distributed resources in electricity power generation and distribution networks |
-
2014
- 2014-11-06 CN CN201410645375.3A patent/CN104376504B/en active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US228154A (en) * | 1880-05-25 | Harvester |
Non-Patent Citations (3)
Title |
---|
"含分布式电源的配电系统可靠性评估方法研究";王浩鸣;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20140515(第05期);参见第4.2.2.1-4.2.2.2节,图4-3 * |
"含分布式电源配电网可靠性评估的点估计法";芦晶晶 等;《电网技术》;20130831;第37卷(第8期);参见第2节 * |
"基于点估计法的有源配电网概率可靠性评估";葛少云 等;《电力系统保护与控制》;20140616;第42卷(第12期);参见第5节,图5 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2723575C1 (en) * | 2019-07-01 | 2020-06-16 | Федеральное государственное бюджетное учреждение "4 Центральный научно-исследовательский институт" Министерства обороны Российской Федерации | Method of estimating probabilities of accidents in items of rocket-and-space equipment using stochastic network models of occurrence and development of emergency situations |
Also Published As
Publication number | Publication date |
---|---|
CN104376504A (en) | 2015-02-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104376504B (en) | A kind of distribution system probabilistic reliability appraisal procedure based on analytic method | |
Hemmati et al. | Day-ahead profit-based reconfigurable microgrid scheduling considering uncertain renewable generation and load demand in the presence of energy storage | |
Orfanos et al. | Transmission expansion planning of systems with increasing wind power integration | |
Eghbali et al. | Stochastic energy management for a renewable energy based microgrid considering battery, hydrogen storage, and demand response | |
Silva et al. | Stochastic assessment of the impact of photovoltaic distributed generation on the power quality indices of distribution networks | |
Liu et al. | A Bayesian learning based scheme for online dynamic security assessment and preventive control | |
CN106972481A (en) | Scale electrically-charging equipment accesses the security quantitative estimation method of active power distribution network | |
Adewuyi et al. | Security-constrained optimal utility-scale solar PV investment planning for weak grids: Short reviews and techno-economic analysis | |
CN104156892A (en) | Active distribution network voltage drop simulation and evaluation method | |
CN103825272A (en) | Reliability determination method for power distribution network with distributed wind power based on analytical method | |
CN107730111A (en) | A kind of distribution voltage risk evaluation model for considering customer charge and new energy access | |
Li et al. | Optimal planning of Electricity–Hydrogen hybrid energy storage system considering demand response in active distribution network | |
Su et al. | An optimized algorithm for optimal power flow based on deep learning | |
Izzatillaev et al. | Short-term load forecasting in grid-connected microgrid | |
Singh et al. | Method for evaluating battery size based on loss of load probability concept for a remote PV system | |
Iweh et al. | Control and optimization of a hybrid solar PV–Hydro power system for off-grid applications using particle swarm optimization (PSO) and differential evolution (DE) | |
Osman et al. | Optimal resilient microgrids formation based on darts game theory approach and emergency demand response program for cyber-physical distribution networks considering natural disasters | |
Bakhtiari et al. | Uncertainty modeling methods for risk-averse planning and operation of stand-alone renewable energy-based microgrids | |
CN104504524A (en) | Reliability assessment method and load curtailing method applied to active distribution network | |
Zhang et al. | Reliability evaluation of high permeability renewable energy distribution network considering energy storage charge and discharge strategy | |
Aghajani et al. | Optimal energy storage sizing and offering strategy for the presence of wind power plant with energy storage in the electricity market | |
Zhang et al. | A review on capacity sizing and operation strategy of grid-connected photovoltaic battery systems | |
Gao et al. | Evaluation on the short-term power supply capacity of an active distribution system based on multiple scenarios considering uncertainties | |
Kafle et al. | Reliability Analysis Techniques in Distribution System: A Comprehensive Review | |
Bagchi et al. | Studying the impacts of incorporating energy storage devices into an aggregated probabilistic model of a virtual power plant |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |