CN107171334B - Event driven cutting load and corrective route switching control method for coordinating and system - Google Patents
Event driven cutting load and corrective route switching control method for coordinating and system Download PDFInfo
- Publication number
- CN107171334B CN107171334B CN201710491138.XA CN201710491138A CN107171334B CN 107171334 B CN107171334 B CN 107171334B CN 201710491138 A CN201710491138 A CN 201710491138A CN 107171334 B CN107171334 B CN 107171334B
- Authority
- CN
- China
- Prior art keywords
- cutting load
- route switching
- corrective
- event driven
- transient
- 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
Classifications
-
- 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
- Y02B—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO BUILDINGS, e.g. HOUSING, HOUSE APPLIANCES OR RELATED END-USER APPLICATIONS
- Y02B70/00—Technologies for an efficient end-user side electric power management and consumption
- Y02B70/30—Systems integrating technologies related to power network operation and communication or information technologies for improving the carbon footprint of the management of residential or tertiary loads, i.e. smart grids as climate change mitigation technology in the buildings sector, including also the last stages of power distribution and the control, monitoring or operating management systems at local level
- Y02B70/3225—Demand response systems, e.g. load shedding, peak shaving
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S20/00—Management or operation of end-user stationary applications or the last stages of power distribution; Controlling, monitoring or operating thereof
- Y04S20/20—End-user application control systems
- Y04S20/222—Demand response systems, e.g. load shedding, peak shaving
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
Abstract
The invention discloses event driven cutting loads and corrective route switching control method for coordinating and system, wherein the method comprising the steps of (1): acquisition actual electric network data build the simulation model of studied power grid；Relevant parameter is set；Step (2): selection forecast accident executes transient emulation and obtains monitoring data, calculate corresponding security margin index and compared with default safety margin critical value when not taking control measure, judges whether transient security has no branch overload to system；Step (3): Optimized model is set up by the forecast accident selected and the parameter of setting, it is iterated solution, calculate cutting load dominant vector and route switching dominant vector, in conjunction with the position of each cutting load node and switching route in power grid, the control program for being directed to current forecast accident is formulated；Step (4): in grid collapses, corresponding control program is executed, and Starting mode is to detect that failure generation then starts immediately；Return step (1) at regular intervals.
Description
Technical field
The invention belongs to power system security stability contorting field more particularly to a kind of event driven cutting loads and correction
Property route switching control method for coordinating and system.
Background technique
Under the driving of energy-saving and emission-reduction national policy, the renewable and clean energy resources such as wind-powered electricity generation, photovoltaic are quickly grown, in electricity
Permeability in net is continuously improved.The power supply of these low inertias largely accesses electric system, by the equivalent inertial time for the system that reduces
Between constant, aggravation power shortage disturbance when frequency decline degree.Since China load consumer center is in energy resources center
Contrary distribution, remote, the high capacity transmission of electric energy become a kind of necessary measure.As a plurality of extra-high voltage AC and DC is defeated
The construction of electric line, being continuously improved by electric ratio for east China area, increases it by potential power disturbance quantity.It is high-power
Vacancy disturbs transient frequency safety, the transient voltage security that can seriously affect electric system, and can cause transmission line of electricity, transformer
Equal electric equipment overloads.If effective control measure can not be taken, accident impact is limited, system safety operation will be seriously threatened.
Common control measure of the cutting load as reply system wide range of frequencies, variation and overload, are power grid securities
The important component of stability contorting.The transient frequency and voltage security of electric system need the transient process after consideration disturbance,
Being related to time scale usually is all second grade；And overload and be related to middle long time scale, generally reach minute grade even hour grade.But
It is that current cutting load measure usually separately studies it, lacks a unified research control framework, be unable to satisfy power train
The requirement for emergent control of uniting, may cause the control measure for a certain security feature, can deteriorate other two aspect instead, lead
Cause the negative effect of control.
It is negative that the 599th command of State Council " disposition of electric power safety accident emergency and regulations of investigating " specifies that stabilized control system is cut
Lotus is equal to breakdown loss load, region cutting load ratio is excessively high or load distribution it is unreasonable will lead to it is even more serious
Incident classification grading and accident are called to account.But current cutting load control measure are mostly based on response driving, such as low-frequency load shedding and
Low-voltage load sheding.At failure initial stage, system frequency, Voltage Drop are unobvious, it is difficult to reach the movement threshold value of cutting load device, this meeting
Deteriorate system running state, and needs to cut off more loads after causing to maintain system safe.
Currently, although corrective route switching can be changed and be closed by breaker actuation switching route, converting network topology
The trend of key transmission cross-section is distributed, thus the loading condition of heavy-haul line or transformer before reducing.Corrective route switching is only
Existing breaker is relied on, hardware investment and maintenance cost are very low, and can cutting load amount after reduction accident.But corrective at present
Route switching is mainly individually controlled, and the association with other control measure (such as cutting load, direct-current emergency are modulated) is lacked
It adjusts.For catastrophe failure, it is possible to occur relying only on the case where corrective route switching can not eliminate overload, and only consider to rectify
Influence of the positivity route switching to overload, does not consider whether it can deteriorate system transient modelling frequency and voltage security, may cause control
The negative effect of system.
In conclusion it is that power system security stability contorting is put forward new requirements that electric network composition transformation and new regulations, which are put into effect,
And challenge, existing system security and stability control technology have certain limitation.Therefore, a kind of electricity net safety stable control is needed
Method processed comprehensively considers system transient modelling frequency security, transient voltage security and electric equipment overload, guarantees power grid security, subtracts
Few load loss reduces control cost.
Summary of the invention
In order to solve the deficiencies in the prior art, the present invention provides a kind of event driven cutting loads and corrective route to throw
Control method for coordinating is cut, control uniformly can be optimized to event driven cutting load and corrective route switching, entirely
Consider to face system transient modelling frequency security, transient voltage security and electric equipment overload, while guaranteeing electricity net safety stable,
Cutting load amount is reduced, control cost is reduced, improves the economy of control measure.
A kind of event driven cutting load and corrective route switching control method for coordinating of the invention, comprising:
Step (1): acquisition actual electric network data build the simulation model of studied power grid；Cutting load, route switching are set
Control parameter, simulation parameter, monitoring data, transient security constrained parameters and overload constrained parameters；
Step (2): selection forecast accident executes transient emulation and obtains monitoring data, calculate when not taking control measure
Corresponding security margin index and compared with default safety margin critical value judges system whether transient security has no branch to pass by
It carries, if it is not, thening follow the steps (3)；If so, terminating this calculating, determines and do not need to take control for current forecast accident
Measure further determines whether that all forecast accidents all handle completion, if so, thening follow the steps (4)；Otherwise, return step
(2)；
Step (3): Optimized model is set up by the forecast accident selected and the parameter of setting, is iterated solution, calculating is cut
Load control system vector sum route switching dominant vector, in conjunction with the position of each cutting load node and switching route in power grid, system
Surely it is directed to the control program of current forecast accident；Further determine whether that all forecast accidents all handle completion, if so, executing
Step (4)；Otherwise, return step (2)；
Step (4): in grid collapses, corresponding control program is executed, and Starting mode is to detect failure
Then start immediately；At regular intervals, it is executed since step (1) again, updates the corresponding control program of each forecast accident.
Further, Optimized model is set up by the forecast accident selected and the parameter of setting in the step (3) are as follows:
The minimum value of event driven cutting load is the sum of the excision load of cutting load node of all determinations；
Wherein, parameter sets up the constraint condition of Optimized model are as follows: system meets transient security respectively and has no branch overload, line
The constraint of road switching and the constraint of cutting load amount.
Further, in the step (3), Optimized model is set up to solve parameter, is broken down into two submodels,
And using the double-deck alternative manner, it is divided into internal layer iteration and external iteration, internal layer carries out event driven cutting load dominant vector
Iteration, outer layer are iterated correcting property route switching dominant vector.
Further, in the step (3), Optimized model is set up to solve parameter, is broken down into two submodels
Respectively submodel one and submodel two；Wherein, submodel first is that consider by cutting load guarantee power grid transient frequency safety with
Transient voltage security, and event driven cutting load is optimized；
Submodel is second is that consider thoroughly to eliminate overload by cutting load by route switching limiting overload degree, and to thing
The driving cutting load of part and corrective route switching optimize.
Further, in the step (3), the specific mistake that parameter sets up Optimized model is solved using the double-deck alternative manner
Journey are as follows:
Step (3-1): parameter initialization, and former Optimized model is decomposed into two submodels；
Step (3-2): a kind of corrective route switching scheme is selected；
Step (3-3): under selected corrective route switching scheme, being iterated solution to two submodels, calculates
Current optimal cutting load dominant vector and corresponding institute's cutting load total amount；
Step (3-4): if having completed to choose corresponding minimum tangential load to the search of all corrective route switchings
The scheme of total amount is as final scheme；Otherwise, return step (3-2).
Further, in the step (3-3), submodel one is solved using linearization technique, use gradually adds up
Method solves submodel two.
Further, the cutting load control section of submodel one and submodel two is successively carried out using iterative method is integrated
Iteration seeks the optimal of event driven cutting load.
Further, in the step (3), the control program of formulation include the load that each cutting load node is cut and
Route switching strategy.
Further, in the step (4), control program Starting mode is that the failure of monitoring starts immediately when occurring, side
After case starting, the control of event driven cutting load cuts off corresponding load immediately, and corrective route switching will be in system dynamic mistake
Cheng Jiben terminates, fluctuates relatively steady rear implementation, in order to avoid bigger disturbance is caused to system due to the variation of network topology.
Further, in the step (4), using off-line calculation, the control strategy of On-line matching is right at regular intervals
Control strategy is recalculated.
The present invention also provides a kind of event driven cutting loads and corrective route switching coordinated control system.
Event driven cutting load and corrective route switching coordinated control system of the invention, including tuning controller,
The tuning controller is given using event driven cutting load described above with corrective route switching control method for coordinating
It realizes.
Compared with prior art, the beneficial effects of the present invention are:
(1) event driven cutting load of the invention and corrective route switching control method for coordinating, pass through event-driven
The cooperation of type cutting load and corrective route switching maintains the safety of electric system steady with less and controllable load loss
It is fixed, reduce cutting load amount, reduce scope of power outage, reduces control cost, improve economy.
(2) present invention is in a unified Optimization Framework, to event driven cutting load and corrective route switching
Coordinated control optimizes, and has taken into account frequency transient security, voltage transient safety and electric equipment overload, has avoided control
Negative effect is more advantageous to the safety and stability of protection electric system.
(3) control program of the invention starts immediately when failure starts, and can contain frequency, voltage transient unstability in time
Trend, avoid the generation of a wide range of power flow transfer, reduce influence of the failure to Operation of Electric Systems, enhancing systemic defence therefore
Hinder the ability of disturbance.
Detailed description of the invention
The accompanying drawings constituting a part of this application is used to provide further understanding of the present application, and the application's shows
Meaning property embodiment and its explanation are not constituted an undue limitation on the present application for explaining the application.
Fig. 1 is event driven cutting load and corrective route switching control method for coordinating flow diagram of the invention；
Fig. 2 is double-layer lap of the invention for method flow schematic diagram；
Fig. 3 is linearization technique flow diagram of the invention；
Fig. 4 is gradually summation flow diagram of the invention；
Fig. 5 is integration iterative method flow diagram of the invention；
Fig. 6 is that example one improves 39 node system structure chart of IEEE；
Fig. 7 is that certain Provincial Electric Power System of example two simplifies geographical wiring diagram.
Specific embodiment
It is noted that following detailed description is all illustrative, it is intended to provide further instruction to the application.Unless another
It indicates, all technical and scientific terms used herein has usual with the application person of an ordinary skill in the technical field
The identical meanings of understanding.
It should be noted that term used herein above is merely to describe specific embodiment, and be not intended to restricted root
According to the illustrative embodiments of the application.As used herein, unless the context clearly indicates otherwise, otherwise singular
Also it is intended to include plural form, additionally, it should be understood that, when in the present specification using term "comprising" and/or " packet
Include " when, indicate existing characteristics, step, operation, device, component and/or their combination.
Fig. 1 is event driven cutting load and corrective route switching control method for coordinating flow diagram of the invention.
As shown in Figure 1, event driven cutting load of the invention and corrective route switching control method for coordinating, comprising:
Step (1): acquisition actual electric network data build the simulation model of studied power grid；Cutting load, route switching are set
Control parameter, simulation parameter, monitoring data, transient security constrained parameters and overload constrained parameters.
Specifically, the electric network data in step (1) includes grid structure, node parameter, line parameter circuit value, power parameter and bears
Lotus parameter.
Wherein, the simulation model of studied power grid can have power system transient simulation in BPA, PSASP and PSS/E etc.
It is built in any software of function.
Setting cutting load control parameter comprises determining that cutting load interstitial content, position and each node maximum can cutting loads
Amount forms cutting load dominant vector and maximum cutting load constrained vector, building cutting load amount constraint.
The number N (N is the positive integer more than or equal to 1) of cutting load node and position, are selected to use from power grid by user
The interstitial content and position for carrying out cutting load determine；Each cutting load node maximum can cutting load amount it is fixed according to the actual situation by user
Justice；Cutting load interstitial content N determines cutting load dominant vector P and maximum cutting load constrained vector PmaxDimension, each cutting load section
The position of point corresponds to it in P and PmaxIn number；Cutting load dominant vector P and maximum cutting load constrained vector PmaxIn variable
It is real number, and
P=[p1,p2,...,pi,...,pN]
Pmax=[p1,max,p2,max,...,pi,max,...,pN,max]
Wherein, p1,p2,...,pi,...,pNThe cutting load amount of respectively first cutting load node, second cutting load
The cutting load amount ... of node, the cutting load amount ... of i-th of cutting load node, the cutting load amount of n-th cutting load node；
p1,max,p2,max,...,pi,max,...,pN,maxThe maximum cutting load amount of respectively first cutting load node, second
The maximum cutting load amount ... of a cutting load node, the maximum cutting load amount ... of i-th of cutting load node, n-th cutting load section
The maximum cutting load amount of point.
In step (1), the constraint of cutting load amount are as follows:
0≤P≤Pmax。
Wherein, comprise determining that can switching number of lines, position and maximum switching route for setting route switching control parameter
Number forms route switching dominant vector, building route switching constraint.
Security of system may further be deteriorated due to cut-offfing branch, and calculation amount and computation complexity can be significantly greatly increased,
For route switching, the investment of specific extension wire is only considered, do not consider cut-offfing for working line.
Can switching number of lines M (M be positive integer) more than or equal to 1 and position, by user selected from power grid can
The number and location of the extension wire to put into operation determine；Maximum switching number of lines lmaxIt is defined according to the actual situation by user；
Can switching number of lines M determine the dimension of route switching dominant vector L, respectively can the position of switching route correspond to its volume in L
Number；Variable in route switching dominant vector L is 0-1 variable, and
L=[l1,l2,...,lj,...,lM]
Wherein, l1,l2,...,lj,...,lMRespectively first switching route cuts control amount, Article 2 switching route
Cut control amount ..., j-th strip switching route cuts control amount ..., and the M articles switching route cuts control amount.
lj=0 and lj=1 respectively indicates route j to cut-off state and investment state；If certain extension wire puts into operation,
Then its state is 1 by 0 variation.
In step (1), the constraint of route switching are as follows:
Step (2): selection forecast accident executes transient emulation and obtains monitoring data, calculate when not taking control measure
Corresponding security margin index and compared with default safety margin critical value judges system whether transient security has no branch to pass by
It carries, if it is not, thening follow the steps (3)；If so, terminating this calculating, determines and do not need to take control for current forecast accident
Measure further determines whether that all forecast accidents all handle completion, if so, thening follow the steps (4)；Otherwise, return step
(2)。
In specific implementation, simulation parameter includes power grid transient emulation total duration T and step delta t；Wherein, T and Δ t are
Parameter preset；
Monitor that data include the frequency of each bus and the apparent energy of voltage trace and each branch during transient emulation
Numerical value；Transient security constrained parameters include frequency two-element list and voltage two-element list, and overload constrained parameters are the volume of each branch
Determine apparent energy.
Specifically, frequency binary sheet form is [fcr, tcr], as frequency transient security constrained parameters, it is used to determination frequency
Transient security, wherein fcrFor the threshold value of bus frequency shift (FS), tcrIt is more than that the maximum of corresponding threshold value can connect to be frequency shift (FS)
By the duration；If the offset of frequency locus is more than f during transient emulationcrTime be no more than tcr, then assert should
Frequency locus safety, conversely, then assert that the frequency locus is dangerous；Frequency two-element list has multiple, only all frequency locus quilts
When all frequency two-element lists are evaluated as safe, corresponding analogue system is just judged as being transient frequency safety.
Specifically, voltage binary sheet form is [vcr, tcr], as voltage transient security constraint parameter, for judging voltage
Transient security, wherein vcrFor the threshold value of busbar voltage offset, tcrIt is more than that the maximum of corresponding threshold value can connect to be variation
By the duration；If the offset of voltage trace is more than v during transient emulationcrTime be no more than tcr, then assert should
Voltage trace safety, conversely, then assert that the voltage trace is dangerous；Voltage two-element list has multiple, only all voltage trace quilts
When all voltage two-element lists are evaluated as safe, corresponding analogue system is just judged as being transient voltage security.
Overload constrained parameters only consider the overload situations of 220kV and above branch in power transmission network, to matching
The overload situations of each branch are not considered temporarily in power grid.
The apparent energy numerical value PC of each branchi,nAs overload constrained parameters, for judging whether branch overloads；If
The apparent energy PC of i-th branch roadiApparent energy numerical value PC corresponding no more than itsi,n, then assert that the branch did not occurred
It carries, conversely, then assert that the branch overloads；When only all branches all do not overload, corresponding analogue system is just judged to
Break is not overload.
Moreover, in this step, security margin index includes transient safe and stable margin index vector ηTAnd branch passes by
Carry margin index ηS, the transient security constrained parameters and mistake of the middle setting of the monitoring data and step (1) obtained by transient emulation
Constrained parameters are carried to calculate；Safety margin critical value includes transient security nargin critical value vector εTWith overload nargin critical value
εS。
Transient safe and stable margin index vector ηTIncluding transient frequency security margin index ηfAnd transient voltage security
Stability margin index ηv, i.e.,
Transient frequency security margin index ηfThe cumulative effect for considering frequency locus offset, can preferably reflect
System transient modelling frequency security situation, and have the advantages that it is linear, smooth and dull, and
ηf=min (ηf,i,j), i=0 ..., Nf, j=0 ..., Mf
Wherein, NfNumber for the bus frequency locus monitored during transient state numerical simulation, MfFor the number of frequency two-element list
Mesh, t are initial time, fNFor system nominal frequency, ηf,i,jFor frequency locus fiIn frequency two-element list [fcr,i, tcr,j] under it is temporary
State frequency security stability margin index；T is power grid transient emulation total duration, is parameter preset.
Transient voltage safe and stable margin index ηvThe cumulative effect for considering voltage trace offset, can preferably reflect
System transient modelling voltage security situation, and have the advantages that it is linear, smooth and dull, and
ηv=min (ηv,i,j), i=0 ..., Nv, j=0 ..., Mv
Wherein, NvNumber for the busbar voltage track monitored during transient state numerical simulation, MvFor the number of voltage two-element list
Mesh, t are initial time, vNFor bus voltage rating, ηv,i,jFor voltage trace viIn voltage two-element list [vcr,i, tcr,j] under it is temporary
State voltage security stability margin index；T is power grid transient emulation total duration, is parameter preset.
Transient security nargin critical value vector εTIncluding transient frequency security margin critical value εfAnd transient voltage security
Stability margin critical value εv, i.e.,
If ηf>εf, then assert that system is that transient frequency is stable；Conversely, then assert system not is that transient frequency is stable.
If ηv>εv, then assert that system is Transient Voltage Stability；Conversely, then assert that system is not Transient Voltage Stability.Only when being
Uniting while meeting transient frequency to stablize with Transient Voltage Stability ability identification system is transient stability, even ηT>εT, then assert and be
System is transient stability；Conversely, then assert that system is not transient stability.
Wherein: branch overloads margin index ηSFor
Wherein, NPCNumber for the branch apparent energy monitored during transient state numerical simulation.
If ηs>εs, then assert that branch overload does not occur for system；Conversely, then assert that branch overload occurs for system.
In this step, if ηT>εTAnd ηs>εs, then assert and be directed to the forecast accident, system meets transient security and has no branch
Load is passed by, does not need to take control measure；Conversely, needing to be implemented step (3), then to take control measure.
Step (3): Optimized model is set up by the forecast accident selected and the parameter of setting, is iterated solution, calculating is cut
Load control system vector sum route switching dominant vector, in conjunction with the position of each cutting load node and switching route in power grid, system
Surely it is directed to the control program of current forecast accident；Further determine whether that all forecast accidents all handle completion, if so, executing
Step (4)；Otherwise, return step (2).
In specific implementation, event driven cutting load and corrective route switching are carried out in a unified frame
Optimization, the Optimized model of establishment is,
Wherein, F is cutting load total amount required by Optimized model.
Specifically, in this step, to solve above-mentioned Optimized model, two submodels are broken down into, and using double-deck
Alternative manner is divided into internal layer iteration and external iteration, and internal layer is iterated event driven cutting load dominant vector, outer layer pair
Correcting property route switching dominant vector is iterated.
For submodel first is that optimizing to event driven cutting load, consideration guarantees power grid transient frequency peace by cutting load
Complete and transient voltage security, and,
Submodel two optimizes event driven cutting load and corrective route switching, and consideration is limited by route switching
Overload degree processed thoroughly eliminates overload by cutting load, and,
Wherein, FTAnd FSIt respectively represents to ensure transient security and eliminating the required cutting load amount of overload, then former Optimized model
Required cutting load amount F is,
F=FT+FS。
As shown in Fig. 2, in step (3), the specific steps of the double-deck alternative manner are as follows:
Step (3-1): parameter initialization, and former Optimized model is decomposed into two submodels；
Step (3-2): a kind of corrective route switching scheme is selected；
Step (3-3): under selected corrective route switching scheme, being iterated solution to two submodels, calculates
Current optimal cutting load dominant vector and corresponding F；
Step (3-4): if having completed the search to all corrective route switchings, corresponding minimum F value is chosen
Scheme is as final scheme；Otherwise, it returns and executes step (3-2).
In step (3-3), the cutting load control section in submodel one and submodel two is iteratively solved, F is soughtTAnd FS
When, it mainly uses linearization technique, gradually summation and integrates iterative method.
The linearization technique is used to solve submodel one, seeks FT, by calculating transient security margin index ηTTrack
Sensitivity, linearisation transient security constraint.
As shown in figure 3, the linearization technique specific steps are as follows:
(3-3-1-1) parameter initialization, settingK=1,；
(3-3-1-2) calculates trace sensitivity matrix, and the transient security in submodel one is constrained linearisation；
In step (3-3-1-2), track needed for calculating is obtained by system transient modelling numerical simulation, corresponding trace sensitivity
It is calculated and is obtained by numerical perturbation method:
Wherein, ηf(P) and ηv(P) margin index when be cutting load amount being P, ηv(P, τi) it be cutting load amount is P, and cutting
Increase perturbation τ at load point iiWhen margin index.
Using trace sensitivity, transient security is constrained and is linearized, i.e., by the transient security constraint condition in submodel one
ηT>εTConversion are as follows:
ηT(PT)+AΔPT≥εT
Submodel one is converted linear optimization model by (3-3-1-3), is solved using linear programming methodIt calculates
Corresponding cutting load dominant vector
In step (3-3-1-3), due to objective function and cutting load the amount constraint of submodel one be it is linear,
After linearisation transient security constraint, linear optimization model can be converted by submodel one, the method that iterative solution can be used.
Model equation when+1 iteration of kth are as follows:
Wherein, unknown variable isThe cutting load dominant vector of+1 iteration of kthIt calculates according to the following formula:
(3-3-1-4) judges whether to meet the condition of convergence, if it is satisfied, then by current cutting load dominant vectorAs
Optimal cutting load dominant vector PT, export result and terminate step calculating；Otherwise, k=k+1, return step (3-3-1-2) are enabled.
In step (3-3-1-4), there are two the conditions of convergence, and need to meet simultaneously, specifically:
(a)ηfOr ηvReach critical value, i.e.,
εf≤ηf≤εf+ΔεfOr εv≤ηv≤εv+Δεv
Wherein, Δ εfWith Δ εvRespectively εfAnd εvAllowable error, and Δ εfWith Δ εvIt should be set as the positive number of very little, with true
The iterative process is protected to converge in a lesser section.
(b) difference of the cutting load amount of iteration is less than pre-set threshold epsilon twice before and afterp, i.e.,
Wherein, εpIt should be set as the positive number of very little, to ensure that the iterative process converges in a lesser section.
The gradually summation is used to solve cutting load control section in submodel two, seeks FS。
As shown in figure 4, the gradually summation, the specific steps are that:
P is arranged in (3-3-2-1) parameter initializationS=0 and load variations step delta P；
In step (3-3-2-1), the setting of load variations step delta P, using variable step technology: the initial stage is used
Biggish Δ P accelerates solution procedure can reduce emulation step number；The improvement of adjoint system overload nargin, is used instead lesser
Δ P reduced cutting load amount to prevent from cutting.
(3-3-2-2) calculates sensitivity of the cutting load amount to overload nargin at each bus, in cutting with peak response
Increase cutting load amount Δ p at load point, and updates PS；
In this step, the sensitivity a of nargin is overloadeds,iCalculating it is as follows:
Wherein, ηS(P) margin index when be cutting load amount being P, ηS(P, τi) it be cutting load amount is P, and in cutting load point i
Place increases perturbation τiWhen margin index.
(3-3-2-3) cuts off current PSAfterwards, it if system is still overloaded, returns and executes step (3-3-2-2)；Otherwise, PSAs
Cutting load vector needed for overloading is eliminated, result is exported and terminates step calculating.
Wherein, iterative method is integrated, above-mentioned linearization technique and gradually summation are integrated, to submodel one and submodel two
Cutting load control section is successively iterated, and seeks F=FT+FSOptimal, rather than FTAnd FSIt is respective optimal.
As shown in figure 5, described integrate iterative method, the specific steps are that:
(3-3-3-1) carries out parameter initialization under selected corrective route switching scheme, and parameter is arranged R=1；
The lower limit of cutting load dominant vector in submodel one is revised as by (3-3-3-2)Pass through above-mentioned linearization technique
Submodel one is solved, the optimal cutting load dominant vector P of current submodel one is calculatedT, and enable
(3-3-3-3) is when excision loadAfterwards, if analogue system is able to maintain transient security and has no branch overload, it is arrangedAnd execute step (3-3-3-4)；Otherwise, submodel two, meter are solved by above-mentioned gradually summation
Calculate the optimal cutting load dominant vector P of current submodel twoS, and enable R=r+1 is returned and is executed step (3-3-3-2)；
(3-3-3-4) selectionAs the final cutting load dominant vector under current corrective route switching scheme, and count
Calculate corresponding F.
In step (3), the control program of formulation includes the load that each cutting load node is cut and route switching plan
Slightly.
Step (4): in grid collapses, corresponding control program is executed, and Starting mode is to detect failure
Then start immediately；At regular intervals, it is executed since step (1) again, updates the corresponding control program of each forecast accident.
In this step, start immediately when control program Starting mode is the failure generation of monitoring, after scheme starting, event
Driving cutting load control cuts off corresponding load immediately, and corrective route switching will terminate substantially in system dynamic course, wave
Dynamic relatively steady rear implementation, in order to avoid bigger disturbance is caused to system due to the variation of network topology.
Using off-line calculation, the control strategy of On-line matching at regular intervals recalculates control strategy.
As shown in fig. 6, example one is to improve 39 node system of IEEE, wherein L1, L2 and L3 are extension wire, specific
Situation is as shown in table 1.
1 example of table, one extension wire situation
It is out of service with generator failure at 30 and 39 two buses, cause the loss of 1250MW generated energy as failure field
Scape.When not considering the cooperation of cutting load and route switching, total cutting load amount is 861.1MW；When applying the present process invention, most
Excellent cutting load amount is 822MW, specific as shown in table 2.Compare two methods, sufficiently demonstrate the method for the present invention cut in reduction it is negative
Superiority in terms of lotus, reduction control cost.
The optimal cutting load scheme of 2 example of table, one the method for the present invention
As shown in fig. 7, example two is certain provincial power network, wherein L1, L2, L3 and L4 are extension wire, and concrete condition is such as
Shown in table 3.
3 example of table, two extension wire situation
With 6 generator trippings simultaneously at A1 bus in system, 1200MW active power shortage is caused, as failure field
Scape.When not considering the cooperation of cutting load and route switching, total cutting load amount is 956.7MW；When applying the present process invention, most
Excellent cutting load amount is 510MW, specific as shown in table 4.Compare two methods, sufficiently demonstrate the method for the present invention cut in reduction it is negative
Superiority in terms of lotus, reduction control cost.
The optimal cutting load scheme of 4 example of table, two the method for the present invention
In addition, the present invention also provides a kind of event driven cutting loads and corrective route switching coordinated control system.
Event driven cutting load and corrective route switching coordinated control system of the invention, including tuning controller,
The tuning controller is given using event driven cutting load as shown in Figure 1 with corrective route switching control method for coordinating
It realizes.
Above-mentioned, although the foregoing specific embodiments of the present invention is described with reference to the accompanying drawings, not protects model to the present invention
The limitation enclosed, those skilled in the art should understand that, based on the technical solutions of the present invention, those skilled in the art are not
Need to make the creative labor the various modifications or changes that can be made still within protection scope of the present invention.
Claims (10)
1. a kind of event driven cutting load and corrective route switching control method for coordinating characterized by comprising
Step (1): acquisition actual electric network data build the simulation model of studied power grid；Cutting load, the control of route switching are set
Parameter, simulation parameter, monitoring data, transient security constrained parameters and overload constrained parameters；
Step (2): selection forecast accident executes transient emulation and obtains monitoring data, calculate corresponding when not taking control measure
Security margin index and compared with default safety margin critical value, judge system whether transient security have no branch overload, if
It is no, then follow the steps (3)；If so, terminating this calculating, determines and does not need to take control measure for current forecast accident,
Further determine whether that all forecast accidents all handle completion, if so, thening follow the steps (4)；Otherwise, return step (2)；
Step (3): Optimized model is set up by the forecast accident selected and the parameter of setting, is iterated solution, calculates cutting load
Dominant vector and route switching dominant vector formulate needle in conjunction with the position of each cutting load node and switching route in power grid
To the control program of current forecast accident；Further determine whether that all forecast accidents all handle completion, if so, thening follow the steps
(4)；Otherwise, return step (2)；
Step (4): in grid collapses, corresponding control program is executed, and Starting mode is to detect that failure generation is then stood
Start；At regular intervals, it is executed since step (1) again, updates the corresponding control program of each forecast accident；
The Optimized model are as follows:
Wherein, F is cutting load total amount required by Optimized model；N is the number of cutting load node；PiIt is controlled for the cutting load of node i
Amount；Pi,maxFor the maximum cutting load amount of node i；ηTFor transient safe and stable margin index vector；εTFace for transient security nargin
Dividing value vector；ηSMargin index is overloaded for branch；εSNargin critical value is overloaded for branch；M be can switching route number；ljFor
Switching route j's cuts control amount, lj=0 and lj=1 respectively indicates route j to cut-off state and investment state；lmaxFor maximum switching
Number of lines.
2. event driven cutting load as described in claim 1 and corrective route switching control method for coordinating, feature exist
In by the forecast accident selected and the parameter of setting establishment Optimized model in the step (3) are as follows:
The minimum value of event driven cutting load is the sum of the excision load of cutting load node of all determinations；
Wherein, parameter sets up the constraint condition of Optimized model are as follows: system meets transient security respectively and has no branch overload, route throwing
Cut constraint and the constraint of cutting load amount.
3. event driven cutting load as claimed in claim 2 and corrective route switching control method for coordinating, feature exist
In setting up Optimized model to solve parameter, be broken down into two submodels, and use double-layer lap generation in the step (3)
Method is divided into internal layer iteration and external iteration, and internal layer is iterated event driven cutting load dominant vector, and outer layer is to correction
Property route switching dominant vector is iterated.
4. event driven cutting load as claimed in claim 3 and corrective route switching control method for coordinating, feature exist
In in the step (3), to solve parameter establishment Optimized model, being broken down into two submodels is respectively submodel one
With submodel two；Wherein, submodel is first is that consider through cutting load guarantee power grid transient frequency safety and transient voltage security, and
Event driven cutting load is optimized；
Submodel is second is that consider thoroughly to eliminate overload by cutting load by route switching limiting overload degree, and drive event
Ejector half cutting load and corrective route switching optimize.
5. event driven cutting load as claimed in claim 4 and corrective route switching control method for coordinating, feature exist
In in the step (3), using the detailed process of double-deck alternative manner solution parameter establishment Optimized model are as follows:
Step (3-1): parameter initialization, and former Optimized model is decomposed into two submodels；
Step (3-2): a kind of corrective route switching scheme is selected；
Step (3-3): under selected corrective route switching scheme, solution is iterated to two submodels, is calculated current
Optimal cutting load dominant vector and institute's cutting load total amount；
Step (3-4): if having completed the search to all corrective route switchings, corresponding minimum tangential load total amount is chosen
Scheme as final scheme；Otherwise, return step (3-2).
6. event driven cutting load as claimed in claim 5 and corrective route switching control method for coordinating, feature exist
In in the step (3-3), submodel one is solved using linearization technique, solves submodel using gradually summation
Two.
7. event driven cutting load as claimed in claim 6 and corrective route switching control method for coordinating, feature exist
In being successively iterated using iterative method is integrated to the cutting load control section of submodel one and submodel two, seek event
Driving cutting load it is optimal.
8. event driven cutting load as described in claim 1 and corrective route switching control method for coordinating, feature exist
In in the step (3), the control program of formulation includes the load and route switching strategy that each cutting load node is cut.
9. event driven cutting load as described in claim 1 and corrective route switching control method for coordinating, feature exist
In,
In the step (4), control program Starting mode is that the failure of monitoring starts immediately when occurring, after scheme starting, event
Driving cutting load control cuts off corresponding load immediately, and corrective route switching will terminate substantially in system dynamic course, wave
Dynamic relatively steady rear implementation, in order to avoid bigger disturbance is caused to system due to the variation of network topology；
Or/and
In the step (4), using off-line calculation, the control strategy of On-line matching at regular intervals carries out control strategy
It recalculates.
10. a kind of event driven cutting load and corrective route switching coordinated control system, which is characterized in that controlled including coordinating
Device processed, the tuning controller is using event driven cutting load as claimed in any one of claims 1-9 wherein and corrective line
Road switching control method for coordinating is achieved.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710491138.XA CN107171334B (en) | 2017-06-23 | 2017-06-23 | Event driven cutting load and corrective route switching control method for coordinating and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710491138.XA CN107171334B (en) | 2017-06-23 | 2017-06-23 | Event driven cutting load and corrective route switching control method for coordinating and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107171334A CN107171334A (en) | 2017-09-15 |
| CN107171334B true CN107171334B (en) | 2019-06-28 |
Family
ID=59820644
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710491138.XA Active CN107171334B (en) | 2017-06-23 | 2017-06-23 | Event driven cutting load and corrective route switching control method for coordinating and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107171334B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111740416B (en) * | 2020-07-08 | 2022-06-10 | 国网河南省电力公司经济技术研究院 | Target cascade analysis-based receiving-end power grid direct-current maximum feed-in quantity evaluation method |
| CN114389265B (en) * | 2022-03-23 | 2022-08-16 | 山东大学 | Power grid emergency load optimization method and system based on proliferation and reduced evolution |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106570791A (en) * | 2016-11-13 | 2017-04-19 | 国网江苏省电力公司连云港供电公司 | Dynamic adjustment method for structure of regional power grid in allusion to main transformer overload |
-
2017
- 2017-06-23 CN CN201710491138.XA patent/CN107171334B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106570791A (en) * | 2016-11-13 | 2017-04-19 | 国网江苏省电力公司连云港供电公司 | Dynamic adjustment method for structure of regional power grid in allusion to main transformer overload |
Non-Patent Citations (1)
| Title |
|---|
| Optimization of the Event-Driven Emergency Load-Shedding Considering Transient Security and Stability Constraints;Xinxu等;《IEEE》;20161031;第2581-2592页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107171334A (en) | 2017-09-15 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108418255B (en) | Planning method and system suitable for extra-high voltage direct current outgoing power grid containing high-permeability new energy | |
| Sansawatt et al. | Smart decentralized control of DG for voltage and thermal constraint management | |
| CN107834567B (en) | Reactive voltage coordination control method for extra-high voltage direct current converter station and near-field power plant | |
| CN105375533B (en) | A kind of independent micro-capacitance sensor operational mode Forecasting Methodology containing wind, light regenerative resource | |
| CN105914753A (en) | Active distribution network (ADN) voltage real-time fuzzy control method | |
| CN106099951A (en) | Load cutting method and device for microgrid type combined generating system | |
| CN105262135A (en) | Wind, light and diesel storage microgrid system containing composite energy storage and coordinated control method during grid connection | |
| CN110896231B (en) | Distributed photovoltaic capacity receiving calculation method and system for power distribution network in poverty alleviation area | |
| CN104934971A (en) | Dynamic section control method based on power flow transfer ratio | |
| CN110994626B (en) | 500-220kV regional power grid automatic voltage control method based on voltage trend prediction | |
| CN107171334B (en) | Event driven cutting load and corrective route switching control method for coordinating and system | |
| Roy et al. | Virtual inertia support in power systems for high penetration of renewables—overview of categorization, comparison, and evaluation of control techniques | |
| CN104638638B (en) | A kind of safety on line stable tendency analysis method for bulk power grid | |
| CN115733188A (en) | Novel distributed new energy bearing capacity assessment method for power system considering multi-source complementary characteristics | |
| CN114024331B (en) | Planning operation measure optimization method for improving direct current bearing capacity | |
| Lei et al. | Look-ahead rolling economic dispatch approach for wind-thermal-bundled power system considering dynamic ramping and flexible load transfer strategy | |
| CN107769218A (en) | A kind of comprehensive coordination control strategy of electric distribution network reactive-voltage | |
| CN104767207A (en) | Voltage stabilization prevention and control method based on reactive storage sensitivity of generators | |
| CN109102196B (en) | Power distribution network diversified load digestion capacity assessment method and system | |
| Aguilera-Gonzalez et al. | MPC Energy Management System for A Grid-Connected Renewable Energy/Battery Hybrid Power Plant | |
| CN110718933A (en) | Multilevel coordinated wind storage isolated network system power balance control strategy | |
| CN106447525B (en) | A kind of power grid dynamic partition policy optimization searching method | |
| Jia et al. | An optimal control strategy for lv distribution network with pv and energy storage units | |
| Vinot et al. | Congestion Avoidance in Low-Voltage Networks by using the Advanced Metering Infrastructure | |
| CN102957147B (en) | Multi-load characteristic-matched under-frequency load shedding optimal control system and multi-load characteristic-matched under-frequency load shedding optimal control method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |