CN105680442B - Consider that the expectation of trend and sensitivity consistency equivalence lacks power supply volume appraisal procedure - Google Patents
Consider that the expectation of trend and sensitivity consistency equivalence lacks power supply volume appraisal procedure Download PDFInfo
- Publication number
- CN105680442B CN105680442B CN201610127267.6A CN201610127267A CN105680442B CN 105680442 B CN105680442 B CN 105680442B CN 201610127267 A CN201610127267 A CN 201610127267A CN 105680442 B CN105680442 B CN 105680442B
- Authority
- CN
- China
- Prior art keywords
- intranet
- node
- network
- trend
- equivalence
- 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
- 238000000034 method Methods 0.000 title claims abstract description 47
- 230000035945 sensitivity Effects 0.000 title claims abstract description 35
- 238000004458 analytical method Methods 0.000 claims abstract description 26
- 230000003068 static effect Effects 0.000 claims abstract description 11
- 238000000342 Monte Carlo simulation Methods 0.000 claims abstract description 9
- 238000010937 topological data analysis Methods 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 7
- 230000005540 biological transmission Effects 0.000 claims description 6
- 230000009977 dual effect Effects 0.000 claims description 5
- 238000005457 optimization Methods 0.000 claims description 5
- 238000002347 injection Methods 0.000 claims description 4
- 239000007924 injection Substances 0.000 claims description 4
- 230000006641 stabilisation Effects 0.000 claims description 4
- 238000011105 stabilization Methods 0.000 claims description 4
- 240000002853 Nelumbo nucifera Species 0.000 claims description 3
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 3
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 3
- 239000000284 extract Substances 0.000 claims description 3
- 238000000205 computational method Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 238000011156 evaluation Methods 0.000 abstract description 4
- 230000008901 benefit Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 230000005611 electricity Effects 0.000 description 4
- 238000003012 network analysis Methods 0.000 description 4
- 238000010248 power generation Methods 0.000 description 4
- 230000000717 retained effect Effects 0.000 description 3
- 239000000243 solution Substances 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 238000012614 Monte-Carlo sampling Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000013316 zoning Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- 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/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06312—Adjustment or analysis of established resource schedule, e.g. resource or task levelling, or dynamic rescheduling
-
- 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
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/80—Management or planning
- Y02P90/82—Energy audits or management systems therefor
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Tourism & Hospitality (AREA)
- Water Supply & Treatment (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Power Engineering (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
A kind of expectation for considering trend and sensitivity consistency equivalence lacks power supply volume appraisal procedure, considers trend and the fast reliability evaluation method of sensitivity consistency equivalence.Using computer, by program, the basic data of network and Intranet reliability assessment parameter are inputted first.Equivalent network is obtained, and calculate equivalent network parameter by the static equivalence method of trend and sensitivity consistency, establishes the topology analysis model of the Intranet reliability assessment containing equivalent network.Intranet system mode is sampled with Monte Carlo method, and the minimum tangential load amount of each sample mode is calculated based on the topology analysis model containing equivalent network, finally statistics obtains Intranet reliability index EENS.
Description
Technical field
The present invention relates to lack power supply volume appraisal procedure for electrical domain, particularly a kind of expectation for equivalent network.
Background technology
As the greatly developing of new energy, the continuous of electricity market promote and the continuous growth of load, electric system by
Gradually develop into the interconnected power grid between layering and zoning and each subregion and being closely connected.Due to influencing each other day between each subnet
Benefit enhancing must effectively examine the influence for interconnecting outer net when the subnet (i.e. Intranet) to being studied carries out analysis decision
Consider;But in some cases, due to technical reason or the problems such as be related to industry confidentiality, electricity cannot be shared between system subnet
Network data causes that complete data exchange can not be carried out in system between each subnet.Therefore, interconnected subnetworks information can be retained again
Equivalent Model with good security performance is widely used in interconnected network analysis calculates.
In the index of reliability assessment calculates, trend and sensitivity are needed for network analysis and minimum tangential load calculating
Important information should fully consider trend and the consistency of sensitivity before and after equivalence when establishing Equivalent Model, can with ensure Intranet
By property index computational accuracy and its validity of operational decisions.
It is existing to consider that equivalent reliability estimation method usually all only considered the consistency of trend before and after equivalence, but do not have
The consistency of sensitivity before and after consideration is equivalent.As in May, 1998 is sent out on IEEE Transactions on Power System
" the Application of an Adequacy Equivalent Method in Bulk Power System of table
A Reliability Evaluation " texts, this article carry out reliability to urban core power grid using abundant intensity equivalence method and comment
Estimate.Abundant intensity equivalence method be it is a kind of simple the check-ins method such as hang, it is by external network equivalent into Injection Current, and this method meets etc.
The consistency of trend before and after value, but it is difficult to ensure that the consistency of sensitivity, and do not account for the rack characteristic of outer net, in system
When changing, it is impossible to the computational accuracy of network analysis and minimum tangential load be effectively ensured, so as to influence the standard of reliability assessment
True property.
EENS (expected energy not supplied) also known as expectation lack power supply volume, refer to system to timing
Between the power generation capacity shortage of section internal cause or power constraint cause the expectation number of workload demand energy curtlaiment.EENS indexs are that system can
It is tightly related by one of property index and the probability, range, duration having a power failure, it is utilized to reflection system reliability Billy and uses
The overload of busbar or power transmission line is more accurate comprehensive.
Invention content
It is an object of the invention to provide a kind of expectations for considering trend and sensitivity consistency equivalence to lack power supply volume assessment
Method, it can be calculated in the case where considering trend and sensitivity consistency it is expected to lack power supply volume appraisal procedure.
It realizes, is as follows the purpose of the present invention is technical solution in this way:
1) interconnected network basic data and reliability assessment parameter are acquired;
2) static equivalence method based on trend and sensitivity consistency is used, the equivalent parameters of equivalent network is calculated, builds
The topology analysis model of vertical equivalent network;
3) it according to gained topology analysis model and constraints in step 2), establishes based on trend, sensitivity and constraint one
The equivalent amendment minimum load shedding model of cause property, and solve;
4) reliability assessment index it is expected that lacking power supply volume EENS calculates;
Internetwork node is divided into three parts, respectively outer net node set E, boundary node set B and Intranet node
Set I;
Interconnected network basic data described in step 1) includes interconnected network topological structure and power equipment parameter, reliably
Property assessment parameter include interior mesh element number and the availability of each element;
The power equipment parameter includes the impedance parameter and susceptance parameter over the ground of whole circuits, the impedance ginseng of transformer
It counts, over the ground admittance parameter and no-load voltage ratio parameter, the admittance parameter over the ground of whole nodes, generator output constraints, line transmission
Constraints;
Equivalent parameters described in step 2) include branch impedanceWith equivalent branch admittance over the groundDuty valueEquivalent generator output SeqGi;
Minimum load shedding model is corrected described in step 3) is
In formula, CkFor Intranet or the cutting load amount of boundary node k, NI,NBRespectively Intranet and boundary node quantity;
Constraints described in step 3) includes:
3-1) trend Constraints of Equilibrium
In formula:I ∈ (I, B, eq), PGiIt is exported for generated power, PDiFor the load active power of node i, QGiFor power generation
Machine is idle output, QDiFor the reactive load power of node i, QcriInjection reactive power, V are compensated for parallel reactiveiFor external node
On high-tension side voltage, VjFor the voltage of external node low-pressure side, GijFor the conductance between external node i and j, BijFor external node
Susceptance between i and j, δijFor the phase angle difference between external node i and j, N be it is equivalent after the whole network number of nodes, NBIt is saved for boundary
Point quantity, NSFor balance nodes;
3-2) Intranet variable bound
kItmin≤kIt≤kItmaxT=1 ..., NT (14)
PIGimin≤PIGi≤PIGimaxI=1 ..., NG (15)
QIGimin≤QIGi≤QIGimaxI=1 ..., NG (16)
QIcrimin≤QIcri≤QIcrimaxI=1 ..., Ncr (17)
In formula (4)~(7), kItFor Intranet on-load regulator transformer LTC no-load voltage ratios, kItmin,kItmaxFor kItBound about
Beam, NTNumber for Intranet on-load regulator transformer LTC;PIGiAnd QIGiActive and idle, P is exported for Intranet generatorIGimin,
PIGimaxAnd QIGimin,QIGimaxIt is P respectivelyIGiAnd QIGiBound constraint, NGNumber for Intranet generator;QIcriFor Intranet simultaneously
Join reactive-load compensation, QIcrimin,QIcrimaxFor QIcriBound constraint, NcrThe number of equipment is compensated for parallel reactive;
3-3) node burden with power abatement amount constrains
0≤Ck≤PDkK=1 ..., NI+NB (18)
So far, it is established and completed based on the amendment minimum load shedding model of trend and sensitivity consistency equivalence, and using pre-
Prim al- dual interior point m ethod is surveyed to solve the minimum load shedding model.
Further, the computational methods of EENS described in step 4) are as follows:
4-1) select Intranet system mode
It is sampled using Monte Carlo method to system mode, if sampling total degree is M times, what is extracted every time is all one complete
Independent state, then the system mode probability extracted each time are as follows:
4-2) topological analysis
After system mode has been selected, topological analysis is carried out to system, the node for having electrical link and branch are divided into
In one subsystem, and the data of each subsystem are sorted out, to form the input data of voltage stabilization Corrective control analysis;If
In the case of system not off-the-line, the whole network is a system, if for the system of off-the-line, gets rid of the network of smaller portions, is retained
Big network;
4-3) total EENS indexs calculate;
In formula:S is the system mode complete or collected works that cannot meet workload demand in given time section;PiState is in for system
The probability of i;CiFor the load power cut down under the conditions of state i;T is 8760h;It is expected that it is MWh/a to lack power supply volume unit.
Further, step 4-2) after the whole network is a system or the big network of reservation, it is further comprising the steps of:
A, trend distribution is calculated, discriminates whether out-of-limit, if out-of-limit, is transferred to step B, step E is transferred to if not out-of-limit;
B, it using adjustment means, discriminates whether out-of-limit, step C is transferred to if out-of-limit, if not out-of-limit be transferred to step E;
C, using optimization algorithm carry out minimum cut lotus scheme processing, setting Intranet it is all can cutting load node total cutting load
Minimum optimization aim is measured, mesh element constrains and network power Constraints of Equilibrium constraints after equivalence, foundation optimization mould in consideration
Type;Then Optimized model is solved using prediction prim al- dual interior point m ethod, judges whether minimum tangential load scheme restrains and scheme is
It is no unique, if being then transferred to step D, if not ignoring this system mode situation and being transferred to step E;
D, the EENS indexs of computing system state X, and it is transferred to step E;
E, judge whether to meet the Monte Carlo method condition of convergence, if satisfied, being then transferred to step F;
F, it is transferred to step 4-3).
By adopting the above-described technical solution, the present invention has the advantage that:
1st, in terms of the Intranet reliability assessment topology analysis model containing equivalent network is established, outer net of the present invention
Equivalence method can not only keep the consistency of trend before and after equivalence, moreover it is possible to sensitivity consistency before and after equivalence be effectively ensured, effectively
Improve the accuracy of the Intranet reliability assessment topology analysis model containing equivalent network;
2nd, compared to the existing reliability estimation method based on on-hook Theory of Equivalence and Ward Theory of Equivalence, the present invention is built
Vertical topology analysis model can be effectively retained power of the outer net between the power, voltage support and boundary node of Intranet and turn
The characteristics such as shifting so that the analyses such as Load flow calculation, minimum tangential load for system mode are more accurate, so as to ensure that reliability refers to
Target validity and accuracy.
Other advantages, target and the feature of the present invention will be illustrated in the following description to a certain extent, and
And to a certain extent, based on will be apparent to those skilled in the art to investigating hereafter, Huo Zheke
To be instructed from the practice of the present invention.The target and other advantages of the present invention can be wanted by following specification and right
Book is sought to realize and obtain.
Description of the drawings
The description of the drawings of the present invention is as follows.
Fig. 1 is the flow diagram of the present invention;
Fig. 2 is step 2) medium value circuit diagram;
Fig. 3 is node standard test system figure.
Specific embodiment
The invention will be further described with reference to the accompanying drawings and examples.
The present invention was using the 5th phase of volume 29 in 2014《IEEE TRANSACTIONS ON POWER SYSTEM》In
“Static Equivalent Method Based on Component Particularity Representation and
The static equivalence method of trend and sensitivity consistency before and after the holding equivalence that Sensitivity Consistency " are proposed
Equivalent network and its parameter are obtained, establishes the topology analysis model of the Intranet reliability assessment containing equivalent network.With Meng Teka
Lip river method is sampled Intranet system mode, calculates and obtains accurate reliability index.
Consider trend and the fast reliability evaluation method of sensitivity consistency equivalence.It is first by program using computer
First input the basic data of network and Intranet reliability assessment parameter.Pass through trend and the static equivalence method of sensitivity consistency
Equivalent network is obtained, and calculates equivalent network parameter, establishes the topology analysis model of the Intranet reliability assessment containing equivalent network.
Intranet system mode is sampled with Monte Carlo method, and each take out is calculated based on the topology analysis model containing equivalent network
The minimum tangential load amount of sample state finally counts and obtains Intranet reliability index EENS.Its specific method step is as follows:
1) basic data and reliability assessment parameter are inputted
Interconnected network basic data is inputted, including interconnected network topological structure, i.e., is respectively saved in sub-area division situation and power grid
The connection relation of point and power equipment parameter, i.e., the impedance parameter and susceptance parameter over the ground of whole circuits, the impedance ginseng of transformer
It counts, over the ground admittance parameter and no-load voltage ratio parameter, the admittance parameter over the ground of whole nodes, generator output constraints, line transmission
Constraints etc..Internetwork node before equivalence is divided into three parts, respectively outer net node set E, boundary node set B
With Intranet node set I.
Reliability assessment parameter is inputted, including interior mesh element number, i.e., total element number and available components number and each member
The availability of part, you can with the availability of element and the degree of unavailability of unavailable element.
2) topology analysis model containing equivalent network is established
Using the existing static equivalence method based on trend and sensitivity consistency, i.e. " Static Equivalent
Method Based on Component Particularity Representation and Sensitivity
Consistency " calculates the equivalent parameters of medium value network of the present invention, including equivalent branch impedanceWith equivalent branch admittance over the groundDuty valueEquivalent generator output SeqGi。
After calculating equivalent parameters, equivalent network figure as shown in Figure 2 can be drawn, so far opening up containing equivalent network
Analysis model is flutterred just to be successfully established.
Branch in figure between 3 nodes and 4 nodes disconnects, and the branch between 8 nodes and 9 nodes disconnects.22 nodes and 29
Node is boundary node;4-8,18-21 node are external node;1-3,9-17,23-28 and 30-57 node are internal node,
Its interior joint 1 is balance nodes.
3) it establishes and corrects minimum load shedding model
It is to include the electricity of boundary node with reference to Intranet based on the topology analysis model containing equivalent network that step 2) step obtains
The constraintss such as pressure, generator and line power are established the amendment minimum based on trend, sensitivity and constraint consistency equivalence and are cut
Load model.
1) object function
Since outer net is substituted by equivalent network, therefore the object function of minimum tangential load amount is represented by after equivalence;
In formula, CkFor Intranet or the cutting load amount of boundary node k, NI,NBRespectively Intranet and boundary node quantity.
2) constraints
A. trend Constraints of Equilibrium
In formula:I ∈ (I, B, eq), PGiIt is exported for generated power, PDiFor the load active power of node i, QGiFor power generation
Machine is idle output, QDiFor the reactive load power of node i, QcriInjection reactive power, V are compensated for parallel reactiveiFor external node
On high-tension side voltage, VjFor the voltage of external node low-pressure side, GijFor the conductance between external node i and j, BijFor external node
Susceptance between i and j, δijFor the phase angle difference between external node i and j, N be it is equivalent after the whole network number of nodes, NBIt is saved for boundary
Point quantity, NSFor balance nodes.
B. Intranet variable bound
kItmin≤kIt≤kItmaxT=1 ..., NT (24)
PIGimin≤PIGi≤PIGimaxI=1 ..., NG (25)
QIGimin≤QIGi≤QIGimaxI=1 ..., NG (26)
QIcrimin≤QIcri≤QIcrimaxI=1 ..., Ncr (27)
In formula (4)~(7), kItFor Intranet on-load regulator transformer LTC no-load voltage ratios, kItmin,kItmaxFor kItBound about
Beam, NTNumber for Intranet on-load regulator transformer LTC;PIGiAnd QIGiActive and idle, P is exported for Intranet generatorIGimin,
PIGimaxAnd QIGimin,QIGimaxIt is P respectivelyIGiAnd QIGiBound constraint, NGNumber for Intranet generator;QIcriFor Intranet simultaneously
Join reactive-load compensation, QIcrimin,QIcrimaxFor QIcriBound constraint, NcrThe number of equipment is compensated for parallel reactive.
C. node burden with power reduction constrains
0≤Ck≤PDkK=1 ..., NI+NB (28)
So far, it is established and completed based on the amendment minimum load shedding model of trend and sensitivity consistency equivalence.The present invention adopts
With prediction prim al- dual interior point m ethod, the minimum load shedding model is solved.
4) reliability assessment index calculates
After the work for completing above 3 steps, reliability assessment index calculating can be carried out, the present invention calculates reliable
Property index be EENS.EENS also known as expectation lack power supply volume, refer to system in the power generation capacity shortage of given time section internal cause or electricity
Net constraint causes the expectation number of workload demand energy curtlaiment.
4-1) select Intranet system mode
The present invention samples to system mode using Monte Carlo method, and basic ideas are a series of with random number sequence generation
Experiment sample.When sample size is sufficiently large, according to central-limit theorem or the law of large numbers, sample average can be used as the mathematics phase
The unbiased esti-mator of prestige.The variance of sample average is a mark of estimated accuracy.If total degree of sampling is M times, extract every time
All it is a completely self-contained state, then the system mode probability extracted each time is as follows:
In the sampling process of Monte Carlo, load is also required to random sampling, though the place of system jam is identical, but
It is the fluctuation due to load, the state extracted every time is also not quite similar.Therefore in the present invention, it is assumed that the state extracted every time all
It is completely independent, then the probability of each system mode is the same.
4-2) topological analysis
After system mode has been selected, it is necessary to topological analysis is carried out to system, task seeks to analysis whole system
Node how many a subsystems are connected by line and transformer, will there is the node of electrical link and branch to be divided into a son
In system, and the related data of each subsystem is sorted out, to form the input data of voltage stabilization Corrective control analysis.It is being
For system not in the case of off-the-line, the whole network is a system.For the system of off-the-line, the network of smaller portions is got rid of, retains big net
Network.
4-3) EENS indexs calculate
EENS indexs are that one of Reliability Index and the probability having a power failure, range, duration are tightly related, are utilized
It come reflect the overload of system reliability Billy busbar or power transmission line more it is accurate comprehensively.Reliability Index EENS's
Calculation formula is as follows:
In formula:S is the system mode complete or collected works that cannot meet workload demand in given time section;PiState is in for system
The probability of i;CiFor the load power cut down under the conditions of state i;T is 8760h.It is expected that it is MWh/a to lack power supply volume unit.
Since EENS indexs are a summation amount, if the two systems state extracted is identical, minimum tangential load amount
Also it is identical, therefore the minimum tangential load amount of the two states can be merged and be cut as the minimum under a system mode automatically
Load is calculated.
Embodiment:
(1) basic data and reliability assessment parameter are inputted
Internetwork node before equivalence is divided into three parts, respectively outer net node set E={ 4-8,18-21 }, side
Boundary node set B={ 22,29 } and Intranet node set I={ 1-3,9-17,23-28,30-57 }.It is shown in input Fig. 3 first
The basic data of system, wherein, interconnected network basic data include topological structure and power equipment parameter, reference capacity 100MW,
Referring primarily to IEEE-57 modular systems.
Secondly, reliability assessment parameter is inputted, the available components number of IEEE57 node modular systems is 70, wherein wrapping
Include interior mesh element number and the number of boundary non-equivalent element.The availability of each available components is as shown in the table:
The availability table of 1 system element of table
Element | Generator | Circuit | Transformer | Shunt capacitance | Parallel reactance |
Availability | 0.008 | 0.01 | 0.01 | 0.008 | 0.008 |
(2) topology analysis model containing equivalent network is established
After the completion of (1) step, in this embodiment, using the existing Static Equivalent side based on trend and sensitivity consistency
Method, i.e. " Static Equivalent Method Based on Component Particularity
Representation and Sensitivity Consistency " calculate the equivalent parameters of medium value network of the present invention, packet
Include equivalent branch impedance Zeq22、Zeq29And ZEq58,59, equivalent branch admittance over the groundWithDuty value PeqL22And PeqL29。
Equivalent parameters are as shown in the table:
Table 2 (p.u) external network equivalent parameter list
(3) it establishes and corrects minimum load shedding model
Based on the topology analysis model containing equivalent network that (2) step obtains, voltage with reference to Intranet including boundary node,
The constraintss such as generator and line power, establish based on trend, sensitivity and constrain consistency equivalence amendment minimum cut it is negative
Lotus model.
Consider trend Constraints of Equilibrium formula (2) and formula (3), Intranet variable bound formula (4), (5), (6) and (7) and
The minimum tangential load amount of system can be calculated according to formula (1) for node burden with power reduction constraints formula (8).
It using new outer net static equivalence method in the network of 57 nodes, is sampled using Monte Carlo method, various states can be calculated
Under minimum tangential load amount, the minimum tangential load amount of which part state is:
3 (unit of table:p.u)
(4) reliability index calculates
After the work for completing above 3 steps, reliability index calculating can be carried out, the reliability that the present invention calculates refers to
It is designated as it is expected to lack power supply volume EENS indexs.
1) Intranet system mode is selected
The present invention samples to system mode using Monte Carlo method, and basic ideas are a series of with random number sequence generation
Experiment sample.When sample size is sufficiently large, according to central-limit theorem or the law of large numbers, sample average can be used as the mathematics phase
The unbiased esti-mator of prestige.Monte Carlo sampling number is set as 15000 times in the present invention, and that extracts every time is all one and is completely independent
State, then the system mode probability extracted each time is:
In the sampling process of Monte Carlo, load is also required to random sampling, though the place of system jam is identical, but
It is the fluctuation due to load, the state extracted every time is also not quite similar.Therefore in the present invention, it is assumed that the state extracted every time all
It is completely independent, then the probability of each system mode is the same.
2) topological analysis
After system mode has been selected, it is necessary to topological analysis is carried out to system, task seeks to analysis whole system
Node how many a subsystems are connected by line and transformer, will there is the node of electrical link and branch to be divided into a son
In system, and the related data of each subsystem is sorted out, to form the input data of voltage stabilization Corrective control analysis.It is being
For system not in the case of off-the-line, the whole network is a system.For the system of off-the-line, the network of smaller portions is got rid of, retains big net
Network.
3) EENS indexs calculate
EENS indexs are that one of Reliability Index and the probability having a power failure, range, duration are tightly related, are utilized
It come reflect the overload of system reliability Billy busbar or power transmission line more it is accurate comprehensively.Reliability Index EENS's
Calculation formula is as follows:
It is as follows to participate in the several method compared:
M0:Original the whole network calculates;
M1:The method of the present invention;
M2:Check-ins method --- the most common equivalence methods in engineering such as the simple extension of outer net constraint are not considered;
M3:Equivalence method based on ward equivalences --- consider the equivalence method of outer net constraint.
Now define two indexs that can reflect external network equivalent precision, absolute error e1And relative error e2, the two indexs
Numerical value it is smaller, illustrate that precision is higher.
Absolute error e1:
e1=| x-xeq|
Relative error e2:
Wherein, x and xeqTrue value and estimated value are represented respectively, and so-called true value refers under simulated conditions, inside and outside using detailed
The reliability index calculated value of refined net model, and estimated value refers in Intranet using detailed network model, outer net is using equivalent
Reliability index calculated value under network model.
It is 1. as shown in table 6 using the EENS indexs that method M1-M3 is calculated.
The EENS indexs that 6 method M0-M3 of table is calculated
Equivalence method | M0 | M1 | M2 | M3 |
EENS results (MWh/a) | 23545 | 22640 | 19288 | 25121 |
It is 2. as shown in table 7 using the EENS scale errors analysis that method M1-M3 is calculated.
The EENS scale errors analysis that 7 method M1-M3 of table is calculated
It was found from experimental result:Use the fast reliability proposed by the present invention based on trend and sensitivity consistency equivalence
After appraisal procedure, the error calculation result of Reliability Index is all very small.
In conclusion the fast reliability evaluation side based on trend and sensitivity consistency equivalence that the method for the present invention proposes
Method computational accuracy be higher than it is existing do not consider sensitivity consistency containing equivalent reliability index calculating method, being capable of mould well
Intend the operating condition of practical outer net, appropriate power and voltage support are provided to Intranet, to ensure the safety and stability of interconnected network
Economical operation.
Finally illustrate, the above embodiments are merely illustrative of the technical solutions of the present invention and it is unrestricted, although with reference to compared with
The present invention is described in detail in good embodiment, it will be understood by those of ordinary skill in the art that, it can be to the skill of the present invention
Art scheme is modified or replaced equivalently, and without departing from the objective and range of the technical program, should all be covered in the present invention
Right in.
Claims (3)
1. a kind of expectation for considering trend and sensitivity consistency equivalence lacks power supply volume appraisal procedure, which is characterized in that specific step
It is rapid as follows:
1) interconnected network basic data and reliability assessment parameter are acquired;
2) static equivalence method based on trend and sensitivity consistency is used, calculates the equivalent parameters of equivalent network, establish etc.
It is worth the topology analysis model of network;
3) it according to gained topology analysis model and constraints in step 2), establishes based on trend, sensitivity and constraint consistency
Equivalent amendment minimum load shedding model, and solve;
4) reliability assessment index it is expected that lacking power supply volume EENS calculates;
Internetwork node is divided into three parts, respectively outer net node set E, boundary node set B and Intranet node set
I;
Interconnected network basic data described in step 1) includes interconnected network topological structure and power equipment parameter, and reliability is commented
Estimate parameter and include interior mesh element number and the availability of each element;
The power equipment parameter includes the impedance parameters of whole circuits and susceptance parameter over the ground, the impedance parameter of transformer,
Admittance parameter and no-load voltage ratio parameter over the ground, the admittance parameter over the ground of whole nodes, generator output constraints, line transmission constraint
Condition;
Equivalent parameters described in step 2) include branch impedanceWith equivalent branch admittance over the groundDuty valueEquivalent generator output SeqGi;
Minimum load shedding model is corrected described in step 3) is
In formula, CkFor Intranet or the cutting load amount of boundary node k, NI,NBRespectively Intranet and boundary node quantity;
Constraints described in step 3) includes:
3-1) trend Constraints of Equilibrium
In formula:I ∈ (I, B, eq), PGiIt is exported for generated power, PDiFor the load active power of node i, QGiFor generator without
Work(exports, QDiFor the reactive load power of node i, QcriInjection reactive power, V are compensated for parallel reactiveiFor external node high pressure
The voltage of side, VjFor the voltage of external node low-pressure side, GijFor the conductance between external node i and j, BijFor external node i and j
Between susceptance, δijFor the phase angle difference between external node i and j, N be it is equivalent after the whole network number of nodes, NBFor boundary node number
Amount, NSFor balance nodes;
3-2) Intranet variable bound
kItmin≤kIt≤kItmaxT=1 ..., NT (4)
PIGimin≤PIGi≤PIGimaxI=1 ..., NG (5)
QIGimin≤QIGi≤QIGimaxI=1 ..., NG (6)
QIcrimin≤QIcri≤QIcrimaxI=1 ..., Ncr (7)
In formula (4)~(7), kItFor Intranet on-load regulator transformer LTC no-load voltage ratios, kItmin,kItmaxFor kItBound constraint, NT
Number for Intranet on-load regulator transformer LTC;PIGiAnd QIGiActive and idle, P is exported for Intranet generatorIGimin,PIGimax
And QIGimin,QIGimaxIt is P respectivelyIGiAnd QIGiBound constraint, NGNumber for Intranet generator;QIcriFor Intranet parallel connection nothing
Work(compensates, QIcrimin,QIcrimaxFor QIcriBound constraint, NcrThe number of equipment is compensated for parallel reactive;
3-3) node burden with power abatement amount constrains
0≤Ck≤PDkK=1 ..., NI+NB (8)
So far, it is established and completed based on the amendment minimum load shedding model of trend and sensitivity consistency equivalence, and is former using prediction
Dual interior point solves the minimum load shedding model.
2. consider that the expectation of trend and sensitivity consistency equivalence lacks power supply volume appraisal procedure as described in claim 1,
It is characterized in that, the computational methods of EENS described in step 4) are as follows:
4-1) select Intranet system mode
It is sampled using Monte Carlo method to system mode, if sampling total degree is M times, that extracts every time is all one and is completely independent
State, then the system mode probability extracted each time is as follows:
4-2) topological analysis
After system mode has been selected, topological analysis is carried out to system, the node for having electrical link and branch are divided into one
In subsystem, and the data of each subsystem are sorted out, to form the input data of voltage stabilization Corrective control analysis;If being
For system not in the case of off-the-line, the whole network is a system, if for the system of off-the-line, gets rid of the network of smaller portions, retains big
Network;
4-3) total EENS indexs calculate;
In formula:S is the system mode complete or collected works that cannot meet workload demand in given time section;PiThe general of state i is in for system
Rate;CiFor the load power cut down under the conditions of state i;T is 8760h;It is expected that it is MWh/a to lack power supply volume unit.
3. consider that the expectation of trend and sensitivity consistency equivalence lacks power supply volume appraisal procedure as claimed in claim 2, it is special
Sign is, step 4-2) in the whole network for a system or after retaining big network, it is further comprising the steps of:
A, trend distribution is calculated, discriminates whether out-of-limit, if out-of-limit, is transferred to step B, step E is transferred to if not out-of-limit;
B, it using adjustment means, discriminates whether out-of-limit, step C is transferred to if out-of-limit, if not out-of-limit be transferred to step E;
C, using optimization algorithm carry out minimum cut lotus scheme processing, setting Intranet it is all can cutting load node total cutting load amount most
It is small for optimization aim, in consideration mesh element constraint and it is equivalent after network power Constraints of Equilibrium constraints, establish Optimized model;So
Afterwards using prediction prim al- dual interior point m ethod to Optimized model solve, judge minimum tangential load scheme whether restrain and scheme whether only
One, if being then transferred to step D, if not ignoring this system mode situation and being transferred to step E;
D, the EENS indexs of computing system state X, and it is transferred to step E;
E, judge whether to meet the Monte Carlo method condition of convergence, if satisfied, being then transferred to step F;
F, it is transferred to step 4-3).
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610127267.6A CN105680442B (en) | 2016-03-07 | 2016-03-07 | Consider that the expectation of trend and sensitivity consistency equivalence lacks power supply volume appraisal procedure |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610127267.6A CN105680442B (en) | 2016-03-07 | 2016-03-07 | Consider that the expectation of trend and sensitivity consistency equivalence lacks power supply volume appraisal procedure |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105680442A CN105680442A (en) | 2016-06-15 |
CN105680442B true CN105680442B (en) | 2018-06-15 |
Family
ID=56306961
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610127267.6A Active CN105680442B (en) | 2016-03-07 | 2016-03-07 | Consider that the expectation of trend and sensitivity consistency equivalence lacks power supply volume appraisal procedure |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105680442B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106887837A (en) * | 2017-02-16 | 2017-06-23 | 广西电网有限责任公司电力科学研究院 | For the substance line fault load cutting method of reliability assessment |
CN107370149B (en) * | 2017-08-25 | 2020-04-14 | 广州供电局有限公司 | Method and system for determining load reduction amount of power system |
CN109066685B (en) * | 2018-08-02 | 2021-10-01 | 国网安徽省电力有限公司 | Line parameter correction method based on parameter sensitivity |
CN109472509B (en) * | 2018-11-22 | 2020-09-29 | 重庆大学 | High-voltage direct-current transmission system reliability interval assessment method |
CN109873421A (en) * | 2019-03-29 | 2019-06-11 | 国网浙江省电力有限公司 | A kind of provincial power network abundant intensity online evaluation method and system |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103646358A (en) * | 2013-12-30 | 2014-03-19 | 重庆大学 | Method for detecting planned maintenance cycle of power grid related to power equipment time-varying failure rate |
CN104715423A (en) * | 2015-03-13 | 2015-06-17 | 国家电网公司 | Method for assessing risk and reliability of power distribution network |
CN105071381A (en) * | 2015-07-28 | 2015-11-18 | 天津大学 | State enumeration reliability evaluation method and device based on influence increment |
CN105226644A (en) * | 2015-09-23 | 2016-01-06 | 重庆大学 | Based on the conforming belt restraining equivalence method of active volume |
-
2016
- 2016-03-07 CN CN201610127267.6A patent/CN105680442B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103646358A (en) * | 2013-12-30 | 2014-03-19 | 重庆大学 | Method for detecting planned maintenance cycle of power grid related to power equipment time-varying failure rate |
CN104715423A (en) * | 2015-03-13 | 2015-06-17 | 国家电网公司 | Method for assessing risk and reliability of power distribution network |
CN105071381A (en) * | 2015-07-28 | 2015-11-18 | 天津大学 | State enumeration reliability evaluation method and device based on influence increment |
CN105226644A (en) * | 2015-09-23 | 2016-01-06 | 重庆大学 | Based on the conforming belt restraining equivalence method of active volume |
Non-Patent Citations (1)
Title |
---|
电力元件和系统的运行可靠性模型研究;杨洋;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20100215(第02期);第C042-122页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105680442A (en) | 2016-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105680442B (en) | Consider that the expectation of trend and sensitivity consistency equivalence lacks power supply volume appraisal procedure | |
Balamurugan et al. | Review of power flow studies on distribution network with distributed generation | |
Su et al. | Optimal PV inverter reactive power control and real power curtailment to improve performance of unbalanced four-wire LV distribution networks | |
CN106329523B (en) | Consider probabilistic active power distribution network intelligence Sofe Switch robust Optimization Modeling method | |
CN103872681A (en) | Online real-time loop closing method based on integration of major network and distribution network | |
US20140244189A1 (en) | System For Achieving Real-Time Monitoring and State Estimation in Power Distribution Networks | |
CN107342586B (en) | Method and device for determining static equivalence of power distribution network | |
CN103236691A (en) | Method of three-phase unbalance load flow calculation based on complex affine mathematical theory | |
CN105656036B (en) | Consider trend and the probability static security analysis method of sensitivity uniformity equivalence | |
CN108649597B (en) | Fault migration method and analysis method for influence of fault on HVDC commutation failure | |
CN110162843B (en) | Real-time simulation modeling method and device for power grid primary system and secondary system | |
Chen et al. | Distribution system state estimation: A survey of some relevant work | |
CN112103950B (en) | Power grid partitioning method based on improved GN splitting algorithm | |
De Albuquerque et al. | Nonlinear analysis on transmission line parameters estimation from noisy phasorial measurements | |
CN105790262B (en) | A kind of method and system of the EENS equivalent calculation based on enumerative technique | |
KR101988670B1 (en) | Device for state estimation of power distribution system | |
Andreoli et al. | Topology validation via simultaneous state & topology estimation with phasor data processing capability | |
CN113139295A (en) | Method and system for estimating comprehensive state of power system | |
CN105226644B (en) | Belt restraining equivalence method based on active volume uniformity | |
CN111834996B (en) | Power grid line loss calculation method and device | |
CN115549093B (en) | Method and system for online modeling and oscillation analysis of new energy power system | |
Parihar et al. | Power flow analysis of balanced radial distribution system with composite load model | |
CN105846447B (en) | Multiple-circuit line Inverter Station dynamic passive compensation configures the determination method of order of priority | |
CN114996635A (en) | Power distribution network parameter determination method, device, equipment and storage medium | |
CN112688332B (en) | Hybrid alternating current-direct current microgrid power flow analysis method and system based on sequence component coordinates |
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 |