CN106532717A - Comprehensive sensitivity analysis-based circuit overload load-shedding coordination optimization method - Google Patents

Comprehensive sensitivity analysis-based circuit overload load-shedding coordination optimization method Download PDF

Info

Publication number
CN106532717A
CN106532717A CN201610978458.3A CN201610978458A CN106532717A CN 106532717 A CN106532717 A CN 106532717A CN 201610978458 A CN201610978458 A CN 201610978458A CN 106532717 A CN106532717 A CN 106532717A
Authority
CN
China
Prior art keywords
node
cutting load
load
circuit overload
shedding
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201610978458.3A
Other languages
Chinese (zh)
Other versions
CN106532717B (en
Inventor
张丹丹
杨冬
蒋哲
武乃虎
李文博
邢鲁华
侯燕文
麻常辉
张荻
张一荻
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Shandong Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
Electric Power Research Institute of State Grid Shandong Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, Electric Power Research Institute of State Grid Shandong Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN201610978458.3A priority Critical patent/CN106532717B/en
Publication of CN106532717A publication Critical patent/CN106532717A/en
Application granted granted Critical
Publication of CN106532717B publication Critical patent/CN106532717B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • H02J3/12—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
    • H02J3/14—Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by switching loads on to, or off from, network, e.g. progressively balanced loading
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas or water supply
    • 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
    • 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
    • Y04S10/00—Systems supporting electrical power generation, transmission or distribution
    • Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
    • 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

Abstract

The invention discloses a comprehensive sensitivity analysis-based circuit overload load-shedding coordination optimization method, which comprises the steps of node filtering, namely building a node comprehensive sensitivity index of considering the branch overload degree according to the sensitively relationship between branch current and node injection current and carrying out alternative load-shedding node filtering through calculation of the node comprehensive sensitivity index; scheme searching, coordinating two objects of the load-shedding economic cost and the node load-shedding balance and building a circuit overload load-shedding coordination optimization control model by taking the minimum weighted value of the economic cost and a node load-shedding balance index as an optimization object; and setting of the weighted value of the load-shedding economic cost and the balance index, namely calculating an optimal circuit overload load-shedding coordination optimization control scheme of the circuit overload load-shedding coordination optimization control model under different weights and finally obtaining a load-shedding node. According to the method disclosed by the invention, the scheme searching efficiency is greatly improved and a dimension disaster in a search optimization process is avoided.

Description

Circuit overload connection cutting load coordination optimizing method based on comprehensive sensitivity analysis
Technical field
The present invention relates to electric power network technique field, and in particular to the circuit overload connection cutting load association based on comprehensive sensitivity analysis Adjust optimization method.
Background technology
The continuous development of trans-regional interconnected power grid, solves energy resources inverse with load center to a certain extent To distribution problem, economical and efficient, the clean environment firendly distributed rationally with operation of power networks across big area resource is realized, but while to electricity The safe and stable operation of net brings unprecedented challenge.
In recent years, many large area blackouts show both at home and abroad：In long-distance and large-capacity power transmission system, circuit overload can Can cause a series of cascading trip accidents, in turn result in large-area power-cuts.Therefore, how to formulate rational circuit overload connection to cut Duty control method, circuit overload is quickly eliminated, be to prevent the chain power outage of large area, ensure power grid security, stable, economic The key of operation.
The circuit overload connection load shedding control method for adopting in real system at present is mainly minute wheel classification cutting load method.Should Method refers to that different according to circuit overload degree, minute wheel level cuts off according to each load excision round and prioritization is preset Load.As the method is when each load excision round and prioritization is preset, relies on related personnel's experience more, lack The theoretical calculation of system, therefore load is often resulted in practical operation owe to cut or cross and cut, and in the method load excision priority Sequence is transshipped just for uniline, when in system, a plurality of circuit transships simultaneously, cannot be to cutting load measure using the method Carry out global optimization.
In a word, solved for uniline overload more than above-mentioned circuit overload connection cutting load method, lacked to a plurality of line Coordination optimization control when road is transshipped simultaneously；And in Practical Project, electric power enterprise, should when overload joining-cutting load measure is taken Pursue financial cost to minimize, consider the harmony of cutting load again, to reduce accident impact, the grading of reduction accident, and show There is method not take in cutting load harmony.
In consideration of it, the application proposes a kind of circuit overload connection cutting load coordination optimization side based on comprehensive sensitivity analysis Method.
The content of the invention
To solve the deficiency that prior art is present, the invention discloses the circuit overload connection based on comprehensive sensitivity analysis is cut Coordinated optimization method, the purpose of the present invention is：For one or more of circuits after failure while overload situations, establish Overall sensitivity index of the node cutting load measure to overload circuit, by the alternative cutting load node of index calculating sifting, it is to avoid Calculating is absorbed in dimension calamity；Method constructs circuit overload connection cutting load coordination optimization Controlling model, is cut with node with financial cost The minimum optimization aim of weighted value of load balancing index, by trend active, reactive balance, node voltage, current limiter ginseng Number is not out-of-limit to wait restriction to include constraints.The method achieve financial cost when overload joining-cutting load is controlled balanced with cutting load Property coordination optimization control, and by arranging different weights, meeting different operations of power networks requires, it is ensured that circuit overload quickly disappears Remove, meet Practical Project demand.
For achieving the above object, concrete scheme of the invention is as follows：
Circuit overload connection cutting load coordination optimizing method based on comprehensive sensitivity analysis, comprises the following steps：
Node is screened：Network node admittance matrix and incidence matrix are formed, it is fixed with incidence matrix according to nodal voltage equation Justice, draws the sensitivity relation between branch current and node Injection Current；
Meter is set up according to the sensitivity relation between branch current and node Injection Current and branch road transships the node of degree Overall sensitivity index, and alternative cutting load node screening is carried out by the calculating of node overall sensitivity index；
Scheme is searched for：Coordinate cutting load financial cost and node cutting load two targets of harmony, with financial cost and node The minimum optimization aim of weighted value of cutting load harmony index, sets up circuit overload connection cutting load coordination optimization Controlling model；
It is different according to requiring to cutting load financial cost and harmony, cutting load financial cost is set with harmonious index Weighted value, calculates the circuit overload connection cutting load and coordinates and optimizes Controlling model under different weights, and optimal route overload connection is cut Coordinated optimization control scheme, finally draws cutting load node.
Further, defined according to nodal voltage equation and incidence matrix, draw branch current and node Injection Current it Between sensitivity relation be：
Wherein,For branch current column vector；Y is Nl×NlRank branch admittance matrix,For the inverse of bus admittance matrix Battle array,For node Injection Current column vector, node is flowed into for just, it is negative to flow out, ATFor incidence matrix transposition.
Further, take cutting load measure ability to be eliminated to the excessively stream of all overload branch roads according to node k define node k Sensitivity index, it is considered to which branch road overload degree is different, defines node overall sensitivity index：
Wherein, δjFor branch road j weighted values, for representing the overload degree of branch road j；Ij、The actual electricity of respectively branch road j Flow valuve and the permission current value upper limit, NsTo transship branch road total number, YjFor branch admittance battle array jth row element, the i.e. self-conductance of branch road j Receive and transadmittance, ATFor incidence matrix transposition,For the inverse matrix of bus admittance matrix.
Further, node overall sensitivity index is normalized calculating, for intuitively expressing all overload branch roads pair The overall sensitivity of a certain node cutting load measure.
Further, when alternative cutting load node screening is carried out, take node overall sensitivity value λkMore than setting value Node includes alternative cutting load set of node.
Further, cutting load measure financial cost is：
Wherein, NkFor cutting load node total number；ΔPki、ciThe i-th type load cutting load amount of respectively node k, excision i-th Type load unit-economy cost coefficient.
Further, to also including the step of being normalized calculating after the calculating of cutting load measure financial cost.
Further, the cutting load harmony index definition is node cutting load rate standard deviation：
Wherein, NkFor cutting load node total number,For k node cutting load rates；For all cutting load node cutting load rates Mean value, PkiFor the i-th type load total amount of node k, Δ PkiFor the i-th type load cutting load amount of node k.
Further, the step of node cutting load rate harmony index being normalized.
Further, circuit overload connection cutting load coordination optimization Controlling model is：
Min f=ω1f1+ω2f2
Wherein, ω1、ω2Respectively cutting load financial cost, the weighted value of node cutting load harmony index, carry out normalizing Change cutting load measure financial cost f after calculating1, it is normalized node cutting load rate harmony index f after calculating2。
Further, the constraints of circuit overload connection cutting load coordination optimization Controlling model includes that trend is active, idle Balance, node voltage and line electricity flow valuve be not out-of-limit.
Beneficial effects of the present invention：
(1) circuit overload connection cutting load control is decomposed into node screening and searches for two steps with scheme by the present invention, by meter Operator node overall sensitivity carries out alternative cutting load node screening, effectively reduces cutting load node searching scope, greatly improves Scheme search efficiency, it is to avoid be absorbed in dimension calamity during Optimum search.
(2) related personnel's experience being relied on circuit overload connection load shedding control method in the past more, often result in load and owe to cut or mistake Cut, when a plurality of circuit of generation transships simultaneously, it is difficult to coordinate optimal control, for this purpose, the present invention is established with financial cost With the circuit overload connection cutting load coordination optimization control mould of the minimum optimization aim of weighted value of node cutting load harmony index Type, by trend active, reactive balance, node voltage, line electricity flow valuve be not out-of-limit etc. limits and includes constraints the model, realizes Coordination optimization control when a plurality of circuit transships simultaneously.
(3), in Practical Project, electric power enterprise should pursue financial cost minimum when overload joining-cutting load is taken, again Consider the harmony of cutting load, reduce accident impact, reduce accident grading and social influence, in the set up model of the present invention, lead to The weighted value of adjust economy cost and both node cutting load harmony indexs is crossed, be can be calculated under different weights, circuit overload Connection cutting load coordination optimization control program, meets the demand of electric power enterprise in Practical Project.
Description of the drawings
Fig. 1 is coordination optimizing method flow chart；
Fig. 2 is somewhere Power grid structure figure.
Specific embodiment：
The present invention is described in detail below in conjunction with the accompanying drawings：
As shown in figure 1, a kind of circuit overload connection cutting load coordination optimizing method based on comprehensive sensitivity analysis, including such as Lower step：
(1) defined with incidence matrix according to nodal voltage equation, derived between branch current and node Injection Current Sensitivity relation；
(2) set up meter and branch road transships the node overall sensitivity index of degree, the index can reflect all overloads Overall sensitivity of the road to a certain node cutting load measure, and alternative cutting load is carried out by calculate node overall sensitivity index Node is screened；
(3) cutting load financial cost and node cutting load two targets of harmony are coordinated, with financial cost and node cutting load The minimum optimization aim of weighted value of harmonious index, sets up circuit overload connection cutting load coordination optimization Controlling model；
(4) cutting load financial cost and harmony are required according to policymaker it is different, setting cutting load financial cost with it is equal The weighted value of weighing apparatus property index, calculates under different weights, optimal route overload joining-cutting load coordination optimization control program.
With regard to the sensitivity relation between branch current and node Injection Current：With the continuous development construction of electrical network, rack Structure is increasingly complicated, and system interior joint is continuously increased, in conducting transmission line overload joining-cutting load scheme optimization solution procedure, if Optimizing is carried out in full range of nodes, may make scheme that dimension calamity is absorbed in searching process, in system some electrically away from Cutting load measure is taken to affect little to eliminating circuit overload from node farther out, therefore, optimization circuit overload connection is cut by the present invention Load scheme is divided into alternative cutting load node screening and scheme searches for two steps, so as to improve Searching efficiency, is prevented effectively from side Dimension calamity is absorbed in case searching process.
By deriving sensitivity relation between computation overload branch current and node Injection Current, pin can be intuitively obtained More sensitive node is eliminated to a certain overload branch road excessively stream, theoretical foundation is provided for the screening of alternative cutting load node.
Between overload branch current and node Injection Current, sensitivity relation derivation is：
Wherein, formula (1) is nodal voltage equation；Formula (2) is defined for incidence matrix；For node Injection Current column vector, Node is flowed into for just, it is negative to flow out；YnFor Nb×NbRank bus admittance matrix, wherein NbFor number of network node；For node voltage Column vector；For branch voltage column vector；ATFor incidence matrix transposition.
Can be derived as follows by formula (1), (2)：
Sensitivity relation between branch current and node Injection Current can be transshipped by formula (3), i.e.,：
Node overall sensitivity index：In Practical Project, failure often results in a plurality of circuit while transshipping, therefore, carrying out When alternative cutting load node is screened, the present invention proposes node overall sensitivity index, for reflecting the measure of node cutting load to institute There is overload overall sensitivity of the branch road to a certain node cutting load measure.
In Practical Project, as circuit overload degree is different, different degrees of power grid accident, social influence, Jing are often resulted in Ji affects to affect with electric power enterprise power supply quality, therefore when a plurality of circuit of generation transships simultaneously, need to preferentially eliminate overload The more serious branch road of degree, it is to avoid even more serious power grid accident occurs.Therefore, the present invention proposes that meter and branch road transship degree Node overall sensitivity index.
Certain node takes cutting load measure to eliminate ability to the excessively stream of all overload branch roads：
Wherein, NsTo transship branch road total number；For the current variation value of branch road j after node k cutting loads；YjFor branch road Admittance battle array jth row element, the i.e. self-admittance of branch road j and transadmittance；For its current variation value after node k cutting loads.
Defining node k sensitivity is：
The present invention proposes that meter and branch road transship the node overall sensitivity index of degree, and computing formula is：
Wherein, δjFor branch road j weighted values, for representing the overload degree of branch road j；Ij、The actual electricity of respectively branch road j Flow valuve and the permission current value upper limit.
Formula (7) interior joint overall sensitivity index is normalized into calculating：
Wherein, λk、λ′kRespectively the overall sensitivity normalized value of node k and actual value；λ′max、λ′minRespectively basis Formula (7) calculates gained node overall sensitivity maxima and minima.
Node overall sensitivity index is normalized into calculating all overload branch roads can be more intuitively expressed to a certain The overall sensitivity of node cutting load measure.When alternative cutting load node screening is carried out, it is proposed that take node overall sensitivity value λkNode more than 0.1 includes alternative cutting load set of node, it is believed that λkNode less than 0.1 does not possess competitiveness.
Circuit overload connection cutting load coordination optimization Controlling model：Circuit overload connection cutting load control belongs to complex nonlinear rule The problem of drawing, in Practical Project, electric power enterprise not only will reduce cutting load as far as possible when circuit overload connection cutting load measure is taken Amount, pursues financial cost and minimizes, and power grid accident grading, social influence and electric power enterprise that cutting load measure causes also is coordinated The aspects such as power supply quality impact, this requires when cutting load control is carried out, will coordinate optimization from the overall situation, should chase after Ask cutting load financial cost to minimize, avoid the excessive unbalanced phenomena of node cutting load as far as possible again, it is to avoid cause some sensitive The high node load excision ratio of degree is too high or even load complete resection phenomenon, causes major accident to grade, so as to look forward to electric power Industry brings harmful effect with society.
Circuit overload connection cutting load coordinates and optimizes Controlling model and is：
Min f=ω1f1+ω2f2 (11)
Vi min≤Vi≤Vi maxI=1,2,3 ..., Nb (15)
ΔQki=σkΔPkiK=1,2,3 ..., Nk (16)
Minimum optimization aim of the formula (11) for model, including cutting load measure financial cost referred to cutting load harmony Mark, the harmonious dimension difference of cutting load measure financial cost and node cutting load rate are excessive, it is impossible to directly compare coordinate it is excellent Change, therefore be normalized calculating, in formula, f1、f2Respectively cutting load measure financial cost, cutting load harmony index are returned One change value；ω1、ω2Respectively f1、f2Weight, value by policymaker according to real system service requirement depending on.NlIt is total for branch road Number；The current value upper limit is allowed for branch road l；Vi max、Vi minRespectively node i voltage upper and lower limit；σkFor node it is idle/it is active Scale factor, can have economize on electricity power factor to calculate and obtain.Pi、QiRespectively node i it is active, without work value；Vi、VjRespectively node I, j magnitude of voltage；θijFor node i, j phase difference of voltage；Gij、BijRespectively branch road lijElectric conductivity value, susceptance value.
In object function, two parts are directed to the active change of node i.e.：Different cutting load measures, node active power are different, According to the definition of object function, different target value is calculated.
Formula (12)-(16) are model constraints, and formula (16) represents that node, while active power is cut off, cuts off node Corresponding reactive power, after taking cutting load measure (change of node active power), becomes according to formula (16) calculate node reactive power Change, for node is active, reactive balance constraint, formula (14) is constrained for the branch current value upper limit, formula (15) is node for formula (12)-(13) Magnitude of voltage is constrained, wherein, the parameter such as the voltage of network each node, branch current can be according to formula (12), (13) by trend Calculating is obtained,, object function is counted with penalty function form model solution up-to-date style (14)-(15) are carried out.
Cutting load financial cost and harmony are required according to policymaker different, cut in setting model object function respectively negative The weighted value ω of lotus financial cost, node cutting load harmony index1、ω2, by calculating, provide under different weights for policymaker Optimal route overload joining-cutting load coordinates and optimizes control program.
In model, cutting load measure financial cost with cutting load harmony index derivation is：
Defining cutting load measure financial cost is：
Wherein, NkFor cutting load node total number；ΔPki、ciThe i-th type load cutting load amount of respectively node k, excision i-th Type load unit-economy cost coefficient.
Cutting load measure financial cost in formula (17) is normalized and is calculated：
Wherein, NbFor network load node total number；PkiFor the i-th type load total amount of node k.
It is node cutting load rate standard deviation to define cutting load harmony index definition：
Wherein,For k node cutting load rates；For all cutting load node cutting load rate mean values.
Formula (19) interior joint cutting load rate harmony index is normalized and is calculated：
Circuit overload connection cutting load coordination optimization control is complicated multiple constraint, a nonlinear optimal problem, therefore, this Inequality constraints in model is counted object function, and self adaptation of the application based on nonlinear weight with penalty function form by invention Particle swarm optimization algorithm is optimized solution.
Real system simulation analysis：By taking the electrical network real system of somewhere as an example, to the validity of institute's extracting method of the present invention with Feasibility is verified.Node cutting load percentage range is set as [0,50], one, two, three load cost coefficient of this area Respectively 1,0.4,0.1, this area's electrical network total load is 5469MW.Certain planning year this area's electric network composition is as shown in Figure 2.Accident There is overload in various degree in circuit 2-6,11-20, F-24 afterwards.
According to alternative cutting load node screening strategy, overall sensitivity calculating is carried out to all nodes in system, and by drop Sequence is arranged, and chooses front N of the overall sensitivity more than 0.1kIndividual node is included alternative cutting load set of node and (thinks that overall sensitivity is little Node in 0.1 does not possess competitiveness), higher front 10 nodes of overall sensitivity are only listed, as shown in table 1.
1 node overall sensitivity of table sorts
Using the adaptive particle swarm optimization Algorithm for Solving model based on nonlinear weight.Particle length is taken for Nk, particle Group's scale is 30, and integer coding is carried out to particle in population, as shown in table 2 below.
2 particle coding structure table of table
z1 z2 z3 … zNk
Wherein, zkFor kth place value, arbitrary integer in [0,50] is taken.If zkArbitrary integer in [1,50] is taken, is represented in node k Upper excision percentage is zkLoad；If zk0 is taken, represents that node k does not take cutting load measure.
In formula (11) object function, the asynchronous optimum results of weight coefficient value are as shown in table 3.
Simulation result under 3 different design conditions of table
By being compared analysis to above-mentioned four kinds of schemes, scheme 1 is represented and is only considered when cutting load measure is taken Cutting load measure economy cost, preferential excision sensitivity is high, the node that economic cost is low, causes cutting load rate difference between node Greatly, part economic cost it is relatively low node loss rate of load condensate be even as high as 40%, increase accident grading, to power grid enterprises bring compared with High policy fine, while causing bad social influence；Scheme 2 is represented when cutting load measure is taken, is compared using node etc. Example excision principle, the method that simulation result shows cause that system cutting load amount is excessive, and system operation is less economical, in real system It is inadvisable；Scheme 3 and scheme 4 using coordinating economic control strategy, both of which can take into account cutting load economy cost with it is balanced Property require, can effective control node cutting load rate, wherein scheme 4 more emphasizes tactful economy, during using scheme 4, system Node maximum cutting load rate is 17%, and contrast scheme 3 understands, scheme 4 can meet node cutting load rate it is less under the premise of, System cutting load amount is substantially reduced, optimal case is obtained.
Above-mentioned simulation result shows that scheme 3,4 can be obviously improved accident grading using coordination optimization control strategy, again may be used To ensure the economic feasibility of scheme.In real system, should be required according to system operation, carry out Scheme Choice.
Although the above-mentioned accompanying drawing that combines is described to the specific embodiment of the present invention, not to present invention protection model The restriction enclosed, one of ordinary skill in the art should be understood that on the basis of technical scheme those skilled in the art are not The various modifications made by needing to pay creative work or deformation are still within protection scope of the present invention.

Claims (10)

1. the circuit overload connection cutting load coordination optimizing method based on comprehensive sensitivity analysis, is characterized in that, comprise the following steps：
Node is screened：Network node admittance matrix and incidence matrix are formed, is defined with incidence matrix according to nodal voltage equation, is obtained Sensitivity relation between out branch electric current and node Injection Current；
Meter is set up according to the sensitivity relation between branch current and node Injection Current and branch road transships the node synthesis of degree Sensitivity index, and alternative cutting load node screening is carried out by the calculating of node overall sensitivity index；
Scheme is searched for：Coordinate cutting load financial cost and node cutting load two targets of harmony, cut with financial cost and node negative The minimum optimization aim of weighted value of lotus harmony index, sets up circuit overload connection cutting load coordination optimization Controlling model；
It is different according to requiring to cutting load financial cost and harmony, set the weight of cutting load financial cost and harmonious index Value, the calculating circuit overload connection cutting load coordinate and optimize Controlling model under different weights, optimal route overload joining-cutting load Coordination optimization control program, finally draws cutting load node.
2. the circuit overload based on comprehensive sensitivity analysis as claimed in claim 1 joins cutting load coordination optimizing method, and which is special Levying is, is defined with incidence matrix according to nodal voltage equation, show that the sensitivity between branch current and node Injection Current is closed It is to be：
I · = YA T Y n - 1 I · n
Wherein,For branch current column vector；Y is Nl×NlRank branch admittance matrix,For the inverse matrix of bus admittance matrix, For node Injection Current column vector, node is flowed into for just, it is negative to flow out, ATFor incidence matrix transposition.
3. the circuit overload based on comprehensive sensitivity analysis as claimed in claim 1 joins cutting load coordination optimizing method, and which is special Levying is, takes cutting load measure to eliminate ability to the excessively stream of all overload branch roads according to node k and defines node k sensitivity indexs, Consider that branch road overload degree is different, define node overall sensitivity index：
λ ′ k = Σ j = 1 N s δ j Y j A T Y n - 1
δ j = ( I j I max j ) 2
Wherein, δjFor branch road j weighted values, for representing the overload degree of branch road j；Ij、Respectively branch road j actual current values With the permission current value upper limit, NsTo transship branch road total number, YjFor branch admittance battle array jth row element, the i.e. self-admittance of branch road j with Transadmittance, ATFor incidence matrix transposition,For the inverse matrix of bus admittance matrix.
4. the circuit overload based on comprehensive sensitivity analysis as claimed in claim 3 joins cutting load coordination optimizing method, and which is special Levying is, node overall sensitivity index is normalized calculating, for intuitively express it is all overload branch roads a certain node is cut it is negative The overall sensitivity of lotus measure.
5. the circuit overload based on comprehensive sensitivity analysis as claimed in claim 4 joins cutting load coordination optimizing method, and which is special Levying is, when carrying out alternative cutting load node and screening, the node for taking node overall sensitivity value more than setting value is included and alternatively cut Load bus collection.
6. the circuit overload based on comprehensive sensitivity analysis as claimed in claim 1 joins cutting load coordination optimizing method, and which is special Levying is, cutting load measure financial cost is：
f 1 ′ = Σ k = 1 N k Σ i = 1 3 c i ΔP k i
Wherein, NkFor cutting load node total number；ΔPki、ciThe i-th type load cutting load amount of respectively node k, the i-th class of excision are born Lotus unit-economy cost coefficient；
Also include the step of being normalized calculating after calculating to cutting load measure financial cost.
7. the circuit overload based on comprehensive sensitivity analysis as claimed in claim 1 joins cutting load coordination optimizing method, and which is special Levying is, the cutting load harmony index definition is node cutting load rate standard deviation：
Wherein, NkFor cutting load node total number,For k node cutting load rates；It is average for all cutting load node cutting load rates Value, PkiFor the i-th type load total amount of node k, Δ PkiFor the i-th type load cutting load amount of node k.
8. the circuit overload based on comprehensive sensitivity analysis as claimed in claim 7 joins cutting load coordination optimizing method, and which is special Levying is, the step of node cutting load rate harmony index is normalized.
9. the circuit overload based on comprehensive sensitivity analysis as claimed in claim 1 joins cutting load coordination optimizing method, and which is special Levying is, circuit overload connection cutting load coordination optimization Controlling model is：
Min f=ω1f1+ω2f2
Wherein, ω1、ω2Respectively cutting load financial cost, the weighted value of node cutting load harmony index, are normalized meter Cutting load measure financial cost f after calculation1, it is normalized node cutting load rate harmony index f after calculating2。
10. the circuit overload based on comprehensive sensitivity analysis as claimed in claim 9 joins cutting load coordination optimizing method, and which is special Levying is, the constraints of circuit overload connection cutting load coordination optimization Controlling model includes that trend is active, reactive balance, node voltage And line electricity flow valuve is not out-of-limit.
CN201610978458.3A 2016-11-04 2016-11-04 Circuit overload based on comprehensive sensitivity analysis joins cutting load coordination optimizing method Active CN106532717B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610978458.3A CN106532717B (en) 2016-11-04 2016-11-04 Circuit overload based on comprehensive sensitivity analysis joins cutting load coordination optimizing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610978458.3A CN106532717B (en) 2016-11-04 2016-11-04 Circuit overload based on comprehensive sensitivity analysis joins cutting load coordination optimizing method

Publications (2)

Publication Number Publication Date
CN106532717A true CN106532717A (en) 2017-03-22
CN106532717B CN106532717B (en) 2018-12-04

Family

ID=58349626

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610978458.3A Active CN106532717B (en) 2016-11-04 2016-11-04 Circuit overload based on comprehensive sensitivity analysis joins cutting load coordination optimizing method

Country Status (1)

Country Link
CN (1) CN106532717B (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106953320A (en) * 2017-03-31 2017-07-14 国网江苏省电力公司 A kind of transmission line of electricity overload controlling method based on UPFC
CN107370149A (en) * 2017-08-25 2017-11-21 广州供电局有限公司 Determine the method and system of power system load reduction
WO2018201703A1 (en) * 2017-05-05 2018-11-08 国电南瑞科技股份有限公司 Adaptive equipment overload emergency control method based on automatic measure sensitivity identification
CN109473991A (en) * 2018-02-12 2019-03-15 国网浙江省电力公司嘉兴供电公司 A kind of area power grid Optimal Load control method considering regional balance
CN112994025A (en) * 2021-02-05 2021-06-18 山东大学 Power grid load removal method and system based on linear programming algorithm
CN113363991A (en) * 2021-07-02 2021-09-07 合肥工业大学 Tidal current control method based on comprehensive sensitivity
CN112994025B (en) * 2021-02-05 2024-04-19 山东大学 Power grid load removal method and system based on linear programming algorithm

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1083023A2 (en) * 1999-09-10 2001-03-14 Yoshiaki Kakino Numerical control apparatus for nc machine tool
KR100772984B1 (en) * 2006-04-21 2007-11-02 고수동 electric socket with breaking power supling for loads
CN101299530A (en) * 2007-11-08 2008-11-05 江苏省电力公司 Intelligent joining-cutting method for large electrified wire netting equipment overloading on-line sensitiveness
CN101702524A (en) * 2009-11-12 2010-05-05 华北电力大学 Generator and control method for load adjustment based on optimization model

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1083023A2 (en) * 1999-09-10 2001-03-14 Yoshiaki Kakino Numerical control apparatus for nc machine tool
KR100772984B1 (en) * 2006-04-21 2007-11-02 고수동 electric socket with breaking power supling for loads
CN101299530A (en) * 2007-11-08 2008-11-05 江苏省电力公司 Intelligent joining-cutting method for large electrified wire netting equipment overloading on-line sensitiveness
CN101702524A (en) * 2009-11-12 2010-05-05 华北电力大学 Generator and control method for load adjustment based on optimization model

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
徐岩等: "基于功率灵敏度的线路过载划区域紧急控制策略", 《电工技术学报》 *
程临燕等: "基于综合灵敏度分析的快速控制算法研究", 《电力自动化设备》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106953320A (en) * 2017-03-31 2017-07-14 国网江苏省电力公司 A kind of transmission line of electricity overload controlling method based on UPFC
CN106953320B (en) * 2017-03-31 2019-05-07 国网江苏省电力公司 A kind of transmission line of electricity overload controlling method based on UPFC
WO2018201703A1 (en) * 2017-05-05 2018-11-08 国电南瑞科技股份有限公司 Adaptive equipment overload emergency control method based on automatic measure sensitivity identification
CN107370149A (en) * 2017-08-25 2017-11-21 广州供电局有限公司 Determine the method and system of power system load reduction
CN107370149B (en) * 2017-08-25 2020-04-14 广州供电局有限公司 Method and system for determining load reduction amount of power system
CN109473991A (en) * 2018-02-12 2019-03-15 国网浙江省电力公司嘉兴供电公司 A kind of area power grid Optimal Load control method considering regional balance
CN109473991B (en) * 2018-02-12 2022-04-19 国网浙江省电力公司嘉兴供电公司 Regional power grid optimal load control method considering regional balance
CN112994025A (en) * 2021-02-05 2021-06-18 山东大学 Power grid load removal method and system based on linear programming algorithm
CN112994025B (en) * 2021-02-05 2024-04-19 山东大学 Power grid load removal method and system based on linear programming algorithm
CN113363991A (en) * 2021-07-02 2021-09-07 合肥工业大学 Tidal current control method based on comprehensive sensitivity

Also Published As

Publication number Publication date
CN106532717B (en) 2018-12-04

Similar Documents

Publication Publication Date Title
CN106532717A (en) Comprehensive sensitivity analysis-based circuit overload load-shedding coordination optimization method
CN103258299B (en) A kind of many direct currents concentrate the receiving end Net Frame of Electric Network optimization method of feed-in
CN104463701B (en) A kind of distribution system and the coordinated planning method of charging electric vehicle network
CN102663515B (en) A kind of method for optimizing of extra-high voltage transformer substation location
CN110969347B (en) Power transmission network structure morphology assessment method
CN103020853B (en) Method for checking short-term trade plan safety
CN107316125A (en) A kind of active distribution network economical operation evaluation method based on economical operation domain
CN110896231B (en) Distributed photovoltaic capacity receiving calculation method and system for power distribution network in poverty alleviation area
CN105514990B (en) The utilization rate of electric transmission line of comprehensive economy and security improves device and method
CN104866919A (en) Multi-target planning method for power grid of wind farms based on improved NSGA-II
CN102738780A (en) Method of optimal configuration of superconducting fault current limiter based on sensitivity technology
CN103870695A (en) Judgment method for voltage level of high power accessing power grid
CN104505821A (en) Power grid operation mode optimizing method for controlling short circuit current level
CN115439000A (en) Power distribution network block division method considering wind-solar-load power uncertainty and correlation
CN103020729A (en) Comprehensive evaluation method for short-circuit current limitation of Gansu power grid
CN104484555A (en) Method for evaluating maximum power supply capability of 220kV self-healing looped network
CN105119279A (en) Distributed power supply planning method and system thereof
CN104484832B (en) The method of assessment 220KV handle net net capability
CN104967121A (en) Large-scale electric power system node load flow computing method
CN103887823A (en) Micro-grid connection position selection method based on fuzzy hierarchical analysis
CN106505576B (en) It is a kind of based on Link-path urban high voltage distribution network load turn for method
CN104410065B (en) A kind of Multiobjective Decision Making Method of receiving end electrical network limiting short-circuit current
CN103560497A (en) Short-circuit current limiting method based on power grid topological adjustment
Dahej et al. Optimal allocation of SVC and TCSC for improving voltage stability and reducing power system losses using hybrid binary genetic algorithm and particle swarm optimization
CN105162143A (en) Method for determining optimal trading capacity between large regional grids

Legal Events

Date Code Title Description
C06 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