CN105654245A - Static voltage stability risk evaluation method of power grid based on load uncertainty modeling - Google Patents
Static voltage stability risk evaluation method of power grid based on load uncertainty modeling Download PDFInfo
- Publication number
- CN105654245A CN105654245A CN201511031154.8A CN201511031154A CN105654245A CN 105654245 A CN105654245 A CN 105654245A CN 201511031154 A CN201511031154 A CN 201511031154A CN 105654245 A CN105654245 A CN 105654245A
- Authority
- CN
- China
- Prior art keywords
- load
- node
- class
- model
- sample
- 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
- 230000003068 static effect Effects 0.000 title claims abstract description 33
- 238000011156 evaluation Methods 0.000 title abstract 10
- 238000000034 method Methods 0.000 claims description 54
- 240000002853 Nelumbo nucifera Species 0.000 claims description 30
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 30
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 30
- 230000012010 growth Effects 0.000 claims description 29
- 238000012502 risk assessment Methods 0.000 claims description 27
- 230000006698 induction Effects 0.000 claims description 15
- 230000005540 biological transmission Effects 0.000 claims description 14
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000004138 cluster model Methods 0.000 claims description 10
- 230000019771 cognition Effects 0.000 claims description 10
- 239000003102 growth factor Substances 0.000 claims description 7
- 238000000342 Monte Carlo simulation Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 230000005284 excitation Effects 0.000 claims description 3
- 238000010248 power generation Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 3
- 241000894007 species Species 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007773 growth pattern Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- 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—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Educational Administration (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention relates to a static voltage stability risk evaluation method of a power grid and particularly relates to a static voltage stability risk evaluation method of the power grid based on load uncertainty modeling, aiming at solving the problem of an existing static voltage stability risk evaluation method that an evaluation result is inaccurate. The static voltage stability risk evaluation method of the power grid based on the load uncertainty modeling is realized by the following steps: (1) establishing a load uncertainty model; (2) establishing a trend operation model; and (3) establishing a risk evaluation index, and recognizing a high-operation-risk region in the power grid by applying the risk evaluation index. The static voltage stability risk evaluation method of the power grid based on the load uncertainty modeling is applicable to static voltage stability risk evaluation of the power grid.
Description
Technical field
The present invention relates to the static electric voltage stability methods of risk assessment of electrical network, specifically a kind of uncertain based on loadThe electrical network static voltage stability methods of risk assessment of modeling.
Background technology
In the actual moving process of electrical network, quiescent voltage collapses the Voltage Instability accident causing can have a strong impact on electrical networkSteady stability operation, causes thus huge economic loss, thereby has a strong impact on national product. Therefore, in order to prevent that voltage from losingThe generation of steady accident, need to carry out risk assessment to the static electric voltage stability of electrical network. Existing static electric voltage stability riskAppraisal procedure mainly comprises following two kinds: the first is certainty appraisal procedure, and the problem of this kind of appraisal procedure existence is: due toCannot take into account the probability distribution of multiple uncertain factor, cause its assessment result to be difficult to objective, to reflect all sidedly electrical network realityBorder operation conditions, causes assessment result inaccurate thus. The second is the probability assessment method that can take into account uncertain factor, thisThe problem of planting appraisal procedure existence is: fail, taking historical load Changing Pattern as basis, to take into account accurately, all sidedly and be present in initiallyUncertainty in node serial number and load growth situation that increasing appears in load consumption, load, and lack and there is physics meaningJustice, can be from the risk indicator of reflection static voltage stability degree at all levels, therefore this kind of probability assessment method can cause equallyAssessment result is inaccurate. Given this, be necessary to invent a kind of brand-new static electric voltage stability methods of risk assessment, existing to solveThere is the inaccurate problem of static electric voltage stability methods of risk assessment assessment result.
Summary of the invention
The present invention, in order to solve the inaccurate problem of existing static electric voltage stability methods of risk assessment assessment result, providesA kind of based on load uncertain modeling electrical network static voltage stability methods of risk assessment.
The present invention adopts following technical scheme to realize: based on the electrical network static voltage stability of the uncertain modeling of loadMethods of risk assessment, the method is to adopt following steps to realize:
1) set up load uncertainty models; Described load uncertainty models comprises: first lotus fluctuation probability Distribution Model,There is the node cluster model of cognition, the load stochastic growth model that increase in load;
2) set up trend operational model: first determine in electrical network and occur according to the node cluster model of cognition of load appearance growthThe node serial number of load growth, afterwards based on monte carlo simulation methodology gather respectively in just lotus fluctuation probability Distribution Model at the beginning ofLoad increment sample in lotus consumption sample and load stochastic growth model, finally takes into account synthetic load spy based on sample valueThe continuous tide computing of property;
3) set up risk assessment index, and use the high operation risk region in risk indicator identification electrical network; Described riskEvaluation index comprises network limit load level risk indicator, the low compressive load risk indicator of node, circuit through-put power accounting windDanger index, line threshold transmission nargin risk indicator.
Described step 1) in,
The method of setting up just lotus fluctuation probability Distribution Model specifically comprises the steps:
1.1) getting clusters number is 2, generates at random one group of initial value that is subordinate to matrix element;
1.2) historical load data are carried out to fuzzy C-means clustering for the first time;
1.3) calculate the interior Classification Index BWP of class between the class of evaluating this cluster result; Specific formula for calculation is as follows:
In formula (1)-(3): BWP (j, i) is the interior Classification Index of class between class; B (j, i) is the infima species of i sample of j classSpacing; W (j, i) is the inter-object distance of i sample of j class; J and k are class mark; I, p and q are specimen number; M is cluster numbersOrder; nkIt is the number of samples in k class; njIt is the number of samples of j class;Be p sample of k class;Be j class iIndividual sample;Be q sample of j class;
1.4) again produce at random one group and be subordinate to matrix element initial value, repeating step 1.2)-1.3), until cluster number of times reachesTill historical load data length 1/2;
1.5) make clusters number add 1, generate at random one group of initial value that is subordinate to matrix element, repeating step 1.2)-1.4), straightTill clusters number reaches the evolution value of historical load data length;
1.6) add up the interior Classification Index of class between whole classes; Select the interior Classification Index maximum of class between class to be BWPoptTime correspondenceCluster result, the period dividing mode as to load power consumption time graph:
In formula (4): BWPoptFor Classification Index in class between the class of corresponding optimum cluster result; N is total sample number; Csize isThe evolution value of load sample total length;
1.7), to every type load power consumption period, calculate the coefficient correlation of all payload node power consumptions of the whole network and form phase relationMatrix number;
1.8) utilize Cholesky method decomposition step 1.7) in every correlation matrix corresponding to type load power consumption period,Obtain obeying the first lotus power consumption sample of multidimensional normal distribution; So far, the modeling of first lotus fluctuation probability Distribution Model is complete.
Load occurs that the node cluster model of cognition increasing is specifically expressed as follows:
In formula (5)-(6): ρxyFor the coefficient correlation of x and y; IxyFor the mutual information of x and y; Cov (x, y) is the association of x and yVariance; D (x) and D (y) are respectively the variance of x, y; P (x) and p (y) are respectively the marginal probability density of x, y; P (x, y) is x and yJoint probability density;
Load stochastic growth model is specifically expressed as follows:
In formula (7)-(9): the reference direction that D0 is load growth; D1 is the actual growing direction of load; (S1,S2,...,SNload) one-dimensional vector of the each payload node apparent energy of serving as reasons composition; Nload is payload node number;Being i carriesThe power factor (1≤i≤Nload) of lotus point; SΔbaseFor the power reference value of system; kLiFor load growth coefficient, kLiValueSet can be by step 1) described method determines, also can be determined by the approximating method of probability distribution.
Described step 2) in,
First, according to given threshold values ρ1、ρ2And I1、I2, determine the node serial number that occurs load growth in electrical network:
ρ1≤ρxy≤ρ2(10);
I1≤Ixy≤I2(11);
In formula (10)-(11): x and y are node serial number; ρ1And ρ2Be respectively lower limit and the upper limit of coefficient correlation; I1And I2Be respectively lower limit and the upper limit of mutual information;
Afterwards, utilize step 1) method introduced, try to achieve first lotus consumption sample in just lotus fluctuation probability Distribution Model andLoad increment sample in load stochastic growth model, determines kLiValue set;
Trend operational model is specifically expressed as follows:
In formula (12)-(13): the sustainable growth factor that λ is load; PLi0、QLi0Be respectively initial meritorious, the nothing of node iMerit power load amount;For the power factor of node i; δijFor the phase angle difference of node i and j voltage; In the time of i=j, Yii=Gii+jBiiFor node self-admittance; In the time of i ≠ j, Yij=Gij+jBijFor node transadmittance; SΔbaseFor the power reference value of system; kLiFor load growth coefficient, kGiFor the exert oneself growth factor relevant with power generation dispatching strategy; Integrated load model is commonly used constant current, perseveranceImpedance, permanent power-type are loaded and the static equivalent model of induction conductivity represents, while taking into account integrated load model, and node admittance squareThe building method of battle array is:
Yij=YIn+YZn+YMn(14);
YIn=In0/Vn0(15);
YZn=1/Zn(16);
YMn=1/(R+jX)(17);
Z2m=(rm+r2/s)2+(xm+x2)2(20);
Tm=l·(α+(1-α)(1-s)p)(21);
In formula (14)-(22): YIn、YZnAnd YMnBe respectively the equivalence of constant current type, constant-impedance type and induction motor loadAdmittance value; S, rm+jxmAnd r2/s+jx2Be respectively revolutional slip, excitation impedance and the secondary equiva lent impedance of induction conductivity; ViFor jointThe voltage magnitude of point i; L, α and p are respectively induction motor load rate, repose resistance square and mechanical load performance index;
Based on step 1) set up load uncertainty models, utilize one group of load sample of the every collection of DSMC,Just the trend operational model that synthetic load is taken into account in application calculates, until gathered whole load samples.
Described step 3) in,
Network limit load level risk indicator is specifically expressed as follows:
In formula (23)-(24): m is total simulation number of times; PLtotal(Mi) and p (Mi) be respectively the system of i kind load scenariosThe probability of ultimate load amount and appearance thereof, the implication of idle parameter and meritorious situation are similar, repeat no more;
The low compressive load risk indicator of node is specifically expressed as follows:
In formula (25)-(26):For i node voltage is the event sets that the whole network is minimum; P (*) is event AijWhen generationThe probability of system crash; Event result res (*) is the system maximum load amount after ultimate load desired value normalized;
Circuit through-put power accounting risk indicator is specifically expressed as follows:
In formula (27)-(28):For the limit transmitted power of circuit l in the whole network thing the highest with the ratio of system load amountPart set; P (*) is event LlkThe probability of system crash when generation; (p_transfer)lk(pcollapse)kBe respectively eventUpper active power and the system limits load capacity transmitting of circuit l when k occurs;
Line threshold transmission nargin risk indicator is specifically expressed as follows:
In formula (29)-(30): p (*) is event LlMThe probability of system crash when generation; (Ls-1) be that circuit transmission limit is abundantDegree; PtFor circuit sending end amount of power transfer; SbaseFor power reference value.
Compared with existing static electric voltage stability methods of risk assessment, of the present invention based on the uncertain modeling of loadElectrical network static voltage stability methods of risk assessment possess following advantage: one, compared with certainty appraisal procedure, institute of the present inventionThe electrical network static voltage stability methods of risk assessment based on the uncertain modeling of load of stating has been taken into account multiple uncertain factorProbability distribution, make thus assessment result can reflect objective, all sidedly the actual operating state of electrical network, thus make assessmentResult is more accurate. Its two, compared with probability assessment method, of the present invention based on load uncertain modeling electrical network quietState voltage stabilization methods of risk assessment, taking historical load Changing Pattern as basis, has been taken into account and has been present in original negative accurately, all sidedlyUncertainty in node serial number and load growth situation that increasing appears in lotus consumption, load, has possessed thus and has had physics meaningJustice, can be from the risk indicator of reflection static voltage stability degree at all levels, thereby make assessment result more accurate.
The present invention efficiently solves the inaccurate problem of existing static electric voltage stability methods of risk assessment assessment result, suitableFor the static electric voltage stability risk assessment of electrical network.
Brief description of the drawings
Fig. 1 is the electrical network static voltage stability methods of risk assessment flow chart based on the uncertain modeling of load.
Fig. 2 is " mutual information-relevant the closing of describing Taiyuan 110kV electrical network the whole network payload node load variations othernessSystem ".
Fig. 3 is the load stochastic growth model of Taiyuan 110kV grid nodes 15,18 and 47.
Fig. 4 is the circuit through-put power situation of change under constant power load model in certain trend computing of Taiyuan 110kV electrical network.
Fig. 5 is the circuit through-put power situation of change under integrated load model in certain trend computing of Taiyuan 110kV electrical network.
Fig. 6 be Taiyuan 110kV electrical network under three kinds of load increases, the core probability of network limit load level value-at-riskDensity Distribution.
Detailed description of the invention
Based on the electrical network static voltage stability methods of risk assessment of the uncertain modeling of load, the method is to adopt following stepRapid realization:
1) set up load uncertainty models; Described load uncertainty models comprises: first lotus fluctuation probability Distribution Model,There is the node cluster model of cognition, the load stochastic growth model that increase in load;
2) set up trend operational model: first determine in electrical network and occur according to the node cluster model of cognition of load appearance growthThe node serial number of load growth, afterwards based on monte carlo simulation methodology gather respectively in just lotus fluctuation probability Distribution Model at the beginning ofLoad increment sample in lotus consumption sample and load stochastic growth model, finally takes into account synthetic load spy based on sample valueThe continuous tide computing of property;
3) set up risk assessment index, and use the high operation risk region in risk indicator identification electrical network; Described riskEvaluation index comprises network limit load level risk indicator, the low compressive load risk indicator of node, circuit through-put power accounting windDanger index, line threshold transmission nargin risk indicator.
Described step 1) in,
The method of setting up just lotus fluctuation probability Distribution Model specifically comprises the steps:
1.1) getting clusters number is 2, generates at random one group of initial value that is subordinate to matrix element;
1.2) historical load data are carried out to fuzzy C-means clustering for the first time;
1.3) calculate the interior Classification Index BWP of class between the class of evaluating this cluster result; Specific formula for calculation is as follows:
In formula (1)-(3): BWP (j, i) is the interior Classification Index of class between class; B (j, i) is the infima species of i sample of j classSpacing; W (j, i) is the inter-object distance of i sample of j class; J and k are class mark; I, p and q are specimen number; M is cluster numbersOrder; nkIt is the number of samples in k class; njIt is the number of samples of j class;Be p sample of k class;Be j class iIndividual sample;Be q sample of j class;
1.4) again produce at random one group and be subordinate to matrix element initial value, repeating step 1.2)-1.3), until cluster number of times reachesTill historical load data length 1/2;
1.5) make clusters number add 1, generate at random one group of initial value that is subordinate to matrix element, repeating step 1.2)-1.4), straightTill clusters number reaches the evolution value of historical load data length;
1.6) add up the interior Classification Index of class between whole classes; Select the interior Classification Index maximum of class between class to be BWPoptTime correspondenceCluster result, the period dividing mode as to load power consumption time graph:
In formula (4): BWPoptFor Classification Index in class between the class of corresponding optimum cluster result; N is total sample number; Csize isThe evolution value of load sample total length;
1.7), to every type load power consumption period, calculate the coefficient correlation of all payload node power consumptions of the whole network and form phase relationMatrix number;
1.8) utilize Cholesky method decomposition step 1.7) in every correlation matrix corresponding to type load power consumption period,Obtain obeying the load power consumption sample of multidimensional normal distribution; So far, the modeling of first lotus fluctuation probability Distribution Model is complete;
Load occurs that the node cluster model of cognition increasing is specifically expressed as follows:
In formula (5)-(6): ρxyFor the coefficient correlation of x and y; IxyFor the mutual information of x and y; Cov (x, y) is the association of x and yVariance; D (x) and D (y) are respectively the variance of x, y; P (x) and p (y) are respectively the marginal probability density of x, y; P (x, y) is x and yJoint probability density;
Load stochastic growth model is specifically expressed as follows:
In formula (7)-(9): the reference direction that D0 is load growth; D1 is the actual growing direction of load; (S1,S2,...,SNload) one-dimensional vector of the each payload node apparent energy of serving as reasons composition; Nload is payload node number;Being i carriesThe power factor (1≤i≤Nload) of lotus point; SΔbaseFor the power reference value of system; kLiFor load growth coefficient, kLiValueSet can be by step 1) described method determines, also can be determined by the approximating method of probability distribution.
Described step 2) in,
First, according to given threshold values ρ1、ρ2And I1、I2, determine the node serial number that occurs load growth in electrical network:
ρ1≤ρxy≤ρ2(10);
I1≤Ixy≤I2(11);
In formula (10)-(11): x and y are node serial number; ρ1And ρ2Be respectively lower limit and the upper limit of coefficient correlation; I1And I2Be respectively lower limit and the upper limit of mutual information;
Afterwards, utilize step 1) method introduced, try to achieve first lotus consumption sample in just lotus fluctuation probability Distribution Model andLoad increment sample in load stochastic growth model, determines kLiValue set;
Trend operational model is specifically expressed as follows:
In formula (12)-(13): the sustainable growth factor that λ is load; PLi0、QLi0Be respectively initial meritorious, the nothing of node iMerit power load amount;For the power factor of node i; δijFor the phase angle difference of node i and j voltage; In the time of i=j, Yii=Gii+jBiiFor node self-admittance; In the time of i ≠ j, Yij=Gij+jBijFor node transadmittance; SΔbaseFor the power reference value of system; kLiForLoad growth coefficient, kGiFor the exert oneself growth factor relevant with power generation dispatching strategy; Integrated load model is commonly used constant current, constant-resistanceAnti-, permanent power-type is loaded and the static equivalent model of induction conductivity represents, while taking into account integrated load model, and bus admittance matrixBuilding method be:
Yij=YIn+YZn+YMn(14);
YIn=In0/Vn0(15);
YZn=1/Zn(16);
YMn=1/(R+jX)(17);
Z2m=(rm+r2/s)2+(xm+x2)2(20);
Tm=l·(α+(1-α)(1-s)p)(21);
In formula (14)-(22): YIn、YZnAnd YMnBe respectively the equivalence of constant current type, constant-impedance type and induction motor loadAdmittance value; S, rm+jxmAnd r2/s+jx2Be respectively revolutional slip, excitation impedance and the secondary equiva lent impedance of induction conductivity; ViFor jointThe voltage magnitude of point i; L, α and p are respectively induction motor load rate, repose resistance square and mechanical load performance index;
Based on step 1) set up load uncertainty models, utilize one group of load sample of the every collection of DSMC,Just the trend operational model that synthetic load is taken into account in application calculates, until gathered whole load samples.
Described step 3) in,
Network limit load level risk indicator is specifically expressed as follows:
In formula (23)-(24): m is total simulation number of times; PLtotal(Mi) and p (Mi) be respectively the system of i kind load scenariosThe probability of ultimate load amount and appearance thereof, the implication of idle parameter and meritorious situation are similar, repeat no more;
The low compressive load risk indicator of node is specifically expressed as follows:
In formula (25)-(26):For i node voltage is the event sets that the whole network is minimum; P (*) is event AijWhen generationThe probability of system crash; Event result res (*) is the system maximum load amount after ultimate load desired value normalized;
Circuit through-put power accounting risk indicator is specifically expressed as follows:
In formula (27)-(28):For the limit transmitted power of circuit l in the whole network and the ratio of system load amount the highestEvent sets; P (*) is event LlkThe probability of system crash when generation; (p_transfer)lk(pcollapse)kBe respectively thingUpper active power and the system limits load capacity transmitting of circuit l when part k occurs;
Line threshold transmission nargin risk indicator is specifically expressed as follows:
In formula (29)-(30): p (*) is event LlMThe probability of system crash when generation; (Ls-1) be that circuit transmission limit is abundantDegree; PtFor circuit sending end amount of power transfer; SbaseFor power reference value.
When concrete enforcement, taking Taiyuan 57 node 110kV electrical networks as example, technical scheme of the present invention is done further specificallyBright:
1) set up load uncertainty models:
By step 1) described flow process, at the beginning of obtaining, lotus fluctuation probability Distribution Model is: the part of first payload node is lotus justSample value is (47.53 ,-46.9 ,-58.5 ,-62.5 ,-50.5 ,-55.7 ,-43.2 ,-58.5), the portion of second payload nodeAt the beginning of point, lotus sample value is (2.2 ,-2.7 ,-3.4 ,-4.0 ,-3.5 ,-2.8 ,-3.1 ,-3.9), and other payload nodes are lotus sample justBe omitted. By step 1) method introduced, the time graph of load power consumption is divided into two time periods, in first time period, theOne and the coefficient correlation of the second payload node sample be 0.77; Second time period, the first and second payload nodes relevantCoefficient is 0.74;
By step 1) described flow process, load is occurred to the node cluster model of cognition increasing is expressed as Fig. 2: abscissa is for describingThe whole network is the mutual information of the node load difference in change opposite sex between two, and ordinate is to describe the whole network node load difference in change opposite sex between twoCoefficient correlation. Artificially the threshold values of given " mutual information-coefficient correlation " is respectively 1.5-0.9,1-0.8, and 0-0, obtains following threeCategory node group, the possibility that same class node load increases is simultaneously larger: first kind node cluster, threshold values 1.5-0.9 to 3-1 itBetween node serial number be 2,6,9,12,13,14,15,16,19,25,29,30,38,41,49,51,52,53,54,56,57; TheTwo category node groups, the node serial number of threshold values between 1-0.8 to 1.5-0.9 is 1,3,5,8,20,27,28,32,35,47,50;The 3rd category node group, the node serial number of threshold values between 0-0 to 1-0.8 is 10,17,18,23,31,33,42,43,44,55;
By step 1) described flow process, the load stochastic growth model obtaining as an example of node 15,18,47 example is expressed as Fig. 3;
2) set up trend operational model, take into account the continuous tide computing of synthetic load characteristic, record each circuit and passDefeated power and the whole network node voltage. In trend computing, make in integrated load model induction conductivity, constant-impedance, constant current and perseveranceThe ratio that power load accounts for total initial load power is respectively 0.5,0.2,0.1,0.2, and induction conductivity model is all chosen domesticTypical case's induction conductivity numerical value calculates (r1=0.04, x1=0.18, r2=0.02, x2=0.12, rm=0.35, xm=3.5, α=0.15, p=2). Statistics only increases first kind node cluster, increases by one or two category node groups and increases by three class joints simultaneously respectivelyTrend operation result in three kinds of situations of point group. Only show in certain continuous tide computing the merit of transmission on circuit 1 and 15 hereinRate size, as Fig. 4-Fig. 5;
3) set up risk assessment index, and based on risk assessment index, continuous tide operation result is carried out to risk assessment.Three kinds of horizontal risk indicator P of ultimate load that load increase is correspondingcollapseBe respectively 2312.3,2434.2,2447.7MW,QcollapseBe respectively 1050.6,1315.2,1400.1MVar, PcollapseCore probability density distribution as shown in Figure 6. Three kinds negativeThe low compressive load risk indicator of node operation result under lotus growth pattern is in table 1; Circuit transmission under three kinds of load increasesPower accounting risk indicator and line threshold transmission nargin risk indicator operation result are in table 2;
The low compressive load value-at-risk of table 1 node
;
Table 2 line power transmission risk
;
Analysis indexes operation result is known, and from causing system generation collapse of voltage angle, 31 and 33 nodes are all risky,And the value-at-risk of 31 nodes is far above 33 nodes; Circuit 1 and 15 is heavy duty elevated track through-put power accounting risk circuit, butThe line threshold transmission nargin risk of circuit 1 is very low, just the opposite with circuit 15.
Claims (4)
1. the electrical network static voltage stability methods of risk assessment based on the uncertain modeling of load, is characterized in that: the partyMethod is to adopt following steps to realize:
1) set up load uncertainty models; Described load uncertainty models comprises: first lotus fluctuation probability Distribution Model, loadThere is the node cluster model of cognition, the load stochastic growth model that increase;
2) set up trend operational model: the node cluster model of cognition that first appearance increases according to load is determined appearance load in electrical networkThe node serial number increasing, gathers respectively the first lotus consumption in just lotus fluctuation probability Distribution Model based on monte carlo simulation methodology afterwardsAmount sample and the load increment sample of load in stochastic growth model, finally take into account synthetic load characteristic based on sample valueContinuous tide computing;
3) set up risk assessment index, and use the high operation risk region in risk indicator identification electrical network; Described risk assessmentIndex comprises that network limit load level risk indicator, the low compressive load risk indicator of node, circuit through-put power accounting risk refer toMark, line threshold transmission nargin risk indicator.
2. the electrical network static voltage stability methods of risk assessment based on the uncertain modeling of load according to claim 1,It is characterized in that: described step 1) in,
The method of setting up just lotus fluctuation probability Distribution Model specifically comprises the steps:
1.1) getting clusters number is 2, generates at random one group of initial value that is subordinate to matrix element;
1.2) historical load data are carried out to fuzzy C-means clustering for the first time;
1.3) calculate the interior Classification Index BWP of class between the class of evaluating this cluster result; Specific formula for calculation is as follows:
In formula (1)-(3): BWP (j, i) is the interior Classification Index of class between class; B (j, i) is the infima species spacing of i sample of j classFrom; W (j, i) is the inter-object distance of i sample of j class; J and k are class mark; I, p and q are specimen number; M is clusters number;nkIt is the number of samples in k class; njIt is the number of samples of j class;Be p sample of k class;Be i sample of j classThis;Be q sample of j class;
1.4) again produce at random one group and be subordinate to matrix element initial value, repeating step 1.2)-1.3), go through until cluster number of times reachesHistory load data length 1/2 till;
1.5) make clusters number add 1, generate at random one group of initial value that is subordinate to matrix element, repeating step 1.2)-1.4), until poly-Till class number reaches the evolution value of historical load data length;
1.6) add up the interior Classification Index of class between whole classes; Select the interior Classification Index maximum of class between class to be BWPoptTime corresponding poly-Class result, the period dividing mode as to load power consumption time graph:
In formula (4): BWPoptFor Classification Index in class between the class of corresponding optimum cluster result; N is total sample number; Csize is loadThe evolution value of sample total length;
1.7), to every type load power consumption period, calculate the coefficient correlation of all payload node power consumptions of the whole network and form coefficient correlation squareBattle array;
1.8) utilize Cholesky method decomposition step 1.7) in every correlation matrix corresponding to type load power consumption period, obtainObey the load power consumption sample of multidimensional normal distribution; So far, the modeling of first lotus fluctuation probability Distribution Model is complete;
Load occurs that the node cluster model of cognition increasing is specifically expressed as follows:
In formula (5)-(6): ρxyFor the coefficient correlation of x and y; IxyFor the mutual information of x and y; Cov (x, y) is the covariance of x and y;D (x) and D (y) are respectively the variance of x, y; P (x) and p (y) are respectively the marginal probability density of x, y; P (x, y) is the connection of x and yClose probability density;
Load stochastic growth model is specifically expressed as follows:
In formula (7)-(9): the reference direction that D0 is load growth; D1 is the actual growing direction of load; (S1,S2,...,SNload) beThe one-dimensional vector being formed by each payload node apparent energy; Nload is payload node number;It is the merit of i the point of loadRate factor (1≤i≤Nload); SΔbaseFor the power reference value of system; kLiFor load growth coefficient, kLiThe set of value can be byStep 1) described method determines, also can be determined by the approximating method of probability distribution.
3. the electrical network static voltage stability methods of risk assessment based on the uncertain modeling of load according to claim 1,It is characterized in that: described step 2) in,
First, according to given threshold values ρ1、ρ2And I1、I2, determine the node serial number that occurs load growth in electrical network:
ρ1≤ρxy≤ρ2(10);
I1≤Ixy≤I2(11);
In formula (10)-(11): x and y are node serial number; ρ1And ρ2Be respectively lower limit and the upper limit of coefficient correlation; I1And I2RespectivelyFor lower limit and the upper limit of mutual information;
Afterwards, utilize step 1) method introduced, try to achieve first lotus consumption sample and load in just lotus fluctuation probability Distribution ModelLoad increment sample in stochastic growth model, determines kLiValue set;
Trend operational model is specifically expressed as follows:
In formula (12)-(13): the sustainable growth factor that λ is load; PLi0、QLi0Be respectively initial meritorious, the idle merit of node iRate load;For the power factor of node i; δijFor the phase angle difference of node i and j voltage; In the time of i=j, Yii=Gii+jBiiFor node self-admittance; In the time of i ≠ j, Yij=Gij+jBijFor node transadmittance; SΔbaseFor the power reference value of system; kLiFor negativeLotus growth factor, kGiFor the exert oneself growth factor relevant with power generation dispatching strategy; Integrated load model is commonly used constant current, constant-resistanceAnti-, permanent power-type is loaded and the static equivalent model of induction conductivity represents, while taking into account integrated load model, and bus admittance matrixBuilding method be:
Yij=YIn+YZn+YMn(14);
YIn=In0/Vn0(15);
YZn=1/Zn(16);
YMn=1/(R+jX)(17);
Z2m=(rm+r2/s)2+(xm+x2)2(20);
Tm=l·(α+(1-α)(1-s)p)(21);
In formula (14)-(22): YIn、YZnAnd YMnBe respectively the equivalent admittance of constant current type, constant-impedance type and induction motor loadValue; S, rm+jxmAnd r2/s+jx2Be respectively revolutional slip, excitation impedance and the secondary equiva lent impedance of induction conductivity; ViFor node iVoltage magnitude; L, α and p are respectively induction motor load rate, repose resistance square and mechanical load performance index;
Based on step 1) set up load uncertainty models, utilize one group of load sample of the every collection of DSMC, just shouldCalculate with the trend operational model of taking into account synthetic load, until gathered whole load samples.
4. the electrical network static voltage stability methods of risk assessment based on the uncertain modeling of load according to claim 1,It is characterized in that: described step 3) in,
Network limit load level risk indicator is specifically expressed as follows:
In formula (23)-(24): m is total simulation number of times; PLtotal(Mi) and p (Mi) be respectively the system limits of i kind load scenariosThe probability of load capacity and appearance thereof, the implication of idle parameter and meritorious situation are similar, repeat no more;
The low compressive load risk indicator of node is specifically expressed as follows:
In formula (25)-(26):For i node voltage is the event sets that the whole network is minimum; P (*) is event AijSystem when generationThe probability of collapse; Event result res (*) is the system maximum load amount after ultimate load desired value normalized;
Circuit through-put power accounting risk indicator is specifically expressed as follows:
In formula (27)-(28):For the limit transmitted power of circuit l in the whole network event set the highest with the ratio of system load amountClose; P (*) is event LlkThe probability of system crash when generation; (p_transfer)lk(pcollapse)kThe event k of being respectively sends outUpper active power and the system limits load capacity transmitting of circuit l when raw;
Line threshold transmission nargin risk indicator is specifically expressed as follows:
In formula (29)-(30): p (*) is event LlMThe probability of system crash when generation; (Ls-1) be circuit transmission limit nargin; PtFor circuit sending end amount of power transfer; SbaseFor power reference value.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511031154.8A CN105654245B (en) | 2015-12-31 | 2015-12-31 | Power grid static voltage stability methods of risk assessment based on negative rules modeling |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511031154.8A CN105654245B (en) | 2015-12-31 | 2015-12-31 | Power grid static voltage stability methods of risk assessment based on negative rules modeling |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105654245A true CN105654245A (en) | 2016-06-08 |
CN105654245B CN105654245B (en) | 2019-11-29 |
Family
ID=56491333
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511031154.8A Active CN105654245B (en) | 2015-12-31 | 2015-12-31 | Power grid static voltage stability methods of risk assessment based on negative rules modeling |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105654245B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106529060A (en) * | 2016-11-15 | 2017-03-22 | 中国电力科学研究院 | Load series modeling method and system |
CN107681691A (en) * | 2017-09-30 | 2018-02-09 | 太原理工大学 | The wind-electricity integration system operation reliability appraisal procedure of meter and uncertain factor |
CN107730111A (en) * | 2017-10-12 | 2018-02-23 | 国网浙江省电力公司绍兴供电公司 | A kind of distribution voltage risk evaluation model for considering customer charge and new energy access |
CN108306303A (en) * | 2018-01-17 | 2018-07-20 | 南方电网科学研究院有限责任公司 | Voltage stability evaluation method considering load increase and new energy output randomness |
CN109492851A (en) * | 2018-09-06 | 2019-03-19 | 国网浙江省电力有限公司经济技术研究院 | One kind being based on the probabilistic rack nargin appraisal procedure of different zones load growth |
CN110705879A (en) * | 2019-09-30 | 2020-01-17 | 国网山东省电力公司滨州供电公司 | Power grid vulnerability assessment method under high-proportion renewable energy access |
CN110932277A (en) * | 2019-12-26 | 2020-03-27 | 广东电网有限责任公司电力科学研究院 | Power grid static voltage stable load margin analysis method, device and equipment |
CN111429027A (en) * | 2020-04-15 | 2020-07-17 | 国网福建省电力有限公司经济技术研究院 | Regional power transmission network operation multidimensional analysis method based on big data |
CN112116235A (en) * | 2020-09-11 | 2020-12-22 | 国网山东省电力公司枣庄供电公司 | Method for evaluating influence of voltage pulse in power grid |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101295875A (en) * | 2008-06-03 | 2008-10-29 | 山东电力研究院 | Method for establishing power transmission network planning model in uncertain surroundings |
CN102509018A (en) * | 2011-11-11 | 2012-06-20 | 山东电力研究院 | System and method for evaluating importance of power system facilities |
CN103279807A (en) * | 2013-05-06 | 2013-09-04 | 国家电网公司 | Static risk assessment method for power grid in severe weather |
CN103366220A (en) * | 2012-04-06 | 2013-10-23 | 华东电力试验研究院有限公司 | Evaluation method of operational risk of electric system |
CN103870700A (en) * | 2014-03-24 | 2014-06-18 | 国家电网公司 | Distribution network static voltage stability probability assessment method based on two-point estimating method |
CN103985066A (en) * | 2014-05-20 | 2014-08-13 | 天津大学 | Static risk assessment method for power system based on hybrid power flow |
CN104659782A (en) * | 2015-03-20 | 2015-05-27 | 太原理工大学 | Power system voltage stability risk assessment method capable of considering load fluctuation limit |
-
2015
- 2015-12-31 CN CN201511031154.8A patent/CN105654245B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101295875A (en) * | 2008-06-03 | 2008-10-29 | 山东电力研究院 | Method for establishing power transmission network planning model in uncertain surroundings |
CN102509018A (en) * | 2011-11-11 | 2012-06-20 | 山东电力研究院 | System and method for evaluating importance of power system facilities |
CN103366220A (en) * | 2012-04-06 | 2013-10-23 | 华东电力试验研究院有限公司 | Evaluation method of operational risk of electric system |
CN103279807A (en) * | 2013-05-06 | 2013-09-04 | 国家电网公司 | Static risk assessment method for power grid in severe weather |
CN103870700A (en) * | 2014-03-24 | 2014-06-18 | 国家电网公司 | Distribution network static voltage stability probability assessment method based on two-point estimating method |
CN103985066A (en) * | 2014-05-20 | 2014-08-13 | 天津大学 | Static risk assessment method for power system based on hybrid power flow |
CN104659782A (en) * | 2015-03-20 | 2015-05-27 | 太原理工大学 | Power system voltage stability risk assessment method capable of considering load fluctuation limit |
Non-Patent Citations (1)
Title |
---|
MING NI等: "Online Risk-Based Security Assessment", 《IEEE POWER ENGINEERING REVIEW》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106529060A (en) * | 2016-11-15 | 2017-03-22 | 中国电力科学研究院 | Load series modeling method and system |
CN107681691A (en) * | 2017-09-30 | 2018-02-09 | 太原理工大学 | The wind-electricity integration system operation reliability appraisal procedure of meter and uncertain factor |
CN107681691B (en) * | 2017-09-30 | 2020-01-10 | 太原理工大学 | Wind power grid-connected system operation reliability assessment method considering uncertainty factors |
CN107730111A (en) * | 2017-10-12 | 2018-02-23 | 国网浙江省电力公司绍兴供电公司 | A kind of distribution voltage risk evaluation model for considering customer charge and new energy access |
CN108306303B (en) * | 2018-01-17 | 2021-06-15 | 南方电网科学研究院有限责任公司 | Voltage stability evaluation method considering load increase and new energy output randomness |
CN108306303A (en) * | 2018-01-17 | 2018-07-20 | 南方电网科学研究院有限责任公司 | Voltage stability evaluation method considering load increase and new energy output randomness |
CN109492851A (en) * | 2018-09-06 | 2019-03-19 | 国网浙江省电力有限公司经济技术研究院 | One kind being based on the probabilistic rack nargin appraisal procedure of different zones load growth |
CN109492851B (en) * | 2018-09-06 | 2021-11-30 | 国网浙江省电力有限公司经济技术研究院 | Grid frame margin evaluation method based on load growth uncertainty of different regions |
CN110705879A (en) * | 2019-09-30 | 2020-01-17 | 国网山东省电力公司滨州供电公司 | Power grid vulnerability assessment method under high-proportion renewable energy access |
CN110932277A (en) * | 2019-12-26 | 2020-03-27 | 广东电网有限责任公司电力科学研究院 | Power grid static voltage stable load margin analysis method, device and equipment |
CN111429027A (en) * | 2020-04-15 | 2020-07-17 | 国网福建省电力有限公司经济技术研究院 | Regional power transmission network operation multidimensional analysis method based on big data |
CN111429027B (en) * | 2020-04-15 | 2023-03-31 | 国网福建省电力有限公司经济技术研究院 | Regional power transmission network operation multidimensional analysis method based on big data |
CN112116235A (en) * | 2020-09-11 | 2020-12-22 | 国网山东省电力公司枣庄供电公司 | Method for evaluating influence of voltage pulse in power grid |
CN112116235B (en) * | 2020-09-11 | 2021-09-14 | 国网山东省电力公司枣庄供电公司 | Method for evaluating influence of voltage pulse in power grid |
Also Published As
Publication number | Publication date |
---|---|
CN105654245B (en) | 2019-11-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105654245A (en) | Static voltage stability risk evaluation method of power grid based on load uncertainty modeling | |
CN102945223B (en) | Method for constructing joint probability distribution function of output of a plurality of wind power plants | |
CN105866725A (en) | Method for fault classification of smart electric meter based on cluster analysis and cloud model | |
CN103279794B (en) | Electric power telecommunication network risk assessment method | |
Li et al. | Copula-ARMA model for multivariate wind speed and its applications in reliability assessment of generating systems | |
CN104410080B (en) | Method for evaluating voltage supporting capability of multi-direct-current feed-in power grid with dynamic reactive compensation | |
CN105678314A (en) | Typical demand-side user screening method based on fuzzy C clustering | |
CN109636009B (en) | Method and system for establishing neural network model for determining line loss of power grid | |
CN106529791A (en) | Evaluation method for evaluating branch importance of power system | |
CN104102832A (en) | Wind power ultrashort-term prediction method based on chaotic time series | |
CN107292502A (en) | A kind of distribution network reliability evaluation method | |
CN103530527A (en) | Wind power probability forecasting method based on numerical weather forecasting ensemble forecasting results | |
CN104020401A (en) | Cloud-model-theory-based method for evaluating insulation thermal ageing states of transformer | |
CN107358542A (en) | A kind of parameter determination method of excitation system Performance Evaluation Model | |
CN106780108A (en) | A kind of distribution transformer state evaluating method based on improvement evidential reasoning fusion | |
CN105005708A (en) | Generalized load characteristic clustering method based on AP clustering algorithm | |
CN106067034A (en) | A kind of distribution network load curve clustering method based on higher dimensional matrix characteristic root | |
CN105843733A (en) | Big data platform performance detection method and device | |
CN105225021A (en) | The optimum choice method of power distribution network project yet to be built | |
CN105406461A (en) | Adaptive dynamic load monitoring method for power distribution network power failure events | |
CN103823968A (en) | Performance evaluation method suitable for multi-region interconnected power grid contact line power control | |
CN103729570A (en) | Power system vibration mode matching method based on matrix perturbation theory | |
CN103400213A (en) | Backbone net rack survivability assessment method based on LDA (Linear Discriminant Analysis) and PCA (Principal Component Analysis) | |
CN109193791B (en) | Wind power convergence tendency state-based quantification method based on improved shape value | |
CN102904252B (en) | Method for solving uncertainty trend of power distribution network with distributed power supply |
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 |